⚛️ general relativity

The quasinormal modes of Kerr-Newman black holes: a separated approach

This paper demonstrates that a recently discovered exact separation of gravitoelectric perturbation equations provides a highly accurate and practical method for computing Kerr-Newman black hole quasinormal modes, yielding results consistent with coupled-system calculations and offering numerical evidence for spin-label isospectrality even in near-extremal configurations.

Original authors: Kyriakos Destounis, Vitor Cardoso, Peter Hintz

Published 2026-10-07

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Kyriakos Destounis, Vitor Cardoso, Peter Hintz

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Black holes are often imagined as silent, static pits in space, but when they are disturbed, they ring like a bell. If a star crashes into a black hole or if two black holes collide, the resulting object does not settle down instantly. Instead, it vibrates, sending out ripples in the fabric of space and time known as gravitational waves. These vibrations have a specific pitch and a specific duration before they fade away. In physics, these characteristic tones are called quasinormal modes. By listening to these tones, scientists can learn the mass and spin of the black hole, much like a musician can identify a violin by the sound of its strings. However, calculating these tones for a black hole that is both spinning and electrically charged is an incredibly difficult mathematical puzzle. For decades, the equations describing such a black hole were so tangled that they could not be solved directly, forcing researchers to rely on approximations or extremely slow computer simulations.

A team of researchers has now found a way to untangle these equations completely. They focused on the Kerr–Newman black hole, a theoretical object that possesses mass, spin, and electric charge. While real black holes in the universe are likely to have negligible electric charge, this theoretical model serves as a crucial testing ground for the laws of gravity and electromagnetism. The researchers discovered a new method to separate the complex equations that describe the black hole's vibrations. Previously, the equations for the gravitational waves and the electromagnetic waves were locked together, making them nearly impossible to solve. The new approach splits the problem into two distinct parts: one that describes how the vibrations behave around the black hole's equator and poles, and another that describes how they move inward toward the event horizon and outward into space. These two parts are linked, but they can be solved independently, which transforms a nearly impossible task into a manageable one.

Using this separated approach, the team built a computer code to calculate the exact frequencies of these vibrations. They tested their method by comparing their results with other high-precision calculations that had solved the original, tangled equations. The match was astonishingly close. For the most fundamental vibrations, the difference between their new method and the old, coupled methods was so small that it was only visible in the tenth decimal place or beyond. Even for more complex, higher-frequency vibrations, the agreement remained within a range of one part in a million. This level of precision confirms that the new separated equations accurately describe the physics of the black hole. It validates the method as a reliable tool for exploring the "ringdown" phase of rotating, charged black holes without the need for the heavy computational cost of solving the full, coupled system.

The researchers also used this new tool to investigate a subtle property of these vibrations called isospectrality. In simpler terms, they asked whether the black hole would ring at the same pitch regardless of how the mathematical description of the vibration was framed. They found that for the fundamental tones, the answer is yes. Whether they described the vibration using one mathematical perspective or its opposite, the resulting frequency was identical to within a tiny fraction of a percent. This numerical evidence suggests a deep symmetry in the way charged, spinning black holes vibrate, even when they are spinning at nearly the maximum possible speed. While this does not yet constitute a formal mathematical proof, it provides strong support for the idea that these symmetries exist in nature.

This work opens the door to a more detailed exploration of black hole physics. Because the new method is efficient and precise, it allows scientists to map out the behavior of black holes across a wide range of spins and charges, including those that are nearly extreme. It also provides a clearer path to understanding how these objects respond to external forces, such as the tidal pull of a nearby star. By offering a practical and accurate way to calculate the ringdown of charged black holes, the study gives astronomers and theorists a sharper instrument to decode the signals from the most violent events in the universe. The ability to separate the equations means that the complex interplay between gravity and electricity can now be studied with a clarity that was previously out of reach, turning a theoretical laboratory into a precise instrument for discovery.

Technical Summary: The Quasinormal Modes of Kerr–Newman Black Holes: A Separated Approach

