Hadronic Scenario for Galactic PeVatron LHAASO J1912+1014u Supported by Fermi-LAT $\gamma$-ray Data and FUGIN CO Data

기술 요약: 페르미-LAT 감마선 데이터와 FUGIN CO 데이터에 의해 지지된 은하계 PeVatron LHAASO J1912+1014u 의 강입자 시나리오

문제 제기
대형 고고도 공기 샤워 관측소 (LHAASO) 는 40 개 이상의 아페타 (sub-PeV) 감마선 천체를 확인했으며, 이는 페타전자볼트 (PeV) 에너지를 달성할 수 있는 우주선 가속기인 은하계 PeVatron 의 잠재적 후보들이다. 천체물리학의 중요한 과제는 이러한 천체들이 양성자 (강입자 PeVatron) 를 가속하는지, 아니면 전자/양전자 (경입자 PeVatron) 를 가속하는지 구분하는 것이다. LHAASO 는 고에너지 관측 범위를 제공하지만, 중간 정도의 각분해능 ($0.3^\circ$–$0.8^\circ$) 으로 인해 종종 방출 메커니즘을 확실히 식별하지 못한다. 특히 TeV 천체 HESS J1912+101 과 공간적으로 일치하는 특정 천체 LHAASO J1912+1014u 는 복잡한 사례로, 이전의 GeV 연구들은 배경 확산 방출과 천체의 형태를 모델링하는 데 어려움을 겪어 가속기의 성질이 해결되지 않은 상태였다.

방법론
저자들은 LHAASO J1912+1014u / HESS J1912+101 영역에 대해 다음과 같은 다중 파장 분석을 수행하였다:

1. 페르미-LAT 데이터: 0.4–409.6 GeV 에너지 범위에서 15 년간 (2008 년 8 월–2023 년 8 월) 관측된 데이터.
2. 성간매질 (ISM) 가스 템플릿: 고해상도 FUGIN CO ($^{12}$CO $J=1-0$) 데이터와 VLA 은하평면 탐사 (VGPS) HI 데이터. 저자들은 문헌에서 제안된 두 가지 특정 시스템 속도 범위, 즉 $V_{LSR} \approx 25$ km s$^{-1}$ (Sano 등 2018) 와 $V_{LSR} \approx 60$ km s$^{-1}$ (Su 등 2017) 에 대한 양성자 열밀도 ($N_p$) 지도를 구성하였다. 중요한 점은 광학적으로 두꺼운 HI 를 보정하여 정확한 열밀도를 유도했다는 것이다.
3. 배경 모델링: 1.6–12.8 GeV 대역에서 큰 잔차를 해결하기 위해 맞춤형 암흑 중성 매질 (DNM) 템플릿과 추가 점천원을 도입하여 표준 페르미-LAT 확산 방출 모델을 개선하였다.
4. 스펙트럼 에너지 분포 (SED) 모델링: Naima 프레임워크를 사용하여 페르미-LAT GeV 데이터와 H.E.S.S., HAWC, LHAASO 의 TeV 측정을 결합한 현상학적 피팅을 수행하였다. 전자 냉각 및 비냉각 경우를 모두 고려하여 전자 - 양성자 플럭스 비율 ($K_{ep}$) 과 가스 밀도 ($n_{gas}$) 매개변수를 스캔하여 경입자 (역콤프턴, 전자 제동복사) 대 강입자 (양성자 - 양성자 충돌) 시나리오를 검증하였다.
5. X 선 제약 조건: 아카이브 Chandra 데이터를 분석하여 X 선 대응체를 탐색하고 확산 싱크로트론 방출에 대한 상한치를 유도하여 경입자 시나리오를 제약하였다.

주요 결과

• GeV 초과 신호 검출: 배경 모델을 정교화함으로써 저자들은 LHAASO/H.E.S.S. 천체와 일치하는 10 GeV 이상의 통계적으로 유의미한 확장된 감마선 초과 신호를 검출하였다. 이 방출은 $\sim 2.1$의 광자 지수를 가진 경한 스펙트럼을 보인다.
• 공간 상관관계: GeV 초과 신호는 시스템 속도가 $\sim 25$ km s$^{-1}$ 및 $\sim 60$ km s$^{-1}$인 분자 구름에 해당하는 성간 가스 템플릿으로 잘 재현된다. 이러한 가스 템플릿을 사용한 피팅 품질은 단순 가우스 모델이나 H.E.S.S. 강도 지도보다 동등하거나 더 우수하다.
• SED 피팅 및 시나리오 평가:
  • 경입자 시나리오: 역콤프턴 (IC) 우세 모델은 인근 펄사 PSR J1913+1011 의 이용 가능한 감속 에너지 ($W_e \sim 10^{50}$ erg) 를 초과하는 과도한 에너지 예산 요구로 인해 배제된다. 전자 제동복사 (eB) 우세 모델은 에너지적으로 가능하지만, 전자 복사 냉각을 고려할 때 고에너지 컷오프를 재현하지 못한다. 필요한 컷오프 에너지 ($\sim 50$ TeV) 는 추정된 나이의 초신성 잔해 (SNR) 에서 전자가 달성할 수 있는 최대 에너지를 초과한다.
  • 강입자 시나리오: 강입자 시나리오 (양성자 - 양성자 충돌) 는 전자 냉각을 포함할 때 ( $E_{cut,e} = 30$ TeV 고정) GeV 에서 TeV 범위의 스펙트럼을 성공적으로 재현한다. 명목상의 경우 ($K_{ep} = 0.01$, $n_{gas} = 10$ cm$^{-3}$) 에서 강입자 성분은 10 TeV 이상에서 방출을 지배한다.
• 유도된 매개변수: 추정된 우주선 양성자 스펙트럼의 스펙트럼 지수는 $\sim 2.2$이다. 1 GeV 이상의 우주선 양성자 총 에너지는 가정된 거리와 가스 속도 연관성에 따라 $(1-5) \times 10^{49}$ erg 로 추정된다.
• X 선 제약 조건: 확산 X 선 방출은 검출되지 않았다. 유도된 95% 확산 X 선 플럭스 상한치는 엄격하다. 일반적인 성간 자기장 조건 ($6$ $\mu$G) 하에서 경입자 시나리오가 예측하는 싱크로트론 방출은 이 상한치를 초과하여 전자 가속 가설을 더욱 배제한다.

의의 및 주장
본 논문은 LHAASO J1912+1014u / HESS J1912+101 시스템이 가장 확률적으로 양성자 PeVatron이라고 주장한다. 저자들은 경한 GeV 스펙트럼, 성간 가스와의 공간 상관관계, TeV 스펙트럼 형태, 그리고 확산 X 선 방출의 부재 등 관측된 다중 파장 데이터에 대해 강입자 시나리오가 가장 일관된 설명을 제공한다고 논증한다.

본 연구는 아페타 천체에서 입자 종류를 식별하기 위해 성간매질 가스 분포와 페르미-LAT 데이터 간의 상세한 비교가 필수적임을 강조한다. 극단적인 가정 (예: 매우 낮은 자기장 또는 특정 펄사 풍선 성운 구성) 하에서는 경입자 시나리오가 공식적으로 배제될 수는 없으나, 에너지학, 냉각 제약 조건, 그리고 X 선 상한치에 의해 강입자 모델이 선호된다. 저자들은 이 천체가 무릎 (knee) 에너지까지의 우주선 양성자를 가속하는 은하계 가속기의 견고한 후보라고 결론지었다.