Redetermination of proton sea distributions
Pith reviewed 2026-06-30 18:13 UTC · model grok-4.3
The pith
NNLO fits to HERA and ATLAS data reveal an excess of anti-up over anti-down quarks in the proton sea for x between 0.01 and 1.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The shapes of light flavor sea quark distributions of the proton are examined directly from two rounds of NNLO global analysis of HERA e±p deep inelastic scattering cross section measurements, termed as HERAshape, and the ATLAS measurement of W±/Z production from pp collision at √s = 7 TeV, termed as ATLASshape. An asymmetric distribution between anti-up (ū) quark and anti-down (d̄) quark is found in both analysis, showing that the anti-up quark distribution exceeds over anti-down quark distribution in the momentum fraction range x∈(10^{-2}, 1) of these partons. The Gottfried Sum Rule is reevaluated from these extracted parton distribution functions and the obtained value differs surprisingl
What carries the argument
NNLO global analysis procedure applied to HERA deep inelastic scattering and ATLAS electroweak boson production data for extracting sea quark parton distributions.
If this is right
- Light sea quarks exhibit ū(x) > d̄(x) asymmetry in the specified x range.
- The Gottfried integral from these PDFs deviates from previous NMC and NuSea results.
- Global parton distribution functions should incorporate this asymmetry in future fits.
- Predictions for processes sensitive to light-flavor sea quarks may change at moderate x.
Where Pith is reading between the lines
- The asymmetry could alter cross-section predictions for other high-energy processes involving sea quarks.
- Independent data from higher-energy colliders could test whether the imbalance persists or evolves with scale.
- Theoretical models of sea-quark generation may need explicit flavor-breaking mechanisms to match this extraction.
Load-bearing premise
The NNLO analysis with the selected HERA and ATLAS data accurately determines the light sea quark distributions without interference from unmodeled effects or data biases.
What would settle it
A new global PDF fit or direct measurement of sea quark distributions that yields ū(x) ≤ d̄(x) across the same x range and restores the prior Gottfried sum rule value would falsify the asymmetry claim.
Figures
read the original abstract
The shapes of light flavor sea quark distributions of the proton are examined directly from two rounds of NNLO global analysis of HERA $e^{\pm}p$ deep inelastic scattering cross section measurements, termed as HERAshape and the ATLAS measurement of $W^{\pm}/Z$ production from $pp$ collision at $\sqrt{s}$ = 7 TeV, termed as ATLASshape. An asymmetric distribution between anti-up ($\ubar$) quark and anti-down ($\dbar$) quark is found in both analysis, showing that the anti-up quark distribution exceeds over anti-down quark distribution in the momentum fraction range $x\in(10^{-2}, 1)$ of these partons. The Gottfried Sum Rule is reevaluated from these extracted parton distribution functions and the obtained value differs surprisingly from that of the NMC and the NuSea Collaborations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript performs two separate NNLO global PDF analyses—one using HERA e±p DIS cross sections (HERAshape) and one using ATLAS W±/Z production data at 7 TeV (ATLASshape)—and reports that both yield ū(x) > d̄(x) for x ∈ (10^{-2}, 1). From the resulting PDFs the Gottfried sum rule is reevaluated and found to differ from the values obtained by the NMC and NuSea collaborations.
Significance. If the asymmetry is shown to be independently constrained by each dataset rather than imposed by parametrization choices, the result would provide a useful cross-check on light-sea distributions and motivate a re-examination of the numerical value of the Gottfried sum rule.
major comments (3)
- [HERAshape analysis description] HERAshape analysis: inclusive HERA e±p structure functions constrain only the charge-weighted sum (4/9)(u + ū) + (1/9)(d + d̄) + …; ū and d̄ enter with identical coefficients, supplying no direct lever arm on ū − d̄. The reported asymmetry therefore requires explicit demonstration (e.g., Δχ² when the asymmetry parameter is fixed to zero) that it is data-driven rather than an artifact of the chosen functional form.
- [Abstract and Section 4] Abstract and results section: the two analyses are presented symmetrically as mutual confirmation, yet the HERAshape result lacks an independent constraint on the asymmetry. The joint claim and the subsequent Gottfried-sum reevaluation therefore rest primarily on the ATLAS W±/Z data; this asymmetry of constraining power should be stated clearly.
- [Gottfried sum rule subsection] Gottfried-sum reevaluation: the numerical value and its uncertainty must be compared directly with the NMC and NuSea measurements, including a quantitative assessment of tension; the present text only states that the value “differs surprisingly.”
minor comments (2)
- [Abstract] The abstract omits any mention of the parametrization form, number of free parameters, data-selection cuts, or error treatment; a one-sentence summary of these elements would improve readability.
- [Introduction] Notation for the sea distributions (ū, d̄) should be defined at first use and kept consistent with standard PDF literature.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate the revisions that will be incorporated.
read point-by-point responses
-
Referee: [HERAshape analysis description] HERAshape analysis: inclusive HERA e±p structure functions constrain only the charge-weighted sum (4/9)(u + ū) + (1/9)(d + d̄) + …; ū and d̄ enter with identical coefficients, supplying no direct lever arm on ū − d̄. The reported asymmetry therefore requires explicit demonstration (e.g., Δχ² when the asymmetry parameter is fixed to zero) that it is data-driven rather than an artifact of the chosen functional form.
Authors: We agree that inclusive HERA e±p data constrain only the charge-weighted combination and therefore supply no direct experimental lever arm on ū − d̄. In the HERAshape analysis the observed asymmetry is generated by the flexible input-scale parametrization and is selected by the global χ² minimization. To demonstrate that the data prefer a non-zero asymmetry, we will add the requested Δχ² comparison (free fit versus asymmetry parameters fixed to zero) in the revised manuscript. revision: yes
-
Referee: [Abstract and Section 4] Abstract and results section: the two analyses are presented symmetrically as mutual confirmation, yet the HERAshape result lacks an independent constraint on the asymmetry. The joint claim and the subsequent Gottfried-sum reevaluation therefore rest primarily on the ATLAS W±/Z data; this asymmetry of constraining power should be stated clearly.
Authors: The referee is correct that the ATLAS W±/Z measurements supply the direct experimental constraint on the sea asymmetry, while the HERAshape result is a consistency check that relies on parametrization freedom. We will revise the abstract and Section 4 to state this difference in constraining power explicitly and to present the two analyses accordingly rather than as fully symmetric confirmations. revision: yes
-
Referee: [Gottfried sum rule subsection] Gottfried-sum reevaluation: the numerical value and its uncertainty must be compared directly with the NMC and NuSea measurements, including a quantitative assessment of tension; the present text only states that the value “differs surprisingly.”
Authors: We will expand the Gottfried-sum subsection to quote the central values and uncertainties from both the NMC and NuSea measurements and to include a quantitative assessment of tension (e.g., pull in units of standard deviation) between our result and each experimental value. revision: yes
Circularity Check
No significant circularity detected; derivation relies on external data fits.
full rationale
The paper performs standard NNLO global PDF fits to independent external datasets (HERA e±p DIS cross sections and ATLAS W±/Z production). The extracted ū > d̄ asymmetry and subsequent Gottfried sum reevaluation are direct outputs of these fits rather than reductions by construction to the paper's own inputs or self-citations. No equations, parametrizations, or claims in the provided text exhibit self-definitional equivalence, fitted inputs renamed as predictions, or load-bearing self-citation chains. The procedure is externally benchmarked against experimental data, satisfying the criteria for a self-contained, non-circular analysis.
Axiom & Free-Parameter Ledger
free parameters (1)
- PDF parametrization parameters
axioms (1)
- domain assumption NNLO perturbative QCD is adequate to describe the HERA and ATLAS data for sea-quark extraction
Forward citations
Cited by 1 Pith paper
-
Probing SU(2) Quark Flavor Asymmetry with W Bosons at RHIC
The STAR W+/W- cross section ratio is minimally affected by fiducial constraints and pT resummation, providing robust discrimination between ar d(x) and ar u(x) PDF models at x ~ 0.1 via L2 sensitivity and reweighting.
Reference graph
Works this paper leans on
-
[1]
(2) If the sea SU(2) symmetry is further assumed, ¯ u(x) = ¯d(x), we obtain the Quark Parton Model (QPM) prediction for GSR, to be 1 /3
reduces to S G = 1 3 + 2 3 ∫ 1 0 [ ¯u(x) −¯d(x) ] dx. (2) If the sea SU(2) symmetry is further assumed, ¯ u(x) = ¯d(x), we obtain the Quark Parton Model (QPM) prediction for GSR, to be 1 /3. The NMC result, 0 . 240 ±0. 016 [ 8], is roughly 10% lower than this value. This discrepancy be- tween the NMC result and QPM prediction of the GSR was considered as ...
-
[2]
Redetermination of proton sea distributions
reads as S G = 1 3 + 10 9 ∫ 1 0 [¯qp(x) −¯qn(x)] dx. (3) arXiv:2605.19918v1 [hep-ph] 19 May 2026 2 Thus, we can interpret the discrepancy between the NMC and QPM predictions, according to Eq. ( 3), as an indication that there are more sea quarks in the neutron than in the pro- ton. The contribution from isospin symmetry breaking to the GSR is expected to ...
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[3]
40 for charged current (CC) interaction (mediated by W ±), with electron beam energy of Ee ≃27
3 ×10−2 ≤x ≤0. 40 for charged current (CC) interaction (mediated by W ±), with electron beam energy of Ee ≃27. 5 GeV and proton beam energies of E p = 920, 820, 575, 460 GeV . The obtained PDFs from these data are free of any hadronic corrections and the assumption of isospin symmetry in the nucleon. Another important feature of the HERA data is that the ...
-
[4]
As for the strange quark factor, fs, it is reported that the ac- tual distribution of the strange quark is not suppressed in the proton [ 20, 23]
with freeing the parameters A, B of x ¯D and x ¯U. As for the strange quark factor, fs, it is reported that the ac- tual distribution of the strange quark is not suppressed in the proton [ 20, 23]. The ratio of ¯ s/ ¯d is measured to be
-
[5]
19+0. 15 −0. 16 with great precision by A TLAS collaboration, us- ing the HERA e±p data jointly with W ±/ Z production cross section data [ 20]. By Eq. ( 5), we have fs = 0.54 in our anal- ysis. By relaxing the restricted parameters in Eq. ( 4), the num- ber of free fit parameters then adds up to 16 from the 14 fit parameters of HERAPDF2.0. The fitting proce...
-
[6]
With the inclusion of the A TLAS W ±/ Z data and the application of the parame- terization form in Eq
The x¯u, x ¯d, x ¯s distributions of the A T- LASshape PDFs are well constrained with smaller uncertain- ties relative to the HERAshape PDFs. With the inclusion of the A TLAS W ±/ Z data and the application of the parame- terization form in Eq. ( 6), x¯u is slightly diminished in A T- LASshape PDFs with respect to both the HERAPDF2.0 and the HERAshape PDF...
-
[7]
The pale green bands in both figure represents the NMC results with 0. 227 ±0. 0016 and 0. 240 ±0. 016 corresponding in two x ranges. The S G obtained from NNPDF3.1, CT18NNLO crosses with that of the NMC within their large uncertainties, with the central values being larger than NMC result in 0 . 004 ≤ 4 4−10 3−10 2−10 1−10 1 x 0 0.2 0.4 0.6 0.8 1 )2(x,Q u...
-
[8]
E. D. Bloom, D. H. Coward, H. C. DeStaebler, J. Drees, G. Miller, L. W . Mo, R. E. Taylor, M. Breidenbach, J. I. Friedman, G. C. Hartmann and H. W . Kendall, Phys. Rev. Lett. 23, 930-934 (1969)
1969
-
[9]
M. Breidenbach, J. I. Friedman, H. W . Kendall, E. D. Bloom, D. H. Coward, H. C. DeStaebler, J. Drees, L. W . Mo and R. E. Taylor, Phys. Rev. Lett. 23, 935-939 (1969) doi:10.1103/PhysRevLett.23.935
-
[10]
Neural network determination of parton distributions: the nonsinglet case
L. Del Debbio et al. [NNPDF], JHEP 03, 039 (2007) [arXiv:hep-ph/0701127 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2007
-
[11]
New parton distribution functions from a global analysis of quantum chromodynamics
S. Dulat, T. J. Hou, J. Gao, M. Guzzi, J. Huston, P . Nadolsky, J. Pumplin, C. Schmidt, D. Stump and C. P . Y uan, Phys. Rev. 6 D 93, no.3, 033006 (2016) [arXiv:1506.07443 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[12]
L. A. Harland-Lang, A. D. Martin, P . Motylinski and R. S. Thorne, Eur. Phys. J. C 75, no.5, 204 (2015) [arXiv:1412.3989 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[13]
A. M. Cooper-Sarkar [ZEUS and H1], PoS EPS-HEP2011, 320 (2011) [arXiv:1112.2107 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[14]
Gottfried, Phys
K. Gottfried, Phys. Rev. Lett. 18, 1174 (1967)
1967
-
[15]
Amaudruzet al
P . Amaudruzet al. [New Muon], Phys. Rev. Lett. 66, 2712- 2715 (1991)
1991
-
[16]
A. Baldit et al. [NA51], Phys. Lett. B 332, 244-250 (1994) doi:10.1016/0370-2693(94)90884-2
-
[17]
K. Ackersta ffet al. [HERMES], Phys. Rev. Lett. 81, 5519- 5523 (1998) doi:10.1103 /PhysRevLett.81.5519 [arXiv:hep- ex/9807013 [hep-ex]]
- [18]
-
[19]
R. S. Towell et al. [NuSea], Phys. Rev. D 64, 052002 (2001) doi:10.1103/PhysRevD.64.052002 [arXiv:hep-ex /0103030 [hep-ex]]
-
[20]
H. Georgi and H. D. Politzer, Phys. Rev. D 14, 1829 (1976) doi:10.1103/PhysRevD.14.1829
-
[21]
B. Badelek and J. Kwiecinski, Phys. Rev. D 50, 4-8 (1994) doi:10.1103/PhysRevD.50.R4 [arXiv:hep-ph /9401314 [hep- ph]]
-
[22]
Ma, Phys
B.-Q. Ma, Phys. Lett. B 274, 111-115 (1992)
1992
-
[23]
B.-Q. Ma, A. Schafer and W . Greiner, Phys. Rev. D47, 51-55 (1993) [arXiv:hep-ph/9211202 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 1993
- [24]
-
[25]
H. Abramowicz et al. [H1 and ZEUS], Eur. Phys. J. C 75, no.12, 580 (2015) [arXiv:1506.06042 [hep-ex]]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[26]
G. Aad et al. [A TLAS], Phys. Rev. D 85, 072004 (2012) doi:10.1103/PhysRevD.85.072004 [arXiv:1109.5141 [hep- ex]]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevd.85.072004 2012
-
[27]
M. Aaboud et al. [A TLAS], Eur. Phys. J. C 77, no.6, 367 (2017) doi:10.1140 /epjc/s10052-017-4911-9 [arXiv:1612.03016 [hep-ex]]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[28]
J. Abdallah et al. [DELPHI], Eur. Phys. J. C 63, 611-623 (2009) [arXiv:0908.1023 [hep-ex]]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[29]
P . M. Nadolsky et al Phys. Rev. D 78, 013004 (2008) [arXiv:0802.0007 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[30]
G. Aad et al. [A TLAS], Phys. Rev. Lett.109, 012001 (2012) [arXiv:1203.4051 [hep-ex]]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[31]
HERAFitter, Open Source QCD Fit Project
S. Alekhin et al., HERAFitter Eur. Phys. J. C 75, no.7, 304 (2015) [arXiv:1410.4412 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[32]
xFitter (formerly HERAFitter) program, 2016, www.xfitter.org
2016
-
[33]
V . N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438- 450 (1972) IPTI-381-71
1972
-
[34]
V . N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 675- 684 (1972)
1972
-
[35]
Y . L. Dokshitzer, Sov. Phys. JETP46, 641-653 (1977)
1977
-
[36]
G. Altarelli and G. Parisi, Nucl. Phys. B 126, 298-318 (1977) doi:10.1016/0550-3213(77)90384-4
-
[37]
QCDNUM: Fast QCD Evolution and Convolution
M. Botje, Comput. Phys. Commun. 182, 490-532 (2011) doi:10.1016/j.cpc.2010.10.020 [arXiv:1005.1481 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1016/j.cpc.2010.10.020 2011
-
[38]
F. James and M. Roos, Comput. Phys. Commun. 10, 343-367 (1975) doi:10.1016/0010-4655(75)90039-9
-
[39]
Uncertainties of predictions from parton distribution functions II: the Hessian method
J. Pumplin, D. Stump, R. Brock, D. Casey, J. Hus- ton, J. Kalk, H. L. Lai and W . K. Tung, Phys. Rev. D 65, 014013 (2001) doi:10.1103 /PhysRevD.65.014013 [arXiv:hep-ph/0101032 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2001
-
[40]
A. D. Martin, W . J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63, 189-285 (2009) doi:10.1140 /epjc/s10052-009- 1072-5 [arXiv:0901.0002 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[41]
J. Gao, L. Harland-Lang and J. Rojo, Phys. Rept. 742, 1-121 (2018) doi:10.1016 /j.physrep.2018.03.002 [arXiv:1709.04922 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[42]
T. Carli, D. Clements, A. Cooper-Sarkar, C. Gwenlan, G. P . Salam, F. Siegert, P . Starovoitov and M. Sutton, Eur. Phys. J. C 66, 503-524 (2010) doi:10.1140 /epjc/s10052-010- 1255-0 [arXiv:0911.2985 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[43]
A. M. Cooper-Sarkar and K. Wichmann, Phys. Rev. D 98, no.1, 014027 (2018) doi:10.1103 /PhysRevD.98.014027 [arXiv:1803.00968 [hep-ex]]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[44]
M. Goncharov et al. [NuTeV], Phys. Rev. D 64, 112006 (2001) doi:10.1103 /PhysRevD.64.112006 [arXiv:hep- ex/0102049 [hep-ex]]
-
[45]
M. Tzanov et al. [NuTeV], Phys. Rev. D 74, 012008 (2006) doi:10.1103/PhysRevD.74.012008 [arXiv:hep-ex /0509010 [hep-ex]]
-
[46]
H. L. Lai, M. Guzzi, J. Huston, Z. Li, P . M. Nadol- sky, J. Pumplin and C. P . Y uan, Phys. Rev. D 82, 074024 (2010) doi:10.1103 /PhysRevD.82.074024 [arXiv:1007.2241 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[47]
O. Samoylov et al. [NOMAD], Nucl. Phys. B 876, 339-375 (2013) doi:10.1016 /j.nuclphysb.2013.08.021 [arXiv:1308.4750 [hep-ex]]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[48]
A. Ruzi and B.-Q. Ma, Phys. Rev. D 104, no.7, 076004 (2021) doi:10.1103 /PhysRevD.104.076004 [arXiv:2109.06854 [hep-ph]]
-
[49]
T. J. Hou, J. Gao, T. J. Hobbs, K. Xie, S. Dulat, M. Guzzi, J. Huston, P . Nadolsky, J. Pumplin and C. Schmidt, et al. Phys. Rev. D 103, no.1, 014013 (2021) doi:10.1103/PhysRevD.103.014013 [arXiv:1912.10053 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevd.103.014013 2021
-
[50]
L. A. Harland-Lang, A. D. Martin, R. Nathvani and R. S. Thorne, Eur. Phys. J. C 79, no.10, 811 (2019) doi:10.1140/epjc/s10052-019-7296-0 [arXiv:1907.02750 [hep-ph]]
-
[51]
Illuminating the photon content of the proton within a global PDF analysis
V . Bertoneet al. [NNPDF], SciPost Phys. 5, no.1, 008 (2018) doi:10.21468/SciPostPhys.5.1.008 [arXiv:1712.07053 [hep- ph]]
work page internal anchor Pith review Pith/arXiv arXiv doi:10.21468/scipostphys.5.1.008 2018
-
[52]
A. V . Manohar, P . Nason, G. P . Salam and G. Zanderighi, JHEP 12, 046 (2017) doi:10.1007 /JHEP12(2017)046 [arXiv:1708.01256 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[53]
S. Alekhin, J. Blümlein and S. Moch, Eur. Phys. J. C 78, no.6, 477 (2018) doi:10.1140 /epjc/s10052-018-5947-1 [arXiv:1803.07537 [hep-ph]]
work page internal anchor Pith review Pith/arXiv arXiv 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.