[论文解读] Non-Spherical Core-Collapse Supernovae II. Late-Time Evolution of Globally Anisotropic Neutrino-Driven Explosions and Implications for SN 1987A

K. Kifonidis, T. Plewa+3 more|ArXiv.org|Nov 11, 2005

Gamma-ray bursts and supernovae参考文献 1被引用 136

一句话总结

本研究对一颗非旋转的15 M☉蓝超巨星核心坍缩超新星进行了二维流体动力学模拟，表明全局各向异性、低阶模态驱动的爆炸（主要由l=2和l=1模态主导）能够再现SN 1987A的关键特征。该模型通过早期激波形变和反向激波相互作用的抑制，实现了高速铁族元素喷射（最高达3300 km/s）、强He/H界面混合（H被混合至500 km/s）以及拉长形喷流各向异性（主轴/次轴≈1.6），支持中微子驱动爆炸机制，无需引入奇异物理。

ABSTRACT

Two-dimensional simulations of strongly anisotropic supernova explosions of a nonrotating 15 solar mass blue supergiant progenitor are presented, which follow the hydrodynamic evolution from times shortly after shock formation until hours later. It is shown that explosions which around the time of shock revival are dominated by low-order unstable modes (i.e. by a superposition of the l=2 and l=1 modes, in which the former is strongest), are consistent with all major observational features of SN 1987A, in contrast to models which show high-order mode perturbations only and were published in earlier work. Among other items, the low-mode models exhibit final iron-group velocities of up to 3300 km/s, strong mixing at the He/H composition interface, with hydrogen being mixed downward in velocity space to only 500 km/s, and a final prolate anisotropy of the ejecta with a major to minor axis ratio of about 1.6. The success of low-mode explosions with an energy of about 2x10**51 erg to reproduce these observed features is based on two effects: the (by 40%) larger initial maximum velocities of metal-rich clumps compared to our high-mode models, and the initial global deformation of the shock. The latter triggers the growth of a strong Richtmyer-Meshkov instability at the He/H interface that results in a global anisotropy of the inner ejecta at late times (i.e. t > 10000 s), although the shock itself has long become spherical by then. The simulations suggest a coherent picture, which explains the observational data of SN 1987A within the framework of the neutrino-driven explosion mechanism using a minimal set of assumptions. It is therefore argued that other paradigms, which are based on (more) controversial physics, may not be required to explain this event. (abbreviated)

研究动机与目标

• 探究全局各向异性、低阶模态驱动的中微子驱动爆炸是否能再现SN 1987A的观测特征。
• 解决先前模型（尤其是高阶模态模拟）无法再现SN 1987A中大尺度混合与高速喷流的问题。
• 考察早期激波形变与初始速度结构在保护高速富金属团块免受反向激波破坏中的作用。
• 评估SN 1987A中观测到的拉长各向异性和向内H混合是否可由早期触发的流体不稳定性解释。
• 检验中微子驱动爆炸范式是否在最小假设下可完全解释SN 1987A的观测数据，而无需引入有争议的物理机制。

提出的方法

• 从核心反弹后20毫秒至反弹后5小时以上，对15 M☉蓝超巨星前身星进行了二维、高分辨率流体动力学模拟。
• 采用改进的中微子加热方案，初始通量较低且衰减较慢，与现代Boltzmann输运模拟一致，以延长爆炸 timescale。
• 在单元界面使用通量分裂格式（AUSM+）计算数值通量，结合迎风重构与界面声速通过临界声速条件计算，以确保激波的锐利解析。
• 通过分析t≈100 s时He/H界面处的涡量沉积，追踪流体不稳定性（特别是Richtmyer-Meshkov不稳定性）的增长。
• 利用勒让德模态分解（l=1和l=2）监测激波形变与喷流各向异性的演化，以量化非球对称性。
• 整合详细的核合成计算，以追踪喷流中铁族元素的速度与混合模式。

实验结果

研究问题

• RQ1 在中微子驱动爆炸中，低阶全局各向异性模态（l=1和l=2）是否能再现SN 1987A中观测到的高速铁族元素喷流（最高约3300 km/s）？
• RQ2 早期激波形变如何通过延迟与反向激波的相互作用，影响高速富金属团块的存活与速度？
• RQ3 He/H界面处的涡量沉积在激波变为球形后，如何在内层喷流中产生晚期各向异性？
• RQ4 为何高阶模态模型无法再现SN 1987A中的观测混合与各向异性？其物理机制差异为何？
• RQ5 中微子驱动爆炸机制本身是否足以完全解释SN 1987A的所有主要观测特征，而无需引入更富争议的物理范式？

主要发现

• 低模态、全局各向异性爆炸可实现高达约3300 km/s的最终铁族元素速度，与SN 1987A的观测一致。
• 激波的初始全局形变导致t≈100 s时在He/H界面处早期沉积涡量，触发强烈的Richtmyer-Meshkov不稳定性，从而驱动晚期各向异性。
• 内层喷流的拉长各向异性在t≈10,000 s时达到主轴/次轴比约1.6，与观测到的非球对称性一致。
• 由于不稳定性驱动的混合，氢在速度空间中被向下混合至约500 km/s的低速区域，与SN 1987A的观测成分结构一致。
• 富金属团块的初始最大速度比高模态模型高出约40%，使其在整个传播过程中保持亚声速，避免了反向激波相互作用导致的能量耗散。
• 模拟表明，仅通过最小假设的中微子驱动爆炸机制，即可完全解释SN 1987A的关键观测特征，表明无需引入更富争议的物理范式。