Wilson-Coefficient C₍Qq₎³,¹ in QCD & SMEFT

• Wilson-Coefficient C₍Qq₎³,¹ is a key term that represents three-loop heavy-flavor corrections in deep-inelastic scattering and a dimension-6 four-quark operator in SMEFT.
• It encapsulates quantum effects from integrating out heavy quarks and incorporates complex analytic structures involving nested harmonic sums and transcendental constants.
• Its precise determination refines theoretical predictions, facilitating improved matching of DIS structure functions and top-quark production observables with experimental data.

The Wilson coefficient $C_{Qq}^{3,1}$ appears in two distinct but related contexts within high-energy QCD and effective field theory analyses: as a key three-loop heavy-flavor correction in deep-inelastic scattering (DIS), and as a dimension-6 four-quark operator coefficient in Standard Model effective field theory (SMEFT), especially in top-quark physics. In both frameworks, $C_{Qq}^{3,1}$ encapsulates the quantum effects of integrating out heavy quarks or incorporating new contact operators and is critical for matching theoretical predictions to experimental precision at three-loop order.

1. Definition and Operator Structure

In the context of deep-inelastic scattering in QCD, $C_{Qq}^{3,1}$ denotes the pure-singlet heavy-flavor Wilson coefficient at three-loop order. Explicitly, it arises in the perturbative expansion of the heavy-flavor contributions to the DIS structure functions such as $F_2(x,Q^2)$. The general form is

$$C_{Qq}^{\rm PS}(N, Q^2/m^2) = \sum_{k=1}^\infty a_s^k\, C_{Qq}^{(k),\rm PS}(N, Q^2/m^2),\quad a_s = \frac{\alpha_s}{4\pi},$$

with the three-loop coefficient written as

$$C_{Qq}^{(3),\rm PS}(N) = a_s^3\bigg\{\sum_{i=1}^6 \mathcal{C}_i P_i^{(3)}(N)\bigg\},$$

where $P_i^{(3)}(N)$ denote linear combinations of nested harmonic sums, rational functions in $N$, and zeta-values up to $\zeta_5$ (Ablinger et al., 2024).

In the SMEFT (Standard Model Effective Field Theory), the notation $C_{Qq}^{3,1}$ refers to the Wilson coefficient of a four-quark operator that modifies processes such as $t$-channel single-top production,

$$\mathcal{O}_{Qq}^{3,1} = \sum_{i=1,2}\Big[(\bar Q\gamma^\mu\tau^Iq_i)(\bar q_i\gamma_\mu\tau^IQ) + \tfrac{1}{6} (\bar Q\gamma^\mu q_i)(\bar q_i\gamma_\mu Q) - \tfrac{1}{6}(\bar Q\gamma^\mu\tau^Iq_i)(\bar q_i\gamma_\mu\tau^IQ)\Big],$$

with $C_{Qq}^{3,1}$ parameterizing its strength relative to the new physics scale $\Lambda$ (Collaboration, 8 Jan 2026).

2. Analytic Structure in Deep-Inelastic Scattering

The full analytic structure of $C_{Qq}^{3,1}$ is determined by both the three-loop operator matrix element $A_{Qq}^{(3),\mathrm{PS}}(N)$ and lower-order massless Wilson coefficients. Its form is expanded as

$$C_{Qq}^{(3,1)}(N, Q^2/m^2) = \sum_{k=0}^3 \ln^k(Q^2/m^2)\,C^{(3,1),(k)}(N),$$

where the coefficients $C^{(3,1),(k)}(N)$ are linear combinations of QCD color factors, harmonic sums up to weight $5-k$, and transcendental constants $\zeta_2,\ldots,\zeta_5$ (Blümlein et al., 2013). For example, the leading cubic log term is

$$C^{(3,1),(3)}(N) = T_F C_F \frac{16}{27} S_1(N) + T_F C_A \frac{4}{27} S_1(N) + T_F^2 \frac{8}{27} S_1(N).$$

The constant piece, $C^{(3,1),(0)}(N)$, is $A_{Qq}^{(3),\mathrm{PS}}(N)$ and contains nested sums up to weight five and products with zeta-values.

Upon Mellin inversion, the result in Bjorken-$x$ space involves harmonic polylogarithms $H_{b_1,\dots,b_k}(x)$ (weight up to five), with rational prefactors in $x$ and $\ln(Q^2/m^2)$ (Ablinger et al., 2024, Ablinger et al., 2014).

3. Color and Transcendental Structure

