New algorithms for Feynman integral reduction and epsilon-factorized differential equations

In this paper, we give a detailed account of the algorithm outlined in [I. Bree ( ϵ Collaboration), companion Letter, .] for Feynman integral reduction and ϵ -factorized 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 on the maximal cut are of a Laurent polynomial form in the regularization parameter ϵ and compatible with a filtration. This step works entirely with rational functions. In a second step, we provide a method to ϵ -factorize 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.