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

Manas Vishal,Scott E. Field,Katie Rink,Sigal Gottlieb,Gaurav Khanna

2024-08-31

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.

General Relativity and Quantum Cosmology,High Energy Astrophysical Phenomena,Computational Physics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to solve the numerical simulation problems of Extreme Mass Ratio Inspirals (EMRIs) systems, especially how to accurately handle the Teukolsky equation with Dirac delta function source terms. The EMRIs system consists of a low - mass celestial body (such as a star or a planet) orbiting a high - mass black hole (such as a supermassive black hole). This system plays an important role in low - frequency gravitational wave signals and is a key observational target for future space - based gravitational wave detectors (such as LISA).

### Specific problems

1. **Accurate handling of the Dirac delta function**:
- Existing time - domain solvers usually approximate the Dirac delta function through regularization methods, for example, using a narrow Gaussian distribution. This method introduces a small - scale length, resulting in systematic errors, especially near the low - mass celestial body.
- This paper proposes a multi - zone Discontinuous Galerkin (DG) method, which can accurately handle the Dirac delta function, thereby achieving global spectral accuracy at the source point location.

2. **Direct access to far - field waveforms**:
- In order to connect the near - field and the future null boundary, this paper adopts the Hyperboloidal Layer Method, which simplifies the outer boundary conditions and enables direct access to far - field waveforms.

3. **Verification by numerical experiments**:
- The author has carried out several numerical experiments, including convergence tests with exact solutions, energy radiation rates of circular orbits, super - convergence properties of the future null boundary, and late - time tail behaviors of scalar fields.
- These experiments verify the effectiveness and high accuracy of the method.

### Method overview

- **Multi - zone Discontinuous Galerkin method**:
- Expand the solution into spherical harmonic functions and convert the Teukolsky equation with a source into a one - dimensional symmetric hyperbolic system.
- Derive numerical fluxes to correctly handle the Dirac delta function and achieve global spectral accuracy.

- **Hyperboloidal Layer method**:
- Use hyperbolic coordinate transformation to simplify the outer boundary conditions and enable direct access to far - field waveforms.

- **Numerical experiments**:
- Include convergence tests with exact solutions, energy radiation rates of circular orbits, super - convergence properties of the future null boundary, and late - time tail behaviors of scalar fields.

### Conclusion

The methods proposed in this paper perform excellently when handling the Teukolsky equation with Dirac delta function source terms, especially achieving high accuracy at the source point location. These methods are of great significance for calculating gravitational self - force effects, because these effects require high accuracy at the location of the Dirac delta function.