The full three-loop coefficient is a sum over six color structures: $C_F C_A^2 T_F$, $C_F^2 C_A T_F$, $C_F^3 T_F$, $n_f C_F C_A T_F^2$, $n_f C_F^2 T_F^2$, and $C_F T_F^3$. Each color channel contributes rational functions in $N$, products of nested harmonic sums, and transcendental constants including all zeta-values up to $\zeta_5$. Example structure for the $C_F C_A^2 T_F$ channel: $$P^{(3)}_{C_A C_A}(N) = A_1(N)\zeta_5 + A_2(N)\zeta_4 + A_3(N)\zeta_3 + \cdots + \sum_{a+b\leq 5} c_{a,b}(N) S_a(N) S_b(N) + \cdots.$$ The transcendental weight of individual terms ranges up to five, matching the maximum weight of harmonic sums and polylogarithms at this order. All rational coefficients are functions of $N$ given in closed analytic form (Ablinger et al., 2024, Ablinger et al., 2014).

4. Role in Deep-Inelastic Scattering Observables

$C_{Qq}^{3,1}$ encodes the three-loop heavy-flavor contribution to the pure-singlet channel in DIS structure functions. In the asymptotic region $Q^2 \gg m^2$, the structure function $F_2(x,Q^2)$ receives contributions of the form

$$F_2(x,Q^2) \supset a_s^3\,C_{Qq}^{(3),\rm PS}(x,Q^2/m^2) \otimes \Sigma(x),$$

where $\Sigma(x)$ is the singlet quark distribution. The coefficient is most important at small $x$, where it can reach corrections at the few-percent level (e.g., at $x=10^{-4}$, $C_{Qq}^{(3,1)}(x)\sim 0.3$ for $\alpha_s = 0.2, Q^2/m^2=50$) (Blümlein et al., 2013), and diminishes for large $x$.

The analytic structure supports theoretical matching with high-precision experimental data, as both small-$x$ and threshold ($x \to 1$) expansions are in agreement with factorization predictions and resummation limits.

5. Calculation Techniques and Renormalization

The derivation of $C_{Qq}^{3,1}$ at three loops involves several advanced techniques:

• Feynman diagram generation via QGRAF.
• Dirac and color algebra projections onto twist-2 operators.
• Integration-by-parts reduction and master integral evaluation using differential equations, moments "guessing," and symbolic summation (Sigma, HarmonicSums).
• Ultraviolet and collinear renormalization in the on-shell/$\overline{\rm MS}$ scheme.
• Consistency checks via fixed-$N$ moments, sum rules, and resummed limits.

Light-quark masses and currents are treated in the $\overline{\rm MS}$ scheme, heavy-quark mass is on-shell, and the strong coupling is $\overline{\rm MS}$. In polarized cases, the Larin $\gamma_5$ prescription is used (Ablinger et al., 2024).

6. Significance in SMEFT and Constraints from Collider Data

Beyond neutral- and charged-current DIS, $C_{Qq}^{3,1}$ plays a central role in SMEFT, encoding new physics in four-quark contact operators that interfere with Standard Model $W$-exchange in $t$-channel single-top production. In the Warsaw basis, all relevant dimension-6 contributions at leading order can be absorbed into $C_{Qq}^{3,1}$, which modifies both inclusive and differential cross-sections.

Recent ATLAS results constrain $C_{Qq}^{3,1}$ using full Run 2 differential data and an EFT approach: $$-0.12\;\text{TeV}^{-2} < C_{Qq}^{3,1} < 0.12\;\text{TeV}^{-2}$$ at 95% CL, with acceptance and selection efficiencies modeled via dedicated MC samples and unfolded with the covariance of the SM measurement. These represent a significant tightening over previous bounds and constitute the first differential t-channel spectra limits in top-quark contact interactions (Collaboration, 8 Jan 2026).

7. Summary Table: Key Features Across Contexts

Context	Definition/Role	Principal Expressions/Observables
QCD/ DIS	3-loop pure-singlet heavy-flavor Wilson coefficient	Linear combinations of nested harmonic sums and zeta values in Mellin N and x space for DIS structure functions ($F_2, F_1$) (Ablinger et al., 2024, Blümlein et al., 2013, Ablinger et al., 2014)
SMEFT/ Top	Wilson coefficient of SU(2)×U(1) invariant four-quark operator	Appears in dimension-6 operator basis affecting single-top production, constrained by LHC data (Collaboration, 8 Jan 2026)

The anomalous dimension structure, color combinatorics, and analytic continuation in $N$ (Mellin space) are common technical themes, ensuring the theoretical consistency and predictive power of $C_{Qq}^{3,1}$ both in precision QCD calculations and in effective field theory searches for physics beyond the Standard Model.