An algorithm for the high-energy expansion of multi-loop diagrams to next-to-leading logarithmic accuracy

We present an algorithm to compute arbitrary multi-loop massive Feynman diagrams in the region where the typical energy scale Q is much larger than the typical mass scale M, i.e., Q≫M, while various different energy and mass parameters may be present. In this region we perform an asymptotic expansion and, using sector decomposition, we extract the leading contributions resulting from ultraviolet and mass singularities, which consist of large logarithms ln(Q2/M2) and 1/ɛ poles in D=4−2ɛ dimensions. To next-to-leading accuracy, at L loops all terms of the form αLɛ−klnj(Q2/M2) with j+k=2L and j+k=2L−1 are taken into account. This algorithm permits, in particular, to compute higher-order next-to-leading logarithmic electroweak corrections for processes involving various kinematical invariants of the order of hundreds of GeV and masses MW∼MZ∼MH∼mt of the order of the electroweak scale, in the approximation where the masses of the light fermions are neglected.