In work published in The European Physical Journal C, theorists Quan-Yun Guo and Dian-Yong Chen, based at Southeast University and Lanzhou University in China, have calculated how the three Pc states could be produced in the process e+ e- → p p̄ J/ψ, where an electron and a positron annihilate into a virtual photon that fragments into an antiproton and a charmed pentaquark, or into a proton, an antiproton, and the charmonium meson J/ψ. Their calculations, built on an effective Lagrangian framework, predict cross sections small enough to challenge experimenters but tantalizingly within reach of the Super Tau-Charm Facility (STCF), a machine designed to operate with center-of-mass energies between 2 and 7 GeV and a peak luminosity of 5.0 × 10^34 cm^-2 s^-1—more than fifty times the event yield of its predecessor, BEPCII.

The three Pc states at the heart of the study were all carved out of the decay Λb^0 → J/ψ K^- p, first measured by LHCb. In 2019, with vastly improved statistics, the collaboration reported a narrow structure at 4312 MeV with a significance of 7.3 sigma, and revealed that the earlier broad Pc(4450) signal was actually two overlapping resonances, Pc(4440) and Pc(4457). The Particle Data Group now lists the Pc(4312) mass as 4311.9 MeV with a width of about 10 MeV, Pc(4440) at 4440 MeV with a width of roughly 21 MeV, and Pc(4457) at 4457.3 MeV with a width of about 6.4 MeV. What makes these numbers so provocative is where they sit: the mass of Pc(4312) lies almost exactly at the threshold for a charmed baryon Σc bound to an anti-charmed meson D̄, while the two heavier states coincide with the fine structure expected of a Σc D̄* molecule with total spin 1/2 and 3/2, respectively.

That near-threshold positioning has pushed a large fraction of the community toward the hadronic molecule interpretation, in which the pentaquark is not a compact bag of five valence quarks but a loosely bound pairing of two color-singlet hadrons—an arrangement more like a deuteron than a proton. One-boson-exchange calculations of the Σc D̄() interaction, quasipotential Bethe-Salpeter equation studies, contact-range effective field theory incorporating heavy-quark spin symmetry, and QCD sum rule analyses have all reached broadly compatible conclusions: Pc(4312) is a Σc D̄ molecule with quantum numbers J^P = 1/2^-, while Pc(4440) and Pc(4457) are Σc D̄ molecules with J^P = 1/2^- and 3/2^-, respectively. Alternative interpretations persist, however, including compact diquark pentaquarks and, most stubbornly, kinematic threshold cusps—sharp peak structures that arise when the opening of a new production channel mimics a resonance without any genuine bound state existing at all.

The new study attacks the problem from the production side, asking not what the states decay into, but how efficiently they can be created. The mechanism is elegant in its simplicity: an electron and positron annihilate into a virtual photon, which through the vector meson dominance mechanism couples to the J/ψ, and the J/ψ in turn couples electromagnetically to the p p̄ Pc vertex. Because the pentaquarks have so far only been reconstructed in J/ψ p final states, the authors focus exclusively on the J/ψ-mediated contribution, fixing the electromagnetic couplings from the strong couplings g_PcψN through vector meson dominance. Those strong couplings, in turn, are pinned down by assuming that the J/ψ p channel accounts for 10 percent of each state’s total width—a branching fraction consistent with several molecular-model predictions, though other frameworks have suggested values anywhere from a few percent to tens of percent, and the results scale linearly with this assumption.

The headline numbers come at a center-of-mass energy of √s = 6 GeV, where the predicted cross sections are 27.8 fb (with an uncertainty band of +42.2/-17.9 fb) for e+ e- → p̄ Pc(4312), 38.7 fb (+58.6/-25.1 fb) for p̄ Pc(4440), and just 1.44 fb (+2.18/-0.93 fb) for p̄ Pc(4457). A femtobarn is an extraordinarily small unit—one fb corresponds to a single event in 10^39 interactions—and these values lie safely below the theoretical upper limit of roughly 0.1 pb established in an earlier independent analysis of the same process. The dominant uncertainty stems from a form-factor parameter Λr, which encodes the internal structure of the hadrons and their off-shell behavior; the authors vary it between 2.0 and 3.0 GeV, a range calibrated against analogous studies of the Λc(2910) and Λc(2940) states, which are themselves interpreted as D*N molecular partners of the pentaquarks.

