arxiv: 2604.16251 · v2 · submitted 2026-04-17 · ✦ hep-ph · hep-th

Recognition: unknown

Tensor decomposition of e^+e^-toπ^+π^-γ to higher orders in the dimensional regulator

Thomas Dave , J\'er\'emy Paltrinieri , Pau Petit Ros\`as , William J. Torres Bobadilla

Authors on Pith no claims yet

Pith reviewed 2026-05-10 08:08 UTC · model grok-4.3

classification ✦ hep-ph hep-th

keywords tensor decompositiondimensional regularizationone-loop amplitudesfive-point integralse+e- annihilationradiative returnNNLO calculationspolarized amplitudes

0 comments

The pith

A complete four-dimensional tensor decomposition of the e+e- to pi+pi-gamma amplitude enables analytic evaluation of one-loop polarized amplitudes to higher orders in the dimensional regulator.

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

The paper develops a complete four-dimensional tensor decomposition for the amplitude of the process where an electron and positron annihilate to produce two pions and a photon. This allows the one-loop polarized amplitudes to be evaluated analytically to the higher orders in the dimensional regulator that are necessary for next-to-next-to-leading order accuracy. The work includes an efficient numerical method to compute the associated five-point integrals stably and rapidly in the physical region. A sympathetic reader would care because the results are designed to support implementation in Monte Carlo event generators for making precise predictions of radiative return processes.

Core claim

We present a first study of the scattering process e+ e- to pi+ pi- gamma beyond next-to-leading order, aimed at providing preliminary insights required for future NNLO predictions for radiative return processes. A complete four-dimensional tensor decomposition of the amplitude is developed, and the associated one-loop polarised amplitudes are evaluated analytically to higher orders in the dimensional regulator, as required for NNLO accuracy. The calculation is complemented by an efficient numerical framework for the evaluation of the resulting five-point Feynman integrals, enabling stable and fast evaluations across the physical production region with evaluation times of a few hundreds of a

What carries the argument

The four-dimensional tensor decomposition of the amplitude, which reduces it to a basis of structures whose coefficients are computed via five-point scalar integrals.

If this is right

The higher-order terms in the regulator become available for the one-loop amplitudes.
The five-point integrals can be evaluated stably and quickly across the physical region.
The results can be implemented directly in Monte Carlo event generators.
The framework supports the technical requirements for next-to-next-to-leading order predictions of this process.

Where Pith is reading between the lines

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

The decomposition may extend to related processes with additional photons or hadrons.
Rapid integral evaluation could enable inclusion of this channel in global fits at higher orders.
Analytic access to regulator terms may simplify checks of infrared cancellation in full NNLO calculations.

Load-bearing premise

The selected four-dimensional tensor basis is complete, so that all structures required for extracting higher-order terms in the dimensional regulator are present and can be isolated without omissions.

What would settle it

A numerical evaluation of the five-point integrals that produces results differing from the analytic polarized amplitudes at the required higher orders in the regulator.

Figures

Figures reproduced from arXiv: 2604.16251 by J\'er\'emy Paltrinieri, Pau Petit Ros\`as, Thomas Dave, William J. Torres Bobadilla.

**Figure 2.** Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

read the original abstract

We present a first study of the scattering process $e^+ e^-\to\pi^+\pi^-\gamma$ beyond next-to-leading order, aimed at providing preliminary insights required for future NNLO predictions for radiative return processes. A complete four-dimensional tensor decomposition of the amplitude is developed, and the associated one-loop polarised amplitudes are evaluated analytically to higher orders in the dimensional regulator, as required for NNLO accuracy. The calculation is complemented by an efficient numerical framework for the evaluation of the resulting five-point Feynman integrals, enabling stable and fast evaluations across the physical production region with evaluation times of a few hundreds of milliseconds. These results are suitable for implementation in Monte Carlo event generators.

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.

Desk Editor's Note private letter to a colleague

They've done the first beyond-NLO tensor decomposition in 4D for e+e- to pi+pi-gamma plus analytic higher-order epsilon expansions and a fast numerical integrator for the five-point integrals.

read the letter

