Pith. sign in

REVIEW 4 minor 61 references

Three-loop variable flavor number scheme matches full asymptotic heavy-flavor DIS corrections at large Q^{2}.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-11 07:36 UTC pith:ZX22NSLX

load-bearing objection Clean three-loop completion of the VFNS with public codes; the RG equivalence holds and the only real limit is the usual asymptotic one.

arxiv 2607.05235 v1 pith:ZX22NSLX submitted 2026-07-06 hep-ph hep-exhep-th

The variable flavor number scheme to three-loop order

classification hep-ph hep-exhep-th PACS 12.38.Bx13.60.Hb14.65.Dw14.65.Fy
keywords variable flavor number schemethree-loop QCDheavy-quark PDFsoperator matrix elementsdeep-inelastic scatteringWilson coefficientscharm-bottom matchingrenormalization group
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper completes the variable flavor number scheme (VFNS) to three-loop order for both single-mass and two-mass heavy quarks. The scheme redefines ordinary light-parton densities by universal massive operator matrix elements and introduces charm and bottom parton distributions. A renormalization-group argument shows that, once power-suppressed mass terms are dropped, the resulting structure functions agree with the fixed-flavor asymptotic heavy-flavor corrections through O(α_s^{3}), the residual difference being O(α_s^{4}). The authors also release fast numerical codes for the associated Wilson coefficients, making the three-loop VFNS immediately usable in global PDF fits and collider analyses.

Core claim

At three-loop order the single- and two-mass VFNS, built from massive operator matrix elements that dress massless parton densities and generate heavy-quark distributions, reproduces the complete non-power-suppressed heavy-flavor corrections to deep-inelastic structure functions at large Q^{2}; a renormalization-group analysis proves the two representations differ only at O(α_s^{4}).

What carries the argument

Massive operator matrix elements (OMEs) that convert N_F-flavor massless PDFs into N_F+1 or N_F+2 PDFs (including heavy-quark densities) while preserving the renormalization-group equation; the resulting VFNS structure functions then match the asymptotic fixed-flavor expressions.

Load-bearing premise

Power corrections of order (m_Q^{2}/Q^{2})^k must be negligible, which the authors take to hold only for Q^{2} greater than roughly ten times the heavy-quark mass squared.

What would settle it

A direct numerical comparison of the three-loop VFNS F_{2} against the full fixed-flavor asymptotic result (including all non-power-suppressed terms) at a common large Q^{2}; any discrepancy larger than the expected O(α_s^{4}) residual would falsify the claimed equivalence.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Global PDF analyses can now evolve and match charm and bottom distributions at full three-loop accuracy.
  • Collider predictions that rely on heavy-quark PDFs (Higgs, top, electroweak bosons) inherit consistent N^{3}LO heavy-flavor matching.
  • Non-singlet extractions of α_s from future EIC data can include complete single- and two-mass three-loop corrections.
  • Fast public codes for charged- and neutral-current three-loop Wilson coefficients become available for phenomenology.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Two-mass OMEs of size ~½ the pure T_F^{2} pieces imply that simultaneous charm-bottom matching cannot be approximated by successive single-mass steps at the percent level.
  • The same OME technology can be reused for other factorization schemes (e.g., ACOT-style) once their renormalization-group consistency is verified.
  • Polarized three-loop VFNS results open a parallel path for spin-dependent PDF fits at N^{3}LO.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

0 major / 4 minor

Summary. The manuscript constructs the variable flavor number scheme (VFNS) to three-loop order for both single-mass and two-mass heavy-quark effects. Massless parton densities are modified by the complete set of three-loop massive operator matrix elements (OMEs), and heavy-quark PDFs for charm and bottom are introduced via the matching relations (1)–(7). A renormalization-group argument (Section 3, Eqs. (12)–(13)) shows that the resulting VFNS representation of F2 agrees with the fixed-flavor asymptotic heavy-flavor corrections up to O(a_s^4), so that all non-power-suppressed contributions are captured at large Q^2. Numerical x-space implementations of the single-mass OMEs and of a series of charged- and neutral-current Wilson coefficients are supplied as ancillary codes.

