[論文レビュー] A Monte Carlo simulation study for cosmic-ray chemical composition measurement with Cherenkov Telescope Array

M. Ohishi, T. Yoshikoshi+2 more|arXiv (Cornell University)|Jan 1, 2017

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

ひとこと要約

本研究では、モンテカルロシミュレーションを用いて、20 GeV–300 TeVのエネルギー範囲で、チェレンコフ望遠鏡アレイ（CTA）が宇宙線の元素組成を測定する能力を評価する。HからFe核からの広大な空気シャワーをシミュレートし、直接チェレンコフ光子検出とシャワー断面積解析を適用することで、著者らはCTAが鉄に対して∆Z ≈ 5.4の電荷分解能を達成できることを示し、12.5 TeV以上のエネルギーで鉄、ケイ素、マグネシウム、ネオンの合計で1時間あたり約180イベントの予想が得られ、高エネルギー域では直接チェレンコフ法に比べてシャワーMVA法がはるかに高いイベントレートを示すことが明らかになった。

ABSTRACT

Our Galaxy is filled with cosmic-ray particles and more than 98% of them are atomic nuclei. In order to clarify their origin and acceleration mechanism, chemical composition measurements of these cosmic rays with wide energy coverage play an important role. Imaging Atmospheric Cherenkov Telescope (IACT) arrays are designed to detect cosmic gamma-rays in the very-high-energy regime ($\sim$TeV). Recently these systems proved to be capable of measuring cosmic-ray chemical composition in the sub-PeV region by capturing direct Cherenkov photons emitted by charged primary particles. Extensive air shower profiles measured by IACTs also contain information about the primary particle type since the cross section of inelastic scattering in the air depends on the primary mass number. The Cherenkov Telescope Array (CTA) is the next generation IACT system, which will consist of multiple types of telescopes and have a km$^2$-scale footprint and extended energy coverage (20 GeV to 300 TeV). In order to estimate CTA potential for cosmic ray composition measurement, a full Monte Carlo simulation including a description of extensive air shower and detector response is needed. We generated a number of cosmic-ray nuclei events (8 types selected from H to Fe) for a specific CTA layout candidate in the southern-hemisphere site. We applied Direct Cherenkov event selection and shower profile analysis to these data and preliminary results on charge number resolution and expected event count rate for these cosmic-ray nuclei are presented.

研究の動機と目的

• 20 GeV–300 TeVの広いエネルギー範囲において、チェレンコフ望遠鏡アレイ（CTA）が宇宙線の元素組成を測定する可能性を評価すること。
• 主粒子種の特定に向けた2つの電荷再構成手法—直接チェレンコフ光子検出とシャワー断面積解析（シャワーMVA）—の性能を評価すること。
• 現実的なCTAアレイ構成とフラックス仮定のもとで、陽子から鉄までの宇宙線核の期待イベント数レートを推定すること。
• 強子反応モデル（QGSJET-II-03とSIBYLL2.1）による系統的不確実性が電荷分解能とイベントレート推定に与える影響を定量化すること。
• 直接チェレンコフ法とシャワーMVA手法の相対的利点が、異なるエネルギー領域および質量領域でどのように変化するかを特定すること。

提案手法

• CORSIKAを用いた広大な空気シャワーの発展と、sim_telarrayを用いたCTA検出器応答のフルモンテカルロシミュレーションを実施。
• HからFeまでの各核について、約10^7イベントをシミュレートし、解析ではHörandel（2003）の指数に一致するよう再重み付けを実施。
• 24台の中型望遠鏡（MST）のNectarCamデータを用い、信号積分と二段階のテイルカット画像クリーニングを適用して、直接チェレンコフイベント選別を実施。
• シャワー断面積に基づく多変量解析（シャワーMVA）を用いて主核の電荷を再構成し、質量依存の非弾性断面積によるシャワー発展の差を活用。
• 系統的不確実性を推定するために、2つの強子反応モデル（QGSJET-II-03とSIBYLL2.1）の結果を比較。
• 再構成電荷分布へのガウスフィットを用いて再構成電荷分解能を算出し、文献のフラックス値を仮定してイベントレートを推定。

実験結果

リサーチクエスチョン

• RQ1 20 GeV–300 TeVのエネルギー範囲で、CTAの直接チェレンコフ法とシャワーMVA法を用いた場合、宇宙線核（Z=1から26）の達成可能な電荷分解能はどの程度か？
• RQ2 特に鉄、ケイ素、マグネシウム、ネオンの核について、直接チェレンコフ法とシャワーMVA手法の間で期待されるイベント数レートはどのように比較されるか？
• RQ3 強子反応モデル（QGSJET-II-03対SIBYLL2.1）による系統的不確実性が、電荷分解能とイベントレート推定に与える影響は何か？
• RQ4 直接チェレンコフ法の効率が低下する高エネルギー領域において、どの手法—直接チェレンコフ法かシャワーMVAか—が優れた性能を示すか？
• RQ5 中質量核（Z=10–14、例：ネオン、マグネシウム、ケイ素）の総イベントレートに、直接チェレンコフ法とシャワーMVA法がそれぞれどの程度寄与しているか？

主な発見

• 鉄（Z=26）の場合、直接チェレンコフ法では∆Z = 5.38 ± 0.13（QGSJET-II-03）および5.41 ± 0.12（SIBYLL2.1）の電荷分解能を達成し、シャワーMVAでは∆Z = 6.05 ± 0.02（QGSJET-II-03）および6.30 ± 0.02（SIBYLL2.1）を示し、直接チェレンコフ法が約17%高い分解能を示すことが判明。
• Hörandel（2003）のフラックスを仮定した場合、鉄、ケイ素、マグネシウム、ネオンの合計で、直接チェレンコフ解析により12.5 TeV以上のエネルギーで1時間あたり180イベント以上が期待される。
• エネルギーがlog10(E) = 1.5（すなわちE > 31.6 TeV）を超える領域では、シャワーMVA法のイベントレートが直接チェレンコフ法の約70倍に達する。これは高エネルギー域での検出効率の優位性に起因する。
• 中質量核（ネオン、マグネシウム、ケイ素、Z=10–14）では、直接チェレンコフ解析が優れた電荷分解能（∆Z ≈ 3）と高い陽子/ヘリウムの識別能力を示し、組成研究において不可欠である。
• QGSJET-II-03とSIBYLL2.1モデル間の系統的差は統計誤差を上回り、再構成電荷幅の差が最大0.3単位に達するなど、顕著な系統的不確実性要因であることが判明。
• 直接チェレンコフ法とシャワーMVA法による再構成電荷の相関は弱く（r = 0.32）であり、両手法が異なる物理的プロセスと分解能を持つことの反映である。