On the astrophysical origin of cosmic rays: Constraining the Ultra-High-Energy Cosmic Ray Horizon through Nearby Galaxy Distributions

기술 요약: 인근 은하 분포를 통한 초고에너지 우주선 지평선의 제약

문제 정의
에너지가 $E \ge 8$ EeV인 초고에너지 우주선(UHECR)의 기원은 입자 물리학의 중심적인 미해결 과제로 남아 있다. 관측 증거는 이들의 외계 기원을 시사하지만, 은하 및 외계 자기장에 의한 상당한 편향으로 인해 이들을 기원지로 추적하는 것은 복잡한 문제이다. 또한, 그라이젠-제타톤-쿠즈민(GZK) 효과는 에너지 의존적 지평을 부과하여 UHECR의 전파 거리를 제한한다. 피에르 오제(Pierre Auger) 관측소에 의해 UHECR 도착 방향에서 대규모 쌍극자 이방성이 검출되었으나, 국부 대규모 구조가 이 이방성에 기여하는 유효한 공간적 척도는 아직 잘 제약되지 않았다. 현재의 연구들은 먼 거리의 희미한 은하들을 과소 표집할 수 있는 크기 제한과 같은 관측 편향으로 인해, UHECR과 그 기원 분포 사이의 진정한 관계를 가릴 위험이 있다.

방법론
본 연구는 UHECR의 전파 지평을 제약하기 위해 UHECR 도착 방향과 인근 은하 분포 사이의 상관관계를 이용한 다각적인 통계적 접근 방식을 채택한다.

1. 데이터 선택:

   • UHECR 데이터: 본 연구는 약 85%의 천구($-90^\circ < \delta < 45^\circ$)를 커버하는 $E \ge 8$ EeV 에너지 대역의 피에르 오제 관측소 전천 플럭스 맵을 활용한다.
   • 은하 카탈로그: 분석에는 약 2,300만 개의 소스로 구성된 균질화된 전천 컴필레이션인 GLADE+ 카탈로그를 사용한다. 거리 및 광도 의존적 선택 편향을 완화하기 위해, 체적 제한 샘플(GLADE+VLS)을 구축하였다. 이 샘플은 $1200 < cz < 10000$ km s$^{-1}$ (광도 거리로 약 18 ~ 150 Mpc에 해당)의 방사 속도 범위 내에서 절대 B-밴드 등급 $M_B < -18$인 은하들을 포함한다. 샘플은 오제 관측 영역과 일치시키기 위해 은하 평면(Zone of Avoidance) 및 적위 $45^\circ$ 이상의 영역을 제외한다.

2. 쌍극자 분석:

   • 저자들은 은하 샘플의 적경에 대해 레이리 분석(Rayleigh analysis)을 수행하여 쌍극자 진폭과 위상을 계산하였다.
   • 이는 은하 쌍극자가 오제 협력단이 보고한 UHECR 쌍극자와 비교하여 어떻게 진화하는지 조사하기 위해 연속적인 방사 속도 슬라이스별로 수행되었다.

3. 각 상관관계:

   • 이산적인 은하 위치와 연속적인 오제 플럭스 맵 사이의 각 상관함수 $w(\theta)$를 TreeCorr 추정량을 사용하여 계산하였다.
   • 분석은 거리에 따라 GLADE+VLS를 세 개의 빈($1200-4000$, $4000-7000$, $7000-10000$ km s$^{-1}$)으로 나누어 계층화하였다.

4. 자기 편향 모델링:

   • 은하 자기장(GMF)에 의한 UHECR 궤적의 퍼짐 현상을 고려하기 위해, 본 연구는 Unger & Farrar (2024) GMF 모델을 사용하였다.
   • 대표적인 강성(rigidity) $R = 5$ EV(가벼운 질량에서 중간 질량 핵을 포함)에 대해 CRPropa 3 프레임워크를 사용하여 역추적(backtracking) 시뮬레이션을 수행하였다.
   • 각 은하는 GMF 모델에 의해 예측된 제곱 각 편향 면적의 역수($w_i \propto \theta_i^{-2}$)에 따라 가중치가 부여되는 가중치 스키마를 적용하였다. 이는 자기적 왜곡이 심한 영역의 소스들을 효과적으로 감쇄시킨다.
   • 강성 값($R=2, 8$ EV)과 대안적인 GMF 모델(JF12, KST24)을 사용하여 결과의 견고성을 확인하였다.

주요 결과

• 쌍극자 정렬: 은하 쌍극자와 오제 UHECR 쌍극자 사이의 각 분리는 가장 가까운 은하 집단($cz \lesssim 4000$ km s$^{-1}$, 즉 $\lesssim 60$ Mpc)에서 최소화된다. 이 범위에서 은하 쌍극자의 진폭 또한 최대가 된다. 거리가 이 임계값을 넘어 증가함에 따라 쌍극자 정렬은 크게 약화되며 진폭은 감소한다.
• 거리 의존적 상관관계: 은하와 UHECR 사이의 각 상관관계 신호는 가장 가까운 은하 빈($cz < 4000$ km s$^{-1}$)에 의해 지배된다. 상관관계 진폭은 중간 및 먼 빈에서 급격히 떨어지며, 이는 관측된 이방성이 주로 국부 대규모 구조에 의해 구동됨을 나타낸다.
• 자기 편향의 영향: 상관관계 분석에 GMF 편향 가중치를 포함하면 가중치를 적용하지 않은 경우보다 상관관계 진폭이 2배 이상 향상된다. 이는 자기 편향이 관측된 이방성을 형성하는 데 결정적인 역할을 하며, 이를 보정함으로써 물질 분포와의 더 강력한 물리적 연관성을 드러낼 수 있음을 시사한다.
• 견고성: 결과는 서로 다른 강성 값($R=2, 5, 8$ EV)과 다양한 GMF 모델(UF23, JF12, KST24)에 대해서도 일관되게 나타나며, 이는 상관관계 신호의 향상이 분석의 견고한 특징임을 확인시켜 준다.

의의 및 결론
본 논문은 $E \ge 8$ EeV에서의 대규모 이방성이 약 $\sim 50-60$ Mpc의 특성 거리 내의 소스들에 의해 주로 구동된다고 결론짓는다. 이 척도는 GZK 에너지 손실 및 자기 편향으로부터 도출된 이론적 기대치와 일치한다.

본 연구는 관측된 UHECR 이방성이 인근 외계 물질 분포의 측정 가능한 흔적을 유지하고 있음을 입증한다. 부피 제한 은하 샘플을 활용하고 은하 자기 편향을 고려함으로써, 저자들은 효과적인 UHECR 소스 지평을 관측적으로 결정하였다. 이러한 결과는 UHECR 플럭스 억제와 이방성이 단순히 소스의 고유한 특성이나 먼 거리 전파의 효과만이 아니라, 국부 우주와 연결되어 있다는 해석을 강화한다. 본 연구는 유사한 상관관계 기술을 특정 은하 집단에 적용함으로써 UHECR 전구체의 천체물리학적 본질을 더욱 정밀하게 규명하기 위한 프레임워크를 구축하였다.