A Nonlinear Endpoint of Charged Horizon Instabilities

We numerically construct asymptotically extremal black holes through the nonlinear evolution of a charged scalar field. Our procedure – which extends the work of Murata-Reall-Tanahashi to include charged scalar dynamics – involves the fine-tuned scattering of wave packets within an initially super-extremal Reissner-Nordstrom spacetime. The resulting extremal solution develops an event horizon along which the energy density diverges and the charge density approaches a constant (i.e., the horizon forms with “hair”). We investigate this behavior from the perspective of critical phenomena in gravitational collapse, giving evidence that dynamical extremal black holes act as universal threshold solutions modulo this family-dependent hair. As in the linear instability of fixed extremal backgrounds, the scalar field decays outside the dynamical extremal horizon. But just inside the horizon, the scalar curvature appears to develop unbounded growth. This implies that near-threshold solutions without a black hole could develop correspondingly large curvatures visible from future null infinity.

💡 Research Summary

This paper presents a comprehensive numerical study of how a charged scalar field evolves non‑linearly on an initially super‑extremal Reissner‑Nordström (RN) background and, through fine‑tuned scattering, gives rise to a dynamical extremal black hole (BH) with “hair.” Building on the earlier work of Murata, Reall, and Tanahashi (MRT), which considered neutral scalar perturbations, the authors extend the analysis to include electromagnetic coupling, thereby probing the enhanced Aretakis‑type instability expected for charged matter.

The authors work in spherical symmetry using double‑null coordinates (U outgoing, V ingoing) with the metric ds² = –2 f dU dV + r² dΩ². The complex scalar field ϕ is decomposed into real and imaginary parts (ξ, Π), and the gauge‑invariant modulus P = |rϕ| is introduced. Maxwell’s equations determine the gauge potential Aμ and the enclosed charge Q, while Einstein’s equations provide evolution equations for the areal radius r(U,V) and the metric function f(U,V). The full coupled Einstein‑Maxwell‑Klein‑Gordon system is solved with a finite‑difference scheme that incorporates adaptive mesh refinement based on the inverse of the UV metric component, ensuring high resolution near any apparent horizon that may form.

Initial data are constructed by pasting a compact, ingoing charged pulse onto a null cone emanating from a chosen point (U₀, V₀). The past of this cone is taken to be exact RN with mass M₀ = 1 and charge Q₀, while the pulse is specified by an amplitude A₀ = 0.01, a width of 20, and a local mass‑to‑charge ratio ˜ω = 1. The coupling constant is set to e Q₀ = 0.6, a regime known to maximize the instability. By varying Q₀ while keeping the pulse profile fixed, the authors perform a bisection search for a critical charge Q* at which the trapped region shrinks to zero size and the inner and outer horizons merge at future null infinity, signalling the formation of a dynamical extremal BH.

The numerical experiments reveal several key phenomena:

1. Critical Scaling – The coordinate V_trap, marking the earliest appearance of a trapped surface, diverges as V_trap ∝ (Q* – Q₀)^{-1/2} when Q₀ approaches Q* from below. This 1/2 exponent matches the prediction of MRT and analytic arguments based on the near‑horizon AdS₂ symmetry of extremal RN.

2. Hair Formation – At the event horizon of the resulting extremal BH, the energy density of the scalar field diverges while the charge density settles to a constant value. This “hair” is a direct consequence of the charged scalar’s ability to deposit charge on the horizon, a feature absent in the neutral case.

3. Aretakis‑Like Growth – Outside the horizon, the scalar field exhibits the familiar non‑decay and gradient blow‑up characteristic of the Aretakis instability. Inside the horizon, however, the Ricci scalar appears to grow without bound, a behavior the authors term “blueshift focusing instability.” This suggests that, for near‑critical solutions on the dispersive side, arbitrarily large curvatures could become visible at future null infinity.

4. Partial Violation of Universality – While the critical solution appears to be “one‑mode unstable” and shares universal scaling laws, the magnitude of the horizon charge density (the hair) and the coefficients governing gradient growth depend on the specific one‑parameter family of initial data used. Thus, the extremal critical solution possesses a family‑dependent hair, breaking the strict universality familiar from scalar collapse.

The paper also discusses the fate of sub‑critical (BH‑forming) versus super‑critical (dispersive) evolutions. In the sub‑critical regime, back‑reaction reduces the charge‑to‑mass ratio, driving the spacetime toward a sub‑extremal BH. In the super‑critical regime, no trapped surfaces form within the simulated domain, yet the interior curvature still shows signs of unbounded growth, hinting at a possible naked singularity formation if the evolution could be continued further.

In conclusion, the authors demonstrate that charged scalar fields can non‑linearly generate extremal black holes with non‑trivial hair, that the associated Aretakis‑type instability persists in the fully back‑reacted setting, and that the critical collapse picture extends to include a family‑dependent hair parameter. The work opens several avenues for future research: extending the analysis to rotating (Kerr) backgrounds, exploring non‑spherical perturbations, and improving numerical stability to follow the interior evolution deeper into the potential singular regime.