Significance. Promotion of the VFNS from two to three loops, including genuine two-mass OMEs, is a concrete advance for precision QCD phenomenology. The RG identity that establishes equivalence with the fixed-flavor asymptotic result is a clean, load-bearing derivation that rests on previously published multi-loop OMEs and massless Wilson coefficients. The release of fast, high-precision numerical codes for the OMEs and Wilson coefficients is a practical strength that immediately enables N3LO analyses at colliders and at the EIC. The explicit restriction to the asymptotic regime (power corrections neglected for Q^2 ≳ 10 m_Q^2) is stated clearly and does not undermine the claim as formulated.

minor comments (4)
  1. Section 2, after Eq. (8): the polarized case is stated to be treated in the Larin scheme, but the matching relations (1)–(7) are written only for the unpolarized densities; a short remark that the polarized analogues follow by the same construction (with the appropriate polarized OMEs) would improve clarity.
  2. Figure 1 caption: the denominators are given as Σ_NF=3 or G_NF=3; it would help the reader to state explicitly that these are the pure massless N_F=3 distributions evaluated at the same scale.
  3. Section 5: the two-mass numerical codes are announced as forthcoming “very soon”; a brief statement of the expected release format (or a pointer to a repository) would make the claim more concrete for users.
  4. References [41] and [43] are cited for the Wilson-coefficient codes; ensuring that the arXiv identifiers or DOIs are final before publication would avoid broken links.

Circularity Check

0 steps flagged

No load-bearing circularity; heavy self-citation of prior OME results is ordinary input, and the three-loop RG equivalence is an independent matching identity.

full rationale

The paper's central claim (Section 3, Eqs. 12-13) is that the renormalization-group operator applied to F2 yields δF2 = F_total - F_VFNS = O(a_s^4), so the three-loop VFNS completely captures the non-power-suppressed heavy-flavor corrections at large Q^{2}. This follows directly from the standard operator-product construction: the matching relations (1)-(7) absorb the already-computed massive OMEs into the PDFs, which are then convolved with the known massless Wilson coefficients. The OMEs themselves (cited throughout Section 2 from the authors' prior multi-loop calculations) are independent Feynman-diagram results whose assumptions do not include the VFNS matching identity; they are therefore legitimate external inputs rather than a self-referential loop. No parameter is fitted and then re-predicted, no uniqueness theorem is imported to forbid alternatives, and no ansatz is smuggled via citation. The asymptotic restriction (power corrections neglected for Q^{2} ≳ 10 m_Q^{2}) is stated explicitly and does not create circularity. The derivation is therefore self-contained against its own inputs; the only residual is ordinary self-citation of computational building blocks, which raises the score by at most one point.

Axiom & Free-Parameter Ledger

0 free parameters · 4 axioms · 0 invented entities

The paper rests on standard QCD renormalization-group machinery, previously computed massive operator matrix elements, and the conventional neglect of power-suppressed mass terms at large Q^2. No free parameters are fitted; no new entities are postulated.

axioms (4)
  • domain assumption The Callan-Symanzik / renormalization-group equation holds for the deep-inelastic structure functions, implying dF2/dt = 0.
    Invoked in Section 3 to prove that the VFNS and fixed-flavor representations differ only by O(a_s^4).
  • domain assumption Power corrections O((m_Q^2/Q^2)^k), k≥1, may be neglected for Q^2 ≫ m_Q^2 (practically Q^2 ≳ 10 m_Q^2 for F2).
    Stated in the Introduction and used throughout to justify the asymptotic VFNS picture.
  • domain assumption All single- and two-mass three-loop operator matrix elements previously computed by the collaboration are correct.
    The entire construction of the VFNS densities (Eqs. 1-7) is built from those OMEs.
  • domain assumption Polarized quantities are treated in the Larin scheme.
    Explicitly adopted in Section 2; requires consistent evolution of polarized PDFs and Wilson coefficients in that scheme.

