Perturbative Renormalization Factors of Bilinear Operators for Massive Wilson Quarks on the Lattice

Renormalization factors for local vector and axial vector currents for the Wilson quark action are perturbatively calculated to one loop order including finite quark masses from the ratio of the on-shell quark matrix elements in the Feynman gauge defined on the lattice and in the continuum. For large quark masses of order unity in lattice units, we find that finite quark mass effects are quite large: one-loop coefficients of the renormalization factors differ by 100% compared to those in the massless limit.