⚛️ phenomenology

ggxy: Fast and flexible NLO QCD corrections to gluon-initiated processes

The paper introduces **ggxy**, a fast and flexible C++ library that calculates Next-to-Leading Order (NLO) QCD corrections for gluon-initiated top-mediated processes like $gg\to HH$ , $gg\to ZH$ , and $gg\to ZZ$ with full top-quark mass dependence, featuring analytical approximations for two-loop virtual amplitudes and an interface to Powheg for parton shower matching.

Original authors: Daniel Stremmer

Published 2026-07-08

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Daniel Stremmer

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 Large Hadron Collider (LHC) as a giant, high-speed particle accelerator that smashes protons together to see what happens. Scientists are getting better at measuring these crashes, but to understand the results, they need a very precise "map" of what the laws of physics predict should happen.

This paper introduces a new, super-fast computer tool called ggxy (pronounced like "giggle-zy" or just "gg-xy") that helps draw that map for specific, complicated crashes involving gluons (the particles that hold atomic nuclei together).

Here is the breakdown of what the authors did, using simple analogies:

1. The Problem: The "Heavy" Calculation

In the world of particle physics, some calculations are like trying to solve a Rubik's cube while blindfolded. Specifically, calculating what happens when two gluons smash together to create a Higgs boson and a Z boson (or two Z bosons) is incredibly hard.

The difficulty comes from the top quark, the heaviest known particle. Because it is so heavy, it acts like a massive anchor in the math. If you try to calculate the crash while ignoring the top quark's weight, the map is wrong. If you try to include its weight exactly, the math becomes so heavy and slow that it takes forever to run on a computer.

2. The Solution: The "Smart Shortcut"

The authors built ggxy, a software library that acts like a GPS with a "smart routing" feature. Instead of trying to calculate every single tiny detail of the crash from scratch (which is too slow), ggxy uses two different "shortcuts" (mathematical approximations) that work perfectly in different parts of the journey:

The "Forward" Shortcut: This works best when the particles are moving in a straight line (low energy).
The "High-Energy" Shortcut: This works best when the particles are moving very fast and far apart (high energy).

By stitching these two shortcuts together, the software can cover the entire journey (the whole "phase space") without ever getting stuck. It's like having a map that uses a city street guide for downtown and a highway guide for the open road, switching between them automatically so you never hit a traffic jam.

3. What Makes It Special?

Full Weight Included: Unlike older tools that had to guess the weight of the top quark, ggxy includes the full, exact weight of the top quark. This makes the predictions much more accurate.
Flexible Settings: You can tell the software to change the "rules" of the game. You can adjust the mass of the top quark or the Z boson to see how the results change. It's like a video game where you can tweak the physics engine to see how the characters behave under different conditions.
Speed and Stability: The authors say this tool is "fast and flexible." It can calculate complex results in about 90 minutes on a single computer core, which is roughly three times slower than a simpler calculation but still incredibly fast compared to doing the math the "hard way."

4. Real-World Application: The "POWHEG" Connection

The paper mentions that ggxy has been connected to another program called POWHEG. Think of ggxy as the architect who draws the perfect blueprint of the crash. POWHEG is the construction crew that takes that blueprint and builds a full simulation, adding in the "debris" and "aftermath" (called parton showers) that happen when the particles fly apart. This allows scientists to simulate the entire event from start to finish, matching what they actually see in the detectors.

5. What Did They Find?

The authors tested their tool against other known results and found they matched perfectly. They also discovered something interesting about the "mass scheme" (how we define the weight of the top quark):

Depending on how you define the top quark's weight, the predicted results can change by 30% to 40%.
This huge difference is mostly driven by the longitudinal polarization of the Z boson (a specific way the Z boson spins or vibrates). It's like realizing that the way a spinning top wobbles changes the entire outcome of a game.

Summary

In short, ggxy is a new, high-precision calculator for particle physicists. It uses clever mathematical shortcuts to quickly and accurately predict what happens when gluons collide to create heavy particles, taking the full weight of the top quark into account. This helps scientists prepare for the High-Luminosity LHC, ensuring they have the right maps to interpret the massive amounts of data they will soon collect.

Technical Summary of "ggxy: Fast and flexible NLO QCD corrections to gluon-initiated processes"

