Black-Hole Quasinormal Modes (QNMs) Definition, and Applications

• Black-Hole quasinormal modes (QNMs) describe complex-frequency solutions to the perturbation equations of black holes, encoding parameters of the background geometry and serving as indicators of stability, modified gravity, and environmental effects.
• Applications of QNMs include testing the no-hair theorem by providing precise frequencies and damping times of black holes; enabling the reconstruction of spatial structure and understanding gravitational wave contributions.
• Notably, QNMs are used in the field of AdS/CFT correspondence, where they map onto poles of retarded boundary correlators, effectively describing transport and relaxation in strongly coupled quantum systems, thus providing a crucial link between different areas of theoretical physics and astrophysics.

Black-hole quasinormal modes (QNMs) are discrete, generally complex-frequency solutions of the linearized perturbation equations of a black-hole spacetime subject to purely ingoing boundary conditions at the future event horizon and purely outgoing conditions at infinity, or to the appropriate reflective conditions in asymptotically anti-de Sitter (AdS) geometries. With the convention $e^{-i\omega t}$, $\omega=\omega_R+i\omega_I$, where $\omega_R$ is the oscillation frequency and stable modes satisfy $\omega_I<0$; the damping time is $\tau=1/|\omega_I|$. QNMs describe the intermediate-time ringdown of perturbed black holes, encode the parameters of the background geometry, and provide a spectral framework for black-hole spectroscopy, stability analysis, modified-gravity tests, environmental studies, nonlinear perturbation theory, and holographic transport (0905.2975).

1. Definition, spectral structure, and physical interpretation

A perturbed black hole is represented schematically by

$$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$

where $g_{\mu\nu}^{(0)}$ is a stationary background such as Schwarzschild, Kerr, Reissner–Nordström, or an asymptotically AdS or de Sitter geometry. After separating temporal and angular dependence, many perturbation problems reduce to a radial equation of Schrödinger type,

$$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$

where $r_*$ is a tortoise coordinate and $V$ is an effective potential.

For asymptotically flat black holes, the QNM conditions are

$\omega=\omega_R+i\omega_I$0

at the future horizon, and

$\omega=\omega_R+i\omega_I$1

at spatial infinity. The first condition excludes radiation emerging from the future horizon; the second excludes radiation incident from infinity. These one-way conditions make the spectral problem non-self-adjoint. QNMs are consequently resonances of an open dissipative system rather than square-integrable normal modes.

The frequency is conventionally labeled by angular indices and an overtone number,

$\omega=\omega_R+i\omega_I$2

where $\omega=\omega_R+i\omega_I$3 is the multipole index, $\omega=\omega_R+i\omega_I$4 is the azimuthal index, and $\omega=\omega_R+i\omega_I$5 is the overtone number. The fundamental mode $\omega=\omega_R+i\omega_I$6 is usually the least damped, whereas higher overtones generally have larger $\omega=\omega_R+i\omega_I$7. Mode frequencies are properties of the background and boundary conditions; excitation amplitudes and phases depend on the initial perturbation or source.

The response to a generic disturbance has three regimes: a prompt response determined by the initial data, an intermediate ringdown dominated by QNM residues, and a late-time tail. In asymptotically flat spacetimes, branch-cut contributions to the Green function produce power-law tails that eventually dominate over exponentially damped QNMs. QNM expansions therefore describe a finite ringdown interval rather than, in general, a complete expansion of arbitrary initial data.

QNMs are poles of the analytically continued Green function. If $\omega=\omega_R+i\omega_I$8 denotes the incoming asymptotic amplitude, the QNM condition can be written

$\omega=\omega_R+i\omega_I$9

The time-domain response contains prompt propagation, residues at QNM poles, and branch-cut contributions. This Green-function interpretation explains why the QNM spectrum is discrete despite the non-Hermitian character of the radial operator.

2. Perturbation equations and canonical black-hole sectors

Schwarzschild perturbations

The Schwarzschild metric is

$\omega_R$0

with tortoise coordinate

$\omega_R$1

For scalar, electromagnetic, and axial gravitational perturbations, a compact potential is

$\omega_R$2

where $\omega_R$3 denotes scalar, electromagnetic, and axial gravitational perturbations. Thus,