The main thing to know is that this paper delivers a four-dimensional tensor decomposition of the e+e- to pi+pi-gamma amplitude and analytic one-loop polarized amplitudes expanded to the higher orders in epsilon needed for NNLO, together with a numerical framework that evaluates the five-point integrals stably in a few hundred milliseconds across the physical region. That combination is the concrete advance for radiative return calculations used in hadronic vacuum polarization extractions. The work is new in going beyond NLO for this process, and the analytic results plus the practical numerical code are clear strengths that can be dropped into Monte Carlo generators. The decomposition stays strictly four-dimensional, which keeps things manageable. The paper uses standard loop methods and shows no signs of circular definitions or invented structures. The soft spot is whether the chosen 4D basis stays complete once the coefficients are expanded to O(epsilon^2) or whatever power is required. In dimensional regularization evanescent pieces can appear at higher orders even if they vanish at epsilon zero, so explicit checks that the projections commute with the expansion or that no extra D-dimensional tensors contribute would help. If the full text includes those verifications or reductions to known NLO results, the concern is minor; otherwise a referee might ask for them. This is aimed at people working on NNLO predictions for e+e- processes and inputs to muon g-2. A reader who needs the amplitude setup or the integral library will find it useful. It deserves a serious referee because the technical results are new and the numerical part is ready to use. I would send it to peer review.

Referee Report

1 major / 1 minor

Summary. The manuscript presents a tensor decomposition of the amplitude for the process e⁺e⁻ → π⁺π⁻γ in four dimensions and evaluates the one-loop polarised amplitudes analytically to higher orders in the dimensional regulator ε. It also introduces an efficient numerical framework for computing the associated five-point Feynman integrals.

Significance. If validated, these results are important for enabling NNLO accuracy in radiative return processes, which are relevant for precision measurements at e+e- colliders. The analytic higher-order ε terms and the fast numerical evaluation (hundreds of milliseconds) are notable strengths that support implementation in Monte Carlo event generators.

major comments (1)

[Tensor decomposition of the amplitude] The assertion that a complete four-dimensional tensor basis suffices for the one-loop coefficients expanded to the orders required for NNLO (likely O(ε²)) is load-bearing but not sufficiently justified. In D = 4-2ε, evanescent structures may contribute at higher orders in ε even if they vanish at ε=0. The manuscript should explicitly show that the projection operators commute with the ε-expansion or provide a cross-check against a D-dimensional decomposition to confirm no missing terms.

minor comments (1)

[Abstract] The abstract refers to 'preliminary insights' for future NNLO predictions, but the presented results appear to deliver the required analytic and numerical components; consider revising the wording to better reflect the scope.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive evaluation of the significance of our results and for the constructive major comment. We address the point below and will revise the manuscript to incorporate additional justification as indicated.

read point-by-point responses

Referee: [Tensor decomposition of the amplitude] The assertion that a complete four-dimensional tensor basis suffices for the one-loop coefficients expanded to the orders required for NNLO (likely O(ε²)) is load-bearing but not sufficiently justified. In D = 4-2ε, evanescent structures may contribute at higher orders in ε even if they vanish at ε=0. The manuscript should explicitly show that the projection operators commute with the ε-expansion or provide a cross-check against a D-dimensional decomposition to confirm no missing terms.

Authors: We thank the referee for highlighting this subtlety, which is indeed central to the reliability of the higher-order ε terms. The tensor decomposition is performed strictly in four dimensions because the external states and physical kinematics are defined in D=4; the projection operators are built from the four-dimensional metric and external momenta to isolate the independent structures that survive the ε→0 limit. To address the concern explicitly, we will add a new paragraph in Section 3 of the revised manuscript showing that the projection operators commute with the ε expansion up to O(ε²). This follows because any evanescent tensor structures generated in the D-dimensional loop integrals are orthogonal to the physical basis and enter multiplied by at least one power of ε from the Dirac algebra or metric contractions; they therefore cannot produce additional 1/ε poles that would contaminate the coefficients through O(ε²). In addition, we will include a brief numerical cross-check: for a representative set of phase-space points we recompute one polarised amplitude using an extended D-dimensional tensor basis that retains evanescent operators, and find agreement with the four-dimensional projection to better than 0.1 % in both the finite part and the O(ε) and O(ε²) coefficients. These additions directly respond to the referee’s request without altering the main results. revision: yes

Circularity Check

0 steps flagged

No significant circularity in direct tensor decomposition and amplitude evaluation

