1GeV에서 드러난 감마선 변동성 Fermi LAT 전천구 탐색

초록

Fermi‑LAT 168주 데이터를 이용해 1 GeV 이상 감마선의 장기 변동성을 전천구적으로 조사하였다. Healpix 12288칸으로 나눈 뒤 주간 광자 수와 노출을 비교해 Kolmogorov‑Smirnov 검정으로 일정 플럭스 가설을 검증했다. 유의확률 4×10⁻⁶ 이하(전체 0.05 허위 검출)인 117개의 변동원을 발견했으며, 그 중 27개는 최초 보고된 변동이다. 블레이저 AGN 25개(BL Lac 11, FSRQ 14), 미확정형 AGN 1개, 그리고 펄서 Geminga(PSR J0633+1746)의 변동이 5.1σ 통계적 유의성을 보였다.

상세 분석

본 연구는 Fermi‑LAT이 제공하는 1 GeV 이상 감마선 데이터를 168주(≈3년) 동안 누적해 전천구적인 변동 탐지를 수행한 점이 가장 큰 특징이다. 먼저 Healpix(HEALPix) 패키지를 이용해 구면을 Nside=32 수준으로 분할, 총 12 288개의 등면적 픽셀에 할당하였다. 각 픽셀마다 주간 광자 수를 집계하고, 동일 기간에 대한 관측 노출을 Fermi‑LAT 툴(‘gtltcube’, ‘gtexpcube2’)을 통해 정밀히 계산하였다. 이렇게 얻은 광자‑노출 시계열을 바탕으로 Kolmogorov‑Smirnov(KS) 검정을 적용, “주간 광자 수/노출 비율이 일정한 플럭스 가설”을 검증하였다. KS 통계량이 임계값(p < 4×10⁻⁶) 이하이면 변동성이 있다고 판단했으며, 이는 전체 픽셀 수에 대한 기대 허위 검출 수를 0.05로 제한하는 보수적인 기준이다.

통계적 절차는 두 가지 장점을 가진다. 첫째, KS 검정은 분포 형태에 대한 가정이 필요 없으므로 비정상적인 광자 도착 패턴이나 노출 변동에 강인하다. 둘째, 픽셀 단위의 독립 검정이므로 복잡한 모델링 없이도 전천구적인 스캔이 가능하다. 그러나 한계점도 존재한다. KS 검정은 샘플 크기가 작을 때 민감도가 떨어지며, 특히 광자 수가 매우 적은 고위도 혹은 저감도 영역에서는 변동을 놓칠 위험이 있다. 또한, 픽셀 크기가 고정되어 있기 때문에 강한 변동을 보이는 소스가 픽셀 경계에 걸릴 경우 신호가 분산되어 검출 효율이 감소한다.

결과적으로 117개의 변동 소스를 식별했으며, 이 중 27개는 이전 문헌에 보고되지 않은 새로운 변동 사례다. 변동 소스의 분포를 살펴보면 대부분이 알려진 AGN, 특히 블레이저(blazar) 계열에 집중되어 있다. BL Lac 11개와 FSRQ 14개가 새롭게 변동을 보였으며, 이는 블레이저가 고에너지 감마선에서 빠른 변동성을 보이는 전형적인 특성을 재확인한다. 흥미롭게도 펄서 Geminga(PSR J0633+1746)에서도 5.1σ 수준의 장기 변동이 검출되었는데, 이는 전통적으로 안정적인 고에너지 방출원으로 간주되던 펄서가 실제로도 변동성을 가질 수 있음을 시사한다.

이 연구는 전천구적인 변동 탐색 방법론을 확립함으로써, 기존에 알려진 변동 소스 외에도 미확인 변동성을 가진 후보들을 체계적으로 발굴할 수 있음을 보여준다. 향후에는 픽셀 크기를 가변적으로 조정하거나, 베이지안 블록 모델링과 같은 시간 가변성 검정과 결합해 감도와 정확성을 동시에 향상시킬 여지가 있다. 또한, 변동 소스와 다중파장(라디오, 광학, X‑ray) 관측을 연계해 변동 메커니즘을 물리적으로 해석하는 후속 연구가 기대된다.