Pith. Sign in

REVIEW

Nuclear parton distributions at next to leading order

High Energy Physics - Phenomenology · Submitted 18 Nov 2003 · 21 references

Pith's verdict: First full NLO extraction of nuclear parton densities, with a real improvement in fit quality and a genuine external check — but the LO-vs-NLO difference is asserted, not quantified.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 5 major objections · 9 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/0311227 (version 1)
Pith IDMCDMLZKT, Pith's permanent record of this paper and its review
CitationD. de Florian, R. Sassot. "Nuclear parton distributions at next to leading order." arXiv:hep-ph/0311227 (2003). https://pith.science/paper/MCDMLZKT
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report
The paper

The paper claims

A next-to-leading-order global fit to nuclear deep-inelastic and Drell-Yan data, built on a convolution parameterization, describes 420 data points far better than earlier leading-order extractions, and shows that the nuclear correction fac

The reading

The paper asks whether nuclear parton densities — the momentum distributions of quarks and gluons inside a nucleon bound in a nucleus — can be extracted at next-to-leading-order accuracy in QCD, the standard already routine for free protons. The authors answer yes, and the enabling move is to write each nuclear density as the free-proton density convolved with a simple weight function, rather than as the free density multiplied by an ad hoc ratio. With three such weight functions (valence, sea, gluon) and nine parameters per nucleus they fit 420 nuclear deep-inelastic and Drell-Yan data points, reaching $300.15$ units of $χ^2$ for 393 degrees of freedom, against roughly 630 for the widely used leading-order set on the same data. The practical payoff is a warning: the nuclear correction factors extracted at leading order and at next-to-leading order differ by amounts comparable to the nuclear effects themselves at low $Q^2$, so mixing leading-order nuclear factors with next-to-leading-order proton densities is not safe.

What carries the argument

The convolution ansatz $f_i^A(x_N,Q_0^2)=int_{x_N}^{A} rac{dy}{y}, W_i(y,A,Z), f_i(x_N/y,Q_0^2)$, where $W_i$ is an effective density of nucleons carrying a fraction $y/A$ of the nucleus longitudinal momentum. It buys three things: the full kinematic range $0<x_N<A$ instead of truncating at $x_N=1$; closed-form Mellin moments $hat{W}_i(N,A,Z)$ so that evolution and next-to-leading-order coefficient functions are multiplications in moment space; and a physical reading of nuclear effects as a rebalancing of momentum fractions between partons. Non-singlet combinations make $hat{W}_v$ scale independent, while sea and gluon weights mix and must be fixed at an initial scale. Momentum and charge conservation pin two of the parameters, leaving nine per nucleus, with a smooth mass-number dependence $ε_i=γ_i+λ_i A^{δ_i}$.

Extended reading notes

Core claim

Two claims carry the paper. First, the convolution parameterization makes a full next-to-leading-order global fit feasible at all: because the weight functions have closed-form Mellin moments, evolution and the hard coefficient functions reduce to ordinary multiplication in moment space, which the earlier multiplicative-ratio parameterizations with many parameters could not support. The resulting fit describes the NMC and SLAC-E139 structure-function ratios and the E772 Drell-Yan rates far better than the leading-order EKS98 set did on the same points, with $χ^2/mathrm{d.o.f.} = 299.91/393$. Second, and more consequential for practice, the leading-order and next-to-leading-order fits, although of comparable quality, give nuclear correction ratios $R_i(x_N,Q^2,A)$ that differ substantially — most in sea quarks and gluons at small $x_N$ in heavy nuclei — with the discrepancy shrinking as $Q^2$ grows. The authors use this to argue that leading-order nuclear factors should not be bolted onto next-to-leading-order proton densities, a common shortcut.

Load-bearing premise

The fit assumes nuclear modification is isospin-invariant and flavour-blind, so one weight function multiplies both up and down valence quarks and another multiplies both anti-up and anti-down sea quarks, making the extracted weights an average that could bias predictions for neutron-rich nuclei.

Editorial extensions

