ARGO‑YBJ가 포착한 테라전자볼트 감마선 원천 MGRO J1908+06의 새로운 모습

초록

ARGO‑YBJ는 4년간의 관측을 통해 MGRO J1908+06를 6.2σ의 통계적 유의성으로 검출하였다. 신호 중심은 펄서 PSR J1907+0602와 일치하며, 2차원 가우시안 모델을 적용한 결과 원천 반경은 σ = 0.49 ± 0.22°로 H.E.S.S.와 일치한다. 1–20 TeV 구간에서 측정된 스펙트럼은 dN/dE = 6.1 ± 1.4 × 10⁻¹³ (E/4 TeV)⁻²·⁵⁴ ± 0.36 ph cm⁻² s⁻¹ TeV⁻¹이며, Milagro와 비슷하지만 H.E.S.S.가 제시한 플럭스보다 2–3배 높다. 장기간에 걸친 안정적인 과잉률은 이 원천이 PSR J1907+0602의 지속적인 펄서 풍선(PWN)임을 강하게 시사한다.

상세 분석

ARGO‑YBJ는 티베트 고원(4300 m)에 설치된 대형 공기샤워 검출기로, 넓은 시야와 연속적인 관측이 가능한 장점이 있다. 이번 연구에서는 2007년부터 2011년까지 약 4년간 수집된 데이터를 이용해 MGRO J1908+06를 분석했으며, 전체 관측 시간은 약 1.5 × 10⁸ 초에 달한다. 감마선 후보 이벤트는 전자기파와 양성자·중성자 배경을 구분하기 위해 ‘핵심-주변’ 매개변수를 활용한 다중 레이어 샤워 패턴을 적용했으며, 최적화된 이벤트 선택 기준을 통해 신호 대 잡음비를 극대화하였다.

통계적 유의성은 6.2σ로, Milagro가 2007년에 보고한 5.3σ와 비교해도 충분히 강력한 검출이다. 위치 측정은 신호 최대점이 RA = 286.8°, Dec = 6.1°(J2000) 근처에 있음을 보여주며, 이는 Fermi‑LAT이 탐지한 펄서 PSR J1907+0602와 겹친다. 원천의 공간적 확장은 2차원 가우시안 모델을 적용해 σ = 0.49 ± 0.22°로 추정했는데, 이는 H.E.S.S.가 보고한 0.34° ± 0.03°와 통계적으로 일치한다.

스펙트럼 분석은 1 TeV에서 20 TeV까지의 에너지 구간에서 수행되었으며, 전력법칙 형태의 피팅 결과는 지수 -2.54 ± 0.36, 정규화 계수 6.1 ± 1.4 × 10⁻¹³ ph cm⁻² s⁻¹ TeV⁻¹이다. 이 스펙트럼은 Milagro가 제시한 (E/10 TeV)⁻²·⁶⁸와 거의 일치하지만, H.E.S.S.가 보고한 (E/1 TeV)⁻²·¹⁰에 비해 2–3배 높은 플럭스를 보인다. 이러한 차이는 관측 방법(대기 전리층 샤워 vs. 대기 전자 광학)과 에너지 재구성 방식의 차이, 그리고 원천의 시간적 변동성 여부를 검토하게 만든다. 그러나 ARGO‑YBJ는 4년간 연속적으로 일정한 과잉률을 기록했으며, 월별·연도별 플럭스 변동이 통계적 오차 범위 내에 머물러 ‘steady’한 방출원임을 뒷받침한다.

펄서 PSR J1907+0602는 스핀다운 파워 2.8 × 10³⁶ erg s⁻¹을 가지고 있으며, 이와 연관된 펄서 풍선(PWN) 모델에 따르면 전자들이 강한 자기장과 복사장에 의해 인버스 컴프턴 및 시냅트론 방출을 일으켜 TeV 감마선을 생성한다. ARGO‑YBJ가 측정한 1 TeV 이상 통합 광도는 약 1.8배 크리브(크래브) 광도와 동등하며, 이는 이 PWN이 매우 효율적인 입자 가속기임을 시사한다.

결론적으로, ARGO‑YBJ의 관측은 MGRO J1908+06가 Milagro와 H.E.S.S.가 각각 보고한 특성을 모두 포괄하는 ‘연속적인’ TeV 방출원임을 확인한다. 위치와 확장성, 스펙트럼 형태, 그리고 장기적인 플럭스 안정성은 이 원천이 PSR J1907+0602와 연관된 강력한 펄서 풍선임을 강하게 뒷받침한다. 향후 CTA와 같은 차세대 고감도 텔레스코프를 통한 정밀 이미지와 시간 변동성 조사, 그리고 다파장(라디오·X‑ray·γ‑ray) 연계 관측이 이 PWN의 입자 가속 메커니즘과 주변 환경을 이해하는 데 핵심적인 역할을 할 것이다.