⚛️ high-energy theory

Three-Loop QCD corrections to the production of a Higgs boson and a Jet

This paper presents the computation of three-loop QCD corrections to the scattering amplitude for Higgs boson and three-parton production in the generalized leading color limit, expressing the results in terms of multiple polylogarithms to enable precise phenomenological predictions for Higgs-plus-jet production at hadron colliders and Higgs decay to three jets at lepton colliders.

Original authors: Xiang Chen, Xin Guan, Bernhard Mistlberger

Published 2026-07-09

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Xiang Chen, Xin Guan, Bernhard Mistlberger

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 a giant, chaotic dance floor where tiny particles called "partons" are constantly bumping into each other. At the center of this dance is the Higgs boson, a heavy, mysterious celebrity that gives other particles their mass. Scientists at the Large Hadron Collider (LHC) are trying to watch this celebrity dance with a partner (a "jet" of particles) to understand the rules of the dance floor perfectly.

To do this, they need to predict exactly how the dance looks using a set of mathematical rules called Quantum Chromodynamics (QCD). The problem is, the dance is so complicated that the rules are incredibly messy. Usually, scientists can only predict the dance steps for the first few moves (one or two loops). But to match the super-precise cameras of the LHC, they need to predict the dance three moves deep (three loops).

The Big Breakthrough
In this paper, the authors, Xiang Chen, Xin Guan, and Bernhard Mistlberger, have finally calculated the dance steps for a Higgs boson interacting with three partons all the way up to the three-loop level. Think of a "loop" as a layer of complexity in the math; adding a third layer is like trying to solve a Rubik's cube while juggling, blindfolded, on a moving train.

They didn't just guess; they computed the exact mathematical "scattering amplitudes" (the probability of the dance happening) for two specific scenarios:

The Decay: A Higgs boson breaking apart into three particles (like a firework exploding into three sparks).
The Collision: Two particles smashing together to create a Higgs boson and a third particle (like a billiard ball hitting two others to create a new, heavy ball).

The "Generalized Leading Color" Shortcut
Here is the tricky part: The math for these three-loop interactions is so incredibly complex that calculating every single possibility would take forever. To get around this, the authors used a clever shortcut called the "generalized leading color limit."

Imagine the particles have "colors" (not red, blue, or green, but a quantum property). In the real world, there are many ways these colors can mix. The authors decided to focus on the most common, dominant color combinations (the "leading" ones) and treat the number of colors as a huge number. They explicitly ruled out the tiny, rare color combinations that are suppressed by the square of the number of colors. They argue that this shortcut is a "highly effective approximation," capturing the main action while ignoring the background noise that would barely change the result (about a 10% difference at most). This allowed them to tackle the non-planar diagrams—the messy, tangled loops that usually break the math—without getting stuck.

The Result: A New Language for the Dance
The authors didn't just get a number; they expressed their results in a special mathematical language called "multiple polylogarithms." You can think of these as a universal translator that turns the messy, tangled math of the three-loop dance into a clean, readable script. This script is now ready for other scientists to use immediately to make predictions for the LHC.

How Sure Are They?
The authors are very confident in their work, but they are careful not to overstate it. They didn't just simulate it on a computer; they derived the results analytically (using pure math). To prove they didn't make a mistake, they ran a battery of checks:

They checked if their math followed the universal rules for how particles behave when they get very close together (collinear limit) or very slow (soft limit).
They compared their results to a different, simpler theory (N=4 super Yang-Mills) and found a perfect match in the most complex parts of the math.
They even used a different numerical method (the "auxiliary mass flow" method) to double-check their numbers, and the two methods agreed perfectly.

What This Means for the Future
This paper doesn't claim to have solved the entire mystery of the Higgs boson. Instead, it provides the essential "ingredients" needed to cook the next level of precision recipes. Specifically, these results are necessary to calculate the production of a Higgs boson with a jet at the next level of precision (N3LO) and to study how the Higgs decays into three jets at future lepton colliders (like the FCC-ee or CEPC).

In short, the authors have built the most detailed map yet of a very specific, very complex corner of the particle dance floor. They haven't mapped the whole universe, but they have cleared the fog over a critical section, allowing physicists to finally see the dance steps with the precision required to test if our current understanding of nature is truly correct.

Technical Summary: Three-Loop QCD Corrections to the Production of a Higgs Boson and a Jet

Problem Statement
The discovery of the Higgs boson has necessitated high-precision comparisons between experimental data from the Large Hadron Collider (LHC) and theoretical predictions to test the Standard Model and search for new physics. The dominant production mechanism for the Higgs boson at the LHC is gluon fusion, mediated by a virtual top-quark loop. To match the experimental precision expected in the High-Luminosity phase, theoretical predictions must be computed to very high orders in the strong coupling constant $\alpha_S$ . While inclusive Higgs production has been calculated up to next-to-next-to-next-to-leading order (N $^3$ LO), differential predictions involving hadronic jets remain challenging. Specifically, extending the calculation of Higgs boson production in association with a jet to N $^3$ LO requires the computation of three-loop scattering amplitudes for a Higgs boson and three partons. This paper addresses the calculation of these specific three-loop amplitudes, which are essential ingredients for such precision phenomenology.

