Zero-Damped Modes: Spectral Phenomena

• Zero-damped modes are spectral branches defined by decay rates approaching zero as extremality is reached in various systems, including black-hole perturbations and mechanical vibrations.
• They reveal distinct spectral bifurcations in nearly extremal regimes by contrasting modes that retain finite damping with those whose decay rates vanish.
• Analytical techniques such as matched asymptotic expansions, WKB analysis, and operator theory are employed to characterize these modes across gravitational, mechanical, and kinetic frameworks.

Searching arXiv for recent and foundational papers on zero-damped modes. Zero-damped modes are spectral branches whose damping rate tends to zero in a distinguished limit. In black-hole perturbation theory, they are quasinormal modes with frequencies $\omega=\omega_R+i\omega_I$ such that $\operatorname{Im}\omega\to0$ as extremality is approached, typically with $\omega_R$ approaching a horizon kinematic bound such as $m\Omega_H$ or $m\Omega_H+q\Phi_H$. In other non-self-adjoint settings, the same designation is used for modes whose spectral parameters become purely imaginary or purely real, so that exponential decay disappears. Across the literature, the term therefore refers less to a single mechanism than to a recurrent spectral phenomenon: dissipation vanishes at a geometric, dynamical, or parameter-tuned threshold (Zimmerman et al., 2015, Jovanovic et al., 2011, Burov, 2021).

1. Spectral meaning and basic distinctions

In the quasinormal-mode setting, zero-damped modes are defined relative to the standard radiation boundary conditions: purely ingoing at the event horizon and purely outgoing at infinity. For Kerr–Newman perturbations written in the radial tortoise coordinate $r_*$, the ingoing horizon behavior is $u\sim\Delta^{-s/2}e^{-ik r_*}$ with $k=\omega-m\Omega_H$, while the outgoing asymptotics at infinity are $u\sim r^s e^{-i\omega r_*}$. Zero-damped modes are those for which the decay rate tends to zero in the extremal limit, equivalently $\omega_I\to0$ as the surface gravity $\operatorname{Im}\omega\to0$0 or $\operatorname{Im}\omega\to0$1. They are contrasted with damped modes, whose decay rates remain finite at extremality (Zimmerman et al., 2015).

The same distinction appears in nearly extremal Kerr, where two families coexist: zero-damping modes and damped modes. Zero-damping modes exist for all allowed values of $\operatorname{Im}\omega\to0$2 and $\operatorname{Im}\omega\to0$3, and their frequencies cluster onto the real axis in the extremal limit. Damped modes have nonzero damping for all black-hole spins; they exist for all counterrotating modes $\operatorname{Im}\omega\to0$4 and, in the eikonal limit, for corotating modes with $\operatorname{Im}\omega\to0$5, where $\operatorname{Im}\omega\to0$6 and $\operatorname{Im}\omega\to0$7 (Yang et al., 2012).

Outside black-hole physics, the spectral meaning is analogous but not identical. In the vibrating-bar problem with viscous boundaries and an internal dashpot, zero damping means that all modal exponents satisfy $\operatorname{Im}\omega\to0$8, so the eigenvalues lie on the imaginary axis; equivalently, the algebraic roots $\operatorname{Im}\omega\to0$9 lie on the unit circle $\omega_R$0 (Jovanovic et al., 2011). In longitudinal beam dynamics, zero-damped modes are discrete van Kampen modes with real tunes above the incoherent spectrum, so Landau damping is lost (Burov, 2021). In weakly collisional stellar systems, the phrase refers to true eigenmodes whose damping rate vanishes in the limit $\omega_R$1, with $\omega_R$2 in the normalized variables (Ng et al., 2021).

A common misconception is to equate zero damping with exact conservation or with instability. The literature does not support that identification. In nearly extremal black holes, zero-damped modes can remain quasinormal and decaying for any nonzero $\omega_R$3, with decay rates merely suppressed by $\omega_R$4 or $\omega_R$5 (Zimmerman et al., 2015). In the bar problem, zero damping can coexist with rigid motion through a multiple eigenvalue at $\omega_R$6 rather than through exponential growth (Jovanovic et al., 2011). In the Sen black hole, the existence of zero-damped modes does not imply an instability of a charged massive scalar field under quasinormal boundary conditions (Kokkotas et al., 2015).

2. Nearly extremal black holes: Kerr, Kerr–Newman, and Reissner–Nordström

For nearly extremal Kerr black holes, matched asymptotic expansions yield the standard near-extremal zero-damping formula

$\omega_R$7

with $\omega_R$8 in $\omega_R$9 units. This implies $m\Omega_H$0 and $m\Omega_H$1 as $m\Omega_H$2. In terms of horizon quantities, the leading scaling can be written as

