Parametrized Black Hole QNMs

• Parametrized black hole QNMs are a systematic framework that expresses small deviations from GR as linear shifts in the frequency spectra of perturbed black holes.
• The approach employs numerical methods, including Leaver’s continued-fraction technique and spectral collocation, to compute corrections in both static and rotating spacetimes.
• Applications to modified gravity and quantum effects yield spectral fingerprints that are crucial for gravitational-wave data analysis and testing alternative theories.

Parametrized black hole quasinormal modes (QNMs) constitute a systematic framework to describe, predict, and constrain the linear oscillation spectra of perturbed black hole spacetimes in the presence of small deviations from general relativity (GR). This formalism translates arbitrary small deformations of the background metric or the perturbation equations—arising from modified gravity, quantum effects, or phenomenological modifications—into explicit shifts of the QNM frequencies, controlled by a set of (in principle) theory-independent parameters. The resulting spectral fingerprints serve as precision probes of strong-field gravity, allowing one to test GR and constrain alternative theories through gravitational-wave observations of ringdown.

1. Formalism for Parametrized QNMs in Static Spacetimes

Parametrized QNM theory for spherically symmetric, nonrotating black holes proceeds by expressing the perturbation equations in a standard Regge–Wheeler (axial) or Zerilli (polar) form,

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

where $dr_*/dr = 1/f(r)$, with $f(r) = 1 - r_H/r$ for Schwarzschild background and $s \in \{\text{axial},\text{polar}\}$. Small deviations from GR are incorporated as corrections to the effective potential,

$$V^{(s)}_\ell(r) = V^{(s),\text{GR}}_\ell(r) + \delta V_s(r)\,,$$

where

$$\delta V_s(r) = \frac{1}{r_H^2} \sum_{j=0}^\infty \alpha_j^{(s)} \left(\frac{r_H}{r}\right)^j.$$

Here, the dimensionless coefficients $\alpha_j^{(s)} \ll 1$ serve as theory-agnostic deformation parameters. The QNM frequencies then admit a linear expansion,

$$\omega_{\ell n}^{(s)} = \omega_{\ell n}^{(s),\text{GR}} + \sum_{j=0}^N \alpha_j^{(s)} e_{\ell n;j}^{(s)}\,,$$

with $e_{\ell n;j}^{(s)}$ numerically determined “basis vectors,” independent of the specific theory once the background is specified. This methodology is codified and implemented using precomputed data tables for $e_j$ (Cardoso et al., 2019).

Parametric deformations generically break the isospectrality between axial and polar modes unless the hierarchy

$dr_*/dr = 1/f(r)$0

is satisfied to all orders.

2. Extensions: Coupled Systems, Quadratic Corrections, and Degeneracies

The formalism generalizes to arbitrary finite sets of coupled wave equations,

$dr_*/dr = 1/f(r)$1

where $dr_*/dr = 1/f(r)$2 is a vector of N master variables, and $dr_*/dr = 1/f(r)$3 is a real symmetric matrix. Generic deviations are expressed via power series in $dr_*/dr = 1/f(r)$4: $dr_*/dr = 1/f(r)$5 with linear and quadratic correction coefficients tabulated by mode. For nondegenerate spectra, the first- and second-order shifts of the QNM frequencies are computed via derivatives of the determinant of a scattering matrix: $dr_*/dr = 1/f(r)$6 with explicit numerical tables for $dr_*/dr = 1/f(r)$7 and $dr_*/dr = 1/f(r)$8 (McManus et al., 2019). For degenerate unperturbed spectra, such as the axial–polar pair in GR, a secular equation must be diagonalized to resolve the linear level splitting.

3. Parametrized QNMs for Rotating Black Holes

The parametrized QNM program extends to spinning (Kerr) black holes by considering small linear deformations of the Teukolsky equation,

$dr_*/dr = 1/f(r)$9

where $f(r) = 1 - r_H/r$0 denotes a monomial basis in $f(r) = 1 - r_H/r$1. The QNM frequency and angular separation constant are expanded as

$f(r) = 1 - r_H/r$2

with $f(r) = 1 - r_H/r$3 (and $f(r) = 1 - r_H/r$4) computed from linear perturbation theory on Leaver’s continued-fraction solution (Cano et al., 2024).

For arbitrary small ($f(r) = 1 - r_H/r$5) modifications, this yields accurate linear corrections. Examples include massive scalar fields, (proxy) Kerr–Newman perturbations, and higher-derivative gravity, with tabulated coefficients available for multiple (n, $f(r) = 1 - r_H/r$6, m) triplets and spin parameters.

4. Methodologies and Computational Implementation

Parametrized QNM computation requires the following elements:

