Analytic results for the two-loop renormalization constants of Feynman gauge QED

The two-loop renormalization constants $Z_m$ and $Z_1\text{\hspace{0.17em}⁣}=\text{\hspace{0.17em}⁣}{Z}_2$ of Feynman gauge quantum electrodynamics (QED) in the on-shell scheme are computed using dimensional regularization. Analytic results are obtained for these constants. In addition, analytic expressions are obtained for their expansion coefficients in the dimensional regularization parameter $\epsilon$ from order $1\mathrm{/}{\epsilon }^2$ through order ${\epsilon }^5$. The high-order expansion coefficients will be needed for future multiloop calculations in QED and will provide a useful check for multiloop calculations in QCD.