[论文解读] About EAS size spectra and primary energy spectra in the knee region

S. Ter–Antonyan, L. S. Haroyan|ArXiv.org|Mar 6, 2000

High-Energy Particle Collisions Research参考文献 1被引用 10

一句话总结

本论文利用KASCADE、AKENO、EAS-TOP和ANI实验的统一EAS大小谱，在10–100 PeV的膝区重建了原初宇宙射线的能量谱和元素组成。通过使用QGSJET和SIBYLL相互作用模型求解逆问题，识别出与刚度相关的谱指数变陡谱，并发现一个明显的质子/中子组分，其微分能量谱为(6.1±0.7)×10⁻¹¹ (E/Eₖ)⁻¹.⁵ (m²·s·sr·TeV)⁻¹，膝能Eₖ = 2030±130 TeV之前，谱指数γ₁ = -1.5±0.6；膝能之后，谱指数γ₂ = -3.1±0.05。

ABSTRACT

Based on the unified analyses of KASCADE, AKENO, EAS-TOP and ANI EAS size spectra, the approximations of energy spectra of different primary nuclei have been found. Calculations were carried out using SIBYLL and QGSJET interaction models in 0.1-100 PeV primary energy range. The results point to existence of both rigidity-dependent steepening energy spectra at R=200-400 TV and an additional proton (neutron) component with differential energy spectrum (6.1+/-0.7)10^{-11}(E/Ek)^{-1.5} (m^2.s.sr.TeV)^{-1} before the knee Ek=2030+/-130 TeV and with power index g2=-3.1+/-0.05 after the knee.

研究动机与目标

• 通过整合不同天顶角和大气深度的数据，统一KASCADE、AKENO、EAS-TOP和ANI实验的EAS大小谱，以解决这些实验之间EAS大小谱的差异。
• 通过测试与刚度相关的谱指数变陡原初能量谱和额外轻组分，探究宇宙射线谱中膝结构的起源。
• 使用逆问题技术，确定0.1–100 PeV能量范围内原初宇宙射线的元素组成和能量谱。
• 通过与多实验EAS数据对比，验证QGSJET和SIBYLL相互作用模型在膝区的适用性。
• 量化在膝能处额外质子（或中子）组分对全粒子通量的贡献。

提出的方法

• 建立积分方程(1)的逆问题框架，将原初能量谱与观测到的EAS大小谱关联，考虑角度和能量相关的探测效率。
• 采用简化假设：探测概率恒定（Pθ ≡ 1），测量误差服从对数正态分布，通过标定因子η对EAS大小进行标准化。
• 将逆问题转化为χ²最小化框架（方程4），以拟合原初核素及额外组分的谱参数。
• 使用QGSJET和SIBYLL模型模拟EAS发展过程，并计算不同原初核素和能量下的核函数Wθ。
• 对四个实验在多个天顶角和大小区间上的EAS大小谱进行联合拟合，以约束原初谱。
• 测试两种模型：与刚度相关的谱指数变陡谱和额外轻组分，后者组成与谱形由最佳拟合确定。

实验结果

研究问题

• RQ1 当使用多实验统一的EAS数据时，膝区（10–100 PeV）原初宇宙射线能量谱的形状是什么？
• RQ2 膝结构是否源于原初谱的刚度相关谱指数变陡，还是需要引入额外的轻组分？
• RQ3 在膝区，额外质子（或中子）组分的谱指数和归一化值是多少？
• RQ4 QGSJET和SIBYLL相互作用模型在不同实验和天顶角下对观测到的EAS大小谱的描述精度如何？
• RQ5 在膝能处，额外组分对全粒子通量的贡献是多少？

主要发现

• 全粒子谱的膝能Eₖ = 2030±130 TeV，在多个实验和模型中保持一致。
• 在膝能之前，必须存在一个具有微分能量谱(6.1±0.7)×10⁻¹¹ (E/Eₖ)⁻¹.⁵ (m²·s·sr·TeV)⁻¹的明显质子（或中子）组分，谱指数γ₁ = -1.5±0.6。
• 在膝能之后，该额外组分具有硬谱指数γ₂ = -3.1±0.05，表明其对膝后通量有显著贡献。
• 在膝能处（2030 TeV），额外质子（中子）组分对全粒子通量的贡献为20±5%。
• 在R ≃ 200–400 TV的刚度范围内，原初谱的刚度相关谱指数变陡得到支持，表明宇宙射线的注入或传播过程存在转变。
• QGSJET相互作用模型在0–37°天顶角范围内，对KASCADE和ANI实验的EAS大小谱描述精度达到2–5%。