Inferring Neutron Star Nuclear Properties from Gravitational-Wave and Gamma-Ray Burst Observations

技術的要約：重力波およびガンマ線バースト観測からの中性子星核の性質の推論

問題提起
キロノバ放射（例：GRB 211211A、GRB 230307A）を伴う長周期ガンマ線バースト（lGRB）の最近の観測は、lGRBが質量星の崩壊のみから生じるという従来のパラダイムに疑問を投げかけている。もしこれらのイベントが連星中性子星（BNS）または中性子星–ブラックホール（NSBH）の合体から生じているのであれば、決定的な問いが生じる。すなわち、バーストの持続時間を支配する物理的特性は何であるか？理論モデルは、合体の総質量が合体残骸の寿命を決定することを示唆している。より軽い合体は、短周期ガンマ線バースト（sGRB）を駆動する長寿命の中性子星（NS）残骸を形成する可能性がある一方で、より重い合体は短寿命の残骸を形成してブラックホール（BH）へと崩壊し、大規模な降着円盤を介してlGRBを駆動する可能性がある。しかし、これらの結果を分ける具体的な質量閾値、およびそれらと基本的な核特性、特にトルマン・オッペンハイマー・ヴォルコフ質量（$M_{\text{TOV}}$）との関係は、観測的に未制約のままである。

手法
著者らは、LIGO-Virgo-KAGRA（LVK）コラボレーションによる重力波（GW）合体率を、観測されたsGRBおよびキロノバを伴うlGRBの発生率と比較することにより、特徴的な質量スケールを推論するためのベイズ・フレームワークを開発した。

1. 特徴的な質量の定義： 本研究では、$M_{\text{TOV}}$に対する3つの決定的な総質量閾値をパラメータ化する：

   • $M_{\text{VL}} = k_1 M_{\text{TOV}}$：極めて長寿命な残骸と長寿命な残骸の境界。
   • $M_{\text{LS}} = k_2 M_{\text{TOV}}$：長寿命な残骸と短寿命な残骸の境界。$M_{\text{VL}} < M_{\text{tot}} \le M_{\text{LS}}$の合体は、sGRBを生成すると仮定される。
   • $M_{\text{PC}} = k_3 M_{\text{TOV}}$：短寿命な残骸と即時崩壊の境界。$M_{\text{LS}} < M_{\text{tot}} \le M_{\text{PC}}$の合体は、lGRBを生成すると仮定される。$M_{\text{tot}} > M_{\text{PC}}$の場合、質量比（$q$）が十分な円盤形成を可能にする場合（BNSでは $1.2 \le q \le 3$、NSBHでは $q \le 3$）にのみlGRBが生成される。

2. データソース：

   • GW発生率： LVKのGWTC-4カタログ（Abac et al. 2025a）およびFullPop-4.0集団モデル（Abac et al. 2025b）は、成分質量ごとの局所的な合体率を提供する。
   • GRB発生率： sGRBの固有発生率（$R_{\text{sGRB}}$）および発生率比 $\eta = R_{\text{lGRB}}/R_{\text{sGRB}}$ は、大きな不確実性を伴う観測的制約として扱われ、正規分布を用いてモデル化される。

3. 推論： 著者らは、ベイズの定理を用いて、パラメータ集合 $\vec{\Theta} = (k_1, k_2, k_3, M_{\text{TOV}})$ の事後分布を推論する。尤度関数は、質量閾値から導出されたGW誘起GRB率と、観測されたGRB率を比較する。$M_{\text{TOV}}$および$M_{\text{TOV}}$における中性子星半径（$R_{\text{TOV}}$）の事前分布は、パルサーおよびGWの観測に基づき、Legred et al. (2021) から抽出される。

主な結果

• 特徴的な質量比（$k_2$）： 解析の結果、中央値として $k_2 = 1.36^{+0.08}_{-0.09}$（68%信頼区間）が得られ、$k_2 > 1.24$ が90%の信頼度で示された。これは、長寿命から短寿命の残骸への遷移が、およそ $1.36 M_{\text{TOV}}$ で起こることを示している。
• 残骸の寿命への示唆： $k_2$ の高い値は、質量が大きな中性子星が、合体直後に崩壊するのではなく、数百ミリ秒という長い期間生存できることを示唆している。
• 相関関係： 本研究は、$k_2$ と $k_3$ の間の強い相関、および $k_2$ と $M_{\text{TOV}}$ の間の負の相関を特定している。この負の相関により、$k_2$ に対する制約を $M_{\text{TOV}}$ の上限に直接マッピングすることが可能となる。
• ロバスト性： $k_2$ の推論は、想定されるsGRB発生率（$30$から$2000 \text{ Gpc}^{-3}\text{yr}^{-1}$）および発生率比（$\eta$が$0.01$から$3$）の広い範囲においてロバストである。質量比の閾値（$q_{\text{low}}$および$q_{\text{high}}$）の変化や、$M_{\text{VL}}$境界の包含・除外による影響は、推論された$k_2$ に最小限の影響しか与えない。

意義および主張
本論文は、GW合体集団をGRBの持続時間と結びつけることにより、中性子星の状態方程式（EoS）を制約するための、新しい独立した手法を確立している。sGRBとlGRBのプロジェニタを分ける特徴的な質量スケール（$1.36 M_{\text{TOV}}$）を特定することで、著者らは中性子星の最大質量に対する直接的な観測的制約を提供している。

著者らは、本研究の結果が、$1.3$から$1.4 M_{\text{TOV}}$ の間の遷移を予測する数値相対論シミュレーション（例：Perna et al. 2025）と一致していると述べている。また、現在の結果は観測の不確実性に対してロバストであるが、将来的なGWおよびGRBの発生率の測定の改善、特にGRBとGWイベントの同時検出が進むことで、特徴的な質量と $M_{\text{TOV}}$ の間の相関がより強固になり、潜在的に $M_{\text{TOV}}$ の不確実性を $\sim 0.1 M_{\odot}$ まで低減できる可能性があると強調している。このアプローチは既存の手法を補完するものであり、中性子星のEoSに対してより厳格な制限を課すための経路を提供するものである。