Modelling Galactic neutrino emission: contributions from massive star clusters and interstellar cosmic rays

기술 요약: 거대 항성 집단 및 성간 우주선으로부터의 은하 중성미자 방출 모델링

문제 제기
우주에서 가장 에너지가 높은 우주선(CR), 특히 PeV 에너지에 도달하는 PeVatron의 기원은 여전히 불확실하다. 초신성 잔해(SNR)가 전통적으로 주요 가속기로 고려되어 왔으나, $\gamma$-선 관측 결과 이들이 PeV 에너지에 도達하는 데 어려움이 있음을 시사하며, 이에 따라 대안적인 소스 클래스의 필요성이 제기되고 있다. 젊은 거대 항성 집단(YMSC)은 충돌하는 항성풍과 다수의 SNR 이벤트로 인해 생성된 강화된 자기 난류 덕분에 유망한 PeVatron 후보로 제안되어 왔다. 동시에, IceCube 관측소는 은하 평면으로부터 고에너지 중성미자 플럭스를 검출하였다. 그러나 현재의 분석은 성간 가스와의 상호작용으로 생성되는 순수 확산 중성미자 방출만을 고려한 템플릿에 의존하고 있으며, 이는 YMSC 인구와 같은 미해결 하드론 가속기의 기여를 명시적으로 해결하지 못하여 이러한 소스들의 특정 시그니처를 가릴 수 있다.

방법론
저자들은 두 가지 구별되는 성분을 통합하여 은하로부터 예상되는 중성미자 방출을 계산하기 위한 새로운 이론적 프레임워크를 제시한다:

1. 합성 YMSC 인구: 로컬 태양 근방($\lesssim 2$ kpc)에서만 YMSC를 관측할 수 있는 한계를 극복하기 위해, 저자들은 은하 YMSC 인구에 대한 100개의 독립적인 실현(realization)을 생성하였다. 이 합성 인구는 관측적 제약 조건에 기반하여 집단 질량($10^3 - 6.3 \times 10^4 M_\odot$), 연령($<30$ Myr), 그리고 공간 분포(나선 팔, 은하 바, Local Spur 추종)를 샘플링한다.
2. 입자 가속 및 수송: 각 합성 클러스터에 대해, 저자들은 항성 인구(주계열성, Wolf-Rayet, 초신성 조상 별)를 시뮬레이션하여 바람과 초신성으로부터 발생하는 집단 역학 광도를 계산한다. 저자들은 결과적으로 발생하는 풍직 블로잉 버블(wind-blown bubbles)을 모델링하여 가스 밀도와 자기장 프로파일을 계산한다.
   • 가속 메커니나즘: 입자 가속은 두 지점, 즉 집단 항성풍 종단 충격파(WTS)와 초신성 잔해(SNR) 충격파에서 모델링된다.
   • 확산 시나리오: 클러스터 버블 내에서의 CR 전파는 세 가지 난류 가정(Kolmogorov, Kraichnan, Bohm 확산) 하에서 계산된다.
   • 중성미자 생성: 중성미자 플럭스는 가속된 양성자와 대상 가스(주로 버블의 밀도가 높은 휩쓸린 껍질 부분) 사이의 하드론 충돌로부터 계산된다.
3. 확산 은하 성분: 모델은 은하 CR 해(sea)와 성간 가스의 상호작용으로부터 발생하는 확산 중성미자 방출 또한 계산한다. 이를 위해 CR 전파 모델(내부 은하에서의 스펙트럼 경화 포함)과 가스 지도(GALPROP 및 대안적인 3D 모델)를 활용한다.
4. 템플릿 생성: 총 예측 플럭스는 YMSC 기여와 CR 유도 확산 성분의 합이다. 이들은 IceCube 데이터와 비교하기 위해 공간적 및 스펙트럼 템플릿으로 합성된다.

주요 기여

