R-process nucleosynthesis from magnetar giant flares in neutron star--white dwarf mergers: A unified picture for peculiar long gamma-ray bursts

技术摘要：中子星-白矮星并合中磁星巨型耀斑引发的 r-过程核合成

问题陈述
一类“异常”的长伽马射线暴（LGRBs），以 GRB 211211A 和 GRB 230307A 为典型，挑战了长暴与短暴的标准分类二分法。这些事件表现出典型长暴（核心坍缩模型）的长持续时间（$T_{90} > 2$ s），但同时具有典型的短伽马射线暴（SGRBs，来自紧凑双星并合）的多组分瞬发辐射特征（前驱、主爆发、延伸发射）以及 X 射线平台期和千新星（kilonova）信号。现有模型面临重大障碍：

1. 黑洞残骸模型： 涉及中子星–黑洞（NS–BH）或中子星–中子星（NS–NS）并合并留下黑洞引擎的模型，难以解释由单一中心引擎产生的三个不同发射阶段（主爆发 [MB]、延伸发射 [EE] 以及 X 射线平台期）以及特定的 MB–EE 凹陷（trough）。
2. 标准中子星–白矮星（NS–WD）并合模型： 虽然 NS–WD 并合可以通过超吸积和磁阻器（propeller）阶段重现瞬发辐射结构，但常规流体动力学模拟表明，由潮汐瓦解的白矮星碎片形成的吸积盘是富质子的。这种环境无法产生稳健的 r-过程重元素，从而导致无法解释这些事件中观察到的千新星特征（特别是 GRB 230307A 中观测到的 r-过程元素的光谱确认）。

研究方法
作者提出了一个基于预先存在的磁星（$B_s \gtrsim 10^{14}$ G）与大质量白矮星（$M_{WD} \gtrsim 1 M_\odot$）并合的统一图景。该方法整合了：

• 吸积动力学： 利用用于恒定熵吸积盘（CENW/CEEW 方案）的流体动力学模拟（Kaltenborn 等人，2023），对超吸积阶段（$\dot{M} \gtrsim 10^{-4} M_\odot \text{s}^{-1}$）及随后的磁阻器阶段进行建模。
• 磁星物理学： 应用发电机机制，在吸积阶段将环向磁场（$B_\phi$）放大至 $10^{15} - 5 \times 10^{16}$ G。这会导致重复的星震（starquakes）和磁星巨型耀斑（GFs）。
• 核合成机制： 采用 Cehula 等人（2024）的框架，其中巨型耀斑的初始硬脉冲驱动膨胀的 $e^\pm$-$\gamma$ 火球。这些火球冲击中子星壳层，通过 $\alpha$-富集冻结（$\alpha$-rich freeze-out）机制抛射重子质量，从而在电子分数较高（$Y_e \gtrsim 0.4$）的情况下合成 r-过程元素。
• 观测建模： 对 GRB 211211A 和 GRB 230307A 的多波长数据（X 射线、光学、红外、无线电）进行联合马尔可夫链蒙特卡洛（MCMC）拟合。该模型结合了标准的外部冲击余辉模型与半解析的磁星增强型千新星形式（用 r-过程加热取代 $^{56}$Ni 加热）。

核心贡献

1. 异常长暴的统一起源： 本文提出了一个自洽的框架，通过单一的 NS–WD 并合情景解释了瞬发辐射结构（MB、EE、凹陷）、千新星信号以及 X 射线平台期。
2. 解决 r-过程缺失问题： 作者发现了一个新型的 NS–WD 并合中的 r-过程核合成通道。他们认为，r-过程元素并非依赖于富质子的白矮星碎片吸积盘，而是通过磁星巨型耀斑驱动的壳层抛射物进行合成。这解决了 NS–WD 并合模型与 GRB 230307A 中观测到的 r-过程元素光谱证据之间的矛盾。
3. “耀斑森林”： 主爆发被重新诠释为不是单次连续的喷流，而是约 800 到 2,000 个连续的巨型耀斑初始脉冲的叠加。这种解释使观测到的最小变率时间尺度（$\tau_{min} \sim 9-16$ ms）与银河系及邻近星系中观测到的单个巨型耀斑脉冲持续时间相吻合。
4. 定量约束： 研究提供了驱动观测到的千新星所需的抛射物质量（$M_{ej} \sim 10^{-5} - 10^{-3} M_\odot$）、膨胀速度和不透明度的具体约束，显示其与流体动力学及核反应网络计算结果一致。

结果

• 瞬发辐射： 模型成功地将 MB 再现为巨型耀斑脉冲的叠加，并将 EE 再现为磁阻器阶段的结果。从吸积到磁阻器的转变自然地解释了 MB 与 EE 之间的观测凹陷。
• 千新星属性： 联合建模得到的抛射物质量为 $M_{ej} \sim 10^{-2} M_\odot$（考虑了束缚效应），其速度与模拟结果一致。对于 GRB 230307A，高不透明度（$\kappa \sim 4.28 \text{ cm}^2 \text{ g}^{-1}$）由第二和第三峰值 r-过程核素的合成在物理上得到了证明。
• 磁星参数： 推断出的并合后磁星表面磁场（$B_{s,post} \sim 4 - 7 \times 10^{14}$ G）和初始自转周期，与并合过程中环向磁场预期的放大效应一致。
• 能量预算： 巨型耀斑脉冲的累积能量比已知的银河系磁星耀斑高出 2–4 个数量级，这意味着其磁场（$B_\phi \sim \text{数} \times 10^{16}$ G）比演化后的磁星显著更强。

意义与主张
本文声称提供了一个“统一图景”，解决了关于 GRB 211211A 和 GRB 230307A 等异常长暴起源的长期争论。通过将瞬发辐射、余辉和千新星都归因于单一的 NS–WD 并合事件（涉及前并合磁星），该模型为以下方面提供了连贯的解释：

• 瞬发辐射的多组分时间结构。
• 由磁星自转能损耗驱动的 X 射线平台期。
• r-过程重元素的产生及由此产生的千新星，克服了以往 NS–WD 并合模型在核合成方面的局限性。

作者指出，尽管该模型是稳健的，但它依赖于“中子星壳层在活跃吸积期间仍能保持足够低温以支持巨型耀斑机制”这一假设，这是作者承认需要在未来工作中进一步研究的一个局限性。