High-energy Multi-messenger Emission from Galaxy Clusters in the Local Universe

기술적 요약: 국부 우주 은하단에서의 고에너지 다중신호 방출

문제 제기
확산 고에너지 중성미자와 감마선 배경의 기원은 천체물리학에서 여전히 해결되지 않은 주요 문제입니다. 활동은하핵 (AGN) 과 항성폭발 은하가 기여하는 것으로 알려져 있지만, 이러한 우주 메신저의 원천으로서 은하단의 구체적인 역할은 완전히 규명되지 않았습니다. 은하단은 뜨거운 은하단 내부 매질 (ICM) 가스 ($\sim 10^8$ K) 와 강력한 자기장 ($\sim 1$ $\mu$G) 을 특징으로 하며, 이론적으로 병합 충격파, 강착 충격파, 난류를 통해 우주선 (CR) 을 가속할 수 있습니다. 이러한 가속된 양성자는 강입자 과정을 통해 ICM 과 상호작용하여 2 차 감마선과 중성미자를 생성할 수 있습니다. 그러나 이전 연구들은 종종 단순화된 구형 대칭 모델에 의존하거나 중앙 은하단 코어에서만 CR 주입이 일어난다고 가정하여, 확산 방출의 공간적 분포와 총 플럭스를 잘못 표현했을 가능성이 있습니다. 또한, 현재 관측의 제한된 각분해능으로 인해 중앙 AGN(예: NGC 1275 또는 M87) 과 같은 점원 방출과 확산 은하단 방출을 구별하는 것은 여전히 어렵습니다.

방법론
본 연구는 국부 우주 ($z \lesssim 0.1$) 의 은하단에서 발생하는 다중신호 방출을 모델링하기 위해 고해상도 자기유체역학 (MHD) 시뮬레이션과 몬테카를로 수송 방법을 결합한 엄격한 수치 프레임워크를 사용합니다.

1. MHD 시뮬레이션 (SLOW): 본 연구는 CosmicFLow 카탈로그에 기반한 관측적으로 제약된 초기 조건을 포함하는 "Simulating the LOcal Web"(SLOW) 시뮬레이션을 활용합니다. 이전 모델들과 달리 SLOW 는 특히 처녀자리, 페르세우스자리, 코마자리 은하단 내에서 가스 밀도, 온도, 자기장 분포에 대한 자기 일관적이고 고해상도의 3 차원 설명을 제공하며, 병합으로 인한 비구형 형태와 하위 구조를 포착합니다.
2. 우주선 주입: 물리적으로 동기화된 주입 모델을 채택하여, CR 주입이 은하단 가스 밀도의 공간적 분포를 따르도록 설정했습니다. 이 접근법은 은하단 전체 부피에 걸친 난류와 충격파와 관련된 고밀도 영역에서의 CR 가속을 고려하며, 주입을 은하단 코어로만 제한하지 않습니다. 총 CR 에너지 예산은 관측적 상한선과 일치하는 비율 $X_{CR} = E_{CR}/E_{cluster} \simeq 0.01$로 매개변수화됩니다.
3. CR 수송 및 2 차 생성: ICM 과 은하간 매질 (IGM) 내 CR(양성자로 가정) 의 전파는 CRPropa 몬테카를로 코드를 사용하여 시뮬레이션됩니다. 이 코드는 대류, 공간 확산, 운동량 확산, 에너지 손실을 고려하여 파커 수송 방정식을 풉니다. 주요 상호작용 채널로는 $10^{17}$ eV 이하에서 지배적인 양성자 - 양성자 ($pp$) 충돌, 광파이온 생성, 그리고 전자기 캐스케이드 (CMB 및 EBL 에 대한 역콤프턴 산란, 싱크로트론 손실) 가 포함됩니다. 2 차 감마선과 중성미자는 이러한 상호작용으로부터 자기 일관적으로 생성됩니다.

주요 기여

