Signatures of $^{56}$Ni Mixing and Neutron-rich Ejecta in Supernovae

技术摘要：$^{56}\text{Ni}$ 混合与富中子抛射物的特征

问题陈述
超新星（SN）光变曲线通常使用单区放射性衰变模型（例如 Arnett 1980, 1982）进行解释，以推断前身星性质，如抛射物质量（$M_{\text{ej}}$）、镍质量（$M_{\text{Ni}}$）和爆炸能量（$E_k$）。这些模型假设存在单一的有效扩散时间尺度，并忽略了放射性物质的空间分布。然而，流体动力学模拟表明，$^{56}\text{Ni}$ 在抛射物深处合成，随后通过瑞利-泰勒不稳定性、对流和湍流向外混合。这种空间分布是一个关键的物理变量，而单区模型忽略了这一点。因此，对于高度混合的抛射物，其光变曲线可能会被误解为具有较低的抛射物质量，或者需要额外的中心引擎功率（例如磁星或回落吸积）来解释快速上升和高亮度。此外，坍缩星（collapsars）中 r-过程富集的可观测性，高度取决于中子富集抛射物相对于镍驱动背景的放置位置和几何结构，而这是一个在球对称模型中经常被简化的因素。

方法论
作者引入了 snmix，这是一个基于 JAX 的 Redback 框架插件实现的半解析多壳层模型。与单区模型不同，snmix 将同向膨胀的抛射物离散化为 $N$ 个具有几何间距速度的壳层。

• 放射性结构： 该模型允许任意的 $^{56}\text{Ni}$ 径向分布（通过顶帽型或“核+尾”剖面的混合分数 $f_{\text{mix}}$ 进行参数化）和 r-过程物质。
• 输运物理： 每个壳层独立演化，追踪热能、绝热损失和辐射扩散。至关重要的是，模型根据每个特定壳层外部的累积柱密度，而非全局平均值，来计算扩散和伽马射线沉积。
• 不透明度与 SED： 模型结合了灰色不透明度和随温度变化的不透明度律。对于 r-过程计算，它包含了镧系元素依赖的不透明度和经过核网络计算校准的解析加热率。
• 几何结构： 虽然主要针对球对称，但该模型扩展到了坍缩星的二分量“盘风”（disk-wind）几何结构，结合了极向 SN 分量和赤道中子富集风，并根据观测视角进行加权。
• 推断： 作者使用该框架生成合成玻利米度光变曲线，并进行贝叶斯推断（使用 NUTS 采样），将单区拟合结果与多壳层真实情况进行对比，分析参数偏差。

核心贡献与结果

1. $^{56}\text{Ni}$ 混合对光变曲线的影响：

   • $^{56}\text{Ni}$ 的向外混合减少了上覆的扩散柱密度和伽马射线捕获柱密度。对于固定的 $M_{\text{Ni}}$、$M_{\text{ej}}$ 和 $E_k$，这会导致更快速、更亮的上升过程以及更陡峭的放射性尾部。
   • 快速且明亮的上升并非低抛射物质量或存在引擎的专属证据；高度混合的放射性模型可以在没有额外动力源的情况下重现这些特征。

2. 单区推断中的偏差：

   • 使用单区模型拟合高度混合的光变曲线可以得到视觉上良好的拟合，但会产生有偏的物理参数。
   • 在固定不透明度时，模型通过推断显著降低的 $M_{\text{ej}}$ 和更高的镍质量分数（$f_{\text{Ni}}$）来吸收缺失的径向结构，而 $M_{\text{Ni}}$ 则保持相对稳定（在固定不透明度下恢复精度约为 $\sim3\%$）。
   • 如果允许不透明度变化，偏差会转向推断出极低的有效不透明度（对于完全混合的情况，$\kappa_{\text{fit}}/\kappa_{\text{input}} \simeq 0.24$）。
   • 这些偏差会传播到对爆炸能量和前身星映射的错误推断中。

3. 案例研究：SN 1998bw：

   • 将 snmix 应用于 SN 1998bw 得到中等混合分数（$f_{\text{mix}} \approx 0.59$）。虽然单区最大似然估计值与后验分布重叠，但在 $M_{\text{ej}}$ 上差异约为 $13\%$，在 $E_k$ 上差异约为 $5\%$，这说明仅靠玻利米度数据无法唯一分解质量、能量、不透明度和混合度。

4. 坍缩星 r-过程特征：

   • r-过程富集的特征并不总是表现为晚期的近红外（NIR）过剩。
   • 放置位置： 埋藏的 r-过程物质最初会因高不透明度而抑制光学光度，只有当富集层变得透明后，NIR 过剩才会出现。
   • 几何结构： 在盘风模型中，如果中子富集物质局限于赤道平面，则轴向 GRB-SNe（极向视角）的 r-过程特征可以忽略不计，因为极向 SN 分量占据主导。信号在离轴观测者处最强，或者当物质被混合到较高纬度（$\gtrsim 30^\circ$）时最为明显。
   • 诊断方法： 色彩演化（$\Delta(r-H)$）是一种有价值的诊断工具，但需要将光学抑制与 H 波段增亮区分开来。成对的光学和近红外监测（10–60 天）比单纯的晚期 H 波段观测点更具信息量。

意义
本文认为，“放射性结构”（加热和不透明度的径向分布）是一个关键且经常被忽视的变量，在超新星解释中至关重要。

• 对于镍驱动的超新星： 它挑战了常规使用中心引擎来解释快速上升的 Ic-BL 型和 GRB-SNe 的做法，指出向外混合是一种可行的替代方案。它警告说，单区拟合可能会提供“良好”的统计拟合，但提供的物理参数却是错误的。
• 对于坍缩星： 它完善了寻找 r-过程富集的策略，超越了简单的“NIR 过剩”启发式方法。它强调，如果抛射物是赤道受限的，那么在轴向 GRB-SNe 中未探测到 r-过程信号并不代表排除了 r-过程的产生。
• 观测策略： 作者得出结论，打破这些简并性需要结合多种观测指标：光球速度、色彩演化、光谱以及晚期玻利米度下降，这些对于区分低抛射物质量、高混合度或引擎功率是必要的。

这项工作强调，尽管 snmix 是一个简化的半解析工具，但它有效地隔离了径向结构的影响，为标准的单区解释提供了必要的修正。