Beyond the Local Void: A data-driven search for the origins of the Amaterasu particle

기술 요약: 로컬 보이드(Local Void)를 넘어: 아마테라스스(Amaterasu) 입자의 기원을 찾기 위한 데이터 기반 탐색

문제 제기
초고에너지 우주선(UHECRs), 특히 $10^{20}$ eV를 초과하는 입자의 기원은 여전히 밝혀지지 않은 상태이다. 기존의 노력들은 UHECR을 스타버스트 은하(SBGs) 및 활동 은하핵(AGN) 목록과 연결해 왔으나, 이러한 모델들은 최고 에너지 수준의 사건들을 설명하는 데 어려움을 겪고 있다. 텔레스코프 어레이(Telescope Array)가 $244 \pm 29$ (통계적) $^{+51}_{-76}$ (계통적) EeV의 에너지로 검출한 "아마테라스스" 입자는 특정한 과제를 제시한다: 이 입자의 도착 방향은 알려진 활성 소스가 없는 저밀도 영역인 로컬 보이드(Local Void)를 가리키고 있다. 1차원 역추적(backtracking) 또는 단순 편향 추정치를 사용한 이전의 분석들은 이 입자가 알려진 활성 은하에서 기원하지 않았음을 시사하며, 이는 외계 물리학(예: 로렌츠 불변성 위반, 초중량 암흑 물질) 또는 일시적 소스(transient sources)에 대한 가설로 이어진다. 복잡한 전파 효과, 자기장 편향 및 재구성 불확실성을 고려하면서, 이러한 극한의 사건들과 일치하는 우주의 부피를 매핑하기 위한 엄격한 데이터 기반 방법론이 필요하다.

방법론
저자들은 개별 UHECR의 기원을 제약하기 위해 근사 베이지안 계산(Approximate Bayesian Computation, ABC)을 활용한 시뮬레이션 기반 추론 프레임워크를 도입한다. 이 방법은 명시적인 우도 함수(likelihood function)를 정의할 필요 없이, 3D 전파 모델링과 베이지안 파라미터 추정을 통합한다.

• 전파 모델링: 프레임워크는 3D 공간에서의 UHECR 전파를 시뮬레이션하기 위해 CRPropa 3를 사용한다. 여기에는 광피온 생성, 광분해, 전자-쌍 생성, 핵 붕괴 및 단열 손실을 포함한 모든 관련 상호작용이 포함된다.
• 자기장:
  • 은하 외 자기장 (EGMF): 25 kpc 격자 상에서 콜모고로프 난류 스펙트럼을 가진 가우시안 무작위 장으로 모델링된다. 파라미터로는 제곱평균제곱근 자기장 세기($B_{rms}$)와 상관 길이($L_c$)가 포함되며, 로그 균등 사전 분포(log-uniform priors)가 할당된다.
  • 은하 자기장 (GMF): 일관된 자기장을 위해 최적 적합된 UF23 기본 모델을, 난류 성분을 위해 JF12 플랑크 튜닝 성분을 사용하여 모델링된다.
• 파라미터 공간: 추론은 6개의 자유 파라미터인 소스 위치(은하 경도 $l$, 위도 $b$, 거리 $D_{src}$), 소스 에너지($E_{src}$), 그리고 EGMF 특성($B_{rms}$, $L_c$)을 대상으로 한다.
  • 소스 조성: 이전 연구와의 비교를 용이하게 하고 로컬 보이드를 벗어난 기원을 탐색하기 위해, 입자는 소스에서 철(iron) 핵이라고 가정한다.
  • 사전 분포 (Priors): $E_{src}$는 멱법칙 스펙트럼($E^{-1}$)을 따른다. 소스 거리는 명목 에너지의 경우 12 Mpc, 저에너지의 경우 15 Mpc까지의 구형 부피 내에서 균등하게 샘플링된다. 소스 방향 사전 분포는 예상되는 GMF 및 EGMF 편향을 고려하여, 철 핵의 역추적 방향을 중심으로 하는 본 미제스-피셔(von Mises-Fisher) 분포를 따른다.
