[论文解读] Study of the muon content of very high-energy EAS measured with the KASCADE-Grande observatory

J. C. Arteaga-Velázquez, W.D. Apel+55 more|arXiv (Cornell University)|Aug 14, 2013

Astrophysics and Cosmic Phenomena被引用 5

一句话总结

本研究利用 KASCADE-Grande 望远镜的观测数据，分析了 10^16–10^18 eV 能量区间的广延大气簇射（EAS）中的μ子含量，将实测的μ子密度分布与吸收系数与 QGSJET II、QGSJET II-04、SIBYLL 2.1 和 EPOS 1.99 宇宙射线强子相互作用模型的预测进行比较。发现实验数据与所有模型之间存在显著偏差，尤其是在μ子密度随大气深度的变化行为上，表明当前高能强子相互作用模型存在缺陷。

ABSTRACT

The KASCADE-Grande detector is an air-shower array devoted to the study of primary cosmic rays with very high-energies (E = 10^{16} - 10^{18} eV). The instrument is composed of different particle detector systems suitable for the detailed study of the properties of Extensive Air Showers (EAS) developed by cosmic rays in the atmosphere. Among the EAS observables studied with the detector, the charged number of particles, the muon content (at different energy thresholds), and the number of electrons are found. By comparing the measurements of these air-shower parameters with the expectations from MC simulations, different hadronic interaction models can be tested at the high-energy regime with the KASCADE-Grande experiment. In this work, the results of a study on the evolution of the muon content of EAS with zenith angle, performed with the KASCADE-Grande instrument, is presented. Measurements are compared with predictions from MC simulations based on the QGSJET II, QGSJET II-04, SIBYLL 2.1 and EPOS 1.99 hadronic interaction models. A mismatch between experiment and simulations is observed. A similar problem is found for the evolution of the lateral distribution function of muons in the atmosphere.

研究动机与目标

• 利用 EAS 中的μ子含量测试主流强子相互作用模型（QGSJET II、QGSJET II-04、SIBYLL 2.1 和 EPOS 1.99）在极高能量（10^16–10^18 eV）下的表现。
• 研究μ子横向密度分布与吸收系数随天顶角的变化，以反映大气深度的影响。
• 识别实验测量与蒙特卡罗模拟之间的系统性偏差，特别是μ子含量与横向分布函数方面。
• 评估组成假设（混合 vs. 纯物质）对数据与模型预测之间偏差的影响。

提出的方法

• 从 KASCADE-Grande 的 1434 天观测数据中选取事件，采用严格筛选条件：触发站数 >11，簇射核心位于 8×10^4 m² 内，核心距离 270–440 m，天顶角 0°–40°。
• 通过 CIC（核心无关修正）方法对总带电粒子数 N_ch 进行修正，以消除大气衰减影响，使数据在参考天顶角 θ_ref = 22° 下可比。
• 通过将固定径向距离 r_f 处的μ子密度拟合至指数衰减律 ρ_μ(r_f) = ρ^0_μ(r_f) exp[-X_0 sec(θ)/λ_μ(r_f)]，提取吸收系数 λ_μ。
• 针对固定 N^CIC_ch 分箱，在五个天顶角区间内分析μ子横向密度分布 ρ_μ(r)，以研究其随大气深度的演化。
• 使用 CORSIKA 搭载 Fluka（低能段）与四种强子模型（高能段）的蒙特卡罗模拟生成参考预测。
• 分析聚焦于 log N_μ ≈ 4.6–5.2 与 N^CIC_ch ∈ [7.04, 7.28] 范围，对应在混合组成假设下 EAS 能量约为 ~10^16.9–17.1 eV。

实验结果

研究问题

• RQ1 KASCADE-Grande 测量的广延大气簇射中，μ子含量与横向密度分布如何随大气深度演化？
• RQ2 QGSJET II、QGSJET II-04、SIBYLL 2.1 与 EPOS 1.99 模型在不同天顶角下对实测μ子吸收系数的预测与实际观测的符合程度如何？
• RQ3 数据与模拟之间的偏差是否依赖于原初宇宙射线组成的假设（如混合 vs. 纯质子或铁）？
• RQ4 在不同簇射能量与核心距离下，实测μ子密度分布是否落在各模型的预测范围内？
• RQ5 当前强子相互作用模型能否准确描述经 CIC 修正后μ子密度随天顶角的演化行为？

主要发现

• 实测的μ子吸收系数 λ_μ 随天顶角的变化行为与所有四种强子相互作用模型（QGSJET II、QGSJET II-04、SIBYLL 2.1、EPOS 1.99）均存在显著差异，表明存在系统性偏差。
• 偏差在较低簇射能量（N^CIC_ch ≈ 7.04–7.28）最为明显，且随能量升高而减小。
• 在径向距离 r_f 介于 220 m 至 400 m 的范围内，实验μ子密度分布始终低于模型预测值，尤其在大天顶角时更为显著。
• 即使在纯质子或纯铁组成假设下，实验数据仍位于模型预测带之外，表明问题并非源于组成不确定性。
• 结果与μ子修正函数无关，因分析直接基于未经修正的μ子密度进行，确认偏差源于模型预测本身。