[論文レビュー] Hyperboloidal discontinuous time-symmetric numerical algorithm with higher order jumps for gravitational self-force computations in the time domain

Lidia J. Gomes Da Silva, Rodrigo Panosso Macedo+4 more|arXiv (Cornell University)|Jun 22, 2023

Radio Astronomy Observations and Technology被引用数 4

ひとこと要約

本稿では、特異性を持つ点粒子における不連続性を考慮した高次跳躍を有する双曲的非連続時間対称数値アルゴリズムを提案し、時間領域における重力自己力計算の高次精度を達成する。不連続ラグランジュおよびエルミート補間を活用することで、粒子付近での特異性にもかかわらず高次精度を実現する。本手法により、放射フラックスの正確な測定と、粒子極限における場および微分の評価が可能となり、極端質量比降下（EMRI）における自己一貫的軌道進化に適していることが示された。

ABSTRACT

Within the next decade the Laser Interferometer Space Antenna (LISA) is due to be launched, providing the opportunity to extract physics from stellar objects and systems, such as extit{Extreme Mass Ratio Inspirals}, (EMRIs) otherwise undetectable to ground based interferometers and Pulsar Timing Arrays (PTA). Unlike previous sources detected by the currently available observational methods, these sources can extit{only} be simulated using an accurate computation of the gravitational self-force. Whereas the field has seen outstanding progress in the frequency domain, metric reconstruction and self-force calculations are still an open challenge in the time domain. Such computations would not only further corroborate frequency domain calculations and models, but also allow for full self-consistent evolution of the orbit under the effect of the self-force. Given we have extit{a priori} information about the local structure of the discontinuity at the particle, we will show how to construct discontinuous spatial and temporal discretisations by operating on discontinuous Lagrange and Hermite interpolation formulae and hence recover higher order accuracy. In this work we demonstrate how this technique in conjunction with well-suited gauge choice (hyperboloidal slicing) and numerical (discontinuous collocation with time symmetric) methods can provide a relatively simple method of lines numerical algorithm to the problem. This is the first of a series of papers studying the behaviour of a point-particle prescribing circular geodesic motion in Schwarzschild in the extit{time domain}. In this work we describe the numerical machinery necessary for these computations and show not only our work is capable of highly accurate flux radiation measurements but it also shows suitability for evaluation of the necessary field and it's derivatives at the particle limit.

研究の動機と目的

• LISAミッション科学にとって不可欠な極端質量比降下（EMRI）における重力自己力計算の、堅牢な時間領域数値フレームワークの構築を目的とする。
• 周波数領域での進展にもかかわらず、計量再構成と粒子位置における場の評価が未解決の時間領域での自己力の正確な計算に挑戦することを目的とする。
• 点粒子における特異性を克服するため、空間的および時間的離散化に不連続性構造に関する事前知識を組み込むこと。
• 双曲的スライシングと時間対称不連続コロケーションを用いた時間方向差分法により、フラックス放射および場の微分測定における高次精度と安定性を示すこと。

提案手法

• 本手法は、未来の光円錐での放射エネルギーを自然に吸収するため、双曲的スライシングを採用し、時間領域シミュレーションにおける長期的安定性を向上させる。
• 空間的および時間的離散化には、粒子における不連続性の局所的構造を捉えるために不連続ラグランジュおよびエルミート補間式を用いる。
• 場の特異的挙動に基づくジャンプ条件を補間フレームワークに組み込むことで、高次精度を回復する。
• 時間対称不連続コロケーション法を、時間方向差分法の定式化に適用し、安定性および収束性を向上させる。
• 本アルゴリズムは、シュワルツシルト時空における円軌道運動を行う点粒子を対象として実装され、EMRI波形の基礎的モデルとして機能する。
• 数値的妥当性は、[39]の周波数領域データと比較して、RWZマスタ関数およびその微分の時間領域結果を検証し、式 (76) を用いて誤差推定値を算出することで確認された。

実験結果

リサーチクエスチョン

• RQ1 点粒子特異性が存在するにもかかわらず、時間領域数値アルゴリズムが重力自己力計算において高次精度を達成できるか。
• RQ2 場が粒子位置でジャンプを示す場合に、不連続な空間的および時間的離散化をどのように構築すれば、精度を保持できるか。
• RQ3 双曲的スライシングは、自己力効果の長期的時間領域シミュレーションにおける安定性および精度をどの程度向上させるか。
• RQ4 本手法は、フラックスおよび場の微分を粒子位置で正確に計算できるか。これにより、自己一貫的軌道進化が可能になるか。
• RQ5 RWZマスタ関数に関して、周波数領域結果と比較して、収束性および数値誤差の観点から、本アルゴリズムの性能はどの程度か。

主な発見

• RWZマスタ関数の内部および外部解において、本アルゴリズムは高次収束を達成し、特定のモードでは数値誤差が $1.7 \times 10^{-9}$ まで低下した。
• フラックス放射測定は非常に高精度であり、最初の $l=5$ モードでは誤差が $[1.7 \times 10^{-9}, 7.9 \times 10^{-8}]$ の範囲に収束した。
• 本手法は、粒子位置におけるRWZ関数の一次時間および径方向微分を正確に計算でき、大多数のモードで誤差が $10^{-8}$ 未満であった。
• $ au=10,000$ における $(l,m)=(2,2)$ モードの位相図では、シンプレクティック構造の保存が観察され、長期的安定性が裏付けられた。
• 数値収束解析により、最適な時間離散化ステップ選択が確認され、誤差のスケーリングは理論的期待と整合的であった。
• 本手法は、自己力の粒子位置における評価と整合性を示し、EMRIモデルにおける自己一貫的軌道進化に不可欠な要件を満たしている。