[論文レビュー] Kilohertz Gravitational Waves From Binary Neutron Star Mergers: Numerical-relativity Informed Postmerger Model

M. Breschi, Sebastiano Bernuzzi+4 more|arXiv (Cornell University)|May 18, 2022

Pulsars and Gravitational Waves Research被引用数 10

ひとこと要約

本稿では、多様な状態方程式をカバーする618件の数値相対論的シミュレーションを用いてキャリブレーションされた、二つの中性子星合体からの合体後重力波信号を解析的にモデル化するNRPMwを紹介する。複素ガウスウェーブレットと状態方程式に依存しない関係式を用いることで、独立した検証セットにおいて0.9を超えるフィッティングファクターを達成し、第三世代検出器（エイントシュタイン・テレスコープ）向けの高精度なベイズ推論を可能にする。

ABSTRACT

We present ${ t NRPMw}$, an analytical model of gravitational-waves from neutron star merger remnants informed using 618 numerical relativity (NR) simulations. ${ t NRPMw}$ is designed in the frequency domain using a combination of complex Gaussian wavelets. The wavelet's parameters are calibrated to equations of state (EOS) insensitive relations from NR data. The NR simulations are computed with 21 EOS (7 of which are finite-temperature microphysical models, and 3 of which contain quark phase transitions or hyperonic degrees of freedom) and span total binary masses $M\in[2.4,3.4]~{ m M}_\odot$, mass ratios up to $q=2$, and (nonprecessing) dimensionless spins magnitudes up to ${0.2}$. The theoretical uncertainties of the EOS-insensitive relations are incorporated in ${ t NRPMw}$ using recalibration parameters that enhance the flexibility and accuracy of the model. ${ t NRPMw}$ is NR-faithful with fitting factors ${\gtrsim}0.9$ computed on an independent validation set of 102 simulations.

研究の動機と目的

• 第三世代検出器向けに、二つの中性子星合体からの合体後重力波信号を高精度で解析的にモデル化し、ベイズ推論を支援すること。
• 従来のモデルを改善するため、有限温度やクォーク・ハドロン相転移モデルを含む広範な状態方程式をカバーすること。
• 618件の数値相対論的シミュレーションから導出された状態方程式に依存しない関係式を用いて、ウェーブレットパラメータをキャリブレーションし、モデルの柔軟性と精度を向上させること。
• 状態方程式に依存しない関係式における理論的不確実性を再キャリブレーションパラメータで表現し、多様な天体物理学的状況においてもモデルの頑健性を高めること。
• 独立した102件のシミュレーションからなる検証セットを用いてモデルの忠実度を検証し、将来的な重力波データ解析における信頼性を保証すること。

提案手法

• 本モデルは、周波数領域における複素ガウスウェーブレットを用いて合体後重力波信号を表現し、効率的かつ高精度なスペクトル分解を実現する。
• ウェーブレットパラメータは、合体後の質量 M ∈ [2.4, 3.4] M⊙、質量比 q まで2.0、スピンの大きさまで0.2をカバーする618件の数値相対論的シミュレーションから導かれた状態方程式に依存しない関係式に合わせてキャリブレーションされる。
• 質量スケーリングされた量（例：Mf₂(R₁.₈/M)）を用いたキャリブレーションにより、異なる二重星系においても普遍性と一般化能力が向上する。
• 状態方程式に依存しない関係式における理論的不確実性は再キャリブレーションパラメータに埋め込まれ、モデルの柔軟性と頑健性が向上する。
• 独立した102件のシミュレーションからなる検証セットにおいて、0.9を超えるフィッティングファクターが得られ、NRの忠実度が確認された。
• インスピレーション・マージャー・フェーズは有効1体枠組みに基づいて構築され、NRPMwは特に合体後の放射をモデル化する。

Figure 1: Representative examples of BNS PM waveforms. The plot shows the plus polarization $h_{+}(t)$ of the time-domain $\ell=m=2$ waveform (solid line) and the instantaneous GW frequency $\omega_{22}(t)$ (dashed line). The NR simulations are from the CoRe database and computed in Refs. Bernuzzi e

実験結果

リサーチクエスチョン

• RQ1 ウェーブレットに基づく解析的モデルは、二つの中性子星合体からのキロヘルツ帯の合体後重力波信号を正確に表現できるか？
• RQ2 f₂(R₁.₈/M) などの状態方程式に依存しない関係式は、多様な状態方程式において合体後のピーク周波数をどれほど正確に予測できるか？
• RQ3 R₁.₄/R₁.₈ の比を合体後スペクトル特徴のキャリブレーションに追加するパラメータとして含めることで、どのような影響を及ぼすか？
• RQ4 状態方程式に依存しない関係式に理論的不確実性を組み込むことで、モデルの忠実度と頑健性はどの程度向上するか？
• RQ5 独立した検証セットにおいて、モデルが0.9を超えるフィッティングファクターを達成できるか。これにより、ベイズ重力波データ解析に適していることが確認できるか？

主な発見

• NRPMw は、独立した102件の数値相対論的シミュレーションからなる検証セットにおいて、0.9を超えるフィッティングファクターを達成し、NRの忠実度が確認された。
• 最も頑健で信頼性の高い準普遍的関係は Mf₂(R₁.₈/M) であり、これは高密度状態方程式の性質と強く相関している。
• 質量スケーリングされた量（例：Mf₂(R₁.₈/M)）を用いたキャリブレーションは、非スケーリングバージョンと比較して普遍性とトレンドの一貫性が向上した。
• R₁.₄/R₁.₈ をフィッティングパラメータとして追加しても χ² はわずかに低下するが、BIC や AIC では優位性が確認されないため、統計的利点は限定的である。
• 再キャリブレーションパラメータを用いて理論的不確実性を効果的にモデルに組み込むことで、精度を損なわず柔軟性が向上した。
• CoReデータベースからのシミュレーションでは、R₁.₄/R₁.₈ を追加しても完全に解消できないずれが観測され、現在の準普遍的関係式の固有の限界が示唆された。

Figure 2: Exemplary case showing the morphology of NRPMw model. Different wavelet components are reported with different colors: ${W}_{\rm fus}$ in blue, ${\tilde{W}}_{\rm bnc}$ in orange, ${\widetilde{W}}_{\rm pul}$ in green, and ${\tilde{W}}_{\rm peak}$ in purple. The top panel shows the time-doma