Modified Teukolsky Equation: Extensions and Insights

• Modified Teukolsky Equation is a family of generalizations of the classic perturbation equation that models spin-weighted fields on Kerr and Schwarzschild backgrounds with additional mass, dispersion, and curvature terms.
• It employs techniques such as potential term additions, parameterized deformations, conformal methods, and Green-function approaches to reorganize the traditional Teukolsky toolkit while preserving separability and curvature-based structures.
• This formulation facilitates precise modeling of deviations from general relativity in quasinormal mode analysis and inspiral waveforms, offering actionable insights into beyond-GR black-hole physics.

Searching arXiv for recent and foundational papers on the modified Teukolsky equation and closely related formalisms. The modified Teukolsky equation is a family of generalizations of the Teukolsky master equation for perturbations of algebraically special black-hole spacetimes. In its classical form, the Teukolsky equation governs spin-weighted fields on Kerr or Schwarzschild backgrounds and is separable into radial and angular equations. In the literature, the term “modified Teukolsky equation” is used for several distinct constructions: phenomenological additions of dispersion or mass terms, decoupled curvature equations on non-Ricci-flat backgrounds, parameterized deformations of the radial and angular sectors, source-driven beyond-GR formalisms for inspirals, and conformally regular versions adapted to hyperboloidal or compactified settings (Chung et al., 2018, Guo et al., 2024, Yu et al., 4 Nov 2025, Gasperin et al., 22 Feb 2026).

1. Classical framework and motivation

The standard Teukolsky framework describes perturbations of a Kerr black hole of mass $M$ and angular momentum $J=aM$ through spin-weighted fields $\psi_s$, with $s=0,\pm1,\pm2$, and the gravitationally relevant outgoing sector at $s=-2$. In vacuum GR, the equation is massless, separable, and admits quasinormal-mode boundary conditions that are purely ingoing at the horizon and purely outgoing at infinity. The separated form uses a radial function, spin-weighted spheroidal harmonics, and a complex frequency $\tilde\omega_{nlm}$ selected by the QNM boundary value problem (Chung et al., 2018).

The need for modification arises whenever one wishes to encode physics not present in vacuum GR. The examples developed in the literature include nonzero graviton mass and alternative dispersion relations, non-Ricci-flat effective geometries, higher-curvature EFT corrections, parity-violating couplings, and formulations in which the conformal factor is treated as a dynamical variable. A recurring design requirement is to preserve as much of the Teukolsky machinery as possible: curvature-based variables, separability, controllable asymptotics, and compatibility with QNM or flux calculations (Guo et al., 2024, Yu et al., 4 Nov 2025, Yang et al., 8 Jun 2026).

2. Principal forms of modification

A first class of modifications adds an explicit potential-like term to the Teukolsky operator. In the massive-graviton phenomenology of “Phenomenological Inclusion of Alternative Dispersion Relations to the Teukolsky Equation and its Application to Bounding the Graviton Mass with Gravitational-wave Measurements” (Chung et al., 2018), the modified equation is

$\mathcal{L}\psi + \mathcal{D}\Sigma\psi = 0,$

and for a massive field $\mathcal{D}=m^2$,

$\mathcal{L}\psi + m^2\Sigma\psi = 0.$

After separation, this produces a radial equation with an additional effective potential

$V^{(1)}(r)=m_g^2\,\frac{r^2\Delta}{(r^2+a^2)^2},$

and the leading QNM shift

$J=aM$0

This construction is explicitly phenomenological: it reproduces the desired asymptotic dispersion while retaining separability, but it does not attempt to implement a complete massive-gravity theory (Chung et al., 2018).

A second class keeps the Kerr background fixed and introduces small, theory-agnostic deformations in both the radial and angular sectors. In “The parameterized quasinormal modes for modified Teukolsky equations” (Yu et al., 4 Nov 2025), the modified master equation is

$J=aM$1

with separability enforced by

$J=aM$2

The result is a radial equation with additive deformation $J=aM$3 and an angular equation with additive deformation $J=aM$4, parameterized by two independent sets of coefficients $J=aM$5 and $J=aM$6 (Yu et al., 4 Nov 2025).

A third class modifies the curvature operator itself to accommodate non-Ricci-flat backgrounds. In “Generic Modified Teukolsky Formalism beyond General Relativity for Spherically Symmetric Cases” (Guo et al., 2024), the key decoupled equation for $J=aM$7 contains an added curvature term $J=aM$8 relative to the GR operator,

$J=aM$9

with an analogous equation for $\psi_s$0. For the $\psi_s$1-symmetric static spherically symmetric subclass, this yields a separated master equation whose extra effective potential is determined by $\psi_s$2 and its derivatives, and which reduces to the standard Schwarzschild Teukolsky equation when $\psi_s$3 (Guo et al., 2024).

A fourth class is covariant rather than separability-driven. “A Fundamental Equation for Gravitational Wave and Its Analogues in Type D Spacetimes” (Li, 2020) introduces the spin-coefficient connection $\psi_s$4 and a unified equation

$\psi_s$5

