Pith. sign in

REVIEW 2 major objections 3 minor 2 cited by

This paper presents the first time-domain implementation of the parametrized modified Teukolsky equation, reproducing frequency-domain quasinormal-mode coefficients to about 2% for low azimuthal numbers while exposing amplitude, phase, and

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-01 02:50 UTC pith:7ILOKJZ4

load-bearing objection First working time-domain solver for the parametrized modified Teukolsky equation, honestly validated against frequency-domain benchmarks; the boundary-reflection worry is overstated because the domain is causally large for the reported times, but the code-release link needs attention. the 2 major comments →

arxiv 2607.25311 v1 pith:7ILOKJZ4 submitted 2026-07-28 gr-qc astro-ph.HEhep-th

Parametrized beyond-Teukolsky framework in the time domain

classification gr-qc astro-ph.HEhep-th
keywords black hole perturbation theoryquasinormal modesmodified Teukolsky equationtime-domain evolutionringdown amplitudes and phaseslate-time tailsbeyond-GR black hole testsKerr spacetime
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper aims to establish that the parametrized beyond-Teukolsky framework—a theory-agnostic way of encoding small deviations from Kerr black-hole dynamics as a power-series deformation of the Teukolsky potential—can be evolved directly in the time domain, not just solved as an eigenvalue problem. It claims that the new time-domain code reproduces the frequency-domain quasinormal-mode coefficients to about 2% for low azimuthal numbers, with higher azimuthal numbers limited by spherical–spheroidal mode mixing, and that it delivers quantities an eigenvalue calculation cannot: the prompt response, ringdown amplitudes and phases, quadratic and mixed coefficient diagnostics, and the onset of the late-time power-law tail. If correct, this gives black-hole spectroscopy a validated, flexible tool for translating parametrized deviations from general relativity into full linear waveforms and for diagnosing when the linear approximation itself breaks down.

Core claim

The central claim is that the parametrized modified Teukolsky equation, with the radial potential augmented by δV(r) = Δ^{-1} Σ_{k=-K}^{4} α^{(k)} (r/r_+)^{k}, can be integrated as a (2+1)-dimensional scattering problem and that the resulting waveforms reproduce the established frequency-domain quasinormal-mode results: for ℓ=2, m∈{-1,0,1} at a=0.25, the extracted linear coefficients d^{(k)}_{ω,nℓm} agree with frequency-domain benchmarks to about 2%, while for m=±2 mode mixing degrades the real part to 5–6%. The same evolutions yield quadratic and mixed coefficients d^{(k1,k2)}_{ω,nℓm} within the linearized model, which the paper proposes as a diagnostic for the validity range of the linear

What carries the argument

The load-bearing object is the modified Teukolsky potential δV(r) = Δ^{-1} Σ_{k=-K}^{4} α^{(k)} (r/r_+)^{k}, a power-series deformation of the standard Kerr potential that is linear in the parameters α^{(k)}; the framework's key identity is the expansion of the quasinormal frequency as ω ≈ ω_0 + Σ_k d^{(k)}_{ω,nℓm} α^{(k)} + ½ Σ_{k1,k2} d^{(k1,k2)}_{ω,nℓm} α^{(k1)}α^{(k2)}, whose coefficients are the objects compared between frequency and time domains. In the time domain, the deformation is inserted into the Teukolsky equation as a radial modification of the angular potential term, and the equation is advanced with a second-order finite-difference scheme on a (r*, θ) grid using ingoing Gauss

Load-bearing premise

The results rest on placing the numerical boundaries at r*=-200M and +400M with ingoing/outgoing advection conditions and trusting that they do not contaminate the computed modes, amplitudes, phases, or tail onsets, even though no domain-size convergence study is reported for the modified potentials and the k=4 deformation is unstable.

What would settle it

