Skip to content

Bosonic Clouds Around Rotating Black Holes

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

A rotating horizon can amplify a bosonic wave only when the wave carries angular momentum with the appropriate sign. A long-lived cloud additionally needs confinement: a mass term or another effective barrier must keep the amplified field from escaping to infinity. A bound eigenmode in the superradiant range can grow; a stationary threshold cloud must also satisfy the horizon and infinity boundary conditions. Two inequalities alone do not construct that eigenmode.

Helpful background. Black-Hole Geometry and Mechanics supplies Kerr horizons; Black-Hole Perturbation Theory supplies separated modes and boundary conditions; Formulations of General Relativity supplies field evolution; Relativistic Matter and Compact Objects supplies bosonic matter models; Computational Relativity supplies nonlinear validation.

Horizon flux and the superradiant inequality

Section titled “Horizon flux and the superradiant inequality”

First take a real positive frequency ω\omega and an integer azimuthal number mm in e−iωt+imϕe^{-i\omega t+im\phi}. In four-dimensional asymptotically flat Kerr, the future horizon is generated by

χ=∂t+ΩH∂ϕ,ΩH=ar+2+a2\chi=\partial_t+\Omega_H\partial_\phi, \qquad \Omega_H=\frac{a}{r_+^2+a^2}

where a=J/Ma=J/M and r+=GM+(GM)2−a2r_+=GM+\sqrt{(GM)^2-a^2}, with physical mass MM, c=1c=1, and ∣a∣<GM|a|\lt GM. Acting on a scalar mode gives

χΦ=−i(ω−mΩH)Φ.\chi\Phi=-i(\omega-m\Omega_H)\Phi.

For a positive-frequency bosonic wave, the horizon energy flux associated with ∂t\partial_t has the schematic sign

E˙H∝ω(ω−mΩH)∣AH∣2.\dot E_H\propto \omega(\omega-m\Omega_H)|\mathcal A_H|^2.

Thus

0<ω<mΩH0\lt \omega\lt m\Omega_H

gives negative Killing-energy flux into the horizon and positive amplification outside. This sign and its scalar, vector, and tensor realizations are derived in Brito, Cardoso, and Pani 2020, §4.4, pp. 50–51, and §4.7.2, pp. 58–59, PDF. The normalization of AH\mathcal A_H varies, but the flux sign does not.

The field extracts angular momentum and energy in the ratio

dEdJ=ωm\frac{dE}{dJ}=\frac{\omega}{m}

for a monochromatic mode. Axisymmetric m=0m=0 modes cannot satisfy the positive-frequency superradiant inequality.

For a scalar of mass μ\mu, a quasibound mode has complex frequency ω=ωR+iωI\omega=\omega_R+i\omega_I, is ingoing at the future horizon, and decays at infinity. The asymptotic exponential is e−kre^{-kr} with k=μ2−ω2k=\sqrt{\mu^2-\omega^2} on the branch Re⁡k>0\operatorname{Re}k\gt0. In the weak-growth or weak-damping regime ∣ωI∣≪ωR|\omega_I|\ll\omega_R, the usual confinement condition is

ωR2<μ2,\omega_R^2\lt \mu^2,

with positive ωR\omega_R. This is the regime in which the monochromatic flux argument applies with ω\omega replaced by ωR\omega_R. The dimensionless gravitational coupling is

α=GMμ\alpha=GM\mu

with ℏ=c=1\hbar=c=1. For α≪1\alpha\ll1, the cloud is hydrogenic:

ωR≃μ[1−α22(n+ℓ+1)2],rcloud∼(n+ℓ+1)2GMα2.\omega_R\simeq\mu \left[ 1-\frac{\alpha^2}{2(n+\ell+1)^2} \right], \qquad r_{\rm cloud}\sim\frac{(n+\ell+1)^2GM}{\alpha^2}.

Here n=0,1,…n=0,1,\ldots counts radial overtones, ℓ\ell is the orbital angular number, and N=n+ℓ+1N=n+\ell+1 is the principal quantum number. The length scale follows by balancing nonrelativistic kinetic energy 1/(2μr2)1/(2\mu r^2) against gravitational binding GMμ/rGM\mu/r: the gravitational Bohr radius is 1/(GMμ2)=GM/α21/(GM\mu^2)=GM/\alpha^2. Excited hydrogenic states are larger by order N2N^2.

The imaginary part ωI\omega_I determines growth or decay. On the bound-state branch it crosses zero at synchronization. Scalar, Proca, and massive tensor theories have different spectra and fastest modes; the tensor case additionally requires its own consistent field equations. “A bosonic cloud” is not one universal rate.

The small-coupling scalar result originates in Detweiler 1980, pp. 2323–2326; the massive-field equations, hydrogenic limit, and unified instability picture are reviewed in Brito, Cardoso, and Pani 2020, §§5.6.2–5.6.5, pp. 105–115, PDF.

At

ωR=mΩH,\omega_R=m\Omega_H,

the horizon flux vanishes. A real-frequency bound solution satisfying regular horizon and decaying infinity conditions is a stationary cloud with respect to the horizon generator. The field itself can have harmonic phase dependence, while its stress tensor is stationary and axisymmetric for a complex field.

Clouds occur on discrete existence lines in Kerr parameter space because both boundary conditions must be satisfied. They are linear configurations: their stress does not yet backreact on the metric. Herdeiro and Radu showed that scalar clouds mark a bifurcation to fully nonlinear stationary black holes with synchronized complex scalar hair Herdeiro and Radu 2014, Eqs. (1)–(4), Figs. 1–2, and pp. 1–5, PDF. Existence of that branch does not prove it is dynamically selected or stable.

Real fields require care. A single real harmonic generally gives an explicitly time-dependent stress tensor, and nonlinear evolutions can radiate rather than settle to the same stationary complex-field branch. Vector clouds and Proca hair are distinct matter models.

  • Set a→0a\to0, so ΩH→0\Omega_H\to0: no positive-frequency mode is superradiant.
  • Set ω=mΩH\omega=m\Omega_H: the horizon flux vanishes.
  • Require both ωR<μ\omega_R\lt \mu and ωR<mΩH\omega_R\lt m\Omega_H for a growing massive bound state.
  • Track dimensions through α=GMμ\alpha=GM\mu; writing MμM\mu silently assumes G=1G=1.
  • Distinguish a linear cloud, a nonlinear stationary hairy solution, and a dynamical endpoint.

Environmental absorption, accretion, companions, plasma effects, or boson self-interactions can compete with the vacuum growth rate. Those effects do not alter the kinematic flux sign, but they can determine whether an astrophysical cloud becomes macroscopic.

A candidate scalar mode has m=1m=1, ωR=0.4/(GM)\omega_R=0.4/(GM), μ=0.45/(GM)\mu=0.45/(GM), and ΩH=0.45/(GM)\Omega_H=0.45/(GM). Do these values meet the weak-growth confinement and superradiance conditions? What changes at ΩH=0.4/(GM)\Omega_H=0.4/(GM)?

Solution

Because ωR<μ\omega_R\lt\mu, the mode is eligible to decay at infinity. Initially ωR<mΩH=0.45/(GM)\omega_R\lt m\Omega_H=0.45/(GM), so the leading monochromatic horizon energy flux is negative. A growing quasibound eigenmode is established only after solving the radial and angular boundary-value problem and finding ωI>0\omega_I\gt0. At ΩH=0.4/(GM)\Omega_H=0.4/(GM) the synchronization condition holds; an actual real-frequency bound eigenfunction then has zero horizon flux and is a threshold cloud.