CTAO/LST-1 探测到的 GRB 221009A 延迟信号能用对回波解释吗？

对回波是延迟的光，不是第二次爆发

从遥远源发出的甚高能伽马射线并不会毫无干扰地穿过星系际空间。它们在途中与微波和红外/光学背景光子相互作用，转化为电子-正电子对。这些电子对随后将宇宙微波背景光子散射到较低的伽马射线能量，重新生成一束次级辐射。如果大尺度结构空洞中的星系际磁场（IGMF）使这些带电粒子对发生偏转，次级伽马射线就会以略微不同的方向、比主信号更晚到达。这一延迟且偏转的成分就是对回波。

延迟正是这一效应可被测量的原因。论文给出了延迟随伽马射线能量、场强 B 和关联长度 lambda_B 变化的表达式，分为两种情形：当 lambda_B 超过电子冷却距离时，延迟与 B 的平方成正比；当 lambda_B 更短时，延迟还随 lambda_B 变化。冷却距离本身依赖于能量，在 LST-1 约 0.3 TeV 的阈值能量处，它约等于 20 kpc。正是这个数字使 20 kpc 成为整个分析中的分界点。

由于回波由能量最高的主光子产生，300 GeV 以上观测到的对回波来自 10 TeV 以上的主伽马射线。这正是 LHAASO 对 GRB 221009A 测量的能段，也是本文能够进行联合拟合的原因。

瞬变爆发把回波与余辉分开

活动星系核一直是搜寻对回波的传统目标，但它们的辐射是持续性的，很难把自身信号与延迟回波区分开。伽马射线暴提供了更干净的几何条件：TeV 信号是瞬变的，因此原则上可以把延迟探测与一个明确的瞬时辐射窗口进行对比。

GRB 221009A 是迄今观测到的最亮伽马暴，LHAASO 探测到其能量超过 10 TeV。LHAASO 在爆发开始后约 3000 s 内记录到该信号，最亮辐射出现在触发后 230 s 至 900 s 之间。该暴的红移 z 约为 0.151。这些经过星系际背景光吸收修正的测量谱，就是作者用来归一化其对回波预测的输入。

CTAO 的 LST-1 望远镜在爆发后 1.3 至 20 天观测了该源。首次观测期间强烈的月光使能量阈值从名义上的约 20 GeV 提高到约 300 GeV。在 1.33 天时，观测在源方向发现了 4.1 sigma 的超出；此后的观测只给出了流量上限。

在元胞状磁场中模拟粒子对级联

作者用截断幂律谱模拟主辐射，对两个时间段内的 LHAASO 数据拟合归一化、谱指数和截断能量。随后他们使用数值程序 CRProp 3.2，将约 10,000 个能量在 1-100 TeV 范围内的主光子注入星系际空间传播，同时追踪伽马射线以及它们产生的电子和正电子。

星系际磁场被处理为元胞状：每个元胞内场强恒定，但元胞之间方向随机。所探索的参数空间覆盖 1e-18 G 至 1e-14 G 的场强，以及 1e-4 Mpc 至 1 Mpc 的元胞尺度。GRB 221009A 被假定为一个点源，在张角半宽 1 度的喷流内辐射，这一角度由 LHAASO 光变曲线中看到的喷流拐折推断而来。

对每一组 B 和 lambda_B 的组合，模型都会预测一条对回波光变曲线和能谱。随后作者对 LHAASO 和 LST-1 的谱数据点进行联合拟合，假定每个能量区间服从高斯统计，并针对固定的关联长度导出作为 B 函数的剖面似然。

拟合倾向于约 1e-16 G 的磁场

剖面似然显示出清晰的最小值，作者将其解读为 LHAASO 和 LST-1 数据给出了 IGMF 参数的测量，而不仅仅是一个下限。B 的 68% 置信范围定义为剖面似然 2 log(L) 相对其最大值下降不超过 1 的区间。

对于大于 20 kpc 的关联长度，得到的估计为 log(B/G) = -16.2 +/- 0.7，对应约 3e-16 G 的磁场，这是作者在讨论中所用的量级意义上的数值。对于小于 20 kpc 的关联长度，该估计要向上修正一个因子，即 (20 kpc / lambda_B) 的平方根。20 kpc 处的转折与那些产生约 0.3 TeV 逆康普顿辐射的粒子对的冷却距离有关，而 0.3 TeV 接近 LST-1 的阈值。

在宇宙学相变遗留的短关联长度磁场这一特定情形下，对 LST-1 信号的对回波解释意味着场强在 1e-14 G 至 1e-12 G 之间，且关联长度预计短于 1 秒差距。作者指出，这些参数与已知的 IGMF 限制一致，包括来自活动星系核伽马射线观测的下限。

源自身余辉也能产生同样信号

论文明确陈述了核心保留意见。把 LHAASO 的 TeV 光变曲线按幂律外推到首次 LST-1 观测时刻，所预测的流量在不确定度范围内与 LST-1 测量一致；把 Fermi/LAT 能谱外推到 0.1-1 TeV 范围也得到类似结论。因此，源自身余辉对 LST-1 信号有可观贡献是可能的。

论文还指出，原则上可以利用信号光变曲线特有的能量依赖关系把对回波与源自身余辉区分开。但对 GRB 221009A 而言这一检验无法进行，因为强烈月光阻碍了爆发后立即进行的切伦科夫观测。作者直接说明，LST-1 的探测不能被视为对回波的证据，只能说对回波可以对所测流量提供量级为一的贡献。

这就是为什么该结果被表述为在一个假设下的测量，而不是一次探测。所引用的场强以对回波解释是那个 1.33 天信号的正确解释为前提。

单个伽马暴无法确立这一现象

作者明确指出，来自单个伽马暴的日尺度延迟信号本身并不构成对回波现象的发现。能加强这一论断的是在多个伽马暴和活动星系核中重复探测到类似信号，且都与同一组 IGMF 参数一致。这样的组合才能提供更有说服力的证据和更可靠的磁场测量。

论文把当前结果定位为对正在建造的下一代成像大气切伦科夫望远镜 CTAO 和 LACT 的方法潜力的展示。作者还提到一篇同期发表的论文，该文从同一 LST-1 信号得出了类似的 IGMF 估计；他们指出自己的分析不同之处在于处理了所估计磁场的不确定度，而不只是给出一个量级数值。

因此对读者而言，关键比较不是 LST-1 信号是否真实，而是它的哪种解释能在未来数据中存活。论文自己设定的门槛是多个源在共享同一组 IGMF 参数下的一致性。

估计值取决于假定的关联长度

场强并不是一个独立于模型的单一数字。只有在关联长度大于 20 kpc 时，它才被报告为 log(B/G) = -16.2 +/- 0.7；低于该尺度时，中心值要向上移动 (20 kpc / lambda_B) 的平方根倍。读者在把这一结果与其他 IGMF 限制比较时，需要确认那些限制假定的是哪个关联长度区间，因为同一批数据在两种情形下意味着不同的场强。

遗留磁场情形是又一个条件分支。1e-14 G 至 1e-12 G 的范围专门适用于宇宙学起源的短关联长度磁场，其允许参数是通过把对回波估计与描述最大湍流处理涡旋尺度的理论带相交得到的。产生较高场强数值的是这一交叠，而不是单独的对回波拟合。

因此本文设定的实际优先事项是：把 LST-1 信号当作一种方法的检验案例，并在读取所引用的场强时同时考虑其假定的关联长度，以及它与源自身余辉之间尚未解决的简并。决定性的下一步是进行被月光阻断的能量依赖测量，并在更多伽马暴中重复。

本文引用的文献