---
title: Single-Mass Heavy-Flavor Contributions
url: https://www.emergentmind.com/topics/single-mass-heavy-flavor-contributions
type: topic
---

# Single-Mass Heavy-Flavor Contributions

Single-mass heavy-flavor contributions are the deep-inelastic scattering (DIS) corrections in which one heavy-quark mass scale is retained at a time—typically charm or bottom—while simultaneous two-mass effects are excluded. In the asymptotic regime \(Q^2 \gg m^2\), these contributions admit an operator-product-expansion description in which heavy-flavor Wilson coefficients factorize into massless Wilson coefficients and massive operator matrix elements (OMEs). This framework governs the three-loop treatment of inclusive neutral-current and charged-current structure functions such as \(F_2\), \(F_L\), \(F_1\), \(g_1\), and \(xF_3\), and it also provides the transition functions used in variable-flavor-number schemes (VFNS) [1007.0375].

## 1. Definition, observables, and kinematic regime

“Single-mass” denotes the situation in which the heavy-flavor sector contains one heavy-quark mass scale only, for example \(m_c\) or \(m_b\), and no simultaneous two-mass effects such as charm-bottom mixing. In the leading-twist treatment of DIS, the coefficient functions are decomposed into a massless part and a heavy-flavor part,
\[
\mathbb{C}_{2,L}^{(i)}\!\left(x,\frac{Q^2}{\mu^2},\frac{m_q^2}{\mu^2}\right)
=
C_{2,L}^{(i)}\left(x,\frac{Q^2}{\mu^2}\right)
+
H_{2,L}^{(i)}\left(x,\frac{Q^2}{\mu^2},\frac{m_q^2}{\mu^2}\right),
\]
so that the heavy contribution appears as a distinct perturbative sector in \(F_2\), \(F_L\), and, in the polarized case, \(g_1\) [2107.09350].

For neutral-current DIS, the measured structure function is decomposed into light- and heavy-flavor pieces,
\[
F_2(x,Q^2)=F_2^{\rm light}(x,Q^2)+F_2^c(x,Q^2)+F_2^b(x,Q^2),
\]
with the heavy terms becoming particularly important at small Bjorken \(x\) [2509.16124]. For charged-current DIS, the single-mass heavy sector has an additional feature: besides heavy-quark pair production, there is also single heavy-flavor excitation, specifically \(s\to c\), which contributes already at parton level in the non-singlet charged-current structure functions [1609.06255].

The entire asymptotic construction is restricted to sufficiently large virtuality. Several papers state that the approximation is numerically reliable for \(F_2\) once \(Q^2/m^2 \gtrsim 10\), while it is not intended for the threshold region \(Q^2 \sim m^2\) [1202.2700]. A recurrent caution in the literature is that asymptotic heavy-flavor expressions capture twist-2 contributions and omit power-suppressed terms \(O(m^2/Q^2)\), so they are high-scale descriptions rather than exact finite-\(Q^2\) results [1007.0375].

## 2. Asymptotic factorization and the role of massive OMEs

The central structural statement is the asymptotic factorization of heavy-flavor Wilson coefficients into massless Wilson coefficients and massive OMEs. In Mellin space, one form used for the neutral-current case is
\[
{\sf H}_{j,(2,L)}(N)=A_{ij}(N)\cdot C_{i,(2,L)}(N), \qquad i,j=q,g,
\]
where \(A_{ij}(N)\) contains the heavy-mass dependence and \(C_{i,(2,L)}(N)\) is massless; the Mellin convolution in \(x\)-space becomes ordinary multiplication in \(N\)-space [1007.0375]. Equivalent generic formulations appear throughout the three-loop literature,
\[
H_i\left(N_F,\frac{Q^2}{m^2}\right)
=
C_i\left(N_F+1,\frac{Q^2}{\mu^2}\right)\otimes
A\left(N_F,\frac{m^2}{\mu^2}\right),
\]
with \(\otimes\) the Bjorken-\(x\) convolution [2407.02006].

The OMEs are local twist-2 matrix elements between on-shell partonic states in the presence of one heavy mass,
\[
A_{ki}^{\sf S,NS}\Bigl(N,\frac{m^2}{\mu^2}\Bigr)
=
\langle i|O_k^{\sf S,NS}|i\rangle_H
=
\delta_{ki}+\sum_{l=1}^{\infty} a_s^l A_{ki}^{(l)}\Bigl(N,\frac{m^2}{\mu^2}\Bigr),
\qquad
a_s=\frac{\alpha_s}{4\pi},
\]
and they encode the full asymptotic mass dependence after renormalization and mass factorization [1007.0375]. For the unpolarized three-loop single-mass case, the relevant OMEs are
\[
A_{qq,Q}^{\rm NS},\quad A_{Qq}^{\rm PS},\quad A_{qq,Q}^{\rm PS},\quad
A_{Qg},\quad A_{qg,Q},\quad A_{gq,Q},\quad A_{gg,Q},
\]
with corresponding polarized OMEs in the helicity-dependent sector [2407.02006].

