On Epsilon Factorized Differential Equations for Elliptic Feynman Integrals

In this paper we develop and demonstrate a method to obtain epsilon factorized differential equations for elliptic Feynman integrals. This method works by choosing an integral basis with the property that the period matrix obtained by integrating the basis over a complete set of integration cycles is diagonal. The method is a generalization of a similar method known to work for polylogarithmic Feynman integrals. We demonstrate the method explicitly for a number of Feynman integral families with an elliptic highest sector.