The Charm Quark Mass to Two-Loop Order

The truncation of the perturbative series at one loop order for the mass renormalization constants remain a significant systematic uncertainty in the determination of heavy quark masses in lattice QCD. We present here a high beta Monte Carlo calculation of the two loop mass renormalization constant for clover-improved fermions near the charm mass in the Fermilab heavy quark formalism. A preliminary value for the charm quark mass in the MS-bar scheme at two loop order is reported.