[论文解读] Fluxes of Atmospheric Leptons at 600 GeV - 60 TeV

D. Chirkin|ArXiv.org|Jul 7, 2004

Astrophysics and Cosmic Phenomena参考文献 2被引用 17

一句话总结

本论文利用 CORSIKA 模拟结合 QGSJET 相互作用模型，对 600 GeV 至 60 TeV 范围内的大气μ子和中微子通量进行了优化参数化。该工作引入了关键能量标度的修正 $ \cos(\theta^*)$ 依赖关系，考虑了μ子的能量损失与衰变效应，并提供了在所有天顶角下均稳定的、能量相关的谱拟合，其精度显著提升，且经受了模拟数据及包括南极冬季大气和磁场效应在内的大气模型验证。

ABSTRACT

The atmospheric muon flux is generated with CORSIKA 6.030 (using QGSJET high-energy interaction model) and fitted with the formula from Gaisser (1990). Spectra and mass distribution of primaries were taken according to the poli-gonato model of Horandel (2003). A convenient parameterization of the cos(th)->cos(th*) correction, originally tabulated in Volkova (1969), is presented. This correction, together with muon energy losses and decay, is shown to significantly improve the fit to the simulated data.

研究动机与目标

• 改进高能区（600 GeV–60 TeV）大气μ子和中微子通量的参数化，超越标准的 $\cos(\theta)$ 近似。
• 通过引入 $\cos(\theta^*)$ 修正，克服平面大气模型的局限性，以考虑大气曲率及π介子和K介子关键能量的天顶角依赖性。
• 将μ子能量损失与衰变效应纳入谱模型，显著提升对 CORSIKA 模拟数据的拟合精度。
• 在多种大气条件（如南极冬季大气）和外部场（如 50 nT 磁场）下验证模型，确保其在实验应用中的鲁棒性。

提出的方法

• 使用 CORSIKA 6.030 结合 QGSJET 高能相互作用模型与 Hörandel (2003) 的宇宙射线能谱模型，模拟大气μ子和中微子通量。
• 应用源自 Volkova (1969) 的 $\cos(\theta^*)$ 修正，扩展通量公式的适用范围至所有天顶角，取代原始的 $\cos(\theta)$ 依赖关系。
• 采用 Gaisser (1990) 公式的改进版本拟合模拟的μ子能谱，通过参数化校正函数引入能量损失与衰变效应。
• 使用 $10^{10}$ 个能量大于 600 GeV 的簇射样本，耗时约 200 GHz-CPU 天，以确保统计收敛。
• 在不同大气剖面（标准美国大气、南极冬季大气）和磁场条件（50 nT）下验证模型，评估拟合参数的敏感性。
• 采用包含能量依赖性归一化 $A$、谱指数 $\gamma$ 及通过 $\cos(\theta^*)$ 表征角度依赖性的解析形式，对μ子、μ子中微子和电子中微子谱进行拟合。

实验结果

研究问题

• RQ1 与标准 $\cos(\theta)$ 模型相比，$\cos(\theta^*)$ 修正在高天顶角下如何提升μ子通量参数化的精度？
• RQ2 在 600 GeV–60 TeV 范围内，μ子能量损失与衰变效应对大气μ子谱形的影响程度如何？
• RQ3 当大气剖面或磁场发生变化时，拟合参数 $A$ 和 $\gamma$ 的稳定性如何？
• RQ4 与简单幂律假设相比，采用 Hörandel (2003) 的刚度依赖性截止模型对推导出的μ子通量谱指数有何影响？
• RQ5 所推导的参数化在不同中微子味和能量范围内，对模拟谱的重现能力如何？

主要发现

• 源自大气曲率与μ子产生几何结构的 $\cos(\theta^*)$ 修正，显著提升了对 CORSIKA 模拟μ子通量在所有天顶角下的拟合效果，尤其在超过 70° 时表现更优。
• 引入μ子能量损失与衰变效应后，谱指数 $\gamma$ 相较于简单拟合降低了约 0.01，修正模型得到的μ子 $\gamma = 2.715 \pm 0.001$。
• 当大气剖面从标准美国大气切换至南极冬季大气时，拟合参数 $A$ 和 $\gamma$ 保持稳健，仅观察到微小变化。
• μ子中微子通量模型的谱指数为 $\gamma = 2.684 \pm 0.001$，归一化为 $A = 0.646 \pm 0.003$，在不同拟合边界下均表现出高度稳定性。
• 电子中微子谱模型表现出更复杂的能量依赖性，其指数 $\zeta(\theta)$ 满足 $\zeta = a + b\log_{10}E$，其中 $a$ 和 $b$ 依赖于 $\cos(\theta)$，且在 $\cos(\theta_l) \geq 0.3$ 时拟合稳定。
• 引入 50 nT 磁场后，μ子散射宽度增加 2–4 倍，但对拟合参数 $A$ 和 $\gamma$ 的影响可忽略不计，证实了模型在真实探测器环境中的鲁棒性。