$\omega_R$4

$\omega_R$5

and the axial gravitational potential is the Regge–Wheeler potential,

$\omega_R$6

Even-parity gravitational perturbations are described by the Zerilli function and potential,

$\omega_R$7

The Regge–Wheeler and Zerilli potentials are related by a supersymmetric-type transformation and are isospectral under the standard Schwarzschild QNM boundary conditions, apart from subtleties involving algebraically special and exceptional frequencies. Isospectrality concerns the frequencies, not equality of the master variables, eigenfunctions, waveforms, or excitation amplitudes (Zhao et al., 2022).

Representative Schwarzschild gravitational frequencies are

$\omega_R$8

with more precise values

$\omega_R$9

and

$\omega_I<0$0

The fundamental scalar $\omega_I<0$1 and electromagnetic $\omega_I<0$2 frequencies are approximately

$\omega_I<0$3

respectively.

In the eikonal regime, $\omega_I<0$4, the Schwarzschild spectrum is governed by the unstable circular photon orbit at $\omega_I<0$5: $\omega_I<0$6 Equivalently,

$\omega_I<0$7

where $\omega_I<0$8 is the photon-orbit frequency and $\omega_I<0$9 its Lyapunov exponent.

Kerr perturbations

For Kerr,

$\tau=1/|\omega_I|$0

and perturbations are described by the Teukolsky formalism. The separated field contains spin-weighted spheroidal harmonics $\tau=1/|\omega_I|$1 and radial Teukolsky functions $\tau=1/|\omega_I|$2. The angular eigenvalue $\tau=1/|\omega_I|$3 and radial QNM frequency must generally be determined simultaneously.

The Kerr horizon condition involves the horizon-frame frequency

$\tau=1/|\omega_I|$4

Near the horizon, the radial field behaves schematically as

$\tau=1/|\omega_I|$5

while at infinity it is outgoing,

$\tau=1/|\omega_I|$6

Kerr QNMs exhibit $\tau=1/|\omega_I|$7-dependent rotational splitting, distinct co-rotating and counter-rotating branches, superradiant behavior controlled by $\tau=1/|\omega_I|$8, and mode clustering near $\tau=1/|\omega_I|$9 as $$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$0. For many co-rotating modes, increasing spin raises the oscillation frequency and reduces the damping rate. The fundamental gravitational $$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$1 mode is the dominant ringdown mode in many mergers.

Other backgrounds

Reissner–Nordström perturbations couple gravitational and electromagnetic degrees of freedom but can be reduced to two master sectors. Schwarzschild–de Sitter perturbations satisfy ingoing conditions at the black-hole horizon and outgoing conditions at the cosmological horizon. Near extremality, the effective potential becomes approximately Pöschl–Teller,

$$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$2

with an analytically tractable QNM spectrum.

In higher dimensions, the Ishibashi–Kodama formalism classifies perturbations into scalar-, vector-, and tensor-type sectors. The tensor sector has no four-dimensional analogue. Black strings and black branes additionally support Gregory–Laflamme instabilities in suitable long-wavelength sectors.

For black branes in AdS, QNMs depend on continuous spatial momentum rather than a discrete angular number. In gauge/gravity duality, bulk QNMs coincide with poles of retarded boundary correlators. Hydrodynamic QNMs include shear diffusion,

$$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$3

sound propagation,

$$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$4

and charge diffusion,

$$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$5

For conformal $$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$6 supersymmetric Yang–Mills theory at strong coupling,

$$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$7

Thus QNMs encode transport coefficients, diffusion, conductivity, sound attenuation, and equilibration in strongly coupled finite-temperature quantum theories (Konoplya et al., 2011).

3. Numerical and analytical computation

No single method is optimal for all backgrounds, multipoles, overtones, and asymptotic regions.

Continued fractions

Leaver’s continued-fraction method factors out the known horizon and infinity behavior in a Frobenius expansion. A typical three-term recurrence is

$$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$8

The QNM condition is an infinite continued fraction,

$$g_{\mu\nu}=g_{\mu\nu}^{(0)}+\delta g_{\mu\nu}, \qquad |\delta g_{\mu\nu}|\ll |g_{\mu\nu}^{(0)}|,$$9

