Chebyshev approximations of Feynman integrals for collider physics
5 affiliations
- 1Theoretical Physics Department, CERN
- 2The University of Edinburgh
- 3Universidade de Lisboa
- 4LIP — Laboratório de Instrumentação e Física Experimental de Partículas
- 5Ghent University
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.
Related
BibTeX
@article{abreu2026chebyshev,
title = {{Chebyshev approximations of Feynman integrals for collider physics}},
author = {Samuel Abreu and Afonso Guerreiro and Ben Page},
journal = {Journal of High Energy Physics},
year = {2026},
month = nov,
doi = {10.1007/jhep11(2026)038},
url = {https://link.springer.com/article/10.1007/JHEP11(2026)038},
note = {CC BY 4.0},
}
RIS
TY - JOUR
TI - Chebyshev approximations of Feynman integrals for collider physics
AU - Samuel Abreu
AU - Afonso Guerreiro
AU - Ben Page
PY - 2026
DA - 2026/11/05
JO - Journal of High Energy Physics
DO - 10.1007/jhep11(2026)038
UR - https://link.springer.com/article/10.1007/JHEP11(2026)038
AB - 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.
SN - 1029-8479
ER -