Top-Yukawa contributions to $pp\to b\bar{b}H$: two-loop leading-colour amplitudes

This paper presents the derivation of two-loop leading-colour scattering amplitudes for bottom-quark pair production in association with a Higgs boson ($pp \to b\bar{b}H$) at the LHC, focusing on top-Yukawa contributions within the massless bottom and heavy-top approximations, and expresses the finite remainder using one-mass pentagon functions with analytically reconstructed rational coefficients.

The Gist

Imagine the Large Hadron Collider (LHC) as the world's most powerful, high-speed particle collider, smashing protons together to recreate the conditions of the early universe. When these protons collide, they sometimes produce a Higgs boson (the particle that gives other particles mass) along with a pair of bottom quarks (heavy cousins of the electron).

This specific event, called $pp \to b\bar{b}H$, is like finding a needle in a haystack. It's rare, messy, and incredibly hard to spot because the "haystack" (background noise from other collisions) is enormous.

This paper is a massive mathematical achievement that helps physicists predict exactly how often this needle should appear, so they can spot it when it actually happens. Here is the breakdown of what the authors did, using simple analogies.

1. The Problem: The "Heavy Top" Ghost

To understand the collision, physicists have to calculate the probability of it happening. This involves looking at the forces at play.

• The Cast: The collision involves bottom quarks and the Higgs. But there's a "ghost" player: the top quark. The top quark is the heaviest particle in the universe. It's so heavy it doesn't even exist as a free particle in the collision; it only appears as a fleeting, virtual loop in the math.
• The Challenge: The authors needed to calculate the "two-loop" contribution. In physics, a "loop" is a complex path a particle takes in a Feynman diagram (a drawing of particle interactions).
  • One-loop is like a simple detour.
  • Two-loops is like a detour that has a detour inside it. It's a tangled knot of math.
• The Difficulty: Calculating these two-loop paths for a process with five particles flying out (two bottom quarks, two gluons, and a Higgs) is like trying to solve a Rubik's Cube while riding a unicycle on a tightrope. The math gets so huge that standard computers can't handle it.

2. The Strategy: The "Heavy Top" Shortcut

The authors used a clever trick called the Heavy Top Limit (HTL).

• The Analogy: Imagine you are trying to calculate the sound of a giant drum (the top quark) hitting a tiny bell (the Higgs). Because the drum is so massive compared to the bell, you don't need to track every vibration of the drum skin. You can pretend the drum is a solid, immovable wall that just transfers energy instantly.
• The Result: This simplifies the math significantly. Instead of tracking the heavy top quark's mass, they treat the interaction as a local "effective" force. This turns a nightmare calculation into a manageable one.

3. The Method: The "Finite Field" Puzzle

Even with the shortcut, the equations were still too complex to solve with standard algebra. The authors used a technique called Finite-Field Arithmetic.

• The Analogy: Imagine you are trying to guess a secret 10-digit phone number. If you try to write down the whole number at once, it's impossible. But, if you ask, "What is the last digit?" and "What is the second-to-last digit?" separately, it's easy.
• The Process:
  1. They broke the massive algebraic equation into tiny pieces.
  2. They solved these pieces using "finite fields" (a type of math where numbers wrap around, like a clock). This keeps the numbers small and prevents the computer from crashing.
  3. They solved the puzzle thousands of times with different "clocks" (different prime numbers).
  4. Finally, they used a reconstruction algorithm to piece all those small answers back together to reveal the giant, complex formula.

4. The Result: A New Map for the LHC

The authors successfully derived the two-loop scattering amplitudes.

• What is an amplitude? Think of it as a "probability map." It tells you exactly how likely a specific collision is to happen.
• The Output: They didn't just write a paper; they built a software library (a C++ tool) that other scientists can use.
• Why it matters:
  • Precision: Before this, predictions for this specific collision had huge "uncertainty margins" (like guessing the weather with a 50% chance of rain). This work reduces that uncertainty, allowing scientists to say, "We expect X amount of these collisions."
  • Discovery: If the LHC sees more or fewer collisions than this new, precise map predicts, it could mean there is new physics hiding there—perhaps a new particle or a new force we haven't discovered yet.
  • Higgs Self-Interaction: This process is crucial for studying how the Higgs boson interacts with itself (the "triple Higgs coupling"), which is key to understanding the stability of the universe.

Summary

In short, these physicists built a super-precise GPS for a very rare particle collision. They used a "heavy object shortcut" to simplify the terrain and a "modular puzzle-solving" technique to calculate the route. Now, experimentalists at the LHC have a much sharper map to find the needle in the haystack, bringing us one step closer to understanding the fundamental rules of our universe.

Technical Summary

1. Problem Statement

The production of a Higgs boson in association with a bottom-quark pair ($pp \to b\bar{b}H$) at the Large Hadron Collider (LHC) is a crucial process for probing the bottom-Yukawa coupling and studying the triple Higgs coupling. The inclusive cross-section for this process receives contributions from three terms proportional to the bottom-Yukawa ($y_b$), top-Yukawa ($y_t$), and their interference ($y_b y_t$).

While the $y_b$-induced contributions have been computed up to Next-to-Next-to-Leading Order (NNLO) in QCD, the $y_t$-induced contributions (where the Higgs couples to a closed top-quark loop) are currently known only at Next-to-Leading Order (NLO). The lack of NNLO precision for the $y_t$ component limits the theoretical accuracy of predictions, with current uncertainties remaining around 45%.