Methodology
The authors compute the three-loop scattering amplitudes for the process of a Higgs boson ( $H$ ) interacting with three partons (gluons or quark-antiquark pairs). The calculation is performed in two kinematic regions:

Decay Region: $H \to ggg$ and $H \to gq\bar{q}$ .
Scattering Region: $gg \to Hg$ , $gq \to Hq$ , and $q\bar{q} \to Hg$ .

The computation utilizes the following technical framework:

Effective Field Theory (EFT): The top-quark loop is integrated out, treating the Higgs-gluon interaction via a point-like effective operator $H G_{\mu\nu}^a G_a^{\mu\nu}$ with a Wilson coefficient $C_0$ . This approximation is valid for the leading contributions.
Generalized Leading Color Limit: To manage the immense complexity of three-loop non-planar diagrams, the calculation is restricted to the generalized leading color limit. In this limit, the number of colors $N_c$ and light quark flavors $n_f$ are treated as large, neglecting terms suppressed by powers of $1/N_c^2$ . This includes terms proportional to $N_c^6, N_c^5 n_f, N_c^4 n_f^2, N_c^3 n_f^3$ for all-gluon external states, and corresponding leading terms for mixed states.
Dimensional Regularization: Calculations are performed in $d = 4 - 2\epsilon$ dimensions using the $\overline{\text{MS}}$ scheme to regulate ultraviolet and infrared singularities.
Diagram Generation and Reduction: Feynman diagrams are generated using qgraf. Gauge-invariant amplitudes are computed via tensor decomposition. The squared amplitudes are reduced to scalar integrals and classified into integral families.
Integration-by-Parts (IBP): The Blade package, utilizing the Laporta algorithm, block-triangular forms, and finite field techniques (via FiniteFlow), is employed to reduce the integrals to a set of master integrals.
Master Integrals: The master integrals are solved using the method of differential equations. The authors relate their integrals to known canonical bases in the literature and derive new canonical sets for previously unidentified families. Boundary conditions are determined using auxiliary mass flow methods (AMFlow) and analytic continuation.
Function Space: The results are expressed in terms of multiple polylogarithms.
Renormalization and Subtraction: Infrared and ultraviolet singularities are removed using universal subtraction factors ( $Z$ ) and renormalization of the coupling and Wilson coefficient, yielding finite remainder amplitudes.

Key Contributions and Results

First Computation: This work presents the first computation of four-particle scattering amplitudes at three-loop order in QCD involving a massive external particle (the Higgs) and non-planar diagrams.
Explicit Results: The authors provide the three-loop scattering matrix elements ( $M^{(3)}_X$ ) for the decay of a Higgs boson into three partons and the production of a Higgs boson with a parton in the scattering of two partons.
Analytic Form: The results are provided as Laurent expansions in the dimensional regulator $\epsilon$ , with coefficients expressed as ratios of polynomials and generalized polylogarithms of the kinematic variables ( $s, t, u$ ).
Finite Interference: The paper provides the finite interfered matrix elements (interference between three-loop and tree-level amplitudes) necessary for N $^3$ LO cross-section calculations. These are made available in ancillary files.
Validation: The results are validated through five independent checks:
1. Consistency with the universal pole structure of amplitudes.
2. Agreement with universal splitting amplitudes in the two-particle collinear limit.
3. Agreement with universal soft currents in the soft-gluon limit.
4. Reproduction of the leading transcendental part corresponding to the $\text{tr}(\phi^2)$ form factor in $\mathcal{N}=4$ Super Yang-Mills theory.
5. Numerical agreement with auxiliary mass flow evaluations in the decay region.

Significance
The paper claims that these results represent a decisive step toward achieving one-percent precision phenomenology at the LHC. The computed amplitudes are the necessary ingredients to:

Compute the N $^3$ LO production cross section for a Higgs boson in association with a hadronic jet at hadron colliders.
Compute the N $^3$ LO decay probability of a Higgs boson into three hadronic jets at future lepton colliders (e.g., FCC-ee, CEPC, ILC).
Serve as a key component for the inclusive Higgs production probability at the next order (N $^4$ LO).

The work demonstrates that the computation of non-planar three-loop amplitudes with massive external particles is feasible, pushing the frontier of perturbative QCD calculations. The availability of these results in terms of multiple polylogarithms makes them ready for immediate use in phenomenological predictions.

More like this