[论文解读] Old Data, New Forensics: The First Second of SN 1987A Neutrino Emission

Shirley Weishi Li, J. F. Beacom+2 more|arXiv (Cornell University)|Jun 13, 2023

Neutrino Physics Research被引用 5

一句话总结

本研究首次对现代1、2和3维核心坍缩超新星模型与SN 1987A中微子数据进行了全面比较，重点关注碰撞后第一秒内的数据。尽管建模方法一致，大多数模型与观测到的低中微子计数存在偏差，且中微子振荡或前身星质量变化均无法缓解这一矛盾，凸显了未来需改进模拟与数据分析，以应对银河系下一次超新星爆发的迫切需求。

ABSTRACT

The next Milky Way supernova will be an epochal event in multi-messenger astronomy, critical to tests of supernovae, neutrinos, and new physics. Realizing this potential depends on having realistic simulations of core collapse. We investigate the neutrino predictions of nearly all modern models (1-, 2-, and 3-d) over the first $\simeq$1 s, making the first detailed comparisons of these models to each other and to the SN 1987A neutrino data. Even with different methods and inputs, the models generally agree with each other. However, even considering the low neutrino counts, the models generally disagree with data. What can cause this? We show that neither neutrino oscillations nor different progenitor masses appear to be a sufficient solution. We outline urgently needed work.

研究动机与目标

• 通过比较现代1、2和3维核心坍缩超新星模拟，评估建模不确定性。
• 通过检验这些模型与实际SN 1987A中微子数据的匹配程度，评估物理不确定性。
• 确定中微子振荡或前身星质量变化是否能缓解模型与数据之间观测到的矛盾。
• 在银河系下一次核心坍缩超新星爆发之前，激发超新星建模与中微子数据分析的创新工作。

提出的方法

• 作者汇编并比较了多种最先进的1、2和3维超新星模拟的中微子通量预测，涵盖不同维度、前身星质量及中微子振荡情景。
• 他们以SN 1987A的KamLAND-II（Kam-II）数据作为基准，聚焦于核心碰撞后约1秒内的数据。
• 在模型-数据比较中，他们基于总事件数和能量谱计算p值，使用源自蒙特卡洛模拟的检验统计量。
• 他们分析了不同中微子振荡情景——正常质量顺序（NH）、倒质量顺序（IH）及味平衡——对预测事件率与谱的影响。
• 比较中包括了0.5秒和1.3秒的截断时间，以评估时间范围选择的影响。
• 他们使用标准振荡公式计算振荡后的电子反中微子通量：$\tilde{f}_{\bar{\nu}_e} = f_{\bar{\nu}_e} P(\bar{\nu}_e \to \bar{\nu}_e) + f_{\bar{\nu}_x} [1 - P(\bar{\nu}_e \to \bar{\nu}_e)]$，其中$P(\bar{\nu}_e \to \bar{\nu}_e) \approx 0.7$（NH）、0（IH）和1/3（平衡）。

Figure 1: Neutrino ( $\bar{\nu}_{e}$ ; others shown in S.M.) luminosity and RMS energy profiles from supernova simulations O’Connor et al. ( 2018 ); Bruenn et al. ( 2013 , 2016 ); O’Connor and Couch ( 2018a ); Summa et al. ( 2016 ); Kotake et al. ( 2018 ); Vartanyan et al. ( 2018 ); Ott et al. ( 201

实验结果

研究问题

• RQ1 现代1、2和3维超新星模拟在核心碰撞后第一秒内是否产生一致的中微子预测？
• RQ2 考虑到观测到的事件数较少，这些模型在多大程度上与实际SN 1987A中微子数据一致？
• RQ3 在不同质量顺序情景下，中微子振荡是否能缓解模型预测与SN 1987A数据之间的差异？
• RQ4 改变前身星质量（例如10–20 M⊙）是否显著影响第一秒内中微子发射的模型-数据一致性？
• RQ5 与真实中微子爆发数据相比，当前超新星模拟中的主要不确定性来源是什么？

主要发现

• 尽管在维度和输入物理方面存在差异，现代1、2和3维超新星模拟在核心碰撞后约1秒内的预测中微子通量方面总体保持一致。
• 模型普遍与SN 1987A数据不一致，当截断时间为1.3秒时，大多数模型的p值低于0.05阈值。
• 截断时间的选择显著影响p值：使用0.5秒可降低统计矛盾，但1.3秒截断揭示了持续存在的模型-数据不一致。
• 中微子振荡无法缓解矛盾；在某些情况下，尽管计数p值改善，但能量谱拟合反而变差，这是由于检验统计量中复杂的相互作用所致。
• 引入振荡会使预测的$¯{\nu}_e$谱向更高能量和更低通量移动，但p值的总体影响主要由计数不匹配主导，导致整体改善效果不明确。
• 本研究识别出迫切需要开展更长时间的模拟和更全面的模型比较，以解决与SN 1987A数据持续存在的差异。

Figure 2: Predicted counts and average energies of supernova models (colors as in Fig. 1 ), compared to each other and to SN 1987A data. Our main calculations are in the inset.