초고에너지 핵 전파 시뮬레이션을 위한 공개 프레임워크 CRPropa 2.0

초록

CRPropa 2.0은 Z < 26인 초고에너지 핵과 그 2차 입자인 감마선·중성미자를, 구조화된 은하간 자기장과 우주배경광자장을 고려해 6 × 10¹⁶ eV ~ A × 10²² eV 구간에서 전파시키는 공개 소프트웨어이다. 핵의 광자 상호작용(광핵분해·광양성자·쌍생성)과 붕괴를 모두 포함하고, 1차와 2차 입자의 다중 메신저 연결을 연구할 수 있다. 모듈식 설계와 샘플 애플리케이션을 통해 GZK 한계, 원소 조성, 2차 감마·중성미자 스펙트럼 등을 손쉽게 탐색한다.

상세 분석

CRPropa 2.0은 초고에너지(UHE) 핵 전파 연구에 필요한 물리와 수치 기법을 하나의 통합 프레임워크에 구현한 최초의 공개 코드 중 하나이다. 핵종 Z < 26, 즉 철까지를 다루며, 핵 질량 A에 따라 에너지 상한을 A × 10²² eV 로 설정함으로써 실제 가속기 모델과 일치하도록 설계되었다. 전파는 1차원(거리만 고려)과 3차원(공간적 자기장 구조 포함) 두 모드로 제공되며, 3D 모드에서는 사용자가 직접 입력한 MHD 시뮬레이션 결과나 파워 스펙트럼을 갖는 난류 자기장을 적용할 수 있다.

광자 상호작용은 네 가지 주요 프로세스로 구현된다. 첫째, 광핵분해(photo‑disintegration)는 TALYS 기반의 최신 파라미터를 사용해 핵의 방출 채널(중성자·양성자·알파 입자 등)을 확률적으로 샘플링한다. 둘째, 광양성자 생산(photo‑pion production)은 SOPHIA와 연동해 고에너지 광자와 핵/양성자 간의 파이온 생성 과정을 정확히 계산한다. 셋째, 전자쌍생성(pair production)은 연속 손실 형태로 구현되어 에너지 감쇠를 빠르게 처리한다. 넷째, 핵 붕괴와 방사성 동위원소의 반감기 역시 시간 의존적으로 추적한다.

코스믹 마이크로파 배경(CMB)뿐 아니라 적외·가시·자외선 배경(EBL) 모델을 다중 선택 가능하게 함으로써, 다양한 우주론적 진화 시나리오를 실험할 수 있다. 또한, 소스 분포, 적분 광도, 스펙트럼 지수, 최대 가속 에너지 등 사용자 정의 파라미터를 입력해 복합적인 천체 물리 모델을 구성한다. 출력은 입자별 궤적, 도착 에너지·방향·시간, 2차 감마·중성미자 스펙트럼 등을 포함하는 HDF5 포맷으로 제공돼, 후처리와 시각화가 용이하다.

샘플 애플리케이션에서는 (1) GZK 한계 거리와 핵 조성 변화, (2) 은하간 자기장의 강도와 구조가 도착 스펙트럼에 미치는 영향, (3) 2차 감마·중성미자 플럭스와 원천 핵 종류 간의 상관관계, (4) 전천구적 관측기(예: Auger, TA)와의 비교를 통한 모델 검증 등을 시연한다. 이러한 사례는 CRPropa 2.0이 다중 메신저 천문학에서 핵·광·중성미자 연계 연구를 수행하기 위한 핵심 도구임을 입증한다.

전반적으로 CRPropa 2.0은 물리적 정확도와 계산 효율성을 동시에 달성했으며, 모듈식 설계와 풍부한 문서·예제 코드를 통해 초고에너지 천체물리학 커뮤니티가 공동으로 코드 개선과 검증에 참여할 수 있는 기반을 제공한다.