REVIEW 3 major objections 6 minor 1 cited by
Next-to-leading order evolution of structure functions without PDFs
T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read This paper derives a closed set of NLO evolution equations for six measurable DIS structure functions, replacing parton distribution functions with the observables themselves.
desk verdict First NLO six-observable physical-basis DGLAP evolution in momentum space, carefully derived but with a load-bearing boundary relation (C1) that needs an independent check before I'd fully trust the kernels. 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 machinery is the six-observable physical basis together with the differential operator $\hat P(x) = x^2 \frac{d^2}{dx^2} - 2x \frac{d}{dx} + 2$, which is the exact inverse of the leading-order gluon coefficient function $C^{(1)}_{FLg}$ in the longitudinal structure function. This operator lets the paper invert the PDF–structure-function relation perturbatively, expressing the gluon and quark distributions as series in $\alpha_s$ of structure functions and their derivatives. Substituting these inverted PDFs into the standard NLO DGLAP equations and using partial integration to move the derivatives onto the structure functions produces the evolution kernels listed in Appendix C; the boundary terms left by partial integration are removed using the relation $fF'_{L}(1) = -0.003\,\frac{\alpha_s}{2\pi} fF'_2(1) + O(\alpha_s^2)$.
What would settle it
Take any NLO PDF set, compute $F_2$ and $F_L$ at several scales, and evaluate the combination $fF'_{L}(1) + 0.003\,\frac{\alpha_s}{2\pi} fF'_2(1)$; if it does not vanish at order $\alpha_s^2$ as a function of $Q^2$, the boundary handling is incorrect, and the physical-basis evolution will differ from conventional PDF evolution by a term of order $\alpha_s$ rather than $\alpha_s^2$.
Extended reading notes
Core claim
The central result is the closed system of equations (32)–(37), accurate to next-to-leading order in $\alpha_s$, for the six structure functions $F_2$, $F_3$, $\Delta F_2^W$, $F_3^{W^-}$, $F_{2c}^{W^-}$, and $F_L$. The paper obtains these equations by inverting the NLO coefficient-function relations to express PDFs as power series in $\alpha_s$ whose coefficients are structure functions and their $x$-derivatives, substituting these inverted PDFs into the conventional DGLAP splitting functions, and partially integrating the differential operator that arises from the inversion. The factorization-scale and scheme dependence cancel between coefficient functions and splitting functions, so the resulting evolution kernels are factorization-scale and scheme independent by construction. When evolved numerically from initial conditions given by a standard NLO PDF set, the physical-basis structure functions track the PDF-evolved ones with differences that are formally higher order in $\alpha_s$; the deviations are typically below ten percent and smallest for the valence-dominated observables.
Load-bearing premise
The construction stands on the boundary relation $fF'_{L}(1) = -0.003\,\frac{\alpha_s}{2\pi} fF'_2(1) + O(\alpha_s^2)$, which eliminates the surface terms generated when the differential operator is partially integrated; if this relation is wrong or misses higher-order terms, the evolution kernels listed in Appendix C no longer reproduce NLO DGLAP evolution.
Editorial extensions
If this is right
- Physical-basis DGLAP evolution is factorization-scale and factorization-scheme independent, so the only remaining scheme choice in a calculation is the renormalization scheme of $\alpha_s$.
- Initial conditions for the physical basis are the same at every perturbative order, in contrast to PDFs, which must be refitted at each order and scheme.
- Because the inverted PDFs are expressed in terms of structure functions, other PDF-dependent cross sections can be rewritten in the physical basis; the paper demonstrates this by computing the non-basis structure function $F_2^{W^+}$ within about twenty percent of the PDF-based value.
- The iterative inversion procedure extends straightforwardly to higher orders in $\alpha_s$, so the approach can be pushed beyond NLO, with the caveat that more nested differential operators appear at each order.
- Extending the framework to a variable-flavour-number scheme and including heavy-quark masses would be needed before the physical basis can serve as a full phenomenology-ready alternative to PDF evolution.
Reading between the lines
- If the boundary relation at $x=1$ is robust, the same partial-integration strategy could be carried to NNLO, but the appearance of up to three nested differential operators at that order makes accurate numerical handling of high derivatives the main technical challenge.
- The negative large-$x$ values of the physical-basis inverted NLO gluon suggest that positivity constraints derived for MS PDFs do not transfer directly to physical-basis PDF counterparts; the paper notes this undercuts a positivity argument made elsewhere for MS PDFs.
- The physical basis offers a concrete way to estimate factorization-scheme uncertainties in LHC cross sections: comparing physical-basis predictions with conventional MS predictions quantifies scheme dependence without needing to construct alternate PDF fits.
- A natural testable extension is to apply the physical-basis construction to Drell-Yan or other hadronic processes, checking that the scheme and scale dependences cancel in those cross sections exactly as they do for DIS structure functions.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper constructs a physical-basis DGLAP evolution at next-to-leading order for six DIS structure functions. Starting from the MS expressions for the structure functions in terms of PDFs, the authors invert these relations iteratively to express PDFs in terms of structure functions, substitute into the NLO DGLAP equations, and obtain a closed set of evolution equations (32)-(37) for the structure functions. After partial integration of the differential operator P-hat, they present explicit evolution kernels in Appendix C. The numerical solution for CT18NLO-based initial conditions is compared with standard PDF evolution, with differences typically at the 10-20% level. The paper also discusses inverted PDFs, sum-rule violations, and the computation of a non-basis structure function. The central claim is that the physical-basis evolution is factorization-scale and scheme independent by construction.
Significance. If correct, this work provides the most extensive NLO physical-basis evolution to date, with potential applications to global fits of DIS data and to cross-section predictions in terms of physical observables. A notable strength is the explicit listing of all evolution kernels and the honest discussion of limitations, including fixed-flavor scheme, massless quarks, and the violation of sum rules by the perturbatively inverted PDFs. The paper also explicitly footnotes a discrepancy with the two-observable limit of Ref. [32], which is commendable but, as detailed below, needs to be resolved. The algebraic derivation is presented in detail and the comparison with CT18 provides a sanity check, but the validation is incomplete due to the asserted boundary relation and the absence of a direct consistency test of the operator and kernel forms.
major comments (3)
- [Appendix C, Eq. (C1)] The boundary relation fF'_L(1) = -0.003 (alpha_s/2pi) fF'_2(1) + O(alpha_s^2) is asserted with only a parenthetical remark that the factor -0.003 multiplies the delta function in C^(2)_FLNS. The text states that the evolution kernels are obtained using this relation, so it is load-bearing for the central result. However, the relation involves the x-derivatives of the structure functions at the endpoint, and it is not shown how the delta-function coefficient alone determines these derivatives; the regular part of C^(2)_FLNS and the O(alpha_s) correction to Xi could contribute to fF'_L(1). Please provide a complete derivation from Eq. (11), including a distributional treatment of the endpoint, and verify the relation numerically on a test structure-function input. Without this, the kernels in Appendix C are not fully justified.
- [Footnote 1] The paper reports an unresolved discrepancy with the two-observable (F2, FL) limit of Ref. [32], attributing it to a possible inconsistency in that reference. Since the two-observable limit is the minimal independent check of the physical-basis construction, and since the same boundary handling is involved, the authors should either demonstrate the inconsistency in Ref. [32] (for example by reproducing its Mellin-space computation) or reconcile the two results. As it stands, the discrepancy leaves the NLO derivation without an independent limiting-case validation.
- [Section III] The numerical solution is presented only via the kernel form (41)-(46), with no independent check that these kernels reproduce the original operator equations (32)-(37). Given the nontrivial partial integration and the use of Eq. (C1), a consistency test comparing dF_i/d log Q^2 computed directly from (32)-(37) (with explicit derivatives of the structure functions) and from the Appendix C kernels for the same input would materially strengthen the validation. Please add such a test, at least for a representative structure function (e.g., F2 or FL) over the x-range used in the figures.
minor comments (6)
- [Appendix A, Eq. (A3)] The notation 'log(1 - x)^2' is ambiguous; please use log^2(1-x) or [log(1-x)]^2 to avoid confusion with log[(1-x)^2].
- [Section III] The numerical method is not described (x-grid, interpolation for derivatives, treatment of plus-distributions and endpoint terms). A brief description would improve reproducibility.
- [Figure 2 caption] The phrase 'colourful curves' is informal; consider 'colored curves' or a more precise description of the lines and markers.
- [Footnote 1] The important caveat about the discrepancy with Ref. [32] is placed in a footnote; it would be better moved to the main text or an appendix, since it bears directly on the validation of the results.
- [Appendix C, Eq. (C17)] The notation (C ⊗ G^(1))_{Fi} is introduced but not fully explained for all Fi in the list; a short definition of how these terms arise after partial integration would aid the reader.
- [General] The paper would benefit from a statement on the availability of the numerical code or input files; as it stands, reproducing the figures requires reimplementation of long expressions.
Circularity Check
No significant circularity: the NLO physical-basis kernels are an algebraic rearrangement of standard MS coefficient functions and DGLAP splitting functions, benchmarked against CT18; self-citations are not load-bearing.
full rationale
The derivation chain starts from standard, externally published NLO MS coefficient functions (Refs. [38,39]) and the conventional DGLAP splitting functions, inverts them in a perturbative expansion (Eqs. (24)-(28), Appendix B), and substitutes into the DGLAP equations to obtain Eqs. (32)-(37). No parameter is fitted to the comparison; the numerical comparison uses CT18NLO NF3 as an external input to set initial conditions and as an independent benchmark. The only self-citations are to the authors' previous LO paper [37] for the LO inverse operator and to [33] for scale-variation consistency; the LO inverse is rederived explicitly (Eqs. (17)-(19), (B1)-(B16)), so the citation is not load-bearing. The boundary relation Eq. (C1), fF'L(1) = -0.003 (alpha_s/2pi) fF'2(1) + O(alpha_s^2), is derived from the delta-function term of the external NLO coefficient function C^(2)_FLNS, not fitted or posited as a prediction; it is a possible correctness risk (as is the unresolved discrepancy with Ref. [32] noted in footnote 1), but not a circular step. The central claim of scheme/scale independence is a consequence of the standard factorization argument, not an input assumed from the authors' prior work. Overall score 1 reflects only the presence of minor, non-load-bearing self-citations.
Assumptions & free parameters
free parameters (1)
- s-quark PDF rescaling factor in Fig. 3 =
0.2635
assumptions (5)
- standard math NLO DGLAP splitting functions in MS scheme (Ref. [39])
- standard math NLO DIS coefficient functions in MS scheme (Refs. [38,39])
- domain assumption The differential operator P-hat is the exact inverse of the LO coefficient function C^(1)_FLg
- ad hoc to paper Boundary relation Eq. (C1) between fF'L(1) and fF'2(1)
- domain assumption Fixed 3-flavour massless quark scheme with s = sbar and no CKM mixing
Cite this review
Pith. "Pith review of Next-to-leading order evolution of structure functions without PDFs." pith.science (2026). https://pith.science/paper/2VVPJ6O6
@misc{pith2026241209589,
author = {Pith},
title = {Pith review of: Next-to-leading order evolution of structure functions without PDFs},
year = {2026},
howpublished = {\url{https://pith.science/paper/2VVPJ6O6}},
note = {Machine review of arXiv:2412.09589}
}
read the original abstract
We formulate and numerically solve the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi~(DGLAP) evolution equations at next-to-leading order in perturbation theory directly for a basis of 6 physical, observable structure functions in deeply inelastic scattering. By expressing the evolution in the physical basis one evades the factorization scale and scheme dependence. Working in terms of observable quantities, rather than parametrizing and fitting unobservable parton distribution functions (PDFs), provides an unambiguous way to confront predictions of perturbative Quantum Chromodynamics with experimental measurements. We compare numerical results for the DGLAP evolution for structure functions in the physical basis to the conventional evolution with PDFs.
Figures
Forward citations
Cited by 1 Pith paper
-
Factorisation schemes for proton PDFs
All common NLO factorisation schemes for proton PDFs are special cases of a single discrete family, and scheme choice can shift LHC jet predictions by up to 20-30%.
Reference graph
Works this paper leans on
-
[32]
M. Hentschinski and M. Stratmann, On the Practical Application of Physical Anomalous Dimensions , arXiv:1311.2825 [hep-ph]
-
[1]
R. Abdul Khalek et. al. , Science Requirements and Detector Concepts for the Electron-Ion Collider: EIC Yellow Report , Nucl. Phys. A 1026 (2022) 122447 [ arXiv:2103.05419 [physics.ins-det]]
arXiv 2022
-
[2]
LHeC, FCC-he Study Group collaboration, P. Agostini et. al. , The Large Hadron-Electron Collider at the HL-LHC , J. Phys. G 48 (2021) no. 11 110501 [ arXiv:2007.14491 [hep-ex]]
arXiv 2021
- [3]
-
[4]
H1, ZEUS collaboration, H. Abramowicz et. al., Combination of measurements of inclusive deep inelastic e±p scattering cross sections and QCD analysis of HERA data , Eur. Phys. J. C 75 (2015) no. 12 580 [ arXiv:1506.06042 [hep-ex]]
arXiv 2015
-
[5]
H1 collaboration, V. Andreev et. al. , Measurement of inclusive ep cross sections at high Q2 at √s = 225 and 252 GeV and of the longitudinal proton structure function FL at HERA , Eur. Phys. J. C 74 (2014) no. 4 2814 [ arXiv:1312.4821 [hep-ex]]
arXiv 2014
-
[6]
Deep inelastic cross-section measurements at large y with the ZEUS detector at HERA
ZEUS collaboration, H. Abramowicz et. al., Deep inelastic cross-section measurements at large y with the ZEUS detector at HERA , Phys. Rev. D 90 (2014) no. 7 072002 [ arXiv:1404.6376 [hep-ex]]
work page Pith review arXiv 2014
-
[7]
J. P. Berge et. al. , A Measurement of Differential Cross-Sections and Nucleon Structure Functions in Charged Current Neutrino Interactions on Iron , Z. Phys. C 49 (1991) 187
work page 1991
Show all 45 references
-
[8]
Onengut et
CHORUS collaboration, G. Onengut et. al. , Measurement of nucleon structure functions in neutrino scattering , Phys. Lett. B 632 (2006) 65
2006
-
[9]
Tzanov et
NuT eVcollaboration, M. Tzanov et. al. , Precise measurement of neutrino and anti-neutrino differential cross sections , Phys. Rev. D 74 (2006) 012008 [ arXiv:hep-ex/0509010]
2006 arXiv
-
[10]
Goncharov et
NuT eVcollaboration, M. Goncharov et. al., Precise Measurement of Dimuon Production Cross-Sections in νµ Fe and ¯νµ Fe Deep Inelastic Scattering at the Tevatron. , Phys. Rev. D 64 (2001) 112006 [ arXiv:hep-ex/0102049]
2001 arXiv
-
[11]
Samoylov et
NOMAD collaboration, O. Samoylov et. al. , A Precision Measurement of Charm Dimuon Production in Neutrino Interactions from the NOMAD Experiment , Nucl. Phys. B 876 (2013) 339 [ arXiv:1308.4750 [hep-ex]]
2013 arXiv
-
[12]
Y. L. Dokshitzer, Calculation of the Structure Functions for Deep Inelastic Scattering and e+e− Annihilation by Perturbation Theory in Quantum Chromodynamics. , Sov. Phys. JETP 46 (1977) 641
1977
-
[13]
V. N. Gribov and L. N. Lipatov, Deep inelastic e p scattering in perturbation theory , Sov. J. Nucl. Phys. 15 (1972) 438
1972
-
[14]
V. N. Gribov and L. N. Lipatov, e+ e- pair annihilation and deep inelastic e p scattering in perturbation theory , Sov. J. Nucl. Phys. 15 (1972) 675. 25
1972
-
[15]
Altarelli and G
G. Altarelli and G. Parisi, Asymptotic Freedom in Parton Language , Nucl. Phys. B 126 (1977) 298
1977
-
[16]
Cacciari and N
M. Cacciari and N. Houdeau, Meaningful characterisation of perturbative theoretical uncertainties , JHEP 09 (2011) 039 [arXiv:1105.5152 [hep-ph]]
2011 arXiv
-
[17]
David and G
A. David and G. Passarino, How well can we guess theoretical uncertainties? , Phys. Lett. B 726 (2013) 266 [arXiv:1307.1843 [hep-ph]]
2013 arXiv
-
[18]
Bagnaschi, M
E. Bagnaschi, M. Cacciari, A. Guffanti and L. Jenniches, An extensive survey of the estimation of uncertainties from missing higher orders in perturbative calculations , JHEP 02 (2015) 133 [ arXiv:1409.5036 [hep-ph]]
2015 arXiv
-
[19]
Abdul Khalek et
NNPDF collaboration, R. Abdul Khalek et. al. , A first determination of parton distributions with theoretical uncertainties, Eur. Phys. J. C (2019) 79:838 [ arXiv:1905.04311 [hep-ph]]
2019 arXiv
-
[20]
Furmanski and R
W. Furmanski and R. Petronzio, Lepton - Hadron Processes Beyond Leading Order in Quantum Chromodynamics , Z. Phys. C 11 (1982) 293
1982
-
[21]
Armesto, T
N. Armesto, T. Lappi, H. M¨ antysaari, H. Paukkunen and M. Tevio, Signatures of gluon saturation from structure-function measurements, Phys. Rev. D 105 (2022) no. 11 114017 [ arXiv:2203.05846 [hep-ph]]
2022 arXiv
-
[22]
J. C. Collins, D. E. Soper and G. F. Sterman, Factorization of Hard Processes in QCD , Adv. Ser. Direct. High Energy Phys. 5 (1989) 1 [ arXiv:hep-ph/0409313]
1989 arXiv
-
[23]
J. Gao, L. Harland-Lang and J. Rojo, The Structure of the Proton in the LHC Precision Era , Phys. Rept. 742 (2018) 1 [arXiv:1709.04922 [hep-ph]]
2018 arXiv
-
[24]
Kovaˇ r ´ ık, P
K. Kovaˇ r ´ ık, P. M. Nadolsky and D. E. Soper,Hadronic structure in high-energy collisions , Rev. Mod. Phys. 92 (2020) no. 4 045003 [ arXiv:1905.06957 [hep-ph]]
2020 arXiv
-
[25]
J. J. Ethier and E. R. Nocera, Parton Distributions in Nucleons and Nuclei , Ann. Rev. Nucl. Part. Sci. 70 (2020) 43 [arXiv:2001.07722 [hep-ph]]
2020 arXiv
-
[26]
Klasen and H
M. Klasen and H. Paukkunen, Nuclear PDFs After the First Decade of LHC Data , arXiv:2311.00450 [hep-ph]
-
[27]
McGowan, T
J. McGowan, T. Cridge, L. A. Harland-Lang and R. S. Thorne, Approximate N3LO parton distribution functions with theoretical uncertainties: MSHT20aN 3LO PDFs , Eur. Phys. J. C 83 (2023) no. 3 185 [ arXiv:2207.04739 [hep-ph]]. [Erratum: Eur.Phys.J.C 83, 302 (2023)]
2023 arXiv
-
[28]
NNPDF collaboration, R. D. Ball et. al. , The path to N 3LO parton distributions , Eur. Phys. J. C 84 (2024) no. 7 659 [arXiv:2402.18635 [hep-ph]]
2024 arXiv
-
[29]
Delorme, A
S. Delorme, A. Kusina, A. Si´ odmok and J. Whitehead, Factorisation schemes for proton PDFs , arXiv:2501.18289 [hep-ph]
-
[30]
Catani, Physical anomalous dimensions at small x , Z
S. Catani, Physical anomalous dimensions at small x , Z. Phys. C 75 (1997) 665 [ arXiv:hep-ph/9609263]
1997 arXiv
-
[31]
Blumlein, V
J. Blumlein, V. Ravindran and W. L. van Neerven, On the Drell-Levy-Yan relation to O(α2 s ), Nucl. Phys. B 586 (2000) 349 [arXiv:hep-ph/0004172]
2000 arXiv
-
[33]
L. A. Harland-Lang and R. S. Thorne, On the Consistent Use of Scale Variations in PDF Fits and Predictions , Eur. Phys. J. C 79 (2019) no. 3 225 [ arXiv:1811.08434 [hep-ph]]
2019 arXiv
-
[34]
Bl¨ umlein and M
J. Bl¨ umlein and M. Saragnese,The N 3LO scheme-invariant QCD evolution of the non-singlet structure functions F NS 2 (x, Q2) and gNS 1 (x, Q2), Phys. Lett. B 820 (2021) 136589 [ arXiv:2107.01293 [hep-ph]]
2021 arXiv
-
[35]
W. L. van Neerven and A. Vogt, NNLO evolution of deep inelastic structure functions: The Nonsinglet case , Nucl. Phys. B 568 (2000) 263 [ arXiv:hep-ph/9907472]
2000 arXiv
-
[36]
Ruiz Arriola, NLO evolution for large scale distances, positivity constraints and the low-energy model of the nucleon , Nucl
E. Ruiz Arriola, NLO evolution for large scale distances, positivity constraints and the low-energy model of the nucleon , Nucl. Phys. A 641 (1998) 461
1998
-
[37]
Lappi, H
T. Lappi, H. M¨ antysaari, H. Paukkunen and M. Tevio,Evolution of structure functions in momentum space , Eur. Phys. J. C 84 (2024) no. 1 84 [ arXiv:2304.06998 [hep-ph]]
2024 arXiv
-
[38]
S. Moch, J. A. M. Vermaseren and A. Vogt, The Longitudinal structure function at the third order , Phys. Lett. B 606 (2005) 123 [ arXiv:hep-ph/0411112]
2005 arXiv
-
[39]
R. K. Ellis, W. J. Stirling and B. R. Webber, QCD and collider physics , vol. 8. Cambridge University Press, 2, 2011
2011
-
[40]
Y. V. Kovchegov and E. Levin, Quantum Chromodynamics at High Energy , vol. 33. Oxford University Press, 2013
2013
-
[41]
T.-J. Hou et. al. , New CTEQ global analysis of quantum chromodynamics with high-precision data from the LHC , Phys. Rev. D 103 (2021) no. 1 014013 [ arXiv:1912.10053 [hep-ph]]
2021 arXiv
-
[42]
Buckley, J
A. Buckley, J. Ferrando, S. Lloyd, K. Nordstr¨ om, B. Page, M. R¨ ufenacht, M. Sch¨ onherr and G. Watt,LHAPDF6: parton density access in the LHC precision era , The European Physical Journal C 75 (mar, 2015)
2015
-
[43]
Candido, S
A. Candido, S. Forte, T. Giani and F. Hekhorn, On the positivity of MS parton distributions , Eur. Phys. J. C 84 (2024) no. 3 335 [ arXiv:2308.00025 [hep-ph]]
2024 arXiv
-
[44]
Barontini, A
A. Barontini, A. Candido, F. Hekhorn, G. Magni and R. Stegeman, An FONLL prescription with coexisting flavor number PDFs, JHEP 10 (2024) 004 [ arXiv:2408.07383 [hep-ph]]
2024 arXiv
-
[45]
Gluck and E
M. Gluck and E. Reya, Deep Inelastic Quantum Chromodynamic Charm Leptoproduction , Phys. Lett. B 83 (1979) 98
1979
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.