Insights on the Gamma-Ray Bursts variability in their cosmological rest frame

技术摘要：伽马射线暴在宇宙学静止参考系中的变率洞察

问题陈述
伽马射线暴（GRBs）表现出极端的时变性，其范围从单脉冲到叠加脉冲的复杂序列。虽然火球模型将这种变率归因于中心引擎活动或磁重联，但 GRB 喷流的具体成分（坡印廷通量还是重子物质）仍不明确，因为这些情景预测了不同的变率特征。一个关键的开放性问题是确定最小变率时间尺度（MVT），其定义为观测到 GRB 轮廓发生显著变化的最短时间尺度。确定 MVT 对于约束中心引擎的维度和活动，以及区分不同理论模型至关重要。以往的研究已对低至数百微秒的 MVT 进行了调查，但需要对数千个处于宇宙学静止参考系中的 GRB 进行全面分析，特别是将 MVT 与光谱参数（如各向同性能量 $E_{iso}$ 和峰值能量 $E_p$）进行关联，以完善这些约束。

方法论
作者分析了 2008 年 7 月至 2021 年 10 月期间由 Fermi 伽马射线暴监测器（GBM）探测到的 3,155 个 GRB。在经过数据完整性过滤后，保留了 3,062 个 GRB（511 个短暴，2,551 个长暴）。选取了 150 个具有测量红移的 GRB 子集（129 个长暴，21 个短暴）用于在宇宙学静止参考系中的分析。

数据处理流程包括：

1. 数据检索： 获取触发探测器的非分箱时间触发事件（TTE）数据文件。
2. 背景减除： 使用 gbm.background 程序拟合并减去爆发前 20–5 秒及爆发后 75–150 秒间隔内的背景计数。
3. MVT 计算： 采用了两种基于 Haar 离散小波变换的方法来确定 MVT：
   • 方法 1： 将 MVT 识别为结构函数（$\sigma_{X,\Delta t}$）偏离线性上升阶段（通过 $\chi^2$ 测试确定）的时间尺度。
   • 方法 2： 将 MVT 定义为上升阶段的线性拟合与随后平坦区域拟合之间的交点（基于 Golkhou 等人 2015 年的研究）。
4. 静止系修正： 对于红移样本，通过除以 $(1+z)$ 将观测到的 MVT（$\Delta t_{min}$）转换为宇宙学静止系 MVT（$\Delta t_{min,z}$）。
5. 相关性分析： 研究了不同能量波段（软能段：20–100 keV；硬能段：100–1000 keV）下 $\Delta t_{min,z}$ 与光谱参数（$E_{iso}$、$E_p$）之间的相关性。利用 Amati 关系同时比较了 $\Delta t_{min,z}$、$E_{iso}$ 和 $E_p$。

主要结果

• MVT 数值： 在观测系中，长暴的 MVT 中位数为 738 ms（方法 1）和 583 ms（方法 2）；短暴的中位数为 77.6 ms 和 53 ms，分别为后两者。在宇宙学静止系中，这些数值分别下降至 236 ms（长暴）和 37.3 ms（短暴）。
• 能量依赖性： 在静止系中证实了“窄度-硬度”关系：在硬能段（100–1000 keV）计算的 MVT 一致短于软能段（20–100 keV），表明时间轮廓在高能下变得更加尖锐。
• 时间膨胀： 通过按红移组对长暴进行分箱显示，其趋势与宇宙学时间膨胀（$\Delta t_{min} \propto 1+z$）相符，尽管受限于 2$\sigma$ 置信区间，统计显著性有限。
• 与各向同性能量的反相关性： 研究发现长暴的 $\Delta t_{min,z}$ 与 $E_{iso}$ 之间存在显著的反相关关系（$r = -0.59$），表明具有更高各向同性能量的爆发表现出更快的变率。即使在控制通量的情况下，该关系依然成立。值得注意的是，一些在观测系持续时间（$T_{90}$）上被归类为“短暴”，但具有“长暴”光谱特征的高红移 GRB，在这一关系中表现得像长暴，这对传统的基于持续时间的分类提出了挑战。
• 峰值能量关系： $\Delta t_{min,z}$ 被发现与峰值能量（$E_{p,z}$）呈反相关。对于短暴，关系遵循 $\Delta t_{min,z} \propto E_{p,z}^{-0.98}$；对于长暴，关系遵循 $\Delta t_{min,z} \propto E_{p,z}^{-0.62}$。
• 总变率： 研究计算了总变率参数（$V$），发现 $V$ 与 MVT 之间没有明显的关联，尽管观察到 $V$ 与 $E_{iso}$ 之间存在轻微的反相关。
• 模拟： 一个包含微秒级脉冲的模拟 GRB 表明，只要信噪比足够，分析工具能够识别低至 3–5 $\mu$s 的变率。

意义与主张
本文声称利用大规模样本对宇宙学静止系中的 GRB 变率进行了稳健的统计表征。其主要贡献在于：

1. 确认物理趋势： 该研究通过比以往工作更大的样本量，验证了 GRB 脉冲的能量依赖性尖锐化现象以及变率时间尺度的时间膨胀效应。
2. 中心引擎约束： 观测到的 MVT 与 $E_{iso}$（以及 $E_p$）之间的反相关性表明，爆发的总能量输出与中心引擎活动的活动时间尺度之间存在物理联系。作者指出，如果内激波模型成立，则 MVT 追踪了中心引擎与 $\gamma$ 射线产生位点之间的距离。
3. 分类挑战： 结果强调了仅基于 $T_{90}$ 的传统长/短 GRB 分类的局限性。具有短持续时间但具有长暴光谱特性（及 MVT-$E_{iso}$ 关系）的高红移爆发表明，替代性的分类标准（如 Amati 关系或光谱硬度）可能具有更强的物理动机。
4. 未来准备： 该工作证明了当前的分析技术有能力解析微秒级变率，为具有改进计时分辨率的未来任务（如 HERMES Pathfinder）奠定了基础，以探测亚毫秒量级，并可能区分磁重联与内激波模型。

作者对这些发现的决定性性质保持谦逊，指出目前针对短暴的统计效力有限，且 MVT-$E_{iso}$ 关系的确切物理起源仍需进一步的理论研究。他们强调，微秒级变率的存在仍需下一代仪器予以证实。