Re-run the same Gaussian wave-packet evolutions for representative deformations (e.g., k=-4, α=0.1 and k=1, α=-3) with the inner boundary moved to r*=-400M and the outer to r*=800M, or with hyperboloidal compactified coordinates, and check whether the extracted d^{(k)}_{ω,nℓm}, ringdown amplitudes and phases, and tail-onset times shift by more than the quoted ~2% (or ~10% for tails); if they do, the validation is compromised. A complementary check is to compute k=3 and k=4 coefficients with a frequency-domain code that includes the 1/r and constant large-distance potential terms.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the central claim holds, the time-domain coefficients independently validate the frequency-domain eigenvalue results for low azimuthal numbers and mark the multipole range (roughly |m|≤1 for ℓ=2 at moderate spin) where direct QNM extraction from these simulations is reliable.
  • The quadratic and mixed coefficients give a practical bound on the deformation parameters—the linear approximation remains trustworthy only while |d^{(k)}| exceeds |d^{(k,k)} α^{(k)}|—which can be used to restrict viable beyond-GR parameters before theory-specific modeling.
  • Ringdown amplitudes and phases are predicted to be strongly affected by the deformation, with a sign pattern explainable by the near-horizon plateau: negative α suppresses transmission and boosts the reflected amplitude, positive α does the opposite; this adds observables beyond mode frequencies.
  • For k<3 deformations, the late-time tail keeps the Kerr power-law exponent (≈ t^{-(2ℓ+3)} within ~10% extraction uncertainty), but the onset of the tail can shift significantly because the deformation changes the damping time of the preceding ringdown.
  • The same time-domain tool can be turned to theory-specific modified-gravity models (higher-derivative gravity, dynamical Chern-Simons), non-separable or frequency-dependent deformations, overtones, superradiance, and instabilities, since only the potential terms in the evolution system need changing.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The most consequential unstated implication is that this turns the parametrized framework into a waveform generator: one could in principle fit observed ringdowns directly for amplitudes, phases, and tail onsets, not just frequencies, which may break degeneracies in parameter estimation that frequency-only analyses suffer.
  • The boundary-truncation caveat suggests a specific test that the paper leaves open: implementing hyperboloidal compactified coordinates (which the authors mention as future work) would settle whether the quoted amplitude, phase, and tail-onset shifts—especially the late-time results the paper itself calls 'not very accurate'—survive with boundaries pushed to null infinity.
  • An interesting extension would be to compare the quadratic coefficients computed inside the linear model with the full second-order modified Teukolsky equation of a concrete theory; this would show whether the linear-model validity criterion overestimates or underestimates the true regime of validity, since the paper notes the true second-order terms can be comparable or larger.
  • Because the mode-mixing degradation for m=±2 is attributed to using spin-weighted spherical harmonics rather than spin-weighted spheroidal harmonics as initial data, a testable improvement would be to iterate the spheroidal basis using the deformed QNM frequency; if the 5–6% errors drop toward the ~2% level, the physical-mode-mixing explanation is confirmed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 3 minor

Summary. The paper presents the first time-domain implementation of the parametrized modified Teukolsky equation introduced in Ref. [39]. The authors evolve Gaussian wave packets on a (2+1)-dimensional grid, extract quasinormal-mode frequencies with the Prony method, and compare the linear deformation coefficients d^{(k)}_{\omega,n\ell m} with frequency-domain Leaver benchmarks. They report ≈2% agreement for low azimuthal numbers, with larger (5–6%) errors for m=±2 that they attribute to physical mode mixing. They also compute quadratic and mixed coefficients within the linearized modified-potential model, analyze ringdown amplitudes and phases, and study the late-time tail exponent and its onset. The central claim is that this time-domain framework provides a validated, theory-agnostic tool that goes beyond eigenvalue calculations, giving access to the full linear evolution and to new diagnostics for black-hole spectroscopy.

Significance. If the validation and new results are reliable, this is a useful complement to the eigenvalue-only framework of Ref. [39]. The authors are careful to quantify uncertainties: the m=0 Kerr baseline matches to ≈0.3%, the TD–FD comparison uses 68% HDI error bars, and the degradation for m=±2 is disclosed and supported by a resolution study. The quadratic-coefficient discussion is appropriately hedged: the authors explicitly note that these are coefficients of the linear modified equation, not complete second-order corrections of a generic theory. However, the paper's most novel claims — amplitude/phase trends and tail-onset shifts — depend on the finite-domain numerical treatment, whose convergence is not demonstrated. The validation against frequency-domain QNM frequencies does not cover these new results, because amplitude, phase, and tail onset are sensitive to the entire propagation, including the outer boundary region.

major comments (2)
  1. [Sec. V and Appendix A] The finite-domain claim is load-bearing for the amplitude/phase and tail-onset results. The domain is r*∈[-200M,400M] with advection boundary conditions (Sec. IIB), and Appendix A asserts without a convergence study that the domain is large enough. For k=3 the modified potential decays as 1/r at large distance, so at r*=400M it is still ~1/400 M^{-1}; for k=4 it tends to a constant. An outgoing advection boundary condition is therefore not asymptotically exact for these modified potentials, and partial reflections can contaminate the very quantities Sec. V uses to draw physical conclusions. The QNM-frequency validation in Sec. III does not test this, since ringing is localized near the potential peak. I request either a domain-size convergence study for the k=3 (and, if possible, k=4) modified potentials, or a clear restriction of the claims to the range where boundary effects are demons
  2. [Sec. VB] The tail analysis is explicitly qualified by the statement 'the code is not very accurate at late times'. The claim that the tail exponent is unchanged but its onset is shifted is based on waveforms such as Fig. 9, whose late-time segments are admitted to be noisy. The onset shifts shown for α^(1)=±3 are outside the linear regime and are not accompanied by a quantified uncertainty. Without an estimate of how numerical noise and boundary reflections affect the extracted onset time, the quantitative statement 'can substantially shift the onset' is not yet supported. A convergence study in domain size and resolution, with a defined criterion for the onset time, would make this result robust.
