quantum field theory - Two-loop regularization

Working out some quantum field theory computations, I have to find out the value of the two-loop Feynman integral

$$I(p)=\int\frac{d^4p_1}{(2\pi)^4}\frac{d^4p_2}{(2\pi)^4}\frac{1}{(p_1^2+m_1^2)(p_2+m_2^2)[(p-p_1-p_2)^2+m_3^2]}.$$ This integral is rather common and so, its value at small $p$ should be already well-knwon. But I was not able to find out it in literature. I would also appreaciate to see all the procedure to get the right value with whatever regularization procedure one likes.

Thanks beforehand.

Answer

In the case of equal masses, there is an analytical solution (of this diagram known by the name "the two loop sunrise diagram" for the obvious reason) in terms of hypergeometric functions given by O.V. Tarasov (equation 4.32). There is also a numerical method given by: Pozzorini and Remiddi.

In the case of unequal masses Müller-Stach, Weinzierl and Zayadeh were able to obtain a second order differential equation for this graph.