The condition selects the minimal solution of the recurrence. Continued fractions are among the most accurate methods for Schwarzschild, Kerr, and Reissner–Nordström QNMs, although convergence becomes difficult for highly damped modes, high overtones, and near-extremal horizons. Nollert-type large-index expansions improve convergence in highly damped regimes.

WKB and solvable potentials

WKB methods expand near the maximum of a single barrier potential. At leading order,

$g_{\mu\nu}^{(0)}$0

with higher-order corrections involving derivatives of the potential at its maximum. WKB is most effective for low overtones, large $g_{\mu\nu}^{(0)}$1, and single-barrier potentials. It becomes unreliable for high overtones, multiple barriers, near-extremal systems, and strongly non-barrier potentials. Sixth- or ninth-order WKB improves fundamental-mode estimates but does not remove the asymptotic nature of the approximation.

The Pöschl–Teller potential provides an exactly solvable approximation in near-extremal Schwarzschild–de Sitter and related geometries. Monodromy and phase-integral methods analytically continue the radial equation into the complex plane and are especially useful for highly damped modes. For Schwarzschild, they yield the asymptotic relation

$g_{\mu\nu}^{(0)}$2

for scalar and gravitational perturbations, whereas electromagnetic perturbations have asymptotic real part tending to zero.

Asymptotic iteration method

The asymptotic iteration method (AIM) starts from

$g_{\mu\nu}^{(0)}$3

Repeated differentiation produces

$g_{\mu\nu}^{(0)}$4

with

$g_{\mu\nu}^{(0)}$5

$g_{\mu\nu}^{(0)}$6

The termination condition is

$g_{\mu\nu}^{(0)}$7

Improved AIM replaces repeated differentiation by Taylor expansions around a fixed point. It avoids Gaussian elimination when the underlying recurrence has more than three terms and can treat Schwarzschild, Reissner–Nordström, Kerr, higher-dimensional, and multiply rotating geometries. Its convergence depends strongly on the expansion point; high overtones, near-extremal horizons, and coupled Kerr radial-angular problems remain difficult (Cho et al., 2011).

Pseudospectral and direct methods

Pseudospectral methods discretize the radial problem at Chebyshev collocation points and convert it into a generalized eigenvalue problem,

$g_{\mu\nu}^{(0)}$8

They are useful for coupled systems, nonstandard geometries, higher-derivative theories, and high overtones. Eddington–Finkelstein formulations can incorporate horizon regularity and asymptotic conditions directly.

Time-domain evolution is useful when separation is unavailable or when the dominant response and stability must be identified directly. The waveform is fitted to damped exponentials,

$g_{\mu\nu}^{(0)}$9

Time-domain fits are affected by prompt response, tails, finite-domain reflections, mode interference, and the rapid decay of high overtones.

Characteristic and Bondi–Sachs formulations

QNMs can also be computed directly in characteristic formulations. In the Bondi–Sachs approach, a fourth-order equation for a metric variable can be reduced to a second-order equation for

$$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$0

A Frobenius expansion then yields a three-term recurrence and a continued-fraction condition. This formulation recovers the standard Schwarzschild spectrum and naturally represents outgoing radiation at null infinity. It also permits a distinct treatment of the algebraically special mode, which satisfies outgoing conditions at both the horizon and infinity (Mongwane et al., 2024).

4. Ringdown spectroscopy and spatial QNM structure

After a binary merger or collapse, the remnant approaches a stationary black hole. The gravitational waveform consists of prompt merger dynamics, an intermediate QNM-dominated ringdown, and, in asymptotically flat spacetimes, a late-time tail. For Kerr, the ringdown strain can be expressed schematically as

$$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$1

The frequencies scale with the remnant mass and dimensionless spin: $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$2 Consequently, measuring a mode frequency and damping time can infer $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$3 and $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$4. Measuring two or more independent modes supplies a test of the Kerr spectrum and the no-hair theorem. The waveform amplitudes, however, depend on the merger history, mass ratio, spins, orientation, precession, and excitation coefficients.

Kerr spheroidal harmonics are not identical to spin-weighted spherical harmonics. They admit an expansion

$$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$5