The renormalization logic is standard but highly constraining. The calculations use \(\overline{\mathrm{MS}}\) for coupling and operator renormalization and on-shell renormalization for the heavy-quark mass; after ultraviolet and mass renormalization, the remaining collinear singularities are removed by mass factorization [1007.0375]. A crucial consequence is that the pole and logarithmic terms of the three-loop OMEs are determined by anomalous dimensions up to three loops, \(\beta\)-function coefficients up to two loops, and lower-order OMEs, while the constant term is the genuinely new three-loop contribution [1007.0375]. This is why OME calculations simultaneously provide independent access to pieces of the three-loop anomalous dimensions.

## 3. Three-loop development and completion of the single-mass sector

The three-loop program developed in stages. Early NNLO asymptotic work established that massive three-loop OMEs account for all heavy-flavor terms except the power corrections \(O(m^2/Q^2)\) and reported fixed Mellin moments for all three-loop massive OMEs. That stage also yielded new all-\(N\) expressions for the \(O(n_f)\) terms, including the constant term \(a_{Qg}^{(3)}\), confirmed the corresponding fermionic pieces of the three-loop anomalous dimensions, and presented the first genuine three-loop ladder-type results for general Mellin moment \(N\) [1007.0375].

A later status report stated that six out of eight relevant three-loop OMEs had been computed,
\[
A_{qq,Q}^{(3), \rm PS},\quad A_{qg,Q}^{(3)},\quad A_{gq}^{(3)},\quad
A_{qq,Q}^{(3), \rm NS},\quad A_{qq,Q}^{(3), \rm NS,\,TR},\quad
A_{Qq}^{(3), \rm PS},
\]
while
\[
A_{gg,Q}^{(3)},\qquad A_{Qg}^{(3)}
\]
were still missing [1409.1804]. In parallel, the pure-singlet channel was advanced by a complete calculation of the asymptotic three-loop heavy-flavor corrections to \(F_2\) in the \(H_{2,q}^{\rm PS}\) sector together with the transition matrix element \(A_{Qq}^{(3),\rm PS}\) and an independent recalculation of the full three-loop pure-singlet anomalous dimension [1409.1135].

The gluonic transition sector was completed by the three-loop calculation of the unpolarized and polarized OMEs
\[
A_{gg,Q}(x,\mu^2), \qquad \Delta A_{gg,Q}(x,\mu^2),
\]
including their Mellin-space and \(x\)-space forms and their analytic continuation from even or odd moments into the complex plane [2211.05462]. A subsequent 2024 overview then stated that all single-mass OMEs and all asymptotic inclusive heavy-flavor Wilson coefficients had been completed to three-loop order for unpolarized and polarized DIS, together with the single-mass VFNS matching relations [2407.02006].

The next step was phenomenological rather than structural: the first numerical three-loop single-mass heavy-flavor corrections to the inclusive structure functions \(F_2(x,Q^2)\) and \(g_1(x,Q^2)\) were reported in 2025, providing a full NNLO numerical realization of the asymptotic single-mass program for both unpolarized and polarized inclusive DIS [2509.16124].

## 4. Analytic structures and computational technology

The analytic function space of single-mass heavy-flavor calculations is considerably richer than the one encountered in lower-order massless DIS. For the \(O(n_f)\) three-loop terms, the Mellin-space results are expressed through harmonic sums up to weight 4, and the basis can be reduced to six basic sums,
\[
S_1,\quad S_{2,1},\quad S_{-2,1},\quad S_{-3,1},\quad S_{2,1,1},\quad S_{-2,1,1},
\]
together with rational functions of \(N\) and zeta values [1007.0375]. Ladder-type diagrams already generate more complicated nested sums, including generalized harmonic sums in intermediate expressions [1007.0375].

The pure-singlet sector introduced a qualitatively new feature in DIS: generalized harmonic sums in the final Mellin-space result for \(A_{Qq}^{(3),\rm PS}\), and corresponding \(x\)-space structures involving not only ordinary harmonic polylogarithms in \(x\) but also functions with argument \(1-2x\) [1409.1135]. More broadly, later summaries emphasized that pure-singlet terms require generalized harmonic sums, whereas non-singlet and some \(A_{gq,Q}\) sectors can be expressed in terms of ordinary harmonic sums or harmonic polylogarithms only [2107.09350].

