PV-reduction of sunset topology with auxiliary vector

The Passarino–Veltman (PV) reduction method has proven to be very useful for the computation of general one-loop integrals. However, not much progress has been made when it is applied to higher loops. Recently, we have improved the PV-reduction method by introducing an auxiliary vector. In this paper, we apply our new method to the simplest two-loop integrals, i.e., the sunset topology. We show how to use differential operators to establish algebraic recursion relations for reduction coefficients. Our algorithm can be easily applied to the reduction of integrals with arbitrary high-rank tensor structures. We demonstrate the efficiency of our algorithm by computing the reduction with the total tensor rank up to four.