[論文レビュー] Event and Apparent Horizon Finders for 3+1 Numerical Relativity

Jonathan Thornburg|ePrints Soton (University of Southampton)|Dec 30, 2005

Astrophysical Phenomena and Observations参考文献 67被引用数 5

ひとこと要約

この論文は、3+1次元数値相対性理論シミュレーションにおけるイベントホライズンおよびありったけホライズンを特定するための数値的アルゴリズムおよびソフトウェアツールについて包括的なレビューを提供している。イベントホライズンは光線の後退的統合または表面の後退的統合により特定され、ありったけホライズンは限界的外方捕獲表面を表す非線形楕円型偏微分方程式（PDE）を解くことで特定される。両者について、軸対称および一般の時空において、効率性、正確性、耐障害性の観点から強調されている。

ABSTRACT

Event and apparent horizons are key diagnostics for the presence and properties of black holes. In this article I review numerical algorithms and codes for finding event and apparent horizons in numerically-computed spacetimes, focusing on calculations done using the 3+1 ADM formalism. There are 3 basic algorithms for finding event horizons, based respectively on integrating null geodesics \emph{forwards} in time, integrating null geodesics \emph{backwards} in time, and integrating null \emph{surfaces} backwards in time. The last of these is generally the most efficient and accurate. There are a large number of apparent-horizon finding algorithms, with differing trade-offs between speed, robustness, accuracy, and ease of programming. In axisymmetry, shooting algorithms work well and are fairly easy to program. In slices with no continuous symmetries, Nakamura et al.'s algorithm and elliptic-PDE algorithms are fast and accurate, but require good initial guesses to converge. In many cases Schnetter's "pretracking" algorithm can greatly improve an elliptic-PDE algorithm's robustness. Flow algorithms are generally quite slow, but can be very robust in their convergence.

研究の動機と目的

• 3+1 ADM形式による時空計算で得られるイベントホライズンおよびありったけホライズンを検出するための数値的アルゴリズムを体系化し、レビューすること。
• 速度、正確性、耐障害性、実装の容易さという観点から、さまざまなホライズン特定手法の長所と短所を特定・比較すること。
• 時空の対称性（例：軸対称性対一般）および計算制約に応じて適切なアルゴリズムを選択するための指針を提供すること。
• 公開可能なホライズン特定ソフトウェアのリスト作成と、数値相対性理論におけるソフトウェア工学のベストプラクティスの強調を通じて、コードの再利用と共有を促進すること。
• 新しい手法が登場するにつれてレビューを更新・拡張できるようにコミュニティの貢献を促すことで、研究者にとっての動的リファレンスとしての役割を果たすこと。

提案手法

• イベントホライズンは、主に後退的統合による3つの主要なアルゴリズム（光線の前向き統合、光線の後退的統合、光の表面の後退的統合）を用いて後処理で特定される。後退的統合による表面の方法が最も効率的かつ正確である。
• ありったけホライズンは、限界的外方捕獲表面（MOTS）条件を解くことで特定され、これは3次元計量、その空間的微分、および外的曲率を含む非線形楕円型偏微分方程式（PDE）である。
• 軸対称性では、シューティング法およびゼロ探索法が効果的で、実装が簡単である。
• 一般の時空では、中村らのスペクトル的統合反復法および楕円型PDEソルバーが高速かつ正確であるが、良好な初期推定値が必要である。
• ホライズンの事前追跡は、楕円型PDEソルバーの収束を改善するために用いられ、試行的ホライズンを真の解へと進化させる。
• 流れアルゴリズムは、速度よりも耐障害性を重視する場合に用いられ、動的進化によってMOTSに収束する。

実験結果

リサーチクエスチョン

• RQ1 数値的に計算された時空におけるイベントホライズンを特定するための最も効率的かつ正確な数値的手法は何か？
• RQ2 さまざまな時空の対称性（軸対称性、一般）において、異なるありったけホライズン特定アルゴリズムの収束性、耐障害性、計算コストはどのように比較できるか？
• RQ3 楕円型PDEに基づくありったけホライズン特定器の収束に、初期推定値の選択が果たす役割は何か？また、どのように改善できるか？
• RQ4 ソフトウェア工学の実践は、数値相対性理論におけるホライズン特定コードの信頼性、再利用性、保守性をどのように向上させられるか？
• RQ5 非線形方程式および常微分方程式（ODE）のための公開可能な数値ライブラリやソルバー（例：netlib）は、どのようにして堅牢なホライズン特定ツールの構築に活用できるか？

主な発見

• 光の表面の後退的統合は、イベントホライズンを特定するための最も効率的かつ正確な手法であり、前向き統合や光線ベースの手法を上回る。
• 楕円型PDEソルバーや中村らのスペクトル法は、一般時空におけるありったけホライズン検出において高い正確性と速度を発揮するが、収束のためには良好な初期推定値が必要である。
• シュナイターの事前追跡アルゴリズムは、楕円型PDEホライズン特定器の耐障害性を顕著に向上させ、解への収束を導く。
• 流れアルゴリズムは収束性に優れているが、計算コストが高いため、速度よりも信頼性が優先される状況に適している。
• 非線形代数方程式を解くために、van Wijngaarden-Dekker-Brent法が推奨される。これは速度と信頼性のバランスが取れているためである。
• LSODE や RKF45 といった堅牢な ODE ソルバーは、ホライズン特定パイプラインにおける時間積分タスクに不可欠であり、netlib を通じて自由に利用可能である。これにより、信頼性と効率性の両立が実現できる。