Dudley–Finley Approximation in Kerr–Newman

• Dudley–Finley approximation is a separable method that freezes part of the coupled gravito-electromagnetic system to yield a Teukolsky-like master equation.
• It enables tractable quasinormal-mode spectroscopy using numerical methods like Leaver’s continued-fraction and asymptotic techniques such as matched expansions.
• The approach is exact for scalar fields and offers high accuracy for small charge or spin, though omitting coupling effects may lead to errors in regimes with strong charge-induced mixing.

Searching arXiv for Dudley-Finley approximation and related Kerr-Newman perturbation papers. The Dudley–Finley approximation is a separable perturbative framework for spin-weighted fields on nonvacuum or charged rotating black-hole backgrounds, most prominently Kerr–Newman, obtained by “freezing” part of the coupled perturbation sector and retaining a single Teukolsky-like master equation. In the Kerr–Newman setting, it replaces the full nonseparable gravito-electromagnetic system by a radial-angular system for a field of spin weight $s$, with exact status for scalar perturbations and approximate status for electromagnetic and gravitational perturbations (Zimmerman et al., 2015, Silva et al., 18 Feb 2025). The approximation has served as a tractable laboratory for quasinormal-mode spectroscopy, nearly extremal asymptotics, WKB analyses, continued-fraction computations, Seiberg–Witten correspondences, parametrized deviations from Kerr ringdown, and extensions to nonvacuum geometries such as hairy Kerr black holes and weakly charged Einstein–Maxwell–dilaton spacetimes (Zimmerman et al., 2015, Cano et al., 2024, Zi et al., 2023, Brito et al., 2018).

1. Definition and formal structure

In Boyer–Lindquist coordinates $\{t,r,\theta,\phi\}$, a Kerr–Newman black hole with mass $M$, specific angular momentum $a$, and charge $Q$ is described by

$$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$

The Dudley–Finley approximation seeks separated solutions of the form

$$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$

where $s\in\{-2,-1,0\}$ denotes the spin weight, $\omega\in\mathbb{C}$ the frequency, and $m\in\mathbb{Z}$ the azimuthal number (Silva et al., 18 Feb 2025).

The angular sector is governed by the spin-weighted spheroidal-harmonic equation. Writing $\{t,r,\theta,\phi\}$0, $\{t,r,\theta,\phi\}$1, and $\{t,r,\theta,\phi\}$2, one has

$\{t,r,\theta,\phi\}$3

with regularity at $\{t,r,\theta,\phi\}$4 selecting a discrete spectrum $\{t,r,\theta,\phi\}$5 (Silva et al., 18 Feb 2025).

The radial Dudley–Finley equation is a deformation of the Teukolsky equation. Defining

$\{t,r,\theta,\phi\}$6

the radial equation reads

$\{t,r,\theta,\phi\}$7

Equivalent forms appear in the literature with closely related sign conventions for $\{t,r,\theta,\phi\}$8 and the potential term (Silva et al., 18 Feb 2025, Zimmerman et al., 2015).

This separable master system is the core mathematical content of the approximation. In the scalar case $\{t,r,\theta,\phi\}$9, it is exact for Kerr–Newman scalar quasinormal modes; for electromagnetic and gravitational sectors it is an approximation because the full perturbation problem is coupled (Zimmerman et al., 2015, Saha et al., 6 Oct 2025).

2. Origin of the approximation and its assumptions

The approximation originates in the observation that the full linearized Kerr–Newman perturbation problem in Newman–Penrose variables is not separable in Boyer–Lindquist coordinates because gravitational and electromagnetic perturbations mix. Schematically, the full system contains coupled equations for the gravitational scalar $M$0 and the electromagnetic scalar $M$1 (Saha et al., 6 Oct 2025).

The Dudley–Finley prescription consists of “freezing” one perturbation sector to its background value. In practice one sets, for example, $M$2 in the gravitational equation or $M$3 in the electromagnetic equation and discards the other field, leaving a single separable equation for a spin-weighted field (Saha et al., 6 Oct 2025). In the language of later applications, the approximation applies when one ignores the Maxwell perturbations or, more generally, matter-fluid perturbations and treats them as frozen (Zi et al., 2023).

Several special limits delimit its status. For $M$4, the equation reduces exactly to the Kerr Teukolsky equation. For $M$5, the scalar test-field problem on Kerr–Newman is exact. No explicit small-$M$6 or small-$M$7 expansion is performed in deriving the approximation; rather, it is a field-freeze ansatz, although one expects better accuracy when $M$8 or $M$9 (Saha et al., 6 Oct 2025). A related characterization appears in parametrized-Teukolsky treatments, where the Dudley–Finley model is viewed as the ordinary Kerr–Teukolsky equation plus a small, charge-induced deformation of the effective potential, exact at $a$0 and complete through $a$1 (Cano et al., 2024).

