IBP reduction coefficients made simple
Janko Boehm, Marcel Wittmann, Zihao Wu, Yingxuan Xu, Yang Zhang
TL;DR
The paper tackles the problem of oversized analytic IBP reduction coefficients in multi-loop Feynman integrals by introducing an improved Leinartas' multivariate partial fraction algorithm implemented in Singular. It leverages UT-basis properties to ensure denominators are symbol letters or polynomials in $D$, enabling substantial size reductions, often by orders of magnitude, when present. The approach is demonstrated across simple and frontier two-loop/topologies, including elliptic cases, with reported compression factors up to around 100x and consistent applicability even without a UT basis. A publicly available Singular library accompanies the work, highlighting practical impact for fast, storage-efficient, and scalable multi-loop amplitude computations. The findings suggest broad utility for analytical simplification and potential integration with finite-field and reconstruction techniques in high-precision collider physics.
Abstract
We present an efficient method to shorten the analytic integration-by-parts (IBP) reduction coefficients of multi-loop Feynman integrals. For our approach, we develop an improved version of Leinartas' multivariate partial fraction algorithm, and provide a modern implementation based on the computer algebra system Singular. Furthermore, We observe that for an integral basis with uniform transcendental (UT) weights, the denominators of IBP reduction coefficients with respect to the UT basis are either symbol letters or polynomials purely in the spacetime dimension $D$. With a UT basis, the partial fraction algorithm is more efficient both with respect to its performance and the size reduction. We show that in complicated examples with existence of a UT basis, the IBP reduction coefficients size can be reduced by a factor of as large as $\sim 100$. We observe that our algorithm also works well for settings without a UT basis.
