⚛️ general relativity

Scalaron Quasinormal and Quasi-Bound Modes of Schwarzschild Black Holes in the Bondi-Sachs Formulation of Metric $f(R)$ Gravity

This paper investigates scalaron quasinormal and quasi-bound modes of Schwarzschild black holes within the Bondi-Sachs formulation of metric $f(R)$ gravity by deriving the master equations, computing spectral branches via Leaver's method, and independently validating these frequencies and damping rates through fourth-order characteristic time evolution and matrix-pencil stability scans.

Original authors: Bishop Mongwane, Didam G. A. Duniya, Seturumane Tema, Sanele Ziqubu, David F. Mota

Published 2026-10-06

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Bishop Mongwane, Didam G. A. Duniya, Seturumane Tema, Sanele Ziqubu, David F. Mota

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

Gravity, as we understand it today, is the story of how mass bends the fabric of space and time. For over a century, Albert Einstein's theory of General Relativity has been the master key to this story, passing every test thrown at it, from the motion of planets in our solar system to the collision of black holes detected by modern observatories. Yet, scientists know this story is likely incomplete. Just as Newton's laws were a special case of Einstein's, Einstein's theory might itself be a special case of something deeper. To explore what lies beyond, researchers study modified theories of gravity, looking for subtle deviations that might explain the universe's accelerating expansion or the nature of dark energy. One such family of theories, known as metric f(R) gravity, suggests that the rules governing gravity change slightly depending on the curvature of space itself. In these theories, gravity is not just the geometry of space but also carries a hidden, extra component—a scalar field—that ripples through the cosmos alongside the usual gravitational waves.

When a black hole is disturbed, it does not simply sit still; it rings like a struck bell, emitting gravitational waves that fade away over time. These vibrations, called quasinormal modes, are unique to the black hole's mass and spin, acting as a fingerprint of the spacetime around it. If the laws of gravity are different from Einstein's predictions, the pitch and decay of this ring should change. Furthermore, if the extra scalar field has mass, it can get trapped in a temporary orbit around the black hole, creating a different kind of vibration known as a quasi-bound state. These states are long-lived and decay very slowly, offering a distinct signature that could be spotted in future observations. Understanding exactly how these black holes ring in these modified theories is essential for interpreting the data from gravitational-wave detectors, which are now listening to the universe with unprecedented sensitivity.

In a recent study, a team of researchers set out to calculate these vibrations for a simple, non-rotating black hole within the framework of metric f(R) gravity. They did not rely on the standard methods used for decades but instead employed a specialized mathematical approach known as the Bondi–Sachs formulation. This method views spacetime not as a static grid but as a series of light-like surfaces moving outward from the black hole, a perspective that is particularly well-suited for tracking how radiation escapes to the far reaches of the universe. By applying this framework to the equations of modified gravity, the team derived a master equation that describes how the extra scalar field behaves around the black hole. They found that while the gravitational waves themselves (the tensor modes) behave exactly as Einstein predicted, the extra scalar field behaves like a massive particle, with its vibrations determined by a single effective mass parameter.

The researchers then tackled the problem of finding the specific frequencies of these vibrations using two independent methods. First, they used a frequency-domain technique, which involves solving the equations directly to find the allowed tones, much like finding the resonant notes of a guitar string. Second, they simulated the black hole's response in the time domain, essentially "plucking" the black hole in a computer model and watching how the waves evolve and fade over time. They compared the results from both methods against each other and against existing calculations from other studies. The two approaches agreed remarkably well, confirming that their new mathematical framework was working correctly. They mapped out how the frequencies of these vibrations change as the mass of the scalar field increases, showing that as the field gets heavier, the vibrations become less damped, meaning the black hole rings for a longer time before falling silent.

A significant portion of their work focused on the tricky task of separating these different types of vibrations from a single computer-generated signal. In a real observation, a detector would see a complex mix of signals: the initial burst, the standard ringing, and the lingering, slow decay of the quasi-bound states. The team demonstrated that by carefully analyzing different time windows of the same simulated signal, they could successfully isolate both the standard ringing and the long-lived quasi-bound states without forcing the computer to assume which one was present. They showed that for certain conditions, both types of vibrations could be extracted from the same data stream, proving that the method is robust enough to handle the complexity of real-world signals.

The study also revealed a surprising simplicity in the results. Despite the complexity of the modified gravity theories, the behavior of the black hole's vibrations depended only on a single combination of numbers: the mass of the black hole multiplied by the mass of the scalar field. This means that many different theories of gravity, which might look very different on paper, would produce identical ringdown signals if they share this specific value. This finding creates a degeneracy, where different underlying theories cannot be distinguished just by looking at the frequency of the ring. However, the researchers noted that this simplicity applies only to the linear, or small, disturbances. If the black hole were disturbed more violently, or if the background space had a different curvature, these distinctions might reappear, offering a path to test these theories in more extreme scenarios.

Ultimately, this work provides a crucial set of benchmarks for the future of gravitational-wave astronomy. By establishing a reliable way to calculate and extract these vibrations within a modified gravity framework, the researchers have given observers a clearer map of what to look for. Their results confirm that the characteristic ringing of a black hole is a powerful tool for testing the laws of gravity, capable of distinguishing between standard Einsteinian gravity and more exotic alternatives. The ability to recover both the fast-decaying standard modes and the slow-decaying quasi-bound states from a single signal suggests that future detectors could potentially uncover the presence of these extra scalar fields, opening a new window into the fundamental nature of the universe.

Technical Summary: Scalaron Quasinormal and Quasi-Bound Modes of Schwarzschild Black Holes in the Bondi–Sachs Formulation of Metric f(R) Gravity

