Within the framework of field-theoretical description of second-order phase transitions via the three-dimensional
O(N)
vector model, accurate predictions for critical exponents can be obtained from
(resummation of) the perturbative series of renormalization-group functions, which are in
turn derived—following Parisi's approach—from the expansions of appropriate field
correlators evaluated at zero external momenta.
Such a technique was fully exploited
30
years ago in two seminal works of Baker, Nickel, Green and Meiron, which led to the knowledge of
the β-function up to the six-loop level; they succeeded in obtaining a precise numerical
evaluation of all needed Feynman amplitudes in momentum space by lowering the
dimensionalities of each integration with a cleverly arranged set of computational
simplifications. In fact, extending this computation is not straightforward, due both to the
factorial proliferation of relevant diagrams and the increasing dimensionality of their
associated integrals; in any case, this task can be reasonably carried on only in the
framework of an automated environment.
On the road towards the creation of such an environment, we here show how a strategy
closely inspired by that of Nickel and co-workers can be stated in algorithmic form, and
successfully implemented on a computer. As an application, we plot the minimized
distributions of residual integrations for the sets of diagrams needed to obtain RG functions
to the full seven-loop level; they represent a good evaluation of the computational effort
which will be required to improve the currently available estimates of critical exponents.