pith. sign in

arxiv: 2606.05300 · v1 · pith:P6NWHW53new · submitted 2026-06-03 · ✦ hep-ph · hep-ex

NNLO+PS Higgs-pair production in MiNNLOPS

Pith reviewed 2026-06-28 05:18 UTC · model grok-4.3

classification ✦ hep-ph hep-ex
keywords Higgs pair productiongluon fusionNNLO QCDparton shower matchingtop quark mass effectstrilinear Higgs couplingMiNNLOPS
0
0 comments X

The pith

NNLO QCD corrections for Higgs pair production are matched to parton showers using top-mass approximations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a matched NNLO plus parton shower prediction for Higgs boson pair production through gluon fusion at hadron colliders. This is done within the MiNNLOPS framework, using approximations for the top quark mass effects at NNLO based on the exact NLO result and the available two-loop amplitude. The approach includes exact contributions for the Born, single-virtual, single-real and double-real terms while approximating the others, and different approximation variants are tested for uncertainty assessment. Such predictions are needed to improve the accuracy of theoretical modeling for di-Higgs production, which is sensitive to the Higgs trilinear self-coupling. The results are validated against fixed-order NNLO calculations and compared to other matched predictions, with phenomenological studies for various decay modes and coupling variations provided.

Core claim

The central claim is that NNLO QCD corrections to Higgs-boson pair production in gluon fusion can be matched to parton showers in the MiNNLOPS framework by incorporating finite top-quark mass effects through approximations based on the exact NLO QCD result and the two-loop amplitude in the full theory, with Born, single-virtual, single-real and double-real contributions treated exactly and the real-virtual and double-virtual corrections approximated in different ways to assess uncertainties.

What carries the argument

The MiNNLOPS matching procedure combined with the scheme for approximating higher-order top-mass dependent corrections using the exact NLO result.

If this is right

  • Validation shows agreement with fixed-order NNLO QCD results.
  • Comparison with GENEVA reveals noticeable differences in some cases.
  • Phenomenological results are presented for different Higgs decay channels.
  • Variations of the trilinear Higgs coupling are explored in the predictions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • These matched predictions can help quantify the theoretical uncertainties in measurements of the Higgs self-coupling at the LHC.
  • Once the full NNLO top-mass dependence becomes available, the approximation uncertainties can be reduced or eliminated.
  • Similar matching techniques could be applied to other loop-induced processes with incomplete higher-order mass corrections.

Load-bearing premise

The approximations used for the real-virtual and double-virtual corrections provide a reliable estimate of the finite top-quark mass effects at NNLO order.

What would settle it

A calculation of the exact NNLO QCD corrections with full top-quark mass dependence, when available, could be compared to these approximate results to check for discrepancies beyond the estimated uncertainties.

read the original abstract

We consider Higgs-boson pair production in gluon fusion at hadron colliders and match next-to-next-to-leading-order (NNLO) QCD corrections to parton showers within the MiNNLO$_{PS}$ framework. Since the full top-quark mass dependence at this order is not available, finite top-quark mass effects are incorporated through approximations based on the exact NLO QCD result, using the available two-loop amplitude in the full theory. Specifically, the Born, single-virtual, single-real and double-real contributions are included exactly, while the real--virtual and double-virtual corrections are approximated. We consider different approximations for the latter to assess the associated uncertainties. We validate our predictions against fixed-order NNLO QCD results and compare with existing NNLO calculations matched to parton shower from GENEVA, where in some cases we find noticeable differences. Finally, we present phenomenological results for different Higgs-decay channels and variations of the trilinear Higgs coupling. Our MiNNLO$_{PS}$ generator for Higgs-boson pair production is available within the POWHEG-BOX-RES framework.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 0 minor

Summary. The manuscript presents an implementation of NNLO QCD corrections for Higgs-boson pair production in gluon fusion, matched to parton showers via the MiNNLOPS method inside the POWHEG-BOX-RES framework. Finite top-quark mass dependence is retained exactly for the Born, single-virtual, single-real and double-real contributions, while real-virtual and double-virtual terms are approximated using the available two-loop amplitude; several variants of the approximation are employed to estimate uncertainties. The predictions are validated against fixed-order NNLO results, compared with GENEVA-matched calculations, and applied to phenomenological studies of decay channels and trilinear-coupling variations.