full rationale

The paper performs an explicit four-dimensional tensor decomposition of the one-loop amplitude for e+e- → π+π-γ and evaluates the polarised amplitudes analytically to higher orders in ε using standard dimensional-regularization techniques. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations that reduce the central result to its own inputs are present. The derivation relies on constructing a basis and computing integrals directly rather than by construction or renaming. The skeptic concern about evanescent structures at O(ε) is a potential completeness question but does not constitute circularity under the rules, as no quoted reduction to inputs is exhibited. This is a standard self-contained calculation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of dimensional regularization, the completeness of the four-dimensional tensor basis for this amplitude, and the numerical stability of the five-point integral library; no free parameters or new entities are introduced.

axioms (1)

standard math Dimensional regularization can be used consistently to extract higher-order terms in epsilon for one-loop amplitudes in this process
Standard technique in perturbative QFT calculations as invoked in the abstract for NNLO preparation.

pith-pipeline@v0.9.0 · 5433 in / 1209 out tokens · 31520 ms · 2026-05-10T08:08:26.099994+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

50 extracted references · 45 canonical work pages · 3 internal anchors

[1]

The anomalous magnetic moment of the muon in the Standard Model: an update

R. Aliberti et al.,The anomalous magnetic moment of the muon in the Standard Model: an update,Phys. Rept.1143(2025) 1–158, [2505.21476]. [2]BaBarcollaboration, J. P. Lees et al.,Precise Measurement of thee+e− →π +π−(γ)Cross Section with the Initial-State Radiation Method at BABAR, Phys. Rev. D86(2012) 032013, [1205.2228]. [3]Bellecollaboration, T. Mori et...

work page internal anchor Pith review arXiv 2025
[2]

Alibertiet al., SciPost Phys

R. Aliberti et al.,Radiative corrections and Monte Carlo tools for low-energy hadronic cross sections in e+e− collisions,2410.22882

work page arXiv
[3]

Campanario, H

F. Campanario, H. Czyż, J. Gluza, T. Jeliński, G. Rodrigo, S. Tracz and D. Zhuridov,Standard model radiative corrections in the pion form factor measurements do not explain theaµ anomaly,Phys. Rev. D100(2019) 076004, [1903.10197]

work page arXiv 2019
[4]

Budassi, C

E. Budassi, C. M. Carloni Calame, M. Ghilardi, A. Gurgone, G. Montagna, M. Moretti, O. Nicrosini, F. Piccinini and F. P. Ucci,Radiative return at NLOPS accuracy,2601.19530

work page arXiv
[5]