minor comments (3)
  1. [General] There are several typographical and formatting issues, e.g., 'Wignerd-functionsd psq ℓm' in Eq. (19) and inconsistent use of r* vs r˚. Please also define all symbols in Eq. (5) at first use and clarify the units in the discussion of d^{(k)}_{\omega,n\ell m}.
  2. [Sec. III, Figs. 4–5] The labels in the figures and captions could be clearer about the specific values of n, ℓ, and a used. In particular, the hierarchy among m values in Fig. 4 is mentioned in the text but not visible from the figure alone; consider adding markers or a legend that encodes the m values consistently.
  3. [Appendix A] The code availability statement refers to Ref. [77], which appears to be the frequency-domain framework. Please state explicitly whether the time-domain code used for this paper is also publicly available, and if so, where.

Circularity Check

0 steps flagged

No significant circularity: the time-domain solver is cross-checked against an independent frequency-domain computation of the same modified operator.

full rationale

The paper's central derivation is a numerical time-domain integration of the modified Teukolsky equation (11) with the potential deformation δV of Eq. (5). The deformation parameters α^(k) are inputs, not outputs, and the time-domain extraction of QNM frequencies is compared with frequency-domain Leaver-method coefficients from Ref. [39] (Sec. III, Figs. 3–5). This is a cross-check of two independent numerical evaluations of the same linear operator; the reported 2% (m=0,±1) and 5–6% (m=±2) disagreements show the comparison is empirical rather than constructional. The quadratic and mixed coefficients (Sec. IV) are explicitly labeled as diagnostics of the linear model: 'the quadratic coefficients are computed in the linear modified Teukolsky equation. The correct quadratic coefficients will, in general, contain additional contributions associated with the second-order modified Teukolsky equation.' The amplitude/phase and tail-onset results (Sec. V) are direct extractions from the simulated waveforms; no parameter is fitted to those quantities. The only author-overlapping reference is Ref. [39], but it is a published independent computation using Leaver's method with a stated code repository, so it is external evidence rather than a load-bearing self-citation. Appendix A's finite-domain boundary discussion is a numerical accuracy concern, not an identity-by-construction issue. No step in the derivation chain reduces to its own input.

Axiom & Free-Parameter Ledger

4 free parameters · 6 axioms · 0 invented entities

The paper's free inputs are the deformation parameters α^(k) and the hand-chosen numerical/initial-data settings. The domain assumptions are the physical mapping of the δV potential family to beyond-GR physics, the small-coupling structure of the modified equation, the correctness of the shared-author frequency-domain benchmark, and the adequacy of standard boundary/initial-data treatments for the modified operator. No new physical entities are introduced: δV is a phenomenological modeling term from Ref. [39], and the quadratic/mixed coefficients are fit outputs of the linear model, explicitly labeled as such.

free parameters (4)
  • α^(k) deformation parameters = scanned by hand: α ∈ [−0.1, 0.1] for linear fits; α = −3 for the tail illustration
    The independent variables of the model (Eq. 5); all results are response functions of them. Ranges are chosen by hand.
  • Gaussian initial data (width, center, extraction point) = σ = 1M, r*_0 = 12M, extraction at r* = 30M, θ = π/2
    The amplitude and phase results (Sec. VA) depend on this hand-chosen wave packet; the paper acknowledges these quantities are initial-data-dependent.
  • Numerical grid and domain = Nr = 8000, Nθ = 32; r* ∈ [−200M, 400M]
    Chosen by hand; a resolution study supports them, but no domain-size convergence test is reported for the modified potentials, on which the tail-onset claims rest.
  • Prony fit windows and α-grids = hundreds of starting times; 5 α values per fit; quadratic α-grid spacing not stated
    Extraction systematics (Sec. III) depend on these choices; the quadratic-fit grid spacing is not fully specified, limiting exact replication.
axioms (6)
  • domain assumption The modified radial Teukolsky equation (Eq. 1) with δV = (1/Δ) Σ_k α^(k)(r/r+)^k is the correct leading linearized perturbation operator for a broad class of beyond-GR theories.
    Taken from Ref. [39] (overlapping authorship); the interpretation of all results as beyond-GR signatures rests on it. Sec. IIA, Eq. (5).
  • domain assumption Small-coupling approximation: the deformation enters only the potential term, leaving the derivative structure of the equation unchanged.
    Sec. IIB: 'the deformation in Eq. (1) only affects the last term in Eq. (11)'; Sec. IV concedes full second-order theory dynamics would add terms of comparable size.
  • standard math Leaver's continued-fraction method gives correct frequency-domain eigenvalues and d^(k) coefficients for the modified radial equation.
    Used as the validation benchmark (Sec. III, Refs. [39,76]); independent of the time-domain code but computed by overlapping authors.
  • standard math The modified Lax-Wendroff scheme remains second-order convergent and stable when δV is added to the potential.
    Appendix A: hyperbolicity is preserved because matrix M is unchanged, but no stability proof or convergence study for the modified equation is given beyond one resolution check.
  • domain assumption The Price-law tail ~ t^−(2ℓ+3) governs late-time decay, and k<3 deformations do not change the exponent.
    Sec. VB; compared with ~10% uncertainty and consistent with Ref. [58]. The tail-onset shift claim depends on identifying the power-law segment in noise-limited data.
  • domain assumption The near-horizon modified field behaves as Δ^−s(r−r+)^{−iσ} with σ in Eq. (A21), so standard QNM-type advection boundary conditions at finite r* remain adequate.
    Appendix A; no domain-size convergence demonstration is provided for modified potentials.