Problem Statement
The High-Luminosity LHC (HL-LHC) demands theoretical predictions with reduced uncertainties to match improved experimental precision. Specifically, calculating Next-to-Leading Order (NLO) QCD corrections for gluon-initiated processes involving top-quark loops, such as $gg \to ZH$ and $gg \to Z^*Z^*$ , presents a significant challenge. The primary bottleneck is the computation of two-loop virtual amplitudes with full top-quark mass dependence. While various expansion methods (large-mass, high-energy, small-mass, $p_T$ ) and numerical approaches have been applied to related processes, a public implementation offering full mass dependence for $gg \to ZH$ at NLO was previously unavailable.

Methodology
The authors address the computational complexity of the two-loop amplitudes by employing a hybrid expansion strategy that combines two distinct kinematic limits to cover the entire phase space:

Forward Expansion: An expansion in the limit where the Mandelstam variable $|t| \ll s, m_t^2$ . This is treated as a Taylor expansion.
High-Energy Expansion: An asymptotic expansion in the limit where $m_t^2 \ll s, |t|$ .

In both limits, the amplitudes are expanded assuming the external masses of the final-state particles ( $q_3^2$ and $q_4^2$ ) are small compared to other scales. This allows the amplitudes to be expressed as a sum of master integrals independent of external masses.

For the forward limit, the expansion is performed in small variables $t$ , $q_3^2$ , and $q_4^2$ .
For the high-energy limit, a Padé approximation is constructed to extend the region of convergence.
The authors also utilize an alternative method for the forward limit involving an expansion in a small variable $\delta q$ introduced via momentum conservation, allowing for deeper mass expansions at the cost of fewer terms in the $t$ expansion.

These analytical approximations are implemented in a new C++ library named ggxy. The library calculates:

One- and two-loop amplitudes.
Hadronic cross sections at NLO QCD.
Real corrections using the Catani-Seymour dipole subtraction scheme, with phase-space integration handled by avhlib.
One-loop 2 $\to$ 3 amplitudes using tools such as Recola, Collier, CutTools, and OneLOop.
Mass scheme conversions (on-shell to $\overline{\text{MS}}$ ) and running masses via CRunDec.

Key Contributions

Public Library (ggxy): The paper presents the first public implementation of $gg \to ZH$ at NLO QCD with full top-quark mass dependence. The library allows users to freely modify input parameters, including the top-quark mass, its renormalization scheme, and the masses of external $Z$ and Higgs bosons.
Process Coverage: The library handles $gg \to ZH$ , $gg \to ZZ$ , and by extension, processes with off-shell decays such as $gg \to \ell^-\ell^+ H$ and $gg \to \nu_\ell \bar{\nu}_\ell H$ .
POWHEG Interface: The results are interfaced with POWHEG, enabling the matching of these NLO calculations to parton showers (e.g., Pythia, Herwig) while preserving spin correlations and off-shell effects of the $Z$ decay.
Polarization Capabilities: The implementation supports the calculation of hadronic cross sections for transversely and longitudinally polarized $Z$ bosons.

Results

Convergence: Numerical tests (Figures 1 and 2) demonstrate that the forward ( $t \to 0$ ) and high-energy expansions agree well in the phase-space region around $p_T \approx 200$ GeV. The forward expansion is reliable at low $p_T$ , while the high-energy expansion is accurate at high $p_T$ . For $gg \to Z^*Z^*$ , both methods yield essentially identical results across the tested $\sqrt{s}$ range.
Validation: The integrated cross sections for $gg \to ZH$ at $\sqrt{s} = 13, 13.6, 14$ TeV agree with existing results from Ref. [35] within statistical uncertainties (Table 1).
Performance: Achieving a precision of 0.1%–0.2% on a single core requires approximately 90 minutes, making the $gg \to ZH$ calculation roughly three times slower than the $gg \to HH$ calculation.
Renormalization Scheme Dependence: The study reveals significant differences (up to 30%–40%) between predictions using the on-shell and $\overline{\text{MS}}$ top-quark mass renormalization schemes. This dependence is particularly pronounced for the longitudinal polarization of the $Z$ boson, which dominates the unpolarized cross section, whereas differences for transverse polarization decrease in the high-energy region.

Significance
The paper claims that ggxy provides a fast, stable, and flexible tool for evaluating two-loop amplitudes and hadronic cross sections for gluon-initiated processes with full top-quark mass effects. By combining analytical approximations in complementary kinematic limits, the authors circumvent the computational bottleneck of exact two-loop calculations. The interface with POWHEG is highlighted as a crucial feature, allowing for the generation of fully exclusive events with parton shower matching for processes like $gg \to ZH$ and $gg \to ZZ$ , which is essential for phenomenological studies at the HL-LHC. The authors note that future versions will extend these capabilities to polarized $gg \to ZH$ production and $gg \to VV$ amplitudes.