The primary challenge in computing the $y_t$ contribution at NNLO is the complexity of the two-loop five-particle scattering amplitudes ($0 \to b\bar{b}ggH$ and $0 \to b\bar{b}\bar{q}qH$) involving an external massive particle (the Higgs). Calculating these amplitudes with massive bottom quarks is extremely difficult. Therefore, the authors aim to compute the two-loop amplitudes in the Heavy-Top-Limit (HTL) effective theory framework, treating the bottom quark as massless. This result serves as the essential missing ingredient to construct the full NNLO prediction via a "massification" procedure (adding mass effects to the massless result).

2. Methodology

The authors employed a sophisticated computational framework combining modern multi-loop techniques:

• Effective Field Theory (EFT): The calculation is performed in the HTL framework where the top quark is integrated out, replacing the top-loop with a local effective operator ($H G_{\mu\nu}^a G^{a\mu\nu}$). The bottom quark is treated as massless.
• Leading-Colour Approximation: To make the two-loop calculation tractable, the authors work in the leading-colour limit ($N_c^2$ and $N_c n_f$ terms), neglecting sub-leading colour terms (expected to be $\sim 10\%$ corrections).
• Feynman Diagram Generation & Decomposition:
  • Diagrams were generated using QGRAF.
  • Colour decomposition was performed to isolate independent partial amplitudes.
  • Tensor reduction was achieved using a four-dimensional projection method to express amplitudes as linear combinations of scalar Feynman integrals.
• Integration-by-Parts (IBP) Reduction:
  • Scalar integrals were reduced to a basis of Master Integrals (MIs) using optimized IBP relations generated by NEATIBP and solved numerically over finite fields.
  • The computation involves non-planar integral families (specifically the $DPzz$ family), which significantly increases complexity compared to standard leading-colour calculations.
• Finite-Field Reconstruction:
  • Instead of manipulating massive symbolic expressions, the authors used finite-field arithmetic (via FINITEFLOW) to evaluate the amplitudes numerically over prime fields.
  • Rational coefficients multiplying the master integrals were reconstructed analytically from these numerical evaluations.
• Kinematic Parametrisation:
  • To ensure rational dependence on kinematics, momentum twistor variables were used to parametrize the external momenta for the five-particle scattering with one massive leg.
• Function Basis:
  • The Master Integrals were expressed in terms of one-mass pentagon functions.
  • The final finite remainders were expressed as linear combinations of monomials built from these pentagon functions and transcendental constants.

3. Key Contributions

• First Two-Loop Amplitudes for $y_t$-induced $b\bar{b}H$: The paper provides the first analytic derivation of the two-loop scattering amplitudes for $pp \to b\bar{b}H$ induced by the top-Yukawa coupling in the HTL framework.
• Handling Non-Planar Topologies: The computation successfully navigated the complexity of non-planar integral families (e.g., $DPzz$) which are typically prohibitive in leading-colour calculations, demonstrating the power of the finite-field reconstruction method.
• Analytic Results in Momentum Twistors: The finite remainders are provided in a compact analytic form using momentum twistor variables and one-mass pentagon functions, facilitating efficient numerical evaluation.
• Public Software Release: The authors released a comprehensive suite of tools, including:
  • MATHEMATICA ancillary files containing the full analytic expressions for finite remainders and pole terms.
  • A C++ library (using PENTAGONFUNCTIONS++ and QD for high precision) for numerical evaluation of hard functions across all partonic channels.

4. Results

• Amplitudes: Analytic expressions for the helicity-dependent leading-colour partial finite remainders were derived for the processes $0 \to b\bar{b}ggH$ and $0 \to b\bar{b}\bar{q}qH$.
• Validation:
  • The one-loop full-colour amplitudes were validated against OPENLOOPS.
  • The singular structure (UV and IR poles) was verified against universal subtraction terms.
  • Renormalization scale dependence was explicitly checked.
  • Ward identities were satisfied for the gluon polarization vectors.
• Numerical Benchmarks: The authors provided numerical benchmarks for the hard functions ($H^{(0)}, H^{(1)}, H^{(2)}$) at a specific physical phase-space point.
  • The two-loop corrections ($H^{(2)}$) were found to be significant (e.g., $\sim 13\%$ for $gg \to b\bar{b}H$ and $\sim 21\%$ for $q\bar{q} \to b\bar{b}H$ relative to the Born term).
• Numerical Stability: A stability analysis revealed that while most phase-space points are stable in double precision, a small fraction (approx. 10% in the $gg$ channel) requires higher precision (quadruple/64-bit) for rational coefficients, and a very small fraction (0.4%) requires higher precision for the pentagon functions. The $b\bar{b}$ channels showed different stability characteristics due to cancellations in the colour decomposition.

5. Significance

• Enabling NNLO Predictions: This work provides the critical two-loop ingredient required to compute the NNLO QCD cross-section for the top-Yukawa contribution to $b\bar{b}H$ production. When combined with existing $y_b$ results and massification techniques, this allows for a full NNLO prediction of the $pp \to b\bar{b}H$ process.
• Reducing Theoretical Uncertainty: Improving the precision of the $y_t$ component from NLO to NNLO is expected to significantly reduce the theoretical scale uncertainty (currently $\sim 45\%$), which is vital for precision Higgs physics and constraints on New Physics (e.g., MSSM).
• Foundation for Higgs+Jets: The scattering processes calculated here are also a subset of $H+2$-jet production in the HTL framework. The techniques and results pave the way for future calculations of $H+4$-gluon amplitudes and full-colour subleading corrections.
• Methodological Benchmark: The successful application of finite-field reconstruction to a two-loop, five-point, one-mass amplitude with non-planar topologies serves as a benchmark for future high-precision QCD calculations at the LHC.