pith-pipeline@v1.3.0-alltime-deepseek · 19718 in / 30339 out tokens · 276468 ms · 2026-08-01T02:50:02.296046+00:00 · methodology

0 comments
read the original abstract

Modifications to General Relativity can significantly alter the perturbative response of black holes, leaving imprints on quasinormal-mode spectra, waveform amplitudes and phases, and late-time tails. The parametrized beyond-Teukolsky framework was introduced to capture possible deviations from Kerr dynamics, but the ringdown has so far only been explored as an eigenvalue problem. We present the first time-domain implementation of this framework and perform (2+1)-dimensional scattering experiments with Gaussian wave packets. This approach provides the full linear evolution of the perturbation, from the initial prompt response through the ringdown and into the late-time regime. Using frequency-domain eigenvalue results as benchmarks, we find excellent agreement for low azimuthal numbers, while higher azimuthal numbers are affected by mode mixing, which limits the precision of the extracted modes. We further use the time-domain waveforms to estimate quadratic and mixed coefficients within the linearized modified-potential model, providing a diagnostic for the regime of validity of the linear approximation. Beyond mode frequencies, we show that the deformation parameters can strongly affect the ringdown amplitude and phase, with trends that can be understood from the near-horizon structure of the modified potential. We also analyze the late- time behavior, finding that near-horizon deformations leave the tail exponent unchanged but can substantially shift the onset of the power-law decay. These results demonstrate that time-domain evolutions provide a complementary and flexible framework for testing parametrized deviations from General Relativity in black-hole perturbation theory.

Figures

Figures reproduced from arXiv: 2607.25311 by Ciro De Simone, Kostas D. Kokkotas, Salvatore Capozziello, Sebastian H. V\"olkel.