• 새로운 중성미자 템플릿: 본 논문은 표준 CR 유도 확산 성분과 함께 YMSC 인구의 방출률을 명시적으로 포함하는 포괄적인 중성위 템플릿을 최초로 제공한다.
• 합성 인구 모델링: 근처의 몇몇 해상된 클러스터에만 의존하는 대신, 전체 은하 YMSC 인구로부터의 누적 중성미자 출력을 통계적으로 견고하게 추정한다.
• 다성분 분석: 클러스터 내부의 "in-situ" 중성미자 생성과 배경 확산 방출을 구분함으로써, 소스 기여를 분리할 수 있게 한다.
• 매개변수 민감도: 입자 확산 체제(Kolmogorov, Kraichnan, Bohm)와 은하 내부의 CR 스펙트럼 경화에 따른 예측 플럭스의 의존성을 탐구한다.

결과

• 개별 소스의 검출 가능성: 개별 YMSC로부터의 예측 플럭스는 세제곱 킬로미터 규모의 검출기(예: IceCube, KM3NeT)로도 간신히 검출 가능한 수준인 것으로 나타났다. 결과적으로, YMSC 인구의 누적 방출은 미해결 확산 성분으로 나타날 것으로 예상된다.
• 플럭스 크기: 모델의 정규화는 기존 Ice-Cube 템플릿의 최적 적합값과 일치한다. 저자들은 클러스터 방출으로부터의 미해결 기여가 무시할 수 없는 수준일 수 있음을 발견했다.
• 확산 의존성:
  • Bohm 확산 시나리오(가장 효율적인 가둠)에서, YMSCs는 약 $\sim 50$ TeV 미만에서 중성미자 플럭스를 지배하며, 만약 스펙트럼 경화가 가정되지 않는다면 Kraichnan 및 Kolmogorov 시나리오에 비해 플럭스를 과대 예측할 가능성이 있다.
  • Kolmogorov 및 Kraichnan 시나리오에서는, CR 유도 확산 성분이 일반적으로 YMSC 성분에 비해 우세하거나 대등한 기여를 하며, 특히 더 높은 에너지에서 그러하다.
• 공간 분포: YMSC 기여는 은하 경도 범위($|l| \le 30^\circ$, 은하 능선)에 집중되어 있는 반면, CR 유도 방출은 더 큰 경도에서 우세하다.
• 스펙트럼 경화: 내부 은하에서의 CR 해의 스펙트럼 경화를 포함하면 Kolmogorov 및 Kraichan 시나리오에서 YMSCs의 상대적 기여가 감소하여, 총 예측 스펙트럼이 $\sim 100$ TeV에서 IceCube 관측 결과와 더 잘 일치하게 된다.
• 일관성: 모델의 예측은 TeV-PeV 범위의 현재 확산 $\gamma$-선 관측과 일치한다.

의의 및 주장
저자들은 자신들의 연구가 은하 중성미자 방출의 이론적 모델링에 있어 필수적인 진화를 제공한다고 주장한다. 순수 확산 템플릿을 넘어 YMSC 인구의 특정한 방출률을 포함함으로써, 이들은 은하 중성미자 방출의 두 성분을 "분리"할 수 있는 도구를 제공한다. 본 논문은 이러한 새로운 템플릿을 미래의 중성미자 데이터셋에 적용함으로써, 극한의 CR 가속 및 중성미자 생성에서 항성 집단의 역할을 규명할 수 있을 것이라고 상정한다. 저자들은 연구 결과가 클러스스로부터의 미해결 기여가 무시할 수 없는 수준임을 시사하며, 소스 기여와 확산 성분의 공동 피팅(joint fitting)이 은하 가속기를 식별하고 우주선의 기원을 해결하는 데 매우 중요하다고 강조한다. 이 작업은 YMSC를 IceCube 신호의 원인으로 확정적으로 식별했다고 주장하는 것이 아니라, 이 가설을 관측 데이터에 대해 테스트할 수 있는 이론적 프레임워크를 제공하는 것이다.