NNLO anomalous dimension matrix for twist-two flavor-singlet operators

Conformal symmetry of QCD is restored at the Wilson-Fisher critical point in noninteger $4-2\epsilon$ space-time dimensions. Correlation functions of multiplicatively renormalizable operators with different anomalous dimensions at the critical point vanish identically. We show that this property allows one to calculate off-diagonal parts of the anomalous dimension matrices for leading-twist operators from a set of two-point correlation functions of gauge-invariant operators which can be evaluated using standard computer algebra techniques. As an illustration, we present the results for the NNLO anomalous dimension matrix for flavor-singlet QCD operators for spin $N\le 8$