Teukolsky Equations in Kerr and Schwarzschild

• Teukolsky Equations are decoupled differential equations that describe massless field perturbations on algebraically special black-hole spacetimes like Kerr and Schwarzschild.
• They separate into radial and angular parts, facilitating the analysis of spin-weighted fields and the computation of quasinormal modes and stability measures.
• Extensions include non-vacuum, Kerr–Newman, and higher-derivative gravity cases, underscoring their versatility in modern gravitational research.

The Teukolsky equations are decoupled equations for extreme-spin components of massless fields on algebraically special black-hole backgrounds, most prominently Kerr and Schwarzschild. In the Kerr setting, the Teukolsky Master Equation governs linear perturbations of scalar, electromagnetic, and gravitational fields and separates into radial and angular ordinary differential equations for spin weight $s$; in Schwarzschild, tensorial spin $\pm1$ and $\pm2$ formulations arise for the extreme null Maxwell and Weyl components (Cano et al., 2024, Masaood, 2020). Across the modern literature, the term also encompasses adjoint equations for Hertz potentials, physical-space formulations tied to Chandrasekhar transformations, and a broad family of Teukolsky-like generalizations on non-vacuum or deformed type D spacetimes (Green et al., 2019, Guo et al., 2023).

1. Definition, field content, and geometric setting

In Kerr spacetime, the Teukolsky formalism gives a separable description of linear perturbations of Kerr for scalar, electromagnetic, and gravitational fields (Cano et al., 2024). In Newman–Penrose language, the radiative variables are the extreme helicity components: for electromagnetism the extreme Maxwell scalars $\phi_0,\phi_2$, and for gravity the extreme Weyl scalars $\psi_0,\psi_4$ or their tensorial counterparts (Giorgi, 2020). In Schwarzschild, the spin $\pm1$ Teukolsky equations are the equations satisfied by the extreme components $\alpha$ and $\underline\alpha$ of the Maxwell field, when expressed with respect to a null frame (Pasqualotto, 2016). For linearized gravity on Schwarzschild, the corresponding extreme curvature components are the $S^2$-tangent symmetric traceless $2$-tensors

$\pm1$0

which satisfy the spin $\pm1$1 and spin $\pm1$2 Teukolsky equations (Masaood, 2020).

A common structural prerequisite is Petrov type D geometry. A general spherically symmetric spacetime is still Petrov type D, so much of the Newman–Penrose/Teukolsky machinery survives even without the vacuum condition (Guo et al., 2023). Likewise, an axisymmetric metric satisfying the Petrov type D property with additional ansätze can be brought to a Kerr-like form deformed by one radial function $\pm1$3, and in the $\pm1$4 sector the homogeneous wave equations admit separation of variables (Guo et al., 2023). This type-D background structure is the main reason that decoupled equations for extreme-spin fields continue to exist well beyond exact Kerr vacuum.

The geometric interpretation has also been reformulated in terms of a preferred weighted connection. On algebraically special conformal spacetimes admitting a shear-free null geodesic congruence, a $\pm1$5-complex-dimensional twistor manifold induces a weighted covariant derivative $\pm1$6, and the scalar wave operators built from $\pm1$7 reproduce the Teukolsky operators for massless fields (Araneda, 2019). This identifies the Teukolsky equations as wave equations with respect to a twistor-induced connection rather than only as Newman–Penrose identities.

2. Separation on Kerr and the master equation

In Boyer–Lindquist coordinates, the Teukolsky Master Equation for a massless field $\pm1$8 of spin weight $\pm1$9 on Kerr is written as

$\pm2$0

with

$\pm2$1

(0903.3617). The separated ansatz

$\pm2$2

reduces the problem to the angular Teukolsky equation

$\pm2$3

and the radial Teukolsky equation

$\pm2$4

where

$\pm2$5

(0903.3617).

The angular equation is the equation for spin-weighted spheroidal harmonics. In the conventions used in the parametrized quasi-normal-mode literature,

$\pm2$6

while the radial Teukolsky equation is written as

$\pm2$7

with

$\pm2$8

(Cano et al., 2024). Quasinormal modes are separated solutions whose radial part is purely ingoing at the event horizon and purely outgoing at infinity, while the angular part is regular on the sphere (Cano et al., 2024).

Several modern developments preserve this separability structure while modifying the operator. A parameterized Kerr framework deforms the radial potential by

$\pm2$9

leading to linear shifts

$\phi_0,\phi_2$0

