고에너지 우주선이 만든 대기·지표 중성자 플럭스 모델링 및 활용 가이드

초록

**
본 논문은 1 GeV부터 1 PeV까지의 고에너지 우주선(HECR)이 대기와 지표면에 생성하는 중성자 플럭스를 몬테카를로 시뮬레이션으로 정량화하고, 다양한 고도에서의 결과를 lookup table 형태로 제공한다. 이를 통해 임의의 우주선 스펙트럼에 대해 중성자 강도를 쉽게 계산할 수 있으며, 지표면 중성자 증강만으로는 원초 입자의 에너지 분포를 역추정하기 어렵다는 결론을 제시한다.

**

상세 분석

**
이 연구는 현재 지구에 도달하는 우주선 스펙트럼이 1 GeV 이하의 저에너지 입자에 의해 지배된다는 점을 출발점으로 삼는다. 그러나 태양 플레어, 초신성, 감마선 폭발 등 급격한 천체 현상은 1 GeV를 훨씬 초과하는 고에너지 입자를 방출하며, 이러한 고에너지 우주선(HECR)이 지구 대기에 충돌하면 광범위한 대기 샤워가 발생한다. 저에너지 우주선이 주로 생성하는 전리층 이온화와 달리, HECR은 깊은 대기층까지 침투해 2차 입자인 뮤온과 중성자를 대량 생산한다. 특히 중성자는 전기적으로 중성이라 대기와 물질에 의해 쉽게 흡수되지 않아, 고고도 비행기 승무원뿐 아니라 지표면까지 도달해 인체에 직접적인 방사선 피폭을 유발한다.

연구팀은 CORSIKA와 GEANT4 기반의 몬테카를로 시뮬레이션 파이프라인을 구축해, 1 GeV, 10 GeV, 100 GeV, 1 TeV, 10 TeV, 100 TeV, 1 PeV 등 7개의 대표적인 기본 입자 에너지에 대해 대기 전단을 모델링하였다. 각 시뮬레이션은 10⁶개의 1차 입자를 투입해 통계적 신뢰도를 확보했으며, 고도별(0 km, 5 km, 10 km, 15 km, 20 km) 중성자 에너지 스펙트럼을 추출했다. 결과는 중성자 플럭스가 기본 입자 에너지와 고도에 따라 비선형적으로 증가함을 보여준다. 예를 들어, 1 PeV 입자는 지표면에서 10⁻⁴ cm⁻²·s⁻¹ 수준의 중성자 플럭스를 생성하는 반면, 1 GeV 입자는 거의 검출되지 않는다. 또한, 고에너지 입자는 중성자 에너지 분포가 넓어 1 MeV 이하부터 수백 MeV까지 포괄한다는 점이 강조된다.

핵심적인 실용적 산출물은 “lookup table”이다. 표는 각 기본 입자 에너지와 고도에 대해 단위 입자당 생성되는 중성자 수(에너지 구간별 포함)를 제공한다. 사용자는 관측된 혹은 가정된 우주선 스펙트럼을 이 표와 컨볼루션(convolution)함으로써, 특정 시점·위치에서 기대되는 중성자 플럭스를 즉시 계산할 수 있다. 이는 과거 대기 방사능 기록(예: 빙핵, 나무 고리)이나 현대 항공·우주 방사선 보호 설계에 직접 적용 가능하다.

하지만 논문은 중요한 역학적 한계도 제시한다. 지표면에서 측정되는 중성자 증강은 기본 입자 에너지 스펙트럼의 형태에 민감하지 않다. 즉, 높은 에너지의 HECR이 소량이라도 존재하면 중성자 플럭스는 급격히 상승하지만, 그 상승 폭만으로는 원천 입자의 에너지 분포를 역추정하기 어렵다. 이는 “역문제(inverse problem)”가 수학적으로 ill‑posed임을 의미한다. 따라서 중성자 관측만으로는 특정 천체 사건(예: 초신성 폭발)의 특성을 정확히 규명하기 힘들며, 뮤온, 감마선, 전리층 전도도 등 다중 관측 지표와 결합해야 한다는 결론에 이른다.

**