Massive Fields and Oscillatory Tails

A field mass moves the relevant branch points from zero frequency to the thresholds $\omega=\pm\mathsf m$. Their interference produces oscillations at approximately the Compton frequency. For a massive scalar on a stationary, spherically symmetric, asymptotically flat subextremal black hole, a fixed spherical harmonic generically has a very-late fixed-radius envelope $(t^*)^{-5/6}$ multiplying a nontrivial oscillatory profile. This precise theorem is narrower than the frequently quoted claim of a universal $t^{-5/6}\sin(\mathsf m t)$ tail.

Required background. PDEs and Green Functions for Gravitation supplies threshold resolvents, stationary phase, and retarded evolution.

Helpful background. Black-Hole Geometry and Mechanics supplies subextremal spherical exteriors; Formulations of General Relativity distinguishes field components; Causality and Global Structure supplies horizon-regular time and observation limits; Black-Hole Perturbation Theory supplies mode and Green-function methods.

Massive dispersion creates two thresholds

On Schwarzschild, write $\Phi=r^{-1}\psi_{\ell m}(t,r)Y_{\ell m}$ for

$$(\Box_g-\mathsf m^2)\Phi=0, \qquad \mathsf m>0.$$

The radial equation is

$$\left[-\partial_t^2+\partial_{r_*}^2 -f\left(\mathsf m^2+\frac{\ell(\ell+1)}{r^2} +\frac{2\mu}{r^3}\right)\right]\psi_{\ell m}=0,$$

with $\mu=GM$ and $f=1-2\mu/r$. At infinity the wavenumber is $k(\omega)=\sqrt{\omega^2-\mathsf m^2}$. It vanishes at $\omega=\pm\mathsf m$, so the outgoing solution and resolvent are multivalued there. Choose the retarded sheet by continuation from $\operatorname{Im}\omega>0$ and a cut joining the thresholds, or equivalent cuts extending from each threshold; the final retarded field is independent of that drawing when the rim orientations are transformed consistently.

Frequencies $\lvert\omega\rvert>\mathsf m$ propagate to infinity, whereas $\lvert\omega\rvert<\mathsf m$ are evanescent there. Attractive curvature also permits quasi-bound poles on the continued resolvent. Those poles generate long-lived exponentially damped oscillations. They are discrete contributions and must not be mistaken for the continuous threshold cut that gives the ultimate inverse-power tail.

Intermediate and asymptotic regimes differ

For a light scalar on Schwarzschild, neglecting strong Coulomb backscatter near the thresholds gives the familiar intermediate fixed-radius behavior

$$\psi_{\ell m}(t,r) \sim C_{\ell m}(r)t^{-(\ell+3/2)} \sin\!\left(\mathsf m t+\delta_\ell\right).$$

This is a regime statement, often realized for $\mu\ll t\ll \mathsf m^{-3}\mu^{-2}$ when $\mathsf m\mu\ll1$, not the final asymptotic theorem. At later time the Coulomb-like $-2\mu\mathsf m^2/r$ part of the far potential makes the threshold phase rapidly varying. A stationary point approaches $\omega=\mathsf m$ with $\mathsf m-\omega_c=O(t^{-2/3})$. Its width and amplitude combine to give the slower $t^{-5/6}$ envelope, independent of fixed $\ell$ at leading power.

The intermediate and asymptotic Schwarzschild regimes, including the role of the long-range curvature potential near threshold, are developed in Koyama and Tomimatsu 2001, §§ II–IV.

The transition can occur extremely late, and quasi-bound ringing can obscure both power laws. A time-domain fit should therefore track the local envelope, instantaneous phase, and several radii; seeing oscillations near $\mathsf m$ does not by itself identify the asymptotic regime.

Precise fixed-harmonic theorem

The sharp result available in the assigned source assumes a stationary, spherically symmetric metric

$$g=-F(r)dt^2+F(r)^{-1}dr^2+r^2d\Omega^2, \qquad F(r)=1-\frac{2\mu}{r}+\frac{D(r)}{r^2},$$

with a nondegenerate outer horizon, positive surface gravity, controlled asymptotic derivatives of $D$, and a positivity condition excluding a zero-frequency resonance. It includes Schwarzschild and every subextremal Reissner–Nordström black hole. Fix one spherical harmonic $L$ and smooth, compactly supported initial data on a spacelike $t^*=0$ slice that may meet the future horizon. For $r$ in any fixed compact interval from the horizon to $R_0$, the theorem gives, for every $\delta>0$,

