[论文解读] 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+3 more|arXiv (Cornell University)|Jul 3, 2023

Pulsars and Gravitational Waves ResearchPhysics and Astronomy被引用 3

一句话总结

本文提出了一种用于克尔时空中标量泰乌克尔斯基方程的多域不连续伽辽金求解器，采用球谐分解与一阶对称双曲形式化，即使在存在奇异狄拉克δ源项的情况下也能实现全局谱精度。该方法实现了极端质量比旋进（EMRI）的指数收敛模拟，能够在未来 null 无穷远处准确获取波形，并正确处理点粒子奇点。

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.

研究动机与目标

• 开发一种用于标量泰乌克尔斯基方程的时域求解器，处理极端质量比旋进（EMRIs）中出现的奇异源项。
• 通过将狄拉克δ函数直接嵌入数值格式中，消除因近似狄拉克δ函数而引入的正则化误差。
• 利用具有适当通量处理的不连续伽辽金方法，在点粒子源位置也实现全局谱精度。
• 通过双曲层方法直接计算未来 null 无穷远处的波形，避免人工边界条件。

提出的方法

• 将标量场在球谐函数上展开，将3+1维泰乌克尔斯基方程简化为在 tortoise 坐标与时间上的1+1维系统。
• 通过引入变量 πₗₘ = -∂ψₗₘ/∂τ 和 φₗₘ = ∂ψₗₘ/∂ρ，将二阶波动方程重写为一阶对称双曲系统。
• 应用多域不连续伽辽金方法，其数值通量特别推导以正确处理狄拉克δ源项。
• 采用双曲层方法将未来 null 无穷远映射到有限边界，从而通过简单的外边界条件直接获取远场波形。
• 该方法确保强双曲性与稳定性，特别处理了右边界（ρ = s）处一个波速为零的情况。
• 比较了两种不同的首阶约化变量选择，发现尽管理论上等价，其中一种在数值上不稳定。

实验结果

研究问题

• RQ1 不连续伽辽金方法是否能在狄拉克δ源项存在的情况下，即使在奇点位置也实现全局谱收敛？
• RQ2 所提出的具有适当通量处理的对称双曲形式是否能消除因近似狄拉克函数而引入的正则化误差？
• RQ3 首阶约化变量的选择如何影响点粒子源存在下的数值稳定性与收敛性？
• RQ4 双曲层方法能否与DG方法有效结合，以在无人工边界条件的情况下直接提取未来 null 无穷远处的准确波形？
• RQ5 标量场的晚期尾行为如何？该格式是否能正确捕捉这种幂律衰减？

主要发现

• 所提出的DG方法由于在数值通量中正确处理了奇异源项，即使在狄拉克δ源位置也实现了全局谱精度。
• 与精确解的收敛性测试验证了指数收敛速率，证实了该方法的高阶精度。
• 该格式正确捕捉了标量场的晚期幂律尾行为，与分析预期一致。
• 双曲层方法成功实现了未来 null 无穷远处波形的直接计算，且边界条件极为简单。
• 一种首阶约化变量的选择导致数值不稳定，约4–5位有效数字后精度损失，而另一种选择则保持了完整的收敛性。
• 利用该方法计算的圆轨道能量辐射率与已知结果一致，证实了物理一致性。