REVIEW 4 major objections 5 minor 3 cited by
Deep-inelastic scattering off helium-3 and tritium cannot be described by free-nucleon quark distributions alone; nucleon off-shell corrections with both isoscalar and isovector components are required.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-02 22:27 UTC pith:U7YRDQCU
load-bearing objection Core result survives scrutiny: A<=3 DIS demands off-shell corrections and MARATHON's KP normalization looks biased; the isovector signal is suggestive but model-dependent. the 4 major comments →
Isospin dependence of nuclear EMC effect from global QCD analysis
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that nucleon off-shell corrections—the modification of the bound nucleon's quark distributions relative to its free-nucleon counterpart—are required to describe inclusive deep-inelastic scattering from A=2 and A=3 nuclei. In the global fit, the 3He/D ratio data rule out the on-shell-only description (chi2_red = 5.10) and are well described once off-shell corrections are introduced (chi2_red = 1.04). The analysis further suggests a nonzero isovector component, visible in the difference of the off-shell u and d distributions at large x, alongside a large isoscalar component. Consequently, the extracted EMC ratios R_D, R_3He, and R_3H at x_B = 0.31 disagree with the values
What carries the argument
The central object is the nuclear structure function expressed as a convolution of on-shell and off-shell nucleon contributions: F2^A = sum_N [ f^on_{N/A} ⊗ F2^N + f^off_{N/A} ⊗ δF2^{N/A} ]. The unknown off-shell correction δF2 is parametrized through off-shell PDFs δq_{N/A}, which are decomposed into an isoscalar part δq0 and isovector parts δu1, δd1 with a prescribed combinatorics (the isovector part vanishes in 3He and doubles in 3H relative to the deuteron). The machinery is a simultaneous Bayesian fit that determines free-nucleon PDFs, off-shell functions, higher-twist corrections, and data normalizations from the same global dataset, so that nuclear effects are inferred rather than ass
Load-bearing premise
The load-bearing premise is the assumed isoscalar/isovector decomposition of the off-shell PDFs: that the isovector part vanishes in 3He, doubles in 3H relative to the deuteron, and that only three fitted functions suffice; if this counting is wrong, the extracted isovector signal and EMC ratios are artifacts of the model.
What would settle it
A precision measurement of the 3He/D and 3H/D ratios with an independent normalization—for example, through tagged deep-inelastic scattering from the spectator nucleon—that disagrees with the fitted off-shell predictions, or a direct measurement of F_L^A showing a nuclear dependence of R = F_L/F_T large enough to shift the structure-function ratios, would falsify the extraction.
If this is right
- The MARATHON data normalization based on the heavy-nucleus model biases the extracted neutron structure function; re-extracted ratios shift accordingly.
- The nuclear EMC effect is shown to persist down to A=2 and A=3, with magnitude tracking the off-shell smearing averages: about −4% for the deuteron, −7% for helium-3, and −9.5% for tritium.
- The d-quark PDF at x ≳ 0.6 gains modest new constraints from the A=3 data, reducing its uncertainties in that region.
- The super-ratio R_{3He/3H} is predicted to fall below unity at large x, contrary to the trend in the normalization model used by the experiment.
- Future A=3 data can discriminate between isoscalar-only and isovector-bearing off-shell scenarios, since the two produce different high-x behavior in the 3H/D ratio.
Where Pith is reading between the lines
- If the isovector off-shell signal holds up, it would imply the medium modification of quark distributions depends on isospin, which is a direct discriminator between mean-field and short-range-correlation mechanisms of the EMC effect—a connection the paper leaves implicit.
- The same isoscalar/isovector decomposition could be applied to 4He, where the balance between the two components differs, offering a testable prediction using existing EMC data on helium-4.
- Because the analysis treats the measured cross-section ratios as structure-function ratios, a future direct measurement of F_L for A=3 nuclei could shift the extracted off-shell functions; the size of that shift is not quantified in the paper.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a global QCD analysis of unpolarized PDFs and nucleon off-shell corrections using proton, deuteron, and A=3 data, including the new MARATHON ^3He/D and ^3H/D cross-section ratios. The authors parametrize off-shell effects through isoscalar and isovector components under a specific spectator-isospin decomposition, fit them simultaneously with PDFs in the JAM Bayesian Monte Carlo framework, and find that off-shell corrections are required: the ^3He/D reduced chi^2 drops from 5.10 to 1.04 when they are included. They further report a 'suggestion' of isovector off-shell effects from a fitted δu_1−δd_1 that is negative at x≳0.5, and they find that the extracted EMC ratios R_D, R_3He, and R_3H deviate from KP-model values at x_B=0.31 by about 3.0σ, 2.3σ, and 2.4σ, respectively. They also argue that the KP-based normalizations applied to MARATHON data introduce a significant bias, with the fitted ^3He/D normalization differing from the KP value by 3.2σ.
Significance. If correct, this would be an important step: it provides a global, data-driven determination of light-nucleus off-shell corrections, tests the MARATHON normalization procedure, and gives quantitative EMC ratios for A=2 and A=3 that can be compared with few-body calculations. The analysis is strengthened by its global data set, Bayesian uncertainty propagation, explicit treatment of normalizations, and checks on higher-twist parametrizations. The need for off-shell corrections to describe ^3He/D appears robust from the goodness-of-fit comparison. However, the isovector signal and the quantitative EMC-ratio deviations from KP are conditional on the assumed isoscalar/isovector decomposition and on a model-selection choice that is not fully justified. Those load-bearing points need to be addressed before the central claims can be regarded as established.
major comments (4)
- [§III and §II.B] The paper chooses scenario (ii) as default even though scenario (i) has a better total MARATHON chi^2 (50 vs 54), with the only stated reason that scenario (ii) better describes the highest three x_B points of ^3H/D. Since the isovector signal x(δu_1−δd_1)<0 is literally the fitted difference of scenario-(ii) functions, this is a post-hoc model-selection rule. Please provide a formal model-selection or cross-validation justification. More importantly, test a more flexible isoscalar-only model (e.g. more parameters for δq_0) and show whether the high-x ^3H/D points can be described without isovector terms. Without such a test, the 'suggestion of isovector off-shell effects' is not a robust data requirement.
- [§II.B, Eqs. (14)–(18)] The isoscalar/isovector decomposition relies on factor-2 combinatorics, the identification of spectator-isospin projections, and the assumption that isovector terms vanish in ^3He. These assumptions are physically motivated but not tested. With only three independent MARATHON ratios constraining three off-shell functions, the extracted signals are strongly shaped by this basis. Please demonstrate sensitivity to these assumptions: allow non-zero isovector contributions in ^3He, vary the combinatorial factors, or use an alternative functional basis for the off-shell corrections, and quantify the shifts in R_D, R_3He, R_3H and δu_1−δd_1. This is also relevant to the abstract's claim of 'without theoretical assumptions about the isospin dependence of nuclear effects,' which is not literally satisfied.
- [§III, paragraph on MARATHON cross-section ratios] The paper equates MARATHON cross-section ratios with F_2 structure-function ratios, neglecting the nuclear dependence of R=σ_L/σ_T. The authors acknowledge that this is an assumption but do not estimate its impact. Since the central EMC-ratio claims are at the few-percent level, and high-x data are most sensitive, please provide a quantitative estimate of the resulting systematic uncertainty (e.g. using model estimates of nuclear R or a conservative uncertainty band) or restrict the conclusion to a kinematic region where this effect is demonstrably negligible.
- [Table II and Fig. 2] The 3.2σ discrepancy between the fitted ^3He/D normalization and the KP-based normalization is presented as evidence that KP normalizations bias the analysis. However, the normalization and the off-shell shape parameters are fitted simultaneously; a more flexible off-shell model could trade a normalization shift for a shape change. Please show that this discrepancy persists under the alternative off-shell parametrizations suggested above, or characterize the correlation between the normalization parameters and the off-shell functional form.
minor comments (5)
- [Fig. 5 caption] The symbol '□' in x(δu_1 □ δd_1) should be a minus sign.
- [Conclusions] The symbols R_ht, R_hd, and R_td appear without definition; presumably R_D, R_3He/D, and R_3H/D. Please use consistent notation.
- [Fig. 9 caption] The lower panel caption says '3He/D and 3He/D'; the second should likely be '3H/D'.
- [Various] Typos: 'constributions' in §II.B; 'T otal' in Table I; 'the shed light' in the Conclusions; 'off-shell isovector δu1−δd1' wording could be clarified.
- [Abstract and §II.B] The phrase 'without theoretical assumptions about the isospin dependence' overstates the situation; the analysis employs a specific isoscalar/isovector decomposition and scenario restrictions. Rephrasing to 'without assuming a specific magnitude/shape of the isospin dependence' would be more accurate.
Circularity Check
No significant circularity: off-shell corrections are fit outputs compared against an external benchmark; the isovector signal is a fitted parameter, not a construction-forced prediction.
full rationale
The closest step to circularity is the extraction of the isovector off-shell signal: δu1−δd1 is a fitted function (Sec. IV.B, Fig. 5) obtained from the same MARATHON data used to argue for it. But this is a normal likelihood-fit output, not a prediction; the data could in principle have returned δu1=δd1=0, and the nonzero result is an inference, not an identity. The preference for scenario (ii) over (i) is explicitly data-driven (Sec. III: 'scenario (ii) was found to better describe the highest 3 data points in x_B for the 3H/D data, suggesting a potential signature of the isovector EMC effect'); this is model selection based on a few data points, which is fragile but not circular. The EMC-ratio deviations from the KP model are comparisons of fit outputs to an external benchmark; the KP model's unity-at-x=0.31 condition is not imposed in the fit, so the disagreement is not by construction. The self-citation to the earlier JAM analysis [32] is confirmatory rather than load-bearing: the current fit includes new 3He/D and 3H/D data and re-derives the result. The acknowledged model dependence (Sec. IV: 'the central values of the nuclear EMC ratios, as well as their uncertainties, are significantly influenced by the modeling of the off-shell functions') is a limitation on robustness, not circularity. No equation in the paper reduces to its input by construction, and no fitted parameter is relabeled as a prediction.
Axiom & Free-Parameter Ledger
free parameters (5)
- Onshell PDF shape parameters (33) =
not quoted in text
- Off-shell function parameters (9) =
not quoted in text
- Higher-twist parameters (10) =
not quoted in text
- Dataset normalizations (22) =
e.g., MARATHON D/p 1.017(4), 3He/D 0.996(6), 3H/D 0.991(5)
- W^2 cut =
3.5 GeV^2
axioms (6)
- domain assumption Impulse approximation and convolution formula Eqs. (1)-(2) for nuclear structure functions.
- domain assumption Exact charge symmetry relations (7)-(8) between mirror nuclei, with no charge-symmetry violation or Coulomb corrections.
- ad hoc to paper Isoscalar/isovector decomposition (14)-(18) with factor-2 counting and zero isovector component in 3He.
- domain assumption MARATHON cross-section ratios equal structure-function ratios.
- domain assumption Off-shell sea-quark and gluon contributions are neglected.
- ad hoc to paper PDF template Eq. (23) and default scenario (ii): symmetric isoscalar and asymmetric isovector off-shell corrections.
read the original abstract
We perform a new global QCD analysis of unpolarized parton distribution functions (PDFs) in the nucleon from proton, deuteron and $A=3$ data, including recent measurements of $^3$He/$D$ and $^3$H/$D$ cross section ratios from the MARATHON experiment at Jefferson Lab. Simultaneously inferring the PDFs and nucleon off-shell corrections allows both to be determined consistently, without theoretical assumptions about the isospin dependence of nuclear effects. The analysis provides strong evidence for the need of nucleon off-shell corrections to describe the $A=3$ data, with large isoscalar and a suggestion of nonzero isovector contributions in $A \leq 3$ nuclei. We find that the extracted EMC ratios of nuclear to nucleon structure functions for $A=2$ and 3 differ from those naively extrapolated from heavy nuclei down to low $A$.
Figures
Forward citations
Cited by 3 Pith papers
-
Mapping data sensitivities in global QCD analysis with linear response and influence functions
A framework based on linear response and influence functions maps data sensitivities in global QCD analyses to show how experiments determine central values, uncertainties, and correlations of non-perturbative functions.
-
Impact of parity-violating deep-inelastic scattering on the weak mixing angle and high-$x$ parton distributions
Simulations of 11 and 22 GeV PVDIS pseudodata indicate potential constraints on low-Q² sin²θ_W and high-x strange quark and d/u PDFs, with strong correlations requiring simultaneous QCD-electroweak analysis.
-
Stability of parton distributions at high $x$: impact of nuclear and power corrections
Global fit finds u and d PDFs stable to x≈0.8 with positive isospin-independent higher-twist corrections and nonzero off-shell nucleon contributions required to describe nuclear data.
Reference graph
Works this paper leans on
-
[1]
isoscalar
= N M xα(1−x) β(1 +γ √x+ηx),(23) wherea={N, α, β, γ, η}is the set of parameters to be inferred. The normalization constant M= B[α+ 2, β+ 1] +γB[α+ 5 2 , β+ 1] +ηB[α+ 3, β+ 1], where B is the Euler beta function, normalizes the function to the second moment. This form allows maximum decorrelation between the fitted normalizations and shape parameters. Our ...
-
[2]
J. J. Aubertet al., Phys. Lett. B123, 275 (1983)
1983
-
[3]
Gomezet al., Phys
J. Gomezet al., Phys. Rev. D49, 4348 (1994)
1994
-
[4]
A. W. Thomas, A. Michels, A. W. Schreiber, and P. A. M. Guichon, Phys. Lett. B233, 43 (1989)
1989
-
[5]
Saito, A
K. Saito, A. Michels, and A. W. Thomas, Phys. Rev. C46, R2149 (1992)
1992
-
[6]
H. Mineo, W. Bentz, N. Ishii, A. W. Thomas, and K. Yazaki, Nucl. Phys.A735, 482 (2004), arXiv:nucl-th/0312097
Pith/arXiv arXiv 2004
-
[7]
I. C. Cloet, W. Bentz, and A. W. Thomas, Phys. Lett. B642, 210 (2006), arXiv:nucl- th/0605061
arXiv 2006
-
[8]
J. R. Smith and G. A. Miller, Phys. Rev. C65, 055206 (2002), arXiv:nucl-th/0202016
Pith/arXiv arXiv 2002
-
[9]
P. A. M. Guichon, J. R. Stone, and A. W. Thomas, Prog. Part. Nucl. Phys.100, 262 (2018), arXiv:1802.08368 [nucl-th]
Pith/arXiv arXiv 2018
-
[10]
Schmookleret al., Nature566, 354 (2019), arXiv:2004.12065 [nucl-ex]
B. Schmookleret al., Nature566, 354 (2019), arXiv:2004.12065 [nucl-ex]
Pith/arXiv arXiv 2019
-
[11]
X. G. Wang, A. W. Thomas, and W. Melnitchouk, Phys. Rev. Lett.125, 262002 (2020), arXiv:2004.03789 [hep-ph]
Pith/arXiv arXiv 2020
-
[12]
E. P. Segarra, J. R. Pybus, F. Hauenstein, D. W. Higinbotham, G. A. Miller, E. Piasetzky, A. Schmidt, M. Strikman, L. B. Weinstein, and O. Hen, Phys. Rev. Res.3, 023240 (2021), arXiv:2006.10249 [hep-ph]
Pith/arXiv arXiv 2021
-
[13]
W. Xing, X. G. Wang, and A. W. Thomas, Phys. Lett. B846, 138195 (2023), 29 arXiv:2305.13666 [hep-ph]
Pith/arXiv arXiv 2023
-
[14]
D. N. Kim, O. Hen, G. A. Miller, E. Piasetzky, M. Strikman, and L. Weinstein, Phys. Rev. C111, 065201 (2025), arXiv:2404.15442 [nucl-th]
Pith/arXiv arXiv 2025
-
[15]
D. F. Geesaman, K. Saito, and A. W. Thomas, Ann. Rev. Nucl. Part. Sci.45, 337 (1995)
1995
-
[16]
P. R. Norton, Rept. Prog. Phys.66, 1253 (2003)
2003
-
[17]
A. W. Thomas, Int. J. Mod. Phys. E27, 1840001 (2019), arXiv:1809.06622 [hep-ph]
Pith/arXiv arXiv 2019
-
[18]
A. J. Tropiano, J. J. Ethier, W. Melnitchouk, and N. Sato, Phys. Rev. C99, 035201 (2019), arXiv:1811.07668 [nucl-th]
Pith/arXiv arXiv 2019
-
[19]
G. V. Dunne and A. W. Thomas, Nucl. Phys.A455, 701 (1986)
1986
-
[20]
S. V. Akulinichev, S. A. Kulagin, and G. M. Vagradov, Phys. Lett. B158, 485 (1985)
1985
-
[21]
R. P. Bickerstaff and A. W. Thomas, J. Phys. G15, 1523 (1989)
1989
-
[22]
W. Melnitchouk, A. W. Schreiber, and A. W. Thomas, Phys. Rev. D49, 1183 (1994), arXiv:nucl-th/9311008
Pith/arXiv arXiv 1994
-
[23]
W. Melnitchouk, A. W. Schreiber, and A. W. Thomas, Phys. Lett. B335, 11 (1994), arXiv:nucl-th/9407007
Pith/arXiv arXiv 1994
-
[24]
S. A. Kulagin, W. Melnitchouk, G. Piller, and W. Weise, Phys. Rev. C52, 932 (1995), arXiv:hep-ph/9504377
Pith/arXiv arXiv 1995
-
[25]
S. A. Kulagin, G. Piller, and W. Weise, Phys. Rev. C50, 1154 (1994), arXiv:nucl-th/9402015
Pith/arXiv arXiv 1994
-
[26]
S. A. Kulagin and R. Petti, Nucl. Phys.A765, 126 (2006), arXiv:hep-ph/0412425
Pith/arXiv arXiv 2006
-
[27]
I. R. Afnan, F. R. P. Bissey, J. Gomez, A. T. Katramatou, W. Melnitchouk, G. G. Petratos, and A. W. Thomas, Phys. Lett. B493, 36 (2000), arXiv:nucl-th/0006003
Pith/arXiv arXiv 2000
-
[28]
I. R. Afnan, F. R. P. Bissey, J. Gomez, A. T. Katramatou, S. Liuti, W. Melnitchouk, G. G. Petratos, and A. W. Thomas, Phys. Rev. C68, 035201 (2003), arXiv:nucl-th/0306054
Pith/arXiv arXiv 2003
-
[29]
E. Pace, G. Salme, S. Scopetta, and A. Kievsky, Phys. Rev. C64, 055203 (2001), arXiv:nucl- th/0109005
arXiv 2001
-
[30]
M. M. Sargsian, S. Simula, and M. I. Strikman, Phys. Rev. C66, 024001 (2002), arXiv:nucl- th/0105052
arXiv 2002
-
[31]
D. Abramset al., Phys. Rev. Lett.128, 132003 (2022), arXiv:2104.05850 [hep-ex]
Pith/arXiv arXiv 2022
-
[32]
D. Abramset al., Phys. Rev. Lett.135, 062502 (2025), arXiv:2410.12099 [nucl-ex]
arXiv 2025
-
[33]
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]. 30
Pith/arXiv arXiv 2021
-
[34]
S. A. Kulagin and R. Petti, Phys. Rev. C82, 054614 (2010), arXiv:1004.3062 [hep-ph]
Pith/arXiv arXiv 2010
-
[35]
L. T. Brady, A. Accardi, T. J. Hobbs, and W. Melnitchouk, Phys. Rev. D84, 074008 (2011), [Erratum: Phys. Rev. D85, 039902 (2012)], arXiv:1108.4734 [hep-ph]
Pith/arXiv arXiv 2011
-
[36]
I. Schienbeinet al., J. Phys. G35, 053101 (2008), arXiv:0709.1775 [hep-ph]
Pith/arXiv arXiv 2008
-
[37]
M. A. G. Aivazis, J. C. Collins, F. I. Olness, and W.-K. Tung, Phys. Rev. D50, 3102 (1994), arXiv:hep-ph/9312319
Pith/arXiv arXiv 1994
-
[38]
E. Moffat, T. C. Rogers, W. Melnitchouk, N. Sato, and F. Steffens, Phys. Rev. D99, 096008 (2019), arXiv:1901.09016 [hep-ph]
Pith/arXiv arXiv 2019
-
[39]
R. L. Jaffe, Nucl. Phys.B229, 205 (1983)
1983
-
[40]
F. E. Close and A. W. Thomas, Phys. Lett. B212, 227 (1988)
1988
-
[41]
N. Sato, W. Melnitchouk, S. E. Kuhn, J. J. Ethier, and A. Accardi, Phys. Rev. D93, 074005 (2016), arXiv:1601.07782 [hep-ph]
Pith/arXiv arXiv 2016
-
[42]
N. Sato, J. J. Ethier, W. Melnitchouk, M. Hirai, S. Kumano, and A. Accardi, Phys. Rev. D 94, 114004 (2016), arXiv:1609.00899 [hep-ph]
Pith/arXiv arXiv 2016
-
[43]
J. J. Ethier, N. Sato, and W. Melnitchouk, Phys. Rev. Lett.119, 132001 (2017), arXiv:1705.05889 [hep-ph]
Pith/arXiv arXiv 2017
-
[44]
N. Sato, C. Andres, J. J. Ethier, and W. Melnitchouk, Phys. Rev. D101, 074020 (2020), arXiv:1905.03788 [hep-ph]
Pith/arXiv arXiv 2020
-
[45]
E. Moffat, W. Melnitchouk, T. C. Rogers, and N. Sato, Phys. Rev. D104, 016015 (2021), arXiv:2101.04664 [hep-ph]
Pith/arXiv arXiv 2021
-
[46]
C. Cocuzza, W. Melnitchouk, A. Metz, and N. Sato, Phys. Rev. D104, 074031 (2021), arXiv:2109.00677 [hep-ph]
Pith/arXiv arXiv 2021
-
[47]
T. Anderson, W. Melnitchouk, and N. Sato, Phys. Rev. D112, 094011 (2025), arXiv:2501.00665 [hep-ph]
arXiv 2025
-
[48]
D. Biswaset al., Phys. Rev. Lett.135, 151902 (2025), arXiv:2409.15236 [hep-ex]
Pith/arXiv arXiv 2025
-
[49]
J. Seelyet al., Phys. Rev. Lett.103, 202301 (2009), arXiv:0904.4448 [nucl-ex]
Pith/arXiv arXiv 2009
-
[50]
Liet al., Nature609, 41 (2022), arXiv:2210.04189 [nucl-ex]
S. Liet al., Nature609, 41 (2022), arXiv:2210.04189 [nucl-ex]
Pith/arXiv arXiv 2022
-
[51]
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]
Pith/arXiv arXiv 2025
-
[52]
S. I. Alekhin, S. A. Kulagin, and R. Petti, Phys. Rev. D96, 054005 (2017), arXiv:1704.00204 [nucl-th]. 31
Pith/arXiv arXiv 2017
-
[53]
Lacombe, B
M. Lacombe, B. Loiseau, R. Vinh Mau, J. Cote, P. Pires, and R. de Tourreil, Phys. Lett. B 101, 139 (1981)
1981
-
[54]
A. Kievsky, E. Pace, G. Salme, and M. Viviani, Phys. Rev. C56, 64 (1997), arXiv:nucl- th/9704050
arXiv 1997
-
[55]
S. I. Alekhin, S. A. Kulagin, and R. Petti, PoSDIS2024, 067 (2025), arXiv:2410.21472 [hep-ph]
Pith/arXiv arXiv 2025
-
[56]
I. Albayraket al., Nucl. Instrum. Meth. A1062, 169190 (2024), arXiv:2402.01904 [physics.ins- det]. 32
Pith/arXiv arXiv 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.