The principal limitation is structural: the approximation does not capture the true coupling of gravitational and electromagnetic perturbations in Kerr–Newman (Zimmerman et al., 2015). This is the central reason it is exact for scalar perturbations, informative but approximate for gravito-electromagnetic spectra, and potentially unreliable in regimes where charge-induced coupling is strong (Cano et al., 2024, Saha et al., 6 Oct 2025).

3. Quasinormal modes, boundary conditions, and computational methods

Quasinormal modes are defined by imposing physically outgoing behavior at infinity and ingoing behavior at the event horizon. In the Dudley–Finley radial problem, near $a$2 one imposes

$a$3

while as $a$4 one imposes purely outgoing waves,

$a$5

These conditions define a discrete set of complex frequencies $a$6 (Silva et al., 18 Feb 2025).

A standard numerical approach is Leaver’s continued-fraction method. One expands a solution satisfying both asymptotic conditions as

$a$7

Substitution into the radial equation yields a three-term recursion

$a$8

and convergence imposes the continued-fraction condition

$a$9

One solves this simultaneously with the corresponding angular continued-fraction equation for $Q$0, using standard root-finding such as Müller’s method (Silva et al., 18 Feb 2025).

Analytic approximations complement continued fractions. In the nearly extremal regime, matched asymptotic expansions introduce the small parameter

$Q$1

and derive zero-damped mode formulas by matching inner and outer hypergeometric solutions (Zimmerman et al., 2015). In the eikonal regime $Q$2, WKB analysis treats the radial potential peak through

$Q$3

with the decay rate determined by the curvature of the peak (Zimmerman et al., 2015). These methods establish the Dudley–Finley approximation as both a numerical and asymptotic framework for Kerr–Newman spectroscopy.

4. Nearly extremal structure: zero-damped and damped branches

One of the principal uses of the Dudley–Finley approximation has been the study of nearly extremal Kerr–Newman black holes. In this regime, the horizons coalesce and the surface gravity becomes small. Zero-damped modes are those whose frequency approaches the corotation frequency $Q$4 with imaginary part of order $Q$5, where

$Q$6

Matched asymptotic analysis yields the first-order formula

$Q$7

or equivalently

$Q$8

with overtone index $Q$9 (Zimmerman et al., 2015).

The WKB treatment reveals two branches in the extremal limit. If the inclination parameter $$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$0 exceeds a critical value $$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$1, the potential peak moves to the horizon and $$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$2, yielding zero-damped modes. If $$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$3, a finite potential peak remains outside the horizon and one obtains damped modes with nonzero decay even at extremality (Zimmerman et al., 2015).

A later reassessment sharpened the boundary between regions containing only zero-damped modes and regions containing both zero-damped and damped modes. In eikonal form,

$$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$4

with

$$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$5

An equivalent criterion defines

$$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$6

and gives

$$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$7

Numerically, the $$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$8-criterion and $$\Delta(r)\equiv r^2-2Mr+a^2+Q^2,\qquad r_{\pm}=M\pm\sqrt{M^2-a^2-Q^2}.$$9-criterion agree extremely well even for $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$0 (Saha et al., 6 Oct 2025).

This branch structure has broader significance. The Dudley–Finley approximation “cleanly exhibits the phenomenon of spectrum bifurcation into zero-damped and damped branches in the nearly extremal limit,” and this feature is shared by the full gravito-electromagnetic spectrum (Zimmerman et al., 2015). A plausible implication is that the approximation isolates the near-extremal kinematics of Kerr–Newman resonances more robustly than it captures all mode-coupling effects.

5. Accuracy, isospectral reformulations, and Seiberg–Witten correspondence

The approximation has been quantitatively tested against more complete Kerr–Newman calculations. A recent comparison across the subextremal spin-charge parameter space found that for the $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$1, $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$2, and $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$3 modes, the agreement is typically within $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$4 and $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$5 for the real and imaginary parts of the frequencies, respectively (Saha et al., 6 Oct 2025). The same study reports that errors grow in the high-charge near-extremal corner, while for small charge or along the diagonal $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$6 they are $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$7 in both real and imaginary parts until near extremality (Saha et al., 6 Oct 2025). In parametrized-Teukolsky language, the quasi-normal frequencies agree with full Kerr–Newman numerics to better than a few percent in the small-charge regime $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$8, whereas for $$\Psi_s(t,r,\theta,\phi)=e^{-i\omega t+im\phi}S_s(\theta)R_s(r),$$9 or when spin and charge strongly couple the model ceases to track mode-mixing effects and higher-order charge corrections (Cano et al., 2024).

A distinct development is the relation between Dudley–Finley quasinormal modes and Seiberg–Witten theory. The radial equation can be mapped to the canonical form

$s\in\{-2,-1,0\}$0

by

$s\in\{-2,-1,0\}$1

Matching to the $s\in\{-2,-1,0\}$2, $s\in\{-2,-1,0\}$3 quantum Seiberg–Witten curve yields a dictionary including

