상호작용 초신성이 만든 초무릎 우주선

초록

밀도 높은 원주 물질(CSM)과 충돌하는 초신성 폭발이 비공명 스트리밍 불안정에 의해 자기장을 증폭시키고, 무거운 원자핵을 포함한 입자를 확산 충격 가속(Diffusive Shock Acceleration)으로 가속한다. 저자들은 IIn형 초신성을 중심으로 Ib, Ic, Ibn, Icn 등 여러 종류를 모델링해, 10¹⁵–10¹⁷ eV 구간의 관측된 초무릎(CR) 플럭스와 무거운 원소 비율을 성공적으로 재현한다는 결론을 제시한다.

상세 분석

이 논문은 핵심적으로 ‘상호작용 초신성(Interacting Supernova, ISN)’이 초무릎(∼10¹⁵–10¹⁷ eV) 구간의 은하계 우주선(Cosmic Ray, CR) 원천이 될 수 있음을 이론적으로 입증한다. 먼저 저자들은 전형적인 핵융합 전 단계에서 대량 질량 손실을 겪는 거대 별이 형성하는 고밀도 원주 물질(CSM)의 물리적 특성을 정량화한다. CSM의 밀도는 질량 손실률(Ṁ)과 풍속(V_w)으로 정의되며, ρ_CSM ∝ Ṁ / (4πR²V_w) 형태를 취한다. 초신성 폭발 직후, 빠른 SN ejecta가 CSM과 충돌해 형성되는 충격파는 초기에는 복사 매개( radiation‑mediated) 상태이지만, 시간이 흐르면서 광학 깊이가 감소해 비복사성(collisonless) 충격으로 전환된다. 이때 비공명 스트리밍 불안정(NRS, Bell instability)이 급격히 자기장을 증폭시켜 B_sat ≈ (12π ξ_CR ρ_CSM V_sh³/ϕ_c)¹ᐟ² 로 표현된다. 여기서 ξ_CR은 충격 램프 압력 중 CR에 할당되는 비율이며, ϕ_c는 에너지 구간 로그(ln(E_max/E_min))에 해당한다. IIn형 초신성에서는 B_sat이 원래 풍 자기장(B_w)보다 3~5배까지 증폭되며, 이는 ε_B ≈ 10⁻³ 수준의 자기장 효율을 의미한다.

가속 가능한 최대 에너지 E_max은 라모어 반경 r_L = E/(ZeB)와 충격 규모를 비교해 추정된다. 저자들은 E_max ≈ Z e (ε_B Ṁ V_w)¹ᐟ² 형태로 도출했으며, Z=2(헬륨) 기준으로 Ṁ ≈ 10⁻² M_⊙ yr⁻¹, V_w ≈ 100 km s⁻¹일 때 E_max ≈ 5 × 10¹⁶ eV가 된다. 그러나 실제 E_max은 핵분열, 광전분해, 팽창 손실 등 파괴 과정과 탈출 시간(t_esc) 사이의 경쟁에 의해 제한된다. 저자들은 ‘escape‑limited’ 모델을 채택해, t_crit에서 탈출 시간이 가장 짧아지는 순간부터 t_f ≈ 300 일까지의 시간 구간을 적분해 CR 스펙트럼을 계산한다.

다양한 SN 유형별 파라미터(Ṁ, V_w, 화학 조성, 발생 비율 f_SN)를 적용해 IIn, Ibn, Icn, Ib, Ic를 각각 모델링했다. IIn은 가장 높은 Ṁ(∼10⁻² M_⊙ yr⁻¹)과 낮은 V_w(∼100 km s⁻¹)를 가져 가장 큰 B_sat과 E_max을 제공한다. 반면 Ib/Ic는 Ṁ ∼ 10⁻³ M_⊙ yr⁻¹, V_w ∼ 1000 km s⁻¹로 자기장 증폭이 거의 일어나지 않아 초무릎 영역에 기여가 제한적이다. 중성자와 중원소 비율은 CSM의 원소 풍부도와 전리 상태에 따라 달라지며, 저자들은 전리된 이온의 전하 Z_eff를 고려해 가속 효율을 조정했다. 결과적으로 IIn 초신성은 10¹⁵–10¹⁷ eV 구간에서 전체 CR 플럭스의 ≈50–70%를 차지하고, 무거운 원소(Fe, Si 등)의 비율이 증가하는 관측된 추세와 일치한다.

또한, 저자들은 전파 관측과 LHAASO의 고에너지 감마선 데이터와 비교해 모델의 일관성을 검증했다. IIn에서 기대되는 고에너지 감마선(π⁰ 붕괴)과 중성미자 플럭스는 현재 감지 한계에 근접하지만, 차세대 중성미자 망원경(KM3NeT 등)에서 확인 가능할 것으로 예측한다. 마지막으로, 전체 은하계 초신성 발생률과 CR 에너지 수지를 비교해, IIn이 차지하는 비율이 전체 초신성 발생률의 약 10% 정도임에도 불구하고, 높은 에너지 효율(ε_CR ∼ 0.05) 덕분에 초무릎 CR에 지배적인 역할을 할 수 있음을 강조한다.