Ultralight Bosons in Astrophysics

• Ultralight bosons are a class of hypothetical, extremely low-mass particles (10⁻²⁰ to 10⁻¹⁰ eV) that form large-scale coherent fields around astrophysical black holes.
• Their superradiant instabilities enable the extraction of black hole spin energy, leading to the formation of boson clouds that emit nearly monochromatic gravitational waves.
• Astrophysical observations, including black hole spin distributions and binary dynamics, constrain boson mass ranges and illuminate their role in dark matter and fundamental physics.

Ultralight bosons are a class of hypothetical, extremely low-mass bosonic particles, including both scalars and vectors, that arise naturally in extensions of the Standard Model, such as the string axiverse. Their masses span a wide range, typically from $10^{-20}$ eV up to $10^{-10}$ eV, with substantial attention given to bosons in the $10^{-19}$ eV to $10^{-11}$ eV range due to their astrophysical and cosmological relevance. Ultralight bosons form large-scale coherent fields with Compton wavelengths comparable to or exceeding the size of astrophysical black holes; as a result, they can give rise to novel strong-gravity phenomena, notably via superradiant instabilities. The amplification of these fields around rapidly spinning black holes can lead to the formation of macroscopic bosonic condensates—typically termed “boson clouds”—whose evolution and gravitational-wave signatures are actively targeted by ground- and space-based gravitational-wave observatories.

1. Superradiant Instability and Boson Cloud Formation

Ultralight bosons exploited in this context can be scalars, such as axions, or vectors (“dark photons”); both satisfy quantized field equations in curved spacetime and are characterized by their mass $m_s$ or $m_V$. Around a Kerr black hole of mass $M$ and spin parameter $a$, the presence of a bosonic field leads to the possibility of extracting the black hole's rotational energy—provided specific conditions are met. The superradiant instability criterion is

$0 < \omega_R < m \Omega_H,$

where $\omega_R$ is the real part of the boson mode frequency, $m$ is the azimuthal number, and $\Omega_H$ is the black hole’s horizon angular velocity. When this condition is satisfied, massive boson modes with Compton wavelength $\lambda \sim \hbar/(m_s c)$ comparable to the gravitational radius of the black hole ($2GM/c^2$) become bound and can exponentially extract angular momentum (and some mass) from the black hole. The instability growth time $\tau_{\text{inst}} = 1/\omega_I$ (with $\omega_I>0$ the imaginary part of the frequency) can be orders of magnitude longer than the black hole’s dynamical time, but for vector bosons in particular, the growth rates can be extremely fast.

The boson cloud or “gravitational atom” is described by a field expansion analogous to hydrogenic orbitals: $$\Psi = \operatorname{Re}\left\{\sum_n \sum_\ell \sum_m e^{-i \omega t + i m \phi} S(\theta) R(r) \right\},$$ with the dominant instability typically in the $\ell = m = 1$ (scalar) or $j = m = 1$ (vector) mode, and the frequencies following a perturbed hydrogenic spectrum, e.g.

$$\omega_{nlm} \simeq m_s - \frac{m_s}{2} \left( \frac{M m_s}{\ell+n+1} \right)^2 + i \, \gamma_\ell (M m_s)^{4\ell+5} \ldots$$

The cloud grows until the superradiant condition is saturated, extracting a significant fraction of the black hole's original spin and up to $\sim$10% of its mass.

2. Gravitational-Wave Signatures from Boson Clouds

The boson cloud’s time-dependent, nonspherical mass distribution acts as a long-lived quadrupole, radiating nearly monochromatic gravitational waves (GWs). The emission frequency is typically set by

$$f_{\text{gw}} \simeq \frac{2\omega_R}{2\pi} \simeq \frac{m_s}{\pi},$$

and the GW strain amplitude (averaged over angles) may be expressed as

$$h \simeq \sqrt{\frac{2}{5\pi}} \, \left( \frac{M}{M_\odot} \right) \left( \frac{M_S}{M} \right) A(\chi, M m_s),$$

where $M_S$ is the cloud mass, $\chi$ is the final dimensionless spin, and $A$ is a numerically computed dimensionless function. The signal is quasi-monochromatic, with small secular spin-up due to the loss of cloud mass and gravitational radiation.

