Pith. sign in

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 →

arxiv 2608.02527 v1 pith:OM2ELYHA submitted 2026-08-03 hep-ph

Two-loop amplitude for tbar{t}W production at hadron colliders in the leading colour approximation

classification hep-ph
keywords two-loop amplitudesttbarW productionleading colourQCD correctionselliptic integralsdifferential equationsfinite-field reconstructionNNLO
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper reports the first exact numerical evaluation of the two-loop QCD amplitude for associated top-antitop-plus-W production, in the leading-colour approximation. This process is an important LHC background and a probe of new physics, and current NNLO predictions rely on an approximation combining the soft-W limit with massification, so an exact result is needed to validate that approximation and to reach differential distributions. The computation reduces 8,959 Feynman integrals to 330 master integrals, expresses the amplitude as a set of special functions that includes elliptic and nested-square-root sectors, and evaluates these functions numerically via differential equations and generalised series expansions while reconstructing rational coefficients from finite-field samples. The author reports that the finite remainder can be evaluated at a phase-space point in about one hour, and that an interpolation grid built from this evaluation has already been used for NNLO inclusive cross-section results.

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.

Watch this falsifier — get emailed when new claim-graph text bears on it.

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

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

  • 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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 3 minor

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)
  1. [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.
  2. [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.
  3. [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)
  1. [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.
  2. [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.
  3. [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

0 steps flagged

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

0 free parameters · 5 axioms · 0 invented entities

The paper introduces no new physical entities. Its assumptions are standard QCD regularisation, the completeness of the IBP reduction, the correctness of the DEs, and two known but unproved conjectures used in the special-function construction. The absence of shipped equations and reduction tables means these assumptions are not independently verifiable from this paper.

axioms (5)
  • domain assumption The 't Hooft-Veltman (tHV) scheme and dimensional regularisation in D=4-2ε provide a consistent regulator for the amplitude.
    Standard in perturbative QCD; invoked in Section 2. The final finite remainder is scheme-dependent and must be matched to the IR subtraction used in [46].
  • domain assumption The 330 master integrals obtained via NeatIBP and FiniteFlow form a complete basis and the IBPs are correct.
    Section 3. The reduction of the 8959 integrals to 330 MIs is a computational step whose correctness is assumed; no IBPs or reduction tables are shown.
  • domain assumption The differential equations for the MIs and for the special functions correctly encode the analytic dependence of the integrals.
    Section 3.1 and 4. The DEs are stated in general form but not given; the solution via generalized series expansions in AMFlow is assumed to produce the correct Laurent coefficients.
  • ad hoc to paper The conjecture that integrands with at most simple poles and constant leading singularity lead to canonical Feynman integrals.
    Section 3.1, eq. (11) and surrounding text. This is an established but unproved conjecture in the field, and the paper explicitly relies on it for constructing the polylogarithmic canonical basis.
  • ad hoc to paper The algebraic-independence conjecture that the only transcendental constants appearing are zeta-values ζ_n.
    Section 4, citing [37]: 'conjectured to be only the zeta-values ζ_n'. This conjecture underlies the pole subtraction and the special-function basis construction.

pith-pipeline@v1.3.0-daily-deepseek · 7791 in / 7439 out tokens · 61454 ms · 2026-08-04T05:20:46.594078+00:00 · methodology

0 comments
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

Figures reproduced from arXiv: 2608.02527 by Mattia Pozzoli.

Figure 1
Figure 1. Figure 1: One-loop (fig. 1a) and two-loop (figs. 1b to 1d) integral families contributing to the leading-colour amplitude. Thin black lines denote massless particles, thick red lines indicate the top quarks and the curly blue line is the 𝑊 boson. In order to regulate the divergences of the Feynman integrals we use dimensional regularisation, working in 𝐷 = 4 − 2𝜀 space-time dimensions. The amplitude admits an expans… view at source ↗
Figure 2
Figure 2. Figure 2: Graphs of the sectors associated with the elliptic curves (figs. 2a to 2c) and with the nested square root (fig. 2d). integrals appearing in the families in figs. 1b to 1d fall under this case, and we construct a basis of canonical MIs for them. However, eq. (11) could reveal also differential forms more complicated than the logarithmic ones, associated for example with an elliptic curve d𝑧 √︁ P4 (𝑧) , wit… view at source ↗

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

44 extracted references · 36 linked inside Pith

  1. [1]

    Buckley, C

    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]

  2. [2]

    J. A. Dror, M. Farina, E. Salvioni, and J. Serra,Strong tW Scattering at the LHC,JHEP01 (2016) 071, [arXiv:1511.03674]

  3. [3]

    Bessidskaia Bylund, F

    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...

  4. [6]

    Buonocore, S

    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]

  5. [7]

    Catani, S

    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]

  6. [8]

    Bärnreuther, M

    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

  7. [9]

    A. A. Penin,Two-loop photonic corrections to massive Bhabha scattering,Nucl. Phys. B734 (2006) 185–202, [hep-ph/0508127]

  8. [10]

    Mitov and S

    A. Mitov and S. Moch,The Singular behavior of massive QCD amplitudes,JHEP05(2007) 001, [hep-ph/0612149]

  9. [11]

    Becher and K

    T. Becher and K. Melnikov,Two-loop QED corrections to Bhabha scattering,JHEP06 (2007) 084, [arXiv:0704.3582]

  10. [12]

    Becchetti, M

    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]

  11. [13]

    Becchetti, D

    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]

  12. [14]

    Barucchi and G

    G. Barucchi and G. Ponzano,Differential equations for one-loop generalized Feynman integrals,J. Math. Phys.14(1973) 396–401

  13. [15]

    A. V. Kotikov,Differential equations method: New technique for massive Feynman diagrams calculation,Phys. Lett. B254(1991) 158–164

  14. [16]

    A. V. Kotikov,Differential equations method: The Calculation of vertex type Feynman diagrams,Phys. Lett. B259(1991) 314–322

  15. [17]

    Gehrmann and E

    T. Gehrmann and E. Remiddi,Differential equations for two loop four point functions,Nucl. Phys. B580(2000) 485–518, [hep-ph/9912329]

  16. [18]

    Z. Bern, L. J. Dixon, and D. A. Kosower,Dimensionally regulated pentagon integrals,Nucl. Phys. B412(1994) 751–816, [hep-ph/9306240]

  17. [19]

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

  18. [20]

    Peraro and L

    T. Peraro and L. Tancredi,Physical projectors for multi-leg helicity amplitudes,JHEP07 (2019) 114, [arXiv:1906.03298]

  19. [21]

    Peraro and L

    T. Peraro and L. Tancredi,Tensor decomposition for bosonic and fermionic scattering amplitudes,Phys. Rev. D103(2021), no. 5 054042, [arXiv:2012.00820]

  20. [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

  21. [23]

    F. V. Tkachov,A Theorem on Analytical Calculability of Four Loop Renormalization Group Functions,Phys. Lett.100B(1981) 65–68

  22. [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

  23. [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]

  24. [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]

  25. [27]

    Peraro,FiniteFlow: multivariate functional reconstruction using finite fields and dataflow graphs,JHEP07(2019) 031, [arXiv:1905.08019]

    T. Peraro,FiniteFlow: multivariate functional reconstruction using finite fields and dataflow graphs,JHEP07(2019) 031, [arXiv:1905.08019]

  26. [28]

    Arkani-Hamed, J

    N. Arkani-Hamed, J. L. Bourjaily, F. Cachazo, and J. Trnka,Local Integrals for Planar Scattering Amplitudes,JHEP06(2012) 125, [arXiv:1012.6032]

  27. [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]

  28. [30]

    Becchetti, C

    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]

  29. [31]

    Caron-Huot, M

    S. Caron-Huot, M. Correia, and M. Giroux,Recursive Landau Analysis,Phys. Rev. Lett.135 (2025), no. 13 131603, [arXiv:2406.05241]

  30. [32]

    Correia, M

    M. Correia, M. Giroux, and S. Mizera,SOFIA: Singularities of Feynman integrals automatized,Comput. Phys. Commun.320(2026) 109970, [arXiv:2503.16601]

  31. [33]

    Badger, M

    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]

  32. [34]

    Gehrmann, J

    T. Gehrmann, J. M. Henn, and N. A. Lo Presti,Pentagon functions for massless planar scattering amplitudes,JHEP10(2018) 103, [arXiv:1807.09812]

  33. [35]

    Chicherin and V

    D. Chicherin and V. Sotnikov,Pentagon Functions for Scattering of Five Massless Particles, JHEP20(2020) 167, [arXiv:2009.07803]

  34. [36]

    Chicherin, V

    D. Chicherin, V. Sotnikov, and S. Zoia,Pentagon functions for one-mass planar scattering amplitudes,JHEP01(2022) 096, [arXiv:2110.10111]

  35. [37]

    Abreu, D

    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]

  36. [38]

    Badger, M

    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

  37. [39]

    Pozzorini and E

    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]

  38. [40]

    Moriello,Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops,JHEP01(2020) 150, [arXiv:1907.13234]

    F. Moriello,Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops,JHEP01(2020) 150, [arXiv:1907.13234]

  39. [41]

    Liu and Y.-Q

    X. Liu and Y.-Q. Ma,AMFlow: A Mathematica package for Feynman integrals computation viaauxiliarymassflow,Comput.Phys.Commun.283(2023)108565,[arXiv:2201.11669]

  40. [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

  41. [43]

    von Manteuffel and R

    A. von Manteuffel and R. M. Schabinger,A novel approach to integration by parts reduction, Phys. Lett. B744(2015) 101–104, [arXiv:1406.4513]

  42. [44]

    Peraro,Scattering amplitudes over finite fields and multivariate functional reconstruction, JHEP12(2016) 030, [arXiv:1608.01902]

    T. Peraro,Scattering amplitudes over finite fields and multivariate functional reconstruction, JHEP12(2016) 030, [arXiv:1608.01902]

  43. [45]

    Badger, H

    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]

  44. [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