Problem Statement
The paper addresses the computation of quasinormal modes (QNMs) for Kerr–Newman (KN) black holes, which describe the characteristic damped oscillations of perturbed spacetimes in Einstein–Maxwell theory. Unlike simpler black hole backgrounds, KN perturbations involve coupled gravitational and electromagnetic degrees of freedom. This coupling renders the spectral problem significantly more challenging than in the Kerr or Reissner–Nordström limits. Previous approaches have relied on approximations (e.g., Dudley–Finley, slow-rotation, or small-charge limits) or fully numerical methods solving the coupled system of equations. The authors aim to investigate the applicability of a recently discovered exact separation method [16] to the KN QNM problem, seeking a tractable description that retains the interaction between gravitational and electromagnetic perturbations while reducing the problem to ordinary differential equations (ODEs).

Methodology
The authors employ a recently derived exact separation of the gravitoelectric perturbation equations. This formulation decouples the radial and angular dependencies into a joint spectral problem for the complex frequency $\omega$ and a separation constant $\lambda$ , linked through two distinct ODEs:

Angular Problem: A matrix angular ODE for complex-valued functions $\chi_1, \chi_2$ depending on the polar angle $\theta$ . The system involves a spin label $\varsigma = \pm 1$ and a separation constant $\lambda$ . The problem is solved by imposing regularity at the angular poles ( $\theta = 0, \pi$ ).
Radial Problem: A scalar second-order radial ODE for the variable $w(r)$ . The boundary conditions require ingoing behavior at the event horizon (after factoring out the physical phase) and outgoing behavior at infinity (defined via analytic continuation rather than simple decay on the real axis).

The numerical implementation combines spectral techniques:

Angular: Galerkin projection.
Radial: Chebyshev collocation.
The code solves the system simultaneously for $\omega$ and $\lambda$ by finding the roots of functional relations $F_{ang}(\omega, \lambda) = 0$ and $F_{rad}(\omega, \lambda) = 0$ . The authors utilize dimensionless variables scaled by the black hole mass $M$ .

Key Contributions and Results

Numerical Benchmarking: The authors benchmark their separated numerical code against independent calculations of the fully coupled gravitoelectric system [13, 15, 17]. The comparison yields absolute differences ( $\delta\omega_n$ ) ranging from $10^{-10}$ to $10^{-5}$ for fundamental modes and first overtones, depending on the black hole parameters and overtone number. This agreement validates the separated formulation as a practical route for computing KN QNMs.
Limiting Cases: The code is tested against known limits:
- Schwarzschild ( $a, Q \to 0$ ): Results agree with matrix solver data to at least 6 digits. The authors identify that the branch label $b$ distinguishes between the gravitational sector ( $b=+1$ , spin $s=\pm 2$ ) and the electromagnetic sector ( $b=-1$ , spin $s=\pm 1$ ).
- Reissner–Nordström ( $a \to 0$ ): Results match published data [21] to all displayed decimal points, confirming the branch selection logic.
- Kerr ( $Q \to 0$ ): Comparisons with Kerr QNM literature [22–24] show agreement around 6 decimal points. The authors note that very small, non-zero charges ( $10^{-4} \lesssim Q/M \lesssim 10^{-2}$ ) are used to avoid singularities in the equations at exactly $Q=0$ .
Spin-Label Isospectrality: The study investigates whether the two spin-label descriptions ( $\varsigma = \pm 1$ ) yield identical fundamental frequencies. The authors find numerical evidence of isospectrality in both the electromagnetic-led and gravitational-led sectors, even for near-extremal configurations. The absolute discrepancies between the $\varsigma = +1$ and $\varsigma = -1$ calculations are extremely small, ranging from $10^{-12}$ to a few $10^{-10}$ . This suggests a consistency in the separated formulation, though the authors note this is numerical evidence rather than a general analytic proof.

Significance and Claims
The paper claims that the separated equations successfully reproduce Kerr–Newman quasinormal frequencies, offering a practical and efficient method for exploring the ringdown of rotating, charged black holes. The formulation allows for independent control of angular and radial resolutions and facilitates the systematic exploration of the black hole parameter space.

The authors emphasize that while the results provide strong numerical evidence for spin-label isospectrality and validate the separated approach against coupled systems, they do not constitute a general analytic proof of parity-isospectrality for all spins. They suggest that the separated equations provide a starting point for future analytical approximations, investigations of limiting regimes, and the study of horizon-regular solutions in the presence of external tidal fields. The work establishes the separated formulation as a viable alternative to solving the full coupled system, particularly for fundamental modes and near-extremal configurations.

More like this