Pseudospectra of Quasinormal Frequencies

• Pseudospectra of Quasinormal Frequencies is a framework that defines and quantifies the sensitivity of black hole quasinormal modes using non-self-adjoint operator theory.
• It utilizes hyperboloidal compactification and high-order Chebyshev discretization to map spectral bulges and predict mode stability or transient instabilities.
• The approach enhances gravitational-wave spectroscopy by uncovering non-modal phenomena and overtone instabilities in both black hole and exotic compact object systems.

Quasinormal frequencies (QNFs) are the discrete complex eigenvalues associated with wave propagation subject to outgoing boundary conditions in black hole or exotic compact object spacetimes. The pseudospectrum of quasinormal frequencies provides a rigorous, operator-theoretic framework for quantifying the sensitivity (stability or instability) of these QNFs to small perturbations—either of the underlying spacetime, of the effective potential, or of the differential operator itself. The notion of the pseudospectrum, rooted in non-self-adjoint spectral theory, exposes non-modal and transient phenomena invisible to standard modal analysis and is now a cornerstone in precision gravitational-wave spectroscopy and the study of both black hole and analog systems.

1. Mathematical Definition and Core Properties

Let $A$ denote a linear, generally non-self-adjoint wave operator whose quasinormal eigenproblem takes the schematic form $(A - \omega^2 I)\psi = 0$ with physically motivated (e.g., outgoing) boundary conditions. The $\epsilon$-pseudospectrum $\Sigma_\epsilon(A)$ is: $$\Sigma_\epsilon(A) = \left\{ \omega \in \mathbb{C} : \|(A - \omega^2 I)^{-1}\| > 1/\epsilon \right\}$$ Equivalently, $\omega \in \Sigma_\epsilon(A)$ if there exists a perturbation $\Delta A$ with $\|\Delta A\| < \epsilon$ such that $\omega^2$ is an eigenvalue of $A + \Delta A$. For eigenproblems in hyperboloidal coordinates and after first-order reduction, one often works with a shifted operator $L$ and spectral parameter $s = -i r_0 \omega$, so that: $$\Sigma_\epsilon(L) = \left\{s : \|(sI - L)^{-1}\| > 1/\epsilon \right\}$$ For non-normal operators encountered in QNM problems, $\Sigma_\epsilon(A)$ typically forms "bulges" or extended regions in the complex plane that can lie far from the true spectrum, signifying a high degree of spectral instability under small perturbations (Siqueira et al., 23 Jan 2025, Jaramillo et al., 2020, Destounis et al., 2023, Cownden et al., 2023).

2. Numerical Computation and Hyperboloidal Approach

The practical computation of QNM pseudospectra proceeds as follows:

• Hyperboloidal compactification: Coordinates such as $t = r_0[\tau - H(\sigma)], \; r = r_0/\sigma$ regularize both the event horizon and future null infinity in a compact domain $\sigma \in [0, 1]$.
• First-order reduction: The wave equation is cast as $L u = s u$ or, for 2D problems (e.g., spinning black holes), $L u = \omega u$ with $L$ a non-self-adjoint differential operator. All outgoing boundary conditions are built into the operator's singular structure and regular solution requirement.
• Spectral discretization: Chebyshev–Lobatto grids are used for high-order accuracy. Differentiation and metric-dependent matrices are assembled, yielding a discretized finite-rank operator $\hat{L}$.
• Energy norm: The operator norm is induced by norms physically motivated by conserved stress-energy or energy flux, not naive $L^2$. For instance,

$$\|u\|^2_E = \tfrac{1}{2}\int_{0}^{1}\left[ w|\partial_\tau \bar{\psi}|^2 + p|\partial_\sigma \bar{\psi}|^2 + q|\bar{\psi}|^2 \right] d\sigma$$

• Pseudospectrum mapping: The smallest generalized singular value $\sigma_\text{min}(sI - \hat{L})$ is computed with respect to the discrete energy norm Gramian. Contours where $\sigma_\text{min} = \epsilon$ are plotted in the complex frequency plane; the level sets trace the $\epsilon$-pseudospectrum (Siqueira et al., 23 Jan 2025, Cai et al., 5 Jan 2025).

3. Physical Interpretation and Spectral Stability

The structure of pseudospectra directly encodes the spectral (in)stability of QNMs:

• Stable spectrum: Concentric, tightly closed pseudospectral contours encircling each QNM indicate that small perturbations ($\epsilon$) produce only $O(\epsilon)$ eigenvalue shifts.
• Unstable spectrum: Extended, open, or merging contours—often seen for overtones—reveal that perturbations of size $\epsilon$ can cause $\gg O(\epsilon)$ QNM displacements.
• Universal features: Overtones are generally much more unstable than the fundamental mode, as revealed by large condition numbers. The fundamental mode typically remains robust up to $O(1)$ perturbations unless special long-range modifications or near-degeneracies arise (Destounis et al., 2023, Destounis et al., 2023, Jaramillo et al., 2020, Destounis et al., 2021).
• Pseudo-resonances and non-modal growth: In some regimes, pseudospectral contours cross into regions with $\operatorname{Im}\omega > 0$ (unstable growth), indicating susceptibility to nonmodal pseudo-resonances and enhanced transient responses, especially in horizonless objects or AdS black brane contexts (Boyanov et al., 2022, Cownden et al., 2023, Arean et al., 2023).

