VERITAS가 밝힌 HESS J1857+026의 다파장 특성

초록

VERITAS는 HESS J1857+026을 30 시간 이상 관측해 확장된 VHE 감마선 원천임을 확인하고, 가우시안 형태(반지름 0.24°)와 전력법 스펙트럼(Γ≈2.46, N₀≈5.5×10⁻¹² TeV⁻¹ s⁻¹ cm⁻²)을 측정했다. 방사선 중심에서 거리별 플럭스 분포를 이용해 전자 확산 길이 θ_d≈1.03°를 도출하고, 거리 5.5 kpc를 가정하면 확산계수 D≈3.5×10²⁸ cm² s⁻¹(≈10⁻¹배)임을 추정했다. 이는 은하 평균보다 현저히 낮으며, 중년 펄사 주변에서 억제된 확산이 일반적임을 뒷받침한다.

상세 분석

VERITAS는 4개의 IACT를 이용해 2008‑2016년 사이에 HESS J1857+026 주변을 30 시간 이상 관측하였다. 데이터는 품질 검증 후 평균 고도 58°, 평균 오프셋 1.21°에서 분석되었으며, EventDisplay와 Gammapy 파이프라인을 통해 이벤트 재구성 및 배경 모델링이 수행되었다. 소스 모델링은 전력법 스펙트럼(지수 Γ)과 대칭 가우시안 공간 분포를 사용했으며, 최적화는 최대우도법으로 진행되었다. 결과는 0.3‑10 TeV 구간에서 15.3σ의 유의성을 보였고, 중심 좌표는 RA 284.37°, Dec 2.72°이며, 확장 반지름은 0.24°±0.015°이다. 스펙트럼은 Γ=2.46±0.06, 정규화 N₀=(5.5±0.43)×10⁻¹² TeV⁻¹ s⁻¹ cm⁻²이며, 10 TeV 이상에서는 상한만 제시된다.

방사선의 방사형 프로파일을 구축해 플럭스가 각도 θ에 따라 f(θ)=f₀ θ/(θ+0.06θ_d) exp(-θ²/θ_d²) 형태를 따른다고 가정하였다. 최적화된 θ_d=1.025°±0.317°는 전자 확산 길이 d=θ_d × l≈98 pc(거리 l=5.5 kpc)와 일치한다. 확산계수는 d=2√(D t) 관계를 이용해 D≈(3.5±2.2)×10²⁸ cm² s⁻¹(펄사 연령 t≈21 kyr)로 추정한다. 이는 Geminga와 같은 펄사 주변에서 보고된 억제된 확산(D≈10²⁶‑10²⁷ cm² s⁻¹)보다 약간 높지만, 은하 평균(D≈10²⁹ cm² s⁻¹)보다 한 단계 낮다. 전자 냉각 시간은 1 μG 정도의 자기장 가정 시 수만 년 수준으로, 펄사 연령과 비슷해 전자들이 충분히 방출·확산될 수 있음을 시사한다.

이러한 결과는 HESS J1857+026이 주로 레프톤(PWN) 메커니즘에 의해 VHE 감마선을 방출한다는 가설을 지지한다. 또한 고에너지(>1 TeV)에서 북쪽 부가 구조가 나타나는 점은 전자 에너지 의존적 확산 혹은 별도의 소스 존재 가능성을 제시한다. 향후 Fermi‑LAT, HAWC, 라디오·X‑ray 데이터와 결합한 다파장 분석이 필요하다.