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$ 的约束差异巨大，从百分之几十到 100% 不等，这取决于对喷流几何结构和内在 BNS 并合率的假设。本研究旨在解决理论 BNS 群体合成模型与观测到的 sGRB 群体之间的张力，具体测试低速率 BNS 模型是否能在不引入非物理喷流比例或不现实喷流几何结构的情况下，重现 Fermi-GBM 的 sGRB 目录。

方法论
作者采用了一个端到端的综合框架，将 64 个不同的 BNS 群体合成模型与来自 Fermi-GBM 仪器（2008 年 7 月至 2025 年 5 月）的 16 年 sGRB 档案联系起来。分析利用了在 MAGGPY 代码中实现的马尔可夫链蒙特卡洛（MCMC）方法来模拟合成的 sGRB 目录，并推断所需的喷流比例（$f_j$）。

1. BNS 群体模型： 研究使用了来自 Iorio 等人（2023）的 16 个基础模型，这些模型是使用 sevn 双星演化代码和 cosmoRate 生成的宇宙并合率。这些模型通过改变共用包层效率（$\alpha_{CE} \in \{0.5, 1, 3, 5\}$）、诞生踢度分布（natal kick distributions）和质量转移稳定性判据，探索了恒星演化的不确定性。由此产生了 64 个独特的 BNS 群体，具有变化的局部并合率（$R_{BNS}(0)$）和红移演化。
2. sGRB 发射模型： 测试了三种不同的喷流几何情景，以评估约束的稳健性：
   • 通用结构化喷流（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 目录的似然度。似然函数结合了一个用于事件率的泊松项，以及一个两样本 Cramér-von Mises (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 的唯一前身。这些模型需要非物理的喷流比例（$f_j > 1$）才能匹配观测到的 sGRB 计数。
• 首选并合率： 局部并合率在 $R_{BNS}(0) \approx 100 \text{ Gpc}^{-3} \text{ yr}^{-1}$ 左右的模型能够以物理上合理的喷流比例（$f_j \approx 0.7 - 0.8$）成功重现观测到的 sGRB 群体。这一速率与当前的 GWTC-4 约束一致。
• 喷流几何张力： 当放宽结构化喷流假设转而采用通用顶帽模型时，只有当喷流核心张角显著变宽（$\theta_c \gtrsim 15^\circ - 20^\circ$）时，低速率 BNS 群体才能与观测结果相协调。这与来自 sGRB 余辉的观测推论产生了张力，后者表明中值核心角 $\langle \theta_c \rangle \approx 6^\circ$。
• 非普遍性： 放宽普遍性假设（允许喷流角度分布）提供的缓解作用有限。虽然对数正态分布的角度可以通过增加宽喷流的比例来降低低速率 BNS 模型所需的其中值 $f_j$，但这要求很大一部分（$>10\%$）的群体拥有宽喷流（$\theta_c > 15^\circ$），这与宽喷流观测稀缺的情况仍然存在张力。
• 延迟时间： 分析确认，受青睐的 BNS 群体（即 $R_{BNS}(0) \approx 100 \text{ Gpc}^{-3} \text{ yr}^{-1}$ 的模型）表现出短延迟时间（$\langle \tau_d \rangle < 1 \text{ Gyr}$），这与近期针对通量限制样本进行选择效应修正后的发现一致。

意义与结论
本文建立了一个纯粹基于电磁学观测的 BNS 局部并合率的稳健下限。通过证明低速率 BNS 模型要么需要非物理的喷流比例，要么需要与余辉数据不一致的喷流几何结构，本研究为双星演化理论提供了关键的多信使约束。

结果倾向于具有标准或高共用包层效率（$\alpha_{CE} \geq 1$）和适中诞生踢度的 BNS 群体合成模型，这些模型自然会产生较高的并合率。相反，具有低抛射效率或高诞生踢度的模型则受到排斥。高喷流产生效率（$f_j \approx 0.8$）的推断表明，在大多数 BNS 并合中，发射相对论性喷流的条件得到了满足，这不支持由于向黑洞快速坍缩而系统性抑制喷流形成的场景。

最终，这项工作强调了结合引力波约束（目前设定了 BNS 率的上界）与 sGRB 群体统计（设定了下界）的力量，从而逐步缩小双中子星演化和喷流物理学的允许参数空间。