(Cano et al., 2024). A later parameterized framework modifies both radial and angular equations through separable deformations $\phi_0,\phi_2$1 and $\phi_0,\phi_2$2, keeping the same Kerr background but altering the master equation itself (Yu et al., 4 Nov 2025). These constructions retain Kerr separability as the organizing principle.

3. Spin reversal, Teukolsky–Starobinsky identities, and exact structures

A central internal symmetry of separated Teukolsky theory is spin reversal. Under separation of variables, spin-reversal for the Teukolsky Master Equation is accomplished through the so-called Teukolsky–Starobinsky identities, with associated angular and radial Teukolsky–Starobinsky constants (Costa et al., 2021). For the separated angular and radial ODEs, repeated application of first-order operators maps spin $\phi_0,\phi_2$3 solutions to spin $\phi_0,\phi_2$4 solutions and conversely (Costa et al., 2021). If the Teukolsky–Starobinsky constant is nonzero, this correspondence is invertible up to multiplication by the constant.

For physical perturbation theory, these identities are indispensable because they tie together the Hertz potential and the observable extreme Newman–Penrose scalar. In canonical quantization on Kerr, the basic pair is

$\phi_0,\phi_2$5

where $\phi_0,\phi_2$6 is a Hertz potential and $\phi_0,\phi_2$7 is an NP scalar satisfying the spin-$\phi_0,\phi_2$8 Teukolsky equation (Iuliano et al., 2023). The canonical commutation relations can be implemented if and only if the Teukolsky–Starobinsky constants are positive, which is the case both for gravitational perturbations and Maxwell fields (Iuliano et al., 2023).

The sign of the radial Teukolsky–Starobinsky constant is therefore not a minor technicality. For spins $\phi_0,\phi_2$9, the radial constant is strictly positive for real separated Kerr modes (Costa et al., 2021). However, contrary to popular belief, these constants can be negative for spin larger than $\psi_0,\psi_4$0 (Costa et al., 2021). The common argument for universal nonnegativity fails because the radial spin-reversal operators are related by complex conjugation, not by an adjointness relation of the type available in the angular problem (Costa et al., 2021). The consequence is a novel form of energy amplification which occurs for non-superradiant frequencies when $\psi_0,\psi_4$1 (Costa et al., 2021). That result corrects a persistent misconception in the separated theory.

The exact analytic structure of the separated equations has also been studied through confluent Heun theory. Both the angular and radial Teukolsky equations on Kerr can be transformed into confluent Heun equations (0903.3617). For $\psi_0,\psi_4$2 and $\psi_0,\psi_4$3, there exists a class of simultaneous angular/radial polynomial solutions in which the polynomiality conditions for the two separated equations coincide, determining the separation constant $\psi_0,\psi_4$4 as a function of $\psi_0,\psi_4$5 while leaving the complex frequency unconstrained; the result is a continuous spectrum in the sense of that special exact family (0903.3617). The paper is explicit that this is not the generic Kerr quasinormal-mode problem but a special exact family of one-way-wave solutions (0903.3617).

4. Physical-space formulations on Schwarzschild

On Schwarzschild, the Teukolsky equations have been reformulated in fully tensorial, physical-space terms. For the Maxwell system, the exterior Schwarzschild spacetime

$\psi_0,\psi_4$6

with

$\psi_0,\psi_4$7

supports a null decomposition of a Maxwell field $\psi_0,\psi_4$8 into extreme components $\psi_0,\psi_4$9 and middle components $\pm1$0 (Pasqualotto, 2016). From the null Maxwell system one derives the tensorial spin $\pm1$1 Teukolsky equations for $\pm1$2 and $\pm1$3: $\pm1$4

$\pm1$5

(Pasqualotto, 2016).

A key structural point is that these equations are not directly amenable to Morawetz estimates because of the bad first-order term

$\pm1$6

whose sign changes at the photon sphere $\pm1$7 (Pasqualotto, 2016). Pasqualotto’s solution is a Chandrasekhar-type differential transformation

$\pm1$8

which converts the Teukolsky variables into fields satisfying the tensorial Fackerell–Ipser equation

$\pm1$9

and similarly for $\alpha$0 (Pasqualotto, 2016). This transformed equation admits standard redshift, energy, Morawetz, and $\alpha$1-weighted estimates, and the transformed quantities vanish on the $\alpha$2 Coulomb sector (Pasqualotto, 2016). The same paper proves inverse polynomial decay for the transformed fields and transfers it back to $\alpha$3, obtaining quantitative decay for the full Maxwell system modulo the stationary Coulomb tail (Pasqualotto, 2016).

