Superconvergent discontinuous Galerkin method for the scalar Teukolsky equation on hyperboloidal domains: Efficient waveform and self-force computation
Manas Vishal, Scott E. Field, Sigal Gottlieb, Jennifer Ryan
{"title":"Superconvergent discontinuous Galerkin method for the scalar Teukolsky equation on hyperboloidal domains: Efficient waveform and self-force computation","authors":"Manas Vishal, Scott E. Field, Sigal Gottlieb, Jennifer Ryan","doi":"10.1007/s10714-025-03435-9","DOIUrl":null,"url":null,"abstract":"<div><p>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 <span>\\(\\propto (\\Delta x)^{N+1}\\)</span>, where <span>\\(\\Delta x\\)</span> is the subdomain size and <span>\\(N\\)</span> is the polynomial degree of the approximation. However, certain DG schemes exhibit superconvergence, where truncation, phase, and dispersion errors can decrease as fast as <span>\\(\\propto (\\Delta x)^{2N+1}\\)</span>. 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 dissipative 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.</p></div>","PeriodicalId":578,"journal":{"name":"General Relativity and Gravitation","volume":"57 7","pages":""},"PeriodicalIF":2.8000,"publicationDate":"2025-06-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"General Relativity and Gravitation","FirstCategoryId":"101","ListUrlMain":"https://link.springer.com/article/10.1007/s10714-025-03435-9","RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"ASTRONOMY & ASTROPHYSICS","Score":null,"Total":0}
引用次数: 0
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 dissipative 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.
期刊介绍:
General Relativity and Gravitation is a journal devoted to all aspects of modern gravitational science, and published under the auspices of the International Society on General Relativity and Gravitation.
It welcomes in particular original articles on the following topics of current research:
Analytical general relativity, including its interface with geometrical analysis
Numerical relativity
Theoretical and observational cosmology
Relativistic astrophysics
Gravitational waves: data analysis, astrophysical sources and detector science
Extensions of general relativity
Supergravity
Gravitational aspects of string theory and its extensions
Quantum gravity: canonical approaches, in particular loop quantum gravity, and path integral approaches, in particular spin foams, Regge calculus and dynamical triangulations
Quantum field theory in curved spacetime
Non-commutative geometry and gravitation
Experimental gravity, in particular tests of general relativity
The journal publishes articles on all theoretical and experimental aspects of modern general relativity and gravitation, as well as book reviews and historical articles of special interest.