$$\begin{aligned} r\Phi_L(t^*,r) ={}&\mathcal A_L[\text{data}]\,\lvert u_H^*(r)\rvert(t^*)^{-5/6}\\ &\times\sum_{q=1}^{\infty}(-\gamma_{\mathsf m})^{q-1}q^{1/3} \cos\!\left(\Psi_q(t^*)+\varphi_H(r)+\theta_L\right) +O_{R_0,L}\!\left((t^*)^{-1+\delta}\right), \end{aligned}$$

where $\lvert\gamma_{\mathsf m}\rvert<1$ and

$$\Psi_q(t^*) =\mathsf m t^* -\frac32(2\pi\mu)^{2/3}\mathsf m q^{2/3}(t^*)^{1/3} +\varphi_- -\frac\pi4+O\!\left((t^*)^{-1/3}\right).$$

Thus the profile contains infinitely many sublinear phase corrections; it is not exactly a single sine. At fixed $r>r_+$, $t^*=t+O(1)$ asymptotically. On the horizon, $t^*$ is an advanced-time coordinate while Schwarzschild $t$ is singular. The amplitude is generic but can vanish on a codimension-one hyperplane of data, in which case the remainder gives faster decay.

These hypotheses and the profile are stated and proved in Pasqualotto, Shlapentokh-Rothman, and Van de Moortel 2024, Theorem I, Remarks 1.1–1.5, and Theorem 3.1, pp. 2–6 and 9–11. As of the evidence cutoff 2026-08-12, arXiv:2303.17767v2, dated 5 January 2024, remains a preprint with no journal reference on its arXiv record. That publication status is separate from the theorem’s explicit mathematical content.

Scope of the five-sixths law

The theorem is for a massive scalar, one fixed harmonic—or a fixed finite sum under its stated version—on a stationary spherical subextremal background, with localized regular data and fixed compact $r$. It does not establish the same profile uniformly at null infinity, for an unrestricted infinite angular sum, on Kerr, for extremal horizons, for massive vector or tensor systems, or for a nonlinear self-gravitating massive field. Earlier numerical and asymptotic evidence in some of those settings is valuable but not the same theorem.

Full nonlinear decay theorems are coordinated with Mathematical Relativity; observational late-time fitting requires a source and detector model. Here, the independent checks are the two threshold locations, the $t^{-5/6}$ stationary phase scaling, the horizon-regular coordinate, and separation of cut and quasi-bound-pole contributions.

Exercise: recover the one-third phase scale

Near $\omega=\mathsf m$, suppose the cut phase has the schematic form

$$S(\delta)=-\mathsf m t+\delta t+A\delta^{-1/2}, \qquad \delta=\mathsf m-\omega>0,\quad A>0.$$

Find the stationary scaling of $\delta_c$ and the size of the sublinear phase correction. Explain why a pure $\sin(\mathsf m t)$ model misses information.

Solution

The stationary equation is

$$S'(\delta)=t-\frac{A}{2}\delta^{-3/2}=0,$$

so $\delta_c=(A/2t)^{2/3}=O(t^{-2/3})$. Both remaining phase terms scale as

$$\delta_ct=O(t^{1/3}), \qquad A\delta_c^{-1/2}=O(t^{1/3}).$$

Hence the oscillation is centered on $\mathsf m t$ but carries a growing sublinear $t^{1/3}$ correction. The rigorous profile also sums harmonics of this scattering phase indexed by $q$. A single sine can reproduce the carrier frequency and envelope while accumulating an order-$t^{1/3}$ phase error.

References

• Koyama, Hiroko, and Akira Tomimatsu. “Asymptotic Tails of Massive Scalar Fields in Schwarzschild Background.” Physical Review D 64 (2001): 044014. DOI. Open PDF
• Pasqualotto, Federico, Yakov Shlapentokh-Rothman, and Maxime Van de Moortel. “The Asymptotics of Massive Fields on Stationary Spherically Symmetric Black Holes for All Angular Momenta.” arXiv preprint arXiv:2303.17767 (2024), version 2, 5 January 2024. arXiv.