[論文レビュー] General Relativistic Numerical Simulation on Coalescing Binary Neutron Stars and Gauge-Invariant Gravitational Wave Extraction

Mari Kawamura, K. Oohara+1 more|arXiv (Cornell University)|Jun 24, 2003

Pulsars and Gravitational Waves Research被引用数 8

ひとこと要約

本論文は、最大スライシングとゲージ不変な重力波抽出を用いた、合体する二重中性子星の3次元一般相対論的数値シミュレーションを提示する。合併後、約25 mHzにおける準正規モード励起が可能であることが示されたが、数値境界における波の反射が精度を制限する。結果は先行研究と定性的に一致するが、非平衡初期条件のため放射エネルギーは低くなる。

ABSTRACT

We are developing 3 dimensional simulation codes for coalescing binary neutron stars. A code using the maximal slicing condition is obtained. To evaluate the gravitational radiation, we implemented a gauge-invariant wave extraction and compared the wave forms with the metric tensors at the wave zone. The energy spectrum of the waves was also evaluated to investigate the possibility that the excitation of the quasi-normal modes of the black hole, which may be formed after the merger of two stars, can be observed.

研究の動機と目的

• 共形スライシングで見られる不安定性を回避するため、最大スライシングを用いた3次元一般相対論的シミュレーションコードの開発。
• 座標選択に依存しない物理的波形を保証するため、ゲージ不変な重力波抽出の実装。
• 放射エネルギースペクトルの評価と、合併後に形成されるブラックホールにおける準正規モード励起の可能性の調査。
• シミュレーション結果を、特に塩田と浦川の研究と比較し、コードの妥当性を検証するとともに、放射エネルギーと角運動量の評価。

提案手法

• 一般相対論の(3+1)形式を採用し、アインシュタイン方程式を12個の発展方程式と4個の制約方程式に分解。
• 長期間の数値安定性を向上させるために、最大スライシング（K=0）を採用し、先行研究で不安定性を示した共形スライシングに代わる。
• アブラハムズらの手法を非軸対称かつ偶関数・奇関数パリティの摂動へ拡張し、ゲージ不変な波の抽出を実装。
• 制約の満たされ方と座標制御を維持するため、ラプス関数α、シフトベクトルβi、共形因子φの楕円型方程式を解く。
• 座標のずれを低減するため、擬似最小歪み条件を採用し、時間経過で∂j(γ̃ij)に約2%の誤差が生じることをモニタリング。
• 475×475×238グリッドを用い、メモリ80 GB、1回の実行あたり約100 CPU時間の計算負荷を要し、主な計算ボトルネックは楕円型方程式の解法。

実験結果

リサーチクエスチョン

• RQ1 最大スライシングを用いた3次元一般相対論的コードは、二重中性子星の合併を正確にシミュレートし、ゲージ不変な重力波形を抽出できるか？
• RQ2 二重中性子星の合併は、形成されるブラックホールの準正規モードを励起するか？その兆候は重力波スペクトルから観測可能か？
• RQ3 非平衡初期条件は、準平衡モデルと比較して放射エネルギーと角運動量にどのような影響を与えるか？
• RQ4 数値境界における波の反射は、重力波信号にどの程度歪みを及けるか？また、その影響を軽減する方法はあるか？

主な発見

• ゲージ不変な波抽出法により、物理的重力波信号が明確に分離され、外部真空中ではゲージ依存モードが抑制された。
• 25 mHz付近に重力波エネルギースペクトルのピークが観測され、2.8 M⊙のシュワルツシルトブラックホールの基本l=2準正規モードと整合的である。
• 放射エネルギーと角運動量は、塩田と浦川のモデルH-1と比較して小さい。これは非平衡初期データが合併を速くするためである。
• ADM質量は1.7 msまで10%以内で保存されたが、合併後は粗いグリッド解像度とブラックホール形成のため保存性が低下した。
• 数値境界での重力波の反射は約2.5 ms以降に顕著になり、境界に近いため反射波が増幅された。
• コードの性能は、ラプス関数α、シフトベクトルβi、共形因子φの楕円型方程式の解法に起因するCPU時間の制限を受けており、今後の改善には高速なソルバーや細かいグリッドが必要である。