지역 감마선 방출량 추정에 숨은 시스템적 오차와 핵강화 인자 재평가

초록

이 논문은 GeV–TeV 범위에서 지역 감마선 방출량을 추정할 때, 원초적 우주선 스펙트럼, 전파 모델, 은하간 헬륨 비율, 그리고 γ‑ray 생성 단면적 등 네 가지 주요 요인이 각각 5%‒50% 수준의 불확실성을 만든다는 점을 보여준다. 특히 최신 프로톤·헬륨 측정에 기반하면 기존에 사용되던 일정한 핵강화 인자(𝜖_M)는 에너지에 따라 크게 변동하므로 실용적이지 않다.

상세 분석

본 연구는 Fermi‑LAT이 제공한 지역 감마선 방출량 측정값과 이론적 예측 사이의 차이를 정량화하기 위해, 네 가지 주요 시스템적 요인을 단계별로 분리하였다. 첫째, 원초적 우주선(프로톤·헬륨) 스펙트럼을 Shikaze, BESS, CREAM, ATIC, PAMELA 등 최신 데이터에 맞춰 여러 가지 피팅 형태로 적용했으며, 특히 PAMELA 데이터에 대한 새로운 파라미터화가 방출량에 10%‒30% 정도의 변동을 야기함을 확인했다. 둘째, 전파 모델은 Donato·et al.이 제시한 최소·중간·최대 세 가지 파라미터 집합(min, med, max)을 사용했으며, 이는 1 kpc 이내에서 우주선 밀도 구배가 최대 10%까지 달라질 수 있음을 보여준다. 셋째, 은하간 헬륨‑수소 비율(n_He/n_H)은 은하 중심에서 0.111, 외곽에서 0.087, 태양 근처에서는 약 0.097으로 추정되는데, 이 변동이 전체 감마선 신호에 약 5% 정도의 영향을 미친다. 넷째, γ‑ray 생성 단면적은 Kamae et al.와 Huang et al. 두 가지 Monte‑Carlo 기반 모델을 비교했으며, 특히 Huang 모델은 50%에 달하는 차이를 보였다. 이러한 네 요인을 모두 고려하면, 이론적 불확실성은 전체적으로 45%‒50%에 이르며, 이는 Fermi‑LAT이 보고한 4%‒30% 통계적 오차와 동등하거나 그보다 크다. 또한, 핵강화 인자 𝜖_M을 에너지와 위치에 따라 고정값(예: 1.84)으로 사용하는 전통적 접근법은, α/프로톤 비율이 고에너지에서 증가한다는 최신 관측 결과와 충돌한다. 실제로 𝜖_M은 1 GeV에서 1.5 TeV까지 30% 이상 변동하므로, 모델링 시 에너지 의존성을 명시적으로 포함해야 한다는 결론에 도달한다.