Additional diagram classes generate cyclotomic harmonic sums and cyclotomic harmonic polylogarithms. These arise from alphabets built from cyclotomic polynomials and extend the standard harmonic-polylogarithm alphabet \(\{1/x,1/(1-x),1/(1+x)\}\) [1202.2700]. The gluonic three-loop sector brings in finite binomial and inverse binomial sums in Mellin space and iterated integrals over square-root-valued alphabets in \(x\)-space, including plus distributions and denominator structures such as \((1\pm x)^{-k}\) with \(k=2,3\) [2211.05462]. The 2025 numerical analysis further noted that \(H_g^{\rm S}\) contains binomial sums, square-root iterated integrals, and \({}_2F_1\)-type solutions in \(x\)-space [2509.16124].

The computational pipeline combines large-scale symbolic reduction with specialized summation technology. Reported tools include QGRAF for diagram generation, FORM for algebraic manipulation, integration-by-parts reduction with Reduze2 implementing Laporta’s algorithm, and master-integral evaluation via generalized hypergeometric functions, Mellin-Barnes representations, hyperlogarithms, and differential equations. The nested sums are handled with Sigma, HarmonicSums, EvaluateMultiSums, SumProduction, and OreSys; later surveys also highlight the method of arbitrarily large Mellin moments and the method of guessing for recurrence reconstruction [1409.1804]. A common misconception addressed explicitly in this literature is that leading logarithms alone are sufficient: both the 2012 and 2024 analyses stress that leading small-\(x\) terms do not approximate the full result well by themselves, because subleading contributions and constant terms are numerically important [1202.2700].

## 5. Phenomenology, polarized observables, and flavor schemes

The phenomenological importance of single-mass heavy-flavor terms is greatest at small \(x\), where the heavy contribution can be a large fraction of the inclusive structure function and its scaling violations differ from the massless case [2509.16124]. For \(F_2\), the asymptotic approximation is stated to be already accurate at the \(\sim 1\%\) level when
\[
\frac{Q^2}{m_Q^2}>10
\]
at two loops, and this same high-\(Q^2\) regime suppresses higher-twist effects [2509.16124]. By contrast, the asymptotic description of \(F_L\) is much less accurate; a survey reports that power corrections remain important and that the pure asymptotic terms describe \(F_L\) only for much larger scales, around \(Q^2/m_q^2 \gtrsim 800\) [2107.09350].

The first complete three-loop numerical analysis for inclusive \(F_2\) used \(N_F=3\), PDFs from Alekhin et al., \(\alpha_s(M_Z^2)=0.1147\), and on-shell masses \(m_c=1.59\), \(m_b=4.78\). It found that, in the small-\(x\) region, the heavy-flavor fraction
\[
R=\frac{F_2^c+F_2^b}{F_2^{\rm light}+F_2^c+F_2^b}
\]
rises from about \(26\%\) at \(Q^2=25^2\) to about \(41\%\) at \(Q^2=1000^2\); around \(x\sim 0.2\) the heavy-flavor correction crosses through zero; and at very large \(x\), NNLO heavy-flavor corrections can be negative and as large as about \(-10\%\), while the total \(F_2\) remains positive [2509.16124]. In the same study, the gluonic coefficient \(H_g^{\rm S}\) dominates the charm contribution, \(H_q^{\rm PS}\) is sizable and negative, \(L_q^{\rm NS}\) dominates at large \(x\), and \(L_g^{\rm S}\) and \(L_g^{\rm PS}\) remain relevant at the percent level [2509.16124]. Earlier pure-singlet analyses likewise found that the heavy-flavor pure-singlet contribution to \(F_2\) is negative in the studied kinematics and significantly larger for charm than for bottom [1409.1135].

For polarized DIS, the same asymptotic single-mass machinery applies, but the three-loop calculation is carried out in the Larin scheme. The 2025 inclusive analysis reports that the gluonic coefficient \(\Delta H_g^{\rm S}\) is dominant also in \(g_1\), and its first moment vanishes through three loops,
\[
\int_0^1 dx\, \Delta H_g^{\rm S}(x,Q^2)=0.
\]
This vanishing first moment drives an oscillatory structure in \(x\)-space [2509.16124]. In the region where \(xg_1\) itself is small, the heavy-flavor ratio can be roughly \(-10\%\), \(-25\%\), and \(-95\%\) for \(Q^2=25^2, 100^2, 1000^2\), respectively; around \(x\sim 10^{-2}\) the ratio changes sign, around \(x\sim 0.3\) it becomes negative again, and yet the full \(g_1\) remains positive [2509.16124]. A 2024 study further notes that, at NNLO, Larin- and \(\overline{\mathrm{MS}}\)-scheme polarized evolution differ noticeably, with quark distributions differing by about \(10\)–\(15\%\) at small \(x\sim 10^{-3}\) and gluon distributions by about \(3\%\), so consistent polarized phenomenology requires PDFs evolved in the Larin scheme [2407.02006].

