RHESSI 위성으로 본 감마선 폭발의 군집 구조

초록

본 연구는 BATSE 자료에서 제시된 감마선 폭발(GRB) 세 번째 군집 존재 여부를 검증하기 위해 RHESSI GRB 카탈로그를 이용하였다. 지속시간과 경도비(Hardness Ratio)를 기준으로 통계적 군집 분석을 수행했으며, 결과적으로 RHESSI 데이터에서는 두 개의 주요 군집(장·단)만이 유의미하게 구분되고, 제3군집은 통계적으로 의미가 없음을 확인하였다.

상세 요약

RHESSI(재조정 고에너지 태양 관측 위성)는 2002년부터 2018년까지 1 MeV 이하의 감마선 및 X선을 관측한 뒤, GRB 탐지를 위한 부가적인 트리거 시스템을 운영하였다. 본 논문은 RHESSI가 기록한 427개의 GRB 사건을 표본으로 삼아, 기존 BATSE 연구에서 제시된 “중간 길이” 혹은 “중간 경도” 군집이 실제로 존재하는지를 검증하고자 했다.

우선, 각 사건에 대해 T90(90% 플루언스가 방출되는 시간) 값을 추출하고, 두 개의 에너지 밴드(25–50 keV와 50–100 keV)에서의 플루언스를 이용해 경도비(HR = F₍50‑100₎/F₍25‑50₎) 를 계산하였다. 데이터 전처리 단계에서는 신호‑대‑잡음비(S/N) > 5인 사건만을 남겨, 통계적 편향을 최소화하였다.

통계적 군집 검증에는 세 가지 방법을 병행했다. 첫째, 로그 변환된 T90 분포에 대해 가우시안 혼합 모델(GMM)을 1~3개의 컴포넌트로 피팅하고, Akaike 정보 기준(AIC)과 Bayesian 정보 기준(BIC)으로 최적 모델을 선정하였다. 둘째, 다변량 K‑means 클러스터링을 적용해 T90와 HR을 동시에 고려한 군집 구조를 탐색하였다. 셋째, 비모수적 방법인 히스토그램 커널 밀도 추정(KDE)와 실루엣 점수를 이용해 군집의 분리도를 정량화하였다.

GMM 결과는 두 개 컴포넌트 모델이 AIC/BIC 모두에서 가장 낮은 값을 보였으며, 세 번째 컴포넌트를 추가할 경우 오히려 기준값이 상승하였다. K‑means는 k = 2일 때 실루엣 점수가 0.62로 가장 높았으며, k = 3으로 확대하면 점수가 0.41로 급격히 감소했다. KDE 분석 역시 T90 로그값이 두 개의 뚜렷한 피크(≈0.3 s와 ≈30 s)를 보였고, 중간 구간에서의 밀도는 현저히 낮아 별도의 군집을 형성하지 못했다.

경도비와의 상관관계를 살펴보면, 장시간(GRB)과 단시간(GRB) 사이에 HR 평균값 차이가 통계적으로 유의미했음(두 표본 t‑검정 p < 0.001)하지만, 중간 길이 구간(2–10 s)에서는 HR 분포가 두 주요 군집의 연속선상에 위치해 별도 군집을 정의하기에 충분한 차이를 보이지 않았다.

결과적으로, RHESSI 데이터는 BATSE에서 보고된 “세 번째 군집”이 관측 편향이나 샘플 크기 차이로 인한 통계적 유희일 가능성을 시사한다. 두 개의 주된 군집(짧고 단단한, 길고 부드러운)만이 물리적으로 의미 있는 구분을 제공한다는 결론에 도달하였다.