REVIEW 4 major objections 6 minor 113 references
Candia-v2: Logarithmic expansions for DGLAP evolution in $x$-space
T0 review · 4 major / 6 minor · reviewed 2026-08-02 · deepseek-v4-flash
Pith's one-line read Candia-v2 solves the DGLAP evolution equations at N3LO accuracy in QCD, directly in x-space.
desk verdict Worth refereeing as a software release, but the N3LO numbers in Table 2 are not backed by any external check. 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 the logarithmic-expansion ansatz with the L_i logs: one log of the alpha_s ratio, one log of a quadratic in alpha_s, one log of a linear root term, and an arctangent coming from the roots of the four-loop beta-function. The coefficient functions a_i(x) are built from convolutions and fixed by the four recursion relations, which generate each higher-order coefficient from lower-index convolutions with the splitting-function kernels. Numerically, convolutions are evaluated by Gauss-Legendre quadrature on a three-segment x-grid, with singular plus-distribution and delta-function pieces factored out and precomputed per active flavor number.
What would settle it
Run Candia-v2 on the Les Houches toy-model input at mu0^2 = 2 GeV^2 with the same aN3LO splitting functions and operator matrix elements, and compare its evolved xg and xuv at mu = 100 GeV, plus the c(-) and b(-) asymmetries, against an independent evolution code using a different discretization, such as Mellin-space with contour inversion. If the difference at the tabulated x points is much larger than the claimed tolerance, or if doubling the Gauss nodes or grid segments materially changes the output, the numerical claim fails even though the analytic recursion may be correct.
Extended reading notes
Core claim
The paper argues that the DGLAP solution can be built in x-space from the ansatz f(x, mu2) = product over i of [sum_n a_i(x)^n/n! L_i^n] convoluted with the initial PDF, with the logarithmic functions L_i depending only on alpha_s, the starting coupling, and the roots of the four-loop beta-function, and with the coefficient functions fixed by recursion relations. Candia-v2 implements this at N3LO: exact four-loop splitting functions in the non-singlet sector, the currently best available approximations in the singlet sector, and three-loop operator matrix elements for heavy-quark threshold matching. The paper presents evolved gluon and valence densities, heavy-flavor asymmetries, and subtrac
Load-bearing premise
The load-bearing premise is that discretizing the convolutions by Gauss-Legendre quadrature on a hand-chosen three-segment grid reproduces the continuum DGLAP solution to the implicit precision, despite there being no grid-refinement or convergence study in the paper and the repository's 'tests' being benchmarking files rather than standard tests.
Editorial extensions
If this is right
- Evolved densities at approximate N3LO with all currently known perturbative ingredients become available in a single open-source package, suitable for LHC and HL-LHC scale predictions.
- Heavy-flavor asymmetries c(-) and b(-), which first appear at NNLO, now include N3LO evolution and matching corrections, and Fig. 2 shows their size.
- Subtraction and residual PDFs for charm and bottom can be generated on demand, giving reusable, process-independent pieces for general-mass variable-flavor-number schemes.
- When exact four-loop singlet splitting functions become available, the recursion structure allows them to be inserted without restructuring the algorithm.
- LHAPDF input and output let the code be embedded in larger fitting and prediction environments as an x-space alternative to Mellin-space evolution.
Reading between the lines
- Editorial: the paper does not report a grid-refinement or quadrature-order convergence study, so a user should verify that the default 100/50/25 grid and 50 Gauss points are sufficient for the x region of interest before trusting the tabulated digits.
- Editorial: because the repository's 'tests' are described as benchmarking and plotting files rather than standard tests, an independent comparison of Candia-v2 with a differently discretized evolution code on identical splitting-function inputs is the direct way to check the N3LO numbers.
- Editorial: the same logarithmic-expansion machinery should carry over to polarized parton densities and to QED corrections, since the recursion relations depend only on the beta-function and splitting-function inputs.
- Editorial: the exact-ansatz recursion cost scales exponentially with perturbative order, so the truncated ansatz is the practical route at higher orders; the code already exposes both, suggesting a path toward N4LO when the required kernels exist.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents Candia-v2, an open-source C++/Fortran reimplementation of the x-space Candia algorithm for DGLAP evolution. It claims to generalize the previous version to (approximate) N3LO accuracy by implementing four-loop non-singlet splitting functions, state-of-the-art approximations for singlet splitting functions, three-loop operator-matrix-element matching via libome, and optional LHAPDF interfaces. The algorithm is based on logarithmic expansions and recursion relations for coefficients in x-space. Results are shown for the Les Houches toy model evolved from 2 GeV^2 to 10^4 GeV^2 in the VFNS, including scale dependence and heavy-quark subtraction/residual PDFs. The paper also provides a detailed program description, API overview, and a list of bug fixes from Candia-v1.
Significance. If the numerical implementation is correct, Candia-v2 would be a useful independent x-space evolution tool at approximate N3LO, with practical features such as subtraction PDFs, LHAPDF I/O, parallelization, and modular design. The paper explicitly uses exact four-loop non-singlet splitting functions, state-of-the-art singlet approximations, and libome for three-loop OMEs, and it openly discusses known v1 bugs. These are strengths that make the package potentially valuable to the community. The main risk is the absence of any quantitative validation of the numerical core, which underlies all reported N3LO results.
major comments (4)
- [Eq. (24)] The recursion relations in Eq. (24) contain the undefined symbol gamma in the denominators of Z21, Z22, Z31, Z32, and Z43. Since these relations are the core of the N3LO implementation, the algorithm is not self-contained as written. Please define gamma, presumably gamma = sqrt(4cbar - bbar^2) as suggested by Eq. (21), or explicitly refer to the corresponding definition in Ref. [6].
- [Sec. 3, Table 2] No external validation is provided for the central N3LO-accuracy claim. The results are compared only across LO/NLO/NNLO/N3LO within the same code. The paper cites the aN3LO benchmarking paper Ref. [75] but does not compare against it or any other evolution code (e.g., PEGASUS, HOPPET, APFEL, EKO). The repository's tests are explicitly 'not standard tests' (Sec. 5.2), and the disclosed v1 bugs (Sec. 5.5) show that implementation errors are a real risk. A direct numerical comparison with at least one independent evolution code is essential.
- [Sec. 5.4.2, Sec. 3] The convolution quadrature and the iteration/truncation of the logarithmic series are not validated. All convolutions use a fixed three-segment grid with 100/50/25 points, 50 Gauss-Legendre nodes, and 4 interpolation points (Sec. 5.4.2). The example uses iterations=10, trunc_idx=10, with no convergence study in these parameters. Since the integrands involve plus-distributions and delta-functions, and the small-x region is highly sensitive, the numerical error could be comparable to the N3LO corrections. Grid-refinement and iteration-convergence tests are needed to show the numerical error is below the target accuracy.
- [Sec. 5.4.3, Sec. 5.5] The implementation of heavy-quark threshold matching inside the solver is not described. The AlphaS class handles n_f changes and libome provides OMEs, but the paper does not explain how the OMEs are applied at each threshold within the recursion framework. Given that v1 had a bug in the matching procedure (Sec. 5.5), the N3LO charm and bottom results in Fig. 2 cannot be assessed. A concise description of the matching step and a threshold-crossing check (e.g., PDF continuity) would address this.
minor comments (6)
- [Eq. (22)] The text 's=n+m+ell+k' uses an undefined k (and ell is not used in the sums). This appears to be a typographical error.
- [Eq. (29)] The notation R(2), R(1), R(0) should be R2, R1, R0 for consistency with Eq. (25).
- [Sec. 4] The citation 'GMVFN schemes currently in use [104? ? ?–113]' contains placeholder question marks and is broken.
- [Abstract / Sec. 6] The abstract states 'N3LO accuracy in QCD', while the conclusions and Sec. 3 correctly use 'approximate N3LO' (aN3LO). Please make the wording consistent, since the singlet splitting functions are approximations.
- [Listing 1] The variables order, iterations, and trunc_idx are declared as const double but represent integer counts; using an unsigned integer type would be clearer.
- [Sec. 5.2] The description of the 'tests' as 'not standard tests' is confusing; consider renaming them to 'benchmark and plot-generation programs'.
Circularity Check
No significant circularity: the code's outputs are computed from externally specified splitting functions, OMEs, and benchmark input PDFs via recursion relations that are stated in the paper and taken as inputs, not fitted to the outputs.
full rationale
The derivation chain in the paper is not circular. The DGLAP equation (1) is combined with the exponential power-series ansatz (17), whose L-functions (18)-(21) are obtained from the algebraic Mellin-space solution, and the coefficient recursions (23a)-(23d) and (31) are stated explicitly and implemented. These recursions are attributed to the authors' earlier Ref. [6], a self-citation, but they are re-used inputs rather than the target result: the paper's claims are about the code reproducing DGLAP evolution, and the recursions are the method being implemented. The actual numerical outputs (Table 2, Figs. 1-3) are generated from external theory inputs — the exact four-loop non-singlet splitting functions [11], state-of-the-art singlet approximations [13-20], three-loop OMEs through libome [21], and the Les Houches toy-model initial conditions [101] — with no parameter fitted to the evolved PDFs. The beta-function roots in Table 1 are computed numerically, not fitted. The paper even explicitly states that the truncated recursion relations 'generate truncated solutions at arbitrary order κ that are equivalent those obtained using the U-matrix method [7, 56, 64],' acknowledging rather than concealing the algorithm's relation to known methods. The possible weaknesses — fixed Gauss-Legendre grid defaults (Sec. 5.4.2), no convergence study, 'tests' described as 'not standard tests' (Sec. 5.2), and the disclosed v1 bugs (Sec. 5.5) — are numerical validation and correctness risks, not circular reasoning. A self-citation to the prior derivation of the algorithm is a normal, legitimate reference for a software paper and does not by itself make the derivation circular.
Assumptions & free parameters
free parameters (2)
- series truncation settings (iterations s, trunc_idx κ) =
10 (default in examples/Listing 1)
- grid/quadrature defaults =
100/50/25 grid nodes; 50 Gauss-Legendre points; 4 interpolation points
assumptions (6)
- domain assumption Collinear factorization: PDFs evolve with scale according to the DGLAP equations, Eq. (1).
- domain assumption The exponential logarithmic ansatz, Eq. (17), with L0-L3 from the roots of the 4-loop β-function (Table 1), spans the exact DGLAP solution in the NS sector.
- domain assumption Envelope approximations (FHM-RUVV, Refs. [13-20]) of the 4-loop singlet splitting functions represent the true singlet P(3) accurately enough for aN3LO evolution.
- domain assumption 3-loop OMEs from libome [21] provide the correct VFNS matching conditions at N3LO.
- domain assumption Les Houches toy-model initial conditions [101] and pole masses mc=µ0, mb=4.5 GeV, mt=175 GeV define the benchmark.
- standard math β-function coefficients β0-β3 are known and correctly implemented (Eqs. (5)-(8)).
Cite this review
Pith. "Pith review of Candia-v2: Logarithmic expansions for DGLAP evolution in $x$-space." pith.science (2026). https://pith.science/paper/FDY5KSOJ
@misc{pith2026260714151,
author = {Pith},
title = {Pith review of: Candia-v2: Logarithmic expansions for DGLAP evolution in $x$-space},
year = {2026},
howpublished = {\url{https://pith.science/paper/FDY5KSOJ}},
note = {Machine review of arXiv:2607.14151}
}
abstract
We present Candia-v2, an open-source software package that generalizes the $x$-space Candia algorithm to next-to-next-to-next-to-leading order (N$^{3}$LO) accuracy in Quantum Chromodynamics (QCD). The code solves the DGLAP evolution equations for unpolarized nucleon parton densities using a highly efficient logarithmic expansion technique that can be systematically extended to all orders in QCD. Candia-v2 supersedes the previous original C implementation of the algorithm with significantly increased efficiency and improved API. The software is publicly available on GitHub under the GPLv3 license.
Figures
Reference graph
Works this paper leans on
-
[6]
C. Hampson, M. Guzzi, DGLAP evolution at N3LO with the candia algorithm, Phys. Rev. D 113 (7) (2026) 074008.arXiv:2512.22667, doi:10.1103/p6vv-8zph
arXiv 2026
-
[101]
M. Dittmar, et al., Working Group I: Parton distributions: Summary report for the HERA LHC Workshop Proceedings (11 2005).arXiv: hep-ph/0511119
arXiv 2005
-
[75]
A. Cooper-Sarkar, T. Cridge, F. Giuli, L. A. Harland-Lang, F. Hekhorn, J. Huston, G. Magni, S. Moch, R. S. Thorne, A Bench- marking of QCD Evolution at ApproximateN3LO(6 2024).arXiv: 2406.16188
arXiv 2024
-
[1]
V. N. Gribov, L. N. Lipatov, e+ e- pair annihilation and deep inelastic e p scattering in perturbation theory, Sov. J. Nucl. Phys. 15 (1972) 675–684
1972
-
[2]
V. N. Gribov, L. N. Lipatov, Deep inelastic e p scattering in perturba- tion theory, Sov. J. Nucl. Phys. 15 (1972) 438–450
1972
-
[3]
L. N. Lipatov, The parton model and perturbation theory, Yad. Fiz. 20 (1974) 181–198
1974
-
[4]
Y. L. Dokshitzer, Calculation of the Structure Functions for Deep In- elastic Scattering and e+ e- Annihilation by Perturbation Theory in Quantum Chromodynamics., Sov. Phys. JETP 46 (1977) 641–653
1977
-
[5]
G. Altarelli, G. Parisi, Asymptotic Freedom in Parton Language, Nucl. Phys. B 126 (1977) 298–318.doi:10.1016/0550-3213(77)90384-4
Show all 113 references
-
[7]
Cafarella, C
A. Cafarella, C. Coriano, M. Guzzi, Nnlo logarithmic expansions and exact solutions of the DGLAP equations from x-space: New algorithms for precision studies at the lhc, Nucl. Phys. B 748 (2006) 253–308. arXiv:hep-ph/0512358,doi:10.1016/j.nuclphysb.2006.05.010
2006 arXiv
-
[8]
Guzzi, QCD Studies at Hadron Colliders and in Deeply Virtual Neutrino Scattering, Other thesis (12 2006).arXiv:hep-ph/0612355
M. Guzzi, QCD Studies at Hadron Colliders and in Deeply Virtual Neutrino Scattering, Other thesis (12 2006).arXiv:hep-ph/0612355
2006 arXiv
-
[9]
Cafarella, C
A. Cafarella, C. Coriano, M. Guzzi, Precision Studies of the NNLO DGLAP Evolution at the LHC with CANDIA, Comput. Phys. Com- mun. 179 (2008) 665–684.arXiv:0803.0462,doi:10.1016/j.cpc. 2008.06.004
2008 arXiv
-
[10]
GitHub,Candia-v2,https://github.com/champso1/candia-v2 (2025)
2025
-
[11]
Gehrmann, A
T. Gehrmann, A. von Manteuffel, V. Sotnikov, T.-Z. Yang, The four- loop non-singlet splitting functions in QCD (4 2026).arXiv:2604. 09534. 36
2026
-
[12]
S. Moch, A. Vogt, Properties and implications of the four-loop non- singlet splitting functions in QCD (5 2026).arXiv:2605.03889
2026 arXiv
-
[13]
Davies, C
J. Davies, C. H. Kom, S. Moch, A. Vogt, Resummation of small-x double logarithms in QCD: inclusive deep-inelastic scattering, JHEP 08 (2022) 135.arXiv:2202.10362,doi:10.1007/JHEP08(2022)135
2022 arXiv
-
[14]
S. Moch, B. Ruijl, T. Ueda, J. A. M. Vermaseren, A. Vogt, Low mo- ments of the four-loop splitting functions in QCD, Phys. Lett. B 825 (2022) 136853.arXiv:2111.15561,doi:10.1016/j.physletb.2021. 136853
2022 arXiv
-
[15]
Falcioni, F
G. Falcioni, F. Herzog, S. Moch, A. Vogt, Four-loop splitting functions in QCD – The quark-quark case, Phys. Lett. B 842 (2023) 137944. arXiv:2302.07593,doi:10.1016/j.physletb.2023.137944
2023 arXiv
-
[16]
Falcioni, F
G. Falcioni, F. Herzog, S. Moch, A. Vogt, Four-loop splitting functions in QCD – The gluon-to-quark case, Phys. Lett. B 846 (2023) 138215. arXiv:2307.04158,doi:10.1016/j.physletb.2023.138215
2023 arXiv
-
[17]
S. Moch, B. Ruijl, T. Ueda, J. Vermaseren, A. Vogt, Additional mo- ments and x-space approximations of four-loop splitting functions in QCD, Phys. Lett. B 849 (2024) 138468.arXiv:2310.05744,doi: 10.1016/j.physletb.2024.138468
2024 arXiv
-
[18]
G.Falcioni, F.Herzog, S.Moch, A.Pelloni, A.Vogt, Four-loopsplitting functions in QCD – The quark-to-gluon case, Phys. Lett. B 856 (2024) 138906.arXiv:2404.09701,doi:10.1016/j.physletb.2024.138906
2024 arXiv
-
[19]
G.Falcioni, F.Herzog, S.Moch, A.Pelloni, A.Vogt, Four-loopsplitting functions in QCD – the gluon-gluon case –, Phys. Lett. B 860 (2025) 139194.arXiv:2410.08089,doi:10.1016/j.physletb.2024.139194
2025 arXiv
-
[20]
Falcioni, F
G. Falcioni, F. Herzog, S. Moch, A. Pelloni, A. Vogt, Additional results on the four-loop flavour-singlet splitting functions in QCD, Phys. Lett. B 875 (2026) 140278.arXiv:2512.10783,doi:10.1016/j.physletb. 2026.140278
2026
-
[21]
Ablinger, A
J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel, C. Schneider, K. Schönwald, The single-mass variable flavor number 37 scheme at three-loop order, JHEP 03 (2026) 248.arXiv:2510.02175, doi:10.1007/JHEP03(2026)248
2026
-
[22]
Ablinger, J
J. Ablinger, J. Blümlein, A. De Freitas, A. von Manteuffel, C. Schnei- der, K. Schönwald, The two-mass contributions to the three-loop mas- sive operator matrix elements eA(3) Qg and∆ eA (3) Qg, JHEP 01 (2026) 111. arXiv:2510.09403,doi:10.1007/JHEP01(2026)111
2026
-
[23]
Ablinger, A
J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel, C. Schneider, K. Schoenwald, The three-loop single-mass heavy flavor corrections to deep-inelastic scattering (7 2024).arXiv:2407.02006
2024 arXiv
-
[24]
Ablinger, A
J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Man- teuffel, C. Schneider, K. Schönwald, The non-first-order-factorizable contributions to the three-loop single-mass operator matrix elements A(3) Qg and∆A (3) Qg, Phys. Lett. B 854 (2024) 138713.arXiv:2403.00513, d...
2024 arXiv
-
[25]
Ablinger, A
J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel, C. Schneider, K. Schönwald, The first–order factorizable contributions to the three–loop massive operator matrix elementsA(3) Qg and∆A (3) Qg, Nucl. Phys. B 999 (2024) 116427.arXiv:2311.00644,doi:10.1016/ ...
2024 arXiv
-
[26]
Bierenbaum, J
I. Bierenbaum, J. Blümlein, A. De Freitas, A. Goedicke, S. Klein, K. Schönwald,O(α 2 s) polarized heavy flavor corrections to deep- inelastic scattering atQ 2 ≫m 2, Nucl. Phys. B 988 (2023) 116114. arXiv:2211.15337,doi:10.1016/j.nuclphysb.2023.116114
2023 arXiv
-
[27]
Ablinger, A
J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. Goedicke, A. von Manteuffel, C. Schneider, K. Schönwald, The unpolarized and polarized single-mass three-loop heavy flavor operator matrix el- ements Agg,Q and∆A gg,Q, JHEP 12 (2022) 134.arXiv:2211.05462, doi:10.1007/JHEP...
2022 arXiv
-
[28]
Behring, J
A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel, K. Schön- wald, C. Schneider, The polarized transition matrix elementAgq(N) of the variable flavor number scheme atO(α 3 s), Nucl. Phys. B 964 (2021) 115331.arXiv:2101.05733,doi:10.1016/j.nuclphysb. 2021.115331. 38
2021 arXiv
-
[29]
Ablinger, J
J. Ablinger, J. Blümlein, A. De Freitas, A. Goedicke, M. Saragnese, C. Schneider, K. Schönwald, The two-mass contribution to the three- loop polarized gluonic operator matrix elementA(3) gg,Q, Nucl. Phys. B 955 (2020) 115059.arXiv:2004.08916,doi:10.1016/j.nuclphysb. 2020.115059
2020 arXiv
-
[30]
Ablinger, A
J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Manteuf- fel, C. Schneider, K. Schönwald, The three-loop single mass polarized pure singlet operator matrix element, Nucl. Phys. B 953 (2020) 114945. arXiv:1912.02536,doi:10.1016/j.nuclphysb.2020.114945
2020 arXiv
-
[31]
Ablinger, J
J. Ablinger, J. Blümlein, A. De Freitas, A. Goedicke, C. Schneider, K. Schönwald, The Two-mass Contribution to the Three-Loop Gluonic Operator Matrix ElementA (3) gg,Q, Nucl. Phys. B 932 (2018) 129–240. arXiv:1804.02226,doi:10.1016/j.nuclphysb.2018.04.023
2018 arXiv
-
[32]
Ablinger, J
J. Ablinger, J. Blümlein, A. De Freitas, C. Schneider, K. Schönwald, The two-mass contribution to the three-loop pure singlet operator ma- trix element, Nucl. Phys. B 927 (2018) 339–367.arXiv:1711.06717, doi:10.1016/j.nuclphysb.2017.12.018
2018 arXiv
-
[33]
Ablinger, J
J. Ablinger, J. Blümlein, A. De Freitas, A. Hasselhuhn, C. Schnei- der, F. Wißbrock, Three Loop Massive Operator Matrix Elements and Asymptotic Wilson Coefficients with Two Different Masses, Nucl. Phys. B 921 (2017) 585–688.arXiv:1705.07030,doi:10.1016/j. nuclphysb.2017.05.017
2017 arXiv
-
[34]
Ablinger, J
J. Ablinger, J. Blümlein, A. De Freitas, A. Hasselhuhn, A. von Man- teuffel, M. Round, C. Schneider, 3-loop MassiveO(T2 F )Contributions to the DIS Operator Matrix ElementAgg, Nucl. Part. Phys. Proc. 258- 259 (2015) 37–40.arXiv:1409.1435,doi:10.1016/j.nuclphysbps. 2015.01.009
2015 arXiv
-
[35]
Ablinger, J
J. Ablinger, J. Blümlein, A. De Freitas, A. Hasselhuhn, A. von Man- teuffel, M. Round, C. Schneider, TheO(α 3 sT 2 F )Contributions to the Gluonic Operator Matrix Element, Nucl. Phys. B 885 (2014) 280–317. arXiv:1405.4259,doi:10.1016/j.nuclphysb.2014.05.028
2014 arXiv
-
[36]
Behring, I
A. Behring, I. Bierenbaum, J. Blümlein, A. De Freitas, S. Klein, F. Wißbrock, The logarithmic contributions to theO(α3 s)asymptotic 39 massive Wilson coefficients and operator matrix elements in deeply inelastic scattering, Eur. Phys. J. C 74 (9) (2014) 3033.arXiv: 1403.6356,d...
2014 arXiv
-
[37]
Ablinger, J
J. Ablinger, J. Blümlein, A. De Freitas, A. Hasselhuhn, A. von Man- teuffel, M. Round, C. Schneider, F. Wißbrock, The transition ma- trix elementA gq(N)of the variable flavor number scheme atO(α 3 s), Nucl. Phys. B 882 (2014) 263.arXiv:1402.0359,doi:10.1016/j. nuclphysb.2014.02.007
2014 arXiv
-
[38]
Ablinger, J
J. Ablinger, J. Blumlein, A. Hasselhuhn, S. Klein, C. Schneider, F. Wissbrock, Massive 3-loop Ladder Diagrams for Quarkonic Local Operator Matrix Elements, Nucl. Phys. B 864 (2012) 52–84.arXiv: 1206.2252,doi:10.1016/j.nuclphysb.2012.06.007
2012 arXiv
-
[39]
Blumlein, A
J. Blumlein, A. Hasselhuhn, S. Klein, C. Schneider, The O(α3 snf T 2 F CA,F )Contributions to the Gluonic Massive Operator Ma- trix Elements, Nucl. Phys. B 866 (2013) 196–211.arXiv:1205.4184, doi:10.1016/j.nuclphysb.2012.09.001
2013 arXiv
-
[40]
Blumlein, A
J. Blumlein, A. De Freitas, W. van Neerven, Two-loop QED Operator Matrix Elements with Massive External Fermion Lines, Nucl. Phys. B 855 (2012) 508–569.arXiv:1107.4638,doi:10.1016/j.nuclphysb. 2011.10.009
2012 arXiv
-
[41]
Ablinger, J
J. Ablinger, J. Blumlein, S. Klein, C. Schneider, F. Wissbrock, The O(α3 s)Massive Operator Matrix Elements ofO(n f )for the Structure FunctionF 2(x, Q2)and Transversity, Nucl. Phys. B 844 (2011) 26–54. arXiv:1008.3347,doi:10.1016/j.nuclphysb.2010.10.021
2011 arXiv
-
[42]
Bierenbaum, J
I. Bierenbaum, J. Blumlein, S. Klein, The Gluonic Operator Matrix Elements atO(α 2 s)for DIS Heavy Flavor Production, Phys. Lett. B 672 (2009) 401–406.arXiv:0901.0669,doi:10.1016/j.physletb. 2009.01.057
2009 arXiv
-
[43]
Bierenbaum, J
I. Bierenbaum, J. Blumlein, S. Klein, C. Schneider, Two-Loop Massive Operator Matrix Elements for Unpolarized Heavy Flavor Production toO(ε), Nucl. Phys. B 803 (2008) 1–41.arXiv:0803.0273,doi:10. 1016/j.nuclphysb.2008.05.016. 40
2008 arXiv
-
[44]
Bierenbaum, J
I. Bierenbaum, J. Blumlein, S. Klein, Two-Loop Massive Operator Ma- trix Elements and Unpolarized Heavy Flavor Production at Asymp- totic ValuesQ 2 ≫m 2, Nucl. Phys. B 780 (2007) 40–75.arXiv: hep-ph/0703285,doi:10.1016/j.nuclphysb.2007.04.030
2007 arXiv
-
[45]
Bierenbaum, J
I. Bierenbaum, J. Blumlein, S. Klein, Calculation of massive 2-loop operator matrix elements with outer gluon lines, Phys. Lett. B 648 (2007) 195–200.arXiv:hep-ph/0702265,doi:10.1016/j.physletb. 2007.02.059
2007 arXiv
-
[46]
Abdul Khalek, et al., Science Requirements and Detector Con- cepts for the Electron-Ion Collider: EIC Yellow Report, Nucl
R. Abdul Khalek, et al., Science Requirements and Detector Con- cepts for the Electron-Ion Collider: EIC Yellow Report, Nucl. Phys. A 1026 (2022) 122447.arXiv:2103.05419,doi:10.1016/j.nuclphysa. 2022.122447
2022 arXiv
-
[47]
Furmanski, R
W. Furmanski, R. Petronzio, Lepton - Hadron Processes Beyond Lead- ing Order in Quantum Chromodynamics, Z. Phys. C 11 (1982) 293. doi:10.1007/BF01578280
1982 doi
-
[48]
Rossi, X SPACE ANALYSIS FOR THE PHOTON STRUCTURE FUNCTIONS IN QCD, Phys
G. Rossi, X SPACE ANALYSIS FOR THE PHOTON STRUCTURE FUNCTIONS IN QCD, Phys. Rev. D 29 (1984) 852.doi:10.1103/ PhysRevD.29.852
1984
-
[49]
J. H. Da Luz Vieira, J. K. Storrow, A Phenomenological analysis of the photon structure function F2 (x, Q**2), Z. Phys. C 51 (1991) 241–258. doi:10.1007/BF01475792
1991 doi
-
[50]
Kumano, J
S. Kumano, J. T. Londergan, A FORTRAN program for numerical so- lution of the Altarelli-Parisi equations by the Laguerre method, Com- put.Phys.Commun.69(1992)373–396.doi:10.1016/0010-4655(92) 90176-Y
1992 doi
-
[51]
Pascaud, F
C. Pascaud, F. Zomer, Fast method to integrate linear integrodifferen- tial equations (6 1994)
1994
-
[52]
Miyama, S
M. Miyama, S. Kumano, Numerical solution of Q**2 evolution equa- tions in a brute force method, Comput. Phys. Commun. 94 (1996) 185– 215.arXiv:hep-ph/9508246,doi:10.1016/0010-4655(96)00013-6
1996 arXiv
-
[53]
Pascaud, F
C. Pascaud, F. Zomer, A Fast and precise method to solve the Altarelli- Parisi equations in x space (12 1996). 41
1996
-
[54]
D. A. Kosower, Evolution of parton distributions, Nucl. Phys. B 506 (1997) 439–467.arXiv:hep-ph/9706213,doi:10.1016/ S0550-3213(97)00526-9
1997 arXiv
-
[55]
Hirai, S
M. Hirai, S. Kumano, M. Miyama, Numerical solution of Q**2 evo- lution equations for polarized structure functions, Comput. Phys. Commun. 108 (1998) 38.arXiv:hep-ph/9707220,doi:10.1016/ S0010-4655(97)00129-X
1998 arXiv
-
[56]
Blumlein, A
J. Blumlein, A. Vogt, The Evolution of unpolarized singlet structure functions at small x, Phys. Rev. D 58 (1998) 014020.arXiv:hep-ph/ 9712546,doi:10.1103/PhysRevD.58.014020
1998 doi
-
[57]
Coriano, C
C. Coriano, C. Savkli, QCD evolution equations: Numerical algorithms from the Laguerre expansion, Comput. Phys. Commun. 118 (1999) 236–258.arXiv:hep-ph/9803336,doi:10.1016/S0010-4655(98) 00158-1
1999 arXiv
-
[58]
Santorelli, E
P. Santorelli, E. Scrimieri, A Semianalytical method to evolve par- ton distributions, Phys. Lett. B 459 (1999) 599–606.arXiv:hep-ph/ 9807572,doi:10.1016/S0370-2693(99)00698-X
1999 doi
-
[59]
P. G. Ratcliffe, A matrix approach to numerical solution of the DGLAP evolution equations, Phys. Rev. D 63 (2001) 116004.arXiv:hep-ph/ 0012376,doi:10.1103/PhysRevD.63.116004
2001 doi
-
[60]
Weinzierl, Fast evolution of parton distributions, Comput
S. Weinzierl, Fast evolution of parton distributions, Comput. Phys. Commun.148(2002)314–326.arXiv:hep-ph/0203112,doi:10.1016/ S0010-4655(02)00584-2
2002 arXiv
-
[61]
Jadach, M
S. Jadach, M. Skrzypek, Solving constrained Markovian evolution in QCD with the help of the non-Markovian Monte Carlo, Comput. Phys. Commun.175(2006)511–527.arXiv:hep-ph/0504263,doi:10.1016/ j.cpc.2006.06.004
2006 arXiv
-
[62]
Jadach, W
S. Jadach, W. Placzek, M. Skrzypek, P. Stephens, Z. Was, Constrained MCforQCDevolutionwithrapidityorderingandminimumkT*, Com- put. Phys. Commun. 180 (2009) 675–698.arXiv:hep-ph/0703281, doi:10.1016/j.cpc.2008.11.014. 42
2009 arXiv
-
[63]
Simonelli, Analytic solutions of the DGLAP evolution and the- oretical uncertainties, Eur
A. Simonelli, Analytic solutions of the DGLAP evolution and the- oretical uncertainties, Eur. Phys. J. C 84 (8) (2024) 867.arXiv: 2401.13663,doi:10.1140/epjc/s10052-024-13169-8
2024 arXiv
-
[64]
Vogt, Efficient evolution of unpolarized and polarized parton dis- tributions with QCD-PEGASUS, Comput
A. Vogt, Efficient evolution of unpolarized and polarized parton dis- tributions with QCD-PEGASUS, Comput. Phys. Commun. 170 (2005) 65–92.arXiv:hep-ph/0408244,doi:10.1016/j.cpc.2005.03.103
2005 arXiv
-
[65]
GitHub,Candia,https://github.com/champso1/candia(2008)
2008
-
[66]
G. P. Salam, J. Rojo, A Higher Order Perturbative Parton Evolu- tionToolkit(HOPPET),Comput.Phys.Commun.180(2009)120–156. arXiv:0804.3755,doi:10.1016/j.cpc.2008.08.010
2009 arXiv
-
[67]
Karlberg, P
A. Karlberg, P. Nason, G. Salam, G. Zanderighi, F. Dreyer, HOPPET v2 release note, Eur. Phys. J. C 86 (2) (2026) 157.doi:10.1140/epjc/ s10052-025-15200-y
2026 doi
-
[68]
Botje, QCDNUM: Fast QCD Evolution and Convolution, Comput
M. Botje, QCDNUM: Fast QCD Evolution and Convolution, Comput. Phys.Commun.182(2011)490–532.arXiv:1005.1481,doi:10.1016/ j.cpc.2010.10.020
2011 arXiv
-
[69]
Bertone, S
V. Bertone, S. Carrazza, J. Rojo, APFEL: A PDF Evolution Library with QED corrections, Comput. Phys. Commun. 185 (2014) 1647–1668. arXiv:1310.1394,doi:10.1016/j.cpc.2014.03.007
2014 arXiv
-
[70]
Bertone, APFEL++: A new PDF evolution library in C++, PoS DIS2017 (2018) 201.arXiv:1708.00911,doi:10.22323/1.297.0201
V. Bertone, APFEL++: A new PDF evolution library in C++, PoS DIS2017 (2018) 201.arXiv:1708.00911,doi:10.22323/1.297.0201
2018 arXiv
-
[71]
Diehl, R
M. Diehl, R. Nagar, F. J. Tackmann, ChiliPDF: Chebyshev interpo- lation for parton distributions, Eur. Phys. J. C 82 (3) (2022) 257. arXiv:2112.09703,doi:10.1140/epjc/s10052-022-10223-1
2022 arXiv
-
[72]
Candido, F
A. Candido, F. Hekhorn, G. Magni, EKO: evolution kernel operators, Eur. Phys. J. C 82 (10) (2022) 976.arXiv:2202.02338,doi:10.1140/ epjc/s10052-022-10878-w
2022 arXiv
-
[73]
Hautmann, H
F. Hautmann, H. Jung, S. T. Monfared, The CCFM uPDF evolution uPDFevolv Version 1.0.00, Eur. Phys. J. C 74 (2014) 3082.arXiv: 1407.5935,doi:10.1140/epjc/s10052-014-3082-1. 43
2014 arXiv
-
[74]
H. Jung, A. Lelek, K. M. Figueroa, S. Taheri Monfared, The Par- ton Branching evolution package uPDFevolv2 (5 2024).arXiv:2405. 20185
2024
-
[76]
R. D. Ball, et al., The path to N3LO parton distributions, Eur. Phys. J. C 84 (7) (2024) 659.arXiv:2402.18635,doi:10.1140/epjc/ s10052-024-12891-7
2024 arXiv
-
[77]
McGowan, T
J. McGowan, T. Cridge, L. A. Harland-Lang, R. S. Thorne, Ap- proximate N3LO parton distribution functions with theoretical un- certainties: MSHT20aN 3LO PDFs, Eur. Phys. J. C 83 (3) (2023) 185, [Erratum: Eur.Phys.J.C 83, 302 (2023)].arXiv:2207.04739, doi:10.1140/epjc/s10052-02...
2023 arXiv
-
[78]
Cridge, L
T. Cridge, L. A. Harland-Lang, R. S. Thorne, Combining QED and approximateN 3LO QCD corrections in a global PDF fit: MSHT20qed_an3lo PDFs, SciPost Phys. 17 (1) (2024) 026.arXiv: 2312.07665,doi:10.21468/SciPostPhys.17.1.026
2024 arXiv
-
[79]
R. D. Ball, A. Barontini, J. Cruz-Martinez, S. Forte, F. Hekhorn, E. R. Nocera, J. Rojo, R. Stegeman, A determination ofαs(mZ)at aN3LOQCD ⊗NLO QED accuracy from a global PDF analysis, Eur. Phys. J. C 85 (9) (2025) 1001.arXiv:2506.13871,doi:10.1140/epjc/ s10052-025-14676-y
2025 arXiv
-
[80]
D. J. Gross, F. Wilczek, Ultraviolet Behavior of Nonabelian Gauge Theories, Phys. Rev. Lett. 30 (1973) 1343–1346.doi:10.1103/ PhysRevLett.30.1343
1973
-
[81]
H. D. Politzer, Asymptotic Freedom: An Approach to Strong Interac- tions, Phys. Rept. 14 (1974) 129–180.doi:10.1016/0370-1573(74) 90014-3
1974 doi
-
[82]
E. G. Floratos, D. A. Ross, C. T. Sachrajda, Higher Order Effects in Asymptotically Free Gauge Theories: The Anomalous Dimen- 44 sions of Wilson Operators, Nucl. Phys. B 129 (1977) 66–88, [Erra- tum: Nucl.Phys.B139, 545–546(1978)].doi:10.1016/0550-3213(77) 90020-7
1977 doi
-
[83]
E. G. Floratos, D. A. Ross, C. T. Sachrajda, Higher Order Effects in Asymptotically Free Gauge Theories. 2. Flavor Singlet Wilson Oper- ators and Coefficient Functions, Nucl. Phys. B 152 (1979) 493–520. doi:10.1016/0550-3213(79)90094-4
1979 doi
-
[84]
Gonzalez-Arroyo, C
A. Gonzalez-Arroyo, C. Lopez, F. J. Yndurain, Second Order Con- tributions to the Structure Functions in Deep Inelastic Scattering. 1. Theoretical Calculations, Nucl. Phys. B 153 (1979) 161–186.doi: 10.1016/0550-3213(79)90596-0
1979 doi
-
[85]
Gonzalez-Arroyo, C
A. Gonzalez-Arroyo, C. Lopez, Second Order Contributions to the Structure Functions in Deep Inelastic Scattering. 3. The Singlet Case, Nucl. Phys. B 166 (1980) 429–459.doi:10.1016/0550-3213(80) 90207-2
1980 doi
-
[86]
Curci, W
G. Curci, W. Furmanski, R. Petronzio, Evolution of Parton Densities Beyond Leading Order: The Nonsinglet Case, Nucl. Phys. B 175 (1980) 27–92.doi:10.1016/0550-3213(80)90003-6
1980 doi
-
[87]
Furmanski, R
W. Furmanski, R. Petronzio, Singlet Parton Densities Beyond Leading Order, Phys. Lett. B 97 (1980) 437–442.doi:10.1016/0370-2693(80) 90636-X
1980 doi
-
[88]
E. G. Floratos, C. Kounnas, R. Lacaze, Higher Order QCD Effects in Inclusive Annihilation and Deep Inelastic Scattering, Nucl. Phys. B 192 (1981) 417–462.doi:10.1016/0550-3213(81)90434-X
1981 doi
-
[89]
Hamberg, W
R. Hamberg, W. L. van Neerven, The Correct renormalization of the gluon operator in a covariant gauge, Nucl. Phys. B 379 (1992) 143–171. doi:10.1016/0550-3213(92)90593-Z
1992 doi
-
[90]
S. Moch, J. A. M. Vermaseren, A. Vogt, The Three loop splitting func- tions in QCD: The Nonsinglet case, Nucl. Phys. B 688 (2004) 101–134. arXiv:hep-ph/0403192,doi:10.1016/j.nuclphysb.2004.03.030. 45
2004 arXiv
-
[91]
A. Vogt, S. Moch, J. A. M. Vermaseren, The Three-loop splitting func- tions in QCD: The Singlet case, Nucl. Phys. B 691 (2004) 129–181. arXiv:hep-ph/0404111,doi:10.1016/j.nuclphysb.2004.04.024
2004 arXiv
-
[92]
J. A. Gracey, Anomalous dimension of nonsinglet Wilson operators at O (1 / N(f)) in deep inelastic scattering, Phys. Lett. B 322 (1994) 141– 146.arXiv:hep-ph/9401214,doi:10.1016/0370-2693(94)90502-9
1994 arXiv
-
[93]
Davies, A
J. Davies, A. Vogt, B. Ruijl, T. Ueda, J. A. M. Vermaseren, Large-nf contributionstothefour-loopsplittingfunctionsinQCD,Nucl.Phys.B 915 (2017) 335–362.arXiv:1610.07477,doi:10.1016/j.nuclphysb. 2016.12.012
2017 arXiv
-
[94]
S. Moch, B. Ruijl, T. Ueda, J. A. M. Vermaseren, A. Vogt, Four-Loop Non-Singlet Splitting Functions in the Planar Limit and Beyond, JHEP 10 (2017) 041.arXiv:1707.08315,doi:10.1007/JHEP10(2017)041
2017 arXiv
-
[95]
Gehrmann, A
T. Gehrmann, A. von Manteuffel, V. Sotnikov, T.-Z. Yang, Complete N 2 f contributions to four-loop pure-singlet splitting functions, JHEP 01 (2024) 029.arXiv:2308.07958,doi:10.1007/JHEP01(2024)029
2024 arXiv
-
[96]
Falcioni, F
G. Falcioni, F. Herzog, S. Moch, J. Vermaseren, A. Vogt, The dou- ble fermionic contribution to the four-loop quark-to-gluon splitting function, Phys. Lett. B 848 (2024) 138351.arXiv:2310.01245,doi: 10.1016/j.physletb.2023.138351
2024 arXiv
-
[97]
Gehrmann, A
T. Gehrmann, A. von Manteuffel, V. Sotnikov, T.-Z. Yang, The NfCF3 contribution to the non-singlet splitting function at four-loop order, Phys. Lett. B 849 (2024) 138427.arXiv:2310.12240,doi:10.1016/ j.physletb.2023.138427
2024 arXiv
-
[98]
B. A. Kniehl, S. Moch, V. N. Velizhanin, A. Vogt, Flavor Nonsinglet Splitting Functions at Four Loops in QCD: Fermionic Contributions, Phys. Rev. Lett. 135 (7) (2025) 071902.arXiv:2505.09381,doi:10. 1103/hkg5-88hr
2025 arXiv
-
[99]
van Ritbergen, J
T. van Ritbergen, J. A. M. Vermaseren, S. A. Larin, The Four loop beta function in quantum chromodynamics, Phys. Lett. B 400 (1997) 379–384.arXiv:hep-ph/9701390,doi:10.1016/S0370-2693(97) 00370-5. 46
1997 arXiv
-
[100]
Czakon, The Four-loop QCD beta-function and anomalous dimen- sions, Nucl
M. Czakon, The Four-loop QCD beta-function and anomalous dimen- sions, Nucl. Phys. B 710 (2005) 485–498.arXiv:hep-ph/0411261, doi:10.1016/j.nuclphysb.2005.01.012
2005 arXiv
-
[102]
GitHub,libome,https://gitlab.com/libome/libome(2025)
2025
-
[103]
M. Buza, Y. Matiounine, J. Smith, R. Migneron, W. L. van Neerven, Heavy quark coefficient functions at asymptotic valuesQ2 ≫m 2, Nucl. Phys. B 472 (1996) 611–658.arXiv:hep-ph/9601302,doi:10.1016/ 0550-3213(96)00228-3
1996 arXiv
-
[104]
M. Buza, Y. Matiounine, J. Smith, W. van Neerven, Charm elec- troproduction viewed in the variable flavor number scheme versus fixed order perturbation theory, Eur. Phys. J. C 1 (1998) 301–320. arXiv:hep-ph/9612398,doi:10.1007/BF01245820
1998 arXiv
-
[105]
M. A. G. Aivazis, F. I. Olness, W.-K. Tung, Leptoproduction of heavy quarks. 1. General formalism and kinematics of charged current and neutral current production processes, Phys. Rev. D 50 (1994) 3085– 3101.arXiv:hep-ph/9312318,doi:10.1103/PhysRevD.50.3085
1994 arXiv
-
[106]
M. A. G. Aivazis, J. C. Collins, F. I. Olness, W.-K. Tung, Leptopro- duction of heavy quarks. 2. A Unified QCD formulation of charged and neutral current processes from fixed target to collider energies, Phys. Rev. D 50 (1994) 3102–3118.arXiv:hep-ph/9312319,doi: 10.1103/PhysRe...
1994 arXiv
-
[107]
W.-K. Tung, S. Kretzer, C. Schmidt, Open heavy flavor production in QCD: Conceptual framework and implementation issues, J. Phys. G 28 (2002) 983–996.arXiv:hep-ph/0110247,doi:10.1088/0954-3899/ 28/5/321
2002 arXiv
-
[108]
Krämer, F
M. Krämer, F. I. Olness, D. E. Soper, Treatment of heavy quarks in deeply inelastic scattering, Phys. Rev. D 62 (2000) 096007.arXiv: hep-ph/0003035,doi:10.1103/PhysRevD.62.096007. 47
2000 arXiv
-
[109]
Guzzi, P
M. Guzzi, P. M. Nadolsky, H.-L. Lai, C. P. Yuan, General-Mass Treat- ment for Deep Inelastic Scattering at Two-Loop Accuracy, Phys. Rev. D 86 (2012) 053005.arXiv:1108.5112,doi:10.1103/PhysRevD.86. 053005
2012 arXiv
-
[110]
Guzzi, P
M. Guzzi, P. Nadolsky, L. Reina, D. Wackeroth, K. Xie, General mass variable flavor number scheme for Z boson production in association with a heavy quark at hadron colliders, Phys. Rev. D 110 (11) (2024) 114030.arXiv:2410.03876,doi:10.1103/PhysRevD.110.114030
2024 arXiv
-
[111]
R. S. Thorne, R. G. Roberts, An Ordered analysis of heavy flavor production in deep inelastic scattering, Phys. Rev. D 57 (1998) 6871– 6898.arXiv:hep-ph/9709442,doi:10.1103/PhysRevD.57.6871
1998 arXiv
-
[112]
R. S. Thorne, A Variable-flavor number scheme for NNLO, Phys. Rev. D 73 (2006) 054019.arXiv:hep-ph/0601245,doi:10.1103/ PhysRevD.73.054019
2006 arXiv
-
[113]
Cacciari, M
M. Cacciari, M. Greco, P. Nason, ThepT spectrum in heavy-flavour hadroproduction., JHEP 05 (1998) 007.arXiv:hep-ph/9803400,doi: 10.1088/1126-6708/1998/05/007. 48
1998 arXiv
Reviewed August 2, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.