where the mixing coefficients preserve $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$6 but couple different spherical $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$7 values. This spheroidal–spherical mixing must be retained in simultaneous multimode fits.

The black-hole cartography program reconstructs the angular shape of a QNM by fitting a fixed QNM time dependence to many spherical-harmonic components while allowing the angular coefficients of the selected mode to vary independently. For linear Kerr modes, the reconstructed shapes agree with the predicted spheroidal harmonics. In one equal-mass nonspinning numerical-relativity simulation, the reconstructed dominant $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$8 mode had spatial mismatch

$$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r,\omega)\right]\Psi=0,$$9

and its contributions to $r_*$0 and $r_*$1 agreed with the corresponding spheroidal–spherical mixing coefficients (Dyer et al., 2024).

Spatial reconstruction also applies to nonlinear quadratic modes. A quadratic mode sourced by parent modes $r_*$2 and $r_*$3 has frequency

$r_*$4

and azimuthal number

$r_*$5

For the dominant self-coupling of the $r_*$6 mode,

$r_*$7

The angular structure of quadratic QNMs is not described by a unique established family of spheroidal harmonics. Candidate prescriptions include products with total spin weight $r_*$8, a spheroidal harmonic at the summed frequency, and spin-weight raising of a direct product. Current numerical-relativity reconstructions find qualitatively similar results for several prescriptions, without decisively selecting one.

The viewing angle affects the detectability of quadratic modes. For a remnant with $r_*$9, the reconstructed dominant quadratic mode has maximum projected amplitude near

$V$0

with values across simulations and prescriptions roughly in the range

$V$1

The strongest quadratic signal therefore need not be observed along the remnant spin axis.

5. Nonlinear, environmental, and modified-background QNMs

Second-order QNMs

At second perturbative order,

$V$2

and the Einstein equations yield

$V$3

The second-order problem is linear in $V$4 but has a source quadratic in the first-order perturbation.

For dominant even-parity Schwarzschild modes $V$5, angular-momentum selection rules generate oscillatory $V$6 modes and a nonoscillatory $V$7 component. If

$V$8

then

$V$9

so the driven quadratic mode has

$\omega=\omega_R+i\omega_I$00

This is a combination frequency, not the frequency of a freely existing homogeneous linear $\omega=\omega_R+i\omega_I$01 QNM. The second-order amplitude can reach approximately $\omega=\omega_R+i\omega_I$02 of the first-order ringdown amplitude for strongly excited binary-black-hole mergers, subject to the excitation, orientation, overtone content, and normalization conventions (0708.0450).

The raw second-order source requires regularization at the horizon and infinity. A shifted master function can transform a source behaving as $\omega=\omega_R+i\omega_I$03 at infinity into one behaving as $\omega=\omega_R+i\omega_I$04. The second-order calculation uses a modified Leaver continued-fraction method, with an inhomogeneous recurrence relation generated by the quadratic source.

Exceptional points and resonant excitation

QNMs form a non-Hermitian spectrum, and mode interactions can occur when a background parameter varies. Near an exceptional point (EP), two eigenvalues and eigenfunctions coalesce. A two-mode model has eigenvalues

$\omega=\omega_R+i\omega_I$05

with EP condition

$\omega=\omega_R+i\omega_I$06

A one-parameter scan generally encounters an avoided crossing near, rather than exactly at, the EP.

QNMs require regularized biorthogonal products rather than ordinary positive norms. If $\omega=\omega_R+i\omega_I$07 denotes a regularized QNM norm, the excitation factor obeys

$\omega=\omega_R+i\omega_I$08

Near an EP, $\omega=\omega_R+i\omega_I$09 becomes small, causing large excitation factors and enhanced sensitivity to perturbations. Frequency trajectories become approximately hyperbolic, whereas excitation-factor trajectories follow an inverse lemniscate-like structure. The resonance profile is described by a quarter-power Lorentzian,

$\omega=\omega_R+i\omega_I$10

High-precision Kerr calculations and double-barrier models exhibit avoided crossings, mode mixing, loops, frequency repulsion, and correlated excitation-factor amplification (Motohashi, 2024).

Deformed black holes and environmental matter

Small deviations from Kerr can be treated perturbatively: $\omega=\omega_R+i\omega_I$11 A contour-regularized QNM inner product gives

