[論文レビュー] Multi-messenger constraints on the neutron-star equation of state and the Hubble constant

Tim Dietrich, M. W. Coughlin+6 more|arXiv (Cornell University)|Feb 26, 2020

Gamma-ray bursts and supernovae被引用数 6

ひとこと要約

本研究では、GW170817、GW190425、AT2017gfo、GRB170817A、およびパルサー観測からのマルチメッセンジャー観測データと、バリオン対称性を用いた有効場理論による核計算を組み合わせ、中性子星の状態方程式を制約し、ハッブル定数を測定した。重力波、キロノーバ光曲線、質量測定を段階的かつベイズ的手法で統合した分析により、1.4 M⊙の中性子星の半径は11.75+0.86−0.81 km、1σ信頼区間でハッブル定数は66.2+4.4−4.2 km Mpc⁻¹ s⁻¹ となった。

ABSTRACT

Observations of neutron-star mergers based on distinct messengers, including gravitational waves and electromagnetic signals, can be used to study the behavior of matter denser than an atomic nucleus, and to measure the expansion rate of the Universe described by the Hubble constant. We perform a joint analysis of the gravitational-wave signal GW170817 with its electromagnetic counterparts AT2017gfo and GRB170817A, and the gravitational-wave signal GW190425, both originating from neutron-star mergers. We combine these with previous measurements of pulsars using X-ray and radio observations, and nuclear-theory computations using chiral effective field theory to constrain the neutron-star equation of state. We find that the radius of a $1.4$ solar mass neutron star is $11.75^{+0.86}_{-0.81}\ m km$ at $90\%$ confidence and the Hubble constant is $66.2^{+4.4}_{-4.2}\ m km \,Mpc^{-1}\, s^{-1}$ at $1\sigma$ uncertainty.

研究の動機と目的

• マルチメッセンジャー観測と核理論を組み合わせることで、中性子星の状態方程式（EOS）の制約を改善すること。
• 二重中性星合体からの重力波および電磁波データを用いてハッブル定数を測定すること。
• 重力波、キロノーバ、ガンマ線バースト、パルサーの質量-半径測定といった天体物理学的観測を統合したベイズ的フレームワークに組み込むこと。
• 波形モデル、分析順序、キロノーバモデルの仮定に対する結果の頑健性を評価すること。
• 低密度領域でバリオン対称性を用いた有効場理論の予測を、高密度領域でモデルに依存しないパラメトリック拡張と統合し、因果律と理論的一致性を保証すること。

提案手法

• 核密度の1.5倍までの密度領域まで、バリオン対称性を用いた有効場理論（EFT）を用いて5,000個の中性星状態方程式（EOS）モデルを構築した。
• 1.5 nsat を超える領域では、因果律を保つために音速平面におけるモデルに依存しないパラメトリック展開を用いてEOSモデルを拡張した。
• 逐次的ベイズ推論を適用し、制約を段階的に統合した：パルサーからの最大中性星質量、GW170817/AT2017gfo、NICERによるPSR J0030+0451の質量-半径測定、IMRPhenomPv2_NRTidalv2を用いたGW170817の重力波波形、AT2017gfoのキロノーバ光曲線、および電磁対応体非検出を伴う/伴わないGW190425のデータ。
• GW170817およびGW190425の重力波形にNRTidalv2波形モデルを適用し、SEOBNRv4Tとの比較により波形モデルの感度をテストした。
• 分析手順の順序を変化させたり、AT2017gfoのキロノーバモデル（I～III）を変更したりすることで、系の不確かさを体系的に伝搬させた。
• 逐次的分析ステージ間で事後分布を伝搬させるためにカーネル密度推定法を用い、一貫性のある事前分布の構築を確保した。

実験結果

リサーチクエスチョン

• RQ1 マルチメッセンジャー観測と核理論を統合した場合、1.4太陽質量の中性星の半径に対する最も精密な制約は何か？
• RQ2 AT2017gfoのキロノーバ光曲線とGW190425の電磁対応体非検出の情報が、中性星状態方程式の制約にどのように影響するか？
• RQ3 異なる重力波形モデル（例：NRTidalv2 と SEOBNRv4T）が、推定される中性星半径およびハッブル定数に与える影響はどの程度か？
• RQ4 ベイズ解析において、天体物理学的および理論的制約を適用する順序が、最終的な制約に与える影響はどの程度か？
• RQ5 核EFTの制約と組み合わせて、GW170817、AT2017gfo、およびパルサー質量データを統合した場合、ハッブル定数はどのように測定されるか？

主な発見

• すべてのマルチメッセンジャー観測データとバリオン対称性EFTの制約を統合した結果、1.4 M⊙の中性星半径は90%信用区間で11.75+0.86−0.81 km に制約された。
• GW170817、AT2017gfo、およびパルサー質量測定の共同解析から、ハッブル定数は1σ信頼区間で66.2+4.4−4.2 km Mpc⁻¹ s⁻¹ に測定された。
• GW170817およびAT2017gfoの情報をGW190425の解析に組み込むことで、質量加重平均の潮汐歪み度 ˜Λ は140.80+144.22−64.73 から117.90+114.60−49.24 に低下し、より小さい半径が支持された。
• 波形モデルの選択に頼らずに結果が安定している：NRTidalv2の代わりにSEOBNRv4Tを用いた場合、半径に及ぼす影響はわずかで、不確かさの範囲内に収まった。
• 分析手順の順序にわずかな影響しか与えず、NICERデータを最初に適用するのではなく最後に適用しても、半径の制約は安定した。
• AT2017gfoのキロノーバモデル（I～III）を変更しても、距離-傾き角の測定値は統計的不確かさの範囲内で一貫しており、結果が統計的要因によって支配されていることが示された。