REVIEW
Scalaron Quasinormal and Quasi-Bound Modes of Schwarzschild Black Holes in the Bondi-Sachs Formulation of Metric $f(R)$ Gravity
General Relativity and Quantum Cosmology · Submitted 5 Oct 2026 · 87 references
Pith's verdict: A careful, honest paper whose Bondi-Sachs scalaron derivation and QNM benchmarks are solid; the real soft spot is the quasi-bound time-domain recovery, where ultra-tight clusters sit six orders of magnitude below a systematic TD–CF offset that no scan explains.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
In metric f(R) gravity, the scalaron around a Schwarzschild black hole is a massive scalar with effective mass $\mu^2 = f'(0)/(3f''(0))$, so the entire linear quasinormal and quasi-bound spectrum is fixed by the dimensionless product…
The reading
What carries the argument
The load-bearing identity is the linearised trace equation $3 f''(0)\,\Box R - f'(0)\,R = 0$, which rearranges into the Klein–Gordon equation $(\Box - \mu^2)R = 0$ with $\mu^2 = f'(0)/(3f''(0))$; this single relation converts the scalaron into a massive scalar on the fixed Schwarzschild background, supplies the closed scalar master equation (25), and makes every computed linear frequency depend on the theory only through $\mu M$. The companion elimination treats the metric-variable subsystem (20b)–(20d) as an algebraic system for $U$ and its derivatives; imposing the compatibility condition $\partial_r U - U_{,r} = 0$ collapses, via the hypersurface relation (20a), into the closed tensor master equation (24) in its GR form, so the spin-2 spectrum is untouched by f(R) at this order. Numerically, the scalar QNM and quasi-bound conditions are the vanishing of the two Wronskians $W_+(\rho)$ and $W_-(\rho)$ on their respective $\kappa$-sheets, enforced through three-term continued-fraction recurrences (the standard series-expansion eigenvalue technique, with a large-overtone termination); the independent time-domain check is a fourth-order Runge–Kutta evolution of the reformulated first-order system on a logarithmic radial grid, with Euler–Maclaurin-corrected trapezoidal integration outward from an inner worldtube placed exponentially close to the horizon and with no spectral branch imposed at the outer boundary.
Extended reading notes
Core claim
The core discovery, on the paper's own terms, is a reduction: linearising the Bondi–Sachs characteristic equations of metric f(R) gravity about a Ricci-flat Schwarzschild background splits the system into a closed massive-scalar problem for the Ricci-scalar perturbation $R$ and a tensor problem that reduces to general relativity's. The trace equation gives $(\Box - \mu^2)R = 0$ with $\mu^2 = f'(0)/(3f''(0))$, so the scalar master equation (25) is exactly the massive scalar equation on Schwarzschild; eliminating the metric variables and using the hypersurface relation yields the tensor master equation (24) in its standard GR form. Because coefficients of $R^3$ and higher never enter at first order, the scalar quasinormal-mode and quasi-bound frequencies depend on the theory only through $\mu M$: at fixed black-hole mass $M$, models sharing $f'(0)/f''(0)$ have identical scalar spectra. The numerically realised spectra are therefore presented as universal massive-scalar Schwarzschild spectra parameterised by $\mu M$; both the outgoing QNM branch (damped oscillations radiating to infinity) and the decaying quasi-bound branch (long-lived, spatially localised states) are computed as Wronskian zero conditions $W_\pm(\rho)=0$ enforced through the continued-fraction method, and a characteristic time evolution with matrix-pencil extraction independently recovers both branches from a single waveform, without imposing either spectral branch at the outer boundary.
Load-bearing premise
The whole spectral equivalence rests on the linearised trace equation $3f''(0)\,\Box R - f'(0)\,R = 0$, which presumes $f$ analytic near $R = 0$ with $f(0) = 0$, $f'(0) \ne 0$ and $f''(0) > 0$, and which discards all $O(\epsilon^2)$ terms — relax any of these conditions and the massive-scalar identification, the $\mu M$ degeneracy, and the GR tensor reduction no longer follow.
Editorial extensions
If this is right
- Any two analytic f(R) models with the same ratio $f'(0)/f''(0)$ predict identical linear scalar quasinormal and quasi-bound spectra on a Schwarzschild background, so linear-order scalaron ringdown data alone cannot distinguish them.
- The tensor master equation reduces exactly to its general-relativistic form, so the spin-2 QNM spectrum of Schwarzschild in metric f(R) gravity is the standard GR spectrum; modified-gravity signatures must show up first in the scalar polarisation.
- Along each overtone branch, increasing $\mu M$ drives the damping toward the quasi-resonant limit $-\operatorname{Im}(2M\omega)\to 0^-$ while leaving the oscillation frequency non-zero, so scalaron modes become arbitrarily long-lived at large effective mass.
- In the small-coupling regime the quasi-bound damping rates tend to the matched-asymptotic hydrogenic formula, giving an analytic anchor for the decaying-branch roots as $\mu M\to 0$.
- A single characteristic waveform at one extraction radius can deliver both the fundamental QNM and quasi-bound complex frequencies from separate fitting windows; the demonstrated $\ell=2$, $\mu M=0.575$ case matches the continued-fraction values to $7.35\times10^{-4}$ (QNM) and $4.46\times10^{-8}$ (QBS) in $2M\omega$.
Reading between the lines
- If the degeneracy holds, scalaron ringdown templates can be shared across entire f(R) model families at linear order; telling individual theories apart would then require sourced perturbations, nonlinear scalaron–tensor coupling, or a non-Ricci-flat background — the directions the authors flag in their concluding remarks.
- The demonstrated single-waveform, single-radius recovery of both mode families implies a reusable extraction protocol for massive fields in null-characteristic codes: tune early windows to the ringdown, late windows to the long-lived mode, and let persistence across window, rank, and decimation variations certify the poles — a procedure that could extend directly to massive vector or Dirac perturb
- The unresolved $\ell=2$, $\mu M=0.25$ quasi-bound decay rate (frequency recovered, damping not) points to a general measurement limit: for extremely long-lived modes a finite waveform constrains the oscillation accurately but the lifetime only weakly, which bounds how finely quasi-bound lifetimes can be probed.
Original abstract
We study scalaron quasinormal modes and quasi-bound states of Schwarzschild black holes directly within the Bondi-Sachs formulation of metric $f(R)$ gravity, combining frequency-domain analyses with characteristic time evolution. Starting from the linearised characteristic equations, we obtain the massive scalar master equation and show that the tensor master equation takes the same form as in general relativity. For the analytic models considered here, the linear scalar spectra depend on the $f(R)$ theory only through the dimensionless combination $\mu M$. We compute both quasinormal and quasi-bound spectral branches using Leaver's continued-fraction method and compare them with independent literature benchmarks. To complement the frequency domain analysis, we evolve the scalaron perturbation equations using a fourth order convergent characteristic scheme. Applying matrix-pencil stability scans to the resulting waveforms, we independently recover fundamental quasinormal frequencies consistent with the continued-fraction results within the measured extraction-stability scales. Finally, for a selected finite-mass case, we demonstrate the extraction of both oscillation frequencies and damping rates of the fundamental quasinormal mode and quasi-bound state from separate time windows of the same waveform at one extraction radius, without imposing either spectral branch at the finite outer boundary.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We study scalaron quasinormal modes and quasi-bound states of Schwarzschild black holes directly within the Bondi-Sachs formulation of metric $f(R)$ gravity, combining frequency-domain analyses with characteristic time evolution. Starting from the linearised characteristic equations, we obtain the massive scalar master equation and show that the tensor master equation takes the same form as in general relativity. For the analytic models considered here, the linear scalar spectra depend on the $f(R)$ theory only through the dimensionless combination $\mu M$. We compute both quasinormal and quasi-bound spectral branches using Leaver's continued-fraction method and compare them with independent literature benchmarks. To complement the frequency domain analysis, we evolve the scalaron perturbation equations using a fourth order convergent characteristic scheme. Applying matrix-pencil stability scans to the resulting waveforms, we independently recover fundamental quasinormal frequencies consistent with the continued-fraction results within the measured extraction-stability scales. Finally, for a selected finite-mass case, we demonstrate the extraction of both oscillation frequencies and damping rates of the fundamental quasinormal mode and quasi-bound state from separate time windows of the same waveform at one extraction radius, without imposing either spectral branch at the finite outer boundary.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We study scalaron quasinormal modes and quasi-bound states of Schwarzschild black holes directly within the Bondi-Sachs formulation of metric $f(R)$ gravity, combining frequency-domain analyses with characteristic time evolution. | contribution 1.00 | 0.95 | 0.12 | asserted 0.98 |
| s1 | Starting from the linearised characteristic equations, we obtain the massive scalar master equation and show that the tensor master equation takes the same form as in general relativity. | method 0.58 | 0.94 | 0.08 | asserted 1.00 |
| s2 | For the analytic models considered here, the linear scalar spectra depend on the $f(R)$ theory only through the dimensionless combination $\mu M$. | result 0.95 | 0.79 | 0.18 | asserted 1.00 |
| s3 | We compute both quasinormal and quasi-bound spectral branches using Leaver's continued-fraction method and compare them with independent literature benchmarks. | method 0.89 | 0.95 | 0.09 | asserted 1.00 |
| s4 | To complement the frequency domain analysis, we evolve the scalaron perturbation equations using a fourth order convergent characteristic scheme. | method 0.99 | 0.96 | 0.03 | asserted 0.98 |
| s5 | Applying matrix-pencil stability scans to the resulting waveforms, we independently recover fundamental quasinormal frequencies consistent with the continued-fraction results within the measured extraction-stability scales. | result 0.88 | 0.96 | 0.15 | asserted 0.95 |
| s6 | Finally, for a selected finite-mass case, we demonstrate the extraction of both oscillation frequencies and damping rates of the fundamental quasinormal mode and quasi-bound state from separate time windows of the same waveform at one extraction radius, without imposing either spectral branch at the finite outer boundary. | result 0.61 | 0.97 | 0.29 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:2610.06370 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.