$\omega=\omega_R+i\omega_I$12

The method remains applicable when the deformation destroys separability, although nearly degenerate modes require a coupled or degenerate treatment (Zimmerman et al., 2014).

A self-gravitating thin disk surrounding Schwarzschild modifies the scalar effective potential at first order in the disk-to-black-hole mass ratio,

$\omega=\omega_R+i\omega_I$13

For physically reasonable disk configurations, the disk lowers both the oscillation frequency and damping magnitude: $\omega=\omega_R+i\omega_I$14 The frequency shifts approximately obey

$\omega=\omega_R+i\omega_I$15

This empirical correlation can serve as an environmental signature within the perturbative single-barrier regime. It fails when the disk is too massive or concentrated and generates multiple potential peaks, trapped modes, or echoes (Chen et al., 2023).

Other modified backgrounds exhibit distinct spectral responses. In a MOG black hole, increasing the MOG parameter at fixed dimensionless mass lowers both real and imaginary parts of the QNM frequency, whereas comparison at fixed effective gravitational mass produces higher real frequencies and nearly GR-like damping (Manfredi et al., 2017). Quantum-corrected Schwarzschild models generally produce lower oscillation frequencies and slower damping in the parameter ranges studied (Saleh et al., 2016). Loop-quantum-gravity-corrected and regular black-hole models can leave fundamental scalar modes close to Schwarzschild while producing large shifts in higher overtones, including oscillatory or spiral-like trajectories near extremality (Gong et al., 2023, Zhang et al., 2024).

6. Stability, limitations, and unresolved problems

A mode is stable when its time dependence decays. With $\omega=\omega_R+i\omega_I$16, this requires

$\omega=\omega_R+i\omega_I$17

A mode with $\omega=\omega_R+i\omega_I$18 signals a linear instability. Stability must be assessed across the full spectrum rather than from the fundamental mode alone.

Schwarzschild and Kerr black holes are linearly stable in the sectors and regimes discussed. Reissner–Nordström and several Schwarzschild–de Sitter and Schwarzschild–AdS systems are likewise stable in their established perturbation sectors. Instabilities arise in other settings: Gregory–Laflamme instabilities affect black strings and branes; superradiance can become unstable when amplification is combined with confinement, as in Kerr–AdS or massive-field bound states; and higher-dimensional or higher-curvature black holes can possess nonspherical and delayed instabilities.

Several limitations constrain the interpretation of QNM spectra. QNM eigenfunctions are not ordinary normalizable states, asymptotic expansions may be incomplete, and late-time tails are not captured by a finite QNM sum. WKB methods are regime-dependent, continued fractions can become delicate near extremality and high damping, AIM depends on an expansion point, and pseudospectral methods require careful control of boundary conditions and spurious eigenvalues. Direct shooting is destabilized by exponentially growing unwanted solutions.

Observationally, high overtones are both potentially informative and difficult to measure. They can have substantial amplitude near the merger peak but decay rapidly and are susceptible to contamination by nonlinear merger dynamics, detector noise, calibration errors, and start-time choices. The interpretation of overtone evidence in GW150914 remains sensitive to these issues. A flexible overtone model can absorb residual merger structure or noise, so robust inference requires stable posterior estimates, multimode consistency, and comparisons across analysis methods (Zhao et al., 2022).

Important open problems include complete decoupling and stability analyses of Kerr–Newman systems, generic higher-dimensional rotating black holes, and higher-curvature theories; nonlinear evolution of black-brane and Gauss–Bonnet instabilities; the ultimate form of late-time nonlinear tails; QNM spectra with environmental and quantum corrections; and observationally resolving multiple modes with sufficient precision to test the no-hair theorem.

Black-hole QNMs consequently form a unified language for open gravitational systems. In astrophysics, their frequencies and damping times describe ringdown and test the Kerr geometry. In nonlinear perturbation theory, quadratic combinations generate predictable combination modes. In modified geometries and environmental systems, QNM shifts probe near-horizon structure and surrounding matter. In characteristic formulations, they connect directly with outgoing radiation at null infinity. In AdS/CFT, they become poles of retarded correlators and encode relaxation and transport in strongly coupled quantum field theories.