[论文解读] SimProp v2r2: a Monte Carlo simulation to compute cosmogenic neutrino fluxes

Roberto Aloisio, Denise Boncioli+4 more|arXiv (Cornell University)|May 6, 2015

Astrophysics and Cosmic Phenomena参考文献 9被引用 6

一句话总结

SimProp v2r2 是一款更新的蒙特卡洛模拟代码，用于模拟宇宙射线在星系际辐射场中的超高端能量宇宙射线传播，全面考虑了在宇宙微波背景（CMB）和星系际背景光（EBL）上光致产生π介子的随机过程，并计算了宇宙成因中微子通量。其结果与Engel–Seckel–Stanev基准结果高度一致，特别是在预测次级中微子能谱方面，仅因忽略了次要的π介子产生通道而存在微小偏差。

ABSTRACT

We present an updated version of the SimProp Monte Carlo code: a simulation scheme to study the propagation of ultra-high-energy cosmic rays through diffuse extragalactic background radiation. The new version of the code presents two important updates: (i) it treats in a full stochastic approach all interaction channels involving ultra-high-energy cosmic rays and (ii) it takes into account the production of secondary (cosmogenic) neutrinos. This new version of SimProp was tested against different simulations code, in particular the production of secondary neutrinos was compared with the fluxes expected in the scenario of Engel, Seckel and Stanev.

研究动机与目标

• 开发一种全面的蒙特卡洛模拟方法，用于超高端能量宇宙射线传播，实现对相互作用的完整随机处理。
• 计算宇宙射线在星系际辐射场中传播过程中产生的宇宙成因中微子通量。
• 通过在宇宙微波背景（CMB）和星系际背景光（EBL）上包含光致产生π介子的随机处理，改进先前版本。
• 将新中微子产生模块与既有的理论通量进行验证，特别是与Engel–Seckel–Stanev情景进行对比。
• 提供一个公开可用、开源的工具，用于模拟超高端能量宇宙射线传播中的次级粒子产生。

提出的方法

• 该代码采用随机方法模拟质子和原子核在CMB和EBL上的光致产生π介子过程，取代了以往的连续近似。
• 利用连续能量损失近似模型能量损失，包括绝热膨胀和电子-正电子对产生。
• 粒子沿红移空间中的直线路径传播，时间与位置仅通过红移z进行追踪。
• 次级粒子（包括中微子）在相互作用点随机生成，并追踪至地球，其中中微子能量因宇宙学红移而按(1+z)缩放。
• 输出包含关于初级注入、次级粒子能量、种类和产生红移的详细事件级数据，以ROOT格式存储。
• 中微子味信息在产生时记录，但模拟中未包含中微子振荡。

实验结果

研究问题

• RQ1 随机蒙特卡洛模拟在多大程度上能准确再现超高端能量宇宙射线相互作用产生的宇宙成因中微子通量？
• RQ2 与仅考虑CMB相比，在CMB和EBL上同时包含光致产生π介子过程，对次级中微子能谱有何影响？
• RQ3 模拟得到的中微子通量与Engel、Seckel和Stanev的基准结果在定量上如何比较？
• RQ4 忽略次要的π介子产生通道对预测的电子反中微子能谱有何影响？
• RQ5 在标准宇宙学假设下，模拟输出与理论预期的一致程度如何？

主要发现

• 模拟得到的宇宙成因中微子通量与Engel–Seckel–Stanev基准结果高度一致，尤其在μ子中微子和反μ子中微子方面。
• SimProp v2r2中的电子反中微子能谱在约10^18 eV处缺少第二个峰，这归因于忽略了次要的π介子产生通道。
• 由于中微子振荡将在地球处使不同味的中微子通量趋于相等，而μ子反中微子主导总反中微子通量，因此该峰的缺失在观测上可忽略不计。
• 该模拟成功实现了对CMB和EBL上光致产生π介子的随机处理，优于早期的连续近似方法。
• 代码的输出结构支持对次级粒子产生的详细分析，包括能量谱、味组成以及红移依赖的产额。
• 结果证实，宇宙成因中微子通量可被该模拟稳健预测，验证了其在未来的UHECR和中微子天体物理学研究中的适用性。