Spiral Density Waves and Torque Balance in the Kerr Geometry

Extreme mass-ratio inspirals (EMRIs) in relativistic accretion discs are a key science target for the upcoming LISA mission. Existing models of disc-EMRI interactions typically rely on crude dynamical friction or Newtonian planetary migration prescriptions, which fail to capture the relativistic fluid response induced by the binary potential. In this work we address this gap by providing the relativistic calculation. We apply standard methods from self-force theory, black hole perturbation theory, and relativistic stellar perturbation theory to perform the full fluid calculation of the relativistic analogue of planetary migration for the first time. We calculate the response of a fluid in the perturbing potential of an EMRI consistently incorporating pressure effects. Using a master enthalpy-like variable and linearised fluid theory, we reconstruct the fluid perturbations and relativistic spiral arm structure for a range of spin values in the Kerr geometry. We conclude by deriving a relativistic torque-balance equation that enables computation and comparison of local torques with advected angular momentum through the disc. This opens a promising route towards establishing torque-balance relations between integrated disc torques arising from fluid perturbations and the forces acting on EMRIs embedded in matter.

💡 Research Summary

This paper presents the first fully relativistic treatment of fluid‑disk perturbations induced by an extreme‑mass‑ratio inspiral (EMRI) orbiting a spinning (Kerr) supermassive black hole, with the aim of providing a rigorous foundation for torque‑balance calculations relevant to future LISA observations. The authors begin by formulating the Einstein equations coupled to a perfect‑fluid stress‑energy tensor and introduce two small expansion parameters: the mass‑ratio q = m_p/M and a dimensionless measure λ of the background disk’s self‑gravity and pressure. The metric is written as a Kerr background plus a perturbation h_{μν}, while the secondary is modeled as a point‑particle source. By expanding to order O(λq) they obtain linearised fluid equations that include both the metric perturbation and the direct force from the particle.

A key innovation is the reduction of the coupled linearised continuity and Euler equations to a single scalar “master enthalpy” equation. This is achieved by assuming the background disk flow respects the Kerr Killing symmetries (stationary, axis‑symmetric) and by imposing an adiabatic relation Δp = c_s² Δe in the fluid’s comoving frame. The Lagrangian displacement ξ^μ is expressed in terms of the velocity perturbation, and the enthalpy‑like variable h^{(1,0)} is defined via p^{(1,1)} = (e + p) h^{(1,0)}. After eliminating the velocity perturbations using an operator Q^{-1}, the authors arrive at a second‑order partial differential equation for h^{(1,0)} that incorporates Kerr frame‑dragging, the radial and polar dependence of the sound speed, and a thermal‑structure term B_μ = −∂_μ log c_s²/(1 + c_s²). This master equation generalises the classic Goldreich‑Treiman‑Papaloizou‑Lin Newtonian density‑wave formalism to full general relativity.

The paper then constructs a stationary background disk solution in Boyer‑Lindquist coordinates, specifying the pressure, density, and sound‑speed profiles consistent with a thin, adiabatic accretion flow around a Kerr black hole. Numerical methods are described in detail: the master equation is cast as an eigenvalue problem for complex frequency ω and azimuthal mode number m, solved using a hybrid spectral‑finite‑difference scheme with appropriate inner‑resonance and outer‑radiation boundary conditions. The operator Q^{-1} is inverted efficiently via matrix decomposition and preconditioning.

Results are presented for a range of black‑hole spins (a = 0–0.99 M). The authors find that higher spin enhances frame‑dragging, leading to tighter, more asymmetric spiral density waves. The morphology of the spirals (pitch angle, wavelength) varies systematically with spin and with the radial temperature gradient encoded in c_s(r). Importantly, the authors derive a relativistic torque‑balance equation that equates the local torque density exerted by the spiral wave on the secondary to the advected angular‑momentum flux through the disk. This balance includes novel terms proportional to the Lagrangian displacement and the enthalpy perturbation, absent in Newtonian treatments. The new formula enables a direct comparison between the torque computed from the fluid response and the angular‑momentum loss inferred from the particle’s inspiral, providing a self‑consistent check on EMRI waveform models that include environmental effects.

In the discussion, the authors emphasise that their framework bridges self‑force theory, black‑hole perturbation theory, and relativistic stellar‑perturbation techniques, offering a unified approach to EMRI‑disk interactions. They argue that incorporating these relativistic torque‑balance corrections will be essential for accurate parameter estimation and tests of general relativity with LISA, especially for sources embedded in dense AGN disks. Limitations are acknowledged: the analysis is linear, assumes a barotropic equation of state, neglects magnetic fields, and treats the disk as geometrically thin. Future work is outlined, including extensions to non‑linear regimes, magnetohydrodynamic disks, and fully three‑dimensional simulations.

Overall, the paper delivers a comprehensive theoretical and numerical toolkit for modelling relativistic spiral density waves and torque exchange in Kerr spacetimes, representing a significant step toward realistic EMRI waveform modelling in realistic astrophysical environments.