⚛️ general relativity

Quasinormal modes of rotating black holes with ultralight scalar fields

This paper presents the first accurate computation of quasinormal-mode frequency shifts in rapidly rotating black holes surrounded by massive, nonminimally coupled scalar fields, enabling the use of existing gravitational-wave signals to search for ultralight dark-matter candidates.

Original authors: Leif Lui, Adrian Ka-Wai Chung

Published 2026-10-14

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Leif Lui, Adrian Ka-Wai Chung

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Dark matter is the invisible substance that holds galaxies together, yet its true nature remains one of the most profound mysteries in physics. While we know it exists because of how it pulls on stars and light, we have never directly detected the particles that make it up. One leading theory suggests that dark matter could be composed of ultralight bosons, a type of particle so light that its quantum wave nature stretches across entire star systems. If these particles exist, they should interact with the most extreme objects in the universe: rotating black holes. When a black hole spins, it can act like a cosmic amplifier, potentially trapping these particles in a cloud around it or altering the way the black hole vibrates after a collision. By listening to the specific tones of these vibrations, known as ringdowns, scientists hope to find the first direct evidence of this hidden cosmic component.

For decades, researchers have focused on how these particles might grow into massive clouds around black holes, a process called superradiance. However, a different possibility has emerged: even without a growing cloud, the mere presence of these particles could subtly change the black hole's behavior if they interact with the curvature of space itself. This interaction, known as nonminimal curvature coupling, means the black hole's own gravity could source a configuration of these particles, leaving a permanent mark on the spacetime around it. Until now, calculating exactly how this changes the black hole's vibrations has been a formidable challenge, particularly for rapidly spinning black holes. The complexity of the equations involved has left a significant gap in our ability to search for these particles using current gravitational-wave detectors.

In a new study, researchers have filled this gap by performing the first highly accurate calculations of how these ultralight particles shift the frequencies of a spinning black hole's vibrations. The team developed a novel computational method that bypasses the need for tedious, manual derivation of complex equations, allowing them to efficiently map out how the vibrations change across a wide range of black hole spins and particle masses. They focused on two specific theoretical models where the particles couple to the geometry of space in different ways, simulating black holes spinning at up to ninety percent of the maximum speed allowed by physics. Their results show that the presence of these particles creates distinct, measurable shifts in the pitch and decay rate of the black hole's ringdown signal. These shifts depend heavily on the mass of the particle and the specific way it interacts with gravity, creating a unique fingerprint that distinguishes it from standard black hole behavior.

The researchers found that these shifts are not random but follow precise patterns. For instance, the way the vibration frequency changes as the black hole spins faster differs depending on whether the particle interaction preserves or flips the symmetry of the gravitational waves. They also discovered a critical threshold where the behavior of the particles changes abruptly. When the particle's mass matches a specific value related to the black hole's spin, the nature of the vibration shifts from a traveling wave to a trapped, decaying wave. This transition creates a sudden jump in the frequency shift, a feature that is finite and predictable rather than chaotic. Crucially, the team found no evidence that these interactions cause the black hole to become unstable or explode; instead, the system settles into a new, stable configuration with slightly altered vibrational tones.

These findings provide the essential blueprint for using gravitational waves to hunt for ultralight dark matter. The researchers applied their new data to a recent gravitational-wave event, estimating what future detectors could reveal. They determined that if a signal matches the predictions of standard physics, it would place strict limits on how strongly these particles can interact with gravity, ruling out certain interaction strengths for particles with masses up to a specific threshold. While a single detection might not be enough to pinpoint the exact mass of the particle, it would be sufficient to constrain the possible theories. As detectors become more sensitive and capture louder signals from more distant collisions, the ability to measure these tiny frequency shifts will improve, potentially turning the ringdown of a black hole into a precise instrument for weighing the invisible universe. This work transforms a theoretical possibility into a concrete search strategy, offering a new path to understanding the fundamental building blocks of reality.

Technical Summary: Quasinormal Modes of Rotating Black Holes with Ultralight Scalar Fields

Problem Statement
The nature of dark matter remains a central open question in fundamental physics. While ultralight bosons are well-motivated candidates, most gravitational-wave searches have focused on superradiant instabilities. A distinct possibility arises when light scalar fields couple nonminimally to spacetime curvature. In such theories, the black hole's curvature itself sources a nonvanishing scalar configuration without relying on superradiant growth. While gravitational perturbations and quasinormal-mode (QNM) spectra have been studied for massive minimally-coupled and massless non-minimally-coupled scalar fields, the case of massive non-minimally-coupled scalar fields remains largely unexplored. This gap hinders the ability to search for ultralight curvature-coupled scalar particles using existing gravitational-wave ringdown signals.