for massless fields with nonzero spin $\psi_s$6 on type D backgrounds. This is best understood as a generalized Teukolsky equation written in a coordinate-free NP-based form (Li, 2020).

A fifth class is conformal. “The linearised conformal Einstein field equations around a Petrov-type D spacetime: the conformal Teukolsky equation” (Gasperin et al., 22 Feb 2026) derives decoupled equations for NP components $\psi_s$7 and $\psi_s$8 of the rescaled Weyl tensor in Friedrich’s CEFE framework. The resulting operator preserves the structural form of the classical Teukolsky equation while remaining regular at the conformal boundary, and it supplies a geometric interpretation of hyperboloidal master variables (Gasperin et al., 22 Feb 2026).

Construction	Representative equation	Representative paper
Alternative dispersion / massive graviton	$\psi_s$9	(Chung et al., 2018)
Parameterized Kerr deformation	$s=0,\pm1,\pm2$0	(Yu et al., 4 Nov 2025)
Non-Ricci-flat SSS formalism	GR Teukolsky operator $s=0,\pm1,\pm2$1	(Guo et al., 2024)
Higher-derivative inspiral formalism	$s=0,\pm1,\pm2$2	(Yang et al., 8 Jun 2026)
Conformal formulation	Conformal Teukolsky equations for $s=0,\pm1,\pm2$3	(Gasperin et al., 22 Feb 2026)

3. Beyond-GR black holes, parity structure, and observables

A major use of modified Teukolsky formalisms is to compute quantities that are inaccessible in vacuum GR Teukolsky theory. In “Love numbers beyond GR from the modified Teukolsky equation” (Cano, 27 Feb 2025), the static spin-$s=0,\pm1,\pm2$4 radial equation is brought to the form

$s=0,\pm1,\pm2$5

which encodes EFT corrections through $s=0,\pm1,\pm2$6 and $s=0,\pm1,\pm2$7. This framework yields electric, magnetic, and “mixing” tidal Love numbers, the last associated to parity-breaking theories and identified there for the first time. It also gives closed-form expressions for running beta functions $s=0,\pm1,\pm2$8 and constant Love numbers $s=0,\pm1,\pm2$9 in the cases $s=-2$0 (Cano, 27 Feb 2025).

Parity structure is treated systematically in “Isospectrality breaking in the Teukolsky formalism” (Li et al., 2023). There, definite-parity modes are defined through parity-conjugated combinations of Weyl scalars, and the modified Teukolsky equation is written schematically as a driven system

$s=-2$1

Using eigenvalue perturbation theory, the paper shows how the degeneracy between even- and odd-parity QNM frequencies is broken in modified gravity, reducing the problem to a $s=-2$2 eigenvalue problem for the first-order frequency shifts (Li et al., 2023).

The same strategy extends to inspirals in non-GR backgrounds. In “Modified Teukolsky Formalism for Extreme Mass-Ratio Inspirals in Higher-Derivative Gravity” (Yang et al., 8 Jun 2026), the formalism is organized as a two-parameter expansion in EFT coupling $s=-2$3 and mass ratio $s=-2$4,

$s=-2$5

and the leading beyond-GR correction satisfies an inhomogeneous GR Teukolsky equation. In a cubic-gravity example around a non-rotating black hole, this is used to compute fluxes to the horizon and null infinity, with the horizon flux enhanced by roughly an order of magnitude relative to GR after factoring out $s=-2$6, while the flux at infinity is slightly reduced (Yang et al., 8 Jun 2026).

A related non-Ricci-flat application appears in “Extreme mass-ratio inspirals and extra dimensions: Insights from modified Teukolsky framework” (Kumar et al., 4 Jul 2025), where a braneworld black hole with tidal charge $s=-2$7 leads to a radial equation

$s=-2$8

with the MTE correction controlled by $s=-2$9. The paper reports significant deviations in dephasing and mismatch values, especially for higher eccentricities, and argues that using the MTE rather than a GR Teukolsky treatment is necessary for accurately modeling such binaries (Kumar et al., 4 Jul 2025).

4. Solution methods and numerical realization

Because the various modified equations preserve different subsets of the Teukolsky structure, the preferred solution method depends on the modification. For the massive-graviton construction, the extra radial potential is treated perturbatively by logarithmic perturbation theory, which yields the leading QNM shift directly from the unperturbed Kerr mode functions and explains why the correction is quadratic in $\tilde\omega_{nlm}$0 (Chung et al., 2018).

For parameterized Kerr deformations, continued-fraction technology remains central. The parameterized framework of (Yu et al., 4 Nov 2025) converts the modified angular and radial recurrences into finite-band matrix problems,

$\tilde\omega_{nlm}$1

and linearizes the spectrum as

$\tilde\omega_{nlm}$2

Those linear coefficients are cross-validated against a full $\tilde\omega_{nlm}$3D hyperboloidal pseudo-spectral solution of the modified master equation, with the two approaches agreeing at the $\tilde\omega_{nlm}$4 level over a broad small-deviation regime (Yu et al., 4 Nov 2025).

