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.

Post-Newtonian templates for phase evolution of spherical extreme mass ratio inspirals

A complete waveform comparison of post-Newtonian and numerical relativity in eccentric orbits

High-order post-Newtonian expansion of the generalized redshift invariant for eccentric-orbit, equatorial extreme-mass-ratio inspirals with a spinning primary

Quasicircular Orbital Parameters for Numerical Relativity Revisited

Asymptotically matched quasi-circular inspiral and transition-to-plunge in the small mass ratio expansion

Post-Newtonian expansion of energy and angular momentum fluxes: inclined spherical orbits about a Kerr black hole

Quasi-Keplerian parametrization for eccentric compact binaries in scalar-tensor theories at second post-Newtonian order and applications

Post-Newtonian expansions of extreme mass ratio inspirals of spinning bodies into Schwarzschild black holes

Tidal heating and torquing of the primary black hole in eccentric-orbit, non-spinning extreme-mass-ratio inspirals to 22PN order

Analytical model of precessing binaries using post-Newtonian theory in the extreme mass-ratio limit I: General Formalism

Post-Newtonian accurate parametric solution to the dynamics of spinning compact binaries in eccentric orbits: The leading order spin-orbit interaction

Optimizing post-Newtonian parameters and fixing the BMS frame for numerical-relativity waveform hybridizations

On the approximate relation between black-hole perturbation theory and numerical relativity

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

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

A Survey of Four Precessing Waveform Models for Binary Black Hole Systems

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

Towards efficient Effective One Body models for generic, non-planar orbits

A note on the conversion of orbital angles for extreme mass ratio inspirals

Post-Newtonian observables for aligned-spin binaries to sixth order in spin from gravitational self-force and Compton amplitudes