[論文レビュー] A Numerical Study of Boson Star Binaries

Bruno C. Mundim|arXiv (Cornell University)|Mar 1, 2010

Pulsars and Gravitational Waves Research参考文献 156被引用数 8

ひとこと要約

本稿では、3+1次元における4階放物型および5階楕円型偏微分方程式の連立系として記述されるアインシュタイン＝クライン＝ゴルドン系を解くために、共形平坦近似（CFA）と最大スライシングを用いたボソン星二重星の数値的検討が行われる。シミュレーションは連続解への収束を達成し、2周期を超える安定した軌道的進化を示し、合体過程においてブラックホール形成の兆候も観測された。これは、コンパクトなボソン星のダイナミクスをシミュレートする上で顕著な前進を示している。

ABSTRACT

This thesis describes a numerical study of binary boson stars within the context of an approximation to general relativity. The approximation we adopt places certain restrictions on the dynamical variables of general relativity (conformal flatness of the 3-metric), and on the time-slicing of the spacetime (maximal slicing). The resulting modeling problem requires the solution of a coupled nonlinear system of 4 hyperbolic, and 5 elliptic partial differential equations (PDEs) in three space dimensions and time. We approximately solve this system as an initial-boundary value problem, using finite difference techniques and well known, computationally efficient numerical algorithms such as the multigrid method in the case of the elliptic equations. Careful attention is paid to the issue of code validation, and a key part of the thesis is the demonstration that, as the basic scale of finite difference discretization is reduced, our numerical code generates results that converge to a solution of the continuum system of PDEs as desired. The thesis concludes with a discussion of results from some initial explorations of the orbital dynamics of boson star binaries. In particular, we describe calculations in which motion of such a binary is followed for more than two orbital periods, which is a significant advance over previous studies. We also present results from computations in which the boson stars merge, and where there is evidence for black hole formation.

研究の動機と目的

• 一般相対性理論の簡略化されたが物理的に妥当な近似において、重力的にコンパクトなボソン星二重星のダイナミクスをモデル化すること。
• ボソン星の進化を支配する結合非線形偏微分方程式系を解くことのできる数値コードの開発と検証を行うこと。
• ボソン星二重星の軌道的進化と合併結果を調査し、ブラックホール形成の可能性を検討すること。
• 有限差分離散化の系統的精錬を通じて、数値的収束性とコードの信頼性を示すこと。

提案手法

• 3次元計量を簡略化し、アインシュタイン方程式の複雑さを低減するため、共形平坦近似（CFA）を採用する。
• 時間ゲージを固定し、進化過程における座標特異点を回避するため、最大スライシングを適用する。
• 3次元空間と時間における4階放物型および5階楕円型偏微分方程式系を、有限差分法（FDA）で離散化する。
• 制約系から生じる楕円型方程式を効率的に解くために、多重グリッド法を用いる。
• 各タイムステップにおける非線形代数的系を解くために、ニュートン・ガウス＝ザイデル反復法（NGS）を実装する。
• グリッド解像度を高めるに従って数値解が連続極限に収束することを示すことにより、コードの妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 CFAおよび最大スライシングに基づく数値コードは、ボソン星二重星の長期的軌道的進化を正確にシミュレートできるか？
• RQ2 インスピレーションおよび合併過程におけるボソン星二重星のダイナミクス的挙動は何か。特に、ありうる視界の形成は？
• RQ3 有限差分グリッドを精錬するに従って、系は数値的収束を示すか。これはシミュレーションの信頼性を裏付けるか？
• RQ4 ボソン星の合併の結果は何か。具体的には、ブラックホール形成の兆候は見られるか？
• RQ5 合併過程中、スカラー場および時空幾何学の性質はどのように変化するか？

主な発見

• 有限差分グリッド間隔を小さくするに従って、数値コードは連続解への一貫した収束を示し、コードの信頼性と正確性が確認された。
• 本シミュレーションは、ボソン星二重星の軌道運動を2周期を超えて安定に追跡できた。これは、先行研究に比べ顕著な向上を示している。
• 合併段階において、高密度コアが形成され、これはありうる視界の形成を示唆しており、ブラックホール形成の兆候である。
• スカラー場および計量成分の進化は、シミュレーション全体を通して安定しており、数値的爆発や不安定性の兆候は観測されなかった。
• 結果から、ボソン星二重星は合併を経て、コンパクトで視界に類似した構造を形成しうること、したがって重力波源としての妥当性が支持されることが示唆された。