All planar three-loop Feynman integrals for the production of two vector bosons at hadron colliders

This paper computes all planar three-loop master integrals required for leading-colour N3LO QCD corrections to the production of two massive or off-shell vector bosons at hadron colliders by constructing pure integral bases, deriving canonical differential equations via finite field techniques, and solving them using generalised power series expansions.

The Gist

Imagine you are trying to predict the exact outcome of a massive, chaotic traffic jam, but instead of cars, you are tracking subatomic particles colliding at nearly the speed of light inside the Large Hadron Collider (LHC).

This paper is essentially a master map created by two physicists, Dhimiter Canko and Mattia Pozzoli, to help us understand what happens when two heavy "vehicles" (vector bosons) are created from the collision of two invisible "drivers" (quarks or gluons).

Here is the breakdown of their work using simple analogies:

1. The Goal: Predicting the Future with Extreme Precision

The LHC is like a high-speed camera taking pictures of the universe. We are now in an era where we don't just want to see that a crash happened; we want to know the exact angle, speed, and debris of every single piece.

To do this, scientists need to calculate the "traffic rules" of the universe (Quantum Chromodynamics, or QCD) with incredible accuracy. They are trying to reach a level of precision called N3LO (Next-to-Next-to-Next-to-Leading Order). Think of it like this:

• Leading Order: A rough sketch of the crash.
• N3LO: A 3D, slow-motion, forensic reconstruction that accounts for every tiny vibration and ripple in the fabric of space-time.

To get this level of detail, you have to calculate complex mathematical shapes called Feynman Integrals. These are the "mathematical blueprints" of the particle collisions.

2. The Problem: The "Three-Loop" Maze

Calculating these blueprints is like trying to solve a maze.

• One-loop: A simple maze with one path. Easy.
• Two-loop: A maze with a few dead ends. Harder, but solvable.
• Three-loop: A massive, multi-dimensional labyrinth with thousands of twists, turns, and hidden traps.

The authors had to solve all the possible "planar" (flat, non-tangled) versions of these three-loop mazes for a specific type of crash: creating two heavy particles. There are nine different types of mazes (families) they had to map out.

3. The Solution: The "Pure" Compass

The biggest challenge with these mazes is that the math gets incredibly messy, filled with confusing variables and "noise" (mathematical divergences).

The authors' breakthrough was finding a "Pure Basis."

• Analogy: Imagine trying to navigate a stormy ocean. The waves (mathematical noise) are huge and chaotic. The authors didn't just try to sail through the waves; they built a special compass (a "pure basis") that cuts through the noise.
• This compass allows them to write down the rules of the maze in a "Canonical" form. This is like translating a chaotic, handwritten note into a clean, standardized instruction manual. It separates the "size" of the problem from the "shape" of the problem, making it solvable.

4. The Alphabet of the Universe

Once they had the clean instructions, they needed to know what "letters" the universe uses to write these rules.

• In math, these are called an Alphabet.
• For simpler problems (two loops), the alphabet was small and stable.
• For these complex three-loop problems, the authors discovered new letters and new square roots that had never been seen before in this specific context.
• Metaphor: It's like discovering that to write a new, complex language, you need not just the standard A-Z, but also some special symbols and accents that didn't exist in the previous version of the language. They had to invent a dictionary for these new symbols.

5. The Calculation: The "Digital River"

Finally, they had to actually calculate the numbers.

• Because the math was so complex (involving those new square roots), they couldn't just write down a single, perfect formula.
• Instead, they used a technique called Generalised Power Series.
• Analogy: Imagine you need to cross a wide river. You can't jump it in one go. Instead, you build a series of stepping stones.
  1. They start at a safe, known point (a "boundary").
  2. They take a tiny step forward, calculating the exact height of the next stone.
  3. They repeat this thousands of times, moving step-by-step across the river of data until they reach the other side (the final answer).
• They did this using powerful computers and special software to ensure the numbers were accurate to 16 decimal places.

Why Does This Matter?

This paper is a foundation stone.

• The authors didn't just solve the math for fun; they built the engine needed to predict what the LHC will see in the future.
• Without these calculations, we can't tell if the Standard Model (our current best theory of physics) is perfect or if there are cracks in it that hint at New Physics (like Dark Matter or extra dimensions).
• By providing these "maps" and "compasses," they allow other scientists to finally calculate the full, ultra-precise predictions for particle collisions.

In short: These two physicists built the ultimate navigation system for the most complex particle collisions we can currently imagine, ensuring that when the next big discovery happens at the LHC, we will have the math ready to understand it.

Technical Summary

Here is a detailed technical summary of the paper "All planar three-loop Feynman integrals for the production of two vector bosons at hadron colliders" by Canko and Pozzoli.

1. Problem Statement

The paper addresses a critical bottleneck in high-energy physics: the computation of Next-to-Next-to-Next-to-Leading Order (N3LO) Quantum Chromodynamics (QCD) corrections for the production of two vector bosons ($pp \to V_1 V_2$, where $V \in \{W, Z, \gamma^*\}$) at the Large Hadron Collider (LHC).

