Superconvergent discontinuous Galerkin method for the scalar Teukolsky equation on hyperboloidal domains: Efficient waveform and self-force computation

IF 2.8 4区 物理与天体物理 Q2 ASTRONOMY & ASTROPHYSICS
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,&nbsp;Scott E. Field,&nbsp;Sigal Gottlieb,&nbsp;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.

双曲域上标量Teukolsky方程的超收敛不连续Galerkin方法:有效波形和自力计算
极端质量比吸气系统的长期演化需要最小的相位和色散误差来精确计算远场波形,而在较小的黑洞附近(以狄拉克三角洲分布建模)进行自力计算则需要高精度。光谱精确的方法,如节点不连续伽辽金(DG)方法,非常适合这些任务。它们的数值误差通常减小为∝(Δx)N+1 \propto (\Delta x)^{N+1},其中Δx \Delta x是子域大小,NN是近似的多项式度。然而,某些DG方案表现出超收敛性,其中截断、相位和色散误差可以以∝(Δx)2N+1 \propto (\Delta x)^{2N+1}的速度减小。超收敛数值解算器在构造上是非常高效和精确的。我们从理论上证明了具有分布源的标量Teukolsky方程的DG格式是超收敛的,并且当与双曲层紧化技术结合使用时,这一性质仍然保持不变。这确保了波形、总能量和角动量通量以及自力计算受益于超收敛。我们通过经验验证了具有不同平滑度的双曲面层紧化族的这种行为。此外,我们还表明,在点粒子位置计算的圆形轨道的耗散自力量也表现出一定程度的超收敛。我们的研究结果强调了数值超收敛对基于DG方法的有效和准确的重力波形模拟的潜在好处。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
求助全文
约1分钟内获得全文 求助全文
来源期刊
General Relativity and Gravitation
General Relativity and Gravitation 物理-天文与天体物理
CiteScore
4.60
自引率
3.60%
发文量
136
审稿时长
3 months
期刊介绍: 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.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术官方微信