[論文レビュー] 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

ひとこと要約

本稿では、Kerr時空におけるスカラーTeukolsky方程式の多領域不連続ガラーキン法を提案する。球面調和関数展開と1階対称双曲型形式を用いることで、特異なDiracデルタ源項に対してもグローバルなスペクトル的精度を達成する。この手法により、正規化に起因する誤差を回避し、極端質量比インスピレーション（EMRI）の指数的収束シミュレーションが可能となり、未来の光円錐無限遠で正確な波形が得られ、点粒子特異性の適切な取り扱いも可能となる。

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.

研究の動機と目的

• 極端質量比インスピレーション（EMRI）に由来する特異な源項を伴う時間領域におけるスカラーTeukolsky方程式の解法を開発すること。
• 正規化に起因する誤差を回避するため、Diracデルタ関数を直接数値スキームに組み込むことにより、近似による誤差を排除すること。
• 不連続ガラーキン法を用い、適切なフラックス処理を施すことにより、点粒子源の位置でもグローバルなスペクトル的精度を達成すること。
• 超仮想層法を用いて未来の光円錐無限遠を有限境界に写像し、人工的境界条件を回避して直接遠方波形を計算すること。

提案手法

• スカラー場を球面調和関数で展開することで、3+1次元Teukolsky方程式を座標τ（トロイズ座標）と時間の1+1次元系に簡略化する。
• 変数πₗₘ = -∂ψₗₘ/∂τおよびφₗₘ = ∂ψₗₘ/∂ρを用いて、2階波動方程式を1階対称双曲型系に再定式化する。
• 多領域不連続ガラーキン法を適用し、数値フラックスを特に設計して、Diracデルタ源項を正しく取り扱う。
• 超仮想層法を用いて、未来の光円錐無限遠を有限境界に写像し、遠方波形への直接アクセスを可能にするとともに、外境界条件を単純化する。
• 強双曲型性と安定性を保証するため、右端境界（ρ = s）における波速の取り扱いに注意を払う。特に、1つの波速が0に近づく点に注意を要する。
• 2種類の1階還元変数の選択を比較し、理論的には同等であるが、一方の選択は数値的不安定性を示すことが判明した。

実験結果

リサーチクエスチョン

• RQ1 不連続ガラーキン法は、特異点付近でもDiracデルタ源項を伴うTeukolsky方程式に対して、グローバルなスペクトル収束を達成できるか？
• RQ2 適切なフラックス処理を施した対称双曲型形式は、δ関数の近似に起因する正規化誤差を完全に排除できるか？
• RQ3 1階還元変数の選択が、点粒子源の存在下での数値的安定性および収束性にどのように影響するか？
• RQ4 超仮想層法をDG法と効果的に組み合わせ、人工的境界条件を用いずに未来の光円錐無限遠での正確な波形を抽出できるか？
• RQ5 スカラー場の遅延尾（late-time tail）の振る舞いはどのようなものか？また、このスキームはそのべき乗則的減衰を正しく捉えているか？

主な発見

• 正しく設計された数値フラックスにより、Diracデルタ源の特異性を適切に処理したため、本提案DG法は特異点付近でもグローバルなスペクトル的精度を達成した。
• 正確な解との収束テストにより、指数的収束率が確認され、本手法の高次精度が裏付けられた。
• 本スキームは、解析的予想と整合的な、スカラー場の遅延尾のべき乗則的減衰を正しく捉えた。
• 超仮想層法により、未来の光円錐無限遠での波形を直接計算可能となり、境界条件が単純化された。
• 1つの1階還元変数の選択では、約4〜5桁の精度を過ぎると数値的不安定性が生じ、収束性が損なわれる一方、他方の選択では完全な収束性が維持された。
• 本手法を用いて円軌道のエネルギー放射率を計算した結果、既知の結果と整合しており、物理的整合性が確認された。