[论文解读] Cosmic Ray Energy Spectrum and Mass Composition from 10^15 to 10^17 eV by Data of the Tunka EAS Cherenkov Array

Н. Буднев, Д. В. Чернов+12 more|arXiv (Cornell University)|Nov 8, 2005

Astrophysics and Cosmic Phenomena被引用 6

一句话总结

本研究利用图尔卡-25切伦科夫阵列的数据，对10^15至10^17 eV能量范围内的宇宙射线能谱和质量组成进行了改进分析。通过基于CORSIKA模拟开发的新横向分布函数（LDF）拟合方法，并应用两种独立方法——LDF陡度与脉冲FWHM——估算簇射最大深度（X_max），作者发现复杂质量组成（p:He:CNO:Fe = 0.25:0.25:0.25:0.25）最能解释观测到的X_max分布，能量分辨率约为15%，X_max分辨率在3×10^15 eV以上约为30 g/cm²。

ABSTRACT

We present results of an improved analysis of the experimental data of the EAS Cherenkov array Tunka-25. A new function to fit the Cherenkov light lateral distribution LDF at core distances from 0 to 350 m has been developed on the base of CORSIKA simulations and applied to the analysis of Tunka data. Two methods to estimate the EAS maximum position have been used. The one is based on the pulse FWHM, the other on the light LDF. We present the primary energy spectrum in the energy range 10^15 - 10^17 eV. The use of the depth of the EAS maximum to determine the mean mass composition is discussed.

研究动机与目标

• 通过改进切伦科夫光横向分布函数（LDF）拟合，提升宇宙射线能量与簇射最大深度的重建精度。
• 利用图尔卡-25阵列评估10^15–10^17 eV能量范围内宇宙射线的初级能谱与质量组成。
• 比较两种独立方法——LDF陡度与脉冲FWHM——在估算簇射最大深度（X_max）方面的表现。
• 评估强子相互作用模型对X_max能量依赖性及质量组成推断的影响。
• 通过基于CORSIKA模拟的详细“实验模型”验证实验重建的准确性。

提出的方法

• 开发了一种新型LDF模型，其幂律指数b与过渡半径R_kn均作为陡度参数P = Q(100)/Q(200)的函数，以更精确拟合模拟的切伦科夫光分布。
• 基于质子与铁初级粒子在0°、15°和25°天顶角下的CORSIKA模拟，得出LDF参数：R_0 = 10^(2.95 - 0.245P) m，R_kn = 155 - 13P m，b = 1.19 + 0.23P。
• “实验模型”生成了10^4–10^5个人工事件，包含真实的阵列几何结构、探测器响应与能谱，用于模拟重建性能与误差。
• 利用E_0 = 400·Q_175^0.95从Q_175（175 m处的光通量）重建初级能量，该公式通过QUEST实验校准。
• 应用两种独立的X_max估计方法：(1) 通过LDF陡度P，使用H_max = 17.63 - 0.0786·(P + 8.916)^2 km；(2) 通过脉冲FWHM，使用ΔX = 1677 + 1006·log10(FWHM) g/cm²（在250 m处）。
• 将模拟的X_max分布通过实验延长率97±3 g/cm²归一化至固定能量（5×10^15 eV），以与数据进行比较。

实验结果

研究问题

• RQ1 图尔卡-25切伦科夫阵列测量的10^15–10^17 eV能量范围内初级宇宙射线能谱为何？
• RQ2 在估算簇射最大深度（X_max）方面，LDF陡度与脉冲FWHM两种独立方法的比较结果如何？
• RQ3 在该能量范围内，基于观测到的X_max分布，宇宙射线最可能的初级质量组成为何？
• RQ4 强子相互作用模型的不确定性如何影响从X_max能量依赖性推断的质量组成？
• RQ5 图尔卡-25阵列在能量与X_max重建中的系统误差与分辨率极限为何？

主要发现

• 从300小时数据中重建了微分能谱，在8×10^14 eV处达到100%探测效率阈值，5×10^14 eV处为50%。
• 初级能量高于3×10^15 eV时，能量分辨率为~15%，核心位置分辨率为~5 m，X_max分辨率为~30 g/cm²。
• 基于FWHM的X_max方法阈值高于LDF陡度方法，但两者在1.5×10^15 eV以上表现出良好一致性。
• 在3–10 PeV能量范围内，归一化至5×10^15 eV的实验X_max分布与假设复杂质量组成（p:He:CNO:Fe = 0.25:0.25:0.25:0.25）的模拟结果最佳吻合。
• 纯质子或纯铁组成在标准差上与数据存在显著偏差，而中间混合组成则与复杂组成模型一致。
• X_max能量依赖性在不同强子相互作用模型下差异显著，QGSJET-01与QGSJET-II分别给出极端边界，凸显了更广能量覆盖对解决模型依赖性问题的必要性。