• Context: As the LHC moves toward the High-Luminosity phase, experimental uncertainties are dropping to the sub-percent level. Theoretical predictions must match this precision to test the Standard Model and search for new physics.
• Specific Challenge: While massless three-loop integrals and single-mass cases have been solved, the planar three-loop four-point integrals with two external massive legs (representing the vector bosons) were incomplete. Specifically, the "pure" master integrals (MIs) required for the generalized leading-colour approximation of the virtual corrections were missing for several topologies and mass configurations (equal-mass and different-mass).

2. Methodology

The authors employed a rigorous, state-of-the-art computational framework combining algebraic geometry, finite-field arithmetic, and differential equations.

• Integral Families: The authors identified that all relevant planar three-loop integrals with two external massive legs and massless internal propagators fall into nine distinct integral families:

  • 2 Reducible Ladder-Box families (RL1, RL2).
  • 3 Ladder-Box families (PL1, PL2, PL3).
  • 4 Tennis-Court families (PT1, PT2, PT3, PT4).
  • These are grouped into two super-families ($F_{123}$ and $F_{132}$) related by momentum permutation.

• Master Integral Reduction:

  • Used Integration-by-Parts (IBP) relations to reduce the infinite set of Feynman integrals to a finite basis of Master Integrals (MIs).
  • Employed the FiniteFlow framework and finite-field techniques to handle the massive algebraic complexity of the IBP systems.
  • Count: The study identified 823 independent MIs for the general (different-mass) case and 523 for the equal-mass case.

• Canonical Differential Equations (DEs):

  • Constructed a basis of "pure" MIs such that the system of DEs with respect to kinematic invariants takes the canonical form (dlog form): $d\vec{I} = \epsilon \, d\tilde{A} \cdot \vec{I}$.
  • This form factorizes the dimensional regulator $\epsilon$, allowing for an iterative solution in terms of iterated integrals.
  • Alphabet Construction: Determined the "alphabet" (the set of logarithmic arguments, or "letters") for the DEs. This involved identifying rational letters and algebraic letters involving square roots.

• Numerical Evaluation:

  • Since the square roots could not be simultaneously rationalized, a fully analytic solution in terms of multiple polylogarithms was not pursued.
  • Instead, the DEs were solved numerically using generalized power series expansions.
  • Boundary conditions were computed using the AMFlow package (Auxiliary Mass Flow method) with high precision (60 digits).
    The solution was propagated to physical phase-space points using the DiffExp package.

3. Key Contributions & Results

• Completeness: The paper provides the complete set of planar three-loop master integrals for two-vector-boson production for both equal-mass and different-mass scenarios. This fills the gaps left by previous works [23, 24], re-evaluating previously known families and computing the missing ones.
• New Analytic Structures (Alphabet):
  • A significant finding is the instability of the alphabet with loop order in the planar sector. Unlike the one-mass case where the alphabet remained stable from two to three loops, this work reveals new algebraic letters and four new square roots ($r_1, r_2, r_3, r_4$) that were absent in two-loop calculations.
  • In the general mass case, the alphabet consists of 37 letters (22 rational, 15 algebraic). In the equal-mass case, it reduces to 22 letters (13 rational, 9 algebraic).
  • The presence of these new square roots prevents a simultaneous rationalization, necessitating the semi-numerical approach.
• Pure Basis Construction: The authors successfully constructed pure bases for all 9 families, including the most complex sector (a 15-MI sector in the $PT_1/PT_2$ families), utilizing a bottom-up approach and Magnus exponential methods to decouple differential equations.
• Validation:
  • Results were cross-checked against independent AMFlow evaluations at specific phase-space points, showing agreement up to 16 significant digits.
  • For the equal-mass ladder-box families, results were compared against existing analytic solutions from Ref. [23], confirming perfect agreement.
• Public Availability: The authors provide ancillary files containing:
  • Definitions of the pure MIs and the canonical DE matrices ($d\tilde{A}$).
  • Boundary values for all families.
  • Mathematica scripts (DiffExpRun.wl, AMFSolver.wl) to evaluate the integrals at any physical point.

4. Significance

• Enabling N3LO Predictions: This work is a prerequisite for computing the three-loop scattering amplitudes for vector-boson pair production in the generalized leading-colour approximation. These amplitudes constitute the purely virtual part of the N3LO QCD corrections, a necessary step for precision phenomenology at the HL-LHC.
• Theoretical Insight: The discovery of new square roots and letters in the planar sector challenges the assumption that planar integrals are structurally simpler or more stable than non-planar ones as loop order increases. It highlights the growing complexity of multi-scale, multi-mass problems in QCD.
• Methodological Benchmark: The successful application of finite-field techniques to construct pure bases for such a large system (823 MIs) and the robust numerical evaluation strategy serve as a blueprint for future high-order calculations in multi-particle processes.

In summary, this paper delivers the essential mathematical building blocks (master integrals) required to push the theoretical precision of diboson production to N3LO, while simultaneously advancing the understanding of the analytic structure of multi-loop Feynman integrals with multiple massive scales.