ggxy: Fast and flexible NLO QCD corrections to gluon-initiated processes
Reviewed by Pith T0 review T1 audit T2 compute T3 formal T4 reserved 2026-07-08 16:30 UTCglm-5.2pith:PLOQPGTVrecord.jsonopen to challenge →
The pith
First public code for gg→ZH at NLO QCD with full top-mass effects
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The key mechanism is the dual-expansion strategy for two-loop amplitudes: the high-energy expansion (asymptotic, valid at large transverse momentum) and the forward-limit or t→0 expansion (Taylor-type, valid at small transverse momentum) are computed analytically to sufficient depth in mass corrections, and their regions of validity overlap around pT ~ 200 GeV for gg→ZH, allowing a seamless switch between them. A Padé approximation extends the convergence region of the high-energy expansion. This combination replaces a full exact two-loop calculation — which is computationally prohibitive — with a fast, analytically grounded approximation that the paper validates against an independent NLO结果
What carries the argument
High-energy expansion (mt² ≪ s, |t|, asymptotic) and forward-limit expansion (|t| ≪ s, mt², Taylor-type) of two-loop master integrals, combined with a Padé approximation for the high-energy branch; Catani-Seymour dipole subtraction for real corrections; POWHEG interface for parton-shower matching; CRunDec for top-mass scheme conversion.
If this is right
- Experimental analyses of ZH production at the HL-LHC can now use NLO QCD predictions with full top-mass dependence and parton-shower matching from a public tool, reducing a source of theoretical systematic uncertainty.
- The 30–40% spread between on-shell and MS-bar top-mass renormalization schemes at NLO — dominated by longitudinal Z polarization — quantifies a leading parametric uncertainty that must be addressed when extracting top-mass-sensitive observables from ZH data.
- The two-loop amplitude results for gg→Z*Z* can be directly reused for gg→Z*γ* and gg→γ*γ*, extending the library's physics reach without new two-loop calculations.
- The dual-expansion strategy, validated here for ZH and ZZ, provides a template for other gluon-initiated processes where exact two-loop amplitudes remain out of reach.
Where Pith is reading between the lines
- The claim of full phase-space coverage rests on the overlap region being well-controlled; if future differential measurements probe kinematic corners where both expansions simultaneously degrade (e.g., intermediate pT with specific rapidity configurations not shown in the convergence plots), the approximation could develop uncontrolled errors not visible in integrated cross-section validation.
- The 30–40% mass-scheme dependence, while large, could itself serve as a diagnostic: if the NLO corrections were dominated by terms that resum cleanly, the scheme dependence would be smaller. The persistence of large scheme dependence suggests substantial higher-order terms that are not yet captured, implying that NNLO corrections (if computed) could shift predictions by amounts comparable to the c
- The factor-of-three runtime increase of ZH over HH suggests the ZH amplitude structure carries additional complexity (extra form factors, polarization channels) that scales unfavorably; this could limit the practicality of parameter scans in data-analysis pipelines unless further optimization is done.
Load-bearing premise
The two-loop amplitudes are approximations, not exact results. The claim that the high-energy and forward-limit expansions together cover the full phase space depends on their overlap region being numerically reliable and on no kinematic configuration existing where both expansions fail at once — a claim supported by convergence plots in specific projections and by agreement of integrated cross sections with one external reference, but not by exhaustive point-by-point testing
What would settle it
A differential observable (a specific bin in pT, rapidity, or invariant mass) at NLO where the ggxy prediction, using the dual-expansion two-loop amplitudes, disagrees with an independent exact two-loop calculation beyond the stated integration uncertainty — particularly in a kinematic region where the overlap between the two expansions is thin or where neither expansion is individually convergent.
Figures
read the original abstract
We present the C++ library ggxy, which can be used for the fast and flexible calculation of partonic and hadronic cross sections to gluon-initiated top-mediated processes such as $gg\to HH$, as well as the newly derived processes $gg\to ZH$ and $gg\to ZZ$, at NLO QCD including full top-quark mass dependence. The two-loop virtual amplitudes are implemented using analytical approximations in different kinematic regions, while all other parts of the calculation are exact. This implementation allows to freely modify all input parameters, such as the top-quark mass and its renormalization scheme and the masses of the external $Z$/Higgs bosons. In addition, ggxy has been interfaced to Powheg, which allows the matching to parton showers.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. These proceedings present the C++ library ggxy for computing NLO QCD corrections to gluon-initiated processes, focusing on the newly implemented gg→ZH process with full top-quark mass dependence. The two-loop virtual amplitudes are approximated using high-energy and forward-limit expansions that, combined, are claimed to cover the full phase space. All other parts of the calculation (one-loop amplitudes, real corrections, subtraction) are treated exactly. The implementation is validated against Ref. [35] at the integrated cross-section level (Table 1) and interfaced to POWHEG for parton-shower matching. The paper also presents results for gg→Z*Z* and studies the top-quark mass scheme dependence for polarized Z bosons.
Significance. The paper provides the first publicly available implementation of gg→ZH at NLO QCD with full top-quark mass dependence, which is a valuable contribution to the phenomenology program for the HL-LHC. The code is publicly available (gitlab.com/ggxy/ggxy-release) and interfaced to POWHEG, enabling parton-shower matching. The study of top-quark mass scheme dependence for polarized Z bosons (Figures 3–4) provides useful physics insight. The two-loop amplitudes rely on expansion methods rather than exact numerical evaluation, which is a legitimate and well-established approach in this field.
major comments (2)
- §3 and §4, Table 1: The validation against Ref. [35] shows agreement at the 0.1% level for integrated cross sections, but no differential-level comparison against an independent calculation is presented. The paper states that 'several cross-checks at the integrated and differential level' have been performed (§4), but only the integrated-level Table 1 is shown. A differential comparison (e.g., dσ/dpT,H or dσ/dMZH) against an exact calculation such as those of Refs. [20–23] would substantially strengthen the claim of full phase-space coverage, since integrated cross sections can mask compensating region-specific errors. This is particularly relevant because the central claim—that the HE and forward-limit expansions together cover the full phase space—depends on no kinematic region where both expansions fail simultaneously. Including at least one differential-level validation plot or table
- §3, Figures 1–2: The convergence of the two expansion methods is demonstrated at √s = 1000 GeV (Fig. 1) and pT = 150 GeV (Fig. 2), but no amplitude-level validation is shown in the near-threshold region 211 < √s < 345 GeV, where only the forward expansion is available (the HE expansion is limited to √s > 2mt ≈ 345 GeV, as acknowledged in §3). In this region, the mass-expansion parameters mZ²/s ~ 0.19 and mH²/s ~ 0.35 are not negligibly small, and |t| can be comparable to s. While the paper notes that 'such small values are never used for the high-energy expansion due to the lower limit on pT,' it would strengthen the manuscript to briefly address whether the forward expansion remains reliable in this near-threshold region, or to cite where this is documented.
minor comments (5)
- §4, Table 1: It is not stated whether Ref. [35] uses exact numerical two-loop amplitudes (as in Refs. [20–23]) or similar expansion methods. This distinction is relevant for assessing how independent the validation is; a brief clarification would help.
- §3: The paper does not quantify a systematic uncertainty from expansion truncation. A brief statement of the expected size of truncation errors, or a reference to where this is quantified in the companion papers, would be useful.
- §1: 'continuos' should be 'continuous'. §4: 'longituidinal' should be 'longitudinal'. §4: 'agree well within in the statistical uncertainties' has a duplicated 'in'.
- Figures 1–2 are taken from Ref. [3] and Figures 3–4 from Ref. [2]. While this is noted, the proceedings would benefit from at least one figure showing original validation results (e.g., a differential comparison) not already available in the companion papers.
- §2: The notation for the expansion orders (m{0,2,4}, t0, HE, etc.) in Figures 1–2 is compact but not fully self-explanatory; a one-line clarification of the labeling convention would improve self-containedness.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments. Both points are well-taken and we will address them in the revised manuscript.
read point-by-point responses
-
Referee: §3 and §4, Table 1: No differential-level comparison against an independent calculation is presented, despite the claim that 'several cross-checks at the integrated and differential level' have been performed. A differential comparison would strengthen the claim of full phase-space coverage.
Authors: The referee is correct that the proceedings only show the integrated-level comparison in Table 1, while the text in §4 states that differential-level cross-checks have been performed. We agree that this is a gap in the presentation. The differential-level validation against the exact numerical results of Ref. [22] (Chen et al., JHEP 08 (2022) 056) is documented in detail in our companion paper Ref. [2] (arXiv:2603.15762), where dσ/dpT,H and dσ/dMZH distributions are compared and agreement at the per-mille level is found across the full kinematic range. However, we acknowledge that this is not visible in the proceedings themselves. In the revised version, we will (i) add a sentence in §4 explicitly referencing where the differential validation can be found, and (ii) if space permits within the proceedings format, include a brief differential comparison plot or table. We note that the concern about compensating region-specific errors is directly addressed by these differential comparisons, which show no such cancellations. revision: yes
-
Referee: §3, Figures 1–2: No amplitude-level validation is shown in the near-threshold region 211 < √s < 345 GeV, where only the forward expansion is available. The mass-expansion parameters mZ²/s ~ 0.19 and mH²/s ~ 0.35 are not negligibly small, and |t| can be comparable to s. The reliability of the forward expansion in this region should be addressed.
Authors: This is a valid concern. We should clarify the situation more explicitly in the manuscript. The forward expansion is a Taylor expansion in the small variables t, q₃², and q₄² simultaneously. In the near-threshold region 211 < √s < 345 GeV, while the mass-expansion parameters mZ²/s and mH²/s are indeed not negligibly small, the key point is that the forward expansion is used in this region only at relatively small pT (since the switch to the HE expansion occurs at larger pT), which corresponds to small |t|/s. The expansion in |t| is the primary expansion parameter, and the mass corrections are included up to sufficiently high order (m⁸ for the external masses, as shown in Figure 1) to ensure convergence. We acknowledge, however, that this reasoning is not spelled out in the proceedings. In the revised version, we will add a brief discussion addressing the reliability of the forward expansion in the near-threshold region, including a reference to Ref. [3] (arXiv:2509.07072) where the convergence properties are studied in more detail. We will also note that the differential-level validation mentioned in our response to the first comment covers this kinematic region and provides an indirect but important check on the amplitude-level accuracy. revision: partial
Circularity Check
No significant circularity found; self-citations are to first-principles amplitude calculations, validated against an independent group.
full rationale
The paper's central claim — that ggxy provides NLO QCD corrections to gg→ZH with full top-mass dependence using high-energy and forward-limit expansions — rests on two pillars: (1) the two-loop amplitude calculations from Ref. [3] (co-authored by Stremmer) and (2) validation against Ref. [35] (an independent group, Campillo Aveleira et al.). Ref. [3] derives the amplitudes from first-principles QCD using systematic Taylor and asymptotic expansions in well-defined kinematic limits, not by fitting to data or to the target observables. The expansions are defined in terms of physical kinematic parameters (s, t, m_t, external masses), not in terms of the cross sections they aim to predict. The Padé approximation (following Ref. [16]) is a standard mathematical resummation technique, not an ansatz fitted to results. Table 1 validates integrated cross sections against Ref. [35], which is an independent collaboration, providing genuine external support. Figures 1–2 (taken from Ref. [3]) demonstrate internal consistency between two independently derived expansion methods in their overlap region — a legitimate cross-check, not circular reasoning. The self-citations to Refs. [1, 2, 3] are normal practice for a software/phenomenology paper building on prior technical work, and the cited calculations are parameter-free QCD computations with externally falsifiable predictions. No step in the derivation chain reduces to its own inputs by construction. The minor self-citation load (Ref. [3] for amplitudes, Refs. [1,2] for the code) is not circular because the amplitude calculations are independent first-principles results, not fitted quantities renamed as predictions.
Axiom & Free-Parameter Ledger
axioms (3)
- domain assumption The high-energy and forward-limit expansions of the two-loop amplitudes, when combined, cover the full physically relevant phase space with sufficient accuracy.
- standard math The Catani-Seymour dipole subtraction scheme correctly handles the infrared singularities in the real corrections for these processes.
- standard math The one-loop 2→3 amplitudes computed by Recola/Collier/CutTools/OneLOop are correct and stable.
Reference graph
Works this paper leans on
-
[1]
{\tt ggxy}: a flexible library to compute gluon-induced cross sections
J. Davies, K. Schönwald, M. Steinhauser and D. Stremmer,ggxy: A flexible library to compute gluon-induced cross sections,Comput. Phys. Commun.320(2026) 109933 [2506.04323]
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [2]
-
[3]
J.Davies,D.Grau,K.Schönwald,M.Steinhauser,D.StremmerandM.Vitti,Two-loopQCD correctionstoZH andoff-shellZbosonpairproduction ingluonfusion,JHEP03(2026)200 [2509.07072]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[4]
A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX
S. Alioli, P. Nason, C. Oleari and E. Re,A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX,JHEP06(2010) 043 [1002.2581]. 6 ggxy: Fast and flexible NLO QCD corrections to gluon-initiated processesDaniel Stremmer
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[5]
A Brief Introduction to PYTHIA 8.1
T. Sjostrand, S. Mrenna and P. Z. Skands,A Brief Introduction to PYTHIA 8.1,Comput. Phys. Commun.178(2008) 852 [0710.3820]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[6]
Bellm et al.,The Physics of Herwig 7,2512.16645
J. Bellm et al.,The Physics of Herwig 7,2512.16645
-
[7]
B. A. Kniehl,Associated Production of Higgs and Z Bosons From Gluon Fusion in Hadron Collisions,Phys. Rev. D42(1990) 2253
work page 1990
-
[8]
D. A. Dicus and C. Kao,Higgs Boson -𝑍0 Production From Gluon Fusion,Phys. Rev. D38 (1988) 1008
work page 1988
-
[9]
D. A. Dicus, C. Kao and W. W. Repko,Gluon Production of Gauge Bosons,Phys. Rev. D36 (1987) 1570
work page 1987
-
[10]
E. W. N. Glover and J. J. van der Bij,Z BOSON PAIR PRODUCTION VIA GLUON FUSION,Nucl. Phys. B321(1989) 561
work page 1989
-
[11]
Gluon-induced Higgs-strahlung at next-to-leading order QCD
L. Altenkamp, S. Dittmaier, R. V. Harlander, H. Rzehak and T. J. E. Zirke,Gluon-induced Higgs-strahlung at next-to-leading order QCD,JHEP02(2013) 078 [1211.5015]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[12]
Production of two Z-bosons in gluon fusion in the heavy top quark approximation
K. Melnikov and M. Dowling,Production of two Z-bosons in gluon fusion in the heavy top quark approximation,Phys. Lett. B744(2015) 43 [1503.01274]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[13]
On top quark mass effects to $gg\to ZH$ at NLO
A. Hasselhuhn, T. Luthe and M. Steinhauser,On top quark mass effects to𝑔𝑔→𝑍 𝐻at NLO,JHEP01(2017) 073 [1611.05881]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[14]
Virtual corrections to $gg\to ZH$ in the high-energy and large-$m_t$ limits
J. Davies, G. Mishima and M. Steinhauser,Virtual corrections to𝑔𝑔→𝑍 𝐻in the high-energy and large-𝑚𝑡 limits,JHEP03(2021) 034 [2011.12314]
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[15]
J. M. Campbell, R. K. Ellis, M. Czakon and S. Kirchner,Two loop correction to interference in𝑔𝑔→𝑍 𝑍,JHEP08(2016) 011 [1605.01380]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[16]
$gg\to ZZ$: analytic two-loop results for the low- and high-energy regions
J. Davies, G. Mishima, M. Steinhauser and D. Wellmann,𝑔𝑔→𝑍 𝑍: analytic two-loop results for the low- and high-energy regions,JHEP04(2020) 024 [2002.05558]
work page internal anchor Pith review Pith/arXiv arXiv 2020
-
[17]
G. Wang, X. Xu, Y. Xu and L. L. Yang,Next-to-leading order corrections for𝑔𝑔→𝑍 𝐻with top quark mass dependence,Phys. Lett. B829(2022) 137087 [2107.08206]
work page internal anchor Pith review Pith/arXiv arXiv 2022
-
[18]
Virtual corrections to $gg\to ZH$ via a transverse momentum expansion
L. Alasfar, G. Degrassi, P. P. Giardino, R. Gröber and M. Vitti,Virtual corrections to 𝑔𝑔→𝑍 𝐻via a transverse momentum expansion,JHEP05(2021) 168 [2103.06225]
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[19]
Virtual QCD corrections to $gg \to ZZ$: top-quark loops from a transverse-momentum expansion
G. Degrassi, R. Gröber and M. Vitti,Virtual QCD corrections to gg→ZZ: top-quark loops from a transverse-momentum expansion,JHEP07(2024) 244 [2404.15113]
work page internal anchor Pith review Pith/arXiv arXiv 2024
-
[20]
L. Chen, G. Heinrich, S. P. Jones, M. Kerner, J. Klappert and J. Schlenk,𝑍 𝐻production in gluon fusion: two-loop amplitudes with full top quark mass dependence,JHEP03(2021) 125 [2011.12325]
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[21]
Two-loop helicity amplitudes for $gg \to ZZ$ with full top-quark mass effects
B. Agarwal, S. P. Jones and A. von Manteuffel,Two-loop helicity amplitudes for𝑔𝑔→𝑍 𝑍 with full top-quark mass effects,JHEP05(2021) 256 [2011.15113]. 7 ggxy: Fast and flexible NLO QCD corrections to gluon-initiated processesDaniel Stremmer
work page internal anchor Pith review Pith/arXiv arXiv 2021
-
[22]
L. Chen, J. Davies, G. Heinrich, S. P. Jones, M. Kerner, G. Mishima et al.,ZH production in gluon fusion at NLO in QCD,JHEP08(2022) 056 [2204.05225]
work page internal anchor Pith review Pith/arXiv arXiv 2022
-
[23]
Complete Next-to-Leading Order QCD Corrections to $ZZ$ Production in Gluon Fusion
B. Agarwal, S. Jones, M. Kerner and A. von Manteuffel,Complete Next-to-Leading Order QCD Corrections to ZZ Production in Gluon Fusion,Phys. Rev. Lett.134(2025) 031901 [2404.05684]
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[24]
The two-loop helicity amplitudes for $gg \to V_1 V_2 \to 4~\mathrm{leptons}$
A. von Manteuffel and L. Tancredi,The two-loop helicity amplitudes for 𝑔𝑔→𝑉 1𝑉2 →4 leptons,JHEP06(2015) 197 [1503.08835]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[25]
Two-loop helicity amplitudes for the production of two off-shell electroweak bosons in gluon fusion
F. Caola, J. M. Henn, K. Melnikov, A. V. Smirnov and V. A. Smirnov,Two-loop helicity amplitudes for the production of two off-shell electroweak bosons in gluon fusion,JHEP06 (2015) 129 [1503.08759]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[26]
RECOLA: REcursive Computation of One-Loop Amplitudes
S. Actis, A. Denner, L. Hofer, J.-N. Lang, A. Scharf and S. Uccirati,RECOLA: REcursive Computation of One-Loop Amplitudes,Comput. Phys. Commun.214(2017) 140 [1605.01090]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[27]
Collier: a fortran-based Complex One-Loop LIbrary in Extended Regularizations
A. Denner, S. Dittmaier and L. Hofer,Collier: a fortran-based Complex One-Loop LIbrary in Extended Regularizations,Comput. Phys. Commun.212(2017) 220 [1604.06792]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[28]
CutTools: a program implementing the OPP reduction method to compute one-loop amplitudes
G. Ossola, C. G. Papadopoulos and R. Pittau,CutTools: A Program implementing the OPP reduction method to compute one-loop amplitudes,JHEP03(2008) 042 [0711.3596]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[29]
OneLOop: for the evaluation of one-loop scalar functions
A. van Hameren,OneLOop: For the evaluation of one-loop scalar functions,Comput. Phys. Commun.182(2011) 2427 [1007.4716]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[30]
A General Algorithm for Calculating Jet Cross Sections in NLO QCD
S. Catani and M. H. Seymour,A General algorithm for calculating jet cross-sections in NLO QCD,Nucl. Phys. B485(1997) 291 [hep-ph/9605323]
work page internal anchor Pith review Pith/arXiv arXiv 1997
-
[31]
PARNI for importance sampling and density estimation
A. van Hameren,PARNI for importance sampling and density estimation,Acta Phys. Polon. B40(2009) 259 [0710.2448]
work page internal anchor Pith review Pith/arXiv arXiv 2009
-
[32]
Kaleu: a general-purpose parton-level phase space generator
A. van Hameren,Kaleu: A General-Purpose Parton-Level Phase Space Generator, 1003.4953
work page internal anchor Pith review Pith/arXiv arXiv
-
[33]
Version 3 of {\tt RunDec} and {\tt CRunDec}
F. Herren and M. Steinhauser,Version 3 of RunDec and CRunDec,Comput. Phys. Commun. 224(2018) 333 [1703.03751]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[34]
H. Frellesvig, D. Tommasini and C. Wever,On the reduction of generalized polylogarithms to Li𝑛 and Li2,2 and on the evaluation thereof,JHEP03(2016) 189 [1601.02649]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[35]
Campillo Aveleira et al.,𝑔𝑔→𝑍 𝐻: updated predictions at NLO QCD, 8, 2025, 2508.09905
B. Campillo Aveleira et al.,𝑔𝑔→𝑍 𝐻: updated predictions at NLO QCD, 8, 2025, 2508.09905. 8
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.