Quasibound states of charged Dirac fields around Einstein-Euler-Heisenberg black holes

Abstract

We study quasibound states of a massive charged Dirac field on the Einstein--Euler--Heisenberg (EEH) black hole, the extension of Reissner--Nordström (RN) by the leading weak-field vacuum-polarization term of QED. We derive the background to first order in the Euler--Heisenberg coupling $a$. We show that it admits a nondegenerate horizon at $Q=M$ and an extremal charge $Q_{\rm ext}=M[1+a/(40M^2)+\mathcal{O}(a^2)]>M$, with horizon data that depend on $\sqrt a$ near extremality. From the separated radial equations, we derive an identity that excludes real-frequency quasibound states on any background of this class, so damped branches cannot become unstable. Since the EEH corrections first enter at $\mathcal{O}(r^{-5})$ in the electrostatic potential, the trapping condition $Mμ^2-qQω_R>0$ and the hydrogenic spectrum coincide with RN, and no fine-structure-order splitting arises. Two-sided shooting with a Muller root finder shows that the real frequencies shift only by near-zone amounts that are linear in $a$ and strongly suppressed with angular momentum. The damping rates respond in a branch-dependent way. At $Q/M=0.99$ the $\ell=0$ ground state becomes slightly longer-lived, while the $(n_r,j,\ell,λ)=(0,\tfrac32,1,-2)$ state decays up to $80\%$ faster for $a/M^2\le1$, and near the trapping threshold the damping rate increases by up to $60\%$. Switching on the metric and potential corrections separately, together with a WKB analysis of the inner barrier, identifies the mechanisms: Coulomb screening always increases leakage, while the metric correction controls the sign and size of the $\ell$-dependent effect. Characteristic time-domain evolutions with Richardson-extrapolated matrix-pencil extraction reproduce the frequency-domain shifts.