Two-loop renormalization of vector, axial-vector, and tensor fermion bilinears on the lattice

We compute the two-loop renormalization functions, in the RI' scheme, of local bilinear quark operators (psi) over bar Gamma psi, where Gamma corresponds to the vector, axial-vector, and tensor Dirac operators, in the lattice formulation of QCD. We consider both the flavor nonsinglet and singlet operators. We use the clover action for fermions and the Wilson action for gluons. Our results are given as a polynomial in c(SW), in terms of both the renormalized and bare coupling constant, in the renormalized Feynman gauge. Finally, we present our results in the (MS) over bar scheme, for easier comparison with calculations in the continuum. The corresponding results, for fermions in an arbitrary representation, together with some special features of superficially divergent integrals, are included in the appendices.