QNMs in Kerr-Newman-de Sitter Black Holes

• The paper demonstrates that QNMs in KNdS spacetimes form a discrete complex spectrum with the real part indicating oscillatory frequencies and the imaginary part quantifying decay rates.
• Methodologies include Heun polynomial quantization, WKB approximations, and Frobenius methods to extract universal spectra independent of the spin of massless fields.
• Insights into regimes like eikonal, near-extremal, and slowly rotating limits underscore the role of QNMs in stability analysis and gravitational-wave parameter estimation.

The quasinormal modes (QNMs) of Kerr-Newman-de Sitter (KNdS) black holes characterize the dissipative oscillations of linear perturbations in spacetime geometries featuring non-zero mass, rotation ($a$), electric charge ($Q$), and a positive cosmological constant ($\Lambda$). As solutions to the separated field equations with specific boundary conditions—purely ingoing at the outer event horizon and purely outgoing at the cosmological horizon—QNMs form a discrete and typically complex spectrum, with the imaginary part quantifying decay rates and the real part encoding oscillation frequencies. Their structure underpins gravitational wave emission, stability analysis, and theoretical issues such as the strong cosmic censorship conjecture.

1. Field Equations, Metric, and Separation of Variables

The four-dimensional Kerr-Newman-de Sitter metric in Boyer–Lindquist coordinates

$$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$

with

$$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$

admits separability of the scalar, spinor, electromagnetic, and gravitational perturbation equations (Li, 25 Dec 2025, Hod, 2018, Churilova et al., 2021). In the context of massless fields of spin $s$, the master equation reduces to angular and radial ODEs. The radial equation for spin-$s$ fields can be cast in general as a Teukolsky-type equation, with regular singularities at the (Cauchy, event, and cosmological) horizons and at infinity (Hortacsu, 2020). Boundary conditions for QNMs are

• $$\psi(r) \sim (r-r_+)^{-\frac{i(\omega - m\Omega_+)}{2\kappa_+}}$$ (ingoing at $r_+$)
• $$\psi(r) \sim (r_c-r)^{-\frac{i(\omega - m\Omega_c)}{2\kappa_c}}$$ (outgoing at $Q$0), where $Q$1 and $Q$2 are the surface gravities and angular velocities at the event and cosmological horizons, respectively.

2. Analytic Spectra and Degeneracy in the General KNdS Geometry

For generic parameters, the quantization of QNMs is linked to the existence of polynomial solutions to the radial Teukolsky equation, which is of Heun type with four regular finite singular points corresponding to the horizons. Imposing truncation of the local Heun series at order $Q$3 (i.e., $Q$4 for the series parameters) discretizes the spectrum (Li, 25 Dec 2025, Hortacsu, 2020). The closed-form expression for the QNM frequencies of all massless spin fields ($Q$5) is

$Q$6

where $Q$7 and $Q$8 ($Q$9) are the surface gravities and angular velocities at the Cauchy and event horizons. This spectrum is universal (i.e., independent of the field spin) for fixed $\Lambda$0, and the radial wavefunction features an internal label $\Lambda$1, indicating an $\Lambda$2-fold degeneracy; different spin perturbations "mimic" one another at the level of observed frequencies (Li, 25 Dec 2025).

3. Limiting Regimes: Eikonal, Near-Extremal, and Slowly Rotating Approximations

Several analytically tractable limits are critical:

• Eikonal (large-$\Lambda$3) regime: The QNMs correspond to the properties of unstable null circular orbits. The frequencies take the asymptotic form $\Lambda$4 where $\Lambda$5 and $\Lambda$6 are the angular velocity and surface gravity at the relevant photon sphere/event horizon. The imaginary part, set by $\Lambda$7, controls the exponential decay rate and is pivotal for stability and strong cosmic censorship (Hod, 2018, Churilova et al., 2021).
• Near-extremal KNdS: As $\Lambda$8, the surface gravity $\Lambda$9, and near-horizon (NH) modes display ultraslow decay, dominating the late ringdown spectrum. In this regime, the QNMs admit a closed analytic form for both bosons and fermions:

$$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$0

where the $$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$1 coefficients and angular eigenvalue $$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$2 generalize the separation constant for the rotating geometry (Churilova et al., 2021).

• Slowly rotating and slowly charged limit: For $$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$3, QNM frequencies split linearly in the azimuthal quantum number $$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$4 with a Zeeman-like effect, lifting the degeneracy of $$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$5. To leading order in semiclassical parameter $$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$6 (large $$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$7),

$$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$8

