Analytic next-to-leading order electroweak corrections to Higgs boson pair production at high energies

This paper presents a complete analytic calculation of next-to-leading order electroweak corrections to gluon-induced Higgs boson pair production in the high-energy limit, demonstrating that these corrections reduce the cross-section by approximately 10% and providing precise numerical results down to low transverse momenta.

The Gist

Imagine the Large Hadron Collider (LHC) as the world's most powerful particle-smashing machine. Its main goal right now is to study the Higgs boson, the particle that gives other particles their mass. But scientists don't just want to see one Higgs boson; they want to see two of them created at the same time (Higgs pair production). This is the "Holy Grail" because it allows them to measure how Higgs bosons interact with themselves, which is crucial for understanding the stability of our universe.

However, predicting exactly how often this happens is incredibly difficult. It's like trying to predict the exact path of a pinball bouncing through a machine filled with millions of moving flippers, springs, and magnets.

Here is what this paper does, explained simply:

1. The Problem: The "Noise" of the Universe

When two protons smash together at the LHC, they usually create Higgs pairs by fusing two gluons (particles that hold atomic nuclei together). This process is messy.

• The "Main Event": The basic calculation (like a simple sketch of the pinball path) is easy.
• The "Corrections": In reality, quantum physics is weird. Virtual particles pop in and out of existence, creating "loops" of activity that change the outcome. The authors of this paper are calculating the Next-to-Leading Order (NLO) Electroweak corrections.
  • Analogy: Imagine you are baking a cake. The "Leading Order" is the recipe (flour, eggs, sugar). The "Electroweak Corrections" are the subtle effects of humidity, the exact temperature of your oven, and the altitude of your kitchen. They are small, but if you want a perfect cake (a precise prediction for the LHC), you must account for them.

2. The Challenge: Too Many Variables

The paper focuses on the "High Energy" limit. This means the particles are moving so fast that their energy is huge compared to their mass.

• The Difficulty: The math involves many different "scales" (masses of the top quark, the W and Z bosons, the Higgs, etc.). Trying to solve the exact equation for all these variables at once is like trying to solve a Rubik's cube that is also on fire and spinning in a hurricane. The equations become so complex that computers can't solve them exactly for every possible scenario.

3. The Solution: The "Zoom Lens" and "Map"

Instead of trying to solve the whole impossible puzzle at once, the authors used a clever strategy: Expansion.

• The Analogy: Imagine you are trying to describe a massive, complex city to someone who has never seen it. Instead of listing every single brick in every building, you create a map.
  • First, you zoom out to see the major highways (the high-energy limit).
  • Then, you zoom in slightly to see the neighborhoods (the mass differences).
  • The authors calculated about 100 different "terms" in this expansion. Think of these as layers of detail on the map.
  • They didn't just get a rough sketch; they got a high-definition, analytic map that works even when you zoom in on specific details (like the transverse momentum of the Higgs boson).

4. The "Magic Trick": Padé Approximation

Even with 100 terms, the math is still a giant, unwieldy beast. To make it usable for real-world predictions, they used a mathematical technique called Padé approximation.

• The Analogy: Imagine you are trying to guess the shape of a curve based on a few dots. You could draw a straight line (too simple), or a wiggly snake (too complex). A Padé approximation is like a smart algorithm that looks at your dots and draws the smoothest, most accurate curve that fits them perfectly, even in areas where you don't have any dots yet.
• This allowed them to turn their massive list of 100 terms into a clean, usable formula that works for a wide range of energies.

5. The Big Discovery: The "Negative Surprise"

After all the hard work, what did they find?

• They discovered that these "electroweak corrections" (the humidity and oven temperature of our cake analogy) are negative.
• Specifically, they reduce the probability of Higgs pair production by about 10% at high energies.
• Why this matters: If scientists didn't know this, they might think the LHC is producing more Higgs pairs than it actually is. By subtracting this 10%, they can get a true picture of what's happening. It's like realizing your scale is off by 10% and adjusting your diet plan accordingly.

Summary

This paper is a massive computational achievement. The authors built a high-precision mathematical map for a very complex particle physics process. They showed that when particles move at incredibly high speeds, the "electroweak" forces act like a brake, reducing the production of Higgs pairs by 10%.

This result is vital because it gives experimentalists at the LHC a more accurate "target" to aim for. Without this map, they might be looking for a treasure in the wrong place, or misinterpreting the size of the treasure they find. Now, they have a much clearer view of the Higgs boson's self-interaction, bringing us one step closer to understanding the fundamental rules of our universe.

Technical Summary

Here is a detailed technical summary of the paper "Analytic next-to-leading order electroweak corrections to Higgs boson pair production at high energies."

1. Problem Statement

The production of Higgs boson pairs ($gg \to HH$) via gluon fusion is the primary channel for probing the Higgs self-coupling at the Large Hadron Collider (LHC). Precise theoretical predictions are essential for interpreting experimental data, particularly in the high-luminosity phase.

• The Gap: While Next-to-Leading Order (NLO) Quantum Chromodynamics (QCD) corrections are well-established, the complete NLO Electroweak (EW) corrections involving top quarks, gauge bosons ($W, Z$), and Higgs bosons in the high-energy regime have been largely inaccessible analytically. Previous studies relied on numerical approaches or expansions in specific limits (e.g., large top mass), which fail to cover the full phase space, especially for boosted Higgs bosons.
• The Challenge: The calculation involves multiple mass scales ($m_t, m_W, m_Z, m_H$) and complex two-loop Feynman diagrams (both 1PI and 1PR), making a full analytic solution extremely difficult.

