Hypergeometric expression for a three-loop vacuum integral

Using the corresponding Mellin–Barnes representation, we derive holonomic hypergeometric system of linear partial differential equations (PDEs) satisfied by Feynman integral of a three-loop vacuum with five propagators. Through the multidimensional residue theorem in dimensional regularization, the scalar integral can be written as the summation of multiple hypergeometric functions, whose convergent regions can be obtained by the Horn’s convergent theory. The numerical continuation of the scalar integral from convergent regions to whole kinematic regions can be accomplished with the finite element methods, when the system of PDEs can be treated as the stationary conditions of a functional under the restrictions.