重新审视具有有限厚度喷射物的伽马射线暴早期余辉：来自XRF 080330和GRB 080710的启示

《Journal of High Energy Astrophysics》：Revisiting early afterglows of gamma-ray bursts with finite-thickness ejecta: Implications from XRF 080330 and GRB 080710

　　摘要：我们重新研究了XRF 080330和GRB 080710在爆发后数千秒内多波段余辉中观测到的无色峰及间断现象的物理起源。通过采用一种同时考虑有限喷射物厚度和广义外部密度分布的数值余辉模型，我们运用贝叶斯推断来估算描述这些事件的模型参数。分析表明，这两种事件中的逐渐上升趋势以及无色时间特征，更自然地可以用具有有限壳层厚度的喷流动态演化来解释，而非偏轴观测效应。对于两次爆发，推断出的喷射物初始径向宽度约为10^13厘米，这意味着中心引擎的活动时间尺度远远长于仅从瞬时伽马射线持续时间所能推断出的时间。综合这些结果可知，当适当考虑有限喷射物厚度时，早期余辉光曲线会受到过渡动态的强烈影响，从而为瞬时阶段与余辉阶段之间的物理联系提供了依据，同时也凸显了在解释早期余辉时简单应用薄壳近似方法的局限性。此外，贝叶斯模型比较结果表明，广义的爆发周围密度分布比传统的均匀或稳态风模型更符合实际，这表明事先将外部密度结构固定为理想化模型可能会掩盖有关恒星前体在爆发前活动的重要信息。

引言：伽马射线暴是电磁辐射最强烈的爆炸性现象之一，普遍认为其源自由致密中心引擎产生的相对论性喷流（例如G?tz等人，2018年；Kumar和Zhang，2015年；Piran，1999年）。在瞬时伽马射线辐射之后，通常会观测到宽频带余辉，这些余辉很可能是由喷流与爆发周围介质相互作用时，非热电子在相对论性外部冲击波作用下产生的同步辐射所形成的（Sari等人，1999年；Sari等人，1998年）。自2004年尼尔·盖赫雷斯斯威夫特天文台投入使用以来，人们能够在爆发触发后几分钟内就开始对早期余辉进行观测。这些观测揭示了不同频率带上多种多样的光曲线行为，而这些行为很难用斯威夫特时代之前建立的标准余辉模型来解释。例如，早期的X射线余辉往往表现出一种典型的行为——最初快速衰减，随后是缓慢衰减阶段，再进入正常衰减阶段（例如Liang等人，2007年；Nousek等人，2006年；Willingale等人，2007年；Zhang等人，2006年、2007年）——而这种行为无法用单一的外部前向冲击波模型来直接解释。非常早期的光学余辉也展现出极大的多样性，在伽马射线暴触发后约500秒时，其亮度可相差三到四个数量级（例如Gao等人，2015年；Kann等人，2010年；Liang等人，2013年；Panaitescu和Vestrand，2011年；Wang等人，2013年）。相当一部分事件在光学上较为暗弱（例如Roming等人，2006年），这与人们对于显著反向冲击波辐射的预期不符（Kobayashi，2000年；Kobayashi和Sari，2000年；Sari和Piran，1995年）。此外，一些观测结果表明，早期的X射线和光学余辉并不一定来自同一发射源（例如Guidorzi等人，2025年；Li等人，2015年；Liang等人，2009年；Liang等人，2008年；Oates等人，2007年；Panaitescu等人，2006年；Sato等人，2007年；Urata等人，2007年），这暗示了它们可能具有不同的物理起源（Ghisellini等人，2007年；Yamazaki，2009年）。

