pith. sign in

arxiv: 2605.31221 · v1 · pith:D6DK65HHnew · submitted 2026-05-29 · ✦ hep-ph · hep-th

CoLoRFulNNLO for color-singlet processes: An update on NNLOCAL

S. Van Thurenhout , V. Del Duca , C. Duhr , L. Fek\'esh\'azy , F. Guadagni , P. Mukherjee , G. Somogyi , F. Tramontano This is my paper

Pith reviewed 2026-06-28 22:00 UTC · model grok-4.3

classification ✦ hep-ph hep-th
keywords NNLO QCDcolor-singlet productionlocal subtraction schemeinfrared factorizationMonte Carlo integrationLHC phenomenologycounterterms
0
0 comments X

The pith

A generic local subtraction scheme based on standard QCD factorization enables analytic NNLO calculations for color-singlet production.

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

The paper updates the NNLOCAL Monte Carlo program that extends the CoLoRFulNNLO subtraction method to color-singlet final states in hadron collisions. It shows that counterterms can be built directly from standard infrared factorization formulas without needing process-by-process adjustments. Analytic integration of those counterterms over unresolved emissions is performed to maintain numerical stability in predictions. The approach is validated by producing fully differential NNLO cross sections for LHC processes. A sympathetic reader would care because this offers a uniform way to reach next-to-next-to-leading order accuracy for key observables such as Higgs or vector-boson production.

Core claim

The construction of the counterterms in our scheme is generic, being based on the standard IR factorization formulae of QCD. Furthermore, the integration of the counterterms over the phase space of unresolved emissions is performed fully analytically, allowing for good control of the numerical stability of our predictions. We validate our method by computing NNLO corrections to fully differential cross sections for the LHC.

What carries the argument

The CoLoRFulNNLO completely local subtraction scheme, whose counterterms are derived from standard QCD infrared factorization formulae and integrated analytically over unresolved parton phase space.

If this is right

  • NNLO corrections to fully differential distributions for any color-singlet final state become accessible with the same subtraction infrastructure.
  • Analytic counterterm integration removes the need for process-dependent numerical tuning to achieve stability.
  • The method applies uniformly across hadron-collider observables without requiring separate counterterm derivations for each process.
  • Validation against LHC data confirms that the local cancellation of infrared poles works at the level of differential cross sections.

Where Pith is reading between the lines

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

  • The same analytic-integration approach could be tested on processes with colored particles by adding the corresponding factorization terms.
  • Full analytic control over counterterms may simplify matching the fixed-order results to parton-shower resummation.
  • Improved numerical stability might reduce uncertainties in global PDF fits that rely on NNLO predictions.

Load-bearing premise

Standard infrared factorization formulae of QCD are enough to construct every needed counterterm for a fully local subtraction scheme in color-singlet production without extra process-specific terms.

What would settle it

A numerical discrepancy between the program's NNLO results and an established benchmark calculation for a color-singlet process such as Drell-Yan or Higgs production at the LHC would show the scheme does not hold.

Figures

Figures reproduced from arXiv: 2605.31221 by C. Duhr, F. Guadagni, F. Tramontano, G. Somogyi, L. Fek\'esh\'azy, P. Mukherjee, S. Van Thurenhout, V. Del Duca.

