Threshold asymptotics and decay for massive Maxwell on subextremal Reissner--Nordström

We study the neutral massive Maxwell (Proca) equation on subextremal Reissner–Nordström exteriors. After spherical-harmonic decomposition, the odd sector is scalar, while the even sector remains a genuinely coupled $2\times 2$ system. Our starting point is that this even system admits an exact asymptotic polarization splitting at spatial infinity. The three resulting channels carry effective angular momenta $\mathrm{\ell }-1$, $\mathrm{\ell }$, and $\mathrm{\ell }+1$, and these are precisely the indices that govern the late-time thresholds. % For each fixed angular momentum we develop a threshold spectral theory for the cut-off resolvent. We prove meromorphic continuation across the massive branch cut, rule out upper-half-plane modes and threshold resonances, and obtain explicit small- and large-Coulomb expansions for the branch-cut jump. Inverting this jump yields polarization-resolved intermediate tails together with the universal very-late $t^{-5\mathrm{/}6}$ branch-cut law. % At the full-field level, high-order angular regularity allows us to sum the modewise leading terms on compact radial sets and obtain a two-regime asymptotic expansion for the radiative branch-cut component of the Proca field, with explicit coefficient fields and quantitative remainders. We also analyze the quasibound resonance branches created by stable timelike trapping, prove residue and reconstruction bounds, and derive a fully self-contained dyadic packet estimate. As a result, the unsplit full Proca field obeys logarithmic compact-region decay, while the radiative branch-cut contribution retains explicit polynomial asymptotics and explicit leading coefficients.