[論文レビュー] High-resolution numerical relativity simulations of spinning binary neutron star mergers

Tim Dietrich, Sebastiano Bernuzzi+2 more|arXiv (Cornell University)|Mar 21, 2018

Pulsars and Gravitational Waves Research被引用数 6

ひとこと要約

本論文は、高分解能数値相対論的シミュレーションを用いて回転する二重中性子星合体の最先端の結果を提示しており、高次の収束法と複数のグリッド解像度を用いて、約12回転の間に位相誤差0.5–1.5ラジアンの波形を達成している。これらのシミュレーションにより、重力波テンプレートの正確なモデリングと準解析的波形モデルの検証が可能となり、多メッセンジャー天文学分野における重要な前進をもたらしている。

ABSTRACT

The recent detection of gravitational waves and electromagnetic counterparts emitted during and after the collision of two neutron stars marks a breakthrough in the field of multi-messenger astronomy. Numerical relativity simulations are the only tool to describe the binary's merger dynamics in the regime when speeds are largest and gravity is strongest. In this work we report state-of-the-art binary neutron star simulations for irrotational (non-spinning) and spinning configurations. The main use of these simulations is to model the gravitational-wave signal. Key numerical requirements are the understanding of the convergence properties of the numerical data and a detailed error budget. The simulations have been performed on different HPC clusters, they use multiple grid resolutions, and are based on eccentricity reduced quasi-circular initial data. We obtain convergent waveforms with phase errors of 0.5-1.5 rad accumulated over approximately 12 orbits to merger. The waveforms have been used for the construction of a phenomenological waveform model which has been applied for the analysis of the recent binary neutron star detection. Additionally, we show that the data can also be used to test other state-of-the-art semi-analytical waveform models.

研究の動機と目的

• 回転する中性子星を伴う二重中性星合体からの高精度な重力波波形の生成を目的とし、数値相対論分野における重要な空白を埋めること。
• 収束テストとリチャードソン補外を厳密に適用することで数値誤差を低減し、約12回転の間に位相不確実性を約1ラジアンにまで抑えること。
• 高精度なシミュレーションデータを用いて、現象論的および効率的1体（EOB）波形モデルの開発とキャリブレーションを可能にすること。
• 波形モデリングにおける仮定（例：スピン・潮汐効果とスピン・スピン効果の分離可能性）の妥当性を、高分解能データの分析を通じて検証すること。
• 将来の重力波探知と多メッセンジャー天文学、特にGW170817に類似するイベントの解析を目的としたベンチマークデータセットの提供。

提案手法

• 複数のHPCクラスタ上で高次収束有限差分法を用い、適応メッシュ細分化（AMR）と複数のグリッド解像度を組み合わせることで、数値収束を保証する。
• 不自然な位相振動を低減し、波形品質を向上させるために、非円形度が0.001未塔の準円形初期データを用いる。
• 自己収束テストにより重力波位相における2次収束が確認され、信頼性の高い誤差推定と連続波形のためのリチャードソン補外が可能になる。
• 波形はNRTidal現象論的波形モデルの構築と、時間領域および周波数領域での比較による既存EOBモデル（SEOBNRv4T、TEOBDResumS）の検証に用いられる。
• Padé近似を周波数領域での潮汐位相補正のモデル化に適用し、NRTidalモデルにおける潮汐効果の閉形式表現を可能にする。
• 有限温度状態方程式（EOS）を含む現実的な物理を組み込み、SuperMUC、JURECA、Marconiを含むシステム上でBAMコードを用いて実行された。

実験結果

リサーチクエスチョン

• RQ1 回転する二重中性星合体の高分解能シミュレーションにおいて、重力波位相の達成可能な精度はどの程度か？
• RQ2 スピンと潮汐効果は波形においてどのように結合するか？また、波形モデリングにおいてそれらを独立に取り扱うことは可能か？
• RQ3 現在の準解析的波形モデル（例：NRTidal、SEOBNRv4T、TEOBDResumS）は、高精度な数値相対論的波形とどの程度一致するか？
• RQ4 収束テストに基づくリチャードソン補外は、複雑で非線形な系からの波形における数値不確実性を信頼性を持って低減できるか？
• RQ5 残留非円形度の波形品質への影響は何か？また、初期データで効果的に低減できるか？

主な発見

• シミュレーションは重力波位相において2次収束を達成しており、信頼性の高い誤差推定とリチャードソン補外による連続波形の作成が可能となり、約12回転の間に位相不確実性を約1ラジアンにまで低減した。
• 位相誤差は約12回転の間に0.5–1.5ラジアンの範囲に評価され、これまでで最も高精度な回転する二重中性星合体のシミュレーションである。
• シミュレーションデータから構築されたNRTidalモデルは、数値波形と~0.1ラジアン以内で一致し、GWデータ解析におけるその有用性が検証された。
• SEOBNRv4TとTEOBDResumSモデルは、最後の数回転においてそれぞれ約0.1ラジアンおよび約0.4ラジアンのずれを示し、改良の余地があることが示された。
• 分析により、スピン・潮汐結合およびスピン・スピン相互作用は現在のシミュレーションでは無視できるほど小さく、波形モデリングにおける効果の分離仮定が支持された。
• シミュレーションには約2500万CPU時間が必要であり、複数のHPCシステムで実行されたことから、現代のスーパーコンピュータ上で大規模かつ高精度な数値相対論が実現可能であることが示された。