Figure 1
Figure 1. Figure 1: Cancellation of kinematic singularities in various IR limits. From left to right, we show [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Rapidity distributions of the Higgs boson with bin widths ∆ [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
read the original abstract

We give an update on the status of NNLOCAL, a parton-level Monte Carlo program implementing the extension of the completely local subtraction scheme CoLoRFulNNLO to color-singlet production in hadron-hadron collisions. The construction of the counterterms in our scheme is generic, being based on the standard IR factorization formulae of QCD. Furthermore, the integration of the counterterms over the phase space of unresolved emissions is performed fully analytically, allowing for good control of the numerical stability of our predictions. We validate our method by computing NNLO corrections to fully differential cross sections for the LHC.

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 provides an update on NNLOCAL, a parton-level Monte Carlo program extending the CoLoRFulNNLO completely local subtraction scheme to color-singlet production in hadron-hadron collisions. It states that counterterms are constructed generically from standard QCD IR factorization formulae, with their integration over unresolved emissions performed fully analytically to ensure numerical stability. Validation is asserted through computation of NNLO corrections to fully differential cross sections for LHC processes.

Significance. If the claims hold, the work advances local subtraction methods for NNLO QCD calculations of color-singlet processes by providing a generic construction and analytic integration that supports numerical stability. The reliance on externally established QCD IR factorization formulae is a strength, as is the emphasis on analytic control over counterterm integration.

major comments (1)
  1. [Abstract] Abstract: The manuscript asserts validation of the method by computing NNLO corrections to fully differential cross sections for the LHC, but supplies no numerical results, error estimates, or comparison data. This prevents assessment of whether the claimed analytic integration and numerical stability are achieved in practice.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their review and the constructive comment on the abstract. We respond to the major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The manuscript asserts validation of the method by computing NNLO corrections to fully differential cross sections for the LHC, but supplies no numerical results, error estimates, or comparison data. This prevents assessment of whether the claimed analytic integration and numerical stability are achieved in practice.

    Authors: We agree that the current abstract claims validation by computation of NNLO corrections without supplying any numerical results, error estimates or comparisons within this manuscript. This manuscript is an update on the implementation status of NNLOCAL, emphasizing the generic construction from QCD IR factorization formulae and the fully analytic integration of counterterms. The explicit numerical validations for LHC processes were presented in our prior publications on the method. To address the referee's concern, we will revise the abstract to remove the claim of validation within this work and instead state that the analytic construction enables stable numerical predictions, with examples available in earlier references. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation rests on external standard QCD IR factorization

full rationale

The paper states that counterterms are constructed generically from the standard IR factorization formulae of QCD and integrated analytically, with validation on LHC processes. This relies on externally established QCD results rather than any self-definition, fitted inputs renamed as predictions, or load-bearing self-citations. No equation or claim reduces by construction to the paper's own inputs; the central method is presented as an application of independent prior knowledge, rendering the derivation self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the direct applicability of standard QCD IR factorization formulae to construct local counterterms for color-singlet processes; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption Standard IR factorization formulae of QCD apply directly to the construction of counterterms in the CoLoRFulNNLO scheme for color-singlet production.
    Explicitly stated in the abstract as the basis for the generic construction of counterterms.

pith-pipeline@v0.9.1-grok · 5663 in / 1233 out tokens · 39855 ms · 2026-06-28T22:00:43.877510+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

36 extracted references · 21 linked inside Pith

  1. [1]

    Nakanishi,General Theory of Infrared Divergence,Progress of Theoretical Physics19(02, 1958) 159

    N. Nakanishi,General Theory of Infrared Divergence,Progress of Theoretical Physics19(02, 1958) 159

    1958

  2. [2]

    Kinoshita,Mass singularities of Feynman amplitudes,J

    T. Kinoshita,Mass singularities of Feynman amplitudes,J. Math. Phys.3(1962) 650

    1962

  3. [3]

    Lee and M

    T.D. Lee and M. Nauenberg,Degenerate Systems and Mass Singularities,Phys. Rev.133(1964) B1549

    1964

  4. [4]

    Frixione, Z

    S. Frixione, Z. Kunszt and A. Signer,Three jet cross-sections to next-to-leading order,Nucl. Phys. B467(1996) 399 [hep-ph/9512328]

    Pith/arXiv arXiv 1996

  5. [5]

    Catani and M.H

    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]

    Pith/arXiv arXiv 1997

  6. [6]

    Heimel, R

    T. Heimel, R. Winterhalder, A. Butter, J. Isaacson, C. Krause, F. Maltoni et al.,MadNIS - Neural multi-channel importance sampling,SciPost Phys.15(2023) 141 [arXiv:2212.06172]

    arXiv 2023

  7. [7]

    De Crescenzo, J.M

    G. De Crescenzo, J.M. Villadamigo, N. Elmer, T. Heimel, T. Plehn, R. Winterhalder et al.,MadNIS at NLO, arXiv:2603.22407

    Pith/arXiv arXiv

  8. [8]

    Fabricius, I

    K. Fabricius, I. Schmitt, G. Kramer and G. Schierholz,Higher Order Perturbative QCD Calculation of Jet Cross-Sections in e+ e- Annihilation,Z. Phys. C11(1981) 315

    1981

  9. [9]

    Ellis, D.A

    R.K. Ellis, D.A. Ross and A.E. Terrano,The Perturbative Calculation of Jet Structure in e+ e- Annihilation,Nucl. Phys. B178(1981) 421

    1981

  10. [10]

    Catani and M

    S. Catani and M. Grazzini,An NNLO subtraction formalism in hadron collisions and its application to Higgs boson production at the LHC,Phys. Rev. Lett.98(2007) 222002 [hep-ph/0703012]. 7

    Pith/arXiv arXiv 2007

  11. [11]

    Gaunt, M

    J. Gaunt, M. Stahlhofen, F.J. Tackmann and J.R. Walsh,N-jettiness Subtractions for NNLO QCD Calculations,JHEP 09(2015) 058 [arXiv:1505.04794]

    Pith/arXiv arXiv 2015

  12. [12]

    Boughezal, J.M

    R. Boughezal, J.M. Campbell, R.K. Ellis, C. Focke, W. Giele, X. Liu et al.,Color singlet production at NNLO in MCFM,Eur. Phys. J. C77(2017) 7 [arXiv:1605.08011]

    Pith/arXiv arXiv 2017

  13. [13]

    Grazzini, S

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

    Pith/arXiv arXiv 2018

  14. [14]

    Gehrmann-De Ridder, T

    A. Gehrmann-De Ridder, T. Gehrmann and E.W.N. Glover,Antenna subtraction at NNLO,JHEP09(2005) 056 [hep-ph/0505111]

    Pith/arXiv arXiv 2005

  15. [15]

    Caola, K

    F. Caola, K. Melnikov and R. R¨ ontsch,Nested soft-collinear subtractions in NNLO QCD computations,Eur. Phys. J. C 77(2017) 248 [arXiv:1702.01352]

    Pith/arXiv arXiv 2017

  16. [16]

    Czakon and D

    M. Czakon and D. Heymes,Four-dimensional formulation of the sector-improved residue subtraction scheme,Nucl. Phys. B890(2014) 152 [arXiv:1408.2500]

    Pith/arXiv arXiv 2014

  17. [17]

    Magnea, E

    L. Magnea, E. Maina, G. Pelliccioli, C. Signorile-Signorile, P. Torrielli and S. Uccirati,Local analytic sector subtraction at NNLO,JHEP12(2018) 107 [arXiv:1806.09570]

    Pith/arXiv arXiv 2018

  18. [18]

    Cacciari, F.A

    M. Cacciari, F.A. Dreyer, A. Karlberg, G.P. Salam and G. Zanderighi,Fully Differential Vector-Boson-Fusion Higgs Production at Next-to-Next-to-Leading Order,Phys. Rev. Lett.115(2015) 082002 [arXiv:1506.02660]

    Pith/arXiv arXiv 2015

  19. [19]

    Herzog,Geometric IR subtraction for final state real radiation,JHEP08(2018) 006 [arXiv:1804.07949]

    F. Herzog,Geometric IR subtraction for final state real radiation,JHEP08(2018) 006 [arXiv:1804.07949]

    Pith/arXiv arXiv 2018

  20. [20]

    Somogyi, Z

    G. Somogyi, Z. Trocsanyi and V. Del Duca,Matching of singly- and doubly-unresolved limits of tree-level QCD squared matrix elements,JHEP06(2005) 024 [hep-ph/0502226]. [21]NNLOJETcollaboration, A. Huss et al.,NNLOJET: a parton-level event generator for jet cross sections at NNLO QCD accuracy,arXiv:2503.22804

    Pith/arXiv arXiv 2005

  21. [21]

    Klein and L

    S.Y. Klein and L. Simon,history: A tool for fully-differential cross sections at next-to-next-to-leading order, arXiv:2603.13729

    Pith/arXiv arXiv

  22. [22]

    Del Duca, C

    V. Del Duca, C. Duhr, L. Fekeshazy, F. Guadagni, P. Mukherjee, G. Somogyi et al.,NNLOCAL: completely local subtractions for color-singlet production in hadron collisions,JHEP05(2024) 151 [arXiv:2412.21028]

    arXiv 2024

  23. [23]

    Del Duca, G

    V. Del Duca, G. Somogyi and F. Tramontano,CoLoRFulNNLO for hadron collisions: regularizing initial-state double real emissions,arXiv:2512.05192

    arXiv

  24. [24]

    Chetyrkin and F.V

    K.G. Chetyrkin and F.V. Tkachov,Integration by parts: The algorithm to calculateβ-functions in 4 loops,Nucl. Phys. B192(1981) 159

    1981

  25. [25]

    Henn,Multiloop integrals in dimensional regularization made simple,Phys

    J.M. Henn,Multiloop integrals in dimensional regularization made simple,Phys. Rev. Lett.110(2013) 251601 [arXiv:1304.1806]

    Pith/arXiv arXiv 2013

  26. [26]

    Guadagni, V

    F. Guadagni, V. Del Duca, C. Duhr, L. Fek´ esh´ azy, P. Mukherjee, G. Somogyi et al.,NNLOCAL: completely local subtractions,PoSEPS-HEP2025(2026) 243 [arXiv:2601.02033]

    arXiv 2026

  27. [27]

    Mukherjee, V

    P. Mukherjee, V. Del Duca, C. Duhr, L. Fek´ esh´ azy, F. Guadagni, G. Somogyi et al.,NNLOCAL: Fully Local Subtractions for Precision Predictions in Hadron Collisions, in17th International Symposium on Radiative Corrections: Applications of Quantum Field Theory to Phenomenology, 2, 2026.arXiv:2602.24257

    arXiv 2026

  28. [28]

    Goncharov,Multiple polylogarithms, cyclotomy and modular complexes,Math

    A.B. Goncharov,Multiple polylogarithms, cyclotomy and modular complexes,Math. Res. Lett.5(1998) 497 [arXiv:1105.2076]

    Pith/arXiv arXiv 1998

  29. [29]

    Van Thurenhout, V

    S. Van Thurenhout, V. Del Duca, C. Duhr, L. Fek´ esh´ azy, F. Guadagni, P. Mukherjee et al.,CoLoRFul for Hadron Collisions: Integrating the Counterterms,Acta Phys. Polon. Supp.18(2025) 6 [arXiv:2412.12750]

    arXiv 2025

  30. [30]

    Fek´ esh´ azy, G

    L. Fek´ esh´ azy, G. Somogyi and S. Van Thurenhout,CoLoRFulNNLO for hadron collisions: integrating the iterated single-unresolved subtraction terms,JHEP05(2026) 035 [arXiv:2512.24403]

    Pith/arXiv arXiv 2026

  31. [31]

    Chargeishvili, L

    B. Chargeishvili, L. Fek´ esh´ azy, G. Somogyi and S. Van Thurenhout,LinApart: Optimizing the univariate partial fraction decomposition,Comput. Phys. Commun.307(2025) 109395 [arXiv:2405.20130]

    arXiv 2025

  32. [32]

    Fek´ esh´ azy and O

    L. Fek´ esh´ azy and O. Schnetz,LinApart2: efficient parallel partial fraction decomposition algorithm for denominators with polynomials of general degree,arXiv:2511.15735

    arXiv

  33. [33]

    Campbell and R.K

    J.M. Campbell and R.K. Ellis,MCFM for the Tevatron and the LHC,Nucl. Phys. B Proc. Suppl.205-206(2010) 10 [arXiv:1007.3492]. 8

    Pith/arXiv arXiv 2010

  34. [34]

    Buckley, J

    A. Buckley, J. Ferrando, S. Lloyd, K. Nordstr¨ om, B. Page, M. R¨ ufenacht et al.,LHAPDF6: parton density access in the LHC precision era,Eur. Phys. J. C75(2015) 132 [arXiv:1412.7420]

    Pith/arXiv arXiv 2015

  35. [35]

    Baglio, C

    J. Baglio, C. Duhr, B. Mistlberger and R. Szafron,Inclusive production cross sections at N 3LO,JHEP12(2022) 066 [arXiv:2209.06138]

    arXiv 2022

  36. [36]

    Beneke and V.A

    M. Beneke and V.A. Smirnov,Asymptotic expansion of Feynman integrals near threshold,Nucl. Phys. B522(1998) 321 [hep-ph/9711391]. 9

    Pith/arXiv arXiv 1998