光变曲线拿到手只是第一步,真正硬核的是搞清楚它背后的物理。研究团队将AT 2020afjz与Kann样本(D. A. Kann et al. 2010汇编的31个具有实测红移且解析出上升阶段的GRB余辉样本)进行了对比。结果发现,AT 2020afjz的形态与样本中的大多数GRB余辉都不一样——它上升得太慢了。
图3:AT 2020afjz(黑色)与GRB 260310A(蓝色方块;Y.-H. Yang et al. 2026)以及D. A. Kann et al. (2010)样本中具有实测上升期的GRB余辉(灰色)的峰值归一化光学光变曲线对比;其中GRB 080710(金色)的静止系上升时标与本文事件最接近。AT 2020afjz的上升形态和峰值时间更接近被解释为离轴起源的慢上升余辉080710和260310A,而非样本主体,支持其相对论余辉起源。
在整个Kann样本中,只有GRB 080710的光变形态与AT 2020afjz最接近——但即便是GRB 080710,持续时间也没有AT 2020afjz这么长。还有另一个有趣的对比对象:GRB 260310A,这是一个异常长的余辉,它紧跟着明亮的Ic型超新星SN 2026fgk出现。AT 2020afjz的时间分辨演化正好介于这两个事件之间,暗示它们可能拥有相似的前身星系统。
有意思的是,GRB 080710和GRB 260310A都不是“标准”的GRB余辉。对GRB 080710的分析(T. Krühler et al. 2009)认为它要么是离轴GRB(off-axis GRB),要么是一个初始洛伦兹因子Γ₀小于100的正轴喷流——正好满足“脏火球”(dirty fireball)的判据。而对GRB 260310A的分析(Y.-H. Yang et al. 2026)则发现它后期余辉涌现、瞬时伽马辐射弱、峰值能量硬,这些特征都支持离轴GRB的起源。既然两个最相似的“亲戚”都指向非标准余辉,那AT 2020afjz是否也是同类?这就需要用物理模型来严格检验了。
VegasAfterglow贝叶斯拟合揭示双解
要对这样一个只有光变曲线、没有多波段数据的暂现源进行物理建模,听起来像是“巧妇难为无米之炊”。但好消息是,AT 2020afjz的光变曲线信噪比足够高、时间分辨率足够好,加上研究团队开发了一套稳健的拟合方法,让这件事变得可行。
建模工具是VegasAfterglow(Y. Wang et al. 2026),这是一个专门用于GRB余辉建模的物理框架。它最大的特色在于把喷流的初始体洛伦兹因子Γ₀作为直接模型参数,并且能够处理反向激波,支持多种喷流几何结构以及灵活的周介介质密度分布。对于观测数据来说,VegasAfterglow的前向激波模型有八个自由参数:初始体洛伦兹因子Γ₀、各向同性能动能量E_K,iso、周介介质密度(均匀ISM的n或恒星风的A★)、喷流核心角θ_c、观测角θ_obs、电子能谱指数p,以及电子和磁场的能量占比ε_e和ε_B。
采用Nautilus(J. U. Lange 2023)进行贝叶斯推断。Nautilus是一种基于重要性嵌套采样的算法,它用深度神经网络来建模等似然边界,相比传统MCMC能更高效地采样复杂的多峰后验分布,同时计算出贝叶斯证据。模型还加入了两个额外的自由参数:爆发时刻t₀和红移z——红移采用了EAZY联合拟合得到的光度红移分布作为先验。
图C3:长暴GRB的周介介质密度分布(均匀ISM介质),数据来自文献中的宽带余辉拟合:A. Panaitescu & P. Kumar (2002,金色)、M. D. Aksulu et al. (2022,蓝色)和S. B. Cenko et al. (2011,红色),共N=24个测量值。黑色曲线是合并样本在对数空间的高斯拟合(μ=0.02,σ=1.46),本文采用该分布(取整为(0.0, 1.45))作为AT 2020afjz余辉拟合的密度先验。表3:AT 2020afjz余辉拟合的先验分布(VegasAfterglow + Nautilus;仅前向激波,禁用侧向扩展)。先验符号说明:U表示线性均匀分布,log10U表示对数均匀分布,N表示正态分布,log10N表示对数正态分布。微观物理和几何先验由观测到的长暴GRB分布所约束。研究团队还测试了不同的喷流几何结构。VegasAfterglow框架里预设了顶帽(top-hat)喷流和高斯(Gaussian)喷流两种几何。对两种几何分别做了拟合后,表4的证据对比显示它们的贝叶斯证据差异极小(ΔlnB < 1),说明只靠TESS的光变曲线无法区分喷流几何。最终团队选用顶帽喷流作为基准模型继续分析。
表4:AT 2020afjz使用VegasAfterglow+Nautilus以4000个活动点建模的喷流几何贝叶斯证据对比。ln B是相对顶帽几何的值,后者是统计上略微优选的模型。仅凭TESS数据无法区分喷流几何。
拟合结果出来后,一个让团队既兴奋又头疼的现象出现了——后验分布是双峰的。用高斯混合模型在(θ_obs, Γ₀, log n)参数空间里拆分,得到了两个截然不同的物理解决方案。表5给出了这两个解的详细参数。
表5:AT 2020afjz全几何顶帽拟合的双群体后验(VegasAfterglow+Nautilus,4000个活动点)。恢复出两个解:正轴低Γ₀模型(M1)和离轴较高Γ₀模型(M2)。用(θ_obs, Γ₀, log n)上的双分量高斯混合来分离这两个解。两种模型对数据的拟合效果相同。数值为后验中位数及16-84百分位不确定度。图4:AT 2020afjz的正轴VegasAfterglow拟合联合图。左:Nautilus先验(灰色)和后验的角图,后验呈双峰分布;高斯混合分解将其拆分为模式1(M1;金色,正轴低Γ₀)和模式2(M2;蓝色,离轴窄喷流,Γ₀正常),各模式的位数和16-84百分位数在对角线上方。两者证据相等(Δln=0),代表Γ₀-几何简并,群体比例主要由M2的宽后验驱动。右上:两个模式中位数模型与数据(黑点)的对比:M1(金色实线)和M2(蓝色虚线),阴影为16-84%区域,下方为残差。爆炸时刻t₀被标出,红色区域为t₀附近Fermi-GBM的覆盖时间段。