pith-pipeline@v1.1.0-grok45 · 16050 in / 2269 out tokens · 20380 ms · 2026-07-11T07:36:41.727437+00:00 · methodology

0 comments
read the original abstract

We describe the variable flavor number scheme to three-loop order, which modifies the massless parton densities by single- and two-mass effects and introduces heavy-quark parton distribution functions for charm and bottom. A renormalization group analysis shows the validity of this picture at large scales $Q^2$, where it resembles the non-power-suppressed heavy-flavor corrections completely. We also provide numerical implementations of a series of charged and neutral current Wilson coefficients.

Figures

Figures reproduced from arXiv: 2607.05235 by A. Behring, A. De Freitas, A. von Manteuffel, C. Schneider, J. Ablinger, J. Bl\"umlein, K. Sch\"onwald.

Figure 1
Figure 1. Figure 1: Left panels: relative heavy-quark contributions due to charm to Δ NS,+ , Σ, 𝐺 and 𝑓𝑐+𝑐¯ . Right panels: relative heavy-quark contributions due to bottom to Δ NS,+ , Σ, 𝐺 and 𝑓𝑏+𝑏¯, both at 𝑄 2 = 100 GeV2 . The massive OMEs have been truncated to the perturbative order 𝑂(𝑎 𝑘 𝑠 ). The denominators of the ratios are Σ NF=3 or 𝐺 NF=3 , respectively. Dashed lines: 𝑂(𝑎𝑠). Dotted lines: contributions to 𝑂(𝑎 2 𝑠 )… view at source ↗
Figure 2
Figure 2. Figure 2: The unpolarized distributions 𝑥[𝑐(𝑥, 𝑄2 ) − 𝑐(𝑥, 𝑄2 )]; from Ref. [7]. 5. Two-mass operator matrix elements The two-mass operator matrix elements were all calculated to three-loop order in the unpolarized and polarized cases. They consist of factorizable and, from three-loop order, also non-factorizable contributions. The emergence of two masses complicates the respective integrals and functional represent… view at source ↗
Figure 3
Figure 3. Figure 3: The contributions of the two-mass OMEs 𝐴 two 𝑔𝑞 and 𝐴 two,PS 𝑄𝑞 normalized to the complete con￾tributions 𝑂(𝑇 2 𝐹 ) in the unpolarized case. Dotted lines: 𝑄 2 = 30GeV2 Full lines: 𝑄 2 = 50GeV2 Dashed lines: 𝑄 2 = 100GeV2 Dash-dotted lines: 𝑄 2 = 1000GeV2 ; from Refs. [33, 43]. implied by the massive renormalization group equation. These schemes describe the deep-inelastic structure functions in the region … view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

