We solve the two-loop sunrise integral with unequal masses systematically to all orders in the dimensional regularisation parameter $\varepsilon$. In order to do so, we transform the system of differential equations for the master integrals to an $\varepsilon$-form. The sunrise integral with unequal masses depends on three kinematical variables. We perform a change of variables to standard coordinates on the moduli space ${ \mathcal M}_{ 1,3}$ of a genus one Riemann surface with three marked points. This gives us the solution as iterated integrals on $\overline{ \mathcal M}_{ 1,3}$. On the hypersurface $\tau=\mbox{ const}$ our result reduces to elliptic polylogarithms. In the equal mass case our result reduces to iterated integrals of modular forms.
Categories: hep-th hep-ph math-ph math.MP