Constraining Binary Neutron Star Populations using Short Gamma-Ray Burst Observations

技術要約：短時間ガンマ線バーストのデータを用いた連星中性子星個体数合成モデルの制約

問題提起
連星中性子星（BNS）の画期的なマルチメッセンジャー検出であるGW170817と、それに伴う短時間ガンマ線バースト（sGRB）GRB 170817Aは、BNS合体がsGRBを生成し得ることを裏付けた。しかし、現在の重力波（GW）観測所によるBNS検出数は限られており、BNS合体のみで観測されたsGRB個体群の全容を説明できるかどうかという疑問が残されている。これらの個体群を結びつける重要な不確実性は、「ジェット分率」（$f_j$）であり、これはBNS合体のうち、合体残骸の噴出物を突き抜けて相対論的ジェットを放出することに成功した割合と定義される。これまでの$f_j$に関する制約は、ジェットの幾何学構造や固有のBNS合体率に関する仮定に応じて、数十パーセントから100%まで幅広く分散している。本研究は、理論的なBNS個体数合成モデルと観測されたsGRB個体群との間の緊張関係に対処し、非物理的なジェット分率や非現実的なジェット幾何学構造を導入することなく、低レートのBNSモデルがFermi-GBMのsGRBカタログを再現できるかどうかを検証する。

手法
著者らは、64種類の異なるBNS個体数モデルを、16年間にわたるFermi-GBM計器（2008年7月〜2025年5月）のsGRBアーカイブに接続する、包括的かつエンドツーエンドのフレームワークを用いている。解析には、合成sGRBカタログをシミュレートし、各BNSモデルに対して必要なジェット分率（$f_j$）を推論するために実装されたMAGGPYコード内のマルコフ連鎖モンテカルロ（MCMC）法を利用する。

1. BNS個体数モデル： 本研究では、sevn連星進化コードとcosmoRate宇宙合体率を用いて生成された、Iorioら（2023）による16の基本モデルを使用する。これらのモデルは、共通外層（common-envelope）効率（$\alpha_{CE} \in \{0.5, 1, 3, 5\}$）、出生時キック分布、および質量移動の安定性基準を変化させることで、恒星進化の不確実性を探索している。これにより、局所的な合体率（$R_{BNS}(0)$）および赤方偏移進化が異なる64個のユニークなBNS個体群が得られる。
2. sGRB放射モデル： ジェット幾何学構造の制約の堅牢性を評価するため、3つの異なるシナリオをテストする：
   • ユニバーサル構造化ジェット（Universal Structured Jet）： GW170817/GRB 170817Aのアフターグローに基づいて較正された、コア開口角 $\theta_c \approx 3.4^\circ$ を持つ構造化角プロファイル。
   • ユニバーサル・トップハット・ジェット（Universal Top-Hat Jet）： 一様な放射円錐と固定されたコア開口角 $\theta_c$ を持つ簡略化されたモデル。
   • 非ユニバーサル・トップハット・ジェット（Non-Universal Top-Hat Jets）： $\theta_c$ がフラット（一様）または対数正規分布に従って個体群間で変化するモデル。
3. 統計的枠組み： MCMCサンプリングは、自由パラメータ（$f_j$、スペクトル指数、および幾何学的パラメータを含む）を調整し、観測されたFermi-GBMカタログを再現する尤度を最大化する。尤度関数は、イベントレートのためのポアソン項と、合成データと実データの観測量（フルエンス、ピークフラックス、ピークエネルギー、および持続時間）の分布形状を比較するための二標本クラメール・フォン・ミーゼス（CvM）検定を組み合わせている。
4. 生存可能性指標： 推論されたジェット分率の事後中央値が $f_j \leq 1$ を満たす場合、そのBNS個体群は「物理的」であるとみなされる。$f_j > 1$ という値は、たとえすべての合体がジェットを放出したとしても、そのBNS個体群が集団としてのsGRBレートを生成するには不十分であることを示している。

主な結果
解析により、BNS個体数およびジェット物理に関する以下の制約が得られた：

• 低レートモデルの排除： ユニバーサル構造化ジェット（GW170817に較正）を仮定した場合、局所的な合体率 $R_{BNS}(0) \lesssim 50 \text{ Gpc}^{-3} \text{ yr}^{-1}$ を予測するBNS個体数モデルは、sGRBの唯一の起源としては事実上排除される。これらのモデルは、観測されたsGRB数を一致させるために、非物理的なジェット分率（$f_j > 1$）を必要とする。
• 好ましい合体率： 局所的な合体率が $R_{BNS}(0) \approx 100 \text{ Gpc}^{-3} \text{ yr}^{-1}$ 付近のモデルは、物理的に妥当なジェット分率（$f_j \approx 0.7 - 0.8$）をもって、観測されたsGRB個体群を正常に再現する。このレートは現在のGWTC-4の制約と一致している。
• ジェット幾何学の緊張関係： ユニバーサル・トップハット・ジェットの仮定に緩和した場合、低レートのBNS個体群は、ジェットのコア開口角が著しく広い（$\theta_c \gtrsim 15^\circ - 20^\circ$）場合にのみ観測と整合する。これは、sGRBのアフターグローからの観測的推論が示す、中央値のコア角 $\langle \theta_c \rangle \approx 6^\circ$ との間に緊張関係を生じさせる。
• 非ユニバーサリティ： ジェット角の分布を許容する（非ユニバーサリティ）ことで得られる緩和効果は限定的である。対数正規分布の角度を用いることで、広いジェットでテール部分を構成することにより、低レートのBNSモデルに必要な中央値 $f_c$ を下げることができるが、これは個体群の相当部分（$>10\%$）が広いジェット（$\theta_c > 15^\circ$）を持つ必要があることを意味しており、そのような広いジェットの観測が極めて稀であることと矛盾する。
• 遅延時間： 解析は、好ましいBNS個体群（$R_{BNS}(0) \approx 100 \text{ Gpc}^{-3} \text{ yr}^{-1}$ のもの）が短い遅延時間（$\langle \tau_d \rangle < 1 \text{ Gyr}$）を示すことを確認しており、これはフラックス制限サンプルにおける選択効果を補正した最近の知見と一致している。

意義および結論
本論文は、電磁波観測（sGRB）のみから導出された、局所的なBNS合体率の堅牢な下限値を確立している。低レートのBNSモデルが、非物理的なジェット分率またはアフターグローのデータと矛盾するジェット幾何学を必要とすることを証明することで、本研究は連星進化理論に対する重要なマルチメッセンジャー制約を提供している。

結果は、標準的または高い共通外層効率（$\alpha_{CE} \geq 1$）と適度な出生時キックを特徴とするBNS個体数合成モデルを支持しており、これらは自然に高い合体率をもたらす。逆に、低い放出効率や高い出生時キックを持つモデルは支持されない。また、高いジェット生成効率（$f_j \approx 0.8$）の推論は、ジェット放出の条件が大多数のBNS合体において満たされていることを示唆しており、ブラックホールへの即時崩壊がジェット形成を系統的に抑制するというシナリオを否定している。

最終的に、本研究は、GWの制約（現在はBNSレートの上限を設定）とsGRBの個体数統計（下限を設定）を組み合わせることで、連星中性子星の進化とジェット物理のための許容パラメータ空間を段階的に縮小できることを強調している。