Recognition: unknown
Stability of parton distributions at high x: impact of nuclear and power corrections
Pith reviewed 2026-05-09 18:55 UTC · model grok-4.3
The pith
The u and d quark PDFs remain stable up to x ≈ 0.8 and describe DIS data down to W² = 3.5 GeV² when nuclear and higher-twist corrections are included.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Within the collinear factorization framework the JAM Bayesian Monte Carlo analysis shows that the u and d quark PDFs and their ratio are relatively stable up to x ≈ 0.8. The distributions continue to describe inclusive DIS data down to W² = 3.5 GeV² and Q² = m_c² once target-mass and higher-twist corrections are fitted. The higher-twist coefficients are positive and largely isospin independent. Nuclear data require nonzero isoscalar and isovector off-shell modifications to the nucleon PDFs, which in turn produce definite predictions for the ratio R_D of deuteron to isoscalar nucleon structure functions.
What carries the argument
The JAM Bayesian Monte Carlo global QCD fit that simultaneously varies the W cut, target-mass corrections, higher-twist coefficients, and nuclear wave functions plus off-shell nucleon PDF modifications for deuteron and A=3 targets.
If this is right
- The extracted u and d distributions can be used with greater confidence in high-x observables such as Drell-Yan production or neutrino scattering.
- Positive, isospin-independent higher-twist terms improve the description of DIS data at moderate Q² without spoiling the leading-twist PDFs.
- The required isoscalar and isovector off-shell corrections give a concrete, testable prediction for R_D that can be checked with future deuteron data.
- Lower W cuts become viable for PDF fits, allowing more inclusive measurements to constrain the large-x region.
Where Pith is reading between the lines
- The same off-shell mechanism may help quantify the EMC effect in heavier nuclei once the A=3 results are extended.
- If the stability persists, global fits could incorporate still lower-Q² data from future facilities to reduce uncertainties at x > 0.8.
- The isospin independence of the higher-twist terms suggests a common dynamical origin that could be tested in separate proton and neutron extractions.
Load-bearing premise
The collinear factorization framework remains valid at the low W and Q values included in the fit, and the chosen nuclear wave functions together with the off-shell PDF modifications capture all relevant nuclear effects without large residual biases.
What would settle it
A new high-precision measurement of the deuteron structure function ratio R_D at x > 0.7 that lies well outside the predicted band from the off-shell contributions, or a clear sign change in the fitted higher-twist coefficients when additional low-W data are added, would falsify the reported stability.
Figures
read the original abstract
We present a comprehensive new global QCD analysis of unpolarized parton distribution functions (PDFs) based upon proton, deuteron and $A\!=\!3$ data, including the latest inclusive deep-inelastic scattering (DIS) measurements from Jefferson Lab at high Bjorken-$x$. Using the JAM Bayesian Monte Carlo framework, we systematically explore the stability of the PDFs with respect to variations in the cuts on the invariant mass $W$ of the DIS final state, the implementation of target mass and higher twist corrections, as well as on the nuclear wave functions for the $A\!=\!2$ and 3 data. We find the $u$ and $d$ quark PDFs (and the $d/u$ ratio) are relatively stable up to $x \approx 0.8$, and able to describe DIS data down to $W^2=3.5$ GeV$^2$ and $Q^2=m_c^2$. Within the collinear factorization framework, the fitted higher twist corrections to DIS are found to be positive, and largely isospin independent. The description of the nuclear data also requires nonzero isoscalar and isovector nucleon off-shell PDF contributions, which gives specific predictions for the ratio, $R_D$, of deuteron to isoscalar nucleon structure functions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a global QCD analysis of unpolarized PDFs using proton, deuteron, and A=3 data, including recent JLab high-x DIS measurements. Employing the JAM Bayesian Monte Carlo framework, the authors systematically vary kinematic cuts on W, implementations of target-mass and higher-twist corrections, and nuclear wave functions for A=2,3. They report that the u and d quark PDFs (and d/u ratio) remain relatively stable up to x≈0.8 and describe the data down to W²=3.5 GeV² and Q²=m_c². Higher-twist corrections are found to be positive and largely isospin-independent, while nonzero isoscalar and isovector off-shell nucleon PDF modifications are required to fit the nuclear data, yielding predictions for the deuteron-to-isoscalar structure function ratio R_D.
Significance. If the separation of leading-twist PDFs from the fitted power corrections and nuclear effects holds, the work would be significant for high-x PDF phenomenology by enabling reliable inclusion of low-W data in global fits. The systematic Bayesian Monte Carlo exploration of model variations (cuts, nuclear wave functions, off-shell terms) is a clear strength, providing better uncertainty quantification than single-fit approaches. This could impact extractions of the d/u ratio at large x and applications involving nuclear corrections in light nuclei.
major comments (2)
- [§4] §4 (results on W² cuts and higher-twist fits): The central claim of PDF stability down to W²=3.5 GeV² rests on the fitted positive, isospin-independent higher-twist corrections fully absorbing deviations from collinear factorization. At these low invariant masses (near the resonance region) and Q²=m_c², resonance contributions and possible target-mass effects beyond the standard implementation could be partially absorbed into the HT coefficients or off-shell parameters, artificially stabilizing the high-x u, d, and d/u distributions. The Bayesian exploration of cuts does not explicitly demonstrate that such degeneracies are ruled out (e.g., via comparison to resonance-subtracted data sets or alternative HT parametrizations).
- [§5] §5 (nuclear corrections and R_D prediction): The specific predictions for R_D are generated from the same fitted isoscalar and isovector off-shell PDF contributions used to describe the A=2,3 data. This reduces R_D to a consistency check within the fitted model rather than an independent test, weakening its support for the overall framework's validity at high x.
minor comments (2)
- [Abstract] Abstract: The phrase 'relatively stable' is used for the PDFs up to x≈0.8; a quantitative statement (e.g., maximum variation in central values or uncertainty bands across W cuts) would improve clarity.
- [Figures] Figure captions (throughout): Several figures comparing PDF bands for different W cuts or nuclear models would benefit from explicit labels indicating the exact parameter variations shown.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below, providing clarifications on our methodology and indicating where revisions will be made.
read point-by-point responses
-
Referee: [§4] §4 (results on W² cuts and higher-twist fits): The central claim of PDF stability down to W²=3.5 GeV² rests on the fitted positive, isospin-independent higher-twist corrections fully absorbing deviations from collinear factorization. At these low invariant masses (near the resonance region) and Q²=m_c², resonance contributions and possible target-mass effects beyond the standard implementation could be partially absorbed into the HT coefficients or off-shell parameters, artificially stabilizing the high-x u, d, and d/u distributions. The Bayesian exploration of cuts does not explicitly demonstrate that such degeneracies are ruled out (e.g., via comparison to resonance-subtracted data sets or alternative HT parametrizations).
Authors: We appreciate the referee's concern about possible absorption of resonance or additional target-mass effects into the higher-twist and off-shell parameters. Our Bayesian Monte Carlo analysis systematically varies the W² cuts and HT implementations, showing that the extracted u and d PDFs (and d/u ratio) remain stable up to x≈0.8 across these choices, with the HT corrections remaining positive and largely isospin-independent. The framework quantifies uncertainties from these variations. We agree that explicit resonance-subtracted comparisons or alternative HT forms would provide stronger evidence against degeneracies and will add a dedicated discussion of this limitation, along with suggestions for future checks, in the revised manuscript. revision: partial
-
Referee: [§5] §5 (nuclear corrections and R_D prediction): The specific predictions for R_D are generated from the same fitted isoscalar and isovector off-shell PDF contributions used to describe the A=2,3 data. This reduces R_D to a consistency check within the fitted model rather than an independent test, weakening its support for the overall framework's validity at high x.
Authors: We acknowledge that the R_D predictions are derived from the off-shell parameters constrained by the A=2 and A=3 data. However, R_D is not an input to the fit but a specific, falsifiable output of the model for a different observable. This provides a non-trivial consistency test of the extracted isoscalar and isovector off-shell modifications. While we agree that fully independent experimental validation would strengthen the case, the current prediction still offers a concrete benchmark for the framework at high x. revision: no
Circularity Check
Fitted off-shell contributions to nuclear data presented as independent predictions for R_D
specific steps
-
fitted input called prediction
[Abstract]
"The description of the nuclear data also requires nonzero isoscalar and isovector nucleon off-shell PDF contributions, which gives specific predictions for the ratio, R_D, of deuteron to isoscalar nucleon structure functions."
Off-shell isoscalar and isovector nucleon PDF contributions are fitted parameters adjusted to reproduce the nuclear (deuteron and A=3) DIS data. R_D is the ratio of deuteron to isoscalar nucleon structure functions, which is directly computed from the same fitted nuclear model. Therefore the 'specific predictions' for R_D are the fitted model's output on the input data by construction, not an independent derivation.
full rationale
The paper performs a global QCD fit using JAM Bayesian Monte Carlo, determining u/d PDFs, higher-twist coefficients, and isoscalar/isovector off-shell nucleon PDF modifications simultaneously from proton, deuteron, and A=3 DIS data. The abstract states that nonzero off-shell contributions are required to describe the nuclear data and 'gives specific predictions for the ratio, R_D'. Since R_D is the deuteron-to-isoscalar-nucleon structure function ratio constructed from the same fitted nuclear model, this step reduces to a fitted output rather than an independent test. The core stability claim for PDFs up to x≈0.8 retains some independence through systematic cut variations and isospin-independent HT findings, preventing a higher score. No self-definitional equations or self-citation chains were identified in the provided text.
Axiom & Free-Parameter Ledger
free parameters (2)
- higher-twist corrections
- isoscalar and isovector off-shell PDF contributions
axioms (1)
- domain assumption Collinear factorization framework holds down to W²=3.5 GeV² and Q²=m_c²
Reference graph
Works this paper leans on
-
[1]
Nucleon structure functions For illustration, we consider the unpolarized inclusive scattering of a charged leptonℓ (with 4-momentumk) from a nucleonNof massM(4-momentumP) to a detected final state leptonℓ ′ (4-momentumk ′) and unobserved hadronsX,ℓN→ℓ ′X. For neutral currents the scattering can proceed through the exchange of a virtual photon or aZboson ...
-
[2]
(0)” and “(1)
Nuclear structure functions For inclusive DIS from nuclei, we follow previous JAM analyses [25] in using the nu- clear impulse (or weak nuclear binding) approximation to relate the structure function of a nucleusAto nucleon on-shell and off-shell structure functions via the convolution frame- work [21, 37–41], F A 2 (xB, Q2) = X N h f (on) N/A ⊗F N 2 +f (...
- [3]
-
[4]
Aaltonenet al., Science376, 170 (2022)
T. Aaltonenet al., Science376, 170 (2022)
2022
- [5]
- [6]
-
[7]
S. Tkachenkoet al., Phys. Rev. C89, 045206 (2014), [Addendum: Phys. Rev. C90, 059901 (2014)], arXiv:1402.2477 [nucl-ex]
-
[8]
Seelyet al., New measurements of the EMC effect in very light nuclei, Phys
J. Seelyet al., Phys. Rev. Lett.103, 202301 (2009), arXiv:0904.4448 [nucl-ex]
-
[9]
D. Biswaset al., Phys. Rev. Lett.135, 151902 (2025), arXiv:2409.15236 [hep-ex]
-
[10]
D. Abramset al., Phys. Rev. Lett.128, 132003 (2022), arXiv:2104.05850 [hep-ex]
-
[11]
D. Abramset al., Phys. Rev. Lett.135, 062502 (2025), arXiv:2410.12099 [nucl-ex]
-
[12]
A. Accardi, M. E. Christy, C. E. Keppel, W. Melnitchouk, P. Monaghan, J. G. Morf´ ın, and J. F. Owens, Phys. Rev. D81, 034016 (2010), arXiv:0911.2254 [hep-ph]
-
[13]
A. Accardi, L. T. Brady, W. Melnitchouk, J. F. Owens, and N. Sato, Phys. Rev. D93, 114017 (2016), arXiv:1602.03154 [hep-ph]
work page Pith review arXiv 2016
-
[14]
M. Cerutti, A. Accardi, I. P. Fernando, S. Li, J. F. Owens, and S. Park, Phys. Rev. D111, 094013 (2025), arXiv:2501.06849 [hep-ph]
- [15]
- [16]
- [17]
- [18]
- [19]
-
[20]
I. Schienbeinet al., J. Phys. G35, 053101 (2008), arXiv:0709.1775 [hep-ph]
-
[21]
Georgi and H
H. Georgi and H. D. Politzer, Phys. Rev. D14, 1829 (1976). 31
1976
-
[22]
C. Cocuzza, W. Melnitchouk, A. Metz, and N. Sato, Phys. Rev. D104, 074031 (2021), arXiv:2109.00677 [hep-ph]
-
[23]
C. Cocuzza, C. E. Keppel, H. Liu, W. Melnitchouk, A. Metz, N. Sato, and A. W. Thomas, Phys. Rev. Lett.127, 242001 (2021), arXiv:2104.06946 [hep-ph]
-
[24]
C. Cocuzza, W. Melnitchouk, A. Metz, and N. Sato, Phys. Rev. D106, L031502 (2022), arXiv:2202.03372 [hep-ph]
-
[25]
T. Anderson, W. Melnitchouk, and N. Sato, Phys. Rev. D112, 094011 (2025), arXiv:2501.00665 [hep-ph]
-
[26]
C. Cocuzza, N. T. Hunt-Smith, W. Melnitchouk, N. Sato, and A. W. Thomas, Phys. Rev. D 112, 114017 (2025), arXiv:2506.13616 [hep-ph]
-
[27]
C. Cocuzza, T. J. Hague, W. Melnitchouk, N. Sato, and A. W. Thomas, arXiv:2602.16589 [hep-ph]
-
[28]
P. A. Zylaet al.(Particle Data Group), PTEP2020, 083C01 (2020)
2020
-
[29]
O. W. Greenberg and D. Bhaumik, Phys. Rev. D4, 2048 (1971)
2048
-
[30]
Nachtmann, Nucl
O. Nachtmann, Nucl. Phys.B63, 237 (1973)
1973
-
[31]
J. C. Collins and D. E. Soper, Nucl. Phys.B194, 445 (1982)
1982
-
[32]
J. C. Collins, D. E. Soper, and G. F. Sterman, Adv. Ser. Direct. High Energy Phys.5, 1 (1989), arXiv:hep-ph/0409313
work page Pith review arXiv 1989
-
[33]
Collins,Foundations of perturbative QCD, Vol
J. Collins,Foundations of perturbative QCD, Vol. 32 (Cambridge University Press, 2013)
2013
-
[34]
E. G. Floratos, C. Kounnas, and R. Lacaze, Nucl. Phys.B192, 417 (1981)
1981
- [35]
- [36]
- [37]
-
[38]
R. L. Jaffe, Nucl. Phys.B229, 205 (1983)
1983
-
[39]
W. Melnitchouk, A. W. Schreiber, and A. W. Thomas, Phys. Rev. D49, 1183 (1994), arXiv:nucl-th/9311008
-
[40]
W. Melnitchouk, A. W. Schreiber, and A. W. Thomas, Phys. Lett. B335, 11 (1994), arXiv:nucl-th/9407007. 32
- [41]
- [42]
-
[43]
A. J. Tropiano, J. J. Ethier, W. Melnitchouk, and N. Sato, Phys. Rev. C99, 035201 (2019), arXiv:1811.07668 [nucl-th]
work page Pith review arXiv 2019
- [44]
- [45]
-
[46]
N. Sato, W. Melnitchouk, S. E. Kuhn, J. J. Ethier, and A. Accardi, Phys. Rev. D93, 074005 (2016), arXiv:1601.07782 [hep-ph]
work page Pith review arXiv 2016
- [47]
-
[48]
J. J. Ethier, N. Sato, and W. Melnitchouk, Phys. Rev. Lett.119, 132001 (2017), arXiv:1705.05889 [hep-ph]
work page Pith review arXiv 2017
- [49]
- [50]
-
[51]
A. C. Benvenutiet al., Phys. Lett. B223, 485 (1989)
1989
-
[52]
M. Arneodoet al., Nucl. Phys.B483, 3 (1997), arXiv:hep-ph/9610231
-
[53]
M. Arneodoet al., Nucl. Phys.B487, 3 (1997), arXiv:hep-ex/9611022
-
[54]
L. W. Whitlow, E. M. Riordan, S. Dasu, S. Rock, and A. Bodek, Phys. Lett. B282, 475 (1992)
1992
-
[55]
H. Abramowiczet al., Eur. Phys. J. C75, 580 (2015), arXiv:1506.06042 [hep-ex]
work page Pith review arXiv 2015
-
[56]
J. C. Webb,Measurement of continuum dimuon production in 800-GeV/C proton nucleon collisions, Ph.D. thesis, New Mexico State U. (2003), arXiv:hep-ex/0301031
work page internal anchor Pith review arXiv 2003
- [57]
-
[58]
Doveet al., Nature590, 561 (2021), arXiv:2103.04024 [hep-ph]
J. Doveet al., Nature590, 561 (2021), arXiv:2103.04024 [hep-ph]
- [59]
-
[60]
T. Aaltonenet al., Phys. Rev. Lett.102, 181801 (2009), arXiv:0901.2169 [hep-ex]
- [61]
- [62]
-
[63]
F. Ringer and W. Vogelsang, Phys. Rev. D91, 094033 (2015), arXiv:1503.07052 [hep-ph]
-
[64]
Chatrchyanet al., JHEP04, 050 (2011), arXiv:1103.3470 [hep-ex]
S. Chatrchyanet al., JHEP04, 050 (2011), arXiv:1103.3470 [hep-ex]
-
[65]
S. Chatrchyanet al., Phys. Rev. Lett.109, 111806 (2012), arXiv:1206.2598 [hep-ex]
-
[66]
S. Chatrchyanet al., Phys. Rev. D90, 032004 (2014), arXiv:1312.6283 [hep-ex]
-
[67]
V. Khachatryanet al., Eur. Phys. J. C76, 469 (2016), arXiv:1603.01803 [hep-ex]
-
[68]
Aaijet al., JHEP12, 079 (2014), arXiv:1408.4354 [hep-ex]
R. Aaijet al., JHEP12, 079 (2014), arXiv:1408.4354 [hep-ex]
-
[69]
Aaijet al., JHEP05, 131 (2016), arXiv:1605.00951 [hep-ex]
R. Aaijet al., JHEP05, 131 (2016), arXiv:1605.00951 [hep-ex]
-
[70]
J. Adamet al., Phys. Rev. D103, 012001 (2021), arXiv:2011.04708 [nucl-ex]
- [71]
-
[72]
A. Abulenciaet al., Phys. Rev. D75, 092006 (2007), [Erratum: Phys. Rev. D75, 119901 (2007)], arXiv:hep-ex/0701051
- [73]
-
[74]
Aadet al., JHEP05, 068 (2014), arXiv:1402.6263 [hep-ex]
G. Aadet al., JHEP05, 068 (2014), arXiv:1402.6263 [hep-ex]
-
[75]
Chatrchyanet al., JHEP02, 013 (2014), arXiv:1310.1138 [hep-ex]
S. Chatrchyanet al., JHEP02, 013 (2014), arXiv:1310.1138 [hep-ex]
- [76]
-
[77]
Lacombe, B
M. Lacombe, B. Loiseau, R. Vinh Mau, J. Cote, P. Pires, and R. de Tourreil, Phys. Lett. B 101, 139 (1981)
1981
-
[78]
A. Kievsky, E. Pace, G. Salme, and M. Viviani, Phys. Rev. C56, 64 (1997), arXiv:nucl- th/9704050
- [79]
-
[80]
R. Machleidt, Phys. Rev. C63, 024001 (2001), arXiv:nucl-th/0006014
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.