Problem Statement
This work investigates the linear perturbations of Schwarzschild black holes within the framework of metric $f(R)$ gravity. Specifically, it focuses on the "scalaron" sector—the additional scalar degree of freedom introduced by replacing the Einstein-Hilbert Lagrangian $R$ with a function $f(R)$ . While the tensor sector of metric $f(R)$ gravity on a Ricci-flat background is known to reduce to General Relativity (GR), the scalar sector obeys a massive Klein-Gordon equation. The authors aim to compute and validate the spectra of quasinormal modes (QNMs) and quasi-bound states (QBSs) for this massive scalar field directly within the Bondi–Sachs (characteristic) formulation of gravity. A key challenge addressed is the rigorous definition and numerical extraction of these modes, particularly the QBSs, which require specific asymptotic decay conditions distinct from the outgoing radiation conditions of QNMs.

Methodology
The study employs a dual approach combining frequency-domain spectral analysis with time-domain characteristic evolution:

Theoretical Formulation:
- Starting from the linearized Bondi–Sachs equations for metric $f(R)$ gravity, the authors derive a closed master equation for the scalar perturbation $R$ .
- They demonstrate that the tensor master equation reduces exactly to the standard GR form on a Schwarzschild background, confirming that the scalaron is the sole carrier of modified gravity signatures in the linear regime.
- The scalar master equation is identified as a massive scalar wave equation with an effective mass $\mu$ , where $\mu^2 = f'(0)/3f''(0)$ . The authors establish that for analytic models on a Ricci-flat background, the linear spectra depend on the specific $f(R)$ theory only through the dimensionless combination $\mu M$ .
Frequency-Domain Analysis:
- The authors utilize Leaver's continued-fraction method to solve the radial master equation.
- They define distinct boundary conditions for QNMs (ingoing at the horizon, outgoing at infinity) and QBSs (ingoing at the horizon, exponentially decaying at infinity).
- These conditions are implemented via branch-dependent Wronskian zeros and corresponding continued-fraction expansions, incorporating Nollert's improvement for high-order convergence.
Time-Domain Evolution:
- A fourth-order convergent characteristic scheme is used to evolve the scalaron perturbation equations on a finite radial domain.
- The inner boundary is treated as a worldtube near the horizon with an ingoing condition, while the outer boundary is a finite radius where no spectral branch (outgoing or decaying) is imposed a priori.
- Waveforms are analyzed using a matrix-pencil stability scan (a Prony-type estimator) to identify poles. This method extracts complex frequencies based on the persistence of modes across variations in fitting windows, matrix ranks, and decimation factors, rather than relying on a single fit.

Key Contributions

Derivation of Master Equations: The paper explicitly derives the scalar master equation from the linearized Bondi–Sachs $f(R)$ system and proves the reduction of the tensor sector to the GR form.
Spectral Degeneracy: It establishes a degeneracy among different analytic $f(R)$ models: at fixed geometric mass $M$ , any two models sharing the same ratio $f'(0)/f''(0)$ (and thus the same $\mu$ ) yield identical linear scalar QNM and QBS spectra.
Numerical Benchmarking: The authors provide a comprehensive set of frequency-domain benchmarks for both QNMs and QBSs across various multipoles ( $\ell$ ) and scalar masses ( $\mu M$ ), validating their continued-fraction implementation against independent literature.
Independent Time-Domain Validation: A significant contribution is the independent recovery of these spectral branches from time-domain evolutions without imposing spectral boundary conditions at the outer edge. This validates the characteristic formulation's ability to capture both outgoing radiation and localized bound states.
Simultaneous Mode Recovery: The paper demonstrates the feasibility of extracting both the fundamental QNM and QBS frequencies from a single finite-radius waveform by utilizing separate time windows (early time for QNMs, late time for QBSs).

Results

Frequency-Domain: The continued-fraction results reproduce literature benchmarks for massless and massive scalar fields with high precision (differences generally $\le 10^{-6}$ ). The study maps the root locus of the fundamental mode and overtones as $\mu M$ increases, showing that increasing mass leads to higher oscillation frequencies and reduced damping (longer lifetimes), approaching a quasi-resonant limit where the imaginary part of the frequency vanishes.
Quasi-Bound States: The computed QBS frequencies agree with published values and asymptotically approach Detweiler's hydrogenic approximation in the small-coupling limit ( $\mu M \ll 1$ ).
Time-Domain: The characteristic evolution exhibits fourth-order convergence. The matrix-pencil analysis successfully recovers the fundamental QNM frequencies consistent with the continued-fraction results within the measured stability scales (MADs).
Robustness: The extracted frequencies remain stable under variations of the inner worldtube location, outer domain size, dissipation strength, and extraction radius.
Simultaneous Extraction: For a selected case ( $\ell=2, \mu M=0.575$ ), the authors successfully recovered both the QNM and QBS complex frequencies from a single waveform, with differences from the frequency-domain benchmarks of $7.35 \times 10^{-4}$ (QNM) and $4.46 \times 10^{-8}$ (QBS).

Significance and Claims
The paper claims to provide a quantitative benchmark set for characteristic evolutions in metric $f(R)$ gravity, validating the Bondi–Sachs framework as a versatile tool for studying linear perturbations. By showing that the scalar spectra depend solely on the effective mass $\mu$ , the work clarifies the specific sector of $f(R)$ theories accessible to linear black hole spectroscopy. The successful independent recovery of both QNMs and QBSs from time-domain data, without imposing spectral conditions at the outer boundary, serves as a robust check on the numerical implementation and the physical interpretation of the modes. The authors modestly note that while the linear spectral problem exhibits degeneracy, higher-order perturbations or non-Ricci-flat backgrounds could break this symmetry and probe higher derivatives of $f(R)$ . The work concludes that the characteristic formulation is well-suited for future extensions to rotating (Kerr) black holes and non-linear regimes.

More like this