A discontinuous Galerkin method for the distributionally-sourced s=0 Teukolsky equation

Extreme Mass Ratio Inspirals (EMRI) will be an important astrophysical source for the upcoming gravitational wave detector Laser Interferometer Space Antenna (LISA). EMRI systems are characterized by mass ratios greater than 10⁴ and will be in the detector's sensitive band for hundreds of thousands of orbital cycles. Numerical relativity codes cannot be used to simulate such systems, so point-particle black hole perturbation theory (ppBHPT) is typically used instead. For perturbations of Kerr, time-domain numerical solvers for the Teukolsky equation have traditionally been based on low-order finite-difference methods. In this talk I will discuss a spectrally-accurate Discontinuous Galerkin method for the distributionally-forced Teukolsky equation. Importantly, the method maintains spectral accuracy even at the location of the delta function. I will present numerical results from the code, including the computation of Price tails and scalar field evolutions from a perturbing black hole in a circular orbit.