Quasinormal Modes in Black Hole Physics

• Quasinormal Modes (QNMs) are dissipative resonances of open dynamical systems, characterized by complex frequencies that encode oscillation and decay rates, with applications in black hole physics and electromagnetic resonators.
• QNMs provide vital information about a system's boundary conditions, exterior potential barriers, and global properties, uniquely different from normal modes in conservative systems.
• They are analyzed through diverse mathematical methods, including perturbation theory, WKB approximations, continued fractions, time-domain evolution, and algebraic approaches, each offering unique insights and computational advantages.

Quasinormal modes (QNMs) are dissipative resonances of open dynamical systems. In black-hole physics, they are linear perturbations satisfying a purely ingoing condition at a future event horizon and a purely outgoing condition at infinity or at the relevant outer asymptotic boundary. Their complex frequencies encode oscillation and decay: with time dependence $e^{-i\omega t}$, $\operatorname{Re}\omega$ is the oscillation frequency and $\operatorname{Im}\omega<0$ denotes temporal damping. QNMs are not ordinary normal modes of a conservative system: their radial or spatial profiles are generally non-square-integrable, and the spectrum is defined through boundary conditions, resonances of meromorphically continued resolvents, or eigenvalues of suitably constructed time-evolution generators.

1. Definition, spectral structure, and physical interpretation

A QNM is a source-free solution of a linearized field equation subject to radiative boundary conditions. For a one-dimensional radial equation of Schrödinger type,

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

the standard asymptotically flat conditions are

$\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$

at the future horizon, and

$\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$

at spatial infinity. The first branch is ingoing into the horizon and the second is outgoing toward infinity. Only discrete complex frequencies satisfy both conditions.

With the convention $e^{-i\omega t}$,

$$e^{-i\omega t}=e^{-i\operatorname{Re}\omega t}e^{\operatorname{Im}\omega t}.$$

Thus $\operatorname{Im}\omega<0$ describes decay, $\operatorname{Im}\omega>0$ describes growth, and $\operatorname{Re}\omega$0 gives an undamped oscillation. An alternative convention, frequently used in the literature, is $\operatorname{Re}\omega$1 with $\operatorname{Re}\omega$2 for damping, or $\operatorname{Re}\omega$3 with $\operatorname{Re}\omega$4 and $\operatorname{Re}\omega$5 for decay. The damping time is approximately $\operatorname{Re}\omega$6.

QNMs depend on the background geometry and on the perturbation sector. Relevant parameters include mass, electric charge, angular momentum, horizon locations and temperatures, inner-horizon structure, field spin and mass, angular momentum number, overtone number, and—in charged-field problems—the field charge. The spectrum therefore provides information about both the exterior potential barrier and the global boundary-value problem.

QNMs are poles of analytically continued Green functions or resolvents. In electromagnetic resonators, for example, the Green tensor has a pole expansion of the form

$\operatorname{Re}\omega$7

where $\operatorname{Re}\omega$8 and $\operatorname{Re}\omega$9 is the decay rate (Kristensen et al., 2019). In rigorous black-hole formulations, QNFs are poles of a meromorphic inverse of a stationary operator or isolated eigenvalues of a time-translation generator. For Kerr–AdS, this construction applies without separation of variables and allows broad elliptic boundary conditions at the timelike conformal boundary (Gannot, 2014). For sub-extremal Kerr, complex scaling at infinity and microlocal radial-point analysis at the horizon yield a cutoff-resolvent definition (Stucker, 2024).

The distinction between a QNM and a normal mode is fundamental. A normal mode is spatially localized, square-integrable, and has a real frequency in a conservative system. A QNM is outgoing and therefore generally has a complex frequency; its spatial profile grows under frequency-domain analytic continuation even though its time-dependent physical signal decays. This behavior occurs in black-hole perturbations, open electromagnetic resonators, and metastable excitations of nonlinear field theories.

2. Perturbation equations and effective potentials

For static, spherically symmetric geometries,

$\operatorname{Im}\omega<0$0

the tortoise coordinate is defined by

$\operatorname{Im}\omega<0$1

Massless neutral scalar perturbations separated as

$\operatorname{Im}\omega<0$2

obey

$\operatorname{Im}\omega<0$3

with

$\operatorname{Im}\omega<0$4

The matter model enters this test-scalar problem through $\operatorname{Im}\omega<0$5 alone. Consequently, distinct theories producing the same lapse function have identical neutral test-scalar spectra (Flachi et al., 2012).

