⚛️ phenomenology

Analytic results for heavy-quark contributions to charged-current DIS at NNLO

This paper presents analytic next-to-next-to-leading-order QCD corrections for heavy-quark production in charged-current deep-inelastic scattering with exact charm mass dependence, expressing two-heavy-quark contributions via Goncharov polylogarithms and three-heavy-quark contributions via elliptic structures and rapidly convergent expansions for efficient numerical evaluation.

Original authors: Fabrizio Caola, Giulio Gambuti, Martin Link

Published 2026-06-23

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Fabrizio Caola, Giulio Gambuti, Martin Link

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

The Big Picture: Zooming in on the Universe's Building Blocks

Imagine you are trying to understand the structure of a brick wall (the proton) by throwing tennis balls (electrons or neutrinos) at it. This is what Deep Inelastic Scattering (DIS) is: a high-energy game of "pin the tail on the donkey" where we smash particles together to see what's inside.

To make sense of the debris flying off, physicists use a map called a Parton Distribution Function (PDF). This map tells us how likely it is to find a specific piece of the wall (like a quark) at a specific spot. However, these maps are only as good as the math used to draw them. If the math is too rough, the map is blurry, and we can't see the details of the wall.

The Problem: The "Heavy" Quark Mystery

Inside the proton, there are light quarks (like up and down) and heavy quarks (like charm).

Light quarks are like feathers; they are easy to ignore when the energy is high.
Heavy quarks are like bowling balls. When the energy of the collision is low or medium, the weight of the bowling ball matters a lot. It changes the trajectory of the debris.

For a long time, physicists had a perfect map for the light feathers (calculated to a very high level of precision, called NNLO). But for the heavy bowling balls (charm quarks), the map was either:

Too simple: It ignored the weight of the ball, which is bad for low-energy collisions.
Too messy: It existed only as a computer number (a "black box" calculation) that couldn't be easily checked or improved.

The Solution: A New, Exact Blueprint

This paper presents a brand new, fully analytic blueprint for how heavy charm quarks behave in these collisions. The authors didn't just guess; they calculated the exact mathematical formula that describes the interaction, keeping the "weight" of the charm quark in the equation at all times.

Think of it like this:

Previous methods were like estimating the flight path of a bowling ball by pretending it was a feather, or by running a simulation that gave you a number but no formula.
This paper gives you the exact physics equation for the bowling ball, valid whether it's moving slowly or fast.

The Technical Hurdles: Navigating a Mathematical Jungle

Calculating this was incredibly hard. The authors had to solve complex equations involving loops of energy that look like tangled spaghetti.

The "Polylogarithm" Forest: Most of the calculations resulted in a type of mathematical function called "Goncharov polylogarithms." The authors managed to tame these functions, turning them into a format that is easy for computers to read and evaluate quickly.
The "Elliptic" Cave: Some parts of the calculation (specifically when three heavy quarks are involved) were so complex they couldn't be described by standard functions. They entered a realm called "elliptic structures." The authors didn't get stuck here; they built a special "tunnel" (using Chen iterated integrals and series expansions) to get through this cave and still get a usable answer.

Why This Matters (According to the Paper)

The authors state that having this exact formula is crucial for two main reasons:

Precision: It allows scientists to create much sharper maps of the proton's interior, specifically for the "strange" quark distribution, which is currently a major source of uncertainty in our understanding of the universe.
Flexibility: Because the result is a clean formula (not just a black-box number), it can be easily plugged into global models to test new theories of physics.

The Validation: Checking the Work

The authors didn't just trust their own math. They ran several "stress tests":

They checked that their new formula matches known results from simpler, older calculations.
They checked that it matches the "massless" version (where the quark has no weight) when the energy is very high.
They compared their results against existing numerical data and found they agreed perfectly.

Summary

In short, this paper is a masterclass in high-precision math. The authors have built a new, exact tool to measure how heavy particles behave inside protons during collisions. They navigated through complex mathematical landscapes (polylogarithms and elliptic curves) to provide a result that is both mathematically beautiful and practically useful for refining our understanding of the fundamental building blocks of matter.

Technical Summary: Analytic Results for Heavy-Quark Contributions to Charged-Current DIS at NNLO

