Similar Fermi-GBM sGRBs to GW/sGRB 170817A in MeV-GeV energies

技術的概要：MeV-GeV エネルギー領域における GW/短ガンマ線バースト 170817A と類似するフェルミ-GBM 短ガンマ線バースト

問題提起
GW/短ガンマ線バースト 170817A の検出は、短ガンマ線バースト（sGRB）が連星中性子星（NS-NS）合体に由来し得ることを確認させた。しかし、この事象は独特の特性を示した：低い光度（$L_{iso} \approx 1.6 \times 10^{47}$ erg s$^{-1}$）と、大きなオフ軸観測角度（$\sim 20^\circ$）で観測された、非熱的ピークに続いて熱的テールが続く明確な二重成分放出構造である。理論モデルは、多くの NS-NS 合体が同様のオフ軸 sGRB を生成する可能性を示唆しているが、大半の事象において距離測定が欠如しているため、歴史的データセットからそれらを特定することは困難である。著者らは、重力波（GW）事象に対する電磁波（EM）対応天体の発生率をより正確に見積もるために、フェルミ-GBM カタログ内の sGRB のうち、sGRB 170817A と特定のスペクトルおよび時間的特徴を共有するものを特定することを目的としている。

手法
本研究は、フェルミ-GBM データセット（2008 年 7 月 17 日 – 2017 年 8 月 17 日）に適用される多段階のフィルタリングおよびクラスタリングアプローチを採用している：

1. 初期選択： 著者らは $T_{90} < 2.05$ 秒の 635 個の sGRB を選択した。
2. スペクトルフィルタリング（ゾーン A）： sGRB 170817A の最初のピークは、特定のパラメータ（$E_{peak} \approx 215$ keV, $\alpha \approx 0.85$）を持つコンプトン化（Comp）モデルに適合することに基づき、サンプルを $E_{peak}$ の $1\sigma$ 上限値（317 keV）以内の 80 個の事象に絞り込んだ。
3. 硬さ比（HR）の計算： 光曲線を 2 つの時間ゾーンに分割した：ゾーン A（$T_{start}^{90}$ から $T_{end}^{50}$）およびゾーン B（$T_{end}^{50}$ から $T_{end}^{90}$）。両ゾーンについて硬さ比（$HR1$および$HR2$）を計算し、10–50 keV 帯に対する 50–300 keV 帯の光子数の比率を表した。
4. 外れ値の除去： マハラノビス距離指標を用いて、80 個の事象サンプルから 2 つの統計的外れ値（sGRB 081101 および sGRB 120814）を除去し、78 個の事象を残した。
5. クラスタリング： $HR1$対$HR2$ のパラメータ空間に対して K-means クラスタリングアルゴリズムを適用した。最適なクラスタ数（$K=6$）は、慣性曲線における最大曲率を検出する Kneedle アルゴリズムを用いて決定された。
6. 最終選択： sGRB 170817A と他の 16 個の事象を含むクラスター 1 をさらに分析した。著者らは 2 つの最終基準を適用した：
   • 光曲線に明確な二重ピークが存在すること。
   • ゾーン B のスペクトル適合が、温度 $kT < 12$ keV の黒体（BB）モデルであること。
7. 発生率の推定： 8 つの選択された事象（そのうち 6 つは赤方偏移が測定されていなかった）について、著者らはアマトイ関係（$E_{peak}$ 対 $E_{iso}$）を用いて疑似赤方偏移（$z_p$）を推定した。これらの距離を用いて事象発生率を計算し、LIGO 観測ラン（O1 から O5）における GW+sGRB 検出数を予測した。

主要な結果

• 類似事象の特定： 本研究は、フェルミ-GBM カタログにおいて、二重ピーク光曲線構造およびスペクトル進化（非熱的放出に続いて熱的放出）の点で sGRB 170817A と酷似する 8 個の sGRB を特定した。特定された事象は以下の通りである：sGRB 150805, sGRB 150101, sGRB 131128, sGRB 131004, sGRB 130808, sGRB 120524, sGRB 090108, および sGRB 081122。
• スペクトル特性： 分析により、これらの事象の最初のゾーン（ゾーン A）はコンプトン化モデルに従い、2 番目のゾーン（ゾーン B）は sGRB 170817A で観測された熱的テールと一致する、$\sim 7.6$ から $12.0$ keV の範囲の温度を持つ黒体モデルでよく適合することが確認された。
• 発生率予測： 導出された光度距離（疑似赤方偏移を含む）を用いて、著者らは LIGO 観測ランにおける GW+sGRB 事象の期待数を計算した：
  • O1 および O2： それぞれ約 0.58 および 1.36 個の事象が期待される。これは O2 期間中の sGRB 170817A の観測と一致する。
  • O3 および O4： 合計期待値は約 6 から 8 個の事象である。
  • O5： 今後の O5 ラン（推定期間 $\sim 884$ 日、距離範囲 240–330 Mpc）への外挿により、著者らは約 7.49 $\pm$ 3.71 個の GW+sGRB 事象を予測している。

意義と主張
本論文は、これら 8 つの「双子」sGRB を特定することが、観測可能な EM 対応天体を生成する NS-NS 合体の発生率を精緻化するための統計的に堅牢なサンプルを提供すると主張している。著者らは、これらの事象の存在が、多くの NS-NS 合体が高傾斜で発生し、検出可能ではあるが軸上バーストとは区別されるオフ軸 sGRB を生成するという仮説を支持すると論じている。

本研究は、単なる光曲線形態だけでなく、非熱的成分から熱的成分へのスペクトル遷移に特化して焦点を当てることで、フェルミ-GBM 協力グループや Swift チームによる以前の探索に対する統計的な補完であることを強調している。著者らは、O4 ラン終了時点では GW+sGRB 事象が約 5 件存在する可能性があり、この数からの逸脱は距離に伴うこれらの事象の進化に関する懸念を提起すると結論付けている。