When the decay of the pentaquarks into J/ψ p is folded in, the full four-body process e+ e- → p p̄ J/ψ acquires a total predicted cross section of 14.7 fb (+22.2/-9.50 fb) at 6 GeV—spanning nearly an order of magnitude across the allowed Λr range, a sobering reminder of the model dependence inherent in such calculations. Crucially, the predictions remain consistent with Belle II’s recent first measurement of this very process, performed via initial-state radiation from threshold up to 7 GeV; no clear structure emerged in the measured p p̄ J/ψ cross section, and the new estimates fall comfortably below the reported upper limits. Among the individual contributions, the Pc(4440)/P̄c(4440) channel dominates, contributing about 6.8 fb at 7 GeV—roughly thirty times the yield from Pc(4457).

The most experimentally valuable prediction concerns the J/ψ p invariant mass spectrum. At √s = 6 GeV, the calculations reveal narrow peaks near 4.31, 4.44, and 4.46 GeV, corresponding to the three pentaquark states, with Pc(4440) providing the strongest signal. At 7 GeV, the picture sharpens further: the prominent structure between 4.42 and 4.46 GeV originates overwhelmingly from Pc(4440) rather than its partner Pc(4457), whose contribution is visible only as a weak signal near 4.45–4.46 GeV. Moreover, the overlap between Pc states and their antiparticle counterparts becomes negligible at this energy, meaning experimenters could disentangle the pentaquark signals without significant interference from P̄c contributions—a clean separation that LHCb’s fixed-target-style decays could never offer.

The study also delivers a diagnostic tool for one of the thorniest open questions in the field: the quantum numbers of the pentaquarks. Helicity-angle distributions of the J/ψ meson turn out to be decisive. For Pc(4312) and Pc(4440), both assigned J^P = 1/2^-, the predicted distributions are isotropic, appearing as flat lines in cos θ. For Pc(4457) with J^P = 3/2^-, the distribution takes on a distinctly concave shape, following an A + B cos²θ pattern with maxima at cos θ = ±1 and a minimum at cos θ = 0. Such a signature, measurable at STCF, would constitute a direct probe of the angular-momentum content of these states—independent of the molecular interpretation that underpins the cross-section estimates.

Translating the numbers into raw event counts, the authors estimate that if the STCF runs 200 days per year at √s = 6 GeV, it could accumulate between roughly 5,000 and 32,000 J/ψ p p̄ events per year—about (1.3^+1.9_-0.8) × 10^4—based on the 10 percent branching fraction assumption. Even under the more conservative 3 percent assumption, the yield remains on the order of 10^3 events annually, sufficient to reconstruct the invariant mass spectra and angular distributions that would confirm or refute the molecular picture. In an era when the LHCb experiment continues to enrich the hidden-charm pentaquark family with strange partners such as Pcs(4459) and Pcs(4380), the prospect of a dedicated e+ e- facility systematically scanning their production properties represents a complementary and potentially decisive line of attack.

The broader stakes extend well beyond a single set of resonances. Hadronic molecules, if confirmed as the explanation for the Pc states, would establish that quantum chromodynamics permits stable bound configurations of two heavy hadrons—insights that feed directly into our understanding of neutron-star matter, where strange and charmed baryon interactions may govern the equation of state at extreme densities. Conversely, if the peaks turn out to be threshold cusps, as some analyses of Pc(4312) and the related Pc(4337) structure have argued, the pentaquark saga would become a textbook demonstration of how final-state interactions can masquerade as resonances. Either outcome would reshape the growing zoo of exotic hadrons that began with the X(3872) in 2003 and now spans tetraquarks, pentaquarks, and doubly heavy states. The new calculations provide exactly the quantitative bridge between theory and experiment that the field needs: concrete, testable cross-section predictions, distinguishable spectral shapes, and angular signatures, all tailored to an experimental program poised to begin. Whether the STCF will ultimately photograph these five-quark shadows in the clean light of electron–positron annihilation now depends on turning femtobarn-scale predictions into measured counts—one collision at a time.