Chebyshev Approximations of Feynman Integrals for Collider Physics
Pith reviewed 2026-07-03 09:37 UTC · model grok-4.3
The pith
Chebyshev polynomial approximations along paths solve canonical differential equations for Feynman integrals with rapid convergence and high efficiency.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By constructing Chebyshev polynomial approximations along a path and using an adaptive method to sample for optimal convergence, the integrals' canonical differential equations can be solved numerically in a stable manner across physical phase space; the resulting implementation requires little case-by-case handling of spurious singularities and performs competitively with state-of-the-art one-fold integral methods for two-loop five-point examples.
What carries the argument
Chebyshev polynomial approximations along a path that exploit analyticity properties to produce rapidly converging representations for numerical evaluation.
If this is right
- Numerical evaluation becomes highly efficient once the polynomial fit is built.
- Stability holds across physical phase space in the tested two-loop five-point cases.
- Little to no manual intervention is needed to handle spurious singularities.
- The framework runs entirely in double-precision arithmetic.
- Performance is competitive with current one-fold integral methods for the examined examples.
Where Pith is reading between the lines
- The same path-based approximation could be applied to integrals at higher loop orders if analyticity continues to hold.
- Automation of path selection might reduce remaining manual choices when moving to new integral families.
- The method could be combined with existing reduction tools to form a more complete pipeline for multi-loop amplitudes.
- Testing on integrals with branch cuts that cross the chosen path would clarify the limits of the current adaptive sampling.
Load-bearing premise
The method depends on Feynman integrals possessing analyticity properties that permit construction of rapidly converging polynomial approximations along suitable paths.
What would settle it
Direct numerical comparison of the Chebyshev-approximated values against independent high-precision results for a known two-loop five-point integral at a physical kinematic point would show whether the claimed stability and accuracy hold.
read the original abstract
We present a novel approach for solving canonical differential equations for Feynman integrals based on an approximation of the integrals with Chebyshev polynomials. By exploiting the analyticity properties of Feynman integrals, the method constructs rapidly converging polynomial approximations along a path, enabling highly efficient numerical evaluation. Moreover, we introduce an adaptive approximation method that dynamically samples to optimise convergence. We implement this framework in double-precision arithmetic and demonstrate its stability across physical phase space using a series of two-loop, five-point examples. Our proof-of-principle implementation proves competitive with state-of-the-art one-fold integral methods, while requiring little to no case-by-case intervention to handle spurious singularities.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a method for numerically evaluating Feynman integrals by constructing Chebyshev polynomial approximations along paths in phase space, exploiting the integrals' analyticity properties to solve canonical differential equations. An adaptive sampling technique optimizes convergence. Implemented in double precision, the approach is tested on two-loop five-point integrals, claiming stability across physical phase space and competitiveness with state-of-the-art one-fold integral methods while requiring minimal case-by-case intervention for spurious singularities.
Significance. If the performance and stability claims are quantitatively validated, the method could serve as an efficient, automatable tool for multi-loop integral evaluations in collider physics. The adaptation of standard Chebyshev techniques to exploit analyticity along paths, combined with adaptive sampling, represents a practical numerical advance that may reduce reliance on specialized handling of singularities in higher-point calculations.
major comments (2)
- [Abstract] Abstract: the central claims that the implementation 'proves competitive with state-of-the-art one-fold integral methods' and demonstrates 'stability across physical phase space' are asserted without any quantitative benchmarks, error estimates, CPU timings, or direct comparisons to existing methods. This information is load-bearing for assessing the performance assertions and must be supplied with explicit data from the results.
- [Numerical examples] Numerical examples section: while two-loop five-point cases are used to illustrate stability, the absence of tabulated error metrics, convergence rates, or side-by-side comparisons with one-fold integral methods prevents independent verification of the competitiveness claim.
minor comments (1)
- [Method] The description of the adaptive sampler would benefit from an explicit statement of the convergence criterion (e.g., a threshold on the Chebyshev coefficient decay) to ensure reproducibility.
Simulated Author's Rebuttal
We thank the referee for the careful review and for highlighting the need for explicit quantitative support of our performance claims. We address each major comment below and will revise the manuscript accordingly.
read point-by-point responses
-
Referee: [Abstract] Abstract: the central claims that the implementation 'proves competitive with state-of-the-art one-fold integral methods' and demonstrates 'stability across physical phase space' are asserted without any quantitative benchmarks, error estimates, CPU timings, or direct comparisons to existing methods. This information is load-bearing for assessing the performance assertions and must be supplied with explicit data from the results.
Authors: We agree that the abstract's performance assertions require explicit supporting data. The current manuscript contains numerical examples in Section 4 that illustrate stability on two-loop five-point integrals, but direct side-by-side timings and error comparisons with one-fold integral methods are not tabulated. In the revised manuscript we will add a dedicated results table (or subsection) reporting maximum absolute errors, relative errors, CPU timings on representative phase-space points, and direct comparisons against existing one-fold integral implementations, with the abstract updated to reference these quantitative results. revision: yes
-
Referee: [Numerical examples] Numerical examples section: while two-loop five-point cases are used to illustrate stability, the absence of tabulated error metrics, convergence rates, or side-by-side comparisons with one-fold integral methods prevents independent verification of the competitiveness claim.
Authors: We accept this assessment. Although the examples demonstrate that the Chebyshev approximations remain stable across physical phase space without case-by-case singularity handling, the section does not include the requested tabulated metrics or comparisons. The revised manuscript will expand the numerical examples section with tables of error metrics versus number of Chebyshev nodes, convergence rates, CPU timings, and direct numerical comparisons to state-of-the-art one-fold integral methods on the same phase-space points. revision: yes
Circularity Check
No significant circularity
full rationale
The paper presents a numerical technique that approximates Feynman integrals via Chebyshev polynomials along paths chosen to exploit their analyticity properties, with an adaptive sampler for convergence. This is a standard construction in numerical analysis applied to the domain of canonical differential equations for integrals; the performance claims rest on explicit implementation and benchmarking against two-loop five-point examples rather than any self-referential definition, fitted parameter renamed as prediction, or load-bearing self-citation chain. No equations or steps in the provided description reduce the claimed result to its own inputs by construction.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Feynman integrals are analytic in appropriate regions of phase space, permitting polynomial approximation along paths.
Reference graph
Works this paper leans on
-
[1]
Heinrich,Collider physics at the precision frontier,Physics Reports922(Aug., 2021)
G. Heinrich,Collider physics at the precision frontier,Physics Reports922(Aug., 2021)
2021
-
[2]
A. Huss, J. Huston, S. Jones, M. Pellen and R. R¨ ontsch,Les houches 2023 – physics at tev colliders: Report on the standard model precision wishlist,SciPost Physics Community Reports(Mar., 2026)
2023
-
[3]
Aliberti, T
R. Aliberti, T. Aoyama, E. Balzani, A. Bashir, G. Benton, J. Bijnens et al.,The anomalous magnetic moment of the muon in the standard model: an update,Physics Reports1143 (Nov., 2025) 1–158. – 26 –
2025
-
[4]
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
-
[5]
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 internal anchor Pith review Pith/arXiv arXiv 2000
- [6]
- [7]
-
[8]
Reduze 2 - Distributed Feynman Integral Reduction
A. von Manteuffel and C. Studerus,Reduze 2 - Distributed Feynman Integral Reduction, 1201.4330
work page internal anchor Pith review Pith/arXiv arXiv
- [9]
-
[10]
R. N. Lee,LiteRed 1.4: a powerful tool for reduction of multiloop integrals,J. Phys. Conf. Ser.523(2014) 012059, [1310.1145]
work page internal anchor Pith review Pith/arXiv arXiv 2014
- [11]
-
[12]
V. A. Smirnov,Feynman integral calculus. Springer, Berlin, Heidelberg, 2006
2006
-
[13]
S. Weinzierl,Feynman Integrals. A Comprehensive Treatment for Students and Researchers. UNITEXT for Physics. Springer, 2022. 10.1007/978-3-030-99558-4
-
[14]
A. V. Kotikov,Differential equations method: New technique for massive Feynman diagrams calculation,Phys. Lett. B254(1991) 158–164
1991
-
[15]
A. V. Kotikov,Differential equation method: The Calculation of N point Feynman diagrams, Phys. Lett. B267(1991) 123–127
1991
-
[16]
Z. Bern, L. J. Dixon and D. A. Kosower,Dimensionally regulated pentagon integrals,Nucl. Phys. B412(1994) 751–816, [hep-ph/9306240]
work page internal anchor Pith review Pith/arXiv arXiv 1994
-
[17]
Remiddi,Differential equations for feynman graph amplitudes,Il Nuovo Cimento A110 (Dec., 1997) 1435–1452
E. Remiddi,Differential equations for feynman graph amplitudes,Il Nuovo Cimento A110 (Dec., 1997) 1435–1452
1997
-
[18]
Gehrmann and E
T. Gehrmann and E. Remiddi,Differential equations for two-loop four-point functions, Nuclear Physics B580(2000) 485–518
2000
-
[19]
J. M. Henn,Multiloop integrals in dimensional regularization made simple,Phys. Rev. Lett. 110(2013) 251601, [1304.1806]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[20]
R. N. Lee, A. V. Smirnov and V. A. Smirnov,Solving differential equations for feynman integrals by expansions near singular points,Journal of High Energy Physics2018(Mar.,
-
[21]
F. Moriello,Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops,JHEP01(2020) 150, [1907.13234]
-
[22]
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]. – 27 –
-
[23]
T. Armadillo, R. Bonciani, S. Devoto, N. Rana and A. Vicini,Evaluation of Feynman integrals with arbitrary complex masses via series expansions,Comput. Phys. Commun.282 (2023) 108545, [2205.03345]
-
[24]
Liu, Y.-Q
X. Liu, Y.-Q. Ma and C.-Y. Wang,A systematic and efficient method to compute multi-loop master integrals,Physics Letters B779(Apr., 2018) 353–357
2018
-
[25]
Liu and Y.-Q
X. Liu and Y.-Q. Ma,Multiloop corrections for collider processes using auxiliary mass flow, Physical Review D105(Mar., 2022)
2022
-
[26]
X. Liu and Y.-Q. Ma,AMFlow: A Mathematica package for Feynman integrals computation via auxiliary mass flow,Comput. Phys. Commun.283(2023) 108565, [2201.11669]
- [27]
-
[28]
Double virtual QCD corrections to $t\bar{t}+$jet production at the LHC
S. Badger, M. Becchetti, C. Brancaccio, M. Czakon, H. B. Hartanto, R. Poncelet et al., Double virtual QCD corrections tot t+jet production at the LHC,JHEP05(2026) 044, [2511.11424]
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [29]
-
[30]
P. Petit Ros` as and W. J. Torres Bobadilla,Fast evaluation of Feynman integrals for Monte Carlo generators,JHEP09(2025) 210, [2507.12548]
-
[31]
Solution of Canonical Differential Equations for Integrals on Arbitrary Geometries
M. Czakon and L. Tancredi,Solution of Canonical Differential Equations for Integrals on Arbitrary Geometries,2606.30354
work page internal anchor Pith review Pith/arXiv arXiv
-
[32]
IterInt: Evaluating iterated integrals via differential equations
G. Baur and C. Duhr,IterInt: Evaluating iterated integrals via differential equations, 2606.02744
work page internal anchor Pith review Pith/arXiv arXiv
- [33]
-
[34]
C. Canuto, M. Y. Hussaini, A. Quarteroni and T. A. Zang,Spectral Methods in Fluid Dynamics. 1988. 10.1007/978-3-642-84108-8
-
[35]
Canuto, M
C. Canuto, M. Hussaini, A. Quarteroni and T. Zang,Spectral Methods: Fundamentals in Single Domains. Scientific Computation. Springer Berlin Heidelberg, 2007
2007
-
[36]
Temperton,On scalar and vector transform methods for global spectral models,Monthly Weather Review119(1991) 1303 – 1307
C. Temperton,On scalar and vector transform methods for global spectral models,Monthly Weather Review119(1991) 1303 – 1307
1991
-
[37]
Grandcl´ ement and J
P. Grandcl´ ement and J. Novak,Spectral methods for numerical relativity,Living Reviews in Relativity12(Jan., 2009)
2009
-
[38]
M. Helmer and F. Tellander,Geometric singularities of Feynman integrals,Phys. Rev. D 112(2025) 065019, [2506.05042]
-
[39]
S. Mizera and S. Telen,Landau discriminants,JHEP08(2022) 200, [2109.08036]
- [40]
- [41]
-
[42]
M. Correia, M. Giroux and S. Mizera,SOFIA: Singularities of Feynman integrals automatized,Comput. Phys. Commun.320(2026) 109970, [2503.16601]
-
[43]
Landau Singularities from Whitney Stratifications
M. Helmer, G. Papathanasiou and F. Tellander,Landau Singularities from Whitney Stratifications,2402.14787
work page internal anchor Pith review Pith/arXiv arXiv
- [44]
-
[45]
T. A. Driscoll, N. Hale and L. N. Trefethen, eds.,Chebfun Guide. Pafnuty Publications, Oxford, 2014
2014
-
[46]
Z. Battles and L. N. Trefethen,An extension of matlab to continuous functions and operators,SIAM Journal on Scientific Computing25(2004) 1743–1770, [https://doi.org/10.1137/S1064827503430126]
-
[47]
J. L. Aurentz and L. N. Trefethen,Chopping a Chebyshev Series,arXiv e-prints(Dec., 2015) arXiv:1512.01803, [1512.01803]
work page internal anchor Pith review Pith/arXiv arXiv 2015
-
[48]
D. Chicherin and V. Sotnikov,Pentagon Functions for Scattering of Five Massless Particles, JHEP20(2020) 167, [2009.07803]
-
[49]
D. Chicherin, V. Sotnikov and S. Zoia,Pentagon functions for one-mass planar scattering amplitudes,JHEP01(2022) 096, [2110.10111]
- [50]
-
[51]
L. N. Trefethen,Approximation Theory and Approximation Practice, Extended Edition. Society for Industrial and Applied Mathematics, Philadelphia, PA, 2019. 10.1137/1.9781611975949
-
[52]
J. Mason and D. C. Handscomb,Chebyshev Polynomials. Chapman and Hall/CRC, 2002. 10.1201/9781420036114
-
[53]
Bernstein,Sur l’ordre de la meilleure approximation des fonctions continues par des polynˆ omes de degr´ e donn´ e,M´ emoires de l’Acad´ emie royale de Belgique4(1912) 1–104
S. Bernstein,Sur l’ordre de la meilleure approximation des fonctions continues par des polynˆ omes de degr´ e donn´ e,M´ emoires de l’Acad´ emie royale de Belgique4(1912) 1–104
1912
-
[54]
Bernstein,Sur la meilleure approximation de|x|par des polynomes de degr´ es donn´ es,Acta Mathematica37(1914) 1–57
S. Bernstein,Sur la meilleure approximation de|x|par des polynomes de degr´ es donn´ es,Acta Mathematica37(1914) 1–57
1914
-
[55]
J. P. Boyd,Chebyshev and Fourier Spectral Methods. Courier Corporation, 2013
2013
-
[56]
C. W. Clenshaw,A note on the summation of chebyshev series,Mathematics of Computation 9(1955) 118–120
1955
-
[57]
Cheney and D
E. Cheney and D. Kincaid,Numerical Mathematics and Computing. Cengage Learning, 2012
2012
-
[58]
W. H. Press, ed.,Numerical Recipes: The Art of Scientific Computing. Cambridge University Press, 3 ed., 2007
2007
-
[59]
Salvatori,The Tropical Geometry of Subtraction Schemes,2406.14606
G. Salvatori,The Tropical Geometry of Subtraction Schemes,2406.14606
-
[60]
M. Giroux, S. Mizera and G. Salvatori,SubTropica,2604.20954
work page internal anchor Pith review Pith/arXiv arXiv
-
[61]
Borinsky,Tropical Monte Carlo quadrature for Feynman integrals,Ann
M. Borinsky,Tropical Monte Carlo quadrature for Feynman integrals,Ann. Inst. H. Poincare D Comb. Phys. Interact.10(2023) 635–685, [2008.12310]
-
[62]
M. Borinsky, H. J. Munch and F. Tellander,Tropical Feynman integration in the Minkowski regime,Comput. Phys. Commun.292(2023) 108874, [2302.08955]. – 29 –
-
[63]
G. Heinrich, S. Jahn, S. P. Jones, M. Kerner, F. Langer, V. Magerya et al.,Expansion by regions with pySecDec,Comput. Phys. Commun.273(2022) 108267, [2108.10807]
-
[64]
Geometric approach to asymptotic expansion of Feynman integrals
A. Pak and A. Smirnov,Geometric approach to asymptotic expansion of Feynman integrals, Eur. Phys. J. C71(2011) 1626, [1011.4863]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[65]
N. Arkani-Hamed, A. Hillman and S. Mizera,Feynman polytopes and the tropical geometry of UV and IR divergences,Phys. Rev. D105(2022) 125013, [2202.12296]
-
[66]
Chebyshev Approximations of Feynman Integrals for Collider Physics
S. Abreu, A. Guerreiro and B. Page,Ancillary files for “Chebyshev Approximations of Feynman Integrals for Collider Physics”, 2026. 10.5281/zenodo.21134108
-
[67]
Byers and C
N. Byers and C. N. Yang,Physical Regions in Invariant Variables for n Particles and the Phase-Space Volume Element,Rev. Mod. Phys.36(1964) 595–609
1964
-
[68]
Analytic form of the two-loop planar five-gluon all-plus-helicity amplitude in QCD
T. Gehrmann, J. M. Henn and N. A. Lo Presti,Analytic form of the two-loop planar five-gluon all-plus-helicity amplitude in QCD,Phys. Rev. Lett.116(2016) 062001, [1511.05409]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[69]
Pentagon functions for massless planar scattering amplitudes
T. Gehrmann, J. M. Henn and N. A. Lo Presti,Pentagon functions for massless planar scattering amplitudes,JHEP10(2018) 103, [1807.09812]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[70]
Differential equations from unitarity cuts: nonplanar hexa-box integrals
S. Abreu, B. Page and M. Zeng,Differential equations from unitarity cuts: nonplanar hexa-box integrals,JHEP01(2019) 006, [1807.11522]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[71]
Analytic result for the nonplanar hexa-box integrals
D. Chicherin, T. Gehrmann, J. M. Henn, N. A. Lo Presti, V. Mitev and P. Wasser,Analytic result for the nonplanar hexa-box integrals,JHEP03(2019) 042, [1809.06240]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[72]
The two-loop five-point amplitude in $\mathcal{N} =4$ super-Yang-Mills theory
S. Abreu, L. J. Dixon, E. Herrmann, B. Page and M. Zeng,The two-loop five-point amplitude inN= 4super-Yang-Mills theory,Phys. Rev. Lett.122(2019) 121603, [1812.08941]
work page internal anchor Pith review Pith/arXiv arXiv 2019
-
[73]
D. Chicherin, T. Gehrmann, J. M. Henn, P. Wasser, Y. Zhang and S. Zoia,All Master Integrals for Three-Jet Production at Next-to-Next-to-Leading Order,Phys. Rev. Lett.123 (2019) 041603, [1812.11160]
- [74]
- [75]
-
[76]
Maheria,Semi- and Fully-Inclusive Phase-Space Integrals at Four Loops
V. Maheria,Semi- and Fully-Inclusive Phase-Space Integrals at Four Loops. PhD thesis, Hamburg U., 2022. [77]https://gitlab.com/pentagon-functions
2022
-
[77]
Kleiss, W
R. Kleiss, W. J. Stirling and S. D. Ellis,A New Monte Carlo Treatment of Multiparticle Phase Space at High-energies,Comput. Phys. Commun.40(1986) 359
1986
- [78]
- [79]
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.