Spin-k Teukolsky Operator in Kerr Perturbations

• Spin-k Teukolsky Operator is a separable second-order operator governing extreme spin-weighted perturbations on Kerr black holes.
• It decomposes into radial and angular equations, enabling precise calculations of quasi-normal modes and gravitational self-force effects.
• Generalizations extend its framework to charged, conformal, and higher-derivative settings, offering broader insights into perturbative dynamics.

The spin-$k$ Teukolsky operator, almost uniformly denoted in the literature by a spin parameter $s$, is the separable second-order operator governing extreme spin-weighted perturbation variables on algebraically special black-hole backgrounds, most prominently Kerr. In its classical form it acts on Newman–Penrose curvature or Maxwell scalars; in Boyer–Lindquist coordinates it separates into radial and angular pieces, converting the perturbation problem into an angular eigenvalue problem and a radial ordinary differential equation. This separability underlies quasi-normal-mode calculations, radiation extraction, self-force formulations, scattering theory, and a wide range of extensions to charged, non-vacuum, conformal, and higher-derivative settings (Spiers, 2024).

1. Definition and basic operator forms

In the modern Kerr literature, the generic operator is written for a spin weight $s$, with physical cases including $s=0,\pm1,\pm2$ in the classical Teukolsky framework, and broader treatments also covering half-integer and higher-spin sectors in related formulations. A compact covariant expression used in GHP language is

