[論文レビュー] Fluxes of Atmospheric Leptons at 600 GeV - 60 TeV

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

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

ひとこと要約

本稿では、QGSJET相互作用モデルを用いたCORSIKAシミュレーションを用いて、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^*)$ の補正を導入すること。
• ミューオンのエネルギー損失および崩壊効果をスペクトルモデルに組み込み、CORSIKAでシミュレートされたデータへのフィット精度を著しく向上させること。
• 実験的応用に耐える堅牢性を確保するため、異なる大気条件（例：南極の冬期大気）および外部磁場（例：50 nT）の下でモデルを検証すること。

提案手法

• CORSIKA 6.030 を用い、QGSJET高エネルギー相互作用モデルと Hörandel (2003) の宇宙線スペクトルモデルを組み合わせて、大気ミューオンおよびニュートリノフラックスをシミュレートした。
• Volkova (1969) に由来する $\cos(\theta^*)$ の補正を適用し、元の $\cos(\theta)$ 依存性を置き換え、全天頂角にわたりフラックス式の有効性を拡張した。
• Gaisser (1990) の公式の修正版を用い、エネルギー損失および崩壊効果をパラメータ化された補正関数で組み込んだ。
• 統計的収束を確保するため、600 GeV 以上のエネルギーを持つ $10^{10}$ 個のシャワーを用い、約 200 GHz-CPU 日の計算時間を要した。
• 標準米国大気および南極の冬期大気の異なる大気プロファイル、および 50 nT の磁場条件下でモデルを検証し、フィットパラメータへの感度を評価した。
• ミューオン、ミューオンニュートリノ、電子ニュートリノのスペクトルを、エネルギー依存の正規化 $A$、スペクトル指数 $\gamma$、および $\cos(\theta^*)$ を用いた角度依存性を持つ解析的形でフィットした。

実験結果

リサーチクエスチョン

• RQ1 標準的な $\cos(\theta)$ モデルと比較して、$\cos(\theta^*)$ の補正は、特に 70° を超える天頂角において、ミューオンフラックスのパラメータ化精度をどの程度向上させるか？
• 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$ のフィットパラメータにはほとんど影響を及ぼさず、実際の検出器環境におけるモデルの堅牢性を確認した。