⚛️ high-energy experiments

The variable flavor number scheme to three-loop order

This paper describes the extension of the variable flavor number scheme to three-loop order, which incorporates mass effects and heavy-quark parton distribution functions, validates this approach at large scales through renormalization group analysis, and provides numerical implementations of relevant Wilson coefficients.

Original authors: J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel, C. Schneider, K. Schönwald

Published 2026-07-07

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: J. Ablinger, A. Behring, J. Blümlein, A. De Freitas, A. von Manteuffel, C. Schneider, K. Schönwald

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is built out of tiny, fundamental Lego bricks. In the world of high-energy physics, these bricks are called quarks and gluons. Sometimes, these bricks are light and easy to move around (like up and down quarks), and sometimes they are heavy and sluggish (like charm and bottom quarks).

This paper is about creating a better instruction manual for how these heavy bricks behave when they are smashed together at incredibly high speeds, like in a particle collider.

Here is the breakdown of their work using simple analogies:

1. The Problem: The "Heavy Suit" Dilemma

When scientists study these particle collisions, they usually use a set of rules called a "scheme" to calculate what happens.

The Old Way (Fixed Flavor): Imagine you are counting the bricks in a box. If the box contains heavy bricks, you have to stop and manually count them every single time you look at the box. It's accurate, but slow and tedious.
The New Way (Variable Flavor Number Scheme - VFNS): The authors propose a smarter way. They say, "If the collision is happening at a very high energy (a very fast speed), the heavy bricks act almost exactly like the light ones." So, instead of treating them as special, heavy exceptions, we can upgrade our instruction manual to include them as regular parts of the team.

2. The Upgrade: From Two Loops to Three Loops

In physics, calculations are often done in "loops." Think of a loop as a layer of detail or a level of precision in a video game.

Two Loops: The previous manuals had two layers of detail. They were good, but the picture was a little blurry.
Three Loops: This paper upgrades the manual to three layers of detail. This is like switching from a standard definition TV to 4K Ultra HD. The picture of how these heavy quarks behave becomes incredibly sharp and precise.

3. The "Massive" Math (Operator Matrix Elements)

To make this upgrade work, the authors had to calculate some very complex math called "Operator Matrix Elements" (OMEs).

The Analogy: Imagine you are baking a cake. You know how to bake a cake with just flour (light quarks). But now you want to add chocolate chips (heavy quarks). You need a new recipe that tells you exactly how the chocolate chips change the texture of the cake.
Single-Mass vs. Two-Mass:
- Single-Mass: They calculated the recipe for adding just one type of heavy chip (like only chocolate).
- Two-Mass: They also calculated the recipe for adding two different types of heavy chips at the same time (like chocolate and peanut butter). This is much harder because the chips interact with each other in complex ways. The authors solved the math for both scenarios up to the highest level of precision (three loops).

4. The Result: A Perfect Match

The authors checked their new manual against the old, heavy-handed way of counting.

The Finding: When the energy is high enough (the "large scales" mentioned in the paper), their new, smart manual matches the old, heavy manual perfectly. The "heavy" quarks are now fully integrated into the system without needing special, clumsy exceptions.
The Limit: They note that if the energy is too low (like a slow-motion collision), the heavy quarks don't act like light ones anymore, and this manual stops working. But for the high-speed collisions they are interested in, it works perfectly.

5. Why This Matters (According to the Paper)

The authors didn't just do the math; they built computer codes (digital tools) that other scientists can use.

They released these tools so other researchers can use them to analyze data from particle colliders.
The Goal: By using these precise tools, scientists can measure the "glue" holding the universe together (the strong force) and the weight of the heavy quarks with much greater accuracy.
Specific Recommendation: The paper suggests that for the most precise measurements of the universe's fundamental forces using future data (like from the Electron-Ion Collider), scientists should use a specific "non-singlet" analysis method that now has all the necessary three-loop corrections available.

In summary: This paper provides the most precise "instruction manual" yet for how heavy particles behave in high-speed collisions. It upgrades the math from a blurry 2D sketch to a sharp 3D model, allowing scientists to predict and measure the behavior of the universe's building blocks with unprecedented accuracy.

Technical Summary: The Variable Flavor Number Scheme to Three-Loop Order

