[论文解读] Long time supernova neutrino simulations

M. Mori|arXiv (Cornell University)|Jan 7, 2022

Astrophysics and Cosmic Phenomena参考文献 56被引用 23

一句话总结

本文提出了一套全面的、端到端的长时间超新星中微子辐射与探测模拟框架，模拟了核心坍缩、爆发及原中子星冷却过程，持续时间达20秒。对于距离10 kpc的银河系超新星，该框架预测在超级神冈探测器中可观测到约1,800个中微子事例，结果与SN 1987A的能量量级一致但事例数较低，凸显了改进通量建模的必要性以匹配观测数据。

ABSTRACT

Long time neutrino simulation data up to 20 seconds with no neutrino oscillation File names are formatted as below: model.A.bary.B.BB_grav.C.CC.dat A is a model number and B.BB is a baryonic mass and C.CC is a gravitational mass. Neutron star masses are measured at 20 seconds. Readme.txt explanation of each file model1.bary.1.27_grav.1.36.dat includes neutrino fluxes and luminosities from z9.6 in the Nakazato format (see below and http://asphwww.ph.noda.tus.ac.jp/snn/guide.pdf). See https://academic.oup.com/ptep/article/2021/2/023E01/6045986 for more details on the z9.6 supernova explosion. model2.bary.1.41_grav.1.31.dat includes neutrino spectra from a supernova which leave a neutron star of 1.41 Msun in baryonic mass and 1.31 in gravitational mass. model3.bary.1.22_grav.1.29.dat includes neutrino spectra from a supernova which leave a neutron star of 1.22 Msun in baryonic mass and 1.29 in gravititonal mass. sn_spectra.py a python script which interprets the nakazato format. If you run the script from command line, the script outputs fluxes and luminosities and average energies for each flavor. Nakazato format The Nakazato format contains time evolutions of neutrino spectra and is a suitable format for estimation of supernova neutrino signals on earth. "T0 E0 E1 dN_nue(T0)/dE1 dN_anue(T0)/dE1 dN_nux(T0)/dE1 dL_nue(T0)/dE1 dL_anue(T0)/dE1 dL_nux(T0)/dE1 E1 E2 dN_nue(T0)/dE2 dN_anue(T0)/dE1 dN_nux(T0)/dE2 dL_nue(T0)/dE2 dL_anue(T0)/dE2 dL_nux(T0)/dE2 . . . E19 E20 dN_nue(T0)/dE20 dN_anue(T0)/dE20 dN_nux(T0)/dE20 dL_nue(T0)/dE20 dL_anue(T0)/dE20 dL_nux(T0)/dE20 T1 E0 E1 dN_nue(T0)/dE1 dN_anue(T1)/dE1 dN_nux(T1)/dE1 dL_nue(T1)/dE1 dL_anue(T1)/dE1 dL_nux(T1)/dE1 . . . ",where Tn [s] is a time measured from the bounce and Ek [MeV] is a neutrino energy.

研究动机与目标

• 开发一个集成的模拟框架，用于超新星中微子辐射与地面探测的端到端建模。
• 将中微子模拟扩展至爆发后第一秒以外，解决SN 1987A中观测到的10秒以上发射现象。
• 实现理论模型与未来高统计量银河系超新星中微子观测之间的详细对比。
• 通过长期、一致的模拟，增进对中微子驱动爆炸机制和原中子星冷却过程的理解。
• 通过预测未来超新星中微子探测的事例率与能谱，支持探测器设计与分析。

提出的方法

• 开发了一维长时间超新星模拟，成功实现爆炸并演化至20秒。
• 将核心坍缩、爆发及原中子星冷却过程整合进一个统一的模拟框架中。
• 采用中微子输运与发射模型，计算所有味态的时间依赖中微子能谱。
• 利用反β衰变（IBD）和电子散射（ES）过程，模拟中微子在超级神冈探测器中的探测。
• 引入中微子振荡（正常与倒质量顺序）及探测器响应函数。
• 通过与SN 1987A观测结果对比，验证模型可靠性与通量标定效果。

实验结果

研究问题

• RQ1 核心坍缩超新星的长时间中微子光曲线与SN 1987A中观测到的10秒以上发射有何异同？
• RQ2 对于距离10 kpc的银河系超新星，20秒内超级神冈探测器中预测的中微事件率是多少？
• RQ3 模拟的中微子能谱与事件分布能否区分不同的中微子振荡质量顺序？
• RQ4 预测的平均能量与事例数与SN 1987A的实测数据相比如何？
• RQ5 理论通量模型需作何改进，才能匹配SN 1987A中观测到的更高事例数？

主要发现

• 模拟成功生成了持续至爆发后20秒的长时序中微子发射轮廓，与SN 1987A中观测到的延长发射一致。
• 对于距离10 kpc的超新星，模型预测在超级神冈探测器中可产生约1,800个中微子事例，主要由反β衰变（IBD）事例主导。
• 预测的平均能量演化与理论预期相符，整体能量量级与SN 1987A具有良好的一致性。
• 模型预测仅产生四个来自中子化爆发的事例，低于SN 1987A中观测到的11个事例，表明需在不改变平均能量的前提下提高总中微子通量。
• 电子散射（ES）事例可在探测器中与IBD事例区分，尤其在计划中的掺钆超级神冈（SK-Gd）探测器中可提供方向信息。
• 该框架实现了超新星物理、中微子发射与探测器响应的一致建模，为未来银河系超新星分析构建了稳健的模拟流程。