Beyond Three Terms: Continued Fractions for Rotating Black Holes in Modified Gravity

TL;DR

This paper extends Leaver's continued-fraction method to handle higher-order recurrence relations in modified gravity black hole perturbations, enabling precise quasinormal mode calculations for non-Kerr spacetimes.

Contribution

We develop a reduction scheme that transforms arbitrary N-term recurrence relations into a three-term form, broadening the applicability of continued fractions in modified gravity perturbation analysis.

Findings

Successfully applied the method to compute quasinormal modes in dynamical Chern-Simons gravity.

Achieved good agreement with independent methods for the fundamental mode.

Extended continued fraction techniques to complex, coupled recurrence relations in non-Kerr spacetimes.

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.