[論文レビュー] 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

ひとこと要約

本研究では、Tunka-25 Cherenkovアレイのデータを用いて、10^15から10^17 eVのエネルギー範囲における宇宙線エネルギースペクトルおよび質量成分の改善された解析を提示する。CORSIKAシミュレーションに基づく新しい横方向分布関数（LDF）フィット手法を開発し、X_max（シャワー最大深度）を推定する二つの独立した手法（LDFの傾きとパルスFWHM）を適用することで、複雑な質量成分（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.

研究の動機と目的

• 洗練されたCherenkov光の横方向分布関数（LDF）フィッティングを用いて、宇宙線エネルギーおよびシャワー最大深度の再構築を改善すること。
• Tunka-25アレイを用いて、10^15–10^17 eVの範囲における宇宙線の一次エネルギースペクトルおよび質量成分を評価すること。
• X_maxを推定する二つの独立した手法（LDFの傾きとパルスFWHM）の比較を行うこと。
• 強子相互作用モデルの不確実性がX_maxのエネルギー依存性および質量成分推定に与える影響を評価すること。
• CORSIKAシミュレーションに基づく詳細な「実験モデル」を用いて、実験的再構築精度の妥当性を検証すること。

提案手法

• Cherenkov光プロファイルにより適したフィットを実現するため、可変のべき乗指数bと遷移半径R_knを導入した新しいLDFモデルを開発した。両者とも、傾きパラメータP = Q(100)/Q(200)の関数として定義された。
• CORSIKAシミュレーションを用いて、0°、15°、25°の天の赤道角における陽子および鉄の一次粒子に対して、R_0 = 10^(2.95 - 0.245P) m、R_kn = 155 - 13P m、b = 1.19 + 0.23P のLDFパラメータを導出。
• 実際のアレイ幾何、検出器応答、エネルギースペクトルを再現する10^4～10^5件の人工イベントを生成した「実験モデル」を構築し、再構築性能と誤差をシミュレート。
• QUEST実験によるキャリブレーションを経て、Q_175（175 mでの光フラックス）からE_0 = 400·Q_175^0.95 を用いて一次エネルギーを再構築。
• 二つの独立した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 Tunka-25 Cherenkovアレイが測定した10^15–10^17 eV範囲における一次宇宙線エネルギースペクトルは何か？
• RQ2 LDFの傾きとパルスFWHMという二つの独立した手法は、X_maxの推定においてどのように比較できるか？
• RQ3 観測されたX_max分布を踏まえると、このエネルギー範囲における宇宙線の最も確率の高い一次質量成分は何か？
• RQ4 強子相互作用モデルの不確実性は、X_maxのエネルギー依存性と質量成分推定にどのように影響を与えるか？
• RQ5 Tunka-25アレイにおけるエネルギーおよびX_max再構築の系統的誤差と分解能限界は何か？

主な発見

• 300時間分のデータから微分エネルギースペクトルが再構築され、検出効率100%の閾値は8×10^14 eV、50%は5×10^14 eVであった。
• 3×10^15 eV以上の一次エネルギーに対して、エネルギー分解能は約15%、コア位置分解能は約5 m、X_max分解能は約30 g/cm²であった。
• FWHMに基づくX_max手法はLDF傾き手法に比べて閾値が高かったが、1.5×10^15 eV以上では良好な一致を示した。
• 3–10 PeVにおける実験的X_max分布（5×10^15 eVに正規化）は、複雑な質量成分（p:He:CNO:Fe = 0.25:0.25:0.25:0.25）を仮定したシミュレーションと最も一致した。
• 純粋な陽子または鉄成分では、データからの標準偏差に顕著なずれが生じたが、中間的な混合成分は複雑な質量成分モデルと整合性のある標準偏差を示した。
• X_maxのエネルギー依存性は強子相互作用モデルによって顕著に異なり、QGSJET-01とQGSJET-IIが極端な境界を示した。これは、モデル依存性を解消するための広いエネルギー範囲のカバーが不可欠であることを示唆している。