KASCADE Grande가 밝힌 10의16승에서 10의18승 전자볼트 우주선 스펙트럼과 구성

초록

KASCADE‑Grande 실험은 10¹⁶ eV부터 10¹⁸ eV까지의 광범위한 에너지 구간에서 대기권 내 에어샤워를 측정해 전체 입자 스펙트럼과 원소 조성을 재구성하였다. 다중 검출기 시스템을 활용해 전하 입자, 전자, 뮤온을 동시에 관측하고, QGSJet‑II와 EPOS 등 최신 강입자 상호작용 모델을 적용해 분석하였다. 결과는 8 × 10¹⁶ eV 부근에서 무거운 원소 흐름 감소에 기인한 ‘무게‑켈프(knee)’ 구조가 존재함을 보여준다.

상세 분석

KASCADE‑Grande는 기존 KASCADE 배열에 700 × 700 m² 면적을 갖는 Grande 어레이를 추가해 총 0.5 km²의 관측 면적을 확보하였다. 이 다중 검출기 시스템은 비가시광학식 섬광 검출기(전하 입자), 전자 전용 계측기, 그리고 얕은 및 깊은 뮤온 검출기로 구성되어, 각각의 샤워 성분을 독립적으로 측정한다. 핵심 관측 파라미터는 전하 입자 수(N_ch)와 뮤온 수(N_μ)이며, 이 두 값의 상관관계는 원초 입자의 질량에 강하게 의존한다.

데이터 처리 과정에서 핵심은 ‘전하‑뮤온 비율(N_μ/N_ch)’과 ‘Y_CIC = log N_μ / log N_ch’와 같은 질량 민감 파라미터를 정의하고, 이를 Monte‑Carlo 시뮬레이션(QGSJet‑II, FLUKA, EPOS)과 비교해 원소 그룹별 스펙트럼을 추출하는 것이다. 시뮬레이션은 CORSIKA 프레임워크를 이용해 다양한 원소(프로톤, 헬륨, 실리콘, 철 등)와 에너지 범위에 대해 에어샤워 전개를 모델링한다.

전체 입자 스펙트럼 재구성 결과, 10¹⁶ eV에서 10¹⁸ eV 사이에 전반적으로 파워‑로우 형태를 보이지만, 8 × 10¹⁶ eV 부근에서 스펙트럼 지수가 γ = −2.95에서 γ = −3.24로 변하는 ‘스몰 브레이크’가 관측되었다. 통계적 유의성은 2σ 이상이며, 이는 강입자 상호작용 모델을 바꾸어도 유지되는 구조이다.

구성 분석에서는 네 가지 독립적인 방법(뮤온‑전하 비율, Y_CUT, 전자‑풍부/전자‑가난 샘플 분리, 그리고 직접 뮤온 측정)을 적용했으며, 모두 무거운 원소(특히 철군)의 플럭스 감소가 전체 스펙트럼 브레이크의 원인임을 일관되게 보여준다. 가벼운 원소(프로톤·헬륨)는 10¹⁷ eV 이상에서도 상대적으로 평탄한 스펙트럼을 유지한다.

시스템atic 불확실성은 주로 검출 효율, 대기 감쇠 보정, 그리고 강입자 모델 차이에 기인한다. 전체 플럭스 불확실성은 약 10‑15 % 수준이며, 질량 비율에 대한 불확실성은 모델에 따라 20‑30 %까지 변동한다. 그럼에도 불구하고, 무거운 원소의 급격한 감소라는 물리적 신호는 모델 의존성을 넘어서는 강인한 결과로 판단된다.

이러한 결과는 ‘켈프(knee)’가 원소별 강성(rigidity) 의존성을 띠며, 무거운 원소의 켈프가 10¹⁶ eV 이하에서 시작해 10¹⁷ eV에 도달함을 시사한다. 따라서 은하계 내 가속 메커니즘이 한계에 다다를 때, 무거운 원소가 먼저 탈출하고 가벼운 원소가 남는 구조적 전이를 관측한 것으로 해석할 수 있다.