[論文レビュー] Coalescing neutron stars -- gravitational waves from polytropic models

M. Ruffert, Markus Rampp+1 more|arXiv (Cornell University)|Nov 7, 1996

Pulsars and Gravitational Waves Research被引用数 5

ひとこと要約

本稿は、3次元流体力学的シミュレーションにおいて多項式状態方程式を用いて、合体する中性子星からの重力波放射を調査する。異なる数値的手法と状態方程式の間で結果を比較したところ、多項式モデルはより複雑な状態方程式と比較して20%以内の精度で重力波シグネチャを再現したが、外層の力学的挙動には顕著な差が生じ、質量損失やディスク形成の予測に影響を及ぼした。

ABSTRACT

The dynamics, time evolution of the mass distribution, and gravitational wave signature of coalescing neutron stars described by polytropes are compared with three simulations published previously: (a) ``Run 2'' of Zhuge et al. (1994), (b) ``Model III'' of Shibata et al. (1992), and (c) ``Model A64'' of Ruffert et al. (1996). We aim at studying the differences due to the use of different numerical methods, different implementations of the gravitational wave backreaction, and different equations of state. Comparison (a) confronts the results of our grid-based PPM scheme with those from an SPH code. We found that due to the lower numerical viscosity of the PPM code, the post-merging oscillations and pulsations could be followed for a longer time and lead to larger secondary and tertiary maxima of the gravitational wave luminosity. In case (b) two grid based codes with the same backreaction formalism but differing hydrodynamic integrators and different numerical resolution are compared. Satisfactory agreement of the amplitude of the gravitational wave luminosity is established, although due to the different initial conditions a small time delay develops in the onset of the dynamical instability. In (c) we find that using a polytropic equation of state instead of the high-density equation of state of Lattimer & Swesty (1991) does not change the overall dynamical evolution of the merger and yields agreement of the gravitational wave signature to within 20% accuracy. However, differences of the structure and evolution of the outer layers of the neutron stars are present, which has important implications for questions like mass loss and disk formation during the merging of binary neutron stars.

研究の動機と目的

• 中性子星合体時の重力波放射に及ぼす異なる数値的手法および状態方程式の影響を調査すること。
• 二重中性子星合体のシミュレーションにおいて、多項式状態方程式を近似として用いる妥当性を評価すること。
• グリッドベースのPPMとSPHコードのシミュレーション結果を比較し、合体後の振動および波の放射度に注目すること。
• 数値的分解能および重力波の逆作用の実装方法が重力波シグネチャに与える影響を評価すること。
• 状態方程式の違いが、二重合体における質量噴出およびディスク形成に与える影響を特定すること。

提案手法

• 3次元グリッドベースの区分的放物線法（PPM）を用いて、多項式状態方程式を伴う二重中性子星合体をシミュレートする。
• 3つの先行シミュレーション（Zhuge et al., 1994年、SPHによるRun 2、Shibata et al., 1992年、グリッドベースのModel III、Ruffert et al., 1996年、グリッドベースのModel A64）と結果を比較する。
• グリッドベースの比較において、一貫した重力波逆作用形式を適用し、流体力学的ソルバーおよび分解能の影響を分離する。
• 同じ数値フレームワークを用いて、多項式状態方程式と高密度のLattimer-Swesty状態方程式を対比する。
• 質量分布および重力波放射度の時間発展を分析し、合体後のダイナミクスの差を評価する。
• シミュレーション間で波の振幅、振動周期、二次的・三次的ピークの差を定量化する。

実験結果

リサーチクエスチョン

• RQ1 異なる数値的手法（PPM 対 SPH）は、合体後の振動および重力波放射度のシミュレーションにどのように影響を与えるか？
• RQ2 流体力学的ソルバーおよび分解能の違いが、グリッドベースのシミュレーションにおける重力波シグネチャに及ぼす影響はどの程度か？
• RQ3 多項式状態方程式は、より現実的な状態方程式の重力波放射をどの程度正確に再現できるか？
• RQ4 状態方程式の違いが、中性子星合体における質量損失およびディスク形成に与える影響は何か？
• RQ5 初期条件および数値的粘性は、重力波放射のタイミングと振幅にどのように影響するか？

主な発見

• PPMコードはSPHコードよりも低い数値的粘性を示し、合体後の振動のシミュレーション期間が長く、重力波放射度の二次的・三次的ピークも大きくなった。
• 同じ逆作用形式を用いた2つのグリッドベースコード間では、重力波放射度の振幅に良好な一致が得られたが、初期条件の違いによりわずかな時間遅れが生じた。
• 多項式状態方程式を用いることで、より複雑なLattimer-Swesty状態方程式と比較して20%以内の精度で重力波シグネチャが再現され、全体的なダイナミクスに良好な一致が得られた。
• 全体的なダイナミクスは類似しているものの、異なる状態方程式下での中性子星の外層の構造および進化に顕著な差が生じ、質量噴出およびディスク形成の予測に影響を及ぼした。
• 結果から、多項式モデルは重力波の主要な特徴を捉えるには十分であるが、核合成および電磁的対応天体に関連する微細なダイナミクスを完全に解明できるとは限らないことが示唆された。
• 本研究は、数値的手法および分解能が、特に合体後の段階において波形の詳細に顕著な影響を及ぼすことを確認した。