REVIEW 3 major objections 5 minor 15 references
nCTEQ PDFs at the LHC: Vector boson production in heavy ion collisions
T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Adding LHC proton-lead W/Z data to a nuclear PDF fit raises the lead strange quark distribution by about 40 percent, and the authors ask whether this shift is physical or a fitting artifact.
desk verdict A candid progress report whose headline strange shift is likely a byproduct of a restricted parameterization; the framework work is solid and worth a look. 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 nCTEQ++, a modular C++ rewrite of the nCTEQ nuclear-PDF global-fitting code. It integrates HOPPET for DGLAP evolution, APPLgrid for precomputed NLO coefficient grids, and MCFM for the W/Z matrix elements, and it exports fits in YAML for analysis. The key technical step is validation that the MCFM-generated ApplGrids are sufficiently PDF-independent that proton and nuclear PDF grids can be interchanged inside the fitting loop. The fit described here uses a deliberately limited parameter set—12 parameters instead of nCTEQ15’s 16—so that the strange PDF, normally fixed, is free to move; this restricted setup is what lets the authors attribute changes in $s(x)$ to the new LHC data.
What would settle it
Perform the same restricted fit with the three W/Z data sets that need roughly 3σ normalization shifts either removed or forced to zero shift; if the strange PDF returns to its nCTEQ15 value and the remaining W/Z $\chi^2$ stays far above 1, the 40% increase is a normalization artifact, not a physical signal.
Extended reading notes
Core claim
On its own terms, the paper's central result is that LHC pPb W±/Z data are a new, independent handle on flavor separation in nuclear parton distributions, and specifically that they prefer an enhanced strange PDF in lead. In a restricted fit that frees only 12 parameters—three for $s+\bar{s}$ and nine for the other fitted flavors—the NLO W/Z data pull $s(x)$ upward by roughly 40% in the $x$ region relevant to vector-boson production, while improving the fit to the W/Z data sets. The same exercise with a 1σ normalization shift cuts that increase roughly in half, which the authors take as a sign that the size of the strange-sea enhancement is entangled with the treatment of data normalization. The paper deliberately stops short of claiming the enhancement is physical; its stated question is whether the shift is dictated by nature or is the fit exploiting $s(x)$ because it is one of the least constrained flavors.
Load-bearing premise
The load-bearing assumption is that the normalization shifts allowed for a few LHC data sets—up to about 3σ—are genuine experimental uncertainties rather than a way to absorb theoretical or nuclear-correction errors.
Editorial extensions
If this is right
- If the result holds, LHC proton-lead W/Z measurements become a usable constraint on the strange quark PDF in nuclei, complementing fixed-target neutrino-dimuon data.
- The extracted strange PDF depends strongly on the assumed normalization uncertainties: without shifts it rises about 40%, with 1σ shifts about half that, so future fits must quote the normalization treatment alongside the PDF.
- Including heavy-ion vector-boson data can improve proton PDFs by reducing the reliance on nuclear corrections applied to heavy-target DIS data.
- With normalization shifts of up to about 3σ for a few data sets, a $\chi^2$/dof near 1 is achievable for all W/Z data, meaning the data are not mutually incompatible once normalizations are free.
Reading between the lines
- If the 40% strange enhancement survives a full 16-parameter fit and the inclusion of ALICE, LHCb, and the second CMS run, it would imply that nuclear shadowing at small $x$ is weaker than nCTEQ15 assumes; the paper does not make this claim, but its own rapidity comparison points in that direction.
- A direct cross-check would be to fit neutrino-DIS dimuon data and LHC pPb W/Z data together; if the strange PDF preferred by the two data sets disagree, the tension would locate the problem in nuclear corrections rather than in the strange flavor itself.
- The near-3σ normalization shifts on a few data sets could be absorbing missing higher-order QCD corrections or underestimated theory uncertainties; under that reading the physical strange enhancement would be closer to the 1σ-shift result, roughly 20%.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings contribution from the nCTEQ collaboration describes a new C++ implementation of the nuclear PDF fitting framework (nCTEQ++) that integrates HOPPET, APPLgrid, and MCFM, and its first application: a preliminary NLO global fit that adds LHC pPb W±/Z vector-boson production data to the existing DIS and Drell-Yan data sets of nCTEQ15. The fit uses 12 parameters, three of which describe the strange-plus-antistrange distribution. The central phenomenological result is that, when no normalization shifts are allowed for the LHC data, the fit raises the lead strange PDF by about 40% at the x values probed by W/Z production; allowing normalization shifts of up to 1σ roughly halves this shift. The paper explicitly asks whether the strange increase is dictated by nature or is an artifact of fitting one of the least constrained flavors, and states that a complete analysis with additional parameters and data is in progress.
Significance. If the reported strange-PDF shift were robust, it would be an interesting hint that LHC pPb W/Z data, together with nuclear corrections, can discriminate strange flavor in a way complementary to neutrino-DIS dimuon production. The paper has real strengths: it introduces a modular, maintainable fitting code (nCTEQ++) with validated APPLgrid-based NLO calculations, and it is admirably candid in identifying the normalization-shift dependence and the question of whether s(x) is simply being exploited as a least-constrained flavor. However, the analysis is explicitly preliminary and the central 40% figure is not accompanied by an uncertainty estimate, is obtained from a deliberately restricted parameter basis, and depends strongly on the normalization treatment. The significance of the claim is therefore conditional: the framework and the questions posed are valuable, but the quantitative conclusion is not yet established.
major comments (3)
- [Sec. 4, Fig. 4] The 40% increase of the lead strange PDF is quoted as a single number with no uncertainty. The paper advertises nCTEQ nPDFs with uncertainties, and the nCTEQ15 baseline has an uncertainty band; the fitted s(x) should be compared with that band at the relevant x values. Without such a comparison, the reader cannot tell whether a 40% shift is statistically significant or consistent with the prior uncertainty, and the abstract and conclusion currently present the number without this qualification.
- [Sec. 4, footnote 3] The restricted parameterization is a load-bearing confound for the central claim. The fit uses 12 parameters (3 for s+sbar and 9 for the remaining flavors), whereas nCTEQ15 uses 16 parameters for the non-strange set. Because the non-strange directions that are frozen here could absorb some of the LHC W/Z data's pull, the 40% strange shift may be an artifact of the reduced basis rather than a property of the data. The paper itself voices this concern in Sec. 4, but it does not provide a test, such as a fit with the full 16-parameter non-strange basis or a fit that frees the gluon and strange parameters separately, to determine how much of the shift survives.
- [Sec. 4, Fig. 3] The no-normalization fit, which produces the 40% strange shift, has an overall chi2/dof of 828/816, but the figure shows that several individual W/Z data sets have unacceptable chi2 values. The paper notes that allowing up to 1-sigma normalization shifts reduces the strange shift by roughly half. Since the headline number is taken from a fit that does not adequately describe the new data, the central claim is not robust to the normalization treatment; at minimum, the paper should give the shift with a systematic uncertainty reflecting the choice of whether and how normalization shifts are allowed.
minor comments (5)
- [Abstract] The abstract uses the word “impact” for the effect of the LHC W/Z data, but the reported chi2 numbers and PDF shifts are results of a fit that includes those data, not predictions made before the fit; a qualifier noting that the result is a preliminary fit would be appropriate in the abstract.
- [Sec. 2] There is a typo in the abstract, “strucure” for “structure,” and in Sec. 2 “Y AML” appears where “YAML” is intended; these should be corrected.
- [Fig. 3 caption] The caption contains the typos “Prelimnary” and “unconstrainted,” and the phrase “w/o Norm Penalty” should be explained, since the body text describes the shifts as constrained by luminosity uncertainties rather than by a penalty term.
- [Sec. 4, status table] The status table lists CMS Z (6235) and CMS_II W± (6232, 6234) as “grids finished not validated,” while the text says several data sets are “fully included”; the paper should specify explicitly which of the listed LHC data sets enter the reported fit and which are only shown as comparisons.
- [Fig. 4] The vertical magenta line is described only as “the central x value for pPb W±/Z production”; stating the x range or value, and the Q scale used for the PDFs, would make the plot more informative.
Circularity Check
No significant circularity: the 40% strange shift is a fit output, not a prediction; the paper explicitly flags the fit-artifact possibility.
full rationale
The paper's central result is a standard PDF re-fit: nCTEQ++ adds LHC pPb W±/Z data to the existing nCTEQ15 global fit and reports the resulting changes in χ² and in the strange PDF. This is not circular because the strange PDF is an output of the fit, not an input used to construct the data or the theory predictions. The χ² improvement from 992/816 to 828/816 is the expected consequence of adding freedom to describe newly included data, and the paper does not present this improvement as an independent prediction. The 40% increase in s(x) at the W/Z x-region is a sensitivity result, explicitly conditional on the chosen 12-parameter fit basis and on the treatment of normalization shifts; the authors themselves ask whether the fit is 'simply exploiting s(x) because that is one of the least constrained flavors.' That caveat is a limitation or model-dependence, not a logical circularity. The derivation chain relies on external inputs (LHC data, MCFM/APPLgrid/HOPPET, nCTEQ15 baseline) and makes a genuine comparison between fits with and without the new data. Self-citations to prior nCTEQ work provide the code and baseline, but the load-bearing quantitative results are computed in this paper. No equation or result is shown to reduce by construction to its own input, and no fitted parameter is renamed as a prediction. The reduced non-strange parameter basis is a legitimate methodological restriction, not a circular step.
Assumptions & free parameters
free parameters (2)
- 12 PDF shape parameters (3 for s+sbar, 9 for {g,uV,dV,ubar+dbar}) =
not reported
- Normalization shifts for LHC W/Z data sets =
up to ~3σ for a few data sets
assumptions (3)
- domain assumption Collinear QCD factorization and DGLAP evolution are valid for pPb W/Z production at NLO.
- domain assumption MCFM APPLgrids are sufficiently PDF-independent to allow interchanging proton and nuclear PDFs.
- ad hoc to paper The nCTEQ15 nuclear correction parameterization provides an adequate baseline, with 12 free parameters sufficient to capture the impact of the new data.
Cite this review
Pith. "Pith review of nCTEQ PDFs at the LHC: Vector boson production in heavy ion collisions." pith.science (2026). https://pith.science/paper/5IO7JHKT
@misc{pith2026190900452,
author = {Pith},
title = {Pith review of: nCTEQ PDFs at the LHC: Vector boson production in heavy ion collisions},
year = {2026},
howpublished = {\url{https://pith.science/paper/5IO7JHKT}},
note = {Machine review of arXiv:1909.00452}
}
read the original abstract
Extraction of the strange quark PDF is a long-standing puzzle. We use the nCTEQ nPDFs with uncertainties to study the impact of the LHC W/Z production data on both the flavor differentiation and nuclear corrections; this complements the information from neutrino-DIS data. As the proton flavor determination is dependent on nuclear corrections (from heavy target DIS, for example), LHC heavy ion measurements can also help improve proton PDFs. We introduce a new implementation of the nCTEQ code (nCTEQ++) based on C++ which has a modular strucure and enables us to easily integrate programs such as HOPPET, APPLgrid, and MCFM. Using ApplGrids generated from MCFM, we use nCTEQ++ to perform a preliminary fit including the pPb LHC W/Z vector boson data.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
nCTEQ15 - Global analysis of nuclear parton distributions with uncertain- ties in the CTEQ framework
K. Kovarik et al. “nCTEQ15 - Global analysis of nuclear parton distributions with uncertain- ties in the CTEQ framework”. In: Phys. Rev. D93.8 (2016), p. 085037. arXiv: 1509.00792 [hep-ph]
arXiv 2016
-
[2]
EPPS16: Nuclear parton distributions with LHC data
Kari J. Eskola et al. “EPPS16: Nuclear parton distributions with LHC data”. In: Eur. Phys. J. C77.3 (2017), p. 163. arXiv: 1612.05741 [hep-ph]
arXiv 2017
-
[3]
Rabah Abdul Khalek, Jacob J. Ethier, and Juan Rojo. “Nuclear parton distributions from lepton-nucleus scattering and the impact of an electron-ion collider”. In: Eur. Phys. J.C79.6 (2019), p. 471. arXiv: 1904.00018 [hep-ph]
arXiv 2019
-
[4]
Strange Quark PDFs and Implications for Drell-Yan Boson Production at the LHC
A. Kusina et al. “Strange Quark PDFs and Implications for Drell-Yan Boson Production at the LHC”. In: Phys. Rev. D85 (2012), p. 094028. arXiv: 1203.1290 [hep-ph]
arXiv 2012
-
[5]
A. M. Cooper-Sarkar and K. Wichmann. “QCD analysis of the ATLAS and CMSW± and Z cross-section measurements and implications for the strange sea density”. In: Phys. Rev. D98.1 (2018), p. 014027. arXiv: 1803.00968 [hep-ex]
arXiv 2018
-
[6]
M. Hirai, S. Kumano, and T. H. Nagai. “Determination of nuclear parton distribution functions and their uncertainties in next-to-leading order”. In: Phys. Rev. C76 (2007), p. 065207. arXiv: 0709.3038 [hep-ph]
arXiv 2007
-
[7]
Global Analysis of Nuclear Parton Distributions
Daniel de Florian et al. “Global Analysis of Nuclear Parton Distributions”. In:Phys. Rev. D85 (2012), p. 074028. arXiv: 1112.6324 [hep-ph]
arXiv 2012
-
[8]
A Higher Order Perturbative Parton Evolution Toolkit (HOPPET)
Gavin P. Salam and Juan Rojo. “A Higher Order Perturbative Parton Evolution Toolkit (HOPPET)”. In: Comput. Phys. Commun. 180 (2009), pp. 120–156. arXiv: 0804.3755 [hep-ph]
arXiv 2009
Show all 15 references
-
[9]
A posteriori inclusion of parton density functions in NLO QCD final- state calculations at hadron colliders: The APPLGRID Project
Tancredi Carli et al. “A posteriori inclusion of parton density functions in NLO QCD final- state calculations at hadron colliders: The APPLGRID Project”. In: Eur. Phys. J. C66 (2010), pp. 503–524. arXiv: 0911.2985 [hep-ph]
2010 arXiv
-
[10]
A Multi-Threaded Version of MCFM
John M. Campbell, R. Keith Ellis, and Walter T. Giele. “A Multi-Threaded Version of MCFM”. In: Eur. Phys. J.C75.6 (2015), p. 246. arXiv: 1503.06182 [physics.comp-ph]
2015 arXiv
-
[11]
Vector boson production in pPb and PbPb collisions at the LHC and its impact on nCTEQ15 PDFs
A. Kusina et al. “Vector boson production in pPb and PbPb collisions at the LHC and its impact on nCTEQ15 PDFs”. In: Eur. Phys. J.C77.7 (2017), p. 488. arXiv: 1610.02925 [nucl-th]
2017 arXiv
-
[12]
PDF Nuclear Corrections for Charged and Neutral Current Processes
I. Schienbein et al. “PDF Nuclear Corrections for Charged and Neutral Current Processes”. In: Phys. Rev. D80 (2009), p. 094004. arXiv: 0907.2357 [hep-ph]
2009 arXiv
-
[13]
Nuclear Corrections in Neutrino-Nucleus DIS and Their Compatibility with Global NPDF Analyses
K. Kovarik et al. “Nuclear Corrections in Neutrino-Nucleus DIS and Their Compatibility with Global NPDF Analyses”. In: Phys. Rev. Lett. 106 (2011), p. 122301. arXiv: 1012.0286 [hep-ph]
2011 arXiv
-
[14]
The Impact of new neutrino DIS and Drell-Yan data on large-x parton distri- butions
J. F. Owens et al. “The Impact of new neutrino DIS and Drell-Yan data on large-x parton distri- butions”. In: Phys. Rev. D75 (2007), p. 054030. arXiv: hep-ph/0702159 [HEP-PH]
2007 arXiv
-
[15]
Review of Particle Physics
M. Tanabashi et al. “Review of Particle Physics”. In: Phys. Rev. D98.3 (2018), p. 030001. 5
2018
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.