Kerr metric perturbations in Lorenz gauge via numerical relativity methods
Abstract
Gravitational self-force calculations are needed to accurately describe the motion of a small body in a spinning (Kerr) black hole spacetime, which is vital to predict gravitational wave signals from asymmetric-mass compact binaries. Accurate calculations for the upcoming LISA mission require a minimum of second perturbative order in the mass ratio, which has thus far only been achieved for nonspinning (Schwarzschild) black holes due to nonlinear complexities of the Einstein equations. In this work, we calculate Kerr metric perturbations with a new computational approach based on numerical relativity methods, extending recent scalar self-force results to full general relativity. We formulate the perturbative Lorenz-gauge Einstein equations in terms of azimuthal m-modes in the frequency domain and solve the resulting elliptic partial differential equations using a spectral discontinuous Galerkin discretization, adaptive mesh refinement, and preconditioned iterative solvers. With our new method, we achieve exponentially convergent first-order gravitational self-force calculations for circular equatorial orbits in Kerr spacetime up to spin |a|/M=0.998 (the Thorne limit). Our framework is designed to carry over to second order by delivering accurate and well-behaved first-order data for nonlinear source construction, and as a solver for the second-order metric perturbations themselves.
Related
- Perturbations and quasinormal modes of black holes in general relativity coupled to nonlinear electrodynamics
- A discontinuous Galerkin method with artificial viscosity geared toward numerical relativity applications
- A Transformer model for precessing binary black-hole merger waveforms calibrated to numerical relativity
- Gravitational waves from core collapse of rotating very-massive stars: 3D numerical relativity computation
- Locating critical solutions in numerical relativity using automatic differentiation
BibTeX
@misc{nishimura2026kerr,
title = {{Kerr metric perturbations in Lorenz gauge via numerical relativity methods}},
author = {Nami Nishimura and Nils L. Vu and Adrià Gispert-Fernández and Cillian Kelly and Lawrence E. Kidder and Benjamin Leather and Thomas Osburn and Jonathan E. Thompson and Samuel D. Upton and Barry Wardell},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.05259},
eprint = {2611.05259},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.05259},
}
RIS
TY - UNPB
TI - Kerr metric perturbations in Lorenz gauge via numerical relativity methods
AU - Nami Nishimura
AU - Nils L. Vu
AU - Adrià Gispert-Fernández
AU - Cillian Kelly
AU - Lawrence E. Kidder
AU - Benjamin Leather
AU - Thomas Osburn
AU - Jonathan E. Thompson
AU - Samuel D. Upton
AU - Barry Wardell
PY - 2026
DA - 2026/11/03
PB - arXiv
DO - 10.48550/arxiv.2611.05259
UR - https://arxiv.org/abs/2611.05259
AB - Gravitational self-force calculations are needed to accurately describe the motion of a small body in a spinning (Kerr) black hole spacetime, which is vital to predict gravitational wave signals from asymmetric-mass compact binaries. Accurate calculations for the upcoming LISA mission require a minimum of second perturbative order in the mass ratio, which has thus far only been achieved for nonspinning (Schwarzschild) black holes due to nonlinear complexities of the Einstein equations. In this work, we calculate Kerr metric perturbations with a new computational approach based on numerical relativity methods, extending recent scalar self-force results to full general relativity. We formulate the perturbative Lorenz-gauge Einstein equations in terms of azimuthal m-modes in the frequency domain and solve the resulting elliptic partial differential equations using a spectral discontinuous Galerkin discretization, adaptive mesh refinement, and preconditioned iterative solvers. With our new method, we achieve exponentially convergent first-order gravitational self-force calculations for circular equatorial orbits in Kerr spacetime up to spin |a|/M=0.998 (the Thorne limit). Our framework is designed to carry over to second order by delivering accurate and well-behaved first-order data for nonlinear source construction, and as a solver for the second-order metric perturbations themselves.
ER -