For spin $\operatorname{Im}\omega<0$6 in ungravity-enhanced black-hole geometries, the effective potentials are

$\operatorname{Im}\omega<0$7

$\operatorname{Im}\omega<0$8

and

$\operatorname{Im}\omega<0$9

The two Dirac potentials are supersymmetric partner potentials. The resulting spectrum resembles that of a Schwarzschild black hole in an effective noninteger dimension, but tensor-unparticle unitarity imposes $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r)\right]\Psi=0,$$0, excluding many formally computed modes that would be allowed in ordinary higher-dimensional Schwarzschild geometries (Lee, 2011).

For gravitational perturbations of Lovelock black holes, tensor, vector, and scalar sectors reduce to

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

At large angular momentum,

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

where $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r)\right]\Psi=0,$$3 depends on the perturbation sector, Lovelock order, and spacetime dimension. The sectors become increasingly similar as the dimension grows, producing approximate isospectrality at sufficiently large dimension (Prasobh et al., 2014).

Charged scalar perturbations contain the electromagnetic coupling through the frequency combination

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

rather than through a simple additive potential. The event-horizon condition becomes

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

while the asymptotic condition remains outgoing. Regular black-hole backgrounds with nonlinear electrodynamics exhibit spectra qualitatively similar to Schwarzschild and Reissner–Nordström geometries, but regularization parameters modify the potential, horizon structure, damping rates, and inner-horizon contributions (Flachi et al., 2012).

Rotating spacetimes generally lead to coupled angular and radial eigenvalue problems. For Kerr, the Teukolsky equation separates into spin-weighted spheroidal harmonics and radial equations. For the rotating C-metric, acceleration makes the angular separation constant frequency-dependent, $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r)\right]\Psi=0,$$6, so the angular and radial continued fractions must be solved simultaneously (Xiong et al., 2023). The separated equations of uniformly accelerated rotating black holes can be transformed into general Heun equations, providing a framework for arbitrary massless spin weights including $$\frac{d^2\Psi}{dr_*^2}+\left[\omega^2-V(r)\right]\Psi=0,$$7 (Chen et al., 2024).

3. Eikonal modes and the geometry of null orbits

In the large-angular-momentum regime, QNMs are controlled by unstable bound null geodesics. For a static black hole, the unstable circular null orbit satisfies

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

Its orbital frequency is

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

and its instability is characterized by a Lyapunov exponent $\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$0. The leading eikonal spectrum is

$\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$1

Hence

$\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$2

The real part measures phase accumulation around the null orbit, whereas the imaginary part measures leakage caused by its instability [140.5334].

For Schwarzschild,

$\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$3

and the eikonal scalar spectrum is

$\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$4

A Penrose limit gives a plane-wave geometry near the relevant null orbit. The transverse perturbation problem becomes an oscillator system: stable transverse directions generate real frequency splittings, while unstable directions generate the imaginary part. For Kerr, generic spherical photon orbits produce periodic transverse dynamics described by Floquet theory; the stable Floquet exponent gives a precession frequency and the unstable one gives a Lyapunov exponent (Fransen, 2023).

In Kerr, the eikonal spectrum of equatorial modes has the form

$\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$5

For generic non-equatorial spherical orbits, the leading real frequency is quantized by an EBK condition and the transverse dynamics is determined by Floquet exponents. The approximation improves for larger angular momentum and lower overtones, and it captures the approach to zero-damped modes for corotating near-extremal Kerr orbits.

Rotating accelerating black holes possess three qualitatively distinct spectral families: photon-sphere modes, acceleration modes, and near-extreme modes. Photon-sphere modes reduce to Kerr modes when acceleration vanishes and are associated with unstable null trajectories. Acceleration modes are controlled by the acceleration horizon and have damping proportional to the Rindler acceleration-horizon surface gravity. Near-extreme modes are associated with the near-horizon region of rapidly rotating black holes and become long-lived as the inner and event horizons approach one another (Xiong et al., 2023, Chen et al., 2024).

4. Analytic, numerical, and algebraic methods

The third-order and sixth-order WKB methods approximate QNMs by expanding about the maximum of a smooth potential barrier. At third order, the schematic quantization condition is

$\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$6

where $\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$7 and $\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$8 are evaluated at the barrier maximum. Sixth-order WKB adds further derivative corrections. WKB is most reliable for low overtones, moderate or large angular momentum, and a single smooth barrier; it becomes less reliable for $\Psi\sim e^{-i\omega r_*},\qquad r_*\to-\infty,$9, high overtones, near-extremal geometries, or multi-peak potentials (Flachi et al., 2012, Prasobh et al., 2014).