在这种背景下，具有高灵敏度软X射线探测能力的爱因斯坦探针卫星有望为伽马射线暴研究打开新的窗口（Yuan等人，2022年）。尤其是那些X射线闪光现象，虽然备受关注（Arefiev等人，2003年；Barraud等人，2003年；Gotthelf等人，1996年；Hamilton等人，1996年；Heise等人，2001年；in ’t Zand等人，2025年；Strohmayer等人，1998年），但其起源仍不确定（Yuan等人，2025年）。从观测上看，X射线闪光具有低光谱峰能量和微弱的伽马射线辐射，通常被视为与典型伽马射线暴相连的连续现象序列的一部分（D’Alessio等人，2006年；Sakamoto等人，2008年；Sakamoto等人，2005年）。了解X射线闪光的起源有助于深入理解典型伽马射线暴的本质。目前已提出了多种关于X射线闪光的理论模型：高红移伽马射线暴（Barraud等人，2003年；Heise等人，2001年）、从偏轴方向观测到的顶帽形喷流（Ioka和Nakamura，2001年；Yamazaki等人，2002年、2003年、2004b年）、宽开口角喷流（Donaghy，2006年；Lamb等人，2005年）、高洛伦兹因子下对比度较小的内部冲击波（Daigne和Mochkovitch，2003年）、失败的伽马射线暴或“脏火球”（Dermer等人，1999年；Dermer和Mitman，2004年；Huang等人，2002年）、以光球为主的火球（Drenkhahn和Spruit，2002年；Mészáros等人，2002年；Ramirez-Ruiz和Lloyd-Ronning，2002年）、结构化或不均匀的喷流（Rossi等人，2002年；Toma等人，2005年；Yamazaki等人，2004a年；Zhang等人，2004年；Zhang和Mészáros，2002年）、来自坍缩星喷流的外围辐射（Zhang等人，2004年），以及偏轴方向的“炮弹状”喷流（Dado等人，2004年）。随着爱因斯坦探针时代中检测到的X射线闪光数量迅速增加（Gao等人，2025年），X射线闪光和伽马射线暴的早期余辉将为这些相互竞争的理论模型提供严格的检验标准。

有一些事件在爆发后约10^3?4秒时，会在多个波段同时出现亮度峰值，形成所谓的无色峰（Liang等人，2013年；Liang等人，2010年；Ror等人，2025年）。这类晚期的无色峰很难用标准余辉框架来解释，因为在标准模型中，10^2秒以后的辐射被认为是由减速的外部冲击波产生的，而特征同步辐射频率的演化则预示着不同观测波段会出现不同时间的彩色峰（Sari和Piran，1999年；Sari等人，1999年；Sari等人，1998年）。XRF 080330和GRB 080710就是具有在10^3秒左右出现无色峰特征的典型事件（Guidorzi等人，2009年；Krühler等人，2009年）。在相关的观测和分析论文中（Guidorzi等人，2009年；Krühler等人，2009年），这两个事件都被定性为偏轴余辉（例如Granot等人，2002年；Granot等人，2005年；Lamb和Kobayashi，2017年；Ramirez-Ruiz等人，2005年；Ryan等人，2020年；Totani和Panaitescu，2002年）。然而，如果X射线闪光本质上是偏轴的典型伽马射线暴，那么将XRF 080330和GRB 080710中的无色峰解释为偏轴余辉，就意味着这两个事件都存在类似X射线闪光的瞬时辐射现象。这一观点与GRB 080710在观测上被归类为典型伽马射线暴的事实难以调和。尽管偏轴观测仍是一种合理的解释，但其他物理因素也可能产生无色的时间特征。对于GRB 080710，Obayashi等人（2024年）指出，偏轴机制并非最有可能的解释，反而提出了一种由喷射物动态驱动的轴向起源。另一种可能性是，非均匀的爆发周围介质和/或有限厚度的喷射物能够改变冲击波的动态，进而以无色方式影响早期余辉光曲线的时间演化（Kobayashi和Sari，2000年；Yi等人，2013年）。此外，最近的数值研究表明，厚壳结构的喷射物会导致与薄壳近似预测不同的冲击波演化过程，从而产生独特的早期余辉行为（Kusafuka和Asano，2025a年、2025b年；Kusafuka等人，2023年）。

