REVIEW 3 major objections 3 minor 44 references
First exact two-loop amplitude for top-antitop-W production at leading colour.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-04 05:20 UTC pith:OM2ELYHA
load-bearing objection Plausible, honest status report; the real check is in the companion NNLO paper, not here. the 3 major comments →
Two-loop amplitude for tbar{t}W production at hadron colliders in the leading colour approximation
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
On the paper's own terms, the discovery is that the leading-colour two-loop amplitude for top-antitop-plus-W boson production is computable exactly, despite the combination of seven kinematic variables, massive internal propagators, three elliptic sectors, and one nested-square-root sector. The author constructs 330 master integrals, derives their differential equations, and shows that a specific choice of elliptic integrals — those with no terms below order ε^4 — lets the first three orders satisfy canonical polylogarithmic equations while elliptic effects enter only at the finite order. The finite remainder is expressed in 298 special functions with rational coefficients, evaluated numeric
What carries the argument
The central object is the system of first-order differential equations satisfied by the 330 master integrals. For the polylogarithmic sectors the connection matrix is built from logarithmic one-forms of algebraic functions of the kinematics; for the five-point elliptic sector, whose elliptic-curve discriminant is a degree-14 irreducible polynomial, a canonical basis is not constructed and the author instead uses a non-canonical form with minimally chosen non-logarithmic one-forms. The evaluation then splits into two parts: the special functions are obtained by solving their sparser, ε-independent differential equations with generalised series expansions, and the rational coefficients are rec
Load-bearing premise
The whole result stands on the correctness and completeness of the unshown differential equations for the 330 master integrals, especially the non-canonical system for the five-point elliptic sector with its degree-14 discriminant, and on the generalised series solutions converging to the true values.
What would settle it
Evaluate one master integral in the degree-14 elliptic sector independently at a fixed rational phase-space point — for example by direct numerical integration of its Feynman-parameter representation — and compare with the value obtained from the paper's differential-equation setup; a mismatch would show the answer is not the true amplitude.
If this is right
- Exact leading-colour NNLO predictions for ttbarW can be made point-by-point, replacing the soft-W/massification approximation in the bulk phase space.
- The finite remainder is available directly for differential distributions, since the pole subtraction is analytic rather than numerical.
- The treatment of the degree-14 elliptic sector shows that exact two-loop five-point amplitudes with massive propagators are numerically feasible at about one hour per point.
- The interpolation grid produced from this evaluation has already fed NNLO inclusive cross-section results, as reported by the author.
Where Pith is reading between the lines
- If this method scales, the same recipe—non-canonical differential equations for hard elliptic sectors plus finite-field reconstruction for coefficients—could be applied to other five-point two-loop processes with massive internal lines, such as ttbarH or four-top production, without first constructing canonical elliptic bases.
- The 'possibly over-complete' special-function basis suggests a smaller irreducible set may exist; identifying it could cut the per-point evaluation time well below one hour.
- Comparing exact results with the soft-W approximation point-by-point would show where the approximation loses accuracy, information the total-cross-section comparison cannot reveal.
- A natural follow-up, not in the paper, would be to release the actual differential equations and special-function basis so other groups can reproduce the numbers independently.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript, a Loops and Legs proceedings contribution, reports the first exact numerical evaluation of the two-loop QCD amplitude for pp -> t tbar W in the generalised leading colour approximation. The author describes the reduction of 8959 Feynman integrals to 330 master integrals, the construction of differential equations in seven kinematic variables, the emergence of three elliptic sectors and one nested-square-root sector, and the evaluation of the finite remainder by solving differential equations for special functions via generalised series expansions and by reconstructing rational coefficients over finite fields, with a reported cost of about one hour per phase-space point. The final result is claimed to be expressed in a possibly over-complete basis of special functions with rational coefficients, and the finite remainder is stated to involve 298 special functions.
Significance. The process is phenomenologically important: an exact two-loop amplitude would remove the soft-W approximation currently used in NNLO predictions and would allow validation of that approximation and of differential distributions. The paper has clear strengths: it is explicit about the conjectural input, it gives useful counts (330 MIs, 298/323 special functions, rational-coefficient prime counts), it uses established tools and companion papers, and it describes analytic pole subtraction and finite-field reconstruction. However, the proceedings contain no numerical value of any form factor or finite remainder, no master integrals, no differential equations, no special-function basis, and no comparison with known limits or approximations. As it stands, the central claim is a methodological status report rather than a checked and reproducible computation. If the underlying computation is correct, this is a significant milestone for the field.
major comments (3)
- [Section 4] The central claim of the paper, namely the numerical evaluation of the two-loop amplitude, is not supported by any numerical output. No value of a form factor or of the finite remainder is given at any phase-space point, and no comparison with the soft-W approximation of [6], with known sub-sectors, or with an independent numerical method is shown. Since the title and abstract assert that the amplitude has been evaluated, please provide at least one explicit phase-space point with rationalised invariants, the corresponding finite remainder or form factors, and a quantitative check (e.g., the soft-W limit or a comparison with the approximate result of [6]). Without such data the reader cannot assess whether the claimed evaluation has actually been performed.
- [Section 3.1, Eq. (13)] The treatment of the five-point elliptic sector is a load-bearing part of the computation, but it is described only schematically. No canonical basis is constructed for that sector, and the non-logarithmic one-forms omega_beta in Eq. (14) are not defined or shown. The paper states that the degree-14 Landau discriminant appears in the denominators and in the discriminant, but does not describe how the generalised series expansions of [39,40] are initialised and continued in this non-canonical setup. Please specify the form of the one-forms, the boundary conditions for the differential equations, and present a validation of the series expansions for the sector in Fig. 2c against a direct numerical integration or another reliable method.
- [Section 4] The analytic pole subtraction relies on the algebraic independence of the polylogarithmic special functions and on the conjecture that the only transcendental constants are zeta-values. The paper also states that the basis of special functions is 'possibly over-complete'. No completeness proof or independent evidence is provided for the finite-remainder basis. Given that a missing function or an incorrect linear relation among the f_k^{(4,*)} would change the finite remainder, please state the evidence that the set of 298 special functions is complete up to the relevant weight and that the elliptic functions f_k^{(4,*)} cannot mix with the polylogarithmic functions in a way that affects the finite part after IR subtraction, or give a precise reference to the companion papers where this is established.
minor comments (3)
- [Section 2, Eq. (5)] The colour expansion is written as N_c^2 A^(N_c^2) + N_c N_f A^(N_c N_f) + N_f^2 A^(N_f^2) + O(N_c). For the numerical values N_c=3, N_f=5 the N_f^2 term is not numerically suppressed relative to N_c^2; the meaning of 'leading colour' in this generalised sense could be clarified for readers not familiar with the 1/N_c expansion at fixed N_f/N_c.
- [Section 3] The reduction setup is described only by the names NeatIBP and FiniteFlow. To make the counts 8959 and 330 reproducible, please state the integral-family definition (or cite the companion paper) and the IBP generation/solving settings, e.g., the seed sectors and the chosen master-integral selection.
- [Throughout] There are minor typographical and formatting issues, such as inconsistent references to figures ('fig.1a' vs 'Fig. 1a') and the use of 'generalised leading colour' without a precise definition. These do not affect the technical content but should be polished.
Circularity Check
No circularity: standard IBPs/DEs/finite-field workflow; self-citations cite independent companion computations, not the conclusion being derived.
full rationale
The derivation chain is a standard reduction: Feynman diagrams are decomposed into 24 tensor structures; the resulting 8959 integrals are reduced by integration-by-parts identities to 330 master integrals; differential equations are derived for these master integrals; and the amplitude is expressed in a basis of special functions whose evaluation is separate from the rational coefficients. Nothing in this chain uses the final finite remainder as an input, and no parameter is fitted to a target prediction. The cited companion papers [12,13,46] are separate computations (one-loop epsilon^2 terms, the two-loop integral families, and a downstream NNLO application); they supply inputs or outputs rather than the argument's conclusion, so overlap in authorship does not make the derivation circular. The paper's own caveat - that it refrains from constructing a canonical basis for the five-point elliptic sector and instead uses the non-canonical form of eq. (13) - is a completeness and rigor limitation, not a circular reduction. Likewise, the absence of released differential equations or numerical benchmark values is a reproducibility concern, not evidence of circularity. No circular step could be located, so score 0 is appropriate.
Axiom & Free-Parameter Ledger
axioms (5)
- domain assumption The 't Hooft-Veltman (tHV) scheme and dimensional regularisation in D=4-2ε provide a consistent regulator for the amplitude.
- domain assumption The 330 master integrals obtained via NeatIBP and FiniteFlow form a complete basis and the IBPs are correct.
- domain assumption The differential equations for the MIs and for the special functions correctly encode the analytic dependence of the integrals.
- ad hoc to paper The conjecture that integrands with at most simple poles and constant leading singularity lead to canonical Feynman integrals.
- ad hoc to paper The algebraic-independence conjecture that the only transcendental constants appearing are zeta-values ζ_n.
read the original abstract
In this contribution I present the first exact calculation of the leading-colour two-loop QCD amplitude for the associated production of a top-anti-top pair and a W boson. I discuss strategies to address the complexity of the computation, which involves complicated analytic structures, such as nested square roots, elliptic functions, and expressions with a high degree of algebraic complexity. The final result is expressed in terms of a set of special functions, which are evaluated using the method of differential equations, and rational coefficients, evaluated via finite field techniques.
Figures
Reference graph
Works this paper leans on
-
[1]
A. Buckley, C. Englert, J. Ferrando, D. J. Miller, L. Moore, M. Russell, and C. D. White, Constraining top quark effective theory in the LHC Run II era,JHEP04(2016) 015, [arXiv:1512.03360]
Pith/arXiv arXiv 2016
-
[2]
J. A. Dror, M. Farina, E. Salvioni, and J. Serra,Strong tW Scattering at the LHC,JHEP01 (2016) 071, [arXiv:1511.03674]
Pith/arXiv arXiv 2016
-
[3]
O. Bessidskaia Bylund, F. Maltoni, I. Tsinikos, E. Vryonidou, and C. Zhang,Probing top quark neutral couplings in the Standard Model Effective Field Theory at NLO in QCD, JHEP05(2016) 052, [arXiv:1601.08193]. [4]CMSCollaboration, A. Tumasyan et al.,Measurement of the cross section of top quark-antiquark pair production in association with a W boson in pro...
Pith/arXiv arXiv 2016
-
[6]
L. Buonocore, S. Devoto, M. Grazzini, S. Kallweit, J. Mazzitelli, L. Rottoli, and C. Savoini, Precise Predictions for the Associated Production of a W Boson with a Top-Antitop Quark Pair at the LHC,Phys. Rev. Lett.131(2023), no. 23 231901, [arXiv:2306.16311]
Pith/arXiv arXiv 2023
-
[7]
S. Catani, S. Devoto, M. Grazzini, S. Kallweit, J. Mazzitelli, and C. Savoini,Higgs Boson Production in Association with a Top-Antitop Quark Pair in Next-to-Next-to-Leading Order QCD,Phys. Rev. Lett.130(2023), no. 11 111902, [arXiv:2210.07846]
Pith/arXiv arXiv 2023
-
[8]
P. Bärnreuther, M. Czakon, and P. Fiedler,Virtual amplitudes and threshold behaviour of hadronic top-quark pair-production cross sections,JHEP02(2014) 078, [arXiv:1312.6279]. 8 Two-loop amplitude for𝑡¯𝑡𝑊productionMattia Pozzoli
Pith/arXiv arXiv 2014
-
[9]
A. A. Penin,Two-loop photonic corrections to massive Bhabha scattering,Nucl. Phys. B734 (2006) 185–202, [hep-ph/0508127]
Pith/arXiv arXiv 2006
-
[10]
A. Mitov and S. Moch,The Singular behavior of massive QCD amplitudes,JHEP05(2007) 001, [hep-ph/0612149]
Pith/arXiv arXiv 2007
-
[11]
T. Becher and K. Melnikov,Two-loop QED corrections to Bhabha scattering,JHEP06 (2007) 084, [arXiv:0704.3582]
Pith/arXiv arXiv 2007
-
[12]
M. Becchetti, M. Delto, S. Ditsch, P. A. Kreer, M. Pozzoli, and L. Tancredi,One-loop QCD corrections to¯ud→𝑡 𝑡𝑊atO 𝜀2 ,JHEP09(2025) 126, [arXiv:2502.14952]
arXiv 2025
-
[13]
M. Becchetti, D. Canko, V. Chestnov, T. Peraro, M. Pozzoli, and S. Zoia,Two-loop Feynman integrals for leading colour𝑡𝑡𝑊production at hadron colliders,JHEP07(2025) 001, [arXiv:2504.13011]
Pith/arXiv arXiv 2025
-
[14]
Barucchi and G
G. Barucchi and G. Ponzano,Differential equations for one-loop generalized Feynman integrals,J. Math. Phys.14(1973) 396–401
1973
-
[15]
A. V. Kotikov,Differential equations method: New technique for massive Feynman diagrams calculation,Phys. Lett. B254(1991) 158–164
1991
-
[16]
A. V. Kotikov,Differential equations method: The Calculation of vertex type Feynman diagrams,Phys. Lett. B259(1991) 314–322
1991
-
[17]
T. Gehrmann and E. Remiddi,Differential equations for two loop four point functions,Nucl. Phys. B580(2000) 485–518, [hep-ph/9912329]
Pith/arXiv arXiv 2000
-
[18]
Z. Bern, L. J. Dixon, and D. A. Kosower,Dimensionally regulated pentagon integrals,Nucl. Phys. B412(1994) 751–816, [hep-ph/9306240]
Pith/arXiv arXiv 1994
-
[19]
J. M. Henn,Multiloop integrals in dimensional regularization made simple,Phys. Rev. Lett. 110(2013) 251601, [arXiv:1304.1806]
Pith/arXiv arXiv 2013
-
[20]
T. Peraro and L. Tancredi,Physical projectors for multi-leg helicity amplitudes,JHEP07 (2019) 114, [arXiv:1906.03298]
Pith/arXiv arXiv 2019
-
[21]
T. Peraro and L. Tancredi,Tensor decomposition for bosonic and fermionic scattering amplitudes,Phys. Rev. D103(2021), no. 5 054042, [arXiv:2012.00820]
Pith/arXiv arXiv 2021
-
[22]
’t Hooft and M
G. ’t Hooft and M. J. G. Veltman,Regularization and Renormalization of Gauge Fields, Nucl. Phys. B44(1972) 189–213
1972
-
[23]
F. V. Tkachov,A Theorem on Analytical Calculability of Four Loop Renormalization Group Functions,Phys. Lett.100B(1981) 65–68
1981
-
[24]
K. G. Chetyrkin and F. V. Tkachov,Integration by Parts: The Algorithm to Calculate beta Functions in 4 Loops,Nucl. Phys. B192(1981) 159–204. 9 Two-loop amplitude for𝑡¯𝑡𝑊productionMattia Pozzoli
1981
-
[25]
Laporta,High precision calculation of multiloop Feynman integrals by difference equations,Int
S. Laporta,High precision calculation of multiloop Feynman integrals by difference equations,Int. J. Mod. Phys.A15(2000) 5087–5159, [hep-ph/0102033]
Pith/arXiv arXiv 2000
-
[26]
Z. Wu, J. Boehm, R. Ma, H. Xu, and Y. Zhang,NeatIBP 1.0, a package generating small-size integration-by-parts relations for Feynman integrals,Comput. Phys. Commun.295(2024) 108999, [arXiv:2305.08783]
Pith/arXiv arXiv 2024
-
[27]
T. Peraro,FiniteFlow: multivariate functional reconstruction using finite fields and dataflow graphs,JHEP07(2019) 031, [arXiv:1905.08019]
Pith/arXiv arXiv 2019
-
[28]
N. Arkani-Hamed, J. L. Bourjaily, F. Cachazo, and J. Trnka,Local Integrals for Planar Scattering Amplitudes,JHEP06(2012) 125, [arXiv:1012.6032]
Pith/arXiv arXiv 2012
-
[29]
J. Henn, B. Mistlberger, V. A. Smirnov, and P. Wasser,Constructing d-log integrands and computing master integrals for three-loop four-particle scattering,JHEP04(2020) 167, [arXiv:2002.09492]
Pith/arXiv arXiv 2020
-
[30]
M. Becchetti, C. Dlapa, and S. Zoia,Canonical differential equations for the elliptic two-loop five-point integral family relevant to tt¯+jet production at leading color,Phys. Rev. D112(2025), no. 3 L031501, [arXiv:2503.03603]
Pith/arXiv arXiv 2025
-
[31]
S. Caron-Huot, M. Correia, and M. Giroux,Recursive Landau Analysis,Phys. Rev. Lett.135 (2025), no. 13 131603, [arXiv:2406.05241]
Pith/arXiv arXiv 2025
-
[32]
M. Correia, M. Giroux, and S. Mizera,SOFIA: Singularities of Feynman integrals automatized,Comput. Phys. Commun.320(2026) 109970, [arXiv:2503.16601]
arXiv 2026
-
[33]
S. Badger, M. Becchetti, N. Giraudo, and S. Zoia,Two-loop integrals for𝑡𝑡+jet production at hadron colliders in the leading colour approximation,JHEP07(2024) 073, [arXiv:2404.12325]
Pith/arXiv arXiv 2024
-
[34]
T. Gehrmann, J. M. Henn, and N. A. Lo Presti,Pentagon functions for massless planar scattering amplitudes,JHEP10(2018) 103, [arXiv:1807.09812]
Pith/arXiv arXiv 2018
-
[35]
D. Chicherin and V. Sotnikov,Pentagon Functions for Scattering of Five Massless Particles, JHEP20(2020) 167, [arXiv:2009.07803]
Pith/arXiv arXiv 2020
-
[36]
D. Chicherin, V. Sotnikov, and S. Zoia,Pentagon functions for one-mass planar scattering amplitudes,JHEP01(2022) 096, [arXiv:2110.10111]
Pith/arXiv arXiv 2022
-
[37]
S. Abreu, D. Chicherin, H. Ita, B. Page, V. Sotnikov, W. Tschernow, and S. Zoia,All Two-Loop Feynman Integrals for Five-Point One-Mass Scattering,Phys. Rev. Lett.132 (2024), no. 14 141601, [arXiv:2306.15431]
Pith/arXiv arXiv 2024
-
[38]
S. Badger, M. Becchetti, C. Brancaccio, H. B. Hartanto, and S. Zoia,Numerical evaluation of two-loop QCD helicity amplitudes for𝑔𝑔→𝑡𝑡𝑔at leading colour,JHEP03(2025) 070, [arXiv:2412.13876]. 10 Two-loop amplitude for𝑡¯𝑡𝑊productionMattia Pozzoli
Pith/arXiv arXiv 2025
-
[39]
S. Pozzorini and E. Remiddi,Precise numerical evaluation of the two loop sunrise graph master integrals in the equal mass case,Comput. Phys. Commun.175(2006) 381–387, [hep-ph/0505041]
Pith/arXiv arXiv 2006
-
[40]
F. Moriello,Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops,JHEP01(2020) 150, [arXiv:1907.13234]
Pith/arXiv arXiv 2020
-
[41]
X. Liu and Y.-Q. Ma,AMFlow: A Mathematica package for Feynman integrals computation viaauxiliarymassflow,Comput.Phys.Commun.283(2023)108565,[arXiv:2201.11669]
Pith/arXiv arXiv 2023
-
[42]
T. Peraro,Analytic multi-loop results using finite fields and dataflow graphs with FiniteFlow, in14th International Symposium on Radiative Corrections: Application of Quantum Field Theory to Phenomenology, 12, 2019.arXiv:1912.03142
Pith/arXiv arXiv 2019
-
[43]
A. von Manteuffel and R. M. Schabinger,A novel approach to integration by parts reduction, Phys. Lett. B744(2015) 101–104, [arXiv:1406.4513]
Pith/arXiv arXiv 2015
-
[44]
T. Peraro,Scattering amplitudes over finite fields and multivariate functional reconstruction, JHEP12(2016) 030, [arXiv:1608.01902]
Pith/arXiv arXiv 2016
-
[45]
S. Badger, H. B. Hartanto, and S. Zoia,Two-Loop QCD Corrections to Wbb¯Production at Hadron Colliders,Phys. Rev. Lett.127(2021), no. 1 012001, [arXiv:2102.02516]
Pith/arXiv arXiv 2021
-
[46]
Becchetti et al.,NNLO QCD predictions for𝑡¯𝑡𝑊production at hadron colliders, arXiv:2606.09503
M. Becchetti et al.,NNLO QCD predictions for𝑡¯𝑡𝑊production at hadron colliders, arXiv:2606.09503. 11
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.