REVIEW 4 major objections 5 minor 20 references
The influence of nuclear short range correlations on sub-threshold particle production in proton-nucleus collisions
T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read Short-range nuclear proton-neutron pairs, via the deuteron-like high-momentum tail, boost above-threshold nucleon-nucleon collisions by up to 1000-fold and reproduce measured p+Nb data.
desk verdict New and plausible SRC-based mechanism for sub-threshold production, but the on-shell treatment in Eq. (2) leaves the size of the claimed enhancement unproven. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the single-nucleon momentum distribution $n(k)=n_0(k)+n_1(k)$ supplied by the spectral-function factorization of Ref. [19], split into a Fermi part $n_0$ and a deuteron-like SRC tail $n_1$ that extends to $k\simeq 5$ fm$^{-1}$. The mechanism is the normalized collision-energy spectrum of Eq. (2), $dN/d\sqrt{s}=\int d^3k\,(dn/d^3k)\,\delta\big(\sqrt{(k^\mu+p^\mu)^2}-\sqrt{s}\big)$, which converts target-nucleon momenta into a distribution of available center-of-mass energies. Integrating this spectrum from each hadron's threshold gives the probability $P(\sqrt{s}\ge\sqrt{s_{\rm thr}})$, and particle yields are obtained by folding the spectrum with $\sigma_{pp\to X}$ and normalizing to the Glauber number of binary collisions.
What would settle it
A decisive test is to measure the $\Xi^-$ (or, better, $\Lambda_c$) excitation function in p+Fe and p+Pb collisions at nominal $\sqrt{s_{NN}}$ values 0.1 to 0.5 GeV below threshold: the SRC prediction from Eq. (2) with the Ref. [19] momentum distribution rises steeply and has a specific per-collision target dependence, so a mismatch in slope or scaling would falsify the mechanism. A complementary check is to recompute the same spectra with the full spectral function $P(k,E)$ including removal energy; if the enhancement disappears, the on-shell treatment at Eq. (2) is the cause.
Extended reading notes
Core claim
The paper's central discovery is that the high-momentum, deuteron-like tail of the nuclear momentum distribution, the signature of short-range correlated $pn$ pairs, acts as a reservoir of high-$\sqrt{s}$ nucleon-nucleon collisions in p+A reactions. In the factorization of Ref. [19], $n(k)=n_0(k)+n_1(k)$, the tail $n_1$ extends to $k\simeq 5$ fm$^{-1}$, far beyond the Fermi momentum, and feeding this distribution into the collision spectrum of Eq. (2) raises the probability that an individual nucleon-nucleon collision sits above an elementary production threshold by up to three orders of magnitude compared with a Fermi gas. Weighting this spectrum by the pp production cross section and normalizing to a Glauber binary-collision count reproduces the measured $\Xi^-$ multiplicity in p+Nb at 61 MeV below threshold. The authors conclude that SRC kinematics, rather than multi-step processes, is the natural explanation for sub-threshold strangeness and charm production in proton-nucleus collisions.
Load-bearing premise
The calculation assumes that a nucleon knocked out of a correlated pair behaves as a free particle whose full momentum is available as kinetic energy, with no energy cost for removing it from the nucleus, and it benchmarks this against a single p+Nb data point only 61 MeV below threshold, where the SRC and Fermi-gas predictions differ by about the experimental uncertainty.
Editorial extensions
If this is right
- Deep sub-threshold yields of $\Lambda$, $\phi$, $\Xi$, $\Omega$, $J/\psi$, and $\Lambda_c$ in p+A collisions are set by SRC kinematics; the heavier the hadron, the farther below its elementary threshold a fixed fraction of collisions remains, reaching about 1.4 GeV below threshold for $\Lambda_c$ versus 0.47 GeV for $\Lambda$ in p+Pb at the 0.1% benchmark.
- Because the SRC tail is dominated by $pn$ pairs and is nearly the same in different nuclei, sub-threshold production becomes a nearly isospin-independent probe, and the target-mass dependence of the yield can separate SRC from multi-step rescattering.
- With a collision-energy reach of $\sqrt{s}=4.92$ GeV for Au+Au, $J/\psi$ mesons can still be produced about 1 GeV below threshold, making sub-threshold charm accessible.
- The measured p+Nb $\Xi^-$ multiplicity at $\sqrt{s_{NN}}=3.18$ GeV, 61 MeV below threshold, is reproduced when SRC are included, whereas a pure Fermi gas alone underpredicts it.
Reading between the lines
- If the same deuteron-like tail governs deep sub-threshold kinematics, production of even heavier exotic states such as $\Xi_{cc}$ or $X/Y/Z$ in p+A collisions should also be enhanced, and their much higher thresholds would make them a sharper test of the mechanism than $\Xi^-$.
- A direct way to probe the paper's on-shell approximation is to repeat the calculation with the full spectral function $P(k,E)$ including removal energy; the difference would quantify how much of the reported enhancement survives when the removal-energy cost is accounted for.
- Because the same correlated-pair tail is invoked to explain the EMC effect, sub-threshold hadron production in p+A offers a complementary, strong-interaction window on short-range correlations with different systematics than electron-scattering measurements.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes that short-range correlated (SRC) proton-neutron pairs in the target nucleus, represented by the deuteron-like tail n1(k) in the single-nucleon momentum distribution n(k)=n0(k)+n1(k), strongly increase the available invariant mass sqrt(s) in p+A collisions. Using the Ciofi degli Atti-Simula parametrization, the authors compute the normalized collision spectrum dN/dsqrt(s) via Eq. (2) and report that the fraction of nucleon-nucleon collisions above the elementary production threshold is enhanced by up to three orders of magnitude relative to a pure Fermi-gas target at deep subthreshold energies. They apply the same spectrum, multiplied by an unspecified pp -> Xi^- X cross section, to estimate the Xi^- multiplicity in p+Nb at sqrt(s_NN)=3.18 GeV and compare with the HADES measurement. The paper argues that SRC kinematics, rather than multi-step processes, dominate deep subthreshold strange and charmed hadron production and that FAIR can test this hypothesis.
Significance. If the mechanism is quantitatively correct, it is an important and simple explanation for sub-threshold strange-baryon production and a strong motivation for sub-threshold charm searches at FAIR. The approach is not circular: the SRC tail is taken from independent electron-scattering data via Ref. [19], and no parameter is fitted to the HADES point, which is a genuine strength. The central quantitative claim, however, is conditional on a spectral-function approximation that is not yet quantified, and the only existing data comparison is too close to threshold to discriminate SRC from ordinary Fermi motion. The paper is therefore a promising contribution whose key magnitude and validation need further substantiation.
major comments (4)
- [Section II, Eq. (2)] Eq. (2) evaluates the invariant mass using a target four-momentum k^mu = (sqrt(m_N^2+k^2), k), placing every target nucleon on the mass shell. The Atti-Simula input is a spectral function P(k,E), and the momentum distribution n(k)=int dE P(k,E) discards the removal-energy dependence. For the SRC tail around |k| = 0.4-1 GeV/c that dominates the deep-subthreshold ratios in Fig. 2, an energy shift of tens to a few hundred MeV is not negligible: it changes sqrt(s) by an amount comparable to the threshold deficits the tail is supposed to overcome. Footnote 3 checks only the mean sqrt(s), which does not constrain the tail. The authors should either use the full spectral function P(k,E), adopt a well-defined off-shell prescription for the target nucleon energy, or provide a quantitative estimate demonstrating that the energy dependence of P(k,E) is negligible for the high-momentum tail. Without this, the magnitude of the claimed enhancement is not established.
- [Section III.B, Fig. 4] The benchmark against HADES is too weak to support the central claim. The data point lies only 61 MeV below the elementary Xi threshold, where the SRC and Fermi curves differ by an amount that the authors themselves describe as 'on a similar magnitude as the experimental uncertainty.' Consequently the agreement in Fig. 4 does not distinguish the SRC mechanism from a standard Fermi-gas treatment and does not validate the deep-subthreshold enhancement that is the paper's main quantitative result. The abstract and conclusion should not present this comparison as confirming SRC without this qualification, and the authors should either add a more constraining comparison or clearly reframe the HADES point as a consistency check only.
- [Section III.B, Fig. 4] The multiplicity calculation is not reproducible as reported. The text states that the sqrt(s) spectrum is multiplied by the production cross section sigma_{pp -> Xi^- X}(sqrt(s)), but this cross section is not specified: no parametrization, reference, or numerical source is given. Since the integral weights the same high-sqrt(s) tail that produces the claimed enhancement, the comparison in Fig. 4 depends directly on this input. Please provide the functional form and normalization of sigma_{pp -> Xi^- X} and, ideally, a sensitivity check with a different parametrization.
- [Section III.A, Fig. 2] The reported enhancement factor of up to 10^3 is presented without any uncertainty or sensitivity analysis. The deep-subthreshold ratios are tail-dominated, and the magnitude of the SRC tail depends on the pair fraction (about 20% from the CLAS data) and on the upper momentum cutoff used in Eq. (2), neither of which is specified in the manuscript. The authors should show how the ratios in Fig. 2 and the multiplicity in Fig. 4 change under reasonable variations of these inputs, for example a +/-20% change in the SRC normalization or a change in the maximum momentum included in the integral.
minor comments (5)
- [Section III.A] In the text and in the expression sqrt(s_NN)|_{P=0.1%}^{SCR}, the abbreviation 'SCR' appears twice; it should be 'SRC'.
- [Introduction] The passage 'Inthisletterwestudytheinfluenceoftheseshortrangecorrelations...' has missing spaces and should be corrected.
- [Section II, Eq. (2)] The notation d^3n(k)/d^3k combined with an integral over d^3k is confusing; since n(k) is already the momentum-space density, the measure should be written as n(k)d^3k, or the differential distribution should be explicitly defined.
- [Section II] The parametrizations of n0(k) and n1(k) from Ref. [19] are not displayed. Including the functional form, or at least the parameter values used for Pb and Fe, would substantially improve reproducibility and would make the tail sensitivity easier to judge.
- [Fig. 4 caption] The vertical axis is labeled 'multiplicity,' but the normalization (per event, per binary collision, or absolute yield) is not stated in the caption; the text mentions a Glauber normalization, but the caption should make this explicit.
Circularity Check
No significant circularity: the SRC enhancement and HADES comparison are computed from external nuclear-structure inputs and an external pp cross section, not from the benchmark itself.
full rationale
Eq. (2) computes the NN center-of-mass spectrum from the independently published Atti-Simula momentum distribution n(k)=n0+n1 (Ref. [19], itself constrained by CLAS lepton-nucleus data) and from the projectile momentum; Eq. (3) integrates that externally specified spectrum. The up-to-10^3 SRC/Fermi ratios in Fig. 2 and the HADES comparison in Fig. 4 follow by folding this spectrum with an external pp->XiX cross section and a Glauber normalization; no parameter is fitted to the HADES point. The only self-citation (Ref. [13], cited for the statement that the spectral-function model 'agrees well with the experimental results') is motivational rather than load-bearing, since n(k) and the threshold integrals are defined by the Atti-Simula parametrization, not by that citation. The on-shell four-momentum choice in Eq. (2) and the unspecified sigma_pp input are correctness and robustness caveats, but they are not cases where an output equals an input by construction. No circular step is exhibited.
Assumptions & free parameters
free parameters (2)
- SRC pair fraction (high-momentum tail normalization) =
approximately 20% of nucleons
- Deuteron-like tail shape of n1(k) =
not given numerically
assumptions (4)
- domain assumption The single-nucleon momentum distribution factorizes into a Fermi part and a deuteron-like SRC part, n(k)=n0(k)+n1(k).
- domain assumption Target nucleons are treated as on-shell with fixed mass m_N in Eq. (2), with removal energy E omitted from the spectral function.
- domain assumption The elementary pp production cross section sigma_pp->Xi^-X is used for pN collisions inside the nucleus.
- domain assumption Multi-step processes and resonance decay channels are neglected as the production mechanism.
Cite this review
Pith. "Pith review of The influence of nuclear short range correlations on sub-threshold particle production in proton-nucleus collisions." pith.science (2026). https://pith.science/paper/S2WZ6W5E
@misc{pith2026250603962,
author = {Pith},
title = {Pith review of: The influence of nuclear short range correlations on sub-threshold particle production in proton-nucleus collisions},
year = {2026},
howpublished = {\url{https://pith.science/paper/S2WZ6W5E}},
note = {Machine review of arXiv:2506.03962}
}
abstract
The apparent production of (multi-)strange baryons and mesons at sub-threshold energies in proton-heavy ion collisions is a consequence of short range correlations (SRC), which have recently been observed in lepton-nucleus scattering. They may enhance the available center of mass energy of individual nucleon-nucleon collisions and allow, therefore, for sub-threshold particle production in proton-nucleus collisions. Calculations demonstrate that SRC enhance the probability for particle production at nominal sub-threshold energies up to a factor of $\times 10^3$ as compared to a simple Fermi gas model. We benchmark the idea by calculating the $\Xi^-$ multiplicity in nominal sub-threshold p+Nb collisions which compare well with the data measured by the HADES collaboration. These findings are of prime relevance for upcoming experiments at the FAIR facility, especially for the study of charmed hadrons.
Figures
Reference graph
Works this paper leans on
-
[19]
C. Ciofi degli Atti and S. Simula, Phys. Rev. C53, 1689 (1996), arXiv:nucl-th/9507024
arXiv 1996
-
[1]
Sub-threshold phi-meson yield in central 58Ni+58Ni collisions
A. Mangiarotti et al. (FOPI), Nucl. Phys. A 714, 89 (2003), arXiv:nucl-ex/0209012
work page Pith review arXiv 2003
-
[2]
Sub-threshold production of $\Sigma$(1385) baryons in Al+Al collisions at 1.9$A$ GeV
X. Lopezet al. (FOPI), Phys. Rev. C76, 052203 (2007), arXiv:0710.5007 [nucl-ex]
work page Pith review arXiv 2007
-
[3]
Deep sub-threshold $\Xi^-$ production in Ar+KCl reactions at 1.76A GeV
G. Agakishiev et al. (HADES), Phys. Rev. Lett. 103, 132301 (2009), arXiv:0907.3582 [nucl-ex]
work page Pith review arXiv 2009
- [4]
-
[5]
Deep sub-threshold {\phi} production and implications for the K+/K- freeze-out in Au+Au collisions
J. Adamczewski-Musch et al. (HADES), Phys. Lett. B 778, 403 (2018), arXiv:1703.08418 [nucl-ex]
work page Pith review arXiv 2018
-
[6]
J. Adamczewski-Musch et al. (HADES), Phys. Lett. B 793, 457 (2019), arXiv:1812.07304 [nucl-ex]
arXiv 2019
-
[7]
Subthreshold Xi- Production in Collisions of p(3.5 GeV)+Nb
G.Agakishiev et al.,Phys.Rev.Lett. 114,212301(2015), arXiv:1501.03894 [nucl-ex]
work page Pith review arXiv 2015
Show all 20 references
-
[8]
Barthet al
R. Barthet al. (KaoS), Phys. Rev. Lett.78, 4007 (1997)
1997
-
[9]
C. T. Sturm et al. (KAOS), Phys. Rev. Lett. 86, 39 (2001), arXiv:nucl-ex/0011001
2001 arXiv
-
[10]
Hartnack, H
C. Hartnack, H. Oeschler, and J. Aichelin, Phys. Rev. Lett. 96, 012302 (2006), arXiv:nucl-th/0506087
2006 arXiv
-
[11]
Hartnack, H
C. Hartnack, H. Oeschler, Y. Leifels, E. L. Bratkovskaya, and J. Aichelin, Phys. Rept. 510, 119 (2012), arXiv:1106.2083 [nucl-th]
2012 arXiv
-
[12]
Steinheimer and M
J. Steinheimer and M. Bleicher, J. Phys. G43, 015104 (2016), arXiv:1503.07305 [nucl-th]
2016 arXiv
-
[13]
T. Song, L. Tolos, J. Wirth, J. Aichelin, and E. Bratkovskaya, Phys. Rev. C 103, 044901 (2021), arXiv:2012.05589 [nucl-th]
2021 arXiv
-
[14]
Navaset al
S. Navaset al. (Particle Data Group), Phys. Rev. D110, 030001 (2024)
2024
-
[15]
K. I. Blomqvistet al., Phys. Lett. B424, 33 (1998)
1998
-
[16]
Duer et al
M. Duer et al. (CLAS), Nature 560, 617 (2018)
2018
-
[17]
Schmookler et al
B. Schmookler et al. (CLAS), Nature 566, 354 (2019), arXiv:2004.12065 [nucl-ex]
2019 arXiv
-
[18]
Duer et al
M. Duer et al. (CLAS), Phys. Rev. Lett. 122, 172502 (2019), arXiv:1810.05343 [nucl-ex]
2019 arXiv
-
[20]
Reichert, J
T. Reichert, J. Steinheimer, and M. Bleicher, Eur. Phys. J. A 61, 104 (2025), arXiv:2503.01613 [nucl-th]
2025 arXiv
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.