[論文レビュー] A fast apparent horizon algorithm

Erik Schnetter|ArXiv.org|Jun 2, 2002

Astrophysical Phenomena and Observations被引用数 7

ひとこと要約

本稿では、テンソル量をデカルト座標表現を用いて非線形楕円型偏微分方程式として解くことで、数値相対性力学における見かけの視界を高速かつ安定して特定するアルゴリズムを提示する。この手法は座標特異性を回避し、任意の分解能と形状の柔軟性を実現でき、単一および複数のブラックホール系、特に合体状態を含む共通の視界を的確に特定する。

ABSTRACT

I present a fast algorithm to find apparent horizons. This algorithm uses an explicit representation of the horizon surface, allowing for arbitrary horizon resolutions and, in principle, shapes. Novel in this approach is that the tensor quantities describing the horizon live directly on the horizon surface, yet are represented using Cartesian coordinate components. This eliminates coordinate singularities, and leads to an efficient implementation. The apparent horizon equation is then solved as a nonlinear elliptic equation with standard methods. I explain in detail the coordinate systems used to store and represent the tensor components of the intermediate quantities, and describe the grid boundary conditions and the treatment of the polar coordinate singularities. Last I give as examples apparent horizons for single and multiple black hole configurations.

研究の動機と目的

• 数値相対性力学の進化コードに統合可能な高速かつ安定した見かけの視界検出法の開発を目的とする。
• 球対称座標における座標特異性の問題を解消するため、テンソル成分をデカルト座標で表現することを目的とする。
• 視界表面の表現における任意の分解能および形状の柔軟性を実現することを目的とする。
• 分離視界の性質を用いてブラックホールの質量およびスピンを正確に推定することを目的とする。
• ブラックホールシミュレーションにおける除去境界条件を支援するため、視界境界を正確に特定することを目的とする。

提案手法

• 見かけの視界は、球座標系における表面をパラメータ表示する関数 r = h(θ, φ) として明示的に表現する。
• 外向き法線や外在的曲率などのテンソル量は、極座標特異性を避けるためにデカルト成分で表現する。
• 視界方程式 H = 0 を、h(θ, φ) に関する非線形楕円型偏微分方程式として定式化し、標準的な数値手法で解く。
• グリッド境界条件および極での座標特異性の取り扱いは、実装において明示的に扱われる。
• 視界表面を定義し、法線ベクトルを計算するために、レベルセット関数 F(r, θ, φ) = r - h(θ, φ) を用いる。
• 本手法は、単一および複数のブラックホール系、特に重ね合わせたKerr–Schildデータを含む構成に対してテストされた。

実験結果

リサーチクエスチョン

• RQ1 球座標における座標特異性を回避できる高速かつ安定な見かけの視界検出法を開発することは可能か？
• RQ2 デカルト座標を用いて、視界表面上のテンソル成分を効率的かつ正確に表現する方法は何か？
• RQ3 標準的な数値手法を用いて、数値的に安定に、楕円型視界方程式をどの程度まで安定して解くことができるか？
• RQ4 二重ブラックホール系において、共通の見かけの視界が形成されたり消失したりする条件は何か？
• RQ5 ADM質量および視界質量推定値における数値誤差は、結合エネルギーおよび放射エネルギーの解釈にどのような影響を及えるか？

主な発見

• アルゴリズムは、重ね合わせたブラックホール構成において共通の見かけの視界を的確に特定でき、1つの合体ブラックホールの存在を示している。
• z = ±2.0 において、共通の視界が形成され、視界質量 M_C が個々の質量の2倍未満であることが判明し、結合エネルギーが存在することを示している。
• z = ±1.75 において、個々の視界は重なり合うが、共通の視界は明確に分離されており、わずかに歪んでいるにとどまる。
• z = ±3.0 において、共通の視界は消失し、視界方程式の残差ノルムは 0.0248 のまま維持され、収束していないが誤差は小さいことが示された。
• 結合エネルギー E_B = M_C - 2M_I は負であり、分離距離が増加するにつれて減少し、物理的期待と整合的である。
• ADM質量の計算では、E_R = M_ADM - M_C に対して物理的に意味のない負の値が得られ、外周境界付近でのグローバル質量推定に数値誤差が生じていることが示唆された。