MAGIC 망원경으로 블레이저 주변 확장 감마선 방출 탐색

초록

본 연구는 MAGIC 단일 망원경을 이용해 블레이저 Mrk 421와 Mrk 501 주변에 존재할 수 있는 감마선 확장(halo) 현상을 조사하였다. 새로운 Random Forest 기반 DISP 재구성 방법으로 각 이벤트의 입사 방향을 정밀하게 추정하고, θ² 분포를 통해 점원천과 확장원천을 구분하였다. 분석 결과 두 대상 모두에서 300 GeV 이상의 에너지 임계에서 4 %–7 % 수준 이하의 확장 방출만이 존재할 가능성을 제시했으며, 현재 데이터로는 확장 감마선 halo를 검출하지 못했다.

상세 분석

이 논문은 VHE(very‑high‑energy) 감마선이 은하간 배경광(CMB, EBL)과 상호작용해 전자‑양전자 쌍을 생성하고, 이후 역컴프턴 산란을 통해 2차 감마선을 방출하는 전자기 캐스케이드 과정을 기반으로, 은하간 자기장(EGMF)의 세기에 따라 관측자에게 재지향되는 ‘halo’ 현상이 나타날 수 있음을 이론적으로 제시한다. 특히 EGMF가 10⁻¹² G 수준으로 매우 약할 경우, 생성된 전자‑양전자 쌍이 크게 산란되지 않아 원래 광자와 거의 같은 방향으로 2차 감마선을 방출하게 되며, 이때 발생하는 감마선은 원천의 점상 이미지보다 넓은 θ² 분포를 만든다.

MAGIC은 17 m 반경의 대형 반사경과 고감도 광센서를 갖춘 IACT(이미징 대기 체렌코프 망원경)로, 현재는 단일 망원경 모드와 85 m 간격의 두 번째 망원경을 추가한 스테레오 모드가 운영 중이다. 본 연구에서는 단일 망원경 데이터를 사용했으며, 기존의 Hillas 파라미터 기반 DISP 추정법을 개선하기 위해 Random Forest(RF) 회귀 모델을 도입하였다. RF DISP은 이미지의 타이밍 그라디언트와 형상 파라미터를 다차원적으로 결합해 입사 방향을 추정함으로써, 전통적인 DISP에 비해 θ² 분포의 폭을 20 %–30 % 정도 축소시켰다. 이는 감마선의 각도 해상도를 크게 향상시켜, 0.1° 수준의 PSF(41 % 포함 반경) 내에서 미세한 확장 신호를 탐지할 수 있게 한다.

시뮬레이션에서는 점원천(100 % point‑like)과 반지름 0.1°, 0.2°, 0.3°인 평탄 원반형 확장원을 혼합한 모델을 구축하고, 각각의 θ² 누적 분포를 비교하였다. 확장 반경이 PSF보다 작으면 점원천과 구분이 어려우며, 반경이 클수록 θ² 분포가 뚜렷하게 넓어지는 특징을 확인했다. 실제 데이터 분석에서는 θ₁, θ₂ 구간을 최적화해 ‘f = N(θ₁<θ<θ₂)/N(0<θ<θ₁)’ 비율을 계산하고, 이를 크랩 성운(점원천)의 동일 비율과 비교하였다.

관측 데이터는 Mrk 501(2008년 4‑5월, ON/OFF 모드, 26 h 유효시간, 평균 플럭스 ≈ 0.15 C.U.)와 Mrk 421(2007‑2009년, wobble 모드, 32 h, 평균 플럭스 ≈ 1.4 C.U.)에 대해 각각 저고도(<30°) 조건만을 선택하였다. θ² 분포는 모두 크랩 성운과 일치했으며, χ²/ndf는 Mrk 421에서 10.7/11, Mrk 501에서 3.8/11로 통계적으로 유의미한 차이를 보이지 않았다. 이를 바탕으로 다양한 확장 반경과 프로파일(β = 1, 0, ‑1, ‑2)에 대해 95 % 신뢰수준의 상한을 산출했으며, Mrk 501에 대해서는 <4 % C.U., Mrk 421에 대해서는 <7 % C.U. 수준의 확장 플럭스만이 허용되었다.

결론적으로, 현재 MAGIC 단일 망원경 데이터와 향상된 RF DISP 분석법으로는 블레이저 주변의 감마선 halo를 검출하지 못했으며, 향후 스테레오 모드와 차세대 CTA(Čerenkov Telescope Array)와 같은 고해상도 장비를 이용하면 더 낮은 플럭스와 작은 반경의 halo 탐지가 가능할 것으로 기대된다.