근처 초신성 잔해가 만든 전자·양전자 급증: PAMELA 이상 현상의 새로운 해석

초록

PAMELA와 Fermi‑LAT이 관측한 양전자 비율 상승과 전자·양전자 총 플럭스 초과를, 인근 초신성 잔해(SNR)에서 2차 전자·양전자가 확산‑충격 가속(DSA) 과정을 겪으며 더 단단한 스펙트럼을 갖게 되는 모델로 설명한다. SNR의 불연속적·확률적 분포에 따른 플럭스 변동을 일반화된 중심극한정리와 안정분포(stable distribution)로 정량화하고, γ‑레이 관측으로부터 얻은 SNR 파라미터를 이용해 자유 파라미터를 고정한다. 모델은 전체 e⁺+e⁻ 스펙트럼을 정확히 재현하고, 독립적인 예측으로 PAMELA 양전자 비율 상승을 재현한다. 또한 고에너지에서 B/C 비율이 상승할 것을 예측해, 향후 AMS‑02·IceCube 관측으로 검증 가능성을 제시한다.

상세 분석

본 논문은 PAMELA가 보고한 양전자 비율 상승과 Fermi‑LAT이 측정한 전자·양전자 총 플럭스가 기존의 1차 전자 가속 모델만으로는 설명되지 않음을 지적한다. 저자들은 초신성 잔해(SNR)를 주요 가속원으로 가정하고, 여기서 2차 전자·양전자가 확산‑충격 가속(DSA) 과정에 참여하면 원래의 부드러운 스펙트럼보다 훨씬 단단한 전력 지수를 갖게 된다고 제안한다.

1. 소스의 불연속성

   • SNR는 수백 파섹 규모이지만, 은하 전체에 퍼진 수천 개의 SNR는 공간·시간적으로 이산적이다. 이를 고려해 각 소스 i가 거리 L_i, 시간 t_i 전후에 방출한 플럭스를 Green 함수 G_disk(E, L_i, t_i) 로 표현하고, 전체 플럭스 J(E)=c/(4π)∑_i G_disk 로 합산한다.
   • L과 t를 독립 확률변수로 두고, L은 r₁≤L<r₂ 구간에서 f_L=2L/(r₂²−r₁²) 로, t는 일정한 발생률을 가정해 f_t=1/t_max (0≤t≤t_max) 로 설정한다. 여기서 t_max=(b₀E_min)⁻¹, b₀는 에너지 손실 계수(b(E)=b₀E²)이다.
   • 중앙극한정리(CLT)가 적용되지 않을 경우, Z=G_disk의 확률밀도 f_Z는 장거리 꼬리를 가지므로 분산이 발산한다. 대신 일반화된 CLT와 안정분포(stable distribution)를 이용해 평균 µ_J와 양쪽 68 %·90 % 신뢰구간을 추정한다. 이는 에너지 증가에 따라 변동폭이 커지는 특성을 만든다.

2. 1차 전자 주입

   • 1차 전자 스펙트럼 R_{e⁻}=R₀ E^{−γ_e} exp(−E/E_cut) 로 가정하고, γ_e≈2.4, E_cut≈20 TeV를 γ‑레이 관측(SNR의 하드 스펙트럼)으로부터 도출한다.
   • 총 주입 에너지 R₀≈1.8×10⁵⁰ GeV⁻¹는 PAMELA 10 GeV 전자 플럭스에 맞추어 결정했으며, 이는 약 7×10⁴⁷ erg에 해당한다.

3. 2차 전자·양전자 가속

   • 전통적으로 2차 전자·양전자는 은하 전파 과정에서만 생성돼 부드러운 스펙트럼을 가진다. 저자들은 SNR 내부에서의 DSA가 2차 입자에도 적용될 경우, 가속된 2차 입자의 스펙트럼이 R_{±}≈R₀ p^{−γ}