⚛️ phenomenology

Two-loop QCD amplitudes for $t\bar{t}W$ production at the LHC in the leading-colour approximation

This paper presents a numerical computation of two-loop QCD scattering amplitudes for $t\bar{t}W$ production at the LHC in the leading-colour approximation, utilizing a hybrid framework that combines numerical evaluation with analytic control to retain exact mass dependence and provide validated hard functions for next-to-next-to-leading-order cross-section calculations.

Original authors: Matteo Becchetti, Dhimiter Canko, Xiang Chen, Vsevolod Chestnov, Maximilian Delto, Sara Ditsch, Tiziano Peraro, Mattia Pozzoli, Chiara Savoini, Simone Zoia

Published 2026-08-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Matteo Becchetti, Dhimiter Canko, Xiang Chen, Vsevolod Chestnov, Maximilian Delto, Sara Ditsch, Tiziano Peraro, Mattia Pozzoli, Chiara Savoini, Simone Zoia

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 as a giant, cosmic LEGO set. For decades, physicists have been trying to figure out exactly how the pieces snap together. They have a rulebook called the Standard Model, which describes the tiny building blocks of matter and the forces that push and pull them. But just like trying to predict the exact path of a falling domino in a room full of wind, the math gets incredibly messy when you try to calculate what happens when these particles collide at super-high speeds.

To get a clear picture, scientists need to calculate these collisions with extreme precision. They start with a simple "first guess," then add layers of corrections to account for the chaos of the quantum world. The first layer is easy, but the second and third layers are like trying to solve a Rubik's cube while blindfolded and spinning on a merry-go-round. This is where "two-loop" calculations come in. Think of a "loop" as a detour a particle takes inside the collision—a temporary ghost particle popping in and out of existence. Calculating these detours is essential because they change the final result. Without them, our predictions are like a weather forecast that ignores the wind; they might be close, but they won't be accurate enough to spot the tiny, hidden signals of new physics hiding in the noise.

This paper tackles one of the messiest, most difficult collisions in the entire LEGO set: the creation of a top quark, an anti-top quark, and a W boson all at once. The top quark is the heavyweight champion of the particle world, and the W boson is a force-carrier that makes things decay. When they are born together, the math becomes a tangled knot of equations involving heavy masses and complex shapes that have never been fully untangled before.

The Paper's Mission: Untying the Knot

In this work, the authors present a new way to solve this mathematical knot. They have successfully computed the "two-loop" corrections for the production of a top-antitop pair and a W boson ( $t\bar{t}W$ ) at the Large Hadron Collider (LHC). This is a massive achievement because the complexity of this specific collision is so high that a traditional, fully written-out math solution (an analytic formula) is currently impossible to find.

Instead of trying to write down the entire solution on a single sheet of paper, the authors used a clever "hybrid" strategy. Imagine trying to describe a complex 3D sculpture. You could try to write a paragraph describing every curve, or you could take a photo of it from every angle and then use a computer to reconstruct the shape. The authors did something similar. They broke the problem down into two parts:

The Structure: They identified the "special functions" (the mathematical shapes) that the solution must be made of. These functions are like the specific curves of the sculpture.
The Numbers: Instead of calculating the exact numbers for every curve by hand, they used a technique called "finite-field evaluation." Think of this as checking the sculpture's shape at thousands of specific points using a digital ruler that never makes a rounding error. By checking these points on a grid of rational numbers, they could reconstruct the exact values of the coefficients (the numbers) without ever needing the full, messy formula.

What They Found and How They Did It

The team calculated the "finite remainder" of the collision. In the world of quantum physics, calculations often blow up with infinite numbers (singularities) that don't make physical sense. The authors showed that their method allows these infinities to cancel out perfectly, leaving behind a clean, finite number that describes the actual probability of the collision happening.

They didn't just do this once; they did it across a massive grid of 224,640 different scenarios (phase-space points). This grid covers all the different ways the particles could fly apart after the crash. To ensure they didn't make a mistake, they built a second, completely different computer program to check their work. This second program used a different mathematical framework (conventional dimensional regularisation) and a different way of solving the equations. When they compared the results from both programs, they matched to five significant digits. This cross-check confirms that their numbers are solid and reliable.

Why It Matters

The results of this paper are already being used to improve the theoretical predictions for the $t\bar{t}W$ process. Currently, experiments at the LHC (like ATLAS and CMS) have measured the rate of these collisions and found them to be slightly higher than the Standard Model predicted. This paper provides the precise "ruler" needed to measure that difference accurately. If the difference remains after applying these new, ultra-precise calculations, it could be a sign of new physics beyond our current understanding. If the difference disappears, it means our Standard Model is even more robust than we thought.

The authors emphasize that while they haven't found a single, beautiful formula that solves the whole problem, their method of combining numerical precision with algebraic control is a powerful new tool. It allows them to handle the "heavy" masses of the top quark and W boson exactly, without making shortcuts that could hide the truth. This approach opens the door to solving other complex collisions that were previously thought to be too difficult to calculate, bringing us one step closer to understanding the fundamental rules of our universe.

Technical Summary: Two-loop QCD amplitudes for $t\bar{t}W$ production at the LHC in the leading-colour approximation

