Pith. sign in

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 →

arxiv 2607.14151 v2 pith:FDY5KSOJ submitted 2026-07-14 hep-ph hep-ex

classification hep-phhep-ex PACS 12.38.-t12.38.Bx
keywords DGLAPevolutionpartondistributionfunctionsN3LOQCDx-spacelogarithmicexpansionsrecursionrelationsheavy-quarkmatchingsubtractionPDFs
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Candia-v2 claims to solve the DGLAP integro-differential equations that govern the scale evolution of proton parton densities directly in Bjorken-x space, at approximate N3LO accuracy in QCD. The method writes the solution as an exponential-type logarithmic power series and determines the unknown coefficient functions through recursion relations, avoiding Mellin-space inversion. The release adds exact four-loop non-singlet splitting functions, state-of-the-art approximations for the four-loop singlet sector, and three-loop operator matrix elements for heavy-quark matching. If the implementation is correct, the package is a production-grade evolution tool for high-precision collider phenomenology and supplies heavy-quark subtraction and residual PDFs for general-mass variable-flavor-number schemes.

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.

Watch

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 extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

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)
  1. [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].
  2. [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.
  3. [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.
  4. [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)
  1. [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.
  2. [Eq. (29)] The notation R(2), R(1), R(0) should be R2, R1, R0 for consistency with Eq. (25).
  3. [Sec. 4] The citation 'GMVFN schemes currently in use [104? ? ?–113]' contains placeholder question marks and is broken.
  4. [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.
  5. [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.
  6. [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

0 steps flagged · score 0.0 of 10

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 2 free parameters · 6 assumptions · 0 invented entities

Everything the numerical answers depend on is imported: the DGLAP framework (Eq. 1), the recursion relations from the authors' own Ref. [6], the exact 4-loop non-singlet splitting functions of Ref. [11], the approximate singlet inputs of Refs. [13-20], the 3-loop OMEs of libome (Refs. [21-45]), the 4-loop β-function (Refs. [99,100]), and the Les Houches toy-model boundary conditions (Ref. [101]). No parameters are fitted to data; the package's hand-chosen numerical settings (grid sizes, Gauss order, series truncation) are the only adjustable inputs and their convergence is not quantified.

free parameters (2)
  • series truncation settings (iterations s, trunc_idx κ) = 10 (default in examples/Listing 1)
    Numerical output depends on the number of terms kept in the exact (s) and truncated (κ) ansätze; the choice 10 is not justified by a convergence study in the paper.
  • grid/quadrature defaults = 100/50/25 grid nodes; 50 Gauss-Legendre points; 4 interpolation points
    Hand-chosen defaults in GridFillerArgs/ConvIntArgs (§5.4.2); no grid-refinement or quadrature-error study is presented, so discretization error is uncontrolled.
assumptions (6)
  • domain assumption Collinear factorization: PDFs evolve with scale according to the DGLAP equations, Eq. (1).
    Standard QCD factorization assumption; the entire package solves 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.
    The ansatz and recursion relations (23a-d) are taken from the authors' Ref. [6]; this paper does not re-derive them. If the ansatz or recursions are incomplete at N3LO, the output is wrong.
  • 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.
    The exact singlet 4-loop calculation is in progress (Refs. [92-98]); the code uses approximate inputs by default (Sec. 1, Sec. 5.4.6).
  • domain assumption 3-loop OMEs from libome [21] provide the correct VFNS matching conditions at N3LO.
    Used for heavy-quark threshold matching (Sec. 3, Sec. 5.4.6); external library input treated as exact.
  • domain assumption Les Houches toy-model initial conditions [101] and pole masses mc=µ0, mb=4.5 GeV, mt=175 GeV define the benchmark.
    Inputs adopted from Ref. [101]; Table 2 results are specific to these inputs.
  • standard math β-function coefficients β0-β3 are known and correctly implemented (Eqs. (5)-(8)).
    From Refs. [99,100]; standard four-loop QCD result.

how reviews work

0 comments
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

Figures reproduced from arXiv: 2607.14151 by the authors.

Figure 1
Figure 1. Results for the Candia-v2 evolution of the gluon and the uv PDFs. 14 [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Results for the Candia-v2 evolution of the c (−) and b (−) heavy flavor combina￾tions. 15 [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 3
Figure 3. Subtraction PDFs for the charm and bottom quark evolved at [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

113 extracted references · 9 canonical work pages

  1. [6]

    Hampson, M

    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

  2. [101]

    Dittmar, et al., Working Group I: Parton distributions: Summary report for the HERA LHC Workshop Proceedings (11 2005).arXiv: hep-ph/0511119

    M. Dittmar, et al., Working Group I: Parton distributions: Summary report for the HERA LHC Workshop Proceedings (11 2005).arXiv: hep-ph/0511119

  3. [75]

    Cooper-Sarkar, T

    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

  4. [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

  5. [2]

    V. N. Gribov, L. N. Lipatov, Deep inelastic e p scattering in perturba- tion theory, Sov. J. Nucl. Phys. 15 (1972) 438–450

  6. [3]

    L. N. Lipatov, The parton model and perturbation theory, Yad. Fiz. 20 (1974) 181–198

  7. [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

  8. [5]

    Altarelli, G

    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
  1. [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

  2. [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

  3. [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

  4. [10]

    GitHub,Candia-v2,https://github.com/champso1/candia-v2 (2025)

  5. [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

  6. [12]

    S. Moch, A. Vogt, Properties and implications of the four-loop non- singlet splitting functions in QCD (5 2026).arXiv:2605.03889

  7. [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

  8. [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

  9. [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

  10. [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

  11. [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

  12. [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

  13. [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

  14. [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

  15. [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

  16. [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

  17. [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

  18. [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...

  19. [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/ ...

  20. [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

  21. [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...

  22. [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

  23. [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

  24. [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

  25. [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

  26. [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

  27. [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

  28. [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

  29. [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

  30. [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...

  31. [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

  32. [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

  33. [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

  34. [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

  35. [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

  36. [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

  37. [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

  38. [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

  39. [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

  40. [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

  41. [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

  42. [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

  43. [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

  44. [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

  45. [51]

    Pascaud, F

    C. Pascaud, F. Zomer, Fast method to integrate linear integrodifferen- tial equations (6 1994)

  46. [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

  47. [53]

    Pascaud, F

    C. Pascaud, F. Zomer, A Fast and precise method to solve the Altarelli- Parisi equations in x space (12 1996). 41

  48. [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

  49. [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

  50. [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

  51. [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

  52. [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

  53. [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

  54. [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

  55. [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

  56. [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

  57. [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

  58. [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

  59. [65]

    GitHub,Candia,https://github.com/champso1/candia(2008)

  60. [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

  61. [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

  62. [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

  63. [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

  64. [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

  65. [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

  66. [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

  67. [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

  68. [74]

    H. Jung, A. Lelek, K. M. Figueroa, S. Taheri Monfared, The Par- ton Branching evolution package uPDFevolv2 (5 2024).arXiv:2405. 20185

  69. [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

  70. [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...

  71. [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

  72. [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

  73. [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

  74. [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

  75. [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

  76. [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

  77. [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

  78. [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

  79. [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

  80. [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

  81. [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

  82. [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

  83. [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

  84. [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

  85. [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

  86. [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

  87. [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

  88. [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

  89. [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

  90. [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

  91. [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

  92. [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

  93. [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

  94. [102]

    GitHub,libome,https://gitlab.com/libome/libome(2025)

  95. [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

  96. [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

  97. [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

  98. [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...

  99. [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

  100. [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

  101. [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

  102. [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

  103. [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

  104. [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

  105. [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

Pith tools

Reviewed August 2, 2026 · model on record in the stance chip above.