Abstract:Gravitational wave signals from extreme mass ratio inspirals are a key target for space-based gravitational wave detectors. These systems are typically modeled as a distributionally-forced Teukolsky equation, where the smaller black hole is treated as a Dirac delta distribution. Time-domain solvers often use regularization approaches that approximate the Dirac distribution that often introduce small length scales and are a source of systematic error, especially near the smaller black hole. We describe a multi-domain discontinuous Galerkin method for solving the distributionally-forced Teukolsky equation that describes scalar fields evolving on a Kerr spacetime. To handle the Dirac delta, we expand the solution in spherical harmonics and recast the sourced Teukolsky equation as a first-order, one-dimensional symmetric hyperbolic system. This allows us to derive the method's numerical flux to correctly account for the Dirac delta. As a result, our method achieves global spectral accuracy even at the source's location. To connect the near field to future null infinity, we use the hyperboloidal layer method, allowing us to supply outer boundary conditions and providing direct access to the far-field waveform. We document several numerical experiments where we test our method, including convergence tests against exact solutions, energy luminosities for circular orbits, the scheme's superconvergence properties at future null infinity, and the late-time tail behavior of the scalar field. We also compare two systems that arise from different choices of the first-order reduction variables, finding that certain choices are numerically problematic in practice. The methods developed here may be beneficial when computing gravitational self-force effects, where the regularization procedure has been developed for the spherical harmonic modes and high accuracy is needed at the Dirac delta's location.

Gravitational wave fluxes on generic orbits in near-extreme Kerr spacetime: higher spin and large eccentricity

Gravitational wave fluxes on generic orbits in near-extreme Kerr spacetime: Higher spin and large eccentricity

Asymptotic gravitational-wave fluxes from a spinning test body on generic orbits around a Kerr black hole

Spinning particles near Kerr black holes: Orbits and gravitational-wave fluxes through the Hamilton-Jacobi formalism

Gravitational waves from extreme-mass-ratio inspirals in the semiclassical gravity spacetime

Gravitational waves from extreme mass ratio inspirals in Kerr-MOG spacetimes

Geometrized effective-one-body formalism for extreme-mass-ratio limits: Generic orbits

Gravitational waveforms from periodic orbits around a quantum-corrected black hole

Spin effects in gravitational waveforms and fluxes for binaries on eccentric orbits to the third post-Newtonian order

Calculating the Gravitational Waves Emitted from High-speed Sources

Gravitational waves for eccentric extreme mass ratio inspirals of self-dual spacetime

Extreme mass-ratio inspiral and waveforms for a spinning body into a Kerr black hole via osculating geodesics and near-identity transformations

Generalized Gibbons-Werner method for stationary spacetimes

Spin wave optics for gravitational waves lensed by a Kerr black hole

Toward exponentially-convergent simulations of extreme-mass-ratio inspirals: A time-domain solver for the scalar Teukolsky equation with singular source terms

Gravitational waves from extreme mass ratio inspirals around a hairy Kerr black hole

Gravitational waves from eccentric extreme mass-ratio inspirals as probes of scalar fields