[論文レビュー] A Way to 3D Numerical Relativity --- Coalescing Binary Neutron Stars ---

Takashi Nakamura, Ken-ichi Oohara|ArXiv.org|Dec 16, 1998

Pulsars and Gravitational Waves Research被引用数 4

ひとこと要約

本論文は、共鳴する二重中性子星の3次元数値相対性理論の基盤的アプローチを提示する。アインシュタイン方程式の(3+1)形式を用い、共形的時間スライシングと擬最小歪みゲージ条件を採用する。初期のテストシミュレーションにより、最終ミリ秒間の合体過程をモデル化するためのフレームワークの実現可能性が裏付けられ、これは重力波の主要な源である。ベクトル並列スーパーコンピュータリソースを活用して、完全相対論的かつ非軸対称的力学の計算的要件に対応する。

ABSTRACT

One of the most promising sources of gravitational radiation is coalescence of binary neutron stars or black holes. In order to study gravitational radiation at the merging phase of coalescing binary neutron stars which is called the last three-milliseconds, full general relativistic simulations are required. Since coalescence of binary stars is a completely non-axisymmetric and 3 dimensional (in space) event, which needs great powers of computers, we have to develop a method of 3D numerical relativity using a vector-parallel supercomputer. As a first step of this final goal, we started simulations using post-Newtonian hydrodynamics including radiation reaction of gravitational waves from 1989. They gave a lot of perspectives on coalescing events and gravitational radiation from them. Next we started to attack 3D, fully relativistic simulations. First, basic equations on the (3+1)-formalism of the Einstein equations are shown. As for gauge conditions, we use conformal time slicing and psuedo minimal distortion condtions at present. For these conditions as well as to solve initial value equations, we should solve some elliptic partial differential equations. It consumes the greatest part of CPU time. We show recent results of test simulations on coalescing binary neutron stars.

研究の動機と目的

• 最終ミリ秒間の合体する二重中性子星をシミュレートできる3次元数値相対性理論フレームワークの開発。これは重力波検出にとって重要な段階である。
• 非軸対称的かつ完全相対論的力学に起因する計算課題に対処し、高性能なベクトル並列スーパーコンピュータを必要とする。
• まずゲージ条件や楕円型方程式ソルバーといった主要な要素をテストすることで、将来の完全一般相対論的シミュレーションの基盤を確立する。
• 合体段階からの重力放射の正確なモデル化を可能にする。これは地上重力波検出器の主な標的である。
• 本格的動的進化の前に、テストシミュレーションを通じて手法の妥当性を検証する。

提案手法

• 3次元数値相対性理論のための中心的数学的フレームワークとして、アインシュタイン方程式の(3+1)形式が採用される。
• 座標特異性の制御と数値的安定性の維持のため、共形的時間スライシングと擬最小歪みゲージ条件が実装される。
• 初期データとゲージ条件のための楕円型偏微分方程式が解かれる。これはシミュレーションの計算負荷が最も高い部分を占める。
• 合体の力学的挙動を理解するための初期ステップとして、重力波放射反動を伴うポストニュートン近似流体力学が使用される。
• 3次元相対論的流体力学の高い計算負荷に対応するため、シミュレーションはベクトル並列スーパーコンピュータへの移植を設計されている。
• 本格的動的進化の前に、数値的手法とゲージ条件の妥当性を検証するためのテストシミュレーションが実施される。

実験結果

リサーチクエスチョン

• RQ1 一般相対性理論の(3+1)形式を、合体する二重中性子星の3次元シミュレーションに効果的に実装する方法は何か？
• RQ2 極めて動的な合体段階において、数値的安定性を維持するのに最も適したゲージ条件は何か？
• RQ3 初期データとゲージ条件を支配する楕円型方程式の解法プロセスは、3次元でどの程度計算的に高負荷であるか？
• RQ4 ポストニュートン近似流体力学と放射反動を用いたテストシミュレーションは、最終合体の力学的挙動について信頼できる洞察を提供できるか？
• RQ5 中性子星合体の完全3次元数値相対性理論を達成する際の主な計算ボトルネックは何か？

主な発見

• (3+1)形式に共形的時間スライシングと擬最小歪みゲージ条件を組み合わせたフレームワークは、二重中性子星の3次元数値相対性理論シミュレーションにおいて実現可能であることが示された。
• 初期データとゲージ条件を支配する楕円型偏微分方程式の解法が、CPU時間の大部分を占め、主要な計算ボトルネックであることが明らかになった。
• テストシミュレーションにより、非軸対称的かつ完全相対論的力学をモデル化するための提案された数値的手法の実現可能性が成功裏に示された。
• ポストニュートン近似流体力学に重力波放射反動を組み合わせた手法により、合体の力学的挙動と波動放出に関する貴重な初期の洞察が得られた。
• フレームワークはスケーラブルであり、将来の本格的動的シミュレーションのためにベクトル並列スーパーコンピュータアーキテクチャを活用するように設計されている。
• 本手法は、二重中性子星合体の最終ミリ秒間の完全一般相対論的シミュレーションのためのしっかりとした基盤を築いた。