the lattice for low moments, [@23], it is the purpose of this article to present the first one loop computations for the moments $n$ $=$ $2$ and $3$. Although lattice computations focus on the Landau gauge, partly as that gauge is less complicated to fix numerically, the Green’s function with the operator insertion is gauge dependent. So in developing another renormalization scheme to overcome one technical problem there are potential numerical errors in measurements from gauge fixing issues. Instead all our results will be in a general linear covariant gauge. Whilst the renormalization of these Wilson operators has already been carried out for the RI${}^\prime$ scheme, [@28; @29; @30], it may appear to be rather simple to follow that procedure for the latest scheme. This is not the case for an elementary reason. This is to do with the fact that the Wilson operators mix under renormalization, [@36]. It is widely accepted that the flavour singlet Wilson operators mix among themselves and the flavour non-singlet ones do not, [@1; @2; @3]. Indeed in the original context of [@1; @2; @3] this is the situation. However, that is only the case for the latter set if the Green’s function containing the 