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

技术摘要：Fermi-GBM 中与 GW/sGRB 170817A 相似的 MeV-GeV 能段短伽马射线暴

问题陈述
GW/sGRB 170817A 的探测证实了短伽马射线暴（sGRBs）可起源于双中子星（NS-NS）并合。然而，该事件表现出独特的特征：低光度（$L_{iso} \approx 1.6 \times 10^{47}$ erg s$^{-1}$）以及在较大离轴视角（$\sim 20^\circ$）下观测到的独特双分量发射结构（非热峰值后跟热尾）。尽管理论模型表明许多 NS-NS 并合可能产生类似的离轴 sGRB，但由于大多数事件缺乏距离测量，在历史数据集中识别它们极具挑战性。作者旨在从 Fermi-GBM 目录中识别出具有与 sGRB 170817A 特定光谱和时序特征的 sGRB，以便更好地估算引力波（GW）事件的电磁（EM）对应体发生率。

方法论
本研究采用多阶段筛选和聚类方法，应用于 Fermi-GBM 数据集（2008 年 7 月 17 日至 2017 年 8 月 17 日）：

1. 初步筛选：作者选取了 635 个 $T_{90} < 2.05$ 秒的 sGRB。
2. 光谱筛选（A 区）：基于 sGRB 170817A 的第一个峰值，该峰值符合具有特定参数的康普顿化（Comp）模型（$E_{peak} \approx 215$ keV, $\alpha \approx 0.85$），样本被缩减至 80 个事件，这些事件的 $E_{peak}$ 落在 $1\sigma$ 上限（317 keV）范围内。
3. 硬度比（HR）计算：将光变曲线划分为两个时间区：A 区（$T_{start}^{90}$ 至 $T_{end}^{50}$）和 B 区（$T_{end}^{50}$ 至 $T_{end}^{90}$）。计算了两个区域的硬度比（$HR1$和$HR2$），代表 50–300 keV 波段与 10–50 keV 波段的光子计数比率。
4. 异常值剔除：使用马氏距离度量，从 80 个事件的样本中剔除了两个统计异常值（sGRB 081101 和 sGRB 120814），剩余 78 个事件。
5. 聚类：将 K-means 聚类算法应用于 $HR1$与$HR2$ 参数空间。利用 Kneedle 算法（检测惯性曲线中的最大曲率）确定了最佳聚类数（$K=6$）。
6. 最终筛选：进一步分析了包含 sGRB 170817A 及其他 16 个事件的聚类 1。作者应用了两个最终标准：
   • 光变曲线中存在明显的双峰。
   • B 区的光谱拟合符合黑体（BB）模型，且温度 $kT < 12$ keV。
7. 发生率估算：对于选定的 8 个事件（其中 6 个缺乏测量红移），作者利用 Amati 关系（$E_{peak}$ 与 $E_{iso}$）估算了伪红移（$z_p$）。利用这些距离计算了事件发生率，并预测了 LIGO 观测运行（$O1$ 至 $O5$）期间的 GW+sGRB 探测数量。

主要结果

• 类似事件的识别：该研究在 Fermi-GBM 目录中识别出8 个 sGRB，它们在双峰光变曲线结构和光谱演化（非热后跟热发射）方面与 sGRB 170817A 高度相似。识别出的事件为：sGRB 150805、sGRB 150101、sGRB 131128、sGRB 131004、sGRB 130808、sGRB 120524、sGRB 090108 和 sGRB 081122。
• 光谱特征：分析证实，这些事件的第一区（A 区）遵循康普顿化模型，而第二区（B 区）则很好地符合黑体模型，温度范围在 $\sim 7.6$ 至 $12.0$ keV 之间，这与 sGRB 170817A 中观测到的热尾一致。
• 事件发生率预测：利用推导出的光度距离（包括伪红移），作者计算了 LIGO 观测运行期间预期的 GW+sGRB 事件数量：
  • $O1$ 与 $O2$：预期事件数分别约为 0.58 和 1.36。这与 $O2$ 期间观测到 sGRB 170817A 相符。
  • $O3$ 与 $O4$：综合预期约为 6 至 8 个事件。
  • $O5$：外推至即将到来的 $O5$ 运行（预计持续时间 $\sim 884$ 天，距离范围 240–330 Mpc），作者预测 GW+sGRB 事件数量约为 7.49 $\pm$ 3.71。

意义与主张
该论文声称，识别出这 8 个“孪生”sGRB 提供了一个统计上稳健的样本，用于细化产生可观测电磁对应体的 NS-NS 并合发生率。作者认为，这些事件的存在支持了以下假设：许多 NS-NS 并合发生在高倾角，产生可探测但不同于轴上爆发的离轴 sGRB。

该研究强调，预测事件发生率的偏差（例如，如果在未来的运行中探测到的事件显著少于或多于预期）可能会挑战当前关于此类事件随距离演化或喷流几何结构的理解。这项工作作为先前搜索（如 Fermi-GBM 合作组和 Swift 团队的搜索）的统计补充，特别关注从非热分量到热分量的光谱转变，而不仅仅是光变曲线形态。作者得出结论，到 $O4$ 运行结束时，可能会有近 5 个 GW+sGRB 事件，若偏离这一数字，将引发人们对这些事件随距离演化的担忧。