Hard-Region Fermion Self-Energy and Fermion--Photon Vertex in Thermal QED through Two Loops

In massless thermal QED in a general covariant $R_{\xi }$ gauge, we compute hard-region contributions to the fermion self-energy and the off-shell fermion–photon vertex at one-loop next-to-leading power and two-loop leading power in the soft-momentum expansion. The zero-temperature counterterms are also included to renormalize these hard amplitudes. Chiral invariance restricts the off-shell two-fermion–$N$-photon vertex to vector (${\gamma }_{\mu }$) and axial-vector (${\gamma }_{\mu }{\gamma }_5$) Dirac structures. The vector part is constrained by the Ward–Takahashi identity (WTI), while the axial part is transverse to the photon momentum. The symmetries and hermiticity of the theory impose definite constraints – including reality, momentum reversal, and fermion-leg exchange properties – which lead to selection rules in the soft expansion: a contribution at power $r$ can be nonzero only when $r+N+1$ is even, with $N=0$ for the self-energy. At the integrand level, we decompose the amplitudes into independent statistical and gauge sectors, and verify the WTI and axial transversality sector by sector. Notably, the gauge-dependent sectors split into mixed metric–longitudinal and fully longitudinal sectors at two loops. The latter vanishes at leading power and reappears at next-to-leading power while satisfying the constraints. We show that these hard-region self-energy corrections do not generate a finite contribution to the fermion damping rate. These results provide hard-region input for mass shifts and inclusive rates, as well as for the construction of the next-to-leading-order fermionic effective Lagrangian.