Methodology
The authors compute the leading-order shifts in the QNM frequencies of rapidly rotating black holes surrounded by massive scalar fields governed by a Lagrangian including Gauss–Bonnet and Pontryagin invariants (covering dynamical Chern–Simons and scalar Gauss–Bonnet gravity).

Key methodological innovations include:

Background Construction: The authors solve the stationary and axisymmetric background equations using a metric ansatz with spectral collocation and radial domain decomposition to handle the multiple radial scales introduced by the finite scalar mass.
Perturbation Scheme: They perturb both the background spacetime and the massive scalar field. The metric perturbation is described by six functions in a Regge–Wheeler-type gauge, while the scalar perturbation incorporates massive-scalar asymptotic behavior.
Numerical Eigenvalue Solver: Unlike previous implementations that rely on explicit symbolic derivation of linearized field equations, the authors employ a newly developed scheme where matrix coefficients are assembled numerically from the field equations at collocation points.
Eigenvalue Perturbation: The problem is formulated as an eigenvalue problem $D(\omega, \zeta)\mathbf{v} = 0$ . By expanding the eigenfrequency and eigenvector in the dimensionless coupling parameter $\zeta$ , the leading-order frequency shift $\omega^{(1)}$ is obtained from a single linear solve (a bordered least-squares problem).
Boundary Conditions: The method consistently treats the transition between propagating (open) and evanescent (closed) scalar channels, defined by the wave number $k = \sqrt{\omega^2 - \mu^2}$ .

Key Results
The study provides the first accurate computation of leading-order frequency shifts for the $(2,2,0)$ , $(2,2,1)$ , and $(3,3,0)$ modes across dimensionless spins ( $a/M$ ) from 0 to 0.9 and reduced Compton wavelengths comparable to or longer than the black-hole mass scale.

Accuracy: For dimensionless spins between 0 and 0.9, numerical errors in frequency shifts are typically $10^{-8} - 10^{-5}$ . Errors reach at most $O(10^{-3})$ near the scalar-mass threshold.
Parity Dependence: The frequency shifts exhibit strong dependence on parity, demonstrating the breaking of axial–polar isospectrality by the curvature coupling.
- In dynamical Chern–Simons (dCS) gravity, the polar-led shift vanishes as $a \to 0$ because the background becomes Schwarzschild and polar perturbations decouple from the pseudoscalar.
- In scalar Gauss–Bonnet (sGB) gravity, scalar hair exists even at $a=0$ , causing shifts in both parities.
Threshold Transitions: A distinct "step" or jump in the frequency shift occurs when $\Re[\omega^{(0)}] = \mu$ . This corresponds to a transition where the scalar channel closes (becoming evanescent). The jump is finite and primarily affects the damping rate. The physical mode hops between branches at this threshold, while the continued branches remain smooth.
Massless Limit: As the scalar mass $\mu \to 0$ $μ \to 0$ , the results continuously recover the massless limit, but with different scaling behaviors:
- Scalar Gauss–Bonnet: The shift scales linearly with mass ( $\omega^{(1)}(\mu) - \omega^{(1)}(0) \propto \mu$ ), attributed to the monopolar nature of the scalar hair.
- Dynamical Chern–Simons: The shift scales quadratically ( $\propto \mu^2$ ), attributed to the dipolar nature of the scalar hair.
Stability: No indication of instability was found among the modes explored within the parameter space.

Significance and Claims
The authors claim that these results establish the necessary ringdown spectra to search for ultralight curvature-coupled scalar particles in gravitational-wave signals.

Observational Prospects: Using a Fisher-information-matrix analysis of a GW250114-like ringdown (assuming LIGO–Virgo–KAGRA O4 sensitivity), the authors forecast that a single event could constrain the coupling parameter to $\zeta \lesssim 0.04–0.09$ (at $2\sigma$ ) for scalar masses $\mu \lesssim 2 \times 10^{-12}$ eV. However, they note that measuring the scalar mass itself would require louder signals or next-generation detectors, as the uncertainty in $\mu$ for a single event is wider than the range of interest.
Methodological Advance: The paper highlights a computational advance in black-hole perturbation theory. By avoiding explicit symbolic derivation and assembling operators numerically, the scheme reduces computational overhead, allowing spectra for hundreds of parameter points to be generated in hours. This facilitates systematic exploration of black-hole spin, scalar mass, and coupling structure.
Future Work: The authors state that their numerical results provide the ingredients for modeling extreme-mass-ratio inspirals around black holes with massive scalar hair. They emphasize that their forecasts are indicative and that a full Bayesian parameter-estimation analysis is required to turn these forecasts into robust constraints.

The paper concludes that these findings open a new avenue for probing ultralight scalar particles directly through black-hole ringdown observations, extending the understanding of gravitational perturbations of rotating black holes carrying scalar matter.

More like this