If this is right

  • Nuclear correction ratios from this fit can be multiplied onto any modern next-to-leading-order proton density set to build nuclear densities, provided both come from the same order in QCD.
  • Predictions for neutral-pion production in deuteron-gold collisions at $sqrt{s_{NN}}=200$ GeV agree with the PHENIX data within experimental errors except at the lowest $p_T$, and although the leading-order and next-to-leading-order cross sections differ by roughly 50 percent, the ratio to proton-proton largely cancels.
  • A variant set with stronger gluon shadowing ($R_g^{Au}=0.75$ at $x_N=0.001$, $Q^2=5$ GeV$^2$) costs about 250 units of $χ^2$, bounding how much small-$x$ gluon shadowing the present data allow.
  • The authors argue that a much stronger reduction of $R_{dAu}$ than their prediction at $p_T$ near 2 GeV would require physics beyond perturbative QCD.
  • Gluon densities remain the least constrained ingredient: among the data used, only the $Q^2$ dependence of $F_2^{Sn}/F_2^{C}$ is sensitive to them.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The paper's own side-by-side comparison of leading-order and next-to-leading-order ratios suggests a broader audit: any phenomenological study that imports a nuclear correction factor should state the order at which it was extracted, and the mismatch size shown here at low $Q^2$ is a rough template for how large that error can be.
  • Because the weight functions are convolutions rather than pointwise ratios, the framework extends naturally to nuclei where data are absent, so long as the mass-number dependence of the nine parameters is trusted; the risk is that this extrapolation does more work in the lead-lead and deuteron-gold predictions than the fit itself constrains.
  • The flavour-blind sea weight is the natural place for future data to bite: neutrino-nucleus scattering and Drell-Yan on neutron-rich targets would test whether an isovector piece is hiding inside the fitted average.
  • The bounded gluon shadowing of the variant set implies that sharper small-$x$ gluon constraints will have to come from future electron-ion colliders or forward hadron measurements rather than from the existing deep-inelastic ratios.
Original abstract

We perform a next to leading order QCD global analysis of nuclear deep inelastic scattering and Drell-Yan data using the convolution approach to parameterize nuclear parton densities. We find both a significant improvement in the agreement with data compared to previous extractions, and substantial differences in the scale dependence of nuclear effects compared to leading order analyses.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 2,294 input tokens, $0.0001 reported.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.07 · human participants 0.01 · controlled comparison 0.33 · randomized assignment 0.01 · states quantitative result 0.08 · compares with prior results 0.97 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

We perform a next to leading order QCD global analysis of nuclear deep inelastic scattering and Drell-Yan data using the convolution approach to parameterize nuclear parton densities. We find both a significant improvement in the agreement with data compared to previous extractions, and substantial differences in the scale dependence of nuclear effects compared to leading order analyses.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We perform a next to leading order QCD global analysis of nuclear deep inelastic scattering and Drell-Yan data using the convolution approach to parameterize nuclear parton densities. method 0.50 0.96 0.08 asserted 0.99
s1 We find both a significant improvement in the agreement with data compared to previous extractions, and substantial differences in the scale dependence of nuclear effects compared to leading order analyses. result 1.00 0.96 0.79 asserted 1.00

What kind of paper the classifier reports this is: simulation study 0.32 · evidence from analysis of existing data 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References
How the literature uses it

73 arXiv papers cite this in 97 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 51 sentences without a confident label

Contests it 1 paper
Uses its method or tool 4 papers
Compares with it 7 papers
Cites it as support 1 paper
Cites it as background 29 papers
  • The fit was done only by eye. * HKM [ref] & HKN [ref] which are LO QCD analyses as well, were the first ones to exploit the automated $χ^2$-minimization procedure and to estimate nPDF uncertainties. * nDS [ref] brought the global analysis of nPDFs to the next-to-leading order level in pQCD.

    A global DGLAP analysis of nuclear PDFs · Introduction · 0.90

  • Since the first of the nPDF sets, EKS98 [ref] , the procedure has been improved by performing the analysis at next-to-leading order (NLO) [ref] and by making uncertainty estimates [ref] in analogy with the free proton case.

    An improved global analysis of nuclear parton distribution functions including RHIC data · Introduction · 0.82

  • Historically, the inclusive pion production at RHIC was the first direct gluon constraint in the fits on nuclear PDFs [ref] .

    EPPS21: A global QCD analysis of nuclear PDFs · Inclusive pion production · 0.97

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.