[论文解读] Test of hadronic interaction models with data from the Pierre Auger Observatory

Ralph Engel|arXiv (Cornell University)|Jun 13, 2007

Particle physics theoretical and experimental studies被引用 29

一句话总结

本研究利用皮埃尔·阿杰尔天文台的荧光探测器与表面探测器联合数据，测试了超高能空气簇射中强子相互作用模型。通过利用电磁簇射剖面的普遍性，独立推断距簇射核心1000米处的μ子密度，发现与QGSJET II.03模型对10^19 eV质子簇射的预测相比，实际观测到的μ子数量多出1.45 ± 0.11（统计）± 0.09（系统）倍，表明当前模型存在显著偏差。

ABSTRACT

The Pierre Auger Observatory allows the measurement of both longitudinal profiles and lateral particle distributions of high-energy showers. The former trace the overall shower development, mainly of the electromagnetic component close to the core where the latter reflect the particle densities in the tail of the shower far away from the core and are sensitive to both the muonic and electromagnetic components. Combining the two complementary measurements, predictions of air shower simulations are tested. In particular the muon component of the tank signals, which is sensitive to hadronic interactions at high energy, is studied with several independent methods. Implications for the simulation of hadronic interactions at ultra-high energy are discussed.

研究动机与目标

• 利用皮埃尔·阿杰尔天文台的独立测量，测试超高能空气簇射中强子相互作用模型的预测。
• 利用电磁簇射剖面的普遍性，不假设原初宇宙射线成分，确定空气簇射中的μ子含量。
• 通过从簇射最大深度推断μ子密度，独立校准表面探测器的能量标度，而不依赖荧光探测器测量。
• 评估多种分析方法（包括恒定强度截断法与混合事件技术）在μ子多重性预测上的一致性。
• 识别观测到的μ子密度与主流强子相互作用模型（如QGSJET II.03和SIBYLL 2.1）预测之间的差异。

提出的方法

• 利用电磁簇射纵向剖面的普遍性，将距簇射核心1000米处的表面探测器信号参数化，其依赖于能量、距簇射最大深度的距离（DG）和天顶角。
• 对表面探测器数据应用恒定强度截断法，假设宇宙射线到达方向各向同性，通过展平dN/dsin²θ分布来提取相对μ子数N_ν^rel。
• 使用CORSIKA和GEANT4的蒙特卡洛模拟来建模簇射发展与探测器响应，包含X_max的涨落和重建分辨率。
• 分析混合事件（同时被荧光探测器和表面探测器触发）以交叉验证μ子信号预测与实际观测信号，利用重建的X_max和荧光能量。
• 比较多种方法的结果——恒定强度截断法、混合事件分析和倾斜簇射分析——以验证一致性并减小系统不确定性。
• 采用μ子信号的参数化形式：S_MC(E,θ,X_max) = S_em(E,θ,DG) + N_ν^rel × S_ν^QGSII,p(10^19 eV,θ,DG)，其中DG = X_ground - X_max。

实验结果

研究问题

• RQ1 观测到的超高能空气簇射中μ子密度在多大程度上与QGSJET II.03和SIBYLL 2.1强子相互作用模型的预测一致？
• RQ2 能否利用电磁簇射剖面的普遍性，在不假设原初成分的前提下独立推断空气簇射中的μ子含量？
• RQ3 从恒定强度截断法和混合事件分析中得出的μ子多重性结果在多大程度上一致？
• RQ4 当通过X_max测量校准μ子含量时，表面探测器的隐含能量标度是多少？
• RQ5 观测到的μ子过量对原初宇宙射线成分和当前强子相互作用模型有效性的含义是什么？

主要发现

• 从恒定强度截断法得出，在10^19 eV能量下，相对μ子数为1.45 ± 0.11（统计）+0.11/-0.09（系统），表明比QGSJET II.03对质子簇射的预测多出45%的μ子。
• 在θ = 38°时，测得的μ子信号为37.5 ± 1.7（统计）+2.1/-2.3（系统）VEM，对应能量标度比荧光探测器重建结果高约30%（E ≈ 1.3E_FD）。
• 混合事件分析得出在E = 1.3E_FD时，N_μ^rel = 1.53 ± 0.05，与恒定强度截断法结果一致，证实了μ子过量的存在。
• 分析未发现垂直簇射与倾斜簇射（60° < θ < 70°）之间μ子含量存在显著差异，支持该方法的稳健性。
• 观测到的μ子密度与QGSJET II.03和SIBYLL 2.1对铁核初级粒子的预测不一致，若假设模型正确，则暗示原初成分重于铁——与实测X_max值矛盾。
• 数据与模拟之间的偏差在定性上与以往空气簇射实验中报告的不一致性相似，凸显了在超高能区模拟强子相互作用时持续存在的挑战。