with $$ds^2 = \frac{\rho\bar\rho}{\Delta_r}\,\Xi^2\bigl(dt - a\sin^2\theta\,d\varphi\bigr)^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,\sin^2\theta \frac{\bigl[a\,dt - (r^2+a^2)\,d\varphi\bigr]^2}{\Xi^2} -\frac{\rho\bar\rho}{\Delta_r}\,dr^2 -\frac{\rho\bar\rho}{\Delta_\theta}\,d\theta^2$$9 and $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$0 determined from the photon sphere geometry (Iantchenko, 2015).

4. Families of Quasinormal Modes and Spectral Hierarchy

A systematic numerical and analytic study identifies three principal families of QNMs for KNdS black holes (Davey et al., 2024):

• de Sitter (dS) modes: Dominant for small black holes or in the $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$1 limit, with frequencies $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$2.
• Photon-sphere (PS) modes: Govern intermediate decay and are associated with wave trapping near the unstable null orbit; asymptotic form $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$3, with $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$4 and $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$5 the orbital frequency and Lyapunov exponent of the photon sphere.
• Near-horizon (NH) modes: Dominant very close to extremality, with slowest decay, $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$6, $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$7.

As parameters are tuned towards extremality, "eigenvalue repulsion" occurs: the imaginary parts of PS and NH modes approach, exchange dominance, and govern the late-time behavior.

5. Exact and Approximate Solution Techniques

The core analytic technique for deducing the QNM spectra is reduction to a Heun equation, truncated by boundary conditions and polynomiality. For specific parameter values or regimes, three primary computational approaches are employed (Hortacsu, 2020, Churilova et al., 2021):

• Heun polynomial quantization: QNM frequencies correspond to truncation conditions of the local Heun series in the radial coordinate, supported by a vanishing determinant that implicitly depends on $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$8.
• WKB and Pöschl-Teller fitting: For large $$\begin{aligned} &\rho = r - ia\cos\theta,\quad \Delta_r = (r^2 + a^2)\left(1 - \tfrac{\Lambda}{3} r^2\right) - 2Mr + Q^2, \ &\Delta_\theta = 1 + \tfrac{\Lambda}{3}a^2 \cos^2\theta,\quad \Xi = 1 + \tfrac{\Lambda}{3}a^2, \end{aligned}$$9 (eikonal), the potential is fit near its maximum to a Pöschl-Teller form, giving approximate analytic QNMs.
• Frobenius/continued-fraction methods: Particularly in the near-extremal regime, the radial differential equation reduces to two-term recursions, enabling analytic extraction of the QNMs for both integer and half-integer spin.

Sample computations, both analytic and numerical, confirm the general statements and support mode stability: all computed $s$0, indicating exponential decay in time (Media et al., 21 Feb 2025, Hortacsu, 2020).

6. Parameter Dependence, Degeneracies, and Physical Implications

The QNM frequencies for KNdS black holes depend algebraically on the physical parameters $s$1 and discrete quantum numbers $s$2, but strikingly not on the spin $s$3 for massless fields; this universal behavior extends to all $s$4 fields (Li, 25 Dec 2025, Churilova et al., 2021). Consequently, diverse perturbing fields (scalar, electromagnetic, Dirac, gravitational) can "mimic" each other's QNM response (frequency degeneracy), though the corresponding radial functions possess extra degeneracy labels (degree $s$5 for overtone $s$6). Observationally, discrimination between field types based solely on the ringdown frequency spectrum is impossible unless further mode structure (e.g., radial profile) is resolved.

In the slowly rotating, bumblebee gravity case (Lorentz-violation parameter $s$7), both $s$8 and $s$9 decrease with increasing $s$0; the potential barrier is suppressed and QNM decay slows (Media et al., 21 Feb 2025). Variation of $s$1, $s$2 and $s$3 systematically shifts the QNM positions, as confirmed by explicit tabular data.

7. Stability, Strong Cosmic Censorship, and Observational Aspects

All calculations and numerical searches indicate mode stability: no QNM with $s$4 is found. In the near-extremal regime, the least-damped NH modes become arbitrarily long-lived, but always decay. Analytically, in the eikonal limit, the decay rate is set by the surface gravity $s$5, ensuring that for neutral massless fields, the fundamental QNM satisfies $s$6, thereby upholding Penrose's strong cosmic censorship conjecture (Hod, 2018, Davey et al., 2024).

A plausible implication is that gravitational-wave detectors can utilize closed-form expressions for QNM frequencies to efficiently extract black hole parameters $s$7 via Bayesian inference, with universality and degeneracy tests constraining possible deviations from standard general relativity. Laboratory analogues can exploit this mimicking to simulate gravitational ringdown with electromagnetic or acoustic systems (Li, 25 Dec 2025).

---

References:

(Hod, 2018, Li, 25 Dec 2025, Churilova et al., 2021, Hortacsu, 2020, Iantchenko, 2015, Davey et al., 2024, Media et al., 21 Feb 2025)