Problem Statement
The paper addresses the theoretical description of deep-inelastic scattering (DIS) structure functions in the presence of heavy quarks (charm and bottom). While Fixed Flavor Number Schemes (FFNS) are accurate at low virtualities, they fail to resum large logarithmic terms at high scales ( $Q^2 \gg m_Q^2$ ). Conversely, standard Variable Flavor Number Schemes (VFNS) account for universal heavy-flavor contributions by modifying massless parton distribution functions (PDFs) and introducing heavy-quark PDFs. However, previous implementations were largely limited to two-loop order or single-mass approximations. The challenge lies in consistently extending the VFNS to three-loop order, incorporating both single-mass (one heavy quark) and two-mass (simultaneous charm and bottom) effects, while ensuring the validity of the renormalization group equation (RGE) and the correct matching at flavor thresholds.

Methodology
The authors utilize the framework of Operator Matrix Elements (OMEs) to connect massive and massless parton densities. The core methodology involves:

Renormalization Group Analysis: The paper demonstrates that the VFNS representation agrees with the Fixed Flavor Number Scheme (FFNS) at large scales $Q^2$ up to three-loop accuracy ( $O(a_s^3)$ ), neglecting power corrections of order $O((m^2/Q^2)^k)$ .
Calculation of OMEs: The study relies on massive OMEs calculated to three-loop order in both unpolarized and polarized cases. These include:
- Single-mass OMEs: Previously calculated in Refs. [22–32].
- Two-mass OMEs: Calculated in Refs. [8, 33–37], comprising both factorizing contributions (products of lower-order terms) and genuine non-factorizing contributions emerging at three loops (e.g., $A_{qq,Q}^{NS,two}$ , $A_{Qq}^{PS,two}$ ).
Mellin Space Relations: The authors derive relations in Mellin space where massless parton densities are transformed into massive densities via convolutions with OMEs ( $A_{ij}$ ). These relations define the evolution of light-quark, gluon, and heavy-quark densities across flavor thresholds ( $N_F \to N_F+1$ or $N_F+2$ ).
Numerical Implementation: The authors developed fast and precise numerical codes in C++ and Fortran. These codes utilize analytic small- and large- $x$ expansions and cubic splines for remainder parts, allowing for efficient use in QCD fitting.

Key Contributions and Results

Three-Loop VFNS Framework: The paper presents the complete structure of the VFNS to three-loop order for both single-mass and two-mass cases. It explicitly details how massless Wilson coefficients are paired with massive OMEs to describe structure functions like $F_2$ , $F_L$ , $xF_3$ , and $g_1$ .
Two-Mass Effects: The work illustrates the emergence of non-factorizing two-mass contributions at the three-loop level. Figure 3 demonstrates that these two-mass OMEs constitute a significant fraction (approximately $O(1/2)$ ) of the complete contributions proportional to the color factor $T_F^2 C_F A$ .
Heavy-Quark PDFs: The study quantifies the impact of heavy-quark corrections on parton distributions. Figure 1 shows the convergence of these corrections from $O(a_s)$ to $O(a_s^3)$ for non-singlet, singlet, and gluon distributions. Additionally, the paper highlights the existence of $f_Q - f_{\bar{Q}}$ distribution functions starting at three-loop order, which are found to be numerically small (Figure 2).
Numerical Tools: The authors provide public numerical codes for:
- Unpolarized and polarized single-mass OMEs to three loops.
- Massless and asymptotic massive Wilson coefficients to three loops.
- Charged and neutral current Wilson coefficients.
  These tools are expressed in terms of harmonic polylogarithms and are designed for integration into QCD fitting frameworks.

Significance and Claims
The paper claims that the proposed VFNS framework provides a consistent and complete description of heavy-flavor corrections at large virtualities, fully respecting the renormalization group equation. The authors emphasize that this framework is essential for high-precision measurements of the strong coupling constant ( $\alpha_s$ ) and heavy-quark masses at hadron colliders, particularly when gluon and flavor-singlet contributions are involved.

Furthermore, the authors advocate for a "strict non-singlet analysis at N3LO" for future high-precision extractions of $\alpha_s$ from world data, including future Electron-Ion Collider (EIC) data. They assert that the necessary theoretical ingredients—massless, single-mass, and two-mass corrections to three loops—are now available to realize this long-standing proposal. The work serves as a foundational step for more precise calculations of processes at hadron colliders and for the consistent treatment of heavy-flavor effects in global QCD analyses.

1. The Problem: The "Heavy Suit" Dilemma

2. The Upgrade: From Two Loops to Three Loops

3. The "Massive" Math (Operator Matrix Elements)

4. The Result: A Perfect Match

5. Why This Matters (According to the Paper)

More like this