Two-loop renormalization of three-quark operators in QCD
  Renormalization of composite three-quark operators in dimensional regularization is complicated by the mixing of physical and unphysical (evanescent) operators. This mixing must be taken into account in a consistent subtraction scheme. In this work we propose a particular scheme that allows one to avoid the necessity of additional finite renormalization and is convenient in QCD applications. As an illustration, we calculate the two-loop anomalous dimensions of local three-quark operators in this scheme. 
S. Kraenkl, A.N. Manashov
hep-ph