For gravitational perturbations on Schwarzschild, an analogous strategy underlies a full scattering theory for the spin $\alpha$4 Teukolsky equations (Masaood, 2020). The spin $\alpha$5 and spin $\alpha$6 equations for $\alpha$7 and $\alpha$8 contain first-order null derivative terms that obstruct a direct Lagrangian energy method (Masaood, 2020). A physical-space Chandrasekhar transformation converts the Teukolsky fields into Regge–Wheeler variables satisfying

$\alpha$9

and the Regge–Wheeler scattering theory is then transplanted back to the Teukolsky fields (Masaood, 2020). The resulting forward and backward maps are unitary Hilbert-space isomorphisms between Cauchy data spaces and radiation-field spaces for both spin sectors (Masaood, 2020). The same paper also treats the Teukolsky–Starobinsky correspondence at the level of scattering data and derives a mixed scattering theorem that isolates the radiative degrees of freedom relevant for linearised Einstein scattering on Schwarzschild (Masaood, 2020).

These Schwarzschild results are important because they show that, even when direct multiplier methods fail at the Teukolsky level, differential transformations can restore robust physical-space estimates and even asymptotic completeness.

5. Extensions beyond vacuum Kerr

A major direction of recent research is the extension of Teukolsky theory beyond vacuum Kerr. On a general spherically symmetric spacetime

$\underline\alpha$0

one obtains a unified Teukolsky-like master equation for massless fields of spin

$\underline\alpha$1

(Guo et al., 2023). After a rescaling $\underline\alpha$2, the unified equation becomes

$\underline\alpha$3

The extra curvature terms involving $\underline\alpha$4 and $\underline\alpha$5 are the clean mathematical signature of the non-vacuum extension (Guo et al., 2023).

For axisymmetric type D non-vacuum backgrounds, a specific class of metrics becomes conformal to Kerr deformed by one radial function $\underline\alpha$6, and in the $\underline\alpha$7 sector the perturbed Weyl-scalar equations admit separation under an explicitly chosen separable tetrad gauge (Guo et al., 2023). The resulting generalized spin $\underline\alpha$8 radial equation is

$\underline\alpha$9

which is literally the Kerr spin-$S^2$0 radial Teukolsky equation with $S^2$1 replaced by a general radial function $S^2$2 (Guo et al., 2023).

That generalization has been extended to all spins $S^2$3 in a deformed Kerr spacetime with metric

$S^2$4

(Nakajima et al., 2024). After a $S^2$5 rescaling and, for spin $S^2$6, the same separable gauge, the unified master equation takes the compact form

$S^2$7

The angular equation remains exactly the spin-weighted spheroidal-harmonic equation of Kerr, while the radial equation is modified by $S^2$8, $S^2$9, and $2$0 (Nakajima et al., 2024). The paper also analyzes the singular-point structure of the generalized radial ODE, showing that simple zeros or poles of $2$1 become regular singular points and $2$2 remains an irregular singular point of rank $2$3 (Nakajima et al., 2024).

Higher-derivative gravity provides another route to generalized Teukolsky equations. The so-called universal Teukolsky equations are identities valid on any four-dimensional spacetime in any theory; theory dependence enters through the Ricci tensor or effective stress tensor on the right-hand side (Cano et al., 2023). When specialized perturbatively to rotating black holes close to Kerr in higher-derivative gravity, they reduce to four coupled equations for $2$4 and their conjugates, which can then be decoupled perturbatively into modified second-order radial equations (Cano et al., 2023). This framework has been applied to six-derivative gravity to compute quasinormal-mode shifts for spinning black holes (Cano et al., 2023).

A further important extension is Kerr–Newman. There, the difficulty is the coupling between spin-$2$5 electromagnetic and spin-$2$6 gravitational perturbations. Because spin-weighted spheroidal harmonics of different spins are not related in the right way by the relevant angular operators, standard mode separation breaks down (Giorgi, 2020). The physical-space resolution is not a single decoupled scalar equation but a coupled hierarchy of gauge-invariant tensorial Teukolsky-type and Regge–Wheeler-type equations suited to multiplier methods (Giorgi, 2020).

6. Sources, nonlinear perturbations, and quantum theory

The classical separability of the homogeneous Teukolsky operator does not by itself solve the sourced problem. For generic extended sources in Kerr, the right-hand side is not naturally separated because background quantities such as

$2$7