$s\in\{-2,-1,0\}$4

$s\in\{-2,-1,0\}$5

$s\in\{-2,-1,0\}$6

A discrete exchange symmetry $s\in\{-2,-1,0\}$7 leaves the spectrum invariant and produces an isospectral reformulation (Silva et al., 18 Feb 2025).

After implementing the symmetry and the gauge transformation

$s\in\{-2,-1,0\}$8

the radial problem becomes

$s\in\{-2,-1,0\}$9

with

$\omega\in\mathbb{C}$0

where $\omega\in\mathbb{C}$1 and $\omega\in\mathbb{C}$2 (Silva et al., 18 Feb 2025). Since this real-potential Regge–Wheeler–like equation has the same quasinormal boundary conditions at $\omega\in\mathbb{C}$3 and $\omega\in\mathbb{C}$4, its complex spectrum is identical to that of the original Dudley–Finley equation (Silva et al., 18 Feb 2025).

The same Seiberg–Witten framework also yields numerical predictions for frequencies through Bohr–Sommerfeld quantization on the $\omega\in\mathbb{C}$5-period of the quantum curve. For the fundamental slowest-damped gravitational mode with $\omega\in\mathbb{C}$6, four representative points in units $\omega\in\mathbb{C}$7 were studied: $\omega\in\mathbb{C}$8, $\omega\in\mathbb{C}$9, $m\in\mathbb{Z}$0, and $m\in\mathbb{Z}$1. At $m\in\mathbb{Z}$2, Leaver’s continued fraction gives

$m\in\mathbb{Z}$3

With the Seiberg–Witten approach, the approximation converges oscillatory toward the continued-fraction value as the instanton order $m\in\mathbb{Z}$4 increases, and for $m\in\mathbb{Z}$5 instantons the result is within $m\in\mathbb{Z}$6 of the continued-fraction answer for all four parameter sets, although convergence slows closer to extremality (Silva et al., 18 Feb 2025).

6. Generalizations and broader applications

The Dudley–Finley idea has been adapted beyond the original Kerr–Newman gravito-electromagnetic setting. In parametrized ringdown studies, it appears as a benchmark example of a small deformation of the Kerr Teukolsky potential. Expanding the charge-induced deformation $m\in\mathbb{Z}$7 in a finite set of $m\in\mathbb{Z}$8 moments, one writes

$m\in\mathbb{Z}$9

and in the Dudley–Finley case only $\{t,r,\theta,\phi\}$00 are nonzero (Cano et al., 2024). This formulation treats the approximation as a prototype for linear-response calculations of quasinormal-mode shifts under small modifications of the Teukolsky potential.

In nonvacuum black-hole spacetimes, the same logic has been used as a practical separability ansatz. For the hairy Kerr black hole, Teukolsky’s derivation fails because the background is not vacuum, but by analogy with Kerr–Newman one obtains a gravitational perturbation equation by replacing $\{t,r,\theta,\phi\}$01 with

$\{t,r,\theta,\phi\}$02

The approximation is justified when matter-fluid perturbations are ignored, and it underlies the computation of energy fluxes, adiabatic inspiral evolution, dephasing, and mismatch for extreme mass ratio inspirals (Zi et al., 2023). The reported results demonstrate that the EMRI waveforms from the HKBH with deviation parameter larger than $\{t,r,\theta,\phi\}$03 and hair charge smaller than $\{t,r,\theta,\phi\}$04 can be discerned by LISA (Zi et al., 2023).

A modified Dudley–Finley scheme also appears in weakly charged Einstein–Maxwell–dilaton black holes. There, the exact polar perturbation problem is a coupled three-field system for electromagnetic, dilaton, and gravitational variables. Freezing the gravitational perturbation reduces the system to a coupled $\{t,r,\theta,\phi\}$05 electromagnetic-dilaton subsystem rather than two fully decoupled matter equations. This preserves the $\{t,r,\theta,\phi\}$06-dependent electromagnetic-scalar mixing responsible for polar-axial isospectrality breaking (Brito et al., 2018). At small charge $\{t,r,\theta,\phi\}$07, the resulting approximate electromagnetic-plus-scalar system reproduces the exact $\{t,r,\theta,\phi\}$08 quasinormal modes to better than $\{t,r,\theta,\phi\}$09, and even up to $\{t,r,\theta,\phi\}$10 the real part of the fundamental electromagnetic modes is accurate to a few percent (Brito et al., 2018).

These generalizations clarify that “Dudley–Finley approximation” denotes less a single fixed equation than a separability strategy: one freezes selected backreacting sectors, keeps a Teukolsky-type master structure where possible, and exploits the resulting tractable radial problem for spectroscopy or waveform generation. This suggests a unifying role for the approximation across charged, weakly nonvacuum, and parametrically deformed black-hole perturbation theory, while preserving its central caveat: omitted couplings remain the dominant source of error whenever they are not parametrically small.