$${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$

with adjoint operator ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$. In the Kinnersley frame this admits a charged Klein–Gordon interpretation through $D_a=\nabla_a+s\Gamma_a$, so that the operator takes the form $g^{ab}D_aD_b-4s^2\Psi_2$ (Iuliano et al., 2023).

A coordinate form widely used in Kerr is the Boyer–Lindquist representation of the modified Teukolsky equation,

$$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$

where setting $\eta=0$ recovers the standard Kerr operator. In Boyer–Lindquist coordinates $s$0, with $s$1, $s$2, and $s$3, the operator contains the characteristic terms

$s$4

together with the standard $s$5, $s$6, $s$7, and spin-coupling contributions. This is the form from which most separated radial and angular equations are derived (Yu et al., 4 Nov 2025).

The operator is tied to algebraic speciality. On Kerr, the Petrov type-D and vacuum structure permit decoupled equations for the extreme NP curvature scalars $s$8 and $s$9, or for the extreme Maxwell scalars in the spin-$s$0 sector. In the Schwarzschild limit, the same formalism reduces to the spin-$s$1 and spin-$s$2 equations studied through Fackerell–Ipser and Regge–Wheeler transformations rather than through rotation-induced spheroidal structure (Pasqualotto, 2016).

2. Separation on Kerr

The operational hallmark of the Teukolsky operator is its decomposition into radial and angular pieces. For the sourced master equation $s$3, one writes

$s$4

and in Boyer–Lindquist coordinates with $s$5 time dependence,

$s$6

$s$7

with $s$8, $s$9, and $s=0,\pm1,\pm2$0. The angular eigenfunctions of $s=0,\pm1,\pm2$1 are the spin-weighted spheroidal harmonics $s=0,\pm1,\pm2$2 (Spiers, 2024).

A separated field is expanded as

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

and, when a source is present,

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

Each mode then satisfies the radial ODE

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

For the homogeneous problem this reduces to the standard radial Teukolsky equation, while the angular equation yields the spheroidal separation constant $s=0,\pm1,\pm2$6 or equivalent shifted eigenvalues used in different conventions (Spiers, 2024).

This separation is not merely formal. It converts a $s=0,\pm1,\pm2$7 dimensional PDE problem into decoupled radial ODEs in the frequency domain, and in time-domain formulations it reduces the problem to coupled $s=0,\pm1,\pm2$8 systems with nearest-neighbor or next-to-nearest-neighbor $s=0,\pm1,\pm2$9-mode coupling. For quasinormal modes, the separated radial solutions are further constrained by the usual boundary conditions: purely ingoing behavior at the horizon and purely outgoing behavior at infinity (Spiers, 2024, Yu et al., 4 Nov 2025).

3. Spin sectors, spin reversal, and operator identities

The classical operator organizes perturbations by spin weight. In Kerr, $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$0 and $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$1 correspond to the gravitational $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$2 sectors, while $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$3 and $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$4 provide the electromagnetic $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$5 sectors. On Schwarzschild, the spin-$${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$6 equations for the extreme Maxwell components $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$7 and $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$8 can be transformed by a physical-space Chandrasekhar-type map into the tensorial Fackerell–Ipser equation, which admits direct Morawetz and $${}_s\mathcal O\,{}_s\psi = \left[g^{ab}(\Theta_a+2sB_a)(\Theta_b+2sB_b)-4s^2\Psi_2\right]{}_s\psi=0,$$9 estimates. That transformation is the spin-${}_s\mathcal O^\dagger={}_{-s}\mathcal O$0 analogue of the spin-${}_s\mathcal O^\dagger={}_{-s}\mathcal O$1 Teukolsky-to-Regge–Wheeler strategy (Pasqualotto, 2016).

For the gravitational ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$2 equations on Schwarzschild, a physical-space Chandrasekhar transformation produces Regge–Wheeler variables ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$3 and ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$4 satisfying a Lagrangian wave equation. This makes it possible to construct a full scattering theory, with unitary forward and backward maps between Cauchy data and scattering data on ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$5 and ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$6. The same analysis also yields a physical-space form of the Teukolsky–Starobinsky identities linking the ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$7 and ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$8 sectors (Masaood, 2020).

Spin reversal is encoded by the Teukolsky–Starobinsky identities and the associated Teukolsky–Starobinsky constants. Under separation of variables, differential chains of order ${}_s\mathcal O^\dagger={}_{-s}\mathcal O$9 map spin $D_a=\nabla_a+s\Gamma_a$0 solutions to spin $D_a=\nabla_a+s\Gamma_a$1 solutions and back. A notable correction to the standard lore is that these constants are not always positive: they can be negative for spin larger than $D_a=\nabla_a+s\Gamma_a$2. The same body of work shows that for $D_a=\nabla_a+s\Gamma_a$3 the positivity properties needed in the canonical quantization scheme do hold for Maxwell and gravitational perturbations, and that positivity is necessary and sufficient there for implementing canonical commutation relations (Costa et al., 2021, Iuliano et al., 2023).

4. Generalizations beyond vacuum Kerr

The most direct extension is the Kerr–Newman case, where gravity and electromagnetism couple already at the linearized level. In that setting the natural scalar operator becomes a generalized Teukolsky operator $D_a=\nabla_a+s\Gamma_a$4, depending not only on the spin $D_a=\nabla_a+s\Gamma_a$5 but also on a conformal type $D_a=\nabla_a+s\Gamma_a$6. In Boyer–Lindquist coordinates it has the schematic form

$D_a=\nabla_a+s\Gamma_a$7

When $D_a=\nabla_a+s\Gamma_a$8 and $D_a=\nabla_a+s\Gamma_a$9, this reduces to the standard Kerr Teukolsky operator. For $g^{ab}D_aD_b-4s^2\Psi_2$0, however, the spin-1 and spin-2 sectors remain coupled, and standard mode decomposition fails because the relevant spheroidal harmonics of different spins are not simply related by the angular operators appearing in the coupled system (Giorgi, 2020).

A distinct non-vacuum extension arises on axisymmetric type-D backgrounds that are conformal to Kerr with a radial deformation $g^{ab}D_aD_b-4s^2\Psi_2$1. In the case where the conformal factor is set to unity, the separated angular equations remain the standard spin-weighted spheroidal harmonic equations, but the radial equations acquire the replacement $g^{ab}D_aD_b-4s^2\Psi_2$2. For the $g^{ab}D_aD_b-4s^2\Psi_2$3 sector, the homogeneous radial equation becomes

$g^{ab}D_aD_b-4s^2\Psi_2$4

which reduces to the standard Teukolsky radial equation when $g^{ab}D_aD_b-4s^2\Psi_2$5 (Guo et al., 2023).

In spherically symmetric non-vacuum type-D spacetimes, a unified Teukolsky-like operator can be written for $g^{ab}D_aD_b-4s^2\Psi_2$6. After a rescaling by $g^{ab}D_aD_b-4s^2\Psi_2$7, the master equation becomes

$g^{ab}D_aD_b-4s^2\Psi_2$8

Here the non-vacuum background enters explicitly through $g^{ab}D_aD_b-4s^2\Psi_2$9, and the separated angular functions become spin-weighted spherical harmonics rather than spheroidal harmonics (Guo et al., 2023).

A conformal reformulation is also available. Linearizing Friedrich’s conformal Einstein field equations around a Petrov type-D background yields a conformal Teukolsky equation for the NP components $$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$0 and $$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$1 of the rescaled Weyl tensor. The resulting operator preserves the classical NP derivative structure but remains regular at the conformal boundary, with the background potential appearing as $$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$2. In the Kerr applications discussed there, this provides a geometric interpretation of hyperboloidal master variables as components of the rescaled Weyl tensor in an unphysical tetrad (Gasperin et al., 22 Feb 2026).

A further abstraction is the “universal Teukolsky equations” program, where Teukolsky-type operators are derived as specific combinations of the differential Bianchi identities and then linearized on higher-derivative or non-vacuum backgrounds. In this formulation, theory dependence enters through effective stress-energy terms or background corrections, while the underlying decoupling strategy remains tied to the type-D structure and to NP/GHP algebra (Cano et al., 2023).

5. Sourced operators and computational frameworks

The homogeneous operator is classical; the sourced operator is a current research frontier. On Kerr, the sourced master equation can be written as

$$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$3

but the source is generically not separable because Kerr background factors such as

$$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$4

mix radial and angular dependence. This becomes especially serious for extended sources in second-order perturbation theory, self-force calculations, and quadratic quasinormal-mode source terms. An analytic source-separation procedure has now been constructed by combining Newman–Penrose, GHP, and Held formalisms, re-expanding the source in spin-weighted spherical harmonics, isolating all residual $$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$5-$$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$6 mixing into a common factor $$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$7, and approximating that factor by a controlled Fourier series in $$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$8. The result is a separated source

$$[(\nabla^a+s\Gamma^a)(\nabla_a+s\Gamma_a)-4s^2\Psi_2]\Psi^{(s)} = (\eta/\lambda^2)\Psi^{(s)},$$9

without performing a sphere integral at every radial grid point (Spiers, 2024).

That same sourced perspective underlies nonlinear Kerr perturbation theory. A sourced linearized Einstein solution can be decomposed into a pure gauge piece, a zero mode, a metric reconstructed from a Debye–Hertz potential $\eta=0$0, and a corrector tensor $\eta=0$1. The potential satisfies a sourced adjoint spin-$\eta=0$2 Teukolsky equation

$\eta=0$3

while the corrector and the source $\eta=0$4 are obtained from decoupled radial ODEs in outgoing Kerr–Newman coordinates. This reduces the tensorial sourced problem to a scalar Teukolsky equation plus radial integrations, and it is designed to iterate order by order in nonlinear perturbation theory (Green et al., 2019).

In the time domain, the point-particle problem has recently been reformulated in comoving, spatially compactified hyperboloidal coordinates. For the scalar testbed and for the $\eta=0$5 Hertz potential on Schwarzschild, the field equations become $\eta=0$6 dimensional PDEs with jump conditions enforced at a fixed grid location. The central numerical observation is that hyperboloidal compactification suppresses the nonphysical growing modes that appear in some uncompactified Teukolsky evolutions, while the comoving coordinate choice makes the jump conditions local and straightforward to impose (Vaswani et al., 24 Jun 2026).

6. Spectral deformations, stability theory, and current directions

One major current direction treats the spin-$\eta=0$7 operator as a parameterized object. In one framework, both the radial and angular equations are deformed by two independent parameter families,

$\eta=0$8

so that

$\eta=0$9

$s$00

The resulting quasinormal-mode problem is solved by coupled determinant conditions for $s$01 and the separation constant, and cross-validated against a $s$02-dimensional hyperboloidal pseudo-spectral calculation (Yu et al., 4 Nov 2025).

A closely related earlier framework parameterized only the radial sector through a $s$03 expansion of the effective potential,

$s$04

and extracted universal linear response coefficients $s$05 and $s$06 for the QNM frequencies and angular separation constants. That program was tested against the massive scalar field, the Dudley–Finley equation, and higher-derivative gravity, and it provides a direct route from a model-dependent $s$07 to linearized QNM shifts (Cano et al., 2024).

Near extremality, the operator simplifies in a different way. In higher-derivative gravity, the near-horizon modified Teukolsky equation for NHEK and near-NHEK geometries still separates, but the radial equation keeps exactly the Kerr/NHEK structure while the angular separation constants are shifted by higher-derivative corrections: $s$08 This is enough to reproduce recent horizon-singularity results for extremal Kerr under certain deformations and to extend them to parity-breaking corrections, while simultaneously constraining the form of the full modified radial potential away from the strict near-horizon limit (Cano et al., 2024).

At the PDE level, the operator now supports full subextremal energy–Morawetz–flux theory even on perturbations of Kerr. In a tensorial non-integrable formalism, the spin-$s$09 Teukolsky system is recast as a wave/transport hierarchy, then regular scalarization via a globally regular triplet of horizontal one-forms reduces the tensorial system to scalarized wave equations with controlled lower-order couplings. This yields energy–Morawetz–flux estimates for perturbations of Kerr compatible with nonlinear stability applications (Ma et al., 24 Mar 2026).

A plausible implication is that the “spin-$s$10 Teukolsky operator” is no longer a single fixed differential operator attached only to vacuum Kerr. In current usage it denotes a family of closely related spin-dependent operators: classical Kerr operators; sourced and adjoint operators; generalized $s$11 operators on Kerr–Newman; non-vacuum and conformal type-D analogues; universal higher-derivative operators; and parameterized deformations for theory-agnostic spectroscopy. Across these settings, the common structural themes remain separability, spin reversal through Teukolsky–Starobinsky identities, and the reduction of perturbative dynamics to radial and angular operator theory (Giorgi, 2020, Cano et al., 2023).