2. Methodology

The authors compute the complete NLO EW corrections to the form factors entering the $gg \to HH$ amplitude using a deep high-energy expansion approach.

A. Theoretical Framework

• Process: $g(q_1)g(q_2) \to H(q_3)H(q_4)$.
• Amplitude Decomposition: The amplitude is decomposed into two Lorentz structures with scalar coefficients $M_1$ and $M_2$, further split into "triangle" ($F_{tri}$) and "box" ($F_{box}$) form factors.
• Renormalization: The calculation uses the on-shell renormalization scheme.
  • Parameters renormalized: $W, Z, H$ masses, top quark mass ($m_t$), and electric charge.
  • Tadpoles: Tadpole contributions are consistently included (Fleischer–Jegerlehner scheme) to ensure finiteness.
  • Scheme: Results are transformed into the $G_\mu$ scheme (replacing $\alpha$ with the Fermi constant $G_F$).

B. Computational Strategy

Due to the complexity of multiple mass scales, the authors employ a nested expansion strategy:

1. Diagram Generation: Using qgraf and tapir, they generate amplitudes and translate them to FORM notation.
2. Mass Expansions:
   • Top Mass Dominance: They expand around the top quark mass ($m_t$), treating $m_t$ as the dominant scale.
   • Mass Differences: They expand in small parameters $\delta_X = 1 - m_X/m_t$ (where $X = H, Z, W$).
   • Higgs Mass: They perform a Taylor expansion in the external Higgs mass ($m_H$).
3. Integral Reduction:
   • The problem is reduced to a set of Master Integrals (MIs).
   • Families: 168 fully massive MIs (28 non-planar) and 161 QCD-like MIs (30 non-planar).
   • Tools: exp for mapping, LiteRed for integrand expansion, Kira and FIRE for reduction to MIs.
4. High-Energy Limit: The MIs are expanded in the limit where Mandelstam variables $s, t \gg m_t^2$. The authors compute approximately 100 expansion terms analytically.
5. Approximation Construction:
   • Padé Approximants: To recover the full dependence on $m_t$ from the truncated expansion, they construct Padé approximants $[n/m](x)$.
   • Sum Acceleration: They use Aitken extrapolation to accelerate the convergence of the $\delta$ and $m_H$ expansions.

3. Key Contributions

• First Complete Analytic Result: This work provides the first full analytic result for NLO EW corrections to $gg \to HH$ involving top quarks, gauge bosons, and Higgs bosons in the high-energy limit.
• Deep Expansion: The authors computed roughly 100 expansion terms, including terms up to $m_t^{108}$, $(m_H^{ext})^4$, and $\delta_X^4$.
• Universal Applicability: The analytic form factors allow for flexible numerical evaluation across a wide phase space, covering energies from the high-energy limit down to relatively low Higgs transverse momenta ($p_T \gtrsim 350$ GeV).
• Public Release: The results are intended for integration into the ggxy framework for computing hadronic cross-sections.

4. Results

• Magnitude of Corrections: The NLO EW corrections at the partonic level are found to be negative, with a magnitude of approximately $-10\%$ in the high-energy limit. This aligns with previous numerical findings but provides a rigorous analytic confirmation.
• Structure of Form Factors:
  • Leading Logarithms: The box form factors exhibit specific logarithmic structures. $F_{box1}$ is dominated by quadratic logarithms ($\ln^2(m_t^2/s)$), while $F_{box2}$ contains cubic logarithms ($\ln^3(m_t^2/s)$).
  • Triangle vs. Box: The triangle form factor is suppressed at high energies by a factor of $m_H^2/s$ due to the $s$-channel propagator.
• Convergence and Uncertainty:
  • The expansion in $\delta$ (mass differences) and $m_H$ converges rapidly.
  • The authors estimate a truncation uncertainty of $\pm 1\%$ for the squared NLO matrix element, which is significantly smaller than the scale uncertainty from NLO QCD corrections.
  • Padé approximants constructed with expansion depths up to $m_t^{108}$ provide stable results, whereas lower depths ($m_t^{32}$) showed instability, particularly in the imaginary parts.

5. Significance

• Phenomenological Impact: These corrections are crucial for precision measurements of the Higgs self-coupling at the High-Luminosity LHC (HL-LHC). A $-10\%$ correction is substantial and must be accounted for to avoid biases in extracting the trilinear Higgs coupling.
• Theoretical Benchmark: The analytic results serve as a critical cross-check for purely numerical approaches (e.g., sector decomposition or Monte Carlo integration) which often struggle with the high-energy regime.
• Future Applications: The methodology and results pave the way for combining high-energy expansions with forward-scattering expansions to cover the entire phase space for EW corrections, mirroring the success achieved previously for QCD corrections.
• High-Energy Behavior: The paper elucidates the high-energy behavior of scattering amplitudes in the full electroweak sector, demonstrating that despite the complexity of the full SM, the high-energy limit yields compact, predictive analytic structures dominated by specific logarithmic terms.

In summary, this paper represents a major advancement in precision Higgs physics, providing the necessary analytic tools to interpret future high-energy LHC data with the required theoretical accuracy.