Iteratively Reduce Auxiliary Scalar Product in Multi-loop Integrals

In this paper, we construct a uniform formula that can iteratively reduce all auxiliary scalar product numerators of arbitrary multi-loop Feynman integrals. Integrals with such numerators commonly appear in Integration-By-Parts (IBP) relations. This formula is constructed with the generalized Sylvester’s determinant identity. Compared to that using only traditional IBP reduction method, the combination of the formula and the traditional IBP method shows a significant speed-up.