[论文解读] The energy spectrum of all-particle cosmic rays around the knee region observed with the Tibet-III air-shower array

M. Amenomori|Sakura: Saitama United Cyber Repository of Academic Resources (Saitama University)|Mar 6, 2008

Astrophysics and Cosmic Phenomena被引用 4

一句话总结

本研究基于西藏-III大气簇射阵列4.5年的数据，给出了10^14至10^17 eV能量区间内全粒子宇宙射线能谱的更新结果，确认了约4 PeV处的能谱膝点。利用QGSJET+HD相互作用模型和改进的蒙特卡罗模拟，测得1 PeV以下的谱指数为-2.68 ± 0.02，4 PeV以上的谱指数为-3.12 ± 0.01，表明在膝点区域能谱的陡降出现显著硬化。

ABSTRACT

We have already reported the first result on the all-particle spectrum around the knee region based on data from 2000 November to 2001 October observed by the Tibet-III air-shower array. In this paper, we present an updated result using data set collected in the period from 2000 November through 2004 October in a wide range over 3 decades between $10^{14}$ eV and $10^{17}$ eV, in which the position of the knee is clearly seen at around 4 PeV. The spectral index is -2.68 $\pm$ 0.02(stat.) below 1PeV, while it is -3.12 $\pm$ 0.01(stat.) above 4 PeV in the case of QGSJET+HD model, and various systematic errors are under study now.

研究动机与目标

• 利用扩展数据集，以高统计精度测量膝点区域的全粒子宇宙射线能谱。
• 改进西藏-III大气簇射阵列的能量标定和簇射大小估计，以实现主粒子能量的精确重建。
• 利用更新的蒙特卡罗模拟和QGSJET+HD相互作用模型，研究约4 PeV附近的谱硬化与膝点结构。
• 评估10^14–10^17 eV能量范围内，相互作用模型、初级成分假设和探测器响应带来的系统误差。
• 通过引入簇射年龄参数s，改进簇射粒子的横向分布函数，提升簇射大小估计的准确性。

提出的方法

• 西藏-III大气簇射阵列包含533个闪烁探测器和0.5 cm厚的铅板，可在10^14 eV以上实现子-0.2°的角分辨率，用于测量簇射大小和到达方向。
• 通过拟合横向粒子密度分布，采用改进的NKG函数估计簇射大小（Ne），其中参数a(s)和b(s)由CORSIKA模拟得出。
• 通过蒙特卡罗模拟标定簇射大小与主粒子能量的关系，得到E0 = 1.88 × (Ne/1000)^0.92 [TeV]（当secθ ≤ 1.1时）。
• 能量分辨率在200 TeV时为36%，在2000 TeV时为17%，当Ne ≥ 10^5时，簇射大小分辨率可达7%。
• 采用QGSJET+HD相互作用模型进行模拟，谱指数通过10^14–10^17 eV能量区间内观测到的Ne–E0关系推导得出。
• 主成分、相互作用模型和探测器响应带来的系统误差正在积极评估中。

实验结果

研究问题

• RQ1 基于扩展的西藏-III数据，全粒子宇宙射线能谱的膝点精确能量位置是什么？
• RQ2 谱指数在膝点区域如何演化，这对宇宙射线加速或传播机制有何启示？
• RQ3 相互作用模型和探测器响应带来的系统误差在多大程度上影响了推导出的能谱？
• RQ4 利用改进的NKG函数，从横向粒子密度分布重建簇射大小的精度如何？
• RQ5 更新的模拟框架（QGSJET+HD）是否提高了观测谱与模拟谱之间的一致性？

主要发现

• 全粒子宇宙射线能谱在约4 PeV处明显表现出膝点，这由谱斜率的显著变化所证实。
• 在1 PeV以下，谱指数为-2.68 ± 0.02（统计误差），表明膝点前区域谱较软。
• 在4 PeV以上，谱指数陡变为-3.12 ± 0.01（统计误差），证实了谱陡降过程中的显著硬化。
• 新能谱在约10^16 eV处的强度比先前结果低约10%，可能归因于改进的蒙特卡罗模拟和更大的数据量。
• 在2000 TeV处17%的能量分辨率和在Ne ≥ 10^5时7%的簇射大小分辨率，支持在膝点区域实现高精度能量重建。
• 主成分、相互作用模型和探测器效应带来的系统误差目前正处于详细调查中。