$m\Omega_H$3

The real part approaches the superradiant bound, while the imaginary part is set by the surface-gravity scale (Yang et al., 2012).

For Kerr–Newman black holes, the same qualitative structure persists. With

$m\Omega_H$4

and near-extremality parameter $m\Omega_H$5, the Dudley–Finley analysis gives, for spin-weighted scalar fields,

$m\Omega_H$6

Hence $m\Omega_H$7 and $m\Omega_H$8 up to $m\Omega_H$9 corrections. In the high-$m\Omega_H+q\Phi_H$0 regime, the WKB analysis yields a threshold $m\Omega_H+q\Phi_H$1 separating a ZDM sector, where the potential peak lies at the horizon, from a DM sector, where the peak remains outside. As charge increases and rotation decreases, $m\Omega_H+q\Phi_H$2 increases, shrinking the zero-damped sector (Zimmerman et al., 2015).

The same bifurcation structure survives in more recent Kerr–Newman analyses of gravitational perturbations in the Dudley–Finley approximation. Near extremality, the parameter space splits into subregions containing either only zero-damped modes or zero-damped modes together with finite-damping modes, and analytic boundaries can be expressed through both an eikonal $m\Omega_H+q\Phi_H$3-criterion and a non-eikonal $m\Omega_H+q\Phi_H$4 criterion. For the benchmark gravitational modes $m\Omega_H+q\Phi_H$5, $m\Omega_H+q\Phi_H$6, and $m\Omega_H+q\Phi_H$7, the Dudley–Finley frequencies agree with fully coupled Kerr–Newman calculations within roughly $m\Omega_H+q\Phi_H$8 in $m\Omega_H+q\Phi_H$9 and $r_*$0 in $r_*$1 across the subextremal parameter space, although the errors grow in the high-charge near-extremal regime (Saha et al., 6 Oct 2025).

In the nonrotating limit, nearly extremal Reissner–Nordström black holes provide a distinctive case. The gravito-electromagnetic zero-damped modes are purely damped,

$r_*$2

with parity sectors $r_*$3 and $r_*$4, $r_*$5, where $r_*$6. Thus $r_*$7 while the damping still scales linearly with $r_*$8. This is a useful corrective to the common expectation that zero-damped modes must approach a nonzero oscillation frequency: in RN they are long-lived but purely damped (Zimmerman et al., 2015).

3. Other black-hole settings: dilatonic and de Sitter-type spacetimes

Zero-damped behavior is not restricted to Kerr or Kerr–Newman. For the rotating dilatonic Sen black hole, quasinormal modes of a massive charged scalar field satisfy the horizon-locking condition $r_*$9 near extremality, which gives

$u\sim\Delta^{-s/2}e^{-ik r_*}$0

In the strict extremal limit this becomes

$u\sim\Delta^{-s/2}e^{-ik r_*}$1

while $u\sim\Delta^{-s/2}e^{-ik r_*}$2. Numerically, the spectrum bifurcates into a branch approaching pure real frequencies and a branch with finite damping, and this bifurcation persists even for not small charge $u\sim\Delta^{-s/2}e^{-ik r_*}$3, unlike the Kerr–Newman behavior described in that paper (Kokkotas et al., 2015).

In the large-coupling regime $u\sim\Delta^{-s/2}e^{-ik r_*}$4, the dominant Sen frequencies admit the asymptotic form

$u\sim\Delta^{-s/2}e^{-ik r_*}$5

so that

$u\sim\Delta^{-s/2}e^{-ik r_*}$6

The paper also reports that no instability is found under quasinormal boundary conditions; superradiant instabilities pertain instead to bound-state boundary conditions (Kokkotas et al., 2015).

A more rigorous existence theory has been developed for de Sitter and de Sitter black-hole spacetimes. In the Laplace-transform convention with time dependence $u\sim\Delta^{-s/2}e^{-ik r_*}$7, zero-damped quasinormal frequencies are defined as functions $u\sim\Delta^{-s/2}e^{-ik r_*}$8 such that each $u\sim\Delta^{-s/2}e^{-ik r_*}$9 is a quasinormal frequency, $k=\omega-m\Omega_H$0 as $k=\omega-m\Omega_H$1, and $k=\omega-m\Omega_H$2 as $k=\omega-m\Omega_H$3 for fixed $k=\omega-m\Omega_H$4. For the conformal Klein–Gordon equation on the spherically symmetric backgrounds treated there, the limiting value is $k=\omega-m\Omega_H$5, so the frequencies cluster near $k=\omega-m\Omega_H$6 (Joykutty, 2021).