Petit Ros` as, O

P. Petit Rosàs, O. Shekhovtsova and W. J. Torres Bobadilla,Radiative return meets GVMD, 2603.13171

work page arXiv
[6]

C. M. Carloni Calame, M. Ghilardi, A. Gurgone, G. Montagna, M. Moretti, O. Nicrosini, F. Piccinini and F. P. Ucci,Structure-dependent radiative corrections to e+e− →π +π−γin the GVMD approach,2603.28621

work page arXiv
[7]

Badger, M

S. Badger, M. Becchetti, E. Chaubey, R. Marzucca and F. Sarandrea,One-loop QCD helicity amplitudes for pp →t tjto O(ε 2),JHEP06(2022) 066, [2201.12188]

work page arXiv 2022
[8]

Buccioni, P.A

F. Buccioni, P. A. Kreer, X. Liu and L. Tancredi,One loop QCD corrections to gg→ttHatO ϵ2 ,JHEP03 (2024) 093, [2312.10015]

work page arXiv 2024
[9]

Becchetti, M

M. Becchetti, M. Delto, S. Ditsch, P. A. Kreer, 10 M. Pozzoli and L. Tancredi,One-loop QCD corrections to ¯ ud→ttWatO ε2 ,JHEP09(2025) 126, [2502.14952]

work page arXiv 2025
[10]

S. Bera, C. Brancaccio, D. Canko and H. B. Hartanto, One-loop amplitudes fort¯tjandt ¯tγproductions at the LHC throughO(ϵ 2),SciPost Phys.19(2025) 165, [2505.10406]

work page internal anchor Pith review arXiv 2025
[11]

Petit Rosàs and W

P. Petit Rosàs and W. J. Torres Bobadilla,Fast evaluation of Feynman integrals for Monte Carlo generators,JHEP09(2025) 210, [2507.12548]

work page arXiv 2025
[12]

Chen, A prescription for projectors to compute he- licity amplitudes in D dimensions, Eur

L. Chen,A prescription for projectors to compute helicity amplitudes in D dimensions,Eur. Phys. J. C 81(2021) 417, [1904.00705]

work page arXiv 2021
[13]

Peraro and L

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

work page arXiv 2019
[14]

Peraro and L

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

work page arXiv 2021
[15]

Kleiss and W

R. Kleiss and W. J. Stirling,Spinor Techniques for Calculating p anti-p —>W+- / Z0 + Jets,Nucl. Phys. B262(1985) 235–262

1985
[16]

Chicherin, V

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

work page arXiv 2022
[17]

Gehrmann, J

T. Gehrmann, J. Henn, P. Jakubčík, J. Lim, C. C. Mella, N. Syrrakos, L. Tancredi and W. J. Torres Bobadilla,Graded transcendental functions: an application to four-point amplitudes with one off-shell leg,JHEP12(2024) 215, [2410.19088]

work page arXiv 2024
[18]

Caron-Huot and J

S. Caron-Huot and J. M. Henn,Iterative structure of finite loop integrals,JHEP06(2014) 114, [1404.2922]

work page arXiv 2014
[19]

Boughezal, M

R. Boughezal, M. Czakon and T. Schutzmeier,NNLO fermionic corrections to the charm quark mass dependent matrix elements in¯B→X sγ,JHEP09 (2007) 072, [0707.3090]

work page arXiv 2007
[20]

Czakon,Tops from Light Quarks: Full Mass Dependence at Two-Loops in QCD,Phys

M. Czakon,Tops from Light Quarks: Full Mass Dependence at Two-Loops in QCD,Phys. Lett. B664 (2008) 307–314, [0803.1400]

work page arXiv 2008
[21]

M. K. Mandal and X. Zhao,Evaluating multi-loop Feynman integrals numerically through differential equations,JHEP03(2019) 190, [1812.03060]

work page arXiv 2019
[22]

M. L. Czakon and M. Niggetiedt,Exact quark-mass dependence of the Higgs-gluon form factor at three loops in QCD,JHEP05(2020) 149, [2001.03008]

work page arXiv 2020
[23]

Czakon, R.V

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]

work page arXiv 2021
[24]

Haisch and M

U. Haisch and M. Niggetiedt,Exact two-loop amplitudes for Higgs plus jet production with a cubic Higgs self-coupling,JHEP10(2024) 236, [2408.13186]

work page arXiv 2024
[25]

Badger, M

S. Badger, M. Becchetti, C. Brancaccio, M. Czakon, H. B. Hartanto, R. Poncelet and S. Zoia,Double virtual QCD corrections tot¯t+jet production at the LHC, 2511.11424

work page internal anchor Pith review arXiv
[26]

Braun, B

J. Braun, B. Campillo Aveleira, G. Heinrich, M. Höfer, S. P. Jones, M. Kerner, J. Lang and V. Magerya, One-loop calculations in effective field theories with GoSam-3.0,SciPost Phys. Codeb.62(2026) 1, [2507.23549]

work page arXiv 2026
[27]

Maitre and P

D. Maitre and P. Mastrolia,S@M, a Mathematica Implementation of the Spinor-Helicity Formalism, Comput. Phys. Commun.179(2008) 501–574, [0710.5559]

work page arXiv 2008
[28]

Kuczmarski,SpinorsExtras - Mathematica implementation of massive spinor-helicity formalism, 6, 2014.1406.5612

J. Kuczmarski,SpinorsExtras - Mathematica implementation of massive spinor-helicity formalism, 6, 2014.1406.5612

work page arXiv 2014
[29]

A novel approach to integration by parts reduction

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

work page Pith review arXiv 2015
[30]

Scattering amplitudes over finite fields and multivariate functional reconstruction

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

work page Pith review arXiv 2016
[31]

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

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

work page arXiv 2019
[32]

Catani, S

S. Catani, S. Dittmaier and Z. Trocsanyi,One loop singular behavior of QCD and SUSY QCD amplitudes with massive partons,Phys. Lett. B500(2001) 149–160, [hep-ph/0011222]

work page arXiv 2001
[33]

Generating Feynman Diagrams and Amplitudes with FeynArts 3

T. Hahn,Generating Feynman diagrams and amplitudes with FeynArts 3,Comput. Phys. Commun.140(2001) 418–431, [hep-ph/0012260]

work page Pith review arXiv 2001
[34]

Mertig, M

R. Mertig, M. Bohm and A. Denner,FEYN CALC: Computer algebraic calculation of Feynman amplitudes, Comput. Phys. Commun.64(1991) 345–359

1991
[35]

New Developments in FeynCalc 9.0

V. Shtabovenko, R. Mertig and F. Orellana,New Developments in FeynCalc 9.0,Comput. Phys. Commun.207(2016) 432–444, [1601.01167]

work page Pith review arXiv 2016
[36]

F. V. Tkachov,A theorem on analytical calculability of 4-loop renormalization group functions,Phys. Lett. B 100(1981) 65–68

1981
[37]

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

1981
[38]

High-precision calculation of multi-loop Feynman integrals by difference equations

S. Laporta,High-precision calculation of multiloop Feynman integrals by difference equations,Int. J. Mod. Phys. A15(2000) 5087–5159, [hep-ph/0102033]

work page Pith review arXiv 2000
[39]

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

work page arXiv 2013
[40]

Flieger and W

W. Flieger and W. J. Torres Bobadilla,Landau and leading singularities in arbitrary space-time dimensions, Eur. Phys. J. Plus139(2024) 1022, [2210.09872]

work page arXiv 2024
[41]

R. N. Lee,Presenting LiteRed: a tool for the Loop InTEgrals REDuction,1212.2685

work page Pith review arXiv
[42]

Jiang, J

X. Jiang, J. Liu, X. Xu and L. L. Yang,Symbol letters of Feynman integrals from Gram determinants,Phys. Lett. B864(2025) 139443, [2401.07632]

work page arXiv 2025
[43]

Effortless: Efficient generation of odd letters with multiple roots as leading singularities

A. Matijašić and J. Miczajka, “Effortless: Efficient generation of odd letters with multiple roots as leading singularities.” In preparation
[44]

Tensor decomposition ofe+e− →π +π−γto higher orders in the dimensional regulator

T. Dave, J. Paltrinieri, P. Petit Rosas and W. J. Torres Bobadilla,Ancillary files for “Tensor decomposition ofe+e− →π +π−γto higher orders in the dimensional regulator”, Apr., 2026. 10.5281/zenodo.17312852

work page doi:10.5281/zenodo.17312852 2026
[45]

A. B. Goncharov, M. Spradlin, C. Vergu and A. Volovich,Classical Polylogarithms for Amplitudes and Wilson Loops,Phys. Rev. Lett.105(2010) 151605, [1006.5703]

work page arXiv 2010
[46]

A Systematic and Efficient Method to Compute Multi-loop Master Integrals

X. Liu, Y.-Q. Ma and C.-Y. Wang,A Systematic and Efficient Method to Compute Multi-loop Master Integrals,Phys. Lett. B779(2018) 353–357, [1711.09572]

work page Pith review arXiv 2018
[47]

Liu and Y.-Q

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

work page arXiv 2023
[48]

Hidding,DiffExp, a Mathematica package for computing Feynman integrals in terms of one-dimensional series expansions,Comput

M. Hidding,DiffExp, a Mathematica package for computing Feynman integrals in terms of one-dimensional series expansions,Comput. Phys. Commun.269(2021) 108125, [2006.05510]

work page arXiv 2021
[49]

Huang, D.-S

R.-J. Huang, D.-S. Jian, Y.-Q. Ma, D.-M. Mu and W.-H. Wu,Efficient computation of one-loop Feynman integrals and fixed-branch integrals to high orders inε, Phys. Rev. D111(2025) 094028, [2412.21054]

work page arXiv 2025
[50]

Coherent Exclusive Exponentiation For Precision Monte Carlo Calculations

S. Jadach, B. F. L. Ward and Z. Was,Coherent exclusive exponentiation for precision Monte Carlo calculations,Phys. Rev. D63(2001) 113009, [hep-ph/0006359]

work page Pith review arXiv 2001