在这项研究中，我们重新分析了XRF 080330和GRB 080710的早期余辉现象，重点探讨了有限喷射物厚度以及广义爆发周围介质密度分布的影响。我们的目标是量化早期余辉光曲线如何限制喷射物的初始条件、爆发周围环境以及外部冲击波模型中的微观物理参数，该模型同时考虑了有限喷射物厚度和幂律形式的爆发周围介质密度分布。尤其由于早期余辉建模不可避免地涉及多种潜在的物理因素以及严重的参数简并性问题，我们采用了贝叶斯框架，通过直接比较后验分布和贝叶斯证据，从而在统一的统计方法中确定哪些因素才是真正必需的（例如偏轴几何结构、有限厚度或爆发周围介质结构）。通过分析这两起具有相似无色时间特征但瞬时辐射特性存在显著差异的事件（Guidorzi等人，2009年；Krühler等人，2009年），我们旨在为早期余辉物理学提供有力的理论依据，同时深入了解恒星前体环境及喷流形成过程。本文的结构如下：第2节介绍由有限厚度喷射物产生的余辉辐射模型；第3节阐述数值计算方法及模型评估手段，我们采用了适用于此研究的Magglow模型（Kusafuka和Asano，2025a年；Kusafuka等人，2025年）；第4节呈现计算结果，第5节则用于讨论。在整篇文章中，我们采用符合一致性宇宙学的模型，其中H0=70kms^-1Mpc^-1，ΩΛ=0.7，ΩM=0.3。XRF 080330的红移为z=1.51，对应的光度距离为dL=3.4×10^28厘米，而GRB 080710的红移为z=0.845，对应的光度距离为dL=1.7×10^28厘米。

片段一：具有有限厚度的喷射物动力学及其产生的余辉辐射
在本节中，我们阐述了由有限喷射物厚度和非均匀爆发周围介质所导致的向前冲击波动力学变化，以及由此产生的余辉辐射，这两者都可能影响早期余辉的演化过程（例如Gill和Granot，2023年；Kusafuka等人，2025年；Sari，1997年；Sari和Piran，1995年；Yi等人，2013年）。遵循Kusafuka等人的研究思路，我们设定一个坐标系，将喷射流产生的中心引擎置于原点R=0处。然后，根据密度分布……

片段二：数值计算与贝叶斯推断
我们使用开源代码Magglow（Kusafuka和Asano，2025a年；Kusafuka等人，2025年）来计算余辉辐射，该代码基于磁弹余辉模型来计算外部冲击波产生的同步辐射。在本研究中，我们忽略了喷射物磁化的作用。一般来说，喷射物的磁化主要会影响流体的加速阶段以及反向冲击波的形成（Kusafuka和Asano……）

片段三：结果
在本节中，我们将展示针对XRF 080330和GRB 080710的贝叶斯推断结果。在结果的理论解读方面，我们采用了最大后验值估计方法，即指在贝叶斯推断过程中所得样本中，使后验概率达到最大的参数组合，后续的分析评估也是基于这一参数值进行的。

片段四：讨论
在本文中，我们运用贝叶斯推断方法，分析了XRF 080330和GRB 080710的多波段余辉数据，所采用的余辉模型既考虑了有限壳层厚度，也考虑了爆发周围介质的幂律密度分布。我们的研究表明，该模型能够再现这两种事件在光学/近红外波段出现的逐渐上升趋势以及无色峰现象。分析结果不支持Guidorzi等人（2009年）和Krühler等人（2009年）提出的偏轴余辉理论，即β=θobs/θj>1的情况。

片段五：结论
在本文中，我们研究了XRF 080330和GRB 080710在爆发后数千秒内多波段余辉中观测到的无色峰的物理起源。我们使用了包含有限喷射物厚度和幂律形式外部密度分布的半解析余辉计算代码Magglow，以及基于MultiNest程序实现的嵌套采样算法的贝叶斯推断方法，对这两个事件进行了模型比较和参数估算。

作者贡献说明：
Kaori Obayashi：撰写——审稿与编辑、撰写——初稿、验证、软件应用、研究工作、正式分析、概念构建。
Ryo Yamazaki：撰写——审稿与编辑、撰写——初稿、验证、监督、项目管理、研究工作、概念构建。
Yo Kusafuka：撰写——审稿与编辑、验证、软件应用、研究工作、概念构建。
Katsuaki Asano：撰写——审稿与编辑、撰写——初稿、验证、监督。

利益冲突声明：
作者声明不存在任何可能影响本文研究结果的已知财务利益或个人关系。

Kaori Obayashi|Ryo Yamazaki|Yo Kusafuka|Katsuaki Asano
日本神奈川县相模原市深之部5-10-1，青山学院大学物理科学系，邮编252-5258