In pure de Sitter space, the exact spectrum contains the arithmetic progression $k=\omega-m\Omega_H$7, $k=\omega-m\Omega_H$8, so $k=\omega-m\Omega_H$9 as $u\sim r^s e^{-i\omega r_*}$0. The same operator-theoretic framework establishes zero-damped quasinormal frequencies for a class of spherically symmetric black holes with a cosmological horizon and, by a conformal transformation, for near-extremal Reissner–Nordström–de Sitter. The phenomenon is stable under smooth decaying potentials; in de Sitter one has $u\sim r^s e^{-i\omega r_*}$1 for sufficiently small $u\sim r^s e^{-i\omega r_*}$2, and under stronger decay one obtains $u\sim r^s e^{-i\omega r_*}$3 (Joykutty, 2021).

4. Non-self-adjoint mechanical realizations

In the vibrating-bar problem with viscous boundary conditions and an internal damper, zero-damped modes arise from a non-self-adjoint spectral problem in which the spectral parameter $u\sim r^s e^{-i\omega r_*}$4 appears in the boundary conditions. The governing equation is

$u\sim r^s e^{-i\omega r_*}$5

with viscous boundary conditions

$u\sim r^s e^{-i\omega r_*}$6

After Laplace transformation, the eigenvalues are zeros of $u\sim r^s e^{-i\omega r_*}$7, and for rational $u\sim r^s e^{-i\omega r_*}$8 the transcendental characteristic equation reduces to the algebraic equation

$u\sim r^s e^{-i\omega r_*}$9

Zero damping means $\omega_I\to0$0 for all $\omega_I\to0$1, equivalently $\omega_I\to0$2 for all algebraic roots (Jovanovic et al., 2011).

The paper gives explicit zero-damping conditions in the symmetric case $\omega_I\to0$3. Without internal damping, zero damping occurs for

$\omega_I\to0$4

With internal damping, one explicit zero-damped case is

$\omega_I\to0$5

for which all eigenvalues lie on the imaginary axis and there is a double pole at zero. A more general family is obtained from $\omega_I\to0$6, $\omega_I\to0$7, and $\omega_I\to0$8, giving a two-dimensional manifold of zero-damping parameter combinations (Jovanovic et al., 2011).

The physical interpretation in that work is an exact balance between dissipation at the boundaries and energy injection by an active internal dashpot. The rigid-motion condition

$\omega_I\to0$9

makes $\operatorname{Im}\omega\to0$00 a double eigenvalue, so the bar exhibits rigid motion at constant velocity in addition to oscillations. If multiplicity at $\operatorname{Im}\omega\to0$01 is higher, the bar can accelerate. Zero damping is therefore compatible with nontrivial Jordan-structure phenomena at the origin rather than with a purely oscillatory semigroup alone (Jovanovic et al., 2011).

The same study also identifies critical cases where modal expansions fail. When one of $\operatorname{Im}\omega\to0$02 equals $\operatorname{Im}\omega\to0$03, the Green function cannot be represented as a sum of simple partial fractions. Transparent boundaries lead to closed-form finite superpositions of traveling waves, while an amplifying boundary with $\operatorname{Im}\omega\to0$04 makes the Laplace-transformed Green function unbounded at infinity and invalidates the eigenmode expansion. Finite-element calculations in zero-damped regimes with an internal damper can produce spurious eigenvalues with nonzero positive real parts that do not converge away under mesh refinement (Jovanovic et al., 2011).

5. Kinetic and collective systems

In longitudinal beam dynamics with repulsive inductive impedance, the linearized Vlasov equation reduces in the weak-space-charge regime to a Hermitian integral equation. For dipole modes,

$\operatorname{Im}\omega\to0$05

and after rescaling one obtains the parameter-less form

$\operatorname{Im}\omega\to0$06

Because the kernel is real and positive for repulsive inductive impedance, the discrete eigenvalues are real and lie above the incoherent spectrum, $\operatorname{Im}\omega\to0$07. The discrete spectrum is infinite for pure inductive impedance with no roll-off, and its accumulation point is the incoherent zero-amplitude synchrotron frequency. In the paper’s interpretation, these are zero-damped modes: notwithstanding RF bucket nonlinearity and potential well distortion, Landau damping is lost (Burov, 2021).

The physical mechanism is that the coherent force raises the collective frequency while depressing the incoherent particle frequencies. Once the coherent eigenvalue moves above the incoherent band, there are no resonant particles and phase mixing no longer damps the mode. A practical implication is that even a tiny coupled-bunch interaction can destabilize the beam. The paper derives analytic coupled-bunch tune shifts for all multipolarities; for dipole modes the perturbative shift is

$\operatorname{Im}\omega\to0$08