Figure 1
Figure 1. Figure 1: FIG. 1. Real and imaginary parts of the modified Teukolsky potential for angular momentum [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Waveforms obtained for [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Relative error [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of the coefficients [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison of the linear coefficients [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Real and imaginary part of the quadratic coefficients [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Real and imaginary part of the mixed coefficients [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Fractional deviation of the amplitude [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Waveforms for [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Comparison of the linear coefficients from frequency and time domain at [PITH_FULL_IMAGE:figures/full_fig_p016_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Color map of the relative error ( [PITH_FULL_IMAGE:figures/full_fig_p016_11.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Discrete symmetries of modified Teukolsky equations

    gr-qc 2026-07 conditional novelty 5.0

    Complex modifications of the Teukolsky potential break m=0 QNM degeneracy and can inject non-physical mode branches when frequency-domain multipole-dependent potentials are evolved in time.

  2. Constraining deviations from the Teukolsky equation with GW250114

    gr-qc 2026-07 conditional novelty 5.0

    GW250114's fundamental ringdown mode bounds theory-agnostic deviations from the Teukolsky equation to be consistent with zero at characteristic scales of 60-100 km.

Reference graph

Works this paper leans on

97 extracted references · 69 linked inside Pith · cited by 2 Pith papers

  1. [1]

    A. G. Abacet al.(LIGO Scientific, Virgo, KAGRA), GW250114: Testing Hawking’s Area Law and the Kerr Nature of Black Holes, Phys. Rev. Lett.135, 111403 (2025), arXiv:2509.08054 [gr-qc]

  2. [2]

    EinsteinWaves

    K.D.K. and S.H.V. acknowledge funding from the European Union’s Horizon MSCA-2022 research and in- novation programme “EinsteinWaves” under grant agree- ment No. 101131233. The authors acknowledge support by the High Performance and Cloud Computing Group at the Zentrum für Datenverarbeitung of the University of Tübingen, the state of Baden-Württemberg thr...

  3. [3]

    A. G. Abacet al.(LIGO Scientific, Virgo, KAGRA), Black Hole Spectroscopy and Tests of General Relativity with GW250114, Phys. Rev. Lett.136, 041403 (2026), arXiv:2509.08099 [gr-qc]

  4. [4]

    Bertiet al., Black hole spectroscopy: from theory to experiment, Class

    E. Bertiet al., Black hole spectroscopy: from theory to experiment, Class. Quant. Grav.43, 123001 (2026), arXiv:2505.23895 [gr-qc]

  5. [5]

    Dreyer, B

    O. Dreyer, B. J. Kelly, B. Krishnan, L. S. Finn, D. Gar- rison, and R. Lopez-Aleman, Black hole spectroscopy: Testing general relativity through gravitational wave ob- servations, Class. Quant. Grav.21, 787 (2004), arXiv:gr- qc/0309007

  6. [6]

    R. P. Kerr, Gravitational field of a spinning mass as an example of algebraically special metrics, Phys. Rev. Lett. 11, 237 (1963)

  7. [7]

    C. V. Vishveshwara, Scattering of Gravitational Radi- ation by a Schwarzschild Black-hole, Nature227, 936 (1970)

  8. [8]

    E.W.Leaver,Spectraldecompositionoftheperturbation responseoftheSchwarzschildgeometry,Phys.Rev.D34, 384 (1986)

  9. [9]

    Andersson, Evolving test fields in a black hole geom- etry, Phys

    N. Andersson, Evolving test fields in a black hole geom- etry, Phys. Rev. D55, 468 (1997), arXiv:gr-qc/9607064

  10. [10]

    K. D. Kokkotas and B. G. Schmidt, Quasinormal modes of stars and black holes, Living Rev. Rel.2, 2 (1999), arXiv:gr-qc/9909058

  11. [11]

    Berti, V

    E. Berti, V. Cardoso, and A. O. Starinets, Quasinormal modes of black holes and black branes, Class. Quant. Grav.26, 163001 (2009), arXiv:0905.2975 [gr-qc]

  12. [12]

    R. H. Price, Nonspherical perturbations of relativistic gravitational collapse. 1. Scalar and gravitational pertur- bations, Phys. Rev. D5, 2419 (1972)

  13. [13]

    Gundlach, R

    C. Gundlach, R. H. Price, and J. Pullin, Late time behav- ior of stellar collapse and explosions: 2. Nonlinear evolu- tion, Phys. Rev. D49, 890 (1994), arXiv:gr-qc/9307010

  14. [14]

    S. A. Teukolsky, Perturbations of a rotating black hole. 1. Fundamental equations for gravitational electromagnetic and neutrino field perturbations, Astrophys. J.185, 635 (1973)

  15. [15]

    J. W. Moffat, Scalar-tensor-vector gravity theory, JCAP 03, 004, arXiv:gr-qc/0506021

  16. [16]

    Capozziello and M

    S. Capozziello and M. De Laurentis, Extended Theories of Gravity, Phys. Rept.509, 167 (2011), arXiv:1108.6266 [gr-qc]

  17. [17]

    Bertiet al., Testing General Relativity with Present and Future Astrophysical Observations, Class

    E. Bertiet al., Testing General Relativity with Present and Future Astrophysical Observations, Class. Quant. Grav.32, 243001 (2015), arXiv:1501.07274 [gr-qc]

  18. [18]

    Baracket al., Black holes, gravitational waves and fundamental physics: a roadmap, Class

    L. Baracket al., Black holes, gravitational waves and fundamental physics: a roadmap, Class. Quant. Grav. 36, 143001 (2019), arXiv:1806.05195 [gr-qc]

  19. [19]

    P. G. S. Fernandes, P. Carrilho, T. Clifton, and D. J. Mulryne, The 4D Einstein–Gauss–Bonnet theory of grav- ity: a review, Class. Quant. Grav.39, 063001 (2022), arXiv:2202.13908 [gr-qc]

  20. [20]

    S. R. Coleman, J. Preskill, and F. Wilczek, Quantum hair on black holes, Nucl. Phys. B378, 175 (1992), arXiv:hep- th/9201059

  21. [21]

    T. P. Sotiriou, Black Holes and Scalar Fields, Class. Quant. Grav.32, 214002 (2015), arXiv:1505.00248 [gr- qc]

  22. [22]

    C. A. R. Herdeiro and E. Radu, Asymptotically flat black holes with scalar hair: a review, Int. J. Mod. Phys. D24, 1542014 (2015), arXiv:1504.08209 [gr-qc]. 12

  23. [23]

    Israel, Event horizons in static vacuum space-times, Phys

    W. Israel, Event horizons in static vacuum space-times, Phys. Rev.164, 1776 (1967)

  24. [24]

    Carter, Axisymmetric Black Hole Has Only Two De- grees of Freedom, Phys

    B. Carter, Axisymmetric Black Hole Has Only Two De- grees of Freedom, Phys. Rev. Lett.26, 331 (1971)

  25. [25]

    S. W. Hawking, Black holes in general relativity, Com- mun. Math. Phys.25, 152 (1972)

  26. [26]

    D. C. Robinson, Uniqueness of the Kerr black hole, Phys. Rev. Lett.34, 905 (1975)

  27. [27]

    Cardoso and L

    V. Cardoso and L. Gualtieri, Perturbations of Schwarzschild black holes in Dynamical Chern-Simons modified gravity, Phys. Rev. D80, 064008 (2009), [Er- ratum: Phys.Rev.D 81, 089903 (2010)], arXiv:0907.5008 [gr-qc]

  28. [28]

    Franchini and S

    N. Franchini and S. H. Völkel, Testing General Rel- ativity with Black Hole Quasi-normal Modes (2024) arXiv:2305.01696 [gr-qc]

  29. [29]

    S. H. Völkel and K. D. Kokkotas, Scalar Fields and Parametrized Spherically Symmetric Black Holes: Can one hear the shape of space-time?, Phys. Rev. D100, 044026 (2019), arXiv:1908.00252 [gr-qc]

  30. [30]

    Pani and A

    P. Pani and A. P. Sanna, Scalar shortcut to beyond-Kerr ringdown tests and their complementarity with black- hole shadow observations, (2026), arXiv:2603.08782 [gr- qc]

  31. [31]

    S. R. Dolan, The Quasinormal Mode Spectrum of a Kerr Black Hole in the Eikonal Limit, Phys. Rev. D82, 104003 (2010), arXiv:1007.5097 [gr-qc]

  32. [32]

    R. A. Konoplya and Z. Stuchlík, Are eikonal quasinor- mal modes linked to the unstable circular null geodesics?, Phys. Lett. B771, 597 (2017), arXiv:1705.05928 [gr-qc]

  33. [33]

    D. Li, P. Wagle, Y. Chen, and N. Yunes, Perturbations of Spinning Black Holes beyond General Relativity: Modi- fiedTeukolskyEquation,Phys.Rev.X13,021029(2023), arXiv:2206.10652 [gr-qc]

  34. [34]

    P. A. Cano, K. Fransen, T. Hertog, and S. Maenaut, Uni- versal Teukolsky equations and black hole perturbations in higher-derivative gravity, Phys. Rev. D108, 024040 (2023), arXiv:2304.02663 [gr-qc]

  35. [35]

    Regge and J

    T. Regge and J. A. Wheeler, Stability of a Schwarzschild singularity, Phys. Rev.108, 1063 (1957)

  36. [36]

    F. J. Zerilli, Gravitational field of a particle falling in a schwarzschild geometry analyzed in tensor harmonics, Phys. Rev. D2, 2141 (1970)

  37. [37]

    McManus, E

    R. McManus, E. Berti, C. F. B. Macedo, M. Kimura, A. Maselli, and V. Cardoso, Parametrized black hole quasinormal ringdown. II. Coupled equations and quadratic corrections for nonrotating black holes, Phys. Rev. D100, 044061 (2019), arXiv:1906.05155 [gr-qc]

  38. [38]

    Cardoso, M

    V. Cardoso, M. Kimura, A. Maselli, E. Berti, C. F. B. Macedo, and R. McManus, Parametrized black hole quasinormal ringdown: Decoupled equations for non- rotating black holes, Phys. Rev. D99, 104077 (2019), arXiv:1901.01265 [gr-qc]

  39. [39]

    Franchini and S

    N. Franchini and S. H. Völkel, Parametrized quasinormal mode framework for non-Schwarzschild metrics, Phys. Rev. D107, 124063 (2023), arXiv:2210.14020 [gr-qc]

  40. [40]

    P. A. Cano, L. Capuano, N. Franchini, S. Maenaut, and S. H. Völkel, Parametrized quasinormal mode frame- work for modified Teukolsky equations, Phys. Rev. D 110, 104007 (2024), [Erratum: Phys.Rev.D 113, 069902 (2026)], arXiv:2407.15947 [gr-qc]

  41. [41]

    Pezzella, K

    L. Pezzella, K. Destounis, A. Maselli, and V. Car- doso, Quasinormal modes of black holes embedded in halos of matter, Phys. Rev. D111, 064026 (2025), arXiv:2412.18651 [gr-qc]

  42. [42]

    T. F. M. Spieksma, V. Cardoso, G. Carullo, M. Della Rocca, and F. Duque, Black Hole Spectroscopy in Environments: Detectability Prospects, Phys. Rev. Lett.134, 081402 (2025), arXiv:2409.05950 [gr-qc]

  43. [43]

    Thomopoulos, S

    S. Thomopoulos, S. H. Völkel, and H. P. Pfeiffer, Ringdown spectroscopy of phenomenologically modi- fied black holes, Phys. Rev. D112, 064054 (2025), arXiv:2504.17848 [gr-qc]

  44. [44]

    Hussain and A

    A. Hussain and A. Zimmerman, Approach to computing spectral shifts for black holes beyond Kerr, Phys. Rev. D 106, 104018 (2022), arXiv:2206.10653 [gr-qc]

  45. [45]

    D. Li, A. Hussain, P. Wagle, Y. Chen, N. Yunes, and A. Zimmerman, Isospectrality breaking in the Teukol- sky formalism, Phys. Rev. D109, 104026 (2024), arXiv:2310.06033 [gr-qc]

  46. [46]

    Wagle, D

    P. Wagle, D. Li, Y. Chen, and N. Yunes, Perturbations of spinning black holes in dynamical Chern-Simons grav- ity: Slow rotation equations, Phys. Rev. D109, 104029 (2024), arXiv:2311.07706 [gr-qc]

  47. [47]

    D. Li, P. Wagle, Y. Chen, and N. Yunes, Perturbations of spinning black holes in dynamical Chern-Simons grav- ity: Slow rotation quasinormal modes, Phys. Rev. D112, 044005 (2025), arXiv:2503.15606 [gr-qc]

  48. [48]

    S. H. Völkel, N. Franchini, and E. Barausse, Theory- agnostic reconstruction of potential and couplings from quasinormal modes, Phys. Rev. D105, 084046 (2022), arXiv:2202.08655 [gr-qc]

  49. [49]

    S. H. Völkel, N. Franchini, E. Barausse, and E. Berti, Constraining modifications of black hole perturbation potentials near the light ring with quasinormal modes, Phys. Rev. D106, 124036 (2022), arXiv:2209.10564 [gr- qc]

  50. [50]

    Hirano, M

    S. Hirano, M. Kimura, M. Yamaguchi, and J. Zhang, Parametrized black hole quasinormal ringdown formal- ism for higher overtones, Phys. Rev. D110, 024015 (2024), arXiv:2404.09672 [gr-qc]

  51. [51]

    Kimura, Note on the parametrized black hole quasi- normal ringdown formalism, Phys

    M. Kimura, Note on the parametrized black hole quasi- normal ringdown formalism, Phys. Rev. D101, 064031 (2020), arXiv:2001.09613 [gr-qc]

  52. [52]

    Krivan, P

    W. Krivan, P. Laguna, P. Papadopoulos, and N. Ander- sson, Dynamics of perturbations of rotating black holes, Phys. Rev. D56, 3395 (1997), arXiv:gr-qc/9702048

  53. [53]

    Pazos-Avalos and C

    E. Pazos-Avalos and C. O. Lousto, Numerical integration of the Teukolsky equation in the time domain, Phys. Rev. D72, 084022 (2005), arXiv:gr-qc/0409065

  54. [54]

    D. D. Doneva, L. G. Collodel, C. J. Krüger, and S. S. Yazadjiev, Black hole scalarization induced by the spin: 2+1 time evolution, Phys. Rev. D102, 104027 (2020), arXiv:2008.07391 [gr-qc]

  55. [55]

    Pedrotti and S

    D. Pedrotti and S. Vagnozzi, Quasinormal modes-shadow correspondence for rotating regular black holes, Phys. Rev. D110, 084075 (2024), arXiv:2404.07589 [gr-qc]

  56. [56]

    De Simone, S

    C. De Simone, S. H. Völkel, K. D. Kokkotas, V. De Falco, and S. Capozziello, Confronting eikonal and post-Kerr methods with numerical evolution of scalar field pertur- bations in spacetimes beyond Kerr, Phys. Rev. D113, 104004 (2026), arXiv:2601.09607 [gr-qc]

  57. [57]

    H. O. Silva, G. Tambalo, K. Glampedakis, K. Yagi, and J. Steinhoff, Quasinormal modes and their excitation be- yond general relativity, Phys. Rev. D110, 024042 (2024), arXiv:2404.11110 [gr-qc]

  58. [58]

    H. O. Silva, G. Tambalo, K. Glampedakis, and K. Yagi, Quasinormal modes and their excitation be- 13 yond general relativity. II. Isospectrality loss in gravi- tational waveforms, Phys. Rev. D113, 084012 (2026), arXiv:2601.13411 [gr-qc]

  59. [59]

    R. F. Rosato and P. Pani, Universality of late-time ringdown tails, Phys. Rev. D112, 024080 (2025), arXiv:2505.08877 [gr-qc]

  60. [60]

    S. A. Teukolsky, Rotating black holes - separable wave equations for gravitational and electromagnetic pertur- bations, Phys. Rev. Lett.29, 1114 (1972)

  61. [61]

    P. A. Cano, L. Capuano, N. Franchini, S. Maenaut, and S. H. Völkel, Higher-derivative corrections to the Kerr quasinormal mode spectrum, Phys. Rev. D110, 124057 (2024), arXiv:2409.04517 [gr-qc]

  62. [62]

    R. Tang, N. Franchini, S. H. Völkel, and E. Berti, Quasi- normal modes of rotating black holes beyond general rel- ativity in the WKB approximation, Phys. Rev. D113, 104052 (2026), arXiv:2512.17786 [gr-qc]

  63. [63]

    P. A. Cano and M. David, Teukolsky equation for near- extremal black holes beyond general relativity: Near- horizon analysis, Phys. Rev. D110, 064067 (2024), arXiv:2407.02017 [gr-qc]

  64. [64]

    C.DeSimone,DiscretesymmetriesofmodifiedTeukolsky equations, in preparation

  65. [65]

    Krivan, P

    W. Krivan, P. Laguna, and P. Papadopoulos, Dynamics of scalar fields in the background of rotating black holes, Phys. Rev. D54, 4728 (1996), arXiv:gr-qc/9606003

  66. [66]

    W. B. Campbell, Tensor and spinor spherical harmonics and the spin-s harmonics y-s-lm(theta, phi), J. Math. Phys.12, 1763 (1971)

  67. [67]

    A. S. Kholyavka, J. A. L. Vega, S. G. Gómez, X. J. Forteza, and S. Datta, Shaping black hole resonances I. Black hole ringdown as a spectral filtering process, (2026), arXiv:2605.24704 [gr-qc]

  68. [68]

    Berti, V

    E. Berti, V. Cardoso, and M. Casals, Eigenvalues and eigenfunctions of spin-weighted spheroidal harmonics in four and higher dimensions, Phys. Rev. D73, 024013 (2006), [Erratum: Phys.Rev.D 73, 109902 (2006)], arXiv:gr-qc/0511111

  69. [69]

    Hod, Mode coupling in rotating gravitational col- lapse of a scalar field, Phys

    S. Hod, Mode coupling in rotating gravitational col- lapse of a scalar field, Phys. Rev. D61, 024033 (2000), arXiv:gr-qc/9902072

  70. [70]

    Baibhav, E

    V. Baibhav, E. Berti, V. Cardoso, and G. Khanna, Black Hole Spectroscopy: Systematic Errors and Ring- down Energy Estimates, Phys. Rev. D97, 044048 (2018), arXiv:1710.02156 [gr-qc]

  71. [71]

    Bhagwat, X

    S. Bhagwat, X. J. Forteza, P. Pani, and V. Ferrari, Ringdown overtones, black hole spectroscopy, and no- hair theorem tests, Phys. Rev. D101, 044033 (2020), arXiv:1910.08708 [gr-qc]

  72. [72]

    Baibhav, M

    V. Baibhav, M. H.-Y. Cheung, E. Berti, V. Cardoso, G. Carullo, R. Cotesta, W. Del Pozzo, and F. Duque, Agnostic black hole spectroscopy: Quasinormal mode content of numerical relativity waveforms and limits of validity of linear perturbation theory, Phys. Rev. D108, 104020 (2023), arXiv:2302.03050 [gr-qc]

  73. [73]

    P. J. Nee, S. H. Völkel, and H. P. Pfeiffer, Role of black hole quasinormal mode overtones for ringdown analysis, Phys. Rev. D108, 044032 (2023), arXiv:2302.06634 [gr- qc]

  74. [74]

    H. Zhu, J. L. Ripley, A. Cárdenas-Avendaño, and F. Pre- torius, Challenges in quasinormal mode extraction: Per- spectivesfromnumericalsolutionstotheTeukolskyequa- tion, Phys. Rev. D109, 044010 (2024), arXiv:2309.13204 [gr-qc]

  75. [75]

    Giesleret al., Overtones and nonlinearities in binary black hole ringdowns, Phys

    M. Giesleret al., Overtones and nonlinearities in binary black hole ringdowns, Phys. Rev. D111, 084041 (2025), arXiv:2411.11269 [gr-qc]

  76. [76]

    S. H. Völkel and A. Dhani, Quantifying systematic biases in black hole spectroscopy, Phys. Rev. D112, 084076 (2025), arXiv:2507.22122 [gr-qc]

  77. [77]

    E. W. Leaver, An Analytic representation for the quasi normal modes of Kerr black holes, Proc. Roy. Soc. Lond. A402, 285 (1985)

  78. [78]

    github, parametrized_qnm_framework (2024)

  79. [79]

    Berti, V

    E. Berti, V. Cardoso, J. A. Gonzalez, and U. Sperhake, Mining information from binary black hole mergers: A Comparison of estimation methods for complex exponen- tials in noise, Phys. Rev. D75, 124017 (2007), arXiv:gr- qc/0701086

  80. [80]

    Berti, V

    E. Berti, V. Cardoso, M. H.-Y. Cheung, F. Di Filippo, F. Duque, P. Martens, and S. Mukohyama, Stability of the fundamental quasinormal mode in time-domain observations against small perturbations, Phys. Rev. D 106, 084011 (2022), arXiv:2205.08547 [gr-qc]

Showing first 80 references.