The same OMEs also determine matching across heavy-flavor thresholds in the VFNS, including the generation of the heavy-quark PDFs \(f_Q+f_{\bar Q}\) [2407.02006]. In this sense, single-mass asymptotic calculations are not only fixed-flavor predictions but also transition functions for an \((N_F+1)\)-flavor description [1007.0375]. A separate NNLO implementation in the S-ACOT-\(\chi\) general-mass scheme illustrates how single-mass heavy-flavor structure functions such as \(F_{2c}\) and \(F_{Lc}\) are assembled from massive and zero-mass ingredients using the rescaling variable
\[
\chi=x\left(1+\frac{4m_h^2}{Q^2}\right),
\]
which enforces threshold kinematics and yields smooth interpolation between fixed-flavor and zero-mass limits [1108.4008].

## 6. Charged-current sector, sum rules, and scope

Single-mass heavy-flavor contributions are not restricted to neutral-current \(F_2\) and \(g_1\). In charged-current non-singlet DIS, the asymptotic three-loop heavy-flavor corrections to
\[
F_L^{W^+-W^-}(x,Q^2), \qquad F_2^{W^+-W^-}(x,Q^2),
\]
and
\[
xF_3^{W^+}(x,Q^2)+xF_3^{W^-}(x,Q^2)
\]
have been calculated for general Mellin moment \(N\) and in \(x\)-space [1609.06255]. These observables contain two distinct heavy-flavor mechanisms: heavy-quark pair production, encoded in \(L\)-type Wilson coefficients, and single heavy-flavor excitation \(s\to c\), encoded in \(H\)-type coefficients [1508.01449].

The charged-current asymptotic results show that charm corrections to \(xF_1^{W^+-W^-}\) and \(F_2^{W^+-W^-}\) are typically in the range \(-8\%\lesssim \text{charm correction}\lesssim 0\%\), with the asymptotic approximation becoming reliable over a broader \(x\)-range as \(Q^2\) increases; at lower \(Q^2\), power corrections remain visible, especially for \(F_L\) and therefore \(F_2\) [1609.06255]. For \(xF_3^{W^+}+xF_3^{W^-}\), the charm effect is up to about \(+3\%\) at small \(x\) and about \(-3\%\) at large \(x\), while the \(O(a_s^3)\) correction relative to the massless three-flavor result is typically at the \(1\)–\(3\%\) level [1508.01449].

The first Mellin moment has a special status. For the Adler sum rule, the charged-current heavy-flavor analysis states that there are no QCD corrections and no quark-mass corrections; in the asymptotic heavy-flavor limit this follows because the relevant non-singlet OMEs vanish at \(N=1\) by fermion-number conservation, and the first moment of the corresponding massless Wilson coefficient also vanishes [1609.06255]. For the Gross–Llewellyn Smith sum rule, the heavy-flavor effect in the asymptotic region reduces essentially to the replacement
\[
N_F\to N_F+1
\]
in the massless coefficient, with CKM weights included [1508.01449]. The same pattern appears for the unpolarized Bjorken sum rule in the charged-current \(F_1\) sector [1609.06255].

The scope of the single-mass literature is also sharply delimited. The main three-loop results are asymptotic, twist-2, and usually formulated in the fixed-flavor-number scheme with \(N_F=3\) light flavors for phenomenological applications [2509.16124]. They are not intended for threshold kinematics, where \(Q^2\sim m^2\), and they exclude simultaneous two-mass effects; a 2024 overview states that two-mass corrections were nearly finished, while the single-mass sector was complete [2407.02006]. A plausible implication is that “single-mass heavy-flavor contributions” should be understood not as a universal finite-\(Q^2\) solution, but as the fully developed high-scale component of heavy-flavor QCD factorization, with direct relevance to precision determinations of \(\alpha_s(M_Z^2)\), heavy-quark masses, twist-2 PDFs, and threshold matching in modern DIS analyses [2509.16124].

Source: https://www.emergentmind.com/topics/single-mass-heavy-flavor-contributions