4. Key Findings from Schwarzschild-like, Kerr, AdS, and Exotic Compact Object Cases

A summary of empirical and semi-analytic results, with focus on physical scenarios:

Scenario	Pseudospectral Feature	Physical Consequence
Schwarzschild, random $\delta V$	Tight contours (stable)	QNMs displaced $O(\epsilon)$
Schwarzschild-like (RZ deformation)	Wildly open contours	Overtone QNMs become unstable, large shifts for $n\gtrsim6$ (Siqueira et al., 23 Jan 2025)
Ad-hoc high-$k$ potential deformations	Irregular, asymmetric bulges	Amplifies overtone instability
Kerr black hole (spin $s=0$)	Highly non-normal $L$, slow convergence for large $|\operatorname{Im}\omega|$	Overtones exhibit instability, depends on norm; higher Sobolev norms tame instability only partially (Cai et al., 5 Jan 2025)
Rotating analog black holes (DBT)	Prograde overtone stabilization at high spin	Prograde overtones inherit stability of fundamental, potential for multi-mode gravitational-wave fits (Paula et al., 31 Mar 2025)
de Sitter black holes ($\ell=0$ zero mode)	Subtle norm singularity	Instability under perturbations recovered when norm correctly handled (Destounis et al., 2023)
Horizonless ECOs	Pseudospectral contours cross real axis	Nonmodal growth, bootstrap instabilities in nonlinear perturbation (Boyanov et al., 2022)

5. Applications and Implications for Gravitational-Wave Astronomy

The pseudospectrum framework modifies core expectations in black hole spectroscopy and ringdown analysis:

• QNM overtone extraction: Overtones are highly sensitive to high-frequency (UV) and/or small-scale ("dirty") perturbations; standard modal fitting may be unreliable without pseudospectral diagnostics.
• Inverse problem complexity: When multiple perturbation sources are active—e.g., spacetime deformation (RZ-parametrized), high-$k$ environmental effects—overtone spectra from distinct physical origins can overlap, hampering identification of new physics solely from QNM spectroscopy (Siqueira et al., 23 Jan 2025).
• Fundamental mode stability: Time-domain ringdown dominated by the fundamental QNM remains robust to $O(\epsilon)$ in a wide range of backgrounds; this underpins the accuracy of mass and spin inference in GW signals under most physical circumstances (Destounis et al., 2023, Destounis et al., 2023).
• Transient and non-modal dynamics: Pseudospectra predict possible large transient responses or pseudo-resonant growth even when modal analysis suggests decay, especially in exotic settings (ECOs, AdS/CFT).
• Holographic context: In strongly coupled duals, QNM pole positions (transport coefficients) can shift dramatically under model perturbations unless hydrodynamic modes are isolated; near-exceptional points (e.g., pole collisions) show marked pseudospectral instability (Garcia-Fariña et al., 2 Feb 2026, Cownden et al., 2023, Arean et al., 2023).

6. Analytic Approaches and Asymptotic Results

• First-order perturbation theory: Exact Rayleigh–Schrödinger-type formulas exist in special cases, e.g., de Sitter static patch with known modes and co-modes, with condition numbers for high-$n$ overtones growing exponentially in $n$ (Warnick, 2024).
• Asymptotic scaling laws: The outer boundary of the pseudospectrum is often logarithmic: $$\operatorname{Im}\omega \sim C_1 + C_2\log(\operatorname{Re}\omega + C_3)$$ for asymptotically flat, polynomial for AdS cases (Destounis et al., 2021, Cownden et al., 2023).
• Universality and gauge independence: Across a wide class of scenarios (Schwarzschild, RNdS, Kerr, AdS), the global structure of pseudospectral bulges and instability regions proves insensitive to coordinate gauge, norm details, or precise operator block structure—establishing the geometric nature of the pseudospectral phenomenon (Destounis et al., 2021, Cownden et al., 2023).

7. Broader Context and Future Directions

Pseudospectra unify the understanding of modal instability, environmental migration, and exotic dynamical growth in a nonperturbative, operator-theoretic manner. Their computation is now standard in the numerical relativity and gravitational-wave modeling toolkit, both as a practical diagnostic and as a theoretical benchmark for the prospects and limitations of precision ringdown measurements in diverse compact object scenarios. Future work will extend these analyses to full Kerr (spin-weighted, mode-coupled), nonlinear and time-dependent perturbations, and coupled multi-field systems, integrating pseudospectral methods directly into waveform modeling and data analysis pipelines for current and next-generation detectors (Siqueira et al., 23 Jan 2025, Paula et al., 31 Mar 2025, Cai et al., 5 Jan 2025, Destounis et al., 2023).