Detection prospects depend on the boson mass and black-hole parameters:

• Advanced LIGO/Virgo: stellar-mass black holes, $m_s \in [2\times10^{-13},10^{-12}]\,\text{eV}$, are the prime targets. In optimistic astrophysical scenarios, LIGO could detect up to $\sim10^4$ resolvable events over 4 years ($m_s \sim 3 \times 10^{-13}$ eV), or observe a stochastic GW background in a similar range.
• LISA: probes massive black holes, $$m_s \in [5\times10^{-19}, 5\times10^{-16}]\,\text{eV}$$, with up to $\sim10^3$ resolvable events over 4 years near $m_s \sim 10^{-17}\,\text{eV}$.

Detection strategies for continuous signals employ both fully coherent and semicoherent search methodologies to optimally extract long-duration, nearly monochromatic signals given their small frequency drift.

3. Imprints on Black Hole Spin Distributions and the Regge Plane

Superradiant extraction of angular momentum by ultralight bosons leads to population-level signatures in black hole mass-spin distributions (“Regge plane”). For a given mass $M$ and boson mass $m_s$, superradiance limits the attainable final spin $J_f$, leaving a “hole” or depletion region at high spins above a calculable threshold. This critical spin is given by

$$J_f = \frac{4 m M_f^3 \omega_R}{m^2 + 4 M_f^2 \omega_R^2},$$

with the associated mass loss given by

$M_f = M_i - (\omega_R/m)(J_i - J_f).$

Population studies utilize hierarchical Bayesian inference over catalogs of detected black hole mergers (e.g., with parameters from GWTC-2) to statistically combine weak spin measurements across many events. By comparing the observed distribution to theoretical exclusion regions as a function of $m_s$, one can:

• Exclude (at high confidence) scalar bosons in the mass interval $[1.3 \times 10^{-13}, 2.7 \times 10^{-13}]$ eV if rapidly spinning black holes are observed in this range (Ng et al., 2020).
• Rule out or confirm the existence of non-interacting bosons in the range $[10^{-13}, 3 \times 10^{-12}]$ eV given hundreds of high-SNR detections (Ng et al., 2019).

The most stringent constraints are driven by individual high-spin outliers and depend on the assumed timescale for superradiant extraction (long inspiral timescales yield stronger constraints; shorter delays reduce their severity).

4. Binary Dynamics and Resonant Phenomena

When a black hole hosting a boson cloud is in a binary system, the presence of the companion modifies the cloud's evolution via tidal perturbations that induce resonant transitions (Landau–Zener type) between growing and decaying modes of the cloud. Key features include:

• Resonant transitions occur when the orbital frequency matches the energy splitting (hyperfine or Bohr transitions) between cloud energy levels, with resonance conditions such as

$$\Omega \simeq \epsilon_h \quad \text{(hyperfine)}, \qquad f_{\text{res}}^{(h)} \simeq \epsilon_h/\pi,$$

where $\epsilon_h$ is the small energy difference between states.

• The companion-induced perturbation can cause rapid depletion ("extinction") of the cloud, sharply reducing or terminating the associated continuous GW signal prior to merger.
• These resonant processes also impart sharp, time-dependent corrections to finite-size parameters of the black hole (e.g., spin-induced quadrupole, $\kappa$, and tidal Love number, $\Lambda$), yielding observable phase features in the binary inspiral waveform.

Enhanced mass quadrupoles carried by the cloud in “gravitational atom” binaries lead to large retrograde apsidal precession in eccentric systems, inducing distinctive modulations either in the GW signal—seen as triplet splittings of each orbital harmonic—or in pulsar-timing parameters if the companion is a pulsar (Su et al., 2021).

5. Multiband and Population-Based Detection Methodologies

Multiband detection strategies leverage the strengths of both space-based (LISA) and ground-based (Cosmic Explorer, Einstein Telescope) interferometers:

