The four-loop non-singlet splitting functions in QCD

This paper presents the first fully analytic calculation of the four-loop non-singlet splitting functions in perturbative QCD, confirming prior partial results and enabling the extraction of four-loop virtual and rapidity anomalous dimensions for logarithmic resummation.

The Gist

Imagine the proton (the core of every atom in your body) not as a solid marble, but as a bustling, chaotic city. Inside this city, there are tiny citizens called quarks and messengers called gluons zooming around, constantly bumping into each other and changing lanes.

To understand how this city behaves when we smash it together at high speeds (like in the Large Hadron Collider), physicists use a map called a Parton Distribution Function (PDF). This map tells us the probability of finding a specific citizen at a specific speed.

But here's the catch: the map isn't static. As you zoom in closer (increasing the resolution), the citizens start interacting more, and the map changes. The rules that govern how this map evolves are called Splitting Functions. Think of these as the "traffic laws" of the subatomic city. They dictate how likely a quark is to split into two, or how a gluon is to turn into a quark.

The Big Achievement: The "Four-Loop" Traffic Law

For decades, physicists have been trying to write down these traffic laws with increasing precision.

• One-loop: The basic rules (like "stop at red lights").
• Two-loop: Adding details (like "yield to pedestrians").
• Three-loop: Getting very specific (like "speed limits on rainy days").

This paper presents the Four-Loop version. This is the "Holy Grail" of precision for a specific type of interaction called the non-singlet case.

What is "Non-Singlet"?
Imagine the city has two types of neighborhoods:

1. The Singlet District: Where everyone mixes freely (quarks and gluons swapping identities constantly). This is messy and hard to calculate.
2. The Non-Singlet District: A more exclusive neighborhood where the citizens keep their identities (quarks stay quarks, anti-quarks stay anti-quarks). This is the "clean" part of the city that the authors focused on.

How They Did It: The "Mathematical Microscope"

The authors didn't just guess these rules; they built a massive mathematical microscope to look at the interactions.

1. The Blueprint (Feynman Diagrams): They generated about 16,000 different blueprints (diagrams) showing every possible way particles could interact. Some of these diagrams were incredibly complex, involving loops of particles that look like tangled spaghetti.
2. The Tangled Knots (Elliptic Geometry): Usually, these knots can be untangled with standard math. But at this level of precision, they found some knots that required a new kind of geometry (called "elliptic") to solve. It's like realizing that to untie a specific knot, you need a tool you've never used before.
3. The Translation (Mellin Moments): They translated the messy, real-world interactions into a language of numbers (Mellin moments) to solve the equations, and then translated the answers back into the "traffic laws" (Splitting Functions).

Why Does This Matter?

You might ask, "Why do we need the 4th level of precision?"

1. Eliminating Guesswork: Previously, scientists had to use "approximations" for the four-loop rules. It was like driving with a map that had some blank spots filled in with "maybe this way." Now, they have the exact, complete formula. There are no more blank spots.
2. The "Virtual" and "Rapidity" Dimensions: The paper also unlocked the exact mathematical form of two hidden quantities (called anomalous dimensions) that are crucial for predicting how particles behave when they are moving very fast or very slow. It's like discovering the exact engine specs of a car that was previously only known by its top speed.
3. Future Experiments: As we build bigger, more powerful particle colliders, the data will be so precise that old, approximate maps will lead to errors. This new, ultra-precise map ensures that when we discover new particles or forces, we aren't just seeing a glitch in our math.

The Bottom Line

This paper is the completion of a massive puzzle. The authors have finally written down the exact, analytical rules for how quarks evolve at the highest level of precision currently possible in our understanding of the universe.

They took a chaotic, tangled mess of 16,000 possibilities, untangled the knots using advanced math, and handed us a perfectly clear, error-free map for the subatomic world. This allows physicists to predict the future of particle collisions with a level of certainty that was previously impossible.

Technical Summary

1. Problem Statement

The evolution of Parton Distribution Functions (PDFs) with the energy scale $Q^2$ is governed by the DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) equations. The kernels of these equations are the splitting functions, $P(x)$, which describe the probability of a parton splitting into two others.

• Current Status: Splitting functions are known analytically up to three loops (NNLO). At four loops (N$^3$LO), only approximate expressions and partial analytical results (e.g., specific color structures or fixed Mellin moments) were available.
• The Gap: The lack of complete analytical expressions for the four-loop non-singlet splitting functions ($P^{\pm, V}_{ns}$) limits the precision of PDF fits and prevents a rigorous study of the asymptotic behavior of these functions at small and large momentum fractions ($x$). Furthermore, the complete analytical form is required to extract the four-loop virtual and rapidity anomalous dimensions necessary for high-order logarithmic resummation (N$^4$LL).

2. Methodology

The authors computed the splitting functions by calculating the anomalous dimensions of non-singlet twist-two operators in the Operator Product Expansion (OPE).

