W28 초신성잔해에서 방출된 확산 양성자에 의해 조명된 분자 구름의 감마선

초록

W28 초신성잔해(SNR)와 인접한 분자구름(MC)에서 관측된 감마선은, SNR 충격전면에서 시간에 따라 방출된 양성자들이 확산하면서 MC에 축적되어 발생한다는 ‘누적 확산 모델’로 설명된다. 저자들은 네 개의 TeV 소스(N, A, B, C)의 스펙트럼을 모델에 맞추어, 동쪽 북동부 소스(N)의 GeV 파워‑브레이크와 남쪽 소스(A, C)의 약한 GeV 방출을 자연스럽게 재현하였다. 소스 A와 C가 SNR 중심으로부터 더 멀리 떨어져 있어 양성자 밀도가 낮아 GeV 감마선이 약하게 나타난다는 결론을 도출한다.

상세 요약

본 논문은 W28 초신성잔해와 그 주변에 존재하는 네 개의 분자구름(소스 N, A, B, C)에서 검출된 감마선 스펙트럼을 ‘누적 확산 모델(Accumulated Diffusion Model)’을 통해 정량적으로 해석한다. 기존 연구에서는 SNR이 폭발 후 일정 시점에 양성자를 방출하고, 이들이 주변 매질에 확산하면서 p‑p 충돌에 의해 감마선을 만든다고 가정했지만, 저자들은 SNR의 전진 충격이 진행되는 전체 시간 동안 지속적으로 양성자가 탈출한다는 점을 강조한다. 이때 탈출 양성자는 에너지 의존적인 확산계수를 갖고, 거리와 시간에 따라 공간적·시간적 분포가 달라진다. 특히, 양성자 탈출률은 SNR의 진화 단계(자유 팽창기 → 라디에이션 단계)와 연관되어, 초기에는 높은 에너지 양성자가 주로 방출되고, 시간이 흐를수록 저에너지 양성자도 점차 방출된다는 물리적 근거를 제시한다.

모델 구현에서는 (1) SNR의 반경 R(t)와 충격속도 V_s(t)를 표준적인 Sedov‑Taylor 해와 관측된 파라미터를 이용해 추정하고, (2) 양성자 탈출 시점 t_esc(E) 를 에너지 의존적인 함수로 설정한다. 탈출 양성자는 3차원 확산 방정식의 해를 이용해 각 MC까지의 전파 시간을 계산하고, (3) 각 MC의 질량·밀도·거리 정보를 바탕으로 p‑p 충돌에 의한 감마선 생성 효율을 구한다. 핵심 파라미터인 확산계수 D(E)=χ D_ISM (E/10 GeV)^δ 에서 χ는 표준 은하계 확산계수 대비 억제 인자이며, δ는 에너지 의존성 지수이다. 저자들은 χ≈0.01, δ≈0.5 정도가 네 개 소스 모두에 일관되게 적용될 수 있음을 보인다.

스펙트럼 피팅 결과, 동쪽 북동부 소스 N은 거리 0.2 kpc 정도로 SNR에 가장 가깝고, 따라서 고에너지 양성자가 풍부하게 충돌해 GeV–TeV 구간에 걸친 부드러운 스펙트럼을 만든다. 특히, 1–10 GeV 구간에서 관측된 스펙트럼 브레이크는 양성자 탈출 시점이 에너지에 따라 달라지는 ‘시간‑에너지’ 효과와, MC 내부에서의 에너지 손실(특히 저에너지 양성자의 감쇠)으로 설명된다. 반면, 남쪽 소스 A와 C는 SNR 중심으로부터 각각 0.5 kpc, 0.6 kpc 정도 떨어져 있어, 양성자 밀도가 약 1/4 수준으로 감소한다. 이로 인해 GeV 대역에서 감마선 플럭스가 거의 검출되지 않으며, TeV 대역에서만 약한 신호가 남는다.

또한, 저자들은 모델이 기존의 ‘단일 탈출·단일 확산’ 시나리오보다 물리적으로 더 타당하다는 점을 강조한다. 실제 SNR는 연속적인 충격 전파와 복합적인 환경(밀도 비균일성, 자기장 구조) 속에서 양성자를 방출하므로, 누적된 확산 양성자 흐름이 주변 MC를 조명한다는 개념이 관측된 복합 스펙트럼을 자연스럽게 재현한다. 이 연구는 SNR‑MC 상호작용을 통한 고에너지 입자 가속·전파 메커니즘을 이해하는 데 중요한 실증적 근거를 제공한다.