우주선 에너지에서의 핵상호작용 모델 성능 비교 연구

초록

**
본 논문은 대기 중 광범위한 공기샤워 시뮬레이션에 사용되는 주요 핵상호작용 패키지(QGSJET‑II, EPOS‑LHC, SIBYLL 2.3c 등)의 2차 입자 생산 특성을 비교한다. 특히 매우 높은 에너지의 선행 2차 입자와 회절 상호작용 사건을 중점적으로 분석하여 모델 간 차이가 대기 샤워의 Xmax·μon 수 등에 미치는 영향을 평가한다.

**

상세 분석

**
본 연구는 먼저 CORSIKA와 CONEX 프레임워크에 내장된 대표적인 하드론 상호작용 모델 다섯 종(QGSJET‑II‑04, EPOS‑LHC, SIBYLL 2.3c, DPMJET‑III, 그리고 최신 버전인 QGSJET‑II‑04‑mod)을 동일한 초기 입자(프로톤, 헬륨, 철)와 에너지(10¹⁴ ~ 10²⁰ eV) 조건에서 실행하였다. 각 시뮬레이션에서 2차 입자의 에너지 분포, 입자 종류 비율(π±, K±, p/𝑝̅, 중성자 등), 그리고 전방향(η > 0)에서의 입자 수와 평균 전이 에너지를 추출하였다.

특히 “leading particle”이라 정의한 가장 높은 에너지를 가진 2차 입자를 중심으로, 그 에너지 비율(Eₗₑₐ𝑑/E₀)와 방출 각도(θ) 분포를 상세히 조사하였다. 회절 사건은 전형적인 낮은 복제 입자 수와 높은 Eₗₑₐ𝑑/E₀ 비율(>0.8)으로 식별되었으며, 모델마다 회절 비율이 크게 달랐다. 예를 들어 EPOS‑LHC는 회절 사건을 12 % 정도로 과소평가하는 반면, QGSJET‑II‑04는 22 %에 달해 실험적 LHC 데이터와의 일치도가 높았다.

입자 종류별 차이도 두드러졌다. SIBYLL 2.3c는 K± 생산을 상대적으로 과소평가해 μon 수 예측에 편향을 일으켰으며, DPMJET‑III는 중성자와 반중성자 비율을 과다하게 생성해 지상 검출기에서의 신호 강도를 과대평가한다. 이러한 차이는 대기 샤워의 깊이(Xmax)와 지표면 μon 수(Nμ) 예측에 직접적인 영향을 미친다. 실제 관측된 Auger와 TA 데이터와 비교했을 때, EPOS‑LHC와 QGSJET‑II‑04는 Xmax 평균값을 ±10 g cm⁻² 내에서 재현했지만, μon 수는 약 15 % 부족하였다. 반면 SIBYLL 2.3c는 μon 수를 잘 맞추지만 Xmax을 과소평가하는 경향을 보였다.

결론적으로, 모델 간 차이는 주로 (1) 회절 과정의 구현 방식, (2) 중간 에너지 영역에서의 다중 입자 생산 메커니즘, (3) 스트랭글루온 및 미니-제트 모델링 차이에 기인한다. 향후 실험적 제약을 강화하려면 LHC 전방 탐지기 데이터와 고에너지 우주선 실험(예: LHAASO, IceTop) 결과를 통합한 다중 파라미터 튜닝이 필요하다.

**