Feynman parameter integration through differential equations
Martijn Hidding, Johann Usovitsch
TL;DR
This paper introduces an automated numerical framework for multi-loop Feynman integrals that iteratively applies Feynman's trick to merge propagators, generating simplified topologies whose master integrals are solved via differential equations in a Feynman parameter. Solutions are expressed as piecewise generalized power series, which are integrated term-by-term to obtain the original integral, avoiding manual boundary-condition specification. Regularization is handled through quasi-finite bases or epsilon-sampling with pole-resolving techniques, and threshold singularities are managed with carefully chosen iδ prescriptions. The approach demonstrates substantial reductions in the number of master integrals and computational effort, showcased by a non-planar double pentagon example and validated against AMFlow, pointing toward an efficient, automatic route for complex Feynman-integral computations.
Abstract
We present a new method for numerically computing generic multi-loop Feynman integrals. The method relies on an iterative application of Feynman's trick for combining two propagators. Each application of Feynman's trick introduces a simplified Feynman integral topology which depends on a Feynman parameter that should be integrated over. For each integral family, we set up a system of differential equations which we solve in terms of a piecewise collection of generalized series expansions in the Feynman parameter. These generalized series expansions can be efficiently integrated term by term, and segment by segment. This approach leads to a fully algorithmic method for computing Feynman integrals from differential equations, which does not require the manual determination of boundary conditions. Furthermore, the most complicated topology that appears in the method often has less master integrals than the original one. We illustrate the strength of our method with a five-point two-loop integral family.
