극단값 분석으로 본 FRB 20250316A의 희귀성

초록

본 논문은 CHIME/FRB 2차 카탈로그의 피크 플럭스와 플루언스 데이터를 일반화 극단값 분포(GEV) 모델에 적용해, 비반복 FRB 20250316A가 얼마나 드문 사건인지 정량화한다. 베이지안 MCMC 추정 결과, 피크 플럭스는 ξ > 0인 Fréchet형으로 1σ 신뢰구간에서 약 600년 회귀주기를 보이며, 플루언스는 전체 샘플에서는 Fréchet형(≈50년)이나, 세 개의 다른 극단값을 제외하면 Weibull형(ξ < 0)으로 상한이 174 Jy ms에 불과해 20250316A는 그 한계를 크게 초과한다는 결론을 얻는다.

상세 분석

본 연구는 CHIME/FRB 2차 카탈로그에 포함된 비반복 FRB들의 피크 플럭스와 플루언스 시계열을 30일 길이의 균등 블록으로 나누어 각 블록에서 최대값을 추출, 블록 최대값 집합을 GEV(Generalized Extreme Value) 분포에 적합하였다. GEV는 위치(µ), 규모(σ), 형태(ξ) 세 파라미터로 정의되며, ξ의 부호에 따라 Fréchet(ξ > 0, 무한 상한, heavy‑tail), Gumbel(ξ ≈ 0, 지수적 감소), Weibull(ξ < 0, 유한 상한, light‑tail) 형태로 구분된다. 베이지안 프레임워크 하에 uniform priors(µ,σ)와 ξ에 대해 평균 0, 표준편차 0.5인 정규 사전분포를 설정하고, emcee 패키지의 affine‑invariant MCMC 샘플러로 사후분포를 탐색하였다. 로그우도 함수는 GEV의 확률밀도함수를 사용해 각 데이터 포인트에 대해 계산했으며, ξ = 0인 경우 Gumbel 한계를 별도 처리하였다.

피크 플럭스에 대한 사후분포는 ξ 평균값 0.32 ± 0.14로, 명확히 양수이며 Fréchet형을 지지한다. 이 경우 1σ(68 %) 신뢰구간에서 회귀주기(T) ≈ 593 년, 2σ(95 %)에서는 약 73 년으로 추정된다. 이는 현재 5년간 관측된 CHIME 데이터에 비해 매우 드문 사건임을 의미한다. 플루언스 분석에서는 전체 샘플에 세 개의 추가 극단값(FRB 20200723B, 20220222B, 20210922C)이 포함될 경우 ξ 평균 0.37 ± 0.18으로 역시 Fréchet형을 보이며, 1σ 회귀주기 14 년, 2σ 51 년으로 나타난다. 그러나 QQ 플롯에서 데이터가 95 % 신뢰밴드를 크게 벗어나는 구간이 존재해 모델 적합도가 떨어진다.

따라서 세 개의 외부 극단값을 제외하고 재분석하면 ξ 평균 − 0.24 ± 0.12, 즉 Weibull형이 가장 가능성이 높다. Weibull형에서는 상한 µ − σ/ξ ≈ 174.5 Jy ms가 계산되며, FRB 20250316A의 플루언스 1.7 kJy ms는 이 상한을 약 10배 초과한다. 사후분포는 ξ가 양·음 모두를 포함하는 폭넓은 형태를 보여, 현재 데이터만으로는 tail 형태를 확정하기 어렵다는 점을 강조한다.

이러한 결과는 FRB 20250316A가 단순히 기존 FRB 집단의 통계적 변동 범위 내에 있지 않으며, 매우 드문 고에너지 물리 과정이나 별도의 하위 집단에 기인할 가능성을 시사한다. 또한 GEV 분석이 관측 기간이 짧은 전파 트랜시언트 연구에 적용될 때, 블록 길이 선택과 극단값 샘플 수가 결과에 큰 영향을 미친다는 methodological 교훈을 제공한다.