Pure inductive impedance has no intrinsic loss-of-Landau-damping threshold; only high-frequency roll-off or intrabeam scattering truncates the discrete spectrum in practice (Burov, 2021).

Weakly collisional stellar systems provide a closely related but mathematically distinct example. With a Lenard–Bernstein collision operator, the linearized Vlasov–Poisson system has a discrete collisional spectrum satisfying a dispersion relation $\operatorname{Im}\omega\to0$09. As $\operatorname{Im}\omega\to0$10, a subset of these eigenfrequencies converges to the collisionless Landau poles of

$\operatorname{Im}\omega\to0$11

while another subset converges to heavily damped Lenard–Bernstein poles. The continuous Case–van Kampen spectrum is thereby eliminated by arbitrarily weak collisions, and the Landau modes become true eigenmodes in the limit of zero collisions (Ng et al., 2021).

In that framework, zero-damped modes arise in two ways. For a Maxwellian, the least-damped mode approaches zero damping at the Jeans threshold $\operatorname{Im}\omega\to0$12: for $\operatorname{Im}\omega\to0$13 it is stable, at $\operatorname{Im}\omega\to0$14 it is marginal, and for $\operatorname{Im}\omega\to0$15 it becomes the Jeans unstable mode. The paper also shows that tiny non-Maxwellian deviations can create weakly damped branches. For example, a two-component distribution

$\operatorname{Im}\omega\to0$16

with very small tail fraction produces new branches with much smaller damping than the Maxwellian spectrum. This is traced to the modification of the resonant slope $\operatorname{Im}\omega\to0$17 near the relevant phase velocity (Ng et al., 2021).

6. Mechanisms, bifurcation, and open problems

Several distinct mechanisms generate zero-damped behavior. In nearly extremal rotating or charged black holes, the decisive structure is near-horizon trapping combined with vanishing surface gravity. Corotating frequencies approach the superradiant threshold, and the leakage rate scales as $\operatorname{Im}\omega\to0$18 or $\operatorname{Im}\omega\to0$19; in WKB terms, the effective-potential peak moves onto the horizon, while damped modes remain associated with an outer photon-orbit barrier (Yang et al., 2012, Zimmerman et al., 2015). In de Sitter-type spherical problems, the rigorous constructions instead use hyperboloidal slicing, Fredholm theory, and horizon-supported co-modes to prove pole families that converge to the origin in the extremal limit (Joykutty, 2021). In the bar problem, zero damping comes from parameter tuning that places all roots on the unit circle (Jovanovic et al., 2011). In the beam and stellar-system problems, it arises when the discrete spectrum detaches from, or replaces, a continuum that normally mediates Landau damping (Burov, 2021, Ng et al., 2021).

These examples suggest a broad spectral pattern: zero-damped modes appear when the dissipative channel remains present but its effective coupling vanishes in a singular limit, or when the spectral support of resonant absorption is removed. The concrete realization varies sharply by discipline—surface gravity in black holes, active/passive dashpot balance in structural vibrations, positive Hermitian kernels in beam dynamics, and weak-collision discretization in stellar kinetics—but the signature is the same: modal decay moves to the boundary of the complex spectrum rather than disappearing by a change of equation class.

Open questions remain most prominently in the Kerr–Newman problem. For full coupled gravito-electromagnetic perturbations, the precise universal near-extremal frequency function is still unknown; the 2015 Kerr–Newman analysis explicitly pointed to “the existence of a simple universal equation for the frequencies of zero-damped gravito-electromagnetic modes of Kerr-Newman black holes, whose precise form remains an open question” (Zimmerman et al., 2015). More recent work in the Dudley–Finley approximation has sharpened the analytic boundaries between ZDM-only and ZDM+DM regions and clarified the relation to near-horizon and photon-sphere families, but it also emphasizes that high-charge near-extremal behavior remains sensitive to the neglected gravito-electromagnetic coupling (Saha et al., 6 Oct 2025).

A final caution concerns terminology. “Zero-damped” does not impose a universal sign or value of the real part. Kerr, Kerr–Newman, and Sen ZDMs typically satisfy $\operatorname{Im}\omega\to0$20 or $\operatorname{Im}\omega\to0$21 (Yang et al., 2012, Kokkotas et al., 2015). Reissner–Nordström gravito-electromagnetic ZDMs are purely damped with $\operatorname{Im}\omega\to0$22 (Zimmerman et al., 2015). In Laplace conventions for de Sitter-type problems, the same phenomenon is formulated as convergence $\operatorname{Im}\omega\to0$23, with $\operatorname{Im}\omega\to0$24 in the spherical cases proved there (Joykutty, 2021). The term therefore denotes a limiting spectral behavior, not a single canonical dispersion law.