[论文解读] Corsika+Herwig Monte Carlo Simulation of Neutrino Induced Atmospheric Air Showers

M. Ambrosio, C. Aramo+5 more|arXiv (Cornell University)|Feb 28, 2003

Astrophysics and Cosmic Phenomena参考文献 2被引用 5

一句话总结

本文提出了一种改进版的CORSIKA蒙特卡罗簇射模拟器，通过与HERWIG事件生成器耦合，首次将中微子作为初级粒子引入。关键结果表明，电子中微子引发的簇射表现出显著增强的纵向剖面，相较于质子簇射；而μ子中微子簇射则被强烈抑制，从而在像皮埃尔·让弗朗索瓦天文台这样的超高能中微子探测实验中具备潜在的味分辨能力。

ABSTRACT

High-energy neutrino astronomy represents an open window both on astrophysical mechanisms of particle acceleration and on fundamental interactions. The possibility of detecting them in large earth-based apparatus, like AUGER, AMANDA, ANTARES, is quite challenging. In view of this, the capability of generating reliable simulations of air showers induced by neutrinos is mandatory in the analysis of experimental data. In this paper we describe preliminary results towards the development of a new version of the Monte Carlo CORSIKA, capable of handling neutrinos too as primary particles. In our approach the first interaction of the primary neutrino is simulated in CORSIKA with a call to the HERWIG event generator.

研究动机与目标

• 开发一种可靠的蒙特卡罗模拟工具，能够模拟超高能中微子引发的大气簇射。
• 扩展CORSIKA簇射模拟器，使其能够将中微子作为初级粒子处理，克服现有版本的局限性。
• 基于簇射发展特征，实现对中微子味（νₑ 和 νₘᵤ）的精确区分。
• 支持下一代中微子与宇宙射线观测站（如皮埃尔·让弗朗索瓦天文台、AMANDA 和 IceCube）的数据分析。
• 研究通过荧光探测器和地面阵列探测器在大天顶角下探测中微子引发簇射的可行性。

提出的方法

• 通过用HERWIG蒙特卡罗事件生成器调用替代标准初级相互作用模型，将中微子相互作用集成到CORSIKA中，用于首次中微子相互作用。
• 利用HERWIG在高能下模拟硬散射过程（例如，νₑ + d → e⁻ + u，νₘᵤ + u → μ⁻ + d），充分发挥其在部分子层次物理中的精确性。
• 将HERWIG输出的部分子末态粒子输入CORSIKA，以完整模拟大气中的级联发展过程，包括强子和电磁簇射演化。
• 在大天顶角（θ = 70°–80°）下模拟簇射，以最大化荧光产额和地面阵列信号，符合上行或掠射事件的典型特征。
• 采用经过修改的CORSIKA版本，通过将初始相互作用作为来自HERWIG的“伪初级”输入来处理中微子引发的簇射。
• 使用能量沉积、纵向剖面和横向分布函数作为关键可观测量，比较中微子与质子引发簇射的形态差异。

实验结果

研究问题

• RQ1 在超高能下，中微子引发的大气簇射的纵向剖面与质子引发的簇射有何不同？
• RQ2 能否基于大气中簇射的发展特征和能量沉积特性，区分电子中微子与μ子中微子簇射？
• RQ3 荧光探测器和地面阵列探测器中，中微子引发簇射的可探测阈值能量是多少？
• RQ4 初级中微子的味如何影响簇射的横向分布和总能量沉积？
• RQ5 中微子的中性流与带电费相互作用在多大程度上改变簇射形态和可探测性？

主要发现

• 在10¹⁵ eV能量下，电子中微子引发的簇射表现出显著增强的纵向剖面，相较于质子簇射，表明νₑ的探测阈值更低。
• μ子中微子引发的簇射在纵向剖面振幅上相对于质子簇射表现出强烈抑制，暗示νₘᵤ的探测阈值更高。
• νₑ簇射的平均能量沉积高于质子簇射，尤其在荧光探测范围，这是由于末态电子带来的增强电磁分量所致。
• 对于νₑ簇射，带电费相互作用产生的簇射剖面更强烈且更集中，而中性流相互作用则能量较低、更弥散。
• 横向分布函数显示，在相同能量和天顶角下，νₑ簇射在簇射核心附近产生的粒子密度高于质子簇射。
• 在10¹⁵ eV时，部分子分布函数可接受范围外的事件比例从0.02%上升至10²¹ eV时的40%，凸显了在极端能量下对PDF外推改进的必要性。