Determining the closed forms of the $O(a_s^3)$ anomalous dimensions and Wilson coefficients from Mellin moments by means of computer algebra

TL;DR

This paper presents a computer algebra method to derive closed-form expressions for 3-loop anomalous dimensions and Wilson coefficients in QCD by solving high-order recurrence relations from Mellin moments.

Contribution

It introduces a novel approach to determine closed forms of complex quantum field theory quantities using recurrence relations and harmonic sums, demonstrated on 3-loop QCD calculations.

Findings

Successfully derived recurrences of high order and degree for 3-loop quantities.

Solved large systems requiring extensive computational resources within feasible time.

Showed the method's potential for high-precision calculations in theoretical physics.

Abstract

Single scale quantities, as anomalous dimensions and hard scattering cross sections, in renormalizable Quantum Field Theories are found to obey difference equations of finite order in Mellin space. It is often easier to calculate fixed moments for these quantities compared to a direct attempt to derive them in terms of harmonic sums and their generalizations involving the Mellin parameter $N$. Starting from a sufficiently large number of given moments, we establish linear recurrence relations of lowest possible order with polynomial coefficients of usually high degree. Then these recurrence equations are solved in terms of d'Alembertian solutions where the involved nested sums are represented in optimal nested depth. Given this representation, it is then an easy task to express the result in terms of harmonic sums. In this process we compactify the result such that no algebraic relations occur among the sums involved. We demonstrate the method for the QCD unpolarized anomalous dimensions and massless Wilson coefficients to 3--loop order treating the contributions for individual color coefficients. For the most complicated subproblem 5114 moments were needed in order to produce a recurrence of order 35 whose coefficients have degrees up to 938. About four months of CPU time were needed to establish and solve the recurrences for the anomalous dimensions and Wilson coefficients on a 2 GHz machine requiring less than 10 GB of memory. No algorithm is known yet to provide such a high number of moments for 3--loop quantities. Yet the method presented shows that it is possible to establish and solve recurrences of rather large order and and degree, occurring in physics problems, uniquely, fast and reliably with computer algebra.