mix radial and angular dependence (Spiers, 2024). A recent analytic framework addresses this by decomposing the source into the angular basis adapted to separability, using Newman–Penrose, GHP, Held formalism, spin-weighted spherical harmonic recoupling, and a truncated Fourier expansion of $2$8 (Spiers, 2024). In a proof-of-concept example for the spin $2$9 equation sourced by the stress-energy tensor of a stationary ideal gas cloud around Kerr, the source is analytically separated into spin-weighted harmonic modes, and for $\pm1$00, $\pm1$01, and Fourier order $\pm1$02, the error in the $\pm1$03 approximation is below $\pm1$04 near the horizon (Spiers, 2024). This makes sourced Teukolsky calculations more tractable in second-order perturbation theory, self-force problems, and nonlinear quasinormal-mode studies (Spiers, 2024).

A complementary nonlinear framework treats higher-order Kerr perturbations by combining a sourced adjoint spin $\pm1$05 Teukolsky equation with transport equations along principal null geodesics (Green et al., 2019). A generic sourced metric perturbation is decomposed as

$\pm1$06

where $\pm1$07 is pure gauge, $\pm1$08 is a Kerr zero mode, $\pm1$09 is a corrector tensor determined by transport equations, and $\pm1$10 satisfies the sourced adjoint Teukolsky equation

$\pm1$11

(Green et al., 2019). In outgoing Kerr–Newman coordinates, the relevant auxiliary equations reduce to ordinary differential equations in $\pm1$12, so the hard part of the sourced Einstein problem is again a separable scalar Teukolsky equation (Green et al., 2019). This same decomposition iterates order by order in nonlinear perturbation theory (Green et al., 2019).

At the quantum level, the Teukolsky formalism supports a canonical quantization for integer spins $\pm1$13 on Kerr (Iuliano et al., 2023). The Teukolsky operator can be written in GHP-covariant form as

$\pm1$14

while the adjoint equation for the Hertz potential is $\pm1$15 (Iuliano et al., 2023). The corresponding Teukolsky action yields a conserved symplectic form, and the Hadamard parametrix for the Teukolsky operator can be constructed by exploiting its charged Klein–Gordon form (Iuliano et al., 2023). The same work proves that the canonical energy of Teukolsky fields equals twice the canonical energy of the reconstructed gravitational perturbation for compactly supported data (Iuliano et al., 2023).

7. Conceptual themes, controversies, and significance

Several recurring themes unify this body of work. First, the Teukolsky equations are powerful because they isolate extreme-spin radiative degrees of freedom in a form compatible with separation of variables or, failing that, with physical-space transformations (Cano et al., 2024, Pasqualotto, 2016). Second, the price of this decoupling is that the resulting equations often carry unfavorable first-order terms, gauge subtleties, or couplings to lower-spin sectors. The Schwarzschild Maxwell and gravitational analyses make explicit that direct Morawetz estimates at the Teukolsky level can fail and that Chandrasekhar-type transformations are the decisive remedy (Pasqualotto, 2016, Masaood, 2020).

Third, several common misconceptions require qualification. The first is that the decoupled equation for a gauge-invariant quantity is itself unique in non-vacuum backgrounds. In fact, for gravitational perturbations on general spherically symmetric non-vacuum spacetimes, different gauge choices produce different decoupled equations for the same gauge-invariant Weyl scalar $\pm1$16; the gauge dependence resides in the source sector, and the ambiguity disappears in vacuum (Guo et al., 2023). The second is that the radial Teukolsky–Starobinsky constants are always nonnegative. That statement is correct for spins up to $\pm1$17 but false for spin larger than $\pm1$18 (Costa et al., 2021). The third is that all Teukolsky generalizations preserve the simple Kerr picture of a single scalar equation. Kerr–Newman and non-vacuum type D settings show that this is generally too optimistic: one often obtains a structured coupled hierarchy rather than a fully decoupled scalar equation (Giorgi, 2020, Guo et al., 2023).

Finally, the significance of the Teukolsky equations lies in their persistence across diverse research programs. They remain central to black-hole ringdown and quasinormal-mode calculations on Kerr and modified Kerr backgrounds (Cano et al., 2024, Yu et al., 4 Nov 2025). They provide the organizing scalar PDE in nonlinear reconstruction schemes (Green et al., 2019). They admit a twistor-geometric reinterpretation that explains the appearance of weighted derivatives, curvature potentials, and hidden symmetries (Araneda, 2019). And in physical-space stability theory, they serve as the curvature-level starting point from which Fackerell–Ipser or Regge–Wheeler equations, decay estimates, and scattering theories are constructed (Pasqualotto, 2016, Masaood, 2020).

Taken together, these developments suggest that “Teukolsky equations” names not only the original Kerr master equation but an entire framework for treating extreme-spin fields on algebraically special black-hole spacetimes: separated or transformed, vacuum or sourced, classical or quantum, exact Kerr or perturbatively deformed.