Significance. If the central matching remains formally accurate under the stated approximations, the work supplies a publicly available generator that extends precision di-Higgs modeling to NNLO+PS level, which is directly relevant for LHC analyses. The explicit uncertainty assessment from multiple approximation schemes and the direct comparison with an independent matching framework (GENEVA) are constructive features.

major comments (1)
  1. [Implementation of finite top-mass effects (as described after the abstract statement of contributions included exactly v] The manuscript states that real-virtual and double-virtual corrections are approximated while Born, single-virtual, single-real and double-real contributions are kept exact. Because the MiNNLOPS construction relies on exact infrared cancellation between the NNLO hard function (including the double-virtual piece) and the real-emission terms, it is not demonstrated that the chosen approximation preserves the pole structure and finite parts to the level required for formal NNLO accuracy of the showered cross section. Fixed-order validation alone does not address this point for the matched prediction.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment of our work and for the detailed major comment. We address the concern regarding the finite top-mass approximation and its implications for formal NNLO accuracy in the MiNNLOPS matching below.

read point-by-point responses
  1. Referee: The manuscript states that real-virtual and double-virtual corrections are approximated while Born, single-virtual, single-real and double-real contributions are kept exact. Because the MiNNLOPS construction relies on exact infrared cancellation between the NNLO hard function (including the double-virtual piece) and the real-emission terms, it is not demonstrated that the chosen approximation preserves the pole structure and finite parts to the level required for formal NNLO accuracy of the showered cross section. Fixed-order validation alone does not address this point for the matched prediction.

    Authors: We appreciate the referee raising this important technical point. The real-virtual and double-virtual contributions are approximated using the available two-loop amplitude computed in the full theory (with exact top-mass dependence). This construction ensures that the infrared pole structure of the NNLO hard function is reproduced exactly, as the poles are fixed by universal factorization and the exact lower-order amplitudes already included in the calculation; only the finite remainders are approximated. Different approximation schemes are used to estimate the associated uncertainty. The fixed-order validation against the full NNLO result confirms that the cancellations hold numerically to high precision within the quoted uncertainties. While a formal all-order proof of NNLO accuracy for the matched prediction would require the complete two-loop amplitude (which is unavailable), the structure of the MiNNLOPS method combined with exact pole preservation supports the accuracy of the showered result. To address the referee's concern explicitly, we will revise the manuscript to include a dedicated discussion of how the pole structure is maintained under the approximation and its consequences for the matched cross section. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper applies the established MiNNLOPS matching framework (external to this work) to Higgs-pair production, using exact Born/single-virtual/single-real/double-real contributions and NLO-based approximations only for the real-virtual and double-virtual terms. No derivation step reduces by construction to a fitted parameter, self-defined observable, or load-bearing self-citation chain; the central predictions remain independently falsifiable against fixed-order NNLO benchmarks and external codes such as GENEVA. The approximations are explicitly described as ad-hoc for uncertainty estimation rather than being tuned to the showered observables themselves.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the MiNNLOPS matching procedure, standard perturbative QCD at NNLO, and the chosen approximations for virtual corrections; no new free parameters, axioms beyond domain standards, or invented entities are introduced.

axioms (1)
  • domain assumption Perturbative QCD expansion at NNLO with parton-shower matching is valid when higher-order mass-dependent terms are approximated from NLO results
    Full top-mass dependence at NNLO is unavailable, so the paper explicitly adopts this approximation strategy.

pith-pipeline@v0.9.1-grok · 5720 in / 1335 out tokens · 31181 ms · 2026-06-28T05:18:16.589700+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

53 extracted references · 25 linked inside Pith

  1. [1]

    Czakon, R

    M. Czakon, R. V. Harlander, J. Klappert and M. Niggetiedt,Exact Top-Quark Mass Dependence in Hadronic Higgs Production,Phys. Rev. Lett.127(2021) 162002, [2105.04436]. [Erratum: Phys.Rev.Lett. 131, 179901 (2023)]

  2. [2]

    Niggetiedt and M

    M. Niggetiedt and M. Wiesemann,Higgs-boson production in the full theory at NNLO+PS, Phys. Lett. B858(2024) 139043, [2407.01354]

  3. [3]

    E. W. N. Glover and J. J. van der Bij,Higgs Boson Pair Production via Gluon Fusion,Nucl. Phys. B309(1988) 282–294

  4. [4]

    Plehn, M

    T. Plehn, M. Spira and P. M. Zerwas,Pair production of neutral Higgs particles in gluon-gluon collisions,Nucl. Phys. B479(1996) 46–64, [hep-ph/9603205]. [Erratum: Nucl. Phys. B 531 (1998) 655–655]

  5. [5]

    Dawson, S

    S. Dawson, S. Dittmaier and M. Spira,Neutral Higgs boson pair production at hadron colliders: QCD corrections,Phys. Rev. D58(1998) 115012, [hep-ph/9805244]

  6. [6]

    Borowka, N

    S. Borowka, N. Greiner, G. Heinrich, S. P. Jones, M. Kerner, J. Schlenk et al.,Higgs Boson Pair Production in Gluon Fusion at Next-to-Leading Order with Full Top-Quark Mass Dependence,Phys. Rev. Lett.117(2016) 012001, [1604.06447]. [Erratum: Phys. Rev. Lett. 117 (2016) 079901]

  7. [7]

    Borowka, N

    S. Borowka, N. Greiner, G. Heinrich, S. P. Jones, M. Kerner, J. Schlenk et al.,Full top quark mass dependence in Higgs boson pair production at NLO,JHEP10(2016) 107, [1608.04798]

  8. [8]

    de Florian and J

    D. de Florian and J. Mazzitelli,Higgs Boson Pair Production at Next-to-Next-to-Leading Order in QCD,Phys. Rev. Lett.111(2013) 201801, [1309.6594]. – 24 –

  9. [9]

    Grigo, K

    J. Grigo, K. Melnikov and M. Steinhauser,Virtual corrections to Higgs boson pair production in the large top quark mass limit,Nucl. Phys. B888(2014) 17–29, [1408.2422]

  10. [10]

    de Florian, M

    D. de Florian, M. Grazzini, C. Hanga, S. Kallweit, J. M. Lindert, P. Maierh¨ ofer et al., Differential Higgs Boson Pair Production at Next-to-Next-to-Leading Order in QCD,JHEP09 (2016) 151, [1606.09519]

  11. [11]

    de Florian and J

    D. de Florian and J. Mazzitelli,Higgs pair production at next-to-next-to-leading logarithmic accuracy at the LHC,JHEP09(2015) 053, [1505.07122]

  12. [12]

    H.-S. Shao, H. T. Li, C. S. Li and J. Wang,Threshold resummation effects in Higgs boson pair production at the LHC,JHEP07(2013) 169, [1301.1245]

  13. [13]

    Davies, G

    J. Davies, G. Mishima, M. Steinhauser and D. Wellmann,Double-Higgs boson production in the high-energy limit: planar master integrals,JHEP03(2018) 048, [1801.09696]

  14. [14]

    L.-B. Chen, H. T. Li, H.-S. Shao and J. Wang,Higgs boson pair production via gluon fusion at N3LO in QCD,Phys. Lett. B803(2020) 135292, [1909.06808]

  15. [15]

    Grazzini, G

    M. Grazzini, G. Heinrich, S. Jones, S. Kallweit, M. Kerner, J. M. Lindert et al.,Higgs boson pair production at NNLO with top quark mass effects,JHEP05(2018) 059, [1803.02463]

  16. [16]

    Frederix, S

    R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, P. Torrielli et al.,Higgs pair production at the LHC with NLO and parton-shower effects,Phys. Lett. B732(2014) 142–149, [1401.7340]

  17. [17]

    Alioli, G

    S. Alioli, G. Billis, A. Broggio, A. Gavardi, S. Kallweit, M. A. Lim et al.,Double Higgs production at NNLO interfaced to parton showers in GENEVA,JHEP06(2023) 205, [2212.10489]

  18. [18]

    Alioli, G

    S. Alioli, G. Marinelli and D. Napoletano,NNLO+PS double Higgs boson production with top-quark mass corrections in GENEVA,JHEP09(2025) 206, [2507.08558]

  19. [19]

    P. F. Monni, P. Nason, E. Re, M. Wiesemann and G. Zanderighi,MiNNLOPS: a new method to match NNLO QCD to parton showers,JHEP05(2020) 143, [1908.06987]

  20. [20]

    Jeˇ zo and P

    T. Jeˇ zo and P. Nason,On the treatment of resonances in next-to-leading order calculations matched to a parton shower,JHEP12(2015) 065, [1509.09071]

  21. [21]

    P. F. Monni, E. Re and M. Wiesemann,MiNNLOPS: optimizing2→1hadronic processes, Eur. Phys. J. C80(2020) 1075, [2006.04133]

  22. [22]

    Lombardi, M

    D. Lombardi, M. Wiesemann and G. Zanderighi,Advancing MiNNLOPS to diboson processes: Zγproduction at NNLO+PS,JHEP06(2021) 095, [2010.10478]

  23. [23]

    Mazzitelli, P

    J. Mazzitelli, P. Nason, E. Re, M. Wiesemann and G. Zanderighi,Next-to-next-to-leading order event generation for top-quark pair production,Phys. Rev. Lett.127(2021) 062001, [2012.14267]

  24. [24]

    Mazzitelli, V

    J. Mazzitelli, V. Sotnikov and M. Wiesemann,Next-to-next-to-leading order event generation forZ-boson production in association with a bottom-quark pair,2404.08598

  25. [25]

    Lombardi, M

    D. Lombardi, M. Wiesemann and G. Zanderighi,W+W− production at NNLO+PS with MiNNLOPS,JHEP11(2021) 230, [2103.12077]

  26. [26]

    Mazzitelli, P

    J. Mazzitelli, P. F. Monni, P. Nason, E. Re, M. Wiesemann and G. Zanderighi,Top-pair production at the LHC with MINNLOPS,JHEP04(2022) 079, [2112.12135]

  27. [27]

    Buonocore, G

    L. Buonocore, G. Koole, D. Lombardi, L. Rottoli, M. Wiesemann and G. Zanderighi,ZZ production at nNNLO+PS with MiNNLOPS,JHEP01(2022) 072, [2108.05337]. – 25 –

  28. [28]

    Lombardi, M

    D. Lombardi, M. Wiesemann and G. Zanderighi,Anomalous couplings in Zγevents at NNLO+PS and improvingνν¯γbackgrounds in dark-matter searches,Phys. Lett. B824 (2022) 136846, [2108.11315]

  29. [29]

    Zanoli, M

    S. Zanoli, M. Chiesa, E. Re, M. Wiesemann and G. Zanderighi,Next-to-next-to-leading order event generation for VH production with H→bbdecay,JHEP07(2022) 008, [2112.04168]

  30. [30]

    Gavardi, C

    A. Gavardi, C. Oleari and E. Re,NNLO+PS Monte Carlo simulation of photon pair production with MiNNLOP S,JHEP09(2022) 061, [2204.12602]

  31. [31]

    Haisch, D

    U. Haisch, D. J. Scott, M. Wiesemann, G. Zanderighi and S. Zanoli,NNLO event generation forpp→Zh→ℓ +ℓ−bbproduction in the SM effective field theory,JHEP07(2022) 054, [2204.00663]

  32. [32]

    J. M. Lindert, D. Lombardi, M. Wiesemann, G. Zanderighi and S. Zanoli,WZ production at NNLO QCD and NLO EW matched to parton showers with MiNNLOPS,JHEP11(2022) 036, [2208.12660]

  33. [33]

    Mazzitelli, A

    J. Mazzitelli, A. Ratti, M. Wiesemann and G. Zanderighi,B-hadron production at the LHC from bottom-quark pair production at NNLO+PS,Phys. Lett. B843(2023) 137991, [2302.01645]

  34. [34]

    Biello, A

    C. Biello, A. Sankar, M. Wiesemann and G. Zanderighi,NNLO+PS predictions for Higgs production through bottom-quark annihilation with MINNLOPS,Eur. Phys. J. C84(2024) 479, [2402.04025]

  35. [35]

    Biello, J

    C. Biello, J. Mazzitelli, A. Sankar, M. Wiesemann and G. Zanderighi,Higgs boson production in association with massive bottom quarks at NNLO+PS,JHEP04(2025) 088, [2412.09510]

  36. [36]

    Biello, C

    C. Biello, C. Savoini, C. Signorile-Signorile and M. Wiesemann,Next-to-next-to-leading order event generation fort¯tHproduction with approximate two-loop amplitude,2603.06143

  37. [37]

    Nason,A New method for combining NLO QCD with shower Monte Carlo algorithms, JHEP11(2004) 040, [hep-ph/0409146]

    P. Nason,A New method for combining NLO QCD with shower Monte Carlo algorithms, JHEP11(2004) 040, [hep-ph/0409146]

  38. [38]

    Frixione, P

    S. Frixione, P. Nason and C. Oleari,Matching NLO QCD computations with Parton Shower simulations: the POWHEG method,JHEP11(2007) 070, [0709.2092]

  39. [39]

    Alioli, P

    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]

  40. [40]

    Cascioli, P

    F. Cascioli, P. Maierh¨ ofer and S. Pozzorini,Scattering Amplitudes with Open Loops,Phys. Rev. Lett.108(2012) 111601, [1111.5206]

  41. [41]

    Buccioni, S

    F. Buccioni, S. Pozzorini and M. F. Zoller,On-the-fly reduction of open loops,Eur. Phys. J. C 78(2018) 70, [1710.11452]

  42. [42]

    Buccioni, J.-N

    F. Buccioni, J.-N. Lang, J. M. Lindert, P. Maierh¨ ofer, S. Pozzorini, H. Zhang et al.,OpenLoops 2,Eur. Phys. J. C79(2019) 866, [1907.13071]

  43. [43]

    Becher and M

    T. Becher and M. Neubert,On the Structure of Infrared Singularities of Gauge-Theory Amplitudes,JHEP06(2009) 081, [0903.1126]. [Erratum: JHEP 11, 024 (2013)]

  44. [44]

    Becher and M

    T. Becher and M. Neubert,Drell-Yan Production at SmallqT, Transverse Parton Distributions and the Collinear Anomaly,Eur. Phys. J. C71(2011) 1665, [1007.4005]

  45. [45]

    Heinrich, S

    G. Heinrich, S. P. Jones, M. Kerner, G. Luisoni and E. Vryonidou,NLO predictions for Higgs – 26 – boson pair production with full top quark mass dependence matched to parton showers,JHEP 08(2017) 088, [1703.09252]

  46. [46]

    Davies, G

    J. Davies, G. Heinrich, S. P. Jones, M. Kerner, G. Mishima, M. Steinhauser et al.,Double Higgs boson production at NLO: combining the exact numerical result and high-energy expansion,JHEP11(2019) 024, [1907.06408]

  47. [47]

    Davies, K

    J. Davies, K. Sch¨ onwald, M. Steinhauser and D. Stremmer,ggxy: A flexible library to compute gluon-induced cross sections,Comput. Phys. Commun.320(2026) 109933, [2506.04323]

  48. [48]

    Sj¨ ostrand, S

    T. Sj¨ ostrand, S. Ask, J. R. Christiansen, R. Corke, N. Desai, P. Ilten et al.,An Introduction to PYTHIA 8.2,Comput. Phys. Commun.191(2015) 159–177, [1410.3012]

  49. [49]

    Grazzini, S

    M. Grazzini, S. Kallweit and M. Wiesemann,Fully differential NNLO computations with MATRIX,Eur. Phys. J. C78(2018) 537, [1711.06631]

  50. [50]

    Aad et al.,Study of Higgs boson pair production in theHH→b bγγ final state with 308 fb−1 of data collected at√s= 13TeV and 13.6 TeV by the ATLAS experiment,Phys

    ATLAScollaboration, G. Aad et al.,Study of Higgs boson pair production in theHH→b bγγ final state with 308 fb−1 of data collected at√s= 13TeV and 13.6 TeV by the ATLAS experiment,Phys. Lett. B876(2026) 140280, [2507.03495]

  51. [51]

    Cacciari, G

    M. Cacciari, G. P. Salam and G. Soyez,The anti-kt jet clustering algorithm,JHEP04(2008) 063, [0802.1189]

  52. [52]

    Aad et al.,Search for nonresonant pair production of Higgs bosons in the bbbb final state in pp collisions at√s = 13TeV with the ATLAS detector,Phys

    ATLAScollaboration, G. Aad et al.,Search for nonresonant pair production of Higgs bosons in the bbbb final state in pp collisions at√s = 13TeV with the ATLAS detector,Phys. Rev. D 108(2023) 052003, [2301.03212]

  53. [53]

    Aad et al.,Highlights of the HL-LHC physics projections by ATLAS and CMS,2504.00672

    ATLAS, CMScollaboration, G. Aad et al.,Highlights of the HL-LHC physics projections by ATLAS and CMS,2504.00672. – 27 –