For source-driven EMRI problems, Green-function methods and metric reconstruction remain effective once the modified source terms are known. In the cubic-gravity EMRI formalism, the radial Teukolsky solution is written through homogeneous GR solutions and a modified source, and the near-horizon singular structure is handled by Taylor expansion plus lower incomplete Gamma functions; the implementation is stabilized numerically by Chebyshev fitting of the GR metric perturbations and recursive differentiation (Yang et al., 8 Jun 2026).

This suggests a broad methodological principle: the most successful modified formalisms do not abandon the classical Teukolsky toolkit, but reorganize it. Separability, continued fractions, Chandrasekhar-type transformations, Green functions, and hyperboloidal compactification reappear in altered roles rather than disappearing altogether. That pattern is explicit in the parameterized, inspiral, and conformal programs (Yu et al., 4 Nov 2025, Gasperin et al., 22 Feb 2026).

5. Terminological distinctions and adjacent constructions

The phrase “modified Teukolsky equation” is narrower than “modified Teukolsky setup,” and several nearby works make that distinction explicit. “Uniform Boundedness for Solutions to the Teukolsky Equation on Schwarzschild from Conservation Laws of Linearised Gravity” (Collingbourne et al., 2023) does not alter the Teukolsky PDE itself; instead, it derives outgoing-cone control of the gauge-invariant Teukolsky quantities from conservation laws of the full linearised Einstein system, using the standard spin $\tilde\omega_{nlm}$5 Teukolsky equations only at a later stage (Collingbourne et al., 2023).

Likewise, “Gravitational-wave echoes from spinning exotic compact objects: numerical waveforms from the Teukolsky equation” (Xin et al., 2021) keeps the Kerr–Teukolsky operator unchanged. What is modified there are the source, via phenomenological trajectories, and the boundary conditions, via a partially reflecting ECO surface related to $\tilde\omega_{nlm}$6 and $\tilde\omega_{nlm}$7 through membrane-paradigm considerations (Xin et al., 2021).

“A Scattering Theory for Linearised Gravity on the Exterior of the Schwarzschild Black Hole I: The Teukolsky Equations” (Masaood, 2020) develops what may be called a modified-Teukolsky viewpoint in a different sense: the physical-space Chandrasekhar transform sends the spin $\tilde\omega_{nlm}$8 Teukolsky equations to Regge–Wheeler equations with real potentials, thereby producing a scattering theory through transformed variables rather than through a new operator on $\tilde\omega_{nlm}$9 or $\mathcal{L}\psi + \mathcal{D}\Sigma\psi = 0,$0 themselves (Masaood, 2020).

Finally, “An alternative to the Teukolsky equation” (Hatsuda, 2020) proposes a new ordinary differential equation conjectured to be exactly isospectral to the radial homogeneous Teukolsky equation. This is not a modification of the Teukolsky operator in the usual beyond-GR sense, but an alternative ODE with the same QNM spectrum, motivated by a hidden symmetry from four-dimensional $\mathcal{L}\psi + \mathcal{D}\Sigma\psi = 0,$1 supersymmetric quantum chromodynamics (Hatsuda, 2020).

6. Research directions and conceptual significance

Several directions recur across the literature. One is to move from model-specific modifications to theory-agnostic parameterizations. The parameterized Kerr program explicitly aims to map ringdown residuals onto a finite set of response coefficients, while the EFT Love-number program extracts observables directly from coefficients in the modified potential (Yu et al., 4 Nov 2025, Cano, 27 Feb 2025). Another is to extend non-rotating results to rotating black holes: this is emphasized in the higher-derivative EMRI formalism, the Love-number computation, and the isospectrality-breaking analysis (Yang et al., 8 Jun 2026, Cano, 27 Feb 2025, Li et al., 2023).

A second direction is to combine strong-field generation effects with propagation effects. The massive-graviton study emphasizes that ringdown probes the strong-field generation of gravitational waves near the horizon, whereas inspiral bounds on the graviton mass are dominated by weak-field propagation over cosmological distances. A plausible implication is that consistent inference in modified gravity will increasingly rely on joint treatments of inspiral, plunge, and ringdown rather than on isolated sectors (Chung et al., 2018).

A third direction is geometric regularization. The conformal Teukolsky equation shows that hyperboloidal master variables can be interpreted as NP components of the rescaled Weyl tensor in the CEFE framework, and the generic SSS formalism shows that one can retain decoupled curvature equations on non-Ricci-flat backgrounds while avoiding the “fake poles” associated with some earlier constructions (Gasperin et al., 22 Feb 2026, Guo et al., 2024).

Taken together, these developments establish the modified Teukolsky equation not as a single canonical PDE but as a research program. Its unifying aim is to preserve the Teukolsky formalism’s curvature-based and often separable structure while extending it to alternative dispersion, non-vacuum effective geometries, higher-curvature corrections, parity-violating dynamics, conformal compactifications, and waveform-generation problems beyond general relativity (Chung et al., 2018, Guo et al., 2024, Yu et al., 4 Nov 2025, Yang et al., 8 Jun 2026).