New algorithms for Feynman integral reduction and $\varepsilon$-factorised differential equations
Iris Bree, Federico Gasparotto, Antonela Matijašić, Pouria Mazloumi, Dmytro Melnichenko, Sebastian Pögel, Toni Teschke, Xing Wang, Stefan Weinzierl, Konglong Wu, Xiaofeng Xu
TL;DR
The paper presents a constructive two-step framework for Feynman integral reduction that first selects master integrals via a geometry-informed Laporta ordering producing $F^{\bullet}$-compatible Laurent-$\varepsilon$ differential equations, and then converts these into an $\varepsilon$-factorised form through a systematic rotation. The first step relies on twisted cohomology and Baikov-represented integrands to define a robust, geometry-agnostic master basis; the second step removes $\varepsilon$-dependence using differential constraints, even when algebraic or transcendental functions appear. The approach is validated across polylogarithmic, elliptic, and Calabi–Yau-type geometries through a spectrum of examples, including one-, two-, three-, and four-loop integrals, demonstrating broad applicability and potential for further mathematical insights. The methodology promises improvements in analytical and numerical computations of perturbative QFT, and opens avenues for rigorous proofs, symmetry-aware extensions in twisted cohomology, and self-dual formulations of $\varepsilon$-factorised systems.
Abstract
In this paper, we give a detailed account of the algorithm outlined in [1] for Feynman integral reduction and $\varepsilon$-factorised differential equations. The algorithm consists of two steps. In the first step, we use a new geometric order relation in the integration-by-parts reduction to obtain a basis of master integrals, whose differential equations are of a Laurent polynomial form in the regularisation parameter $\varepsilon$ and compatible with a filtration. This step works entirely with rational functions. In a second step, we provide a method to $\varepsilon$-factorise the aforementioned Laurent differential equations. The second step may introduce algebraic and transcendental functions. We illustrate the versatility of the algorithm by applying it to different examples with a wide range of complexity.