• LISA first identifies and precisely characterizes a BBH inspiral, tightly constraining the masses and spins of the component black holes.
• Ground-based instruments then perform a targeted search for a continuous GW signal from a post-merger boson cloud, with cloud properties predicted from the observed pre-merger binary (Ng et al., 2020).
• The joint observation allows constraints (or possible detections) of boson masses in the band $25 \lesssim \mu/(10^{-15}\,\text{eV}) \lesssim 500$.

Continuous-wave directed searches exploit robust, computationally tractable methods such as hidden Markov model (HMM) tracking, which can accommodate parameter uncertainties and the mildly time-dependent frequency evolution due to slow cloud mass loss (Isi et al., 2018). Analysis pipelines often focus on recently formed, well-localized merger remnants, as the GW signal is expected to be strongest and the cloud youngest in these systems.

High-sensitivity stochastic searches are performed using both stochastic backgrounds and individually resolvable events, with third-generation detectors capable of probing boson masses up to $\sim2 \times 10^{-11}$ eV (Yuan et al., 2021).

6. Astrophysical and Cosmological Constraints

Beyond GW signals, ultralight bosons leave indirect signatures in other astrophysical observables:

• The Event Horizon Telescope's imaging of Sgr A*'s black hole shadow, combined with spin estimates, constrains ULB masses based on absence of significant spin-down—ruling out parameter slices for scalar, vector, and tensor bosons (Saha et al., 2022).
• High-precision astrometry and spectroscopy of stellar orbits (notably S2) near the Galactic Center limit the size and density of possible fuzzy dark matter solitonic cores, translating to bounds $m_a < 1 \times 10^{-19}\,\text{eV}$ at 95% C.L. in combination with cosmological constraints (Monica et al., 2022).
• Modeling of nonlinear structure growth with ultralight bosons (described by Schrödinger–Poisson equations) reveals distinctive small-scale power suppression and smoothing of bumpy features in the power spectrum, consistent with properties of “fuzzy dark matter” (Medellin-Gonzalez et al., 2020).

7. Broader Theoretical and Phenomenological Implications

Ultralight bosons are realized in diverse theoretical settings, including global symmetry breaking (yielding pseudo-Goldstones), minimal extensions producing hidden abelian or non-abelian gauge bosons, and composite models. Their masses are naturally protected by symmetries, rendering $\mu \sim 10^{-20}$–$10^{-10}$ eV technically natural. Gravitational interactions enable the formation of “boson stars” or “Proca stars” and induce “hairy” non-Kerr black holes with observable consequences for gravitational-wave astronomy, strong-gravity tests, and cosmic structure formation (Freitas et al., 2021).

Detection of ultralight bosons would have consequences for dark matter, the nature of quantum gravity, astrophysical compact object populations, and potentially lead to deep insights into the particle content of the universe. Conversely, the absence of evidence in relevant regimes already sets significant limits on both dark matter candidates and new physics scenarios.

---

Table: Observational Probes and Constraints on Ultralight Boson Mass

Probe/Method	Mass Range (eV)	Notable Features / Outcome
GW Direct Search (LIGO)	$2\times10^{-13}$–$10^{-12}$	Up to $10^4$ resolvable events in 4 yrs (Brito et al., 2017)
GW Direct Search (LISA)	$5\times10^{-19}$–$5\times10^{-16}$	Up to $10^3$ events, spins in $10^3$–$10^7 M_\odot$ BHs
Stochastic GW Background (3G)	$7\times10^{-14}$–$2\times10^{-11}$	Higher $m$-modes crucial, extended mass window (Yuan et al., 2021)
Black Hole Spin Demographics	$1.3\times10^{-13}$–$2.7\times10^{-13}$	Constraints via absence of “spin holes” (Ng et al., 2020)
EHT: Sgr A* Spin	Model-dependent	Excludes various boson masses by spin nonzero (Saha et al., 2022)
S2 star astrometry	$<10^{-19}$	Upper bound via orbital stability (Monica et al., 2022)
Structure formation	$10^{-23}$–$10^{-20}$	Power suppression, core-halo phenomenology (Medellin-Gonzalez et al., 2020)

---

The multi-pronged astrophysical and gravitational-wave search program, underpinned by the detailed theoretical modeling of superradiance, cloud evolution, and wave emission, makes the ultralight boson sector an active and robust frontier in fundamental physics.