• Operator Matrix Elements (OME): The calculation focuses on the renormalization of non-singlet operators $O_{ns}$ inserted between off-shell quark states. The anomalous dimension $\gamma_{ns}$ is extracted from the single $\epsilon$-pole of the renormalization constant $Z_{ns}$ in dimensional regularization ($d = 4 - 2\epsilon$).
• Diagram Generation: They generated approximately 16,000 Feynman diagrams using Qgraf, selecting those contributing to the non-singlet anomalous dimensions. These diagrams involve up to four gluons and two external quarks.
• Algebraic Reduction:
  • Lorentz, Dirac, and color algebra were performed using FORM and Color.h.
  • To handle the symbolic spin $n$ in the operator $(\Delta \cdot D)^{n-1}$, they introduced an auxiliary parameter $t$ to convert powers of momenta into linear propagators.
• Integral Reduction:
  • The resulting scalar integrals were reduced to a set of Master Integrals (MIs) using Integration-by-Parts (IBP) identities.
  • The reduction utilized Reduze 2, Finred, and finite field sampling techniques.
  • The system involved $\sim 3$ million unreduced integrals, reduced to $\sim 6,000$ master integrals across 549 top-level sectors.
• Differential Equations (DEs):
  • Differential equations with respect to the auxiliary parameter $t$ were derived for $\sim 300$ top-level sectors.
  • Solutions were obtained as Laurent series in $\epsilon$ with coefficients being Taylor series in $t$.
  • Elliptic Geometry: A significant technical challenge was the appearance of integrals related to elliptic geometries (specifically a 4$\times$4 coupled system) at the four-loop level, which had not been encountered in previous three-loop calculations. The authors successfully decoupled these using $d$-log forms and elliptic integrals.
• Reconstruction:
  • The four-loop bare OMEs were computed up to $\mathcal{O}(4000)$ moments.
  • Using the HarmonicSums package, the authors reconstructed the anomalous dimensions $\gamma^{(3)}_{ns}$ in $n$-space as linear combinations of harmonic sums.
  • An inverse Mellin transformation was applied to obtain the splitting functions $P(x)$ in $x$-space.

3. Key Contributions

• First Full Analytical Result: This is the first derivation of fully analytic expressions for all three non-singlet splitting functions ($P^+_{ns}, P^-_{ns}, P^V_{ns}$) at the four-loop order.
• New Color Structures: The authors computed previously unknown contributions, including:
  • Non-fermionic terms for $P^{\pm}_{ns}$.
  • The $n_f$ (number of quark flavors) term for $P^s_{ns}$.
  • These correspond to complex color structures involving $C_A^3 C_F$, $C_A^2 C_F^2$, $C_A C_F^3$, $C_F^4$, and terms with $d_{abcd}$ tensors.
• Elliptic Integrals: The work demonstrates that while elliptic integrals appear in the intermediate steps of the calculation (specifically in the homogeneous blocks of the differential equations), they cancel out in the final renormalized anomalous dimensions, which can be expressed entirely in terms of harmonic sums.
• Extraction of Anomalous Dimensions: The results allow for the first time the extraction of the analytic form of the four-loop virtual anomalous dimension ($B_4$) and the four-loop rapidity anomalous dimension ($\gamma_{rap}$).

4. Key Results

• Analytic Form: The splitting functions are expressed in terms of harmonic polylogarithms (up to weight 6) and transcendental constants (up to weight 7).
• Asymptotic Behavior ($x \to 1$):
  • The authors derived the exact coefficient $B_4$ (the four-loop virtual anomalous dimension), which enters the resummation of soft-gluon logarithms.
  • The result confirms the all-order conjecture relating the coefficients of the splitting function near $x=1$ to the cusp anomalous dimension and the $\beta$-function.
• Asymptotic Behavior ($x \to 0$):
  • The paper provides the explicit leading logarithmic behavior ($\log^k x$) for $x \to 0$.
  • The results show that the approximations used in previous N$^3$LO PDF fits are reliable within their quoted uncertainties, though small deviations exist in the low-$x$ region for $P^+_{ns}$.
• Numerical Validation: The analytical results were compared against fixed-moment calculations and existing approximate codes (e.g., from Ref. [19]), showing full agreement.
• Rapidity Anomalous Dimension: A complete analytic expression for the four-loop rapidity anomalous dimension is provided in the appendix, replacing previous numerical-only determinations.

5. Significance

• Precision Physics: The availability of exact four-loop non-singlet splitting functions eliminates approximation uncertainties in N$^3$LO PDF fits, crucial for high-precision Standard Model tests and New Physics searches at the LHC and future colliders.
• Resummation: The extraction of the analytic four-loop virtual and rapidity anomalous dimensions enables N$^4$LL (Next-to-Next-to-Next-to-Leading Logarithmic) resummation of soft and collinear logarithms, significantly improving the theoretical description of processes like Higgs production and Drell-Yan scattering.
• Theoretical Insight: The successful handling of elliptic integrals in the intermediate steps and their subsequent cancellation provides deep insight into the structure of QCD amplitudes at high loop orders, suggesting that the final physical observables may retain a simpler structure (harmonic sums) than the intermediate mathematical objects.
• Future Applications: The authors provide simple power-logarithmic approximations with accuracy better than $10^{-6}$ across the full $x$ range, facilitating immediate implementation in global PDF fitting frameworks.

In summary, this work represents a major milestone in perturbative QCD, closing the gap between partial four-loop results and the complete analytical understanding required for the next generation of precision phenomenology.