61 extracted references · 49 linked inside Pith

  1. [1]

    Callan, Jr.,Broken scale invariance in scalar field theory, Phys

    C.G. Callan, Jr.,Broken scale invariance in scalar field theory, Phys. Rev. D2(1970) 1541–1547

  2. [2]

    Symanzik,Small distance behavior in field theory and power counting, Commun

    K. Symanzik,Small distance behavior in field theory and power counting, Commun. Math. Phys.18 (1970) 227–246. 6Analyticresultsforthefullphasespacewereobtainedat𝑂(𝛼 2 𝑠 )inRefs.[9,21,54,55]andcanbeusedtodetermine the scale𝑄2 0 estimating the agreement with the asymptotic representations at the per-cent level. 7 The VFNS to three-loop orderJ. Blümlein

  3. [3]

    M. Buza, Y. Matiounine, J. Smith and W.L. van Neerven,Charm electroproduction viewed in the variable flavor number scheme versus fixed order perturbation theory, Eur. Phys. J. C1(1998) 301– 320 [hep-ph/9612398]

  4. [4]

    Blümlein, A

    J. Blümlein, A. De Freitas, C. Schneider and K. Schönwald,The Variable Flavor Number Scheme at Next-to-Leading Order, Phys. Lett. B782(2018) 362–366 [arXiv:1804.03129 [hep-ph]]

  5. [5]

    Bierenbaum, J

    I. Bierenbaum, J. Blümlein, A. De Freitas, A. Goedicke, S. Klein and K. Schönwald,𝑂(𝛼2 𝑠) polarized heavy flavor corrections to deep-inelastic scattering at𝑄2 ≫𝑚 2, Nucl. Phys. B988(2023) 116114 [arXiv:2211.15337 [hep-ph]]

  6. [6]

    Ablinger, A

    J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel, C. Schneider and K. Schön- wald,The single-mass variable flavor number scheme at three-loop order, JHEP03(2026) 248 [arXiv:2510.02175 [hep-ph]]

  7. [7]

    A.Behring, J.Blümlein, A.DeFreitas, A.vonManteuffel, C.SchneiderandK.Schönwald,Theheavy quark-antiquark asymmetry in the variable flavor number scheme, Phys. Lett. B876(2026) 140411 [arXiv:2512.13508 [hep-ph]]

  8. [8]

    Ablinger, J

    J. Ablinger, J. Blümlein, A. De Freitas, A. Hasselhuhn, C. Schneider and F. Wißbrock,Three Loop Massive Operator Matrix Elements and Asymptotic Wilson Coefficients with Two Different Masses, Nucl. Phys. B921(2017) 585–688 [arXiv:1705.07030 [hep-ph]]

  9. [9]

    M. Buza, Y. Matiounine, J. Smith, R. Migneron and W.L. van Neerven,Heavy quark coefficient functions at asymptotic values𝑄2 ≫𝑚 2, Nucl. Phys. B472(1996) 611–658 [hep-ph/9601302]

  10. [10]

    M.A.G.Aivazis,J.C.Collins,F.I.OlnessandW.K.Tung,Leptoproductionofheavyquarks.2.AUnified QCD formulation of charged and neutral current processes from fixed target to collider energiesPhys. Rev. D50(1994) 3102–3118 [hep-ph/9312319]

  11. [11]

    Collins,Hard scattering factorization with heavy quarks: A General treatment, Phys

    J.C. Collins,Hard scattering factorization with heavy quarks: A General treatment, Phys. Rev. D58 (1998) 094002 [hep-ph/9806259]

  12. [12]

    Cacciari, M

    M. Cacciari, M. Greco and P. Nason,The𝑝𝑇 spectrum in heavy-flavour hadroproduction, JHEP05 (1998) 007 [hep-ph/9803400]

  13. [13]

    R. S. Thorne and R. G. Roberts,An Ordered analysis of heavy flavor production in deep inelastic scattering, Phys. Rev. D57(1998) 6871–6898 [hep-ph/9709442]

  14. [14]

    S.Forte,E.Laenen,P.NasonandJ.Rojo,Heavyquarksindeep-inelasticscattering,Nucl.Phys.B834 (2010) 116–162 [arXiv:1001.2312 [hep-ph]]

  15. [15]

    Bierenbaum, J

    I. Bierenbaum, J. Blümlein and S. Klein,Mellin Moments of the𝑂(𝛼3 𝑠 )Heavy Flavor Contributions to unpolarized Deep-Inelastic Scattering at𝑄2 ≫𝑚 2 and Anomalous Dimensions, Nucl. Phys. B820 (2009) 417–482 [arXiv:0904.3563 [hep-ph]]

  16. [16]

    E.Witten,HeavyQuarkContributionstoDeepInelasticScattering,Nucl.Phys.B104(1976)445–476

  17. [17]

    Watson,Spin Spin Asymmetries In Inclusive Muon Proton Charm Production, Z

    A.D. Watson,Spin Spin Asymmetries In Inclusive Muon Proton Charm Production, Z. Phys. C12 (1982) 123–125

  18. [18]

    M. Buza, Y. Matiounine, J. Smith and W.L. van Neerven,𝑂(𝛼2 𝑠 )corrections to polarized heavy flavor production at𝑄2 ≫𝑚 2.Nucl. Phys. B485(1997) 420–456 [hep-ph/9608342]. 8 The VFNS to three-loop orderJ. Blümlein

  19. [19]

    I.Bierenbaum,J.BlümleinandS.Klein,Two-LoopMassiveOperatorMatrixElementsandUnpolarized Heavy Flavor Production at Asymptotic Values𝑄2 ≫𝑚 2 Nucl. Phys. B780(2007) 40–75 [hep- ph/0703285]

  20. [20]

    Bierenbaum, J

    I. Bierenbaum, J. Blümlein and S. Klein,The Gluonic Operator Matrix Elements at𝑂(𝛼2 𝑠 )for DIS Heavy Flavor Production, Phys. Lett. B672(2009) 401–406 [arXiv:0901.0669 [hep-ph]]

  21. [21]

    J.Blümlein,G.FalcioniandA.DeFreitas,TheComplete𝑂(𝛼 2 𝑠 )Non-SingletHeavyFlavorCorrections totheStructureFunctions𝑔 𝑒 𝑝 1,2 (𝑥, 𝑄 2),𝐹 𝑒 𝑝 1,2,𝐿 (𝑥, 𝑄 2),𝐹 𝜈(¯𝜈) 1,2,3 (𝑥, 𝑄 2)andtheAssociatedSumRules,Nucl. Phys. B910(2016) 568–617 [arXiv:1605.05541 [hep-ph]]

  22. [22]

    J.Ablinger,J.Blümlein,S.Klein,C.SchneiderandF.Wißbrock,The𝑂(𝛼 3 𝑠 )MassiveOperatorMatrix Elements of𝑂(𝑁 𝐹)for the Structure Function𝐹2 (𝑥, 𝑄 2)and Transversity, Nucl. Phys. B844(2011) 26–54 [arXiv:1008.3347 [hep-ph]]

  23. [23]

    Behring, I

    A. Behring, I. Bierenbaum, J. Blümlein, A. De Freitas, S. Klein and F. Wißbrock,The logarithmic contributions to the𝑂(𝛼 3 𝑠 )asymptotic massive Wilson coefficients and operator matrix elements in deeply inelastic scatteringEur. Phys. J. C74(2014) 9, 3033 [arXiv:1403.6356 [hep-ph]]

  24. [24]

    Ablinger, J

    J. Ablinger, J. Blümlein, A. De Freitas, A. Hasselhuhn, A. von Manteuffel, M. Round, C. Schneider and F. Wißbrock,The Transition Matrix Element𝐴𝑔𝑞 (𝑁)of the Variable Flavor Number Scheme at 𝑂(𝛼 3 𝑠 ), Nucl. Phys. B882(2014) 263–288 [arXiv:1402.0359 [hep-ph]]

  25. [25]

    Ablinger, A

    J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. Hasselhuhn, A. von Manteuffel, M. Round, C. Schneider, and F. Wißbrock,The 3-Loop Non-Singlet Heavy Flavor Contributions and Anomalous Dimensions for the Structure Function𝐹2 (𝑥, 𝑄 2)and Transversity, Nucl. Phys. B886(2014) 733–823 [arXiv:1406.4654 [hep-ph]]

  26. [26]

    Ablinger, A

    J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel and C. Schneider,The 3-loop puresingletheavyflavorcontributionstothestructurefunction𝐹 2 (𝑥, 𝑄 2)andtheanomalousdimension, Nucl. Phys. B890(2014) 48–151 [arXiv:1409.1135 [hep-ph]]

  27. [27]

    J.Ablinger,A.Behring,J.Blümlein,A.DeFreitas,A.vonManteuffel,C.SchneiderandK.Schönwald, The three-loop single mass polarized pure singlet operator matrix element, Nucl. Phys. B953(2020) 114945 [arXiv:1912.02536 [hep-ph]]

  28. [28]

    Blümlein, A

    J. Blümlein, A. De Freitas, M. Saragnese, C. Schneider and K. Schönwald,Logarithmic contributions to the polarized𝑂(𝛼3 𝑠 )asymptotic massive Wilson coefficients and operator matrix elements in deeply inelastic scattering, Phys. Rev. D104(2021) no.3, 034030 [arXiv:2105.09572 [hep-ph]]

  29. [29]

    Behring, J

    A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel, K. Schönwald and C. Schneider,The polarized transition matrix element𝐴𝑔𝑞 (𝑁)of the variable flavor number scheme at𝑂(𝛼 3 𝑠 ), Nucl. Phys. B964(2021) 115331 [arXiv:2101.05733 [hep-ph]]

  30. [30]

    Ablinger, A

    J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. Goedicke, A. von Manteuffel, C. Schneider andK.Schönwald,Theunpolarizedandpolarizedsingle-massthree-loopheavyflavoroperatormatrix elements𝐴 𝑔𝑔,𝑄 andΔ𝐴 𝑔𝑔,𝑄, JHEP12(2022) 134 [arXiv:2211.05462 [hep-ph]]

  31. [31]

    J.Ablinger,A.Behring,J.Blümlein,A.DeFreitas,A.vonManteuffel,C.SchneiderandK.Schönwald, Thefirst–orderfactorizablecontributionstothethree–loopmassiveoperatormatrixelements𝐴 (3) 𝑄𝑔 and Δ𝐴 (3) 𝑄𝑔, Nucl. Phys. B999(2024) 116427 [arXiv:2311.00644 [hep-ph]]. 9 The VFNS to three-loop orderJ. Blümlein

  32. [32]

    J.Ablinger,A.Behring,J.Blümlein,A.DeFreitas,A.vonManteuffel,C.SchneiderandK.Schönwald, The non-first-order-factorizable contributions to the three-loop single-mass operator matrix elements 𝐴(3) 𝑄𝑔 andΔ𝐴 (3) 𝑄𝑔.Phys. Lett. B854(2024) 138713 [arXiv:2403.00513 [hep-ph]]

  33. [33]

    J.Ablinger,J.Blümlein,A.DeFreitas,C.SchneiderandK.Schönwald,Thetwo-masscontributiontothe three-looppuresingletoperatormatrixelement,Nucl.Phys.B927(2018)339–367[arXiv:1711.06717 [hep-ph]]

  34. [34]

    Ablinger, J

    J. Ablinger, J. Blümlein, A. De Freitas, A. Goedicke, C. Schneider and K. Schönwald,The Two-mass Contribution to the Three-Loop Gluonic Operator Matrix Element𝐴(3) 𝑔𝑔,𝑄, Nucl. Phys. B932(2018) 129–240 [arXiv:1804.02226 [hep-ph]]

  35. [35]

    Ablinger, J

    J. Ablinger, J. Blümlein, A. De Freitas, M. Saragnese, C. Schneider and K. Schönwald,The three-loop polarized pure singlet operator matrix element with two different masses, Nucl. Phys. B952(2020) 114916 [arXiv:1911.11630 [hep-ph]]

  36. [36]

    Ablinger, J

    J. Ablinger, J. Blümlein, A. De Freitas, A. Goedicke, M. Saragnese, C. Schneider and K. Schönwald, The two-mass contribution to the three-loop polarized gluonic operator matrix element𝐴(3) 𝑔𝑔,𝑄, Nucl. Phys. B955(2020) 115059 [arXiv:2004.08916 [hep-ph]]

  37. [37]

    Ablinger, J

    J. Ablinger, J. Blümlein, A. De Freitas, A. von Manteuffel, C. Schneider and K. Schönwald,The two-mass contributions to the three-loop massive operator matrix elementse𝐴(3) 𝑄𝑔 andΔ e𝐴 (3) 𝑄𝑔, JHEP01 (2026) 111 [arXiv:2510.09403 [hep-ph]]

  38. [38]

    S.A.Larin,TheRenormalizationoftheaxialanomalyindimensionalregularization,Phys.Lett.B303 (1993) 113–118 [hep-ph/9302240]

  39. [39]

    J.BlümleinandM.Saragnese,Next-to-next-to-leadingorderevolutionofpolarizedpartondensitiesin the Larin scheme, Phys. Rev. D110(2024) no.3, 034006 [arXiv:2405.17252 [hep-ph]]

  40. [40]

    J.Blümlein,P.Marquard,C.SchneiderandK.Schönwald,Themasslessthree-loopWilsoncoefficients for the deep-inelastic structure functions F2, F𝐿, xF3 and g1, JHEP11(2022) 156 [arXiv:2208.14325 [hep-ph]]

  41. [41]

    J.Ablinger,A.Behring,J.Blümlein,A.DeFreitas,A.vonManteuffel,C.SchneiderandK.Schönwald, Thethree-loopsingle-massheavy-flavorcorrectionstothestructurefunctions𝐹 2(𝑥, 𝑄 2)and𝑔 1(𝑥, 𝑄 2), Phys. Lett. B (2026) in print [arXiv:2509.16124 [hep-ph]]

  42. [42]

    J.A.M.Vermaseren,A.VogtandS.Moch,TheThird-orderQCDcorrectionstodeep-inelasticscattering by photon exchange, Nucl. Phys. B724(2005) 3–182 [hep-ph/0504242]

  43. [43]

    J.Ablinger,A.Behring,J.Blümlein,A.DeFreitas,A.vonManteuffel,C.SchneiderandK.Schönwald, The complete three-loop unpolarized and polarized massive operator matrix elements and asymptotic Wilson coefficients, PoS (RADCOR2025) 076 [arXiv:2602.10334 [hep-ph]]

  44. [44]

    Zhu,The N3LO twist-2 matching of helicity TMDs and SIDIS q∗ spectrum, JHEP02(2026) 115 [arXiv:2509.01655v4 [hep-ph] 07 Feb 2026]

    Y.J. Zhu,The N3LO twist-2 matching of helicity TMDs and SIDIS q∗ spectrum, JHEP02(2026) 115 [arXiv:2509.01655v4 [hep-ph] 07 Feb 2026]

  45. [45]

    Behring, J

    A. Behring, J. Blümlein, A. De Freitas, A. Hasselhuhn, A. von Manteuffel and C. Schneider,O(𝛼3 𝑠) heavy flavor contributions to the charged current structure function𝑥𝐹3 (𝑥, 𝑄 2)at large momentum transferPhys. Rev. D92(2015) no.11, 114005 [arXiv:1508.01449 [hep-ph]]. 10 The VFNS to three-loop orderJ. Blümlein

  46. [46]

    Behring, J

    A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel and C. Schneider,The 3-Loop Non-Singlet Heavy Flavor Contributions to the Structure Function𝑔1(𝑥, 𝑄 2)at Large Momentum Transfer, Nucl. Phys. B897(2015) 612–644 [arXiv:1504.08217 [hep-ph]]

  47. [47]

    Behring, J

    A. Behring, J. Blümlein, G. Falcioni, A. De Freitas, A. von Manteuffel and C. Schneider,Asymp- totic 3-loop heavy flavor corrections to the charged current structure functions𝐹𝑊 + −𝑊 − 𝐿 (𝑥, 𝑄 2)and 𝐹𝑊 + −𝑊 − 2 (𝑥, 𝑄 2), Phys. Rev. D94(2016) no.11, 114006 [arXiv:1609.06255 [hep-ph]]

  48. [48]

    Blümlein, A

    J. Blümlein, A. Hasselhuhn and T. Pfoh,The𝑂(𝛼2 𝑠 )heavy quark corrections to charged current deep- inelastic scattering at large virtualities, Nucl. Phys. B881(2014) 1-41 [arXiv:1401.4352 [hep-ph]]

  49. [49]

    Caola, G

    F. Caola, G. Gambuti and M. Link,Analytic results for heavy-quark contributions to charged-current DIS at NNLO, [arXiv:2606.21702 [hep-ph]]

  50. [50]

    Gehrmann and E

    T. Gehrmann and E. Remiddi,Numerical evaluation of harmonic polylogarithms, Comput. Phys. Commun.141(2001) 296–312 [hep-ph/0107173]

  51. [51]

    J.Ablinger,J.Blümlein,M.RoundandC.Schneider,NumericalImplementationofHarmonicPolylog- arithms to Weight w = 8, Comput. Phys. Commun.240(2019) 189–201 [arXiv:1809.07084 [hep-ph]]

  52. [52]

    Ruijl, T

    B. Ruijl, T. Ueda and J. Vermaseren,FORM version 4.2[arXiv:1707.06453 [hep-ph]]

  53. [53]

    Blümlein and C

    J. Blümlein and C. Schneider,Analytic computing methods for precision calculations in quantum field theory, Int. J. Mod. Phys. A33(2018) no.17, 1830015 [arXiv:1809.02889 [hep-ph]]

  54. [54]

    J.Blümlein,A.DeFreitas,C.G.RaabandK.Schönwald,Theunpolarizedtwo-loopmassivepuresinglet Wilson coefficients for deep-inelastic scattering, Nucl. Phys. B945(2019) 114659 [arXiv:1903.06155 [hep-ph]]

  55. [55]

    Blümlein, C.G

    J. Blümlein, C.G. Raab and K. Schönwald,The Polarized Two-Loop Massive Pure Singlet Wilson Coefficient for Deep-Inelastic Scattering, Nucl. Phys. B948(2019) 114736 [arXiv:1904.08911 [hep- ph]]

  56. [56]

    d’Enterriaet al.,The strong coupling constant: state of the art and the decade ahead, J

    D. d’Enterriaet al.,The strong coupling constant: state of the art and the decade ahead, J. Phys. G51 (2024) no.9, 090501 [arXiv:2203.08271 [hep-ph]]

  57. [57]

    S.Alekhin,J.Blümlein,K.Daum,K.LipkaandS.Moch,Precisecharm-quarkmassfromdeep-inelastic scattering, Phys. Lett. B720(2013) 172–176 [arXiv:1212.2355 [hep-ph]]

  58. [58]

    D. Boer, M. Diehl, R. Milner, R. Venugopalan, W. Vogelsang, D. Kaplan, H. Montgomery, S. Vigdor, A. Accardi and E. C. Aschenauer,et al. Gluons and the quark sea at high energies: Distributions, polarization, tomography, [arXiv:1108.1713 [nucl-th]]

  59. [59]

    J.BlümleinandM.Saragnese,TheN 3LOscheme-invariantQCDevolutionofthenon-singletstructure functions𝐹 NS 2 (𝑥, 𝑄 2)and𝑔 NS 1 (𝑥, 𝑄 2), Phys. Lett. B820(2021) 136589 [arXiv:2107.01293 [hep-ph]]

  60. [60]

    Gehrmann, A

    T. Gehrmann, A. von Manteuffel, V. Sotnikov and T.Z. Yang,The four-loop non-singlet splitting functions in QCD, [arXiv:2604.09534 [hep-ph]]

  61. [61]

    Guyot, A

    C. Guyot, A. Milsztajn, A. Oraou, A. Staudte, K.M. Teichert, and M. Virchaux,A new fixed-target experiment for precise tests of QCD, CERN/SPSC 89/56, March 14 1989. 11