[论文解读] Latest Results from the Air Shower Simulation Programs CORSIKA and CONEX

T. Pierog, R. Engel+3 more|ArXiv.org|Feb 9, 2008

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

一句话总结

本文使用 CORSIKA 和 CONEX 对广延空气簇射进行了更新的模拟，引入了新的强子相互作用模型 epos 1.6 和 qgsjet II-3，结果表明 epos 模型预测的地表 muon 产额显著更高——比 qgsjet II-3 高出 30–40%，这对像 Pierre Auger 观测站这样的超高能宇宙射线实验的能量和成分重建产生影响。

ABSTRACT

Interpretation of EAS measurements strongly depends on detailed air shower simulations. The uncertainty in the prediction of shower observables for different primary particles and energies is currently dominated by differences between hadronic interaction models. The new models QGSJETII-3 and EPOS 1.6, which reproduce all major results of existing accelerator data (including detailed data of RHIC experiments for EPOS), have been implemented in the air shower simulation programs CORSIKA and CONEX. We show predictions of these new models and compare them with those from older models such as QGSJET01 or SIBYLL. Results for important air shower observables are discussed in detail.

研究动机与目标

• 评估新的强子相互作用模型 epos 1.6 和 qgsjet II-3 对 CORSIKA 和 CONEX 中广延空气簇射（EAS）模拟的影响。
• 评估不同强子模型对关键可观测量（如 muon 密度、簇射最大深度 Xmax 和阻滞能量转换因子）的影响。
• 通过整合基于近期加速器数据（包括 RHIC 结果）校准的模型，提升超高能宇宙射线实验模拟的准确性。
• 通过改进 CORSIKA 的追踪算法和倾斜深度计算，实现对倾斜和向上运动簇射的模拟。
• 通过采用现代构建工具（autoconf/automake）和改进与外部代码（如 FLUKA 和 HERWIG）的接口，更新模拟基础设施。

提出的方法

• 使用统一的强子相互作用建模框架，将 epos 1.6 和 qgsjet II-3 集成到 CORSIKA 和 CONEX 中。
• 实现改进的粒子追踪算法，支持所有天顶角（包括 89°）的倾斜/弯曲深度选项。
• 基于 Kokoulin 与 Bogdanov 的工作，引入 Sternheimer 密度修正，用于 muon 电离能量损失。
• 在 CORSIKA 中使用蒙特卡洛方法模拟完整的簇射级联过程，而在 CONEX 中通过蒙特卡洛采样求解级联方程。
• 应用更新的构建系统（autoconf/automake），以简化 CORSIKA 和 CONEX 的编译与配置流程。
• 使用质子和铁核初级粒子，对比不同模型下横向切伦科夫信号分布、Xmax 和阻滞能量转换因子。

实验结果

研究问题

• RQ1 新的强子相互作用模型 epos 1.6 和 qgsjet II-3 如何影响广延空气簇射中地表 muon 含量的预测？
• RQ2 这些模型对表面探测器（如 Pierre Auger 观测站）中切伦科夫光横向分布的影响是什么？
• RQ3 在超高能范围内，这些模型在预测质子和铁核初级粒子簇射最大深度（Xmax）方面有何差异？
• RQ4 新模型在多大程度上改变了从阻滞能量到总初级能量的转换因子？
• RQ5 更新后的模拟工具（CORSIKA 6.611 和 CONEX）在倾斜和向上运动簇射模拟中如何提升准确性和灵活性？

主要发现

• 对于 10^19 eV 的质子簇射，epos 1.6 模型预测在地表 1 km 处的 muon 密度比 qgsjet II-3 高出 30–40%，这是由于增强的重子-反重子对产生。
• 在 10^19 eV 能量下，epos 1.6 的阻滞能量转换因子比其他模型高出最多 3.5%，表明更多能量保留在 muon 和中微子中。
• 对于质子引发的簇射，epos 1.6 模型预测的地表 muon 数量多于使用 qgsjet II-3 模拟的铁核引发簇射，这挑战了关于 muon 含量与成分的传统假设。
• 在高达 10^19 eV 的范围内，质子和铁核簇射的平均 Xmax 在各模型间保持在 ±20 g/cm² 以内，但 epos 模型在此范围内表现出略高的展宽率。
• 在低于 10^18 eV 时，基于 Xmax 的分析显示 qgsjet II-3 预测的初级成分更重，而高于 10^18 eV 时，这一趋势反转，这是由于高多重性行为的差异所致。
• 更新后的 CORSIKA 6.611 和 CONEX 版本在极端 89° 倾斜簇射中也表现出纵向能量沉积剖面的一致性，验证了跨程序的一致性。