[論文レビュー] Variational Principles in General Relativity

Brian D. Baker|ArXiv.org|May 26, 2002

Pulsars and Gravitational Waves Research参考文献 35被引用数 4

ひとこと要約

本論文は、一般相対性理論における定常的で軸対称的かつ渦なしの二重中性子星系に対して、3+1形式と表面積分を用いてエネルギー・質量の一貫性を確保する変分原理を定式化する。非線形適応多重グリッド法を用いて自己無矛盾な解を導出し、数値的妥当性を後期ニュートン近似理論と照合することで、高精度な重力波および多極モーメントの予測を達成した。

ABSTRACT

This thesis details an effort to generate astrophysically interesting solutions to the two-body problem in General Relativity. The thesis consists of two main parts. The first part presents an analytical variational principle for describing binary neutron stars undergoing irrotational fluid flow. The variational principle is a powerful tool for generating accurate estimates of orbital parameters for neutron stars in quasi-equilibrium orbits. The second part of the thesis details the numerical application of the variational principle by solving the initial value problem for binary black holes in quasi- equilibrium circular orbits. The analysis draws from the novel ``puncture'' method of describing the black holes, and relies on nonlinear adaptive multigrid techniques for generating numerical results. We arrive at two important conclusions. First, the analytical variational principle describing binary neutron stars in irrotational motion provides a road map for future numerical simulations, and also lends credence to previous simulations by other authors. Second, the numerical application and description of binary black holes in quasi-equilibrium circular orbits simplifies the analyses of previous authors, and allows for the imposition of realistic boundary data in simulations with relatively high grid densities. Both the variational principle and its application are used to generate accurate estimates of the orbital parameters.

研究の動機と目的

• 一般相対性理論における定常的で軸対称的かつ渦なしの二重中性子星系に対して、変分原理を構築すること。
• ADM質量とエネルギー保存則を正しく保つために、表面積分と境界条件を組み込むこと。
• 非線形適応多重グリッドスolv erを実装し、結合されたアインシュタイン方程式を数値的に解くこと。
• 数値的結果を後期ニュートン近似と照合し、多極モーメントの収束性を検証すること。
• 高精度な重力波および質量モーメント予測が可能な、自己無矛盾なフレームワークを提供すること。

提案手法

• 時空を3+1分解することで、アインシュタイン方程式を初期値問題と発展問題に分離する。
• 対称トレースレステンソルとターネの一般化解を用いて、流体の運動と計量の摂動をモデル化する。
• 変分原理から生じる楕円型方程式を解くために、非線形適応多重グリッドアルゴリズムを実装する。
• グリッドの細分化と補間・縮小作用素（例：symrstrct, rstrct）を用いて、マルチスケール解を扱う。
• ガラーキン有限差分法とロビン境界条件を用いて、偏微分方程式の数値解を安定化する。
• 計算グリッド上での数値的統合によりADM質量と多極モーメントを計算し、グローバル保存則の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 一般相対性理論において、定常的で軸対称的かつ渦なしの二重中性子星系に対して、一貫性のある変分原理をどのように定式化できるか？
• RQ2 表面積分と境界項は、正しいADM質量とエネルギー保存則を保つために果たす役割は何か？
• RQ3 非線形適応多重グリッド法は、コンパクトな二重星系のアインシュタイン方程式を解く際の収束性と精度をどのように向上させるか？
• RQ4 数値的結果は、 lapse 関数および重力波成分について、後期ニュートン近似とどの程度一致するか？
• RQ5 この手法は、平衡状態における多極モーメントおよび質量エネルギー分布を高精度に計算できるか？

主な発見

• 変分原理は、正しい表面および境界条件を満たす自己無矛盾な二重中性子星の平衡状態を効果的に生成した。
• 数値的結果は、lapse 関数および重力波振幅について、後期ニュートン近似理論と良好に一致し、手法の正確性が妥当であることが検証された。
• 適応多重グリッドソルバーは効率的に収束し、複数のグリッドレベルで残差低減と誤差制御が確認された。
• 多極モーメントとADM質量は高精度に計算され、グローバルエネルギーおよび運動量保存則の妥当性が確認された。
• 対称トレースレステンソルと表面積分項の実装により、一般相対性理論の制約と整合性が保たれた。
• コードは、強力場領域でさえも最小限の数値的アーティファクトを伴って、重力ポテンシャルおよび計量成分を的確に計算できた。