• 추론 과정 (ABC):
  1. 사전 분포로부터 파라미터를 샘플링한다.
  2. 샘플링된 소스 위치에서 $10^6$개의 입자를 등방적으로 방출하여 지구까지 전파시킨다.
  3. 시뮬레이션된 사건을 관측된 아마테라스스 데이터(도착 방향 및 재구성된 에너지)와 거리 메트릭을 사용하여 비교한다.
  4. 시뮬레이션된 이벤트가 관측된 에너지의 $3\sigma$ 범위 내에 들어오면 파라미터 세트를 수락한다.
  5. 수락된 세트는 일치 정도에 따라 가중치가 부여되어 사후 분포(posterior distributions)를 생성한다.
  6. 두 가지 에너지 시나리오가 분석된다: 명목 케이스($E_{det} = 244$ EeV)와 낮은 계통적 하한값($E_{det} = 168$ EeV).

주요 결과
이 프레임워크를 아마테라스스 입자에 적용한 결과는 다음과 같다:

• 소스 후보: 사후 분포는 이전 분석에서 발견된 것보다 더 넓은 범위의 인근 소스 후보들을 드러낸다. 로컬 보이드가 여전히 관심 영역으로 남아 있으나, 결과는 로컬 그룹(Local Group) 또는 인근 은하 내의 소스 가능성을 배제하지 않는다.
• 질량 의존성: 수락된 시뮬레이션은 무거운 핵(질량수 $A > 28$)에 의해 지배된다. 명목 에너지 케이스의 경우, 수락된 이벤트의 77.9%가 무거운 조성을 나타내는 반면, 중간 질량은 16.2%, 가벼운 질량은 5.9%였다.
• 거리 제약:
  • 명목 에너지($244$ EeV)의 경우, 가장 확률이 높은 소스 거리는 $\sim 3$ Mpc이다. 무거운 핵의 경우 10–12 Mpc 사이에서 이차적인 피크가 나타난다.
  • 낮은 에너지 케이스($168$ EeV)의 경우, 가장 확률이 높은 거리는 가벼운/중간 질량에 대해 $\sim 4$ Mpc로 이동하며, 무거운 핵의 경우 더 긴 꼬리를 가진다.
  • 도착 방향은 가장 엄격한 선택 기준으로 작용하여, 큰 에너지 불확실성에도 불구하고 유효 거리 지평(distance horizon)을 제한한다.
• 편향 분석: 연구는 은하 자기장(GMF)이 가볍거나 중간 질량의 입자 편향을 지배하는 반면, 은하 외 자기장(EGEF)은 무거운 입자나 더 먼 거리, 혹은 더 낮은 에너지에서 유의미하게 기여함을 정량화한다.
• 강성 (Rigidity): 선호되는 소스 강성은 명목 에너지에 대해 $R \approx 10^{19.3}$ V이며, 낮은 에너지 케이스에 대해서는 $R \approx 10^{18.85}$ V로 나타난다.

의의 및 주장
본 논문은 "최고 에너지 영역의 우주선에 대한 시뮬레이션 기반 분석"을 위한 기초적인 프레임워크를 제공한다고 주장한다. 이 논문의 주요 기여는 다음과 같이 개별 UHECR 관측을 직접 해석 가능한 소스 제약으로 변환하는 능력에 있다:

1. 에너지와 도착 방향을 결합하여 제약함으로써, 그들의 상관관계와 모델링 불확실성을 자연스럽게 통합한다.
2. 자기장, 소스 조성 및 재구성 오류를 한계 제거(marginalizing over) 함으로써, 가능한 소스 위치에 대한 3D 사후 분포를 생성한다.
3. "로컬 보이드" 독점성에 대한 도전: 무거운 핵을 이용한 표준 전파 모델이 로컬 보이드 외부의 소스와도 일치하는 결과를 낼 수 있음을 입증함으로써, 저자들은 외계 물리학적 설명이 반드시 필수적인 것은 아니라는 점을 시사한다. 다만, 이를 완전히 배제하는 것은 아니다.

저자들은 현재의 구현이 계산적 제약(예: 특정 GMF 모델의 선택 및 고정된 철 조성)에 의해 제한된다는 점을 언급하며 연구 범위를 겸손하게 설정하였으며, 향후 연구에서는 더 넓은 범위의 GMF 모델을 포함하고 입자 유형을 자유 파라미터로 다루는 것을 목표로 한다고 밝혔다.