Continued-fraction methods impose local Frobenius behavior at the relevant horizons and convert the resulting recurrence into a minimal-solution condition. They are central to Kerr calculations and to the rotating C-metric, where angular and radial continued fractions determine $\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$0 and $\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$1 simultaneously. Leaver’s continued-fraction method provides the zeroth-order Kerr spectrum in perturbative analyses of weakly charged Kerr–Newman spacetimes (Mark et al., 2014).

Direct integration evolves solutions initialized with the required local behavior at each boundary and imposes a vanishing Wronskian-like determinant at an interior matching point. This method provides an independent check on continued fractions but can become unstable near extremality or when the initial frequency guess is poor. Heun-function shooting supplies another boundary-matching method for accelerating black holes (Chen et al., 2024).

Time-domain evolution identifies QNMs from the damped oscillatory phase of a signal. In near-extremal RN–AdS calculations, horizon-penetrating spacelike foliations remain regular as the surface gravity tends to zero and permit initial data nontrivial at the future horizon. The time-domain signal may contain an initial transient, a QNM ringing phase, and a late-time tail. This formulation reveals a transition in which a purely damped mode becomes less damped than oscillatory modes near extremality; at exact extremality, the exponential tail is replaced by power-law decay (Ficek et al., 2023).

Algebraic methods are available in special geometries. Higher-spin de Sitter QNMs can be constructed as lowest-weight modules of $\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$2. Massive fields produce two towers with frequencies

$\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$3

For massless higher-spin fields, gauge symmetry removes the naive second massive tower and replaces it with a gauge-invariant boundary-curvature tower. The resulting QNM character agrees with the Harish–Chandra character of the corresponding de Sitter representation (Sun, 2020).

Perturbative eigenvalue methods treat a deformed Kerr metric as

$\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$4

For a non-Hermitian QNM problem, ordinary Hilbert-space perturbation theory fails because QNM wavefunctions are not square-integrable. A contour-defined bilinear form gives the first-order shift

$\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$5

This method applies to weakly deformed Kerr geometries even when the deformation destroys separability. For weakly charged Kerr–Newman black holes, gravitational and electromagnetic sectors couple at order $\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$6, but the induced field in the initially absent sector affects the frequency only at order $\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$7; the first-order frequency shift is determined by the diagonal operator correction (Zimmerman et al., 2014, Mark et al., 2014).

5. Rigorous definitions, non-Hermiticity, and boundary conditions

The non-Hermitian character of QNM problems requires more than a formal outgoing-wave prescription. In electromagnetic resonators, right QNMs are outgoing while adjoint QNMs are incoming under the analytically continued radiation condition. The appropriate product is an unconjugated bilinear form with a surface term:

$\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$8

The surface term compensates for the exponential spatial growth of the outgoing QNM. Equivalent normalization schemes use perfectly matched layers, analytic regularization, resonant-state expansions, or Green-tensor residues (Kristensen et al., 2019).

For asymptotically flat black holes, outgoing infinity is an irregular singular endpoint. In sub-extremal Kerr, asymptotically hyperboloidal foliations bring future null infinity to a finite boundary. The rescaled field $\Psi\sim e^{+i\omega r_*},\qquad r_*\to+\infty,$9 has a finite radiation field there. A Gevrey-regular Hilbert space controls infinitely many derivatives with factorial weights, selecting the appropriate outgoing branch. The time-evolution semigroup on this space has a generator $e^{-i\omega t}$0; regularity QNFs are isolated eigenvalues of $e^{-i\omega t}$1, while scattering resonances are poles of the meromorphically continued cutoff resolvent (Gajic et al., 2024).

A related rigorous construction for extremal Reissner–Nordström requires Gevrey regularity at both the horizon and null infinity. Ordinary Sobolev spaces are insufficient because extremality removes the red-shift effect: smooth mode solutions can exist throughout the left half-plane and do not define a discrete spectrum. The QNM set is discrete in the sector

$e^{-i\omega t}$2

and scattering resonances are contained in the regularity QNF set. Conserved quantities at the extremal horizon and null infinity generate polynomial tails in addition to exponentially damped resonant contributions (Gajic et al., 2019).

For Kerr–AdS, the conformal boundary is timelike rather than outgoing. Boundary conditions must therefore be specified explicitly. Massive scalar fields have asymptotic behaviors

$e^{-i\omega t}$3