Problem Statement
Deep Inelastic Scattering (DIS) is a primary tool for determining Parton Distribution Functions (PDFs), which are critical for precision Standard Model measurements and Beyond Standard Model searches at the Large Hadron Collider (LHC). While global PDF fits require high-precision theoretical inputs, current determinations of heavy-flavor contributions in Charged-Current (CC) DIS are limited. Specifically, the Next-to-Next-to-Leading Order (NNLO) QCD corrections for heavy-quark (charm) production in CC DIS have previously been available only in numerical form. This numerical nature hinders their efficient integration into global PDF fits and complicates systematic studies of mass effects and power-suppressed contributions, particularly in the intermediate and low $Q^2$ regions ( $Q^2 \lesssim 100 \text{ GeV}^2$ ) relevant for neutrino DIS data. Existing analytic results are largely restricted to the asymptotic limit ( $Q^2 \gg m_c^2$ ) or the massless case.

Methodology
The authors present a fully analytic computation of the NNLO QCD corrections to heavy-flavor CC DIS, retaining the exact dependence on the charm quark mass ( $m_c$ ). The calculation employs the following technical framework:

Process Definition: The study covers the inclusive CC DIS process $l + P \to l' + X$ , mediated by a $W$ boson. The calculation is performed in a mixed renormalization scheme: the proton contains three active light flavors ( $n_f=3$ ), while the strong coupling $\alpha_s$ is renormalized with four flavors ( $n_f+n_h=4$ ) to account for virtual and real heavy-quark emissions.
Reverse Unitarity and IBP: The calculation utilizes the reverse-unitarity framework to treat real-emission contributions as cuts of loop integrals. Integration-by-parts (IBP) identities are used to reduce the resulting Feynman integrals to a minimal set of Master Integrals (MIs).
Differential Equations and Canonical Bases: The MIs are solved using the method of differential equations. The authors construct a canonical basis for the integrals, transforming the system into an $\epsilon$ $ϵ$ -factorized form.
- Polylogarithmic Sector: For contributions involving up to two final-state heavy quarks, the solutions are expressed in terms of manifestly real Goncharov Polylogarithms (GPLs). The authors handle square-root singularities by rationalizing the kinematic variables and mapping the results to GPLs.
- Elliptic Sector: The contribution involving three final-state heavy quarks (RR3) involves elliptic structures (specifically a three-mass sunrise graph). These are expressed in terms of Chen iterated integrals.
Threshold and Boundary Conditions: To fix boundary constants, the authors perform deep Frobenius expansions around the threshold limit ( $y_r \to 0$ ) and the massless limit ( $z \to 0$ ). They utilize high-precision numerical evaluations (via AMFlow) near the threshold to match against these expansions, fixing all integration constants analytically.
Renormalization: The calculation includes UV renormalization of the charm mass (on-shell scheme) and the strong coupling ( $\overline{\text{MS}}$ scheme), as well as collinear renormalization of the PDFs. Special attention is paid to the Larin scheme for handling $\gamma_5$ in dimensional regularization, including the necessary finite renormalization of the axial current.

Key Contributions and Results

Complete Analytic Coefficient Functions: The paper provides the first analytic expressions for the complete partonic coefficient functions for the structure functions $F_2$ , $F_L$ , and $F_3$ at NNLO, retaining exact mass dependence. This covers quark and gluon channels with up to three heavy quarks in the final state.
Elliptic Structures: The work provides a general representation for the three-heavy-quark contribution (RR3) in terms of Chen iterated integrals. Furthermore, it offers rapidly convergent series expansions valid in the perturbative region $Q \gtrsim 5 \text{ GeV}$ , enabling fast and flexible numerical evaluation.
Numerical Validation: The results are validated against:
- Known exact results at lower orders (LO and NLO).
- Massless NNLO coefficient functions available in the literature.
- Leading-power asymptotic expansions in the limit $Q^2 \gg m_c^2$ .
- Existing numerical NNLO results (though a high-accuracy one-to-one comparison was limited by the unavailability of source code for previous numerical results).
Implementation: The authors have implemented these coefficient functions in a FORTRAN code interfaced with LHAPDF. They demonstrate the flexibility of the results by computing hadronic structure functions and comparing them with asymptotic expansions, showing that exact mass effects are significant at low $Q^2$ but converge to the asymptotic limit at high $Q^2$ .

Significance and Claims
The authors claim that their results provide the missing ingredients for a full NNLO analysis of CC DIS that retains exact charm-quark mass dependence. By moving from numerical to analytic results, they enable:

More robust and efficient incorporation of heavy-flavor effects into global PDF fits.
Systematic studies of mass effects and power-suppressed contributions in the intermediate and low $Q^2$ regions, which are crucial for neutrino physics.
A reliable description of charm production where mass effects cannot be approximated by simple rescaling or asymptotic expansions.

The paper concludes that while the current work focuses on charm production, the framework is extendable to cases involving bottom quarks or non-inclusive scenarios (e.g., fragmentation), though these are left for future investigation. The results are provided in computer-readable formats to facilitate their use by the community.