중성자별 합병 중력파와 베타 붕괴 연관성 부정

2017년 8월 17일, 고감도 레이저 간섭계인 LIGO와 VIRGO가 최초로 이중 중성자별 합병에서 발생한 중력파 신호 GW170817을 포착하였다. 이 사건 직후, 중성미자 플럭스가 급증할 가능성이 제기되었으며, 일부 연구자는 이러한 중성미자 증가가 β‑붕괴율에 영향을 미칠 수 있다고 주장했다. 특히, Fishbach 등(2018)은 32Si와 36Cl 두 동위원소의 붕괴율 변동이 5시간 구간 내에서 2.5σ 수준의 피어슨 상관을 보였다고 보고했지만, 해당 분석이 임의 구간 선택에 대한 보정을 누락했다는 비판이 제기되었다. 본 논문은 이러한 주장에 대한 독립적인 검증을 목표로, 네덜란드 암스테르담에 위치한 니케프 연구소에서 44Ti, 60Co, 137Cs 세 방사성 소스를 동시에 측정하였다. 각각의 소스는 약 0.8 kBq(44Ti), 0.5 kBq(60Co), 0.8 kBq(137Cs) 수준이며, NaI(Tl) 검출기 4쌍(각각 76 mm × 76 mm)으로 감시되었다. 검출기들은 5 cm 두께의 납 차폐를 받아 외부 방사능으로부터 보호되었으며, 네 번째 검출기 쌍은 환경 배경을 실시간 모니터링하였다. 데이터 수집은 1시간 단위로 이루어졌으며, 각 시간 구간에 대해 전자기 스펙트럼의 전 흡수 피크를 가우시안 형태로 피팅하고, GEANT4 기반 시뮬레이션을 통해 얻은 컴프톤 백그라운드와 환경 방사능을 동시에 고려하였다. 측정 기간이 짧아 반감기에 따른 지수 감소는 무시할 수 있었으며, 따라서 기대값과 표준편차는 단순히 관측된 평균과 표본 분산으로 정의하였다. 통계 분석은 두 소스 간 피어슨 상관계수 r을 계산하고, 해당 r에 대한 양측 p값을 구하는 전통적인 방법을 적용하였다. 초기 로컬 검정에서는 44Ti‑60Co 조합에서 r≈0.95에 해당하는 2.4σ(p=0.015), 44Ti‑137Cs에서 1.3σ(p=0.19), 60Co‑137Cs에서 1.8σ(p=0.07) 정도의 상관이 관측되었다. 그러나 이러한 로컬 유의성은 “Look‑Elsewhere 효과”를 고려하지 않은 결과이므로, 실제 의미를 평가하기 위해 추가 보정이 필요했다. 보정을 위해 저자들은 3시간, 4시간, … 8시간까지 총 5개의 시간 구간을 정의하고, 각 구간마다 모든 소스 조합에 대해 피어슨 검정을 수행하였다. 각 시뮬레이션 트라이얼은 무효 가설 하에서 포아송 분포를 따르는 카운트를 생성했으며, 30 000회의 몬테카를로 실험을 통해 극단적인 p값(가장 작은 p값)의 분포를 구축하였다. 이 분포를 기반으로 관측된 최소 p값이 무효 가설 하에서 얼마나 드물게 나타나는지를 전역 p값으로 환산하였다. 시간 구간 선택에 대한 보정(식 2)과 동시에 세 소스를 동시에 모니터링한 점을 고려한 다중 비교 보정(식 3)을 적용한 결과, 전역 유의성은 44Ti‑60Co 조합에서 최대 1.4σ(p≈0.15) 수준으로 감소하였다. 이는 통계적 우연에 의한 변동으로 해석될 수 있으며, 귀무가설을 기각하기엔 충분치 않다. 실험 환경 변수(온도, 전압, 자기장, 라돈 농도 등)는 측정 기간 내내 1 % 이내의 변동만을 보였으며, 이들 변동이 카운트에 미치는 영향은 2 × 10⁻⁵ 이하로 추정되었다. 따라서 관측된 작은 상관도 실험적 시스템atics에 의해 설명될 가능성은 낮다. 결론적으로, GW170817 중력파 사건이 지구상의 β‑붕괴율에 미치는 영향을 현재의 정밀도와 통계적 방법으로는 검출할 수 없으며, 이전에 보고된 2.5σ 상관은 시간 구간 선택과 다중 검정 보정을 고려하지 않아 과대평가된 결과임을 확인하였다. 저자들은 향후 유사한 연구에서 사전 모델 정의와 블라인드 분석을 통해 인간 편향을 최소화할 것을 권고한다.