Bilinear quark operators in the RI/SMOM scheme at three loops

We consider the renormalization of the matrix elements of the bilinear quark operators ψ bar ψ, ψ bar γ_μ ψ, and ψ bar σ_μν ψ at next-to-next-to-next-to-leading order in QCD perturbation theory at the symmetric subtraction point. This allows us to obtain conversion factors between the MS ‾ scheme and the regularization invariant symmetric momentum subtraction (RI/SMOM) scheme. The obtained results can be used to reduce the errors in determinations of quark masses from lattice QCD simulations. The results are given in Landau gauge.