• 현실적인 환경 모델링: SLOW MHD 시뮬레이션의 통합을 통해 이전의 분석적 또는 준분석적 연구들의 단순화된 구형 대칭 가정을 넘어, 현실적인 비구형 은하단 환경에서의 CR 전파 모델링이 가능해졌습니다.
• 확장된 주입 체계: 가스 밀도 프로파일을 추적하는 공간적으로 확장된 CR 주입 모델의 채택은 중앙만 주입한다고 가정한 이전 작업들에 비해 중요한 방법론적 개선을 나타냅니다. 이는 코어 내 유효 CR 가둠 시간과 상호작용률을 감소시켜 더 확산된 방출 프로파일을 유도합니다.
• 종합적 플럭스 추정: 본 연구는 개별 은하단 (처녀자리, 페르세우스자리, 코마자리) 의 플럭스와 500 Mpc 공변 거리 내의 전체 거대 은하단 집단 ($M \gtrsim 5 \times 10^{13} M_\odot$) 으로부터의 누적 확산 플럭스를 모두 계산했습니다.

결과

• 개별 은하단: 처녀자리, 페르세우스자리, 코마자리 은하단에 대한 예측된 확산 감마선 플럭스는 LHAASO 협회가 보고한 현재 상한선보다 훨씬 낮습니다. 예측된 플럭스는 또한 페르세우스자리 중앙원천 NGC 1275 에서 MAGIC 이 관측한 감마선 방출보다 낮으며, MAGIC 이 점원을 감지하는 반면 본 모델은 확산 방출을 예측한다는 기대와 일치합니다.
• 누적 플럭스: 500 Mpc 내 모든 은하단으로부터의 통합 감마선 플럭스는 Fermi-LAT 가 측정한 확산 감마선 배경보다 현저히 낮습니다. 연구는 $\sim 100$ GeV 미만의 에너지에서 다른 원천들 (블레이자, 정렬되지 않은 AGN, 항성폭발 은하) 이 배경을 지배한다고 지적합니다.
• 중성미자 방출: 국부 은하단 집단으로부터 예측된 모든 맛깔 중성미자 플럭스는 IceCube 가 관측한 확산 플럭스보다 훨씬 낮지만, 특히 은하단 ($M > 10^{14} M_\odot$) 에 대한 기존 IceCube 상한선과 비교 가능합니다. 처녀자리 은하단으로부터 예측된 중성미자 플럭스는 KM3NeT/ARCA 검출기의 예상 감도와 비교 가능한 것으로 나타났습니다.
• 스펙트럼 특성: 예측된 플럭스의 스펙트럼 지수 ($\alpha \approx 2.0 - 2.5$) 는 병합 및 강착 충격파에 의한 가속 기대치와 일치합니다. 결과는 질량이 낮은 은하단이 $\sim 10$ GeV 감마선 방출에 효율적으로 기여하는 반면, 거대 시스템은 우수한 CR 가둠 능력으로 인해 고에너지 ($\gtrsim 10$ TeV) 중성미자 플럭스를 지배함을 나타냅니다.

의의 및 주장
본 논문은 그 결과가 LHAASO 에 의한 유의미한 확산 방출의 비검출과 IceCube 에 의해 설정된 상한선을 포함한 현재 관측적 제약과 일치한다고 주장합니다. 저자들은 종종 더 높은 플럭스를 가정했던 이전 모델 예측과 IceCube 상한선 사이의 긴장 관계가 공간적으로 확장된 주입 모델의 채택과 고해상도 SLOW 시뮬레이션의 사용으로 인해 현저히 완화된다고 제시합니다.

본 연구는 은하단이 전체 확산 고에너지 중성미자 배경을 지배할 가능성은 낮지만, 그 중 무시할 수 없는 부분을 차지할 수 있다고 결론지었습니다. 이 작업은 KM3NeT 와 IceCube-Gen2 를 이용한 처녀자리 은하단으로부터의 확산 중성미자 탐색과 CTA 및 LHAASO 를 이용한 확산 감마선 탐색을 독려하는 이론적 프레임워크를 제공합니다. 저자들은 그들의 발견이 차세대 다중신호 관측소의 전망과 부합하며, 대규모 구조 내 강입자 CR 가속 과정에 대한 시기적절한 검증을 제공한다고 강조합니다.