Recognition: unknown
Beyond Three Terms: Continued Fractions for Rotating Black Holes in Modified Gravity
Pith reviewed 2026-05-10 04:04 UTC · model grok-4.3
The pith
A reduction scheme converts arbitrary N-term recurrence relations into three-term form, extending continued fractions to black-hole perturbations in modified gravity.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors present a general reduction scheme that maps arbitrary scalar and matrix N-term recurrence relations to equivalent three-term relations. In the concrete application to slowly rotating black holes in dynamical Chern-Simons gravity, the polar sector yields a 16-term decoupled scalar recurrence and the axial sector a 12-term coupled matrix recurrence. After the reduction is applied, both systems are solved with continued fractions. The resulting frequencies for the fundamental (ℓ,m)=(2,2) mode agree with independent results obtained from eigenvalue-perturbation and metric/spectral methods across the parameter range examined.
What carries the argument
The general reduction scheme that transforms N-term scalar and matrix recurrence relations into three-term form while retaining the information required for continued-fraction convergence.
If this is right
- Quasinormal modes of slowly rotating black holes in dynamical Chern-Simons gravity become computable with continued fractions after the reduction is performed.
- Frequencies for the fundamental (2,2) mode match independent calculations from eigenvalue-perturbation and metric/spectral methods over the studied parameter range.
- The same reduction procedure applies to any perturbation problem in modified gravity that generates higher-order or coupled recurrences.
- Precision ringdown calculations become feasible for testing gravity with current and future gravitational-wave observations.
Where Pith is reading between the lines
- The reduction may extend to perturbation equations in other modified-gravity theories whose recurrences exceed three terms.
- Numerical stability of the reduced relations should be checked in regimes of rapid rotation or strong coupling not covered by the slow-rotation examples.
- Existing Leaver implementations could be reused after the reduction step, lowering the barrier to new calculations.
- The approach may simplify mode calculations for non-Kerr backgrounds once the original recurrences are obtained.
Load-bearing premise
The reduction from N-term to three-term relations preserves the accuracy, convergence, and numerical stability of the continued-fraction method.
What would settle it
A direct numerical comparison in which the quasinormal frequencies obtained from the reduced three-term continued-fraction equations diverge from those found by solving the original higher-order recurrence relations would falsify the scheme.
Figures
read the original abstract
Black-hole ringdown offers a clean probe of strong gravity, but one of its most accurate tools--Leaver's continued-fraction method--requires a three-term recurrence relation. Beyond general relativity, and more generally in non-Kerr spacetimes, Frobenius expansions of the perturbation equations generically produce higher-order recurrence relations and, often, couplings among the series coefficients, obstructing a direct application of Leaver's method. Here we develop a general reduction scheme that maps arbitrary scalar and matrix $N$-term recurrence relations to a three-term form, thereby extending continued fractions to a broad class of perturbation problems in modified gravity. As an application, we compute the quasinormal-mode spectrum of slowly-rotating black holes in dynamical Chern-Simons gravity, where the polar sector yields a 16-term, decoupled, scalar recurrence relation, and the axial sector yields a 12-term, coupled, matrix recurrence relation. After applying our reduction scheme, both systems can be solved with continued fractions. For the fundamental $(\ell,m)=(2,2)$ mode, our results agree well with independent calculations based on eigenvalue-perturbation and metric/spectral methods across the parameter range studied. This framework provides a robust and practical route to precision ringdown calculations beyond the standard three-term setting and supports tests of gravity with current and future gravitational-wave observations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a general reduction scheme to map arbitrary scalar and matrix N-term recurrence relations (arising from Frobenius expansions of perturbation equations) to equivalent three-term recurrences, thereby extending Leaver's continued-fraction method to quasinormal-mode problems in modified gravity. It applies the scheme to slowly rotating black holes in dynamical Chern-Simons gravity, reducing a 16-term decoupled scalar recurrence in the polar sector and a 12-term coupled matrix recurrence in the axial sector, and reports that the resulting continued-fraction solutions for the fundamental (ℓ,m)=(2,2) mode agree with independent eigenvalue-perturbation and spectral calculations across the studied parameter range.
Significance. If the reduction is shown to preserve exact equivalence and convergence properties, the work would meaningfully extend a high-precision tool to the broader class of non-Kerr and higher-order recurrence problems common in modified gravity, supporting future ringdown-based tests of gravity with LIGO/Virgo and next-generation detectors. The explicit treatment of both scalar and matrix cases, together with the reported numerical agreement for the (2,2) mode, provides a concrete demonstration of practicality.
major comments (2)
- [§3] §3 (reduction scheme): The mapping is asserted to preserve the exact quasinormal-mode condition for both scalar and matrix recurrences, yet the manuscript provides no explicit algebraic verification that the characteristic equation for the frequency remains unchanged after reduction (particularly for the 12-term axial matrix case). Without this, it is unclear whether the continued-fraction root-finding procedure solves precisely the same transcendental equation as the original N-term system.
- [§4.2] §4.2 (application to dCS gravity): Agreement is shown only for the fundamental (2,2) mode. The manuscript does not report convergence rates, the number of continued-fraction approximants required, or any comparison of the asymptotic growth of the new recurrence coefficients versus the original ones. This leaves the skeptic's concern unaddressed: the transformed coefficients are rational functions of the original ones and may alter the radius of convergence or numerical stability for overtones or other parameter values.
minor comments (2)
- [Abstract] The abstract states agreement 'across the parameter range studied' without specifying the range of the Chern-Simons coupling or spin; this should be stated explicitly in the text and figure captions.
- [§3] Notation for the reduced recurrence coefficients (e.g., the explicit form of the new a_n, b_n, c_n in terms of the original N-term coefficients) should be collected in a single table or appendix for readability.
Simulated Author's Rebuttal
We thank the referee for the positive summary and for identifying points that will strengthen the manuscript. We address the two major comments below and will revise accordingly.
read point-by-point responses
-
Referee: [§3] §3 (reduction scheme): The mapping is asserted to preserve the exact quasinormal-mode condition for both scalar and matrix recurrences, yet the manuscript provides no explicit algebraic verification that the characteristic equation for the frequency remains unchanged after reduction (particularly for the 12-term axial matrix case). Without this, it is unclear whether the continued-fraction root-finding procedure solves precisely the same transcendental equation as the original N-term system.
Authors: The reduction is constructed via a finite sequence of linear equivalence transformations (row operations and substitutions) that map the original N-term system onto an equivalent three-term system without altering its solution space or the associated characteristic equation for the frequency. For the scalar case this is immediate from the recursive elimination of higher-order coefficients. For the coupled matrix case the same holds because the block reductions preserve the determinant condition that defines the quasinormal-mode roots. We acknowledge that an explicit algebraic check was not written out in the original text. In the revised manuscript we will add a short appendix that carries out this verification step by step for both the 16-term scalar and 12-term matrix recurrences, confirming that the continued-fraction condition is identical to the original N-term condition. revision: yes
-
Referee: [§4.2] §4.2 (application to dCS gravity): Agreement is shown only for the fundamental (2,2) mode. The manuscript does not report convergence rates, the number of continued-fraction approximants required, or any comparison of the asymptotic growth of the new recurrence coefficients versus the original ones. This leaves the skeptic's concern unaddressed: the transformed coefficients are rational functions of the original ones and may alter the radius of convergence or numerical stability for overtones or other parameter values.
Authors: We agree that additional quantitative information on convergence would be valuable. The revised version will include (i) a table or figure showing the number of continued-fraction approximants needed to reach a prescribed tolerance for the (2,2) mode across the studied range of the Chern-Simons coupling, and (ii) a brief comparison of the leading asymptotic coefficients before and after reduction. Because the transformation is algebraic and finite, the dominant large-n behavior of the recurrence coefficients is unchanged; the radius of convergence of the underlying Frobenius series is therefore inherited. We will add a short paragraph addressing numerical stability. The present work focuses on the fundamental mode as a first demonstration; the scheme itself applies unchanged to overtones, but we have not yet performed those calculations and therefore do not claim results for them here. revision: partial
Circularity Check
No significant circularity; reduction scheme is a direct mathematical construction
full rationale
The paper's central contribution is an explicit algorithmic reduction that converts arbitrary N-term scalar or matrix recurrences into equivalent three-term recurrences. This mapping is defined by the paper's own transformation rules and is then applied to the 16-term polar and 12-term axial recurrences arising in dynamical Chern-Simons gravity. The resulting three-term coefficients are used to compute quasinormal modes via Leaver's continued-fraction method, with numerical results cross-checked against independent eigenvalue-perturbation and spectral methods. No load-bearing step reduces to a fitted parameter, a self-citation that itself lacks external verification, or a renaming of a known result; the derivation remains self-contained once the reduction algorithm is stated.
Axiom & Free-Parameter Ledger
axioms (1)
- standard math Standard properties of recurrence relations and Frobenius series expansions hold for the perturbation equations in modified gravity.
Forward citations
Cited by 1 Pith paper
-
Non-radial pulsations of gravitationally coupled two-fluid neutron stars in general relativity
A new relativistic linear perturbation formalism is derived for polar modes in gravitationally coupled two-fluid neutron stars, with numerical computation of representative f- and p-mode spectra classified by dominant...
Reference graph
Works this paper leans on
-
[1]
We have evaluated the fundamental (Table I) and the first overtone modes for the ℓ = 2 case (Table II)
Polar sector The QNM frequencies for the polar sector are shown below. We have evaluated the fundamental (Table I) and the first overtone modes for the ℓ = 2 case (Table II). For each value of the black hole spin, we study three different dCS coupling constant values. Observe that the m = 0 polar frequencies do not change with a or α. This happens because...
-
[2]
For each value of the black-hole spin, we once again study 3 different dCS coupling constant values
Axial sector Similarly to the polar case, we provide the QNM frequencies for the axial sector for the fundamental mode (Table III) and the first overtone modes (Table IV) in the ℓ = 2 case. For each value of the black-hole spin, we once again study 3 different dCS coupling constant values. a/M α/M 2 m=−2m=−1m= 0m= 1m= 2 0.00 0.00 .37367, .08896 .37367, .0...
-
[3]
B. P. Abbottet al.(LIGO Scientific, Virgo), Phys. Rev. Lett.116, 061102 (2016), arXiv:1602.03837 [gr-qc]
work page internal anchor Pith review arXiv 2016
-
[4]
Tests of General Relativity with GWTC-3
R. Abbottet al.(LIGO Scientific, VIRGO, KAGRA), Phys. Rev. D112, 084080 (2025), arXiv:2112.06861 [gr- qc]
work page internal anchor Pith review arXiv 2025
-
[5]
C. M. Will, Living Rev. Rel.17, 4 (2014), arXiv:1403.7377 [gr-qc]
work page internal anchor Pith review arXiv 2014
-
[6]
Yunes, X
N. Yunes, X. Siemens, and K. Yagi, Living Rev. Rel.28, 3 (2025)
2025
-
[7]
Theoretical Physics Implications of the Binary Black-Hole Mergers GW150914 and GW151226
N. Yunes, K. Yagi, and F. Pretorius, Phys. Rev. D94, 084002 (2016), arXiv:1603.08955 [gr-qc]
work page Pith review arXiv 2016
- [8]
- [9]
-
[10]
Birefringence tests of gravity with multi- messenger binaries,
M. Lagos, L. Jenks, M. Isi, K. Hotokezaka, B. D. Met- zger, E. Burns, W. M. Farr, S. Perkins, K. W. K. Wong, and N. Yunes, Phys. Rev. D109, 124003 (2024), 19 arXiv:2402.05316 [gr-qc]
-
[11]
Y. Xie, G. Narayan, and N. Yunes, (2025), arXiv:2510.02515 [gr-qc]
work page internal anchor Pith review arXiv 2025
-
[12]
Black Hole Based Tests of General Relativity
K. Yagi and L. C. Stein, Class. Quant. Grav.33, 054001 (2016), arXiv:1602.02413 [gr-qc]
work page Pith review arXiv 2016
-
[13]
Testing the black hole no-hair hypothesis,
V. Cardoso and L. Gualtieri, Class. Quant. Grav.33, 174001 (2016), arXiv:1607.03133 [gr-qc]
- [14]
-
[15]
Islam, (2021), arXiv:2111.00111 [gr-qc]
T. Islam, (2021), arXiv:2111.00111 [gr-qc]
- [16]
- [17]
- [18]
-
[19]
K. Schumacher, N. Yunes, and K. Yagi, Phys. Rev. D 108, 104038 (2023), arXiv:2308.05589 [gr-qc]
-
[20]
B. P. Abbottet al.(LIGO Scientific, Virgo), Phys. Rev. Lett.116, 221101 (2016), [Erratum: Phys.Rev.Lett. 121, 129902 (2018)], arXiv:1602.03841 [gr-qc]
work page Pith review arXiv 2016
- [21]
-
[22]
S. Bhagwat, X. J. Forteza, P. Pani, and V. Ferrari, Phys. Rev. D101, 044033 (2020), arXiv:1910.08708 [gr-qc]
- [23]
-
[24]
J. Calder´ on Bustillo, P. D. Lasky, and E. Thrane, Phys. Rev. D103, 024041 (2021), arXiv:2010.01857 [gr-qc]
-
[25]
Quasinormal modes of black holes and black branes
E. Berti, V. Cardoso, and A. O. Starinets, Class. Quant. Grav.26, 163001 (2009), arXiv:0905.2975 [gr-qc]
work page internal anchor Pith review arXiv 2009
- [26]
-
[27]
Gupta, (2025), arXiv:2511.15890 [gr-qc]
A. Gupta, (2025), arXiv:2511.15890 [gr-qc]
-
[28]
S. A. Teukolsky, Astrophys. J.185, 635 (1973)
1973
-
[29]
W. H. Press and S. A. Teukolsky, Astrophys. J.185, 649 (1973)
1973
-
[30]
E. W. Leaver, Proc. Roy. Soc. Lond. A402, 285 (1985)
1985
-
[31]
Testing General Relativity with Present and Future Astrophysical Observations
E. Bertiet al., Class. Quant. Grav.32, 243001 (2015), arXiv:1501.07274 [gr-qc]
work page internal anchor Pith review arXiv 2015
-
[32]
S. Alexander and N. Yunes, Phys. Rept.480, 1 (2009), arXiv:0907.2562 [hep-th]
-
[33]
A. Maselli, S. Yi, L. Pierini, V. Vellucci, L. Reali, L. Gualtieri, and E. Berti, Phys. Rev. D109, 064060 (2024), arXiv:2311.14803 [gr-qc]
-
[34]
Extreme Gravity Tests with Gravitational Waves from Compact Binary Coalescences: (I) Inspiral-Merger
E. Berti, K. Yagi, and N. Yunes, Gen. Rel. Grav.50, 46 (2018), arXiv:1801.03208 [gr-qc]
work page Pith review arXiv 2018
- [35]
- [36]
- [37]
- [38]
- [39]
-
[40]
A. Hussain and A. Zimmerman, Phys. Rev. D106, 104018 (2022), arXiv:2206.10653 [gr-qc]
- [41]
- [42]
- [43]
- [44]
- [45]
- [46]
- [47]
-
[48]
M. Srivastava, Y. Chen, and S. Shankaranarayanan, Phys. Rev. D104, 064034 (2021), arXiv:2106.06209 [gr-qc]
-
[49]
T. Alapati and S. Shankaranarayanan, Phys. Rev. D112, 024053 (2025), arXiv:2506.03600 [gr-qc]
- [50]
-
[51]
S. J. B. Husken, T. van der Steen, S. Maenaut, K. K.-H. Lam, M. D. Jockwer, A. K.-W. Chung, T. Hertog, T. G. F. Li, and N. Yunes, (2026), arXiv:2604.02214 [gr-qc]
work page internal anchor Pith review Pith/arXiv arXiv 2026
- [52]
- [53]
-
[54]
L. Pierini and L. Gualtieri, Phys. Rev. D103, 124017 (2021), arXiv:2103.09870 [gr-qc]
-
[55]
L. Pierini and L. Gualtieri, Phys. Rev. D106, 104009 (2022), arXiv:2207.11267 [gr-qc]
-
[56]
Maggiore,Gravitational Waves
M. Maggiore,Gravitational Waves. Vol. 2: Astrophysics and Cosmology(Oxford University Press, 2018)
2018
-
[57]
Chandrasekhar,The mathematical theory of black holes (1985)
S. Chandrasekhar,The mathematical theory of black holes (1985)
1985
-
[58]
Regge and J
T. Regge and J. A. Wheeler, Phys. Rev.108, 1063 (1957)
1957
-
[59]
F. J. Zerilli, Phys. Rev. Lett.24, 737 (1970)
1970
-
[60]
Moncrief, Annals Phys.88, 323 (1974)
V. Moncrief, Annals Phys.88, 323 (1974)
1974
- [61]
- [62]
- [63]
-
[64]
H. Roussille, D. Langlois, and K. Noui, JCAP01, 054, arXiv:2312.11986 [gr-qc]
- [65]
-
[66]
E. W. Leaver, Phys. Rev. D41, 2986 (1990)
1990
-
[67]
Quasinormal modes of maximally charged black holes
H. Onozawa, T. Mishima, T. Okamura, and H. Ishihara, Phys. Rev. D53, 7033 (1996), arXiv:gr-qc/9603021
work page Pith review arXiv 1996
- [68]
- [69]
-
[70]
Pani, Advanced Methods in Black-Hole Perturba- tion Theory , Int
P. Pani, Int. J. Mod. Phys. A28, 1340018 (2013), arXiv:1305.6759 [gr-qc]
-
[71]
R. Jackiw and S. Y. Pi, Phys. Rev. D68, 104012 (2003), arXiv:gr-qc/0308071
-
[72]
Karikos, J
G. Karikos, J. A. Saes, and P. Wagle, Matrix- valued-continued-fractions-method-for-dcs-gravity, https://github.com/Giorgos-Karikos/Matrix- 20 valued-continued-fractions-method-for-dCS- gravity/tree/main(2026)
2026
-
[73]
C. W. Misner, K. S. Thorne, and J. A. Wheeler,Gravita- tion(W. H. Freeman, San Francisco, 1973)
1973
-
[74]
A. Zimmerman and Z. Mark, Phys. Rev. D93, 044033 (2016), [Erratum: Phys.Rev.D 93, 089905 (2016)], arXiv:1512.02247 [gr-qc]
-
[75]
S. Alexander, G. Gabadadze, L. Jenks, and N. Yunes, Phys. Rev. D104, 064033 (2021), arXiv:2104.00019 [hep- th]
-
[76]
T. Delsate, D. Hilditch, and H. Witek, Phys. Rev. D91, 024027 (2015), arXiv:1407.6727 [gr-qc]
-
[77]
H. Motohashi and T. Suyama, Phys. Rev. D84, 084041 (2011), arXiv:1107.3705 [gr-qc]
-
[78]
H. Motohashi and T. Suyama, Phys. Rev. D85, 044054 (2012), arXiv:1110.6241 [gr-qc]
-
[79]
S. Alexander, H. Bernardo, and N. Yunes, Phys. Rev. D 112, 104042 (2025), arXiv:2506.14889 [hep-th]
-
[80]
Dynamical Chern-Simons Modified Gravity I: Spinning Black Holes in the Slow-Rotation Approximation
N. Yunes and F. Pretorius, Phys. Rev. D79, 084043 (2009), arXiv:0902.4669 [gr-qc]
work page Pith review arXiv 2009
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.