1. Specify the background metric and the form of linearized deviations in terms of the chosen parametric basis.
2. Numerically obtain the unperturbed QNM spectrum ($f(r) = 1 - r_H/r$7) and mode functions (via Leaver's method, direct integration, or spectral collocation).
3. Introduce each deformation monomial, solve the perturbed master equation, and extract the linear shift ($f(r) = 1 - r_H/r$8 or equivalent) by root-finding. Recursive relations are available for generating higher-order coefficients from a minimal set of basis terms (Kimura, 2020).
4. For multi-channel or coupled systems, the determinant of the constructed scattering matrix encodes the spectral condition.
5. Asymptotic formulas and connection with analytic techniques (eikonal, WKB, Leaver) underpin the high-overtone regime and inform the convergence of the linear expansion (Miyachi et al., 21 Dec 2025).

The process is fully algorithmic, with basis coefficient tables provided in public repositories for rapid QNM evaluation in arbitrary small-deviation models (Cardoso et al., 2019, McManus et al., 2019).

5. Representative Applications to Modified Gravity and Quantum Corrections

Parametrized QNM theory has been exploited for multiple modified- and quantum-gravity models:

• Loop Quantum Gravity: QNMs of self-dual LQG black holes (Santos et al., 2015) and holonomy-corrected Schwarzschild backgrounds (Moreira et al., 2023) depend controllably on polymeric and area-gap parameters, manifest as shifted frequencies and dampings. The WKB and continued-fraction methods benchmark and cross-validate the parametric expansions.
• Generalized Uncertainty Principle and Quintessence: GUP-modified Schwarzschild black holes surrounded by quintessence fluids feature deformation parameter dependencies that can be captured semi-analytically with Mashhoon/Pöschl–Teller and high-order WKB expansions, with observable trends in both real and imaginary QNM parts (Karmakar et al., 2022).
• Spherically Symmetric Metrics—Continued-Fraction Formalism: General asymptotically flat spherically symmetric metrics admit a continued-fraction expansion, with QNMs well approximated by the lowest-order parameter coefficients, supporting a phenomenologically hierarchical approach to gravitational ringdown constraints (Konoplya et al., 2022).
• Quadrupole Moment and Multipole Deformations: In the eikonal regime, QNM frequencies for static or rotating spacetimes with arbitrary quadrupole and higher mass moments can be analytically linked to the light-ring location and Lyapunov exponent, yielding explicit dependence on the quadrupole parameter (Allahyari et al., 2018).
• Near-Horizon Universality: For charged scalar fields in general spacetimes, the near-horizon QNM spectrum—a “tower” of modes with frequencies $f(r) = 1 - r_H/r$9—demonstrates the universal dependence on horizon radii and surface gravity, independent of further asymptotic structure (Tangphati et al., 2018).

6. Asymptotics, High-overtone Instabilities, and Observational Implications

In classical GR, the QNM spectrum for black holes displays characteristic convergence of $s \in \{\text{axial},\text{polar}\}$0 for large overtone index $s \in \{\text{axial},\text{polar}\}$1, particularly at high damping. Parametric deviations introduce new analytic structures:

• For deformations entering at order $s \in \{\text{axial},\text{polar}\}$2, $s \in \{\text{axial},\text{polar}\}$3, the real part $s \in \{\text{axial},\text{polar}\}$4 typically diverges as $s \in \{\text{axial},\text{polar}\}$5, either logarithmically (for tuned cubic deformations) or as a power law ($s \in \{\text{axial},\text{polar}\}$6) (Miyachi et al., 21 Dec 2025, Moreira et al., 2023).
• This breakdown of the “GR plateau” in $s \in \{\text{axial},\text{polar}\}$7 is robust evidence for strong-field beyond-GR signatures, distinct from corrections to low-lying (fundamental or first overtone) modes.
• The existence of near-GR “quasi-bound” states and algebraically special (noncontinuous) modes in perturbed potentials highlights the need to scan the full complex-frequency plane when interpreting deviations in observational data (Cardoso et al., 2019).
• For moderate background changes (radius, redshift, higher multipole coefficients), QNMs are well constrained by only a few leading parametric coefficients. For nonmoderate geometries (rapid metric variations or strong echoing), higher-order coefficients become important, but current observational data already rule out drastic departures from the GR ringdown template (Konoplya et al., 2022).

7. Impact on Gravitational Wave Data Analysis and Constraints

Parametrized QNM frameworks have been directly incorporated into gravitational-wave ringdown pipelines, enabling simultaneous parameter estimation of mass, spin, and beyond-GR coefficients. The “degeneracy web” between quantum/coupling parameters and intrinsic black hole parameters can significantly affect posterior distributions, quantifying systematic uncertainties due to the underlying theory (Chen et al., 31 Oct 2025). Constraints derived from observed frequencies and damping rates of dominant and overtone modes permit model-agnostic exclusion of large classes of non-GR or quantum-corrected solutions.

Parametrized QNM theory thus provides a unified, quantitative, and extensible toolbox for probing the nature of compact-object spacetimes with precision gravitational-wave spectroscopy, connecting analytic, numerical, and observational facets of fundamental gravity research (Cardoso et al., 2019, McManus et al., 2019, Cano et al., 2024, Miyachi et al., 21 Dec 2025, Konoplya et al., 2022).