AGN와 초고에너지 우주선의 상관관계 통계분석

초록

본 논문은 파리오거 증강 관측소(PAO)와 AGASA에서 측정된 초고에너지 우주선(UHECR)과 100 Mpc 이내의 활동은성핵(AGN) 사이의 위치 상관을 검증하기 위해, 2차원 방위 분포를 1차원으로 축소한 뒤 Kolmogorov‑Smirnov(KS) 검정을 적용하는 새로운 통계 방법을 제시한다. 결과는 PAO 데이터가 순수한 등방성 혹은 전적으로 AGN 기원이라는 두 극단 가설을 모두 배제하지만, 적절한 등방성 배경 혹은 큰 스메어링을 포함하면 AGN 기원 모델이 타당함을 보여준다. 반면 AGASA 데이터에서는 유의미한 상관을 찾지 못한다.

상세 분석

이 연구는 두 단계의 통계 절차를 도입한다. 첫 번째는 관측된 UHECR의 방위각(α, δ) 분포를 적도좌표계에서 원점(관측소)까지의 각거리와 방위각 두 개의 1차원 히스토그램으로 변환하는 ‘거리‑방위 변환’이다. 이렇게 하면 복잡한 구면 위의 비등방성 구조를 1차원 누적분포 함수(CDF)로 압축할 수 있다. 두 번째 단계는 변환된 CDF와 모델이 예측하는 CDF 사이의 최대 차이값 D를 구해 KS 검정의 p‑값을 산출한다. 모델은 두 파라미터, 즉 AGN 기원 비율 f와 스메어링 각도 σ를 갖는 혼합 분포로 정의된다. f = 1이면 모든 UHECR가 선택된 AGN에서 방출된 것으로 가정하고, σ는 각 UHECR가 AGN에서 출발해 자기장 등에 의해 흐트러지는 각도이다. f < 1이면 나머지 (1 − f) 비율은 완전 등방성 배경으로 처리한다.

PAO 데이터(27건, E > 5.7 × 10¹⁹ eV)에 대해 f = 1, σ ≈ 0°인 순수 AGN 모델은 KS p‑값이 0.001 이하로, 관측과 크게 불일치함을 보인다. 반대로 f = 0, 즉 완전 등방성 모델도 p‑값이 0.02 이하로 배제된다. 그러나 f ≈ 0.5 ~ 0.7, σ ≈ 10°~30° 범위에서 p‑값이 0.2 ~ 0.5로 상승하여, 적당한 등방성 혼합과 중간 정도 스메어링을 포함하면 모델이 데이터를 충분히 설명한다는 결론에 도달한다. 특히 40 ~ 60 Mpc 거리 구간에 있는 AGN만을 선택하면 σ ≈ 5°~15° 정도의 작은 스메어링만으로도 p‑값이 0.4 이상이 되어, 별도의 등방성 성분 없이도 좋은 적합을 보인다.

AGASA 데이터(57건, E > 4 × 10¹⁹ eV)에서는 동일한 분석을 적용했을 때, 모든 파라미터 조합에서 p‑값이 0.05 이하에 머물러, AGN와의 상관관계가 통계적으로 유의미하지 않음을 확인한다. 이는 관측 영역(남반구)과 에너지 임계값 차이, 혹은 실험별 시스템오차가 결과에 영향을 미쳤을 가능성을 시사한다.

이 논문은 2차원 구면 데이터의 1차원 변환을 통해 KS 검정을 적용하는 방법론을 제시함으로써, 복잡한 천체물리학적 상관관계를 보다 직관적이고 계산 효율적으로 평가할 수 있음을 보여준다. 또한, AGN가 UHECR의 잠재적 근원이라는 가설을 완전히 배제하지는 않으며, 등방성 배경과 자기장에 의한 스메어링을 동시에 고려해야 함을 강조한다.