[논문 리뷰] Observational constraints on cosmic-ray escape from UHE accelerators

Quentin Luce, Sullivan Marafico외 3명|arXiv (Cornell University)|2022. 07. 17.

Astrophysics and Cosmic Phenomena인용 수 1

한 줄 요약

이 논문은 페르디나우르 관측소의 데이터를 바탕으로 초고에너지 우주선(UHECR) 방출에 대한 현상학적 모델을 시험하며, 전입자 스펙트럼, 암반 이하의 양성자 스펙트럼, 그리고 쇼워 최대 깊이 분포를 동시에 맞추고자 한다. 결과적으로, 은하간 물질 밀도와 일치하는 공간적 소스 분포가 데이터를 잘 재현하며, 핵종에 대해는 경량 스펙트럼 지수(γ ≈ 0.27)가 필요하지만, 탈출하는 양성자에 대해서는 상당히 부드러운 지수(γ ≈ 3.54)가 요구된다. 이는 극단적인 가정에도 불구하고 소스 내 상호작용이 강력한 메커니즘임을 지지한다.

ABSTRACT

Interactions of ultra-high energy cosmic rays (UHECRs) accelerated in specific astrophysical environments have been shown to shape the energy production rate of nuclei differently from that of the secondary neutrons escaping from the confinement zone. Here, we aim at testing a generic scenario of in-source interactions through a phenomenological modeling of the flux and composition of UHECRs. We fit a model in which nucleons and nuclei follow different particle energy distributions to the all-particle energy spectrum, proton spectrum below the ankle energy and distributions of maximum shower depths above this energy, as inferred at the Pierre Auger Observatory. We obtain that the data can be reproduced using a spatial distribution of sources that follows the density of extragalactic matter on both local and large scales, providing hence a realistic set of constraints for the emission mechanisms in cosmic accelerators, for their energetics and for the abundances of elements at escape from their environments. While the quasi mono-elemental increase of the cosmic-ray mass number observed on Earth from ${\simeq}\: 2\:$EeV up to the highest energies calls for nuclei accelerated with a hard spectral index, the inferred flux of protons down to ${\simeq}\: 0.6\:$ EeV is shown to require for this population a spectral index significantly softer than that generally obtained up to now. We demonstrate that modeling UHECR data across the ankle substantiate the conjecture of in-source interactions in a robust statistical framework, although pushing the mechanism to the extreme.

연구 동기 및 목표

• 페르디나우르 관측소의 관측 제약 조건을 이용해 UHECR 가속 환경에서의 일반적인 소스 내 상호작용 시나리오를 시험하기 위해.
• 관측된 에너지 스펙트럼과 조성 변화—특히 약 ~2 EeV 이상에서 질량 수의 증가—가 소스 수준의 상호작용과 탈출 과정에 의해 설명될 수 있는지 판단하기 위해.
• 통합 통계 모델을 사용해 UHECR 가속기에서 탈출할 때의 스펙트럼 지수, 소스 빛의 세기, 원소 조성을 제약하기 위해.
• 모든 입자 스펙트럼, 암반 이하의 양성자 스펙트럼, Xmax 분포를 동시에 맞추는 방식으로 소스 내 상호작용 모델의 강건성을 평가하기 위해.

제안 방법

• 핵종과 양성자 각각 다른 스펙트럼 지수 γp와 γA를 가진 소스에서의 에너지 분포를 가진 현상학적 모델을 구성한다.
• 2차 중성자를 생성하는 소스 내 상호작용을 포함하며, 이들은 자유롭게 탈출하여 붕괴해 관측된 양성자 플럭스에 기여한다.
• 세 가지 핵심 관측량에 대해 모델을 적합한다: 전입자 에너지 스펙트럼, 암반 이하의 양성자 스펙트럼(log10(E/eV) < 19.0), 여러 에너지 영역에 걸친 쇼워 최대 깊이(Xmax) 분포.
• 에어 쇼워 발달과 Xmax 분포 계산을 위해 EPOS-LHC 및 Sibyll2.3c 강입자 상호작용 모델을 사용한다.
• deviance 통계량을 사용한 χ²-최우도 적합을 수행하며, 플럭스(J), 양성자 스펙트럼(Jp), Xmax 자료의 개별 기여를 고려하여 총 deviance를 최소화함으로써 최적의 매개변수를 도출한다.
• 소스의 공간 분포는 대규모 및 국지적 은하간 물질 분포를 따르도록 모델링하여 천체물리학적으로 현실적인 제약 조건을 제공한다.

실험 결과

연구 질문

• RQ1 단일 모델에서 소스 내 상호작용을 고려할 때, 전입자 스펙트럼, 암반 이하의 양성자 스펙트럼, Xmax 분포를 동시에 재현할 수 있는가?
• RQ2 관측된 데이터를 적합하기 위해 탈출 시 양성자와 핵종의 스펙트럼 지수는 무엇이 필요한가? 그리고 전통적인 가속 모델과 비교해보면 어떻게 되는가?
• RQ3 은하간 물질 밀도와 일치하는 소스 분포가 UHECR 데이터에 통계적으로 강건한 적합을 제공하는가?
• RQ4 약 ~2 EeV 이상에서 관측된 질량 수의 준일원소적 증가 현상은 경량 핵종 지수와 부드러운 양성자 지수와 일치하는가?
• RQ5 데이터는 UHECR 소스의 빛의 세기와 탈출 시 원소 조성에 대해 어떤 제약을 가하는가?

주요 결과

• 소스가 국지적 및 대규모에서 은하간 물질 분포를 따를 경우 데이터가 잘 재현되며, 이는 소스 기하학에 대한 현실적인 천체물리학적 제약 조건을 제공한다.
• 탈출하는 양성자에 대한 추정된 스펙트럼 지수는 일반적으로 가정되는 것보다 상당히 부드럽다(γp ≈ 3.54), 반면 핵종에 대해서는 경량 지수(γA ≈ 0.27)가 필요하다.
• 모델은 에너지 증가에 따른 질량 수 증가를 성공적으로 재현하며, 이는 소스 내 상호작용과 차별화된 탈출이 조성 변화를 이끄는 데 핵심적임을 시사한다.
• 양성자와 핵종에 대해 별개의 스펙트럼 지수(γp ≠ γA)를 允허할 경우, 공통 지수를 사용한 경우에 비해 deviance가 130 단위 감소하여 통계적으로 유의미한 향상임을 나타낸다.
• EPOS-LHC 대신 Sibyll2.3c를 사용할 경우 최적의 매개변수는 유사하지만 조성 deviance가 약 20 단위 증가하여 강입자 상호작용 모델에 민감함을 시사한다.
• 최적의 모델은 소스 내 상호작용에서 생성된 2차 중성자가 암반 이하에서 양성자 플럭스에 상당히 기여함을 시사하며, 이는 경량 주입에도 불구하고 부드러운 양성자 스펙트럼을 설명할 수 있다.