[论文解读] Non-spherical Core Collapse Supernovae I. Neutrino-Driven Convection, Rayleigh-Taylor Instabilities, and the Formation and Propagation of Metal Clumps

K. Kifonidis, T. Plewa+2 more|CERN Bulletin|Feb 12, 2003

Gamma-ray bursts and supernovae参考文献 81被引用 192

一句话总结

本文通过二维流体动力学模拟研究了大质量恒星核心塌缩超新星，结果表明中微子驱动的对流和瑞利-泰勒不稳定性在激波恢复后约几分钟内将各向异性、低熵的富金属壳层破碎为团块。这些团块在约300秒后与流场解耦，以3500–5500 km/s的速度作弹道运动；与II型模型相比，Ib型模型更符合观测到的团块速度，这归因于前身星包层结构及反向激波效应的差异。

ABSTRACT

Two-dimensional simulations of a Type II and a Type Ib-like supernova explosion are presented that encompass shock revival by neutrino heating, neutrino-driven convection, explosive nucleosynthesis, the growth of Rayleigh-Taylor instabilities, and the propagation of newly formed metal clumps through the exploding star. In both cases we find very high Ni56 velocities of 17000 km/s shortly after shock-revival, and a complete fragmentation of the progenitor's metal core within the first few minutes after core bounce, due to the growth of Rayleigh-Taylor instabilities at the Si/O and (C+O)/He composition interfaces. This leads to the formation of high-velocity, metal-rich clumps which eventually decouple from the flow and move ballistically through the ejecta. Maximum final metal velocities of 3500-5500 km/s and 1200 km/s are obtained for the Type Ib model and the Type II model, respectively. The low maximum metal velocities in the Type II model, which are significantly smaller than those observed in SN 1987A, are due to the massive hydrogen envelope of the progenitor. The envelope forces the supernova shock to decelerate strongly, leaving behind a reverse shock below the He/H interface, which interacts with the clumps and slows them down significantly. This reverse shock is absent in the Type Ib-like model. The latter is in fairly good agreement with observations of Type Ib supernovae. In addition, in this case the pattern of clump formation in the ejecta is correlated with the convective pattern prevailing during the shock-revival phase. This might be used to deduce observational constraints for the dynamics during this early phase of the evolution, and the role of neutrino heating in initiating the explosion.

研究动机与目标

• 研究中微子驱动的对流和瑞利-泰勒不稳定性在核心塌缩超新星爆发期间对喷流各向异性和金属团块形成的影响。
• 在相同的前身星框架下，比较II型和Ib型爆炸的流体动力学特性与混合行为。
• 确定早期多维流体动力学过程如何影响最终金属团块的速度及其对合成光谱和光曲线的启示。
• 探讨喷流中金属团块的结构是否可用于推断激波恢复期间中微子驱动对流的动力学特征。
• 评估反向激波和前身星包层结构对团块减速和最终速度的影响，特别是与观测到的超新星遗迹的关系。

提出的方法

• 对15 M⊙的蓝超巨星前身星进行二维流体动力学模拟，包含中微子加热、激波恢复和核合成过程。
• 采用Woosley和Bruenn提供的前身星和爆发后模型作为初始条件，以确保物理解析一致性。
• 追踪硅/氧和（碳+氧）/氦界面处由中微子驱动对流引发的瑞利-泰勒不稳定性增长。
• 模拟金属富集团块在喷流中的形成、传播与减速过程，考虑阻力与压力梯度的影响。
• 通过从同一前身星中移除氢包层，比较II型与Ib型类似模型，以隔离包层效应。
• 分析团块速度与混合模式的演化，特别关注从流体耦合到弹道解耦的转变过程。

实验结果

研究问题

• RQ1 中微子驱动的对流和瑞利-泰勒不稳定性如何影响核心塌缩超新星中富金属团块的形成与破碎？
• RQ2 前身星包层结构在决定激波恢复后金属团块的最终速度与减速行为中起什么作用？
• RQ3 为何II型与Ib型类似超新星尽管初始爆炸机制相似，却表现出不同的最终金属团块速度？
• RQ4 喷流中团块的空间分布能否为激波恢复期间中微子驱动对流的动力学提供观测约束？
• RQ5 早期反向激波相互作用在多大程度上改变了金属团块的最终速度，这对合成光谱建模有何影响？

主要发现

• 富金属团块在激波恢复后约20秒内，由中微子驱动对流引发的瑞利-泰勒不稳定性，使各向异性、低熵的铁群壳层发生破碎。
• 团块在爆发后约300秒与流场解耦，开始弹道运动，II型与Ib型类似模型中的最大速度均达约3500–5500 km/s。
• 在II型模型中，团块速度在约1500秒后显著下降，这是由于在致密氢包层中受到强烈阻力及反向激波作用，至爆发后20,000秒时降至约1200 km/s。
• Ib型类似模型与观测到的Ib型超新星（如SN 1987A和SN 1987F）的团块速度更吻合，其速度在3500–5500 km/s范围内保持稳定。
• Ib型模型中最终的团块速度分布与激波恢复期间的对流结构相关，表明其可能成为早期爆发动力学的潜在观测诊断工具。
• 反向激波引起的团块减速是解释观测与模型间速度差异的关键因素，尤其在II型事件中，必须在光谱建模中予以考虑。