Superconvergent Discontinuous Galerkin Method for the Scalar Teukolsky Equation on Hyperboloidal Domains: Efficient Waveform and Self-Force Computation

Abstract

The long-time evolution of extreme mass-ratio inspiral systems requires minimal phase and dispersion errors to accurately compute far-field waveforms, while high accuracy is essential near the smaller black hole (modeled as a Dirac delta distribution) for self-force computations. Spectrally accurate methods, such as nodal discontinuous Galerkin (DG) methods, are well suited for these tasks. Their numerical errors typically decrease as $\propto (\Delta x{)}^{N+1}$, where $\Delta x$ is the subdomain size and $N$ is the polynomial degree of the approximation. However, certain DG schemes exhibit superconvergence, where truncation, phase, and dispersion errors can decrease as fast as $\propto (\Delta x{)}^{2N+1}$. Superconvergent numerical solvers are, by construction, extremely efficient and accurate. We theoretically demonstrate that our DG scheme for the scalar Teukolsky equation with a distributional source is superconvergent, and this property is retained when combined with the hyperboloidal layer compactification technique. This ensures that waveforms, total energy and angular-momentum fluxes, and self-force computations benefit from superconvergence. We empirically verify this behavior across a family of hyperboloidal layer compactifications with varying degrees of smoothness. Additionally, we show that self-force quantities for circular orbits, computed at the point particle's location, also exhibit a certain degree of superconvergence. Our results underscore the potential benefits of numerical superconvergence for efficient and accurate gravitational waveform simulations based on DG methods.

arXiv comment14 pages, 7 figures, and 2 tables. Accepted to the General Relativity and Gravitation journal's topical collection on hyperboloidal foliations. Version 2 includes an extended discussion of long-lived "junk tails" at future null infinity and additional dissipative self-force computations

References in corpus (28)

• The motion of point particles in curved spacetime
• Intermediate and Extreme Mass-Ratio Inspirals -- Astrophysics, Science Applications and Detection using LISA
• Self-force and radiation reaction in general relativity
• Gravitational self force in extreme mass-ratio inspirals
• The evolution of circular, non-equatorial orbits of Kerr black holes due to gravitational-wave emission
• A consequence of the gravitational self-force for circular orbits of the Schwarzschild geometry
• Pseudospectrum and black hole quasi-normal mode (in)stability
• Hyperboloidal foliations and scri-fixing
• Black hole perturbation theory and gravitational self-force
• Perspective on gravitational self-force analyses
• Binary black hole coalescence in the large-mass-ratio limit: the hyperboloidal layer method and waveforms at null infinity
• Spectral decomposition of black-hole perturbations on hyperboloidal slices
• Hyperboloidal layers for hyperbolic equations on unbounded domains
• Self force on a scalar charge in Kerr spacetime: circular equatorial orbits
• Numerical solution of the 2+1 Teukolsky equation on a hyperboloidal and horizon penetrating foliation of Kerr and application to late-time decays
• Self-consistent orbital evolution of a particle around a Schwarzschild black hole
• Discontinuous Galerkin method for computing gravitational waveforms from extreme mass ratio binaries
• Asymptotics of Schwarzschild black hole perturbations
• An Efficient Pseudospectral Method for the Computation of the Self-force on a Charged Particle: Circular Geodesics around a Schwarzschild Black Hole
• Rapid Evaluation of Radiation Boundary Kernels for Time-domain Wave Propagation on Blackholes
• Analytic structure of radiation boundary kernels for blackhole perturbations
• Persistent junk solutions in time-domain modeling of extreme mass ratio binaries
• Gravitational radiation reaction
• Fast evaluation of far-field signals for time-domain wave propagation
• Fast evaluation of asymptotic waveforms from gravitational perturbations
• Self-force on a scalar charge in a circular orbit about a Reissner-Nordström black hole
• Toward exponentially-convergent simulations of extreme-mass-ratio inspirals: A time-domain solver for the scalar Teukolsky equation with singular source terms
• Discontinuous Galerkin method for linear wave equations involving derivatives of the Dirac delta distribution

Cited by in corpus (1)

• Topical Collection-Hyperboloidal Foliations in the Era of Gravitational-Wave Astronomy: From Mathematical Relativity to Astrophysics