[논문 리뷰] A Monte Carlo simulation study for cosmic-ray chemical composition measurement with Cherenkov Telescope Array

M. Ohishi, T. Yoshikoshi외 2명|arXiv (Cornell University)|2017. 01. 01.

Astrophysics and Cosmic Phenomena참고 문헌 1인용 수 2

한 줄 요약

이 연구는 20 GeV–300 TeV 범위에서 우주선의 화학적 조성을 측정할 수 있는 케렌코프 망원경 어레이(CrTA)의 능력을 평가하기 위해 몬테카를로 시뮬레이션을 사용한다. 수소에서 철까지의 핵에서 발생하는 광범위한 공기 샤워를 시뮬레이션하고 직접 케렌코프 광자 탐지 및 샤워 프로파일 분석을 적용함으로써, 저자들은 CTA가 철에 대해 ∆Z ≈ 5.4의 전하 해상도를 달성할 수 있음을 보여주며, 12.5 TeV 이상에서 철, 규소, 마그네슘, neon의 경우 시계열당 약 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의 넓은 에너지 범위에서 케렌코프 망원경 어레이(CrTA)가 우주선의 화학적 조성을 측정할 잠재력을 평가한다.
• 주입 입자 종류를 식별하기 위한 두 가지 전하 재구성 방법—직접 케렌코프 광자 탐지 및 샤워 프로파일 분석(샤워 MVA)—의 성능을 평가한다.
• 실제 CTA 어레이 구성과 통계적 복사율 가정 하에 수소에서 철까지의 우주선 핵에 대한 기대 사건 수를 추정한다.
• 고에너지 상호작용 모델(QGSJET-II-03 및 SIBYLL2.1)의 체계적 오차가 전하 해상도와 사건 레이트 추정에 미치는 영향을 정량화한다.
• 직접 케렌코프와 샤워 MVA 기술의 상대적 이점이 다양한 에너지 및 질량 영역에서 어떻게 나타나는지 규명한다.

제안 방법

• CORSIKA를 사용해 광범위한 공기 샤워의 발달을, sim_telarray를 사용해 CTA 검출기 반응을 전면 몬테카를로 시뮬레이션 수행.
• 스펙트럼 지수 −2.0로 H에서 Fe까지의 핵 각각에 대해 약 10^7건의 사건을 시뮬레이션하고 분석에서 Hörandel(2003)의 지수와 일치하도록 재가중.
• 24개의 중형 망원경(MST)에서의 NectarCam 자료를 바탕으로 신호 통합 및 이중 수준 꼬리 컷 이미지 정제를 통해 직접 케렌코프 사건 선별 적용.
• 핵 질량에 따른 비탄성 단면적의 차이를 활용해 샤워 발달의 차이를 이용해 주입 전하를 재구성하는 샤워 프로파일 기반 다변량 분석(샤워 MVA) 사용.
• 체계적 오차를 평가하기 위해 두 개의 고에너지 상호작용 모델(QGSJET-II-03 및 SIBYLL2.1) 간 결과 비교.
• 재구성된 전하 분포에 대해 정규분포 피팅을 통해 재구성 전하 해상도 계산 및 문헌 복사율 값 가정 하에 사건 레이트 추정.

실험 결과

연구 질문

• RQ1 20 GeV–300 TeV 범위에서 CTA의 직접 케렌코프 및 샤워 MVA 방법을 사용해 우주선 핵의 전하 해상도(Z=1에서 26까지)를 얼마나 달성할 수 있는가?
• RQ2 특히 철, 규소, 마그네슘, neon에 대해 직접 케렌코프와 샤워 MVA 기술 간 기대 사건 수 레이트는 어떻게 비교되는가?
• RQ3 고에너지 상호작용 모델(QGSJET-II-03 대비 SIBYLL2.1)의 체계적 오차가 전하 해상도와 사건 레이트 추정에 어떤 영향을 미치는가?
• RQ4 직접 케렌코프 탐지 효율이 감소하는 고에너지 영역에서 직접 케렌코프와 샤워 MVA 중 어느 방법이 더 뛰어난 성능을 보이는가?
• RQ5 중간 질량 핵(예: Ne, Mg, Si, Z=10–14)에 대해 직접 케렌코프와 샤워 MVA가 총 사건 레이트에 기여하는 상대적 기여도는 어떠한가?

주요 결과

• 철(Z=26)에 대해 직접 케렌코프 방법은 QGSJET-II-03 기준 ∆Z = 5.38 ± 0.13, SIBYLL2.1 기준 ∆Z = 5.41 ± 0.12의 해상도를 달성하며, 샤워 MVA는 QGSJET-II-03 기준 ∆Z = 6.05 ± 0.02, SIBYLL2.1 기준 ∆Z = 6.30 ± 0.02를 기록한다. 이는 직접 케렌코프 방식이 약 17% 향상된 해상도를 제공함을 의미한다.
• Hörandel(2003) 복사율을 가정할 경우, 직접 케렌코프 분석을 통해 철, 규소, 마그네슘, neon의 합산 사건 수가 12.5 TeV 이상에서 시간당 180건 이상을 초과한다.
• 에너지가 log10(E) = 1.5 이상(즉, E > 31.6 TeV)일 경우, 샤워 MVA는 직접 케렌코프 방식 대비 약 70배 높은 사건 레이트를 제공한다. 이는 고에너지에서의 탐지 효율이 더 뛰어나기 때문이다.
• 중간 질량 핵(예: Ne, Mg, Si, Z=10–14)에 대해 직접 케렌코프 분석은 ∆Z ≈ 3의 뛰어난 전하 해상도와 높은 수소/Helium 거부 능력을 제공하므로 조성 연구에 필수적이다.
• QGSJET-II-03과 SIBYLL2.1 모델 간 체계적 차이가 통계 오차를 초월하며, 재구성된 전하 폭의 차이가 최대 0.3 단위까지 발생하여 중요한 체계적 오차 원천임을 나타낸다.
• 직접 케렌코프와 샤워 MVA 방법 간 재구성 전하 간 상관계수가 약 r = 0.32로 낮아, 두 방법 간 물리적 과정과 해상도의 상이함을 반영한다.