Problem Statement
This work addresses the computation of two-loop Quantum Chromodynamics (QCD) scattering amplitudes for the production of a top–antitop quark pair in association with a $W$ boson ( $t\bar{t}W$ ) at the Large Hadron Collider (LHC). This process is phenomenologically critical due to its distinctive multilepton signatures, its role as an irreducible background for Higgs and four-top production, and the persistent tension between experimental measurements and Standard Model predictions.

Theoretical predictions for this process at next-to-next-to-leading order (NNLO) are hindered by significant technical challenges. The process involves a $2 \to 3$ kinematic configuration with multiple independent mass scales (top-quark mass $m_t$ and $W$ -boson mass $m_W$ ) and massive internal particles. These features lead to:

Algebraic Complexity: A large number of independent kinematic variables and intricate rational coefficients.
Analytic Complexity: The presence of massive internal lines introduces higher-genus geometries (elliptic curves) and nested square roots, which prevent the use of standard analytic methods developed for massless processes.
While partial results exist for related processes like $t\bar{t}H$ and $t\bar{t}j$ , a fully analytic two-loop calculation for $t\bar{t}W$ with exact mass dependence remains out of reach with current techniques.

Methodology
The authors adopt a hybrid computational strategy that combines numerical evaluation with rigorous algebraic and analytic control, avoiding a fully analytic reconstruction of the final result while retaining exact treatment of singularities. The calculation is performed in the 't Hooft-Veltman (tHV) scheme, with an independent cross-check in the Conventional Dimensional Regularisation (CDR) scheme.

Amplitude Decomposition and Projection: The bare amplitude is decomposed into a basis of 24 Lorentz tensor structures. Using physical projectors, the authors extract scalar form factors. The integrand generation involves 722 Feynman diagrams, reduced to 210 in the generalised leading-colour approximation (LCA).
Integral Reduction: The scalar loop integrals are reduced to a set of 330 Master Integrals (MIs) using Integration-by-Parts (IBP) relations. This reduction is performed numerically over finite fields using the FiniteFlow framework to handle the high algebraic complexity.
Special Function Representation: Instead of solving the differential equations (DEs) for the MIs in canonical form (which is difficult due to elliptic curves), the authors express the $\epsilon$ $ϵ$ -expansion of the MIs in terms of a set of "special functions." These functions are constructed to satisfy specific conditions that allow for the exact cancellation of Ultraviolet (UV) and Infrared (IR) poles.
- Most functions are algebraically independent iterated integrals.
- A small subset encoding the most complex analytic structures (elliptic curves) is treated via numerical $\mathbb{Q}$ -linear independence checks using the PSLQ algorithm.
Rational Coefficient Reconstruction: Rather than deriving analytic expressions for the rational coefficients multiplying these special functions, the authors reconstruct their exact rational values point-by-point. This is achieved by evaluating the amplitudes over finite fields ( $\mathbb{Z}_p$ ) on a grid of rationalised phase-space points and applying the Chinese Remainder Theorem and Wang's rational reconstruction algorithm.
Numerical Evaluation: The special functions are evaluated by numerically solving their systems of differential equations using power-series expansions (via AMFlow). The authors employ a "replica method" to estimate numerical uncertainties, ensuring high precision across the phase space.
Cross-Check: A second, independent calculation is performed in the CDR scheme. This approach bypasses tensor projectors, fixes masses and the dimensional regulator $\epsilon$ to numerical values early in the process, and uses a block-triangular system for integral reduction solved with floating-point arithmetic.

Key Results

Finite Remainders: The authors successfully computed the two-loop finite remainder for the $t\bar{t}W$ process in the leading-colour approximation, retaining exact dependence on $m_t$ and $m_W$ .
Hard Functions: The colour- and polarisation-summed hard functions ( $H^{(0)}$ , $H^{(1)}$ , $H^{(2)}$ ) were evaluated on a grid of 224,640 phase-space points. These values are provided in ancillary files and benchmarked in the paper.
Validation:
- The analytic cancellation of UV and IR poles was verified, confirming consistency with the universal structure of singularities.
- The one-loop results were cross-checked against existing literature.
- The two-loop results were cross-checked against the independent CDR calculation, showing agreement within five significant digits.
Efficiency: The hybrid approach reduced the number of required prime fields for rational reconstruction by approximately 30% compared to reconstructing the $\epsilon$ -expanded coefficients of the MIs directly.

Significance and Claims
The paper claims that this work provides the first exact two-loop QCD amplitude for $t\bar{t}W$ production retaining full kinematic dependence on massive particles, a result previously unattainable via fully analytic methods.

The significance of the work lies in:

Enabling NNLO Predictions: The computed hard functions have already been used to obtain the double-virtual contribution to the NNLO QCD cross section for $t\bar{t}W$ production, as referenced in the paper.
Methodological Advancement: The study demonstrates that a hybrid strategy—combining finite-field rational reconstruction with special-function numerical evaluation—can successfully tackle multi-scale, massive two-loop amplitudes where fully analytic solutions are currently impossible.
Robustness: The agreement between two distinct computational frameworks (tHV with special functions vs. CDR with floating-point reduction) provides a stringent validation of the results.

The authors conclude that this hybrid approach, potentially combined with improved interpolation techniques, offers a viable path toward NNLO predictions for a range of complex scattering processes that remain beyond the reach of traditional analytic calculations.

More like this