When $e^{-i\omega t}$4, finite-energy conditions select the admissible branch. For $e^{-i\omega t}$5, both branches are admissible and an elliptic boundary operator involving weighted traces is required. Different elliptic boundary conditions generally produce different QNM spectra. The QNFs are poles of a meromorphic inverse of the stationary Klein–Gordon operator, with finite-rank residues (Gannot, 2014).

Complex scaling provides an alternative rigorous treatment of sub-extremal Kerr. The spatial contour is deformed only at large radius, where outgoing waves become square-integrable. Microlocal radial-point estimates impose the horizon-regular condition, and normally hyperbolic trapping controls high-energy propagation. The cutoff resolvent extends meromorphically on a logarithmic cover of the frequency plane; its finite-rank poles define Kerr QNMs. The construction excludes the upper half-plane and establishes a high-energy resonance-free strip (Stucker, 2024).

6. Spectral families, stability, and physical applications

Black-hole spectra often contain distinct mode families with different geometric origins. In rotating accelerating black holes, photon-sphere modes are associated with the unstable null orbit, acceleration modes with the acceleration horizon, and near-extreme modes with the near-horizon region. These families can undergo eigenvalue repulsion: when acceleration and near-extreme branches approach one another, they avoid crossing, exchange character, and reconnect. Family labels and overtone assignments are therefore not globally invariant (Xiong et al., 2023).

Near extremality, the least damped mode need not be oscillatory. In RN–AdS, purely damped modes can overtake oscillatory modes and dominate late-time evolution. As extremality is reached, these frequencies approach zero and the limiting decay becomes polynomial rather than exponential. This is a nonuniform limit: every fixed subextremal black hole has exponentially decaying QNMs, while the exactly extremal geometry supports horizon-related power-law behavior (Ficek et al., 2023).

QNMs can also diagnose modified geometries. Ungravity black holes mimic Schwarzschild geometries in a fractional effective dimension $e^{-i\omega t}$6, but the tensor-unparticle unitarity condition $e^{-i\omega t}$7 removes modes that would exist in ordinary higher-dimensional black holes. The allowed-mode pattern therefore supplies a proposed spectral distinction between ungravity and extra dimensions (Lee, 2011).

Near-horizon modifications illustrate a limitation of interpreting the early ringdown solely through the exact QNM spectrum. Replacing the horizon by a second barrier produces a dense set of weakly damped cavity modes rather than the original Schwarzschild poles. Nevertheless, before a signal can travel to the modification and return, a coherent sum over these new modes reconstructs the ordinary GR ringdown. The changed spectrum becomes visible later through echoes or other delayed structure (Mirbabayi, 2018).

QNMs also occur outside black-hole perturbation theory. In open electromagnetic resonators, a few normalized QNMs can model near-field enhancement, scattering, Purcell factors, and temporal coupled-mode dynamics. The quality factor of an isolated electromagnetic resonance is

$e^{-i\omega t}$8

Because QNM fields diverge spatially under outgoing analytic continuation, modal expansions have restricted regions of convergence and should not be evaluated naively in the far field. Green-tensor residues, surface-equivalence methods, PMLs, and non-QNM continuum contributions may be required (Kristensen et al., 2019).

In nonlinear kink–antikink scattering, a narrow QNM can temporarily store energy and return part of it to translational motion. Unlike a true internal mode, it radiates during the interval between collisions. Increasing its decay rate closes resonance windows and raises the critical collision velocity. A damped collective-coordinate model,

$e^{-i\omega t}$9

captures this qualitative mechanism (Dorey et al., 2017).

The interpretation of QNM stability must remain sector- and approximation-dependent. Negative imaginary parts for sampled scalar, electromagnetic, or gravitational modes establish decay only for the specified backgrounds, perturbations, and parameter ranges. They do not by themselves prove nonlinear stability or stability in every perturbation sector. Generic Lovelock theories can possess ultraviolet instabilities even though the restricted Lovelock family studied has barrier-like potentials and damped modes (Prasobh et al., 2014). Similarly, weak-charge perturbation theory finds no unstable Kerr–Newman modes in the examined regime but does not settle the full nearly extremal spectrum (Mark et al., 2014).

Modern QNM theory therefore combines boundary-value analysis, scattering theory, microlocal propagation, asymptotic geometry, numerical spectral methods, and representation theory. The common structure is a discrete set of resonant complex frequencies whose real parts describe coherent oscillation and whose imaginary parts describe leakage through horizons, radiation, absorption, or other open channels.