[논문 리뷰] Constraining models for the origin of ultra-high-energy cosmic rays with a novel combined analysis of arrival directions, spectrum, and composition data measured at the Pierre Auger Observatory

Pierre Auger Collaboration, Abdul Halim, A.외 373명|arXiv (Cornell University)|2023. 01. 01.

Astrophysics and Cosmic Phenomena인용 수 4

한 줄 요약

이 논문은 초고에너지우주선(UHECR) 기원 모델을 제약하기 위해 피에르 아제르 관측소에서의 도착 방향, 에너지 스펙트럼, 조성 데이터를 새로운 통합 분석 방식으로 제시한다. 여러 데이터 스트림을 통합한 통일된 가능도 프레임워크를 통해, 별 형성률(SFR)-유사한 원천 분포가 관측과 가장 일치하는 것으로 밝혀졌으며, 이는 이전의 단일 데이터 분석에 비해 모델 제약 조건을 크게 향상시킨다.

ABSTRACT

The combined fit of the measured energy spectrum and shower maximum depth distributions of ultra-high-energy cosmic rays is known to constrain the parameters of astrophysical models with homogeneous source distributions. Studies of the distribution of the cosmic-ray arrival directions show a better agreement with models in which a fraction of the flux is non-isotropic and associated with the nearby radio galaxy Centaurus A or with catalogs such as that of starburst galaxies. Here, we present a novel combination of both analyses by a simultaneous fit of arrival directions, energy spectrum, and composition data measured at the Pierre Auger Observatory. We find that a model containing a flux contribution from the starburst galaxy catalog of around 20% at 40 EeV with a magnetic field blurring of around $20^\circ$ for a rigidity of 10 EV provides a fair simultaneous description of all three observables. The starburst galaxy model is favored with a significance of $4.5σ$ (considering experimental systematic effects) compared to a reference model with only homogeneously distributed background sources. By investigating a scenario with Centaurus A as a single source in combination with the homogeneous background, we confirm that this region of the sky provides the dominant contribution to the observed anisotropy signal. Models containing a catalog of jetted active galactic nuclei whose flux scales with the $γ$-ray emission are, however, disfavored as they cannot adequately describe the measured arrival directions.

연구 동기 및 목표

• 피에르 아제르 관측소에서의 도착 방향, 에너지 스펙트럼, 조성 데이터를 통합하는 통일된 통계적 프레임워크를 개발하기 위해.
• 단일이고 일관된 가능도 접근법을 사용하여 초고에너지우주선(UHECR) 기원에 대한 경쟁적 천체물리 모델을 테스트하기 위해.
• 평탄한, SFR-유사한, 또는 AGN-유사한 세 가지 원천 분포 모델 중에서 관측된 UHECR 데이터를 가장 잘 설명하는 모델을 특정하기 위해.
• 실험적 체계적 불확실성(에너지 스케일, Xmax 재구성, 조성 등)이 모델 선택 및 파rameter 추정에 미치는 영향을 정량화하기 위해.
• 스펙트럼 지수, 절단 반경, 조성 진화 등의 UHECR 원천 특성에 대한 제약 조건을 향상시키기 위해.

제안 방법

• 도착 방향 집합, 에너지 스펙트럼 형태, 그리고 에너지에 따른 쇼어의 최대 깊이(Xmax)를 동시에 모델링하는 연합 가능도 함수를 구성한다.
• 모든 원천 분포 모델과 UHECR 특성의 전반적인 매개변수 공간을 탐색하기 위해 베이지안 프레임워크와 마르코프 체인 몬테카를로(MCMC) 샘플링을 사용한다.
• 원천 분포는 평탄한, 별 형성률(SFR)-유사한, AGN-유사한 것으로 모델링되며, 적색편이에 따라 변화하는 주입률을 가진다.
• 에너지 스케일, Xmax 재구성, 조성에 대한 체계적 불확실성은 가우시안 사전분포를 가진 부수적 매개변수로 포함된다.
• 파워-법칙 스펙트럼에 절단이 포함되며, 조성 진화는 수소(H), He, C, Ne, Si, Fe 핵의 에너지에 따라 변화하는 분율로 모델링된다.
• 모델 비교는 베이지안 정보 기준(BIC)과 deviance 정보 기준(DIC)을 사용하며, 적합도 평가에는 로그 가능도 및 deviance 통계량을 사용한다.

실험 결과

연구 질문

• RQ1 피에르 아제르에서의 통합 UHECR 데이터를 바탕으로, 평탄한, SFR-유사한, 또는 AGN-유사한 원천 분포 모델 중 어느 것이 가장 잘 설명하는가?
• RQ2 에너지 스케일 및 Xmax 재구성에 대한 체계적 불확실성이 추정된 매개변수와 모델 선택에 어떤 영향을 미치는가?
• RQ3 방향, 스펙트럼, 조성의 통합 분석 하에서 UHECR의 최적 적합 스펙트럼 지수와 절단 반경은 무엇인가?
• RQ4 조성 데이터를 포함함으로써, 스펙트럼 전용 또는 도착 방향 전용 분석에 비해 UHECR 원천 집단에 대한 제약 조건은 어떻게 향상되는가?
• RQ5 스펙트럼 및 조성 데이터와 결합된 도착 방향의 관측된 군집 현상은 어떤 의미를 갖는가?

주요 결과

• SFR-유사한 원천 분포 모델이 데이터에 가장 잘 맞으며, 평탄한 모델 대비 베이지안 정보 기준(BIC)에서 20.4 단위, AGN-유사 모델 대비 12.7 단위 향상되었다.
• SFR-유사 모델의 최적 적합 스펙트럼 지수는 γ = −3.09이며, 90% 신뢰구간은 [−3.23, −3.09]로, 경직된 스펙트럼과 일치한다.
• 절단 반경 Rcut/V는 log10 단위로 18.10 ± 0.02로 제약되며, 이는 원천 분포가 너무 빽빽하거나 너무 희박하지 않음을 시사한다.
• 분석에 체계적 불확실성을 포함함으로써 deviance 정보 기준(DIC)이 크게 감소하였으며, SFR-유사 모델의 DIC는 134.1로 평탄한 모델의 160.5보다 낮다.
• 조성 분석 결과, 고에너지에서 무거운 핵이 강하게 선호되며, 10^20 eV 에서 철 계열 분율이 약 3.8 × 10−2로 증가하여 고에너지에서 혼합 조성과 일치한다.
• 통합 분석은 원천 주입률과 조성 진화에 대한 불확실성을 감소시켰으며, SFR-유사 모델은 테스트한 모든 모델 중에서 가장 높은 로그 가능도(−274.6)와 가장 낮은 deviance(197.3)를 보였다.