A Multiwavelength Interpretation of HESS J1857+026 Emission Using the Fermi-LAT, VERITAS, and HAWC Observatories

기술 요약: HESS J1857+026 방출에 대한 다파장 해석

문제 제기
HESS J1857+026은 에너지 넘치는 펄서 PSR J1856+0245와 일치하는, 확장된 미확인 초고에너지(VHE) $\gamma$-선원이다. 이 소스는 펄서 풍 성운(PWN)의 후보이지만, 그 다파장 특성은 여전히 복잡하다. 이전 관측들은 GeV와 TeV 방출 영역 사이의 공간적 불일치와 확인된 X선 대응체의 부재를 드러냈으며, 이는 지배적인 방출 메커니즘에 대한 의문을 제기한다. 구체적으로, 이 방출이 순수하게 렙토닉(PWN으로부터의 역 컴프턴 산란)인지, 순수하게 하드론(분자 가스와 상호작용하는 초신성 잔해로부터의 파이온 붕괴)인지, 혹은 두 가지의 조합인지가 불분명하다. 더욱이, 진화된 PWN의 즉각적인 환경 내에서의 입자 확산 특성은 일반적인 성간 매질(ISM)과 비교하여 여전히 잘 규명되지 않은 상태이다.

방법론
본 연구는 15년간의 Fermi-LAT 데이터(300 MeV–2 TeV), VERITAS 데이터(0.3–10 TeV), 그리고 새로운 HAWC 데이터(1–316 TeV)를 사용하여 HESS J1857+026에 대한 포괄적인 다기기 분석을 제시한다. 분석은 다음 단계들을 거쳐 진행된다:

1. 데이터 분석:

   • Fermi-LAT: Pass 8 SOURCE 클래스 데이터에 대해 빈(bin) 가능도 분석(binned likelihood analysis)을 수행하였다. 관심 영역(ROI)은 4FGL-DR3 카탈로그와 확산 템플릿을 사용하여 모델링되었다. 공간적 확장 테스트는 점원(point-source), 방사형 디스크, 그리고 방사형 가우시안 템플릿을 비교하였다.
   • VERITAS: Gammapy 프레임워크를 사용하여 30시간의 노출 데이터에 대해 3D 빈 최대 가능도 분석을 수행하였으며, 소스를 전력 법칙(power law)과 대칭 가우시안 공간 형상으로 모델링하였다.
   • HAWC: Pass 5 데이터셋(2860일)에 대해 순방향 폴딩(forward-folding) 최대 가능도 분석을 적용하였다. 연구에서는 여러 공간 템플릿(가우시안 vs. 확산 기반)과 스펙트럼 형상(전력 법칙, 로그-파라볼라, 지수 컷오프가 있는 전력 법칙)을 테스트하였으며, 최적의 모델을 선택하기 위해 베이지안 및 아카이케 정보 기준(BIC/AIC)을 사용하였다.

2. 복사 모델링:

   • 정적(Time-Independent) 모델: 순수 하드론(SNR) 및 렙토-하드론(PWN/SNR) 시나리오를 테스트하기 위해 NAIMA 패키지를 사용하여 광대역 스펙트럼 에너지 분포(SED)를 피팅하였다.
   • 동적(Time-Dependent) 모델: PWN의 팽창과 냉각을 시뮬레이션하기 위해 진화 모델(Gelfand et al. 2009와 유사)을 채택하였다. 마르코프 연쇄 몬테카를로(MCMC) 절차를 통해 시스템 연령, 자기장 강도, 입자 주입 스펙트럼과 같은 물리적 매개변수를 제약하였다.

3. 확산 분석:

   • 세 가지 기기 모두에서 방사형 표면 밝기 프로파일을 추출하였다. 이 프로파일들은 입자 확산 계수 $D(E)$와 ISM 대비 확산 억제 정도를 규명하기 위해 분석적 확산 기반 모델(식 16)로 피팅되었다.

주요 결과

• 스펙트럼 및 공간적 특성:

  • Fermi-LAT 데이터는 반지름이 $0.38^\circ \pm 0.04^\circ$인 방사형 가우시안에 의해 가장 잘 설명되는 확장된 소스를 보여주며, 스펙트럼 지수는 $\Gamma \approx 2.07$이다. 5 GeV 이하에서는 더 부드러운 지수(\Gamma \approx 2.63}), 5 GeV 이상에서는 더 단단한 지수($\Gamma \approx 1.99$)를 갖는 스펙트럼 브레이크가 관찰된다.
  • VERITAS는 $15.2\sigma$의 유의성으로 소스를 검출하였으며, 스펙트럼 지수 $\Gamma = 2.48 \pm 0.06$과 가우시안 폭 $0.236^\circ \pm 0.017^\circ$를 측정하였다.
  • HAWC 데이터는 소스가 $\sim 37$ TeV까지 확장됨을 보여준다. 최적의 공간 모델은 $1 \mu$G의 자기장에 대해 단순 가우시안보다 확산 템플릿을 선호하며, 컷오프 에너지 $E_c \approx 14$ TeV를 갖는다.
  • 소스의 크기는 에너지 증가에 따라 감소하며(GeV 에너지에서의 $0.38^\circ$에서 HAWC 에너지에서의 $0.26^\circ$까지), 이는 진화된 PWN와 일치하는 에너지 의존적 모폴로지를 나타낸다.

• 복사 기원:

  • 렙토-하드론 vs. 하드론: 저에너지($<10$ GeV) 성분에 대해서는 순수 하드론 모델을 완전히 배제할 수 없으나, 렙토-하드론 모델이 광대역 데이터에 대해 더 우수한 적합도를 제공한다. 검출된 SNR 껍질의 부재와 진화 모델에서 도출된 낮은 주변 가스 밀도($n_0 \lesssim 0.13$ cm$^{-3}$)는 TeV 방출에 대한 지배적인 하드론 기원을 부정한다.
  • PWN 진화: 진화 모델은 이 시스템이 연령 $\tau = [16, 21]$ kyr, 낮은 자기장 강도 $B = [0.4, 1.6] \mu$G를 가진 진화된 PWN임을 시사한다. 이러한 매개변수들은 (싱크로트론 냉각으로 인한) 희미한 X선 방출과 에너지 의존적 모폴로지를 설명한다. 또한, 낮은 자기장은 TeV 방출이 국부 광장(photon fields)에 의한 역 컴프턴 산란에 의해 지배된다는 점을 뒷받침한다.

• 입자 확산:

  • HESS J1857+026 주변의 확산 계수는 평균적인 은하 ISM에 비해 억제되어 있는 것으로 밝혀졌다. 50 TeV에서 확산 계수는 $D_{PWN} \sim 2.14 \times 10^{28}$ cm$^2$ s$^{-1}$ ($B=1 \mu$G일 때)이며, 이는 ISM 값($D_{ISM} \sim 1.2 \times 10^{30}$ cm$^2$ s$^{-1}$)보다 약 2개 차수(order) 낮다.
  • 억제는 리지디티 의존성(spectral index $\delta$)보다는 주로 정규화($D_0$)에서 나타나며, $\delta$는 Kolmogorov ($\delta=1/3$) 또는 Kraichnan ($\delta=1/2$) 난류 모델과 일치한다. 이는 소스 근처의 국부적으로 강화된 자기 난류 시나리오를 지지한다.

의의 및 주장
본 논문은 결합된 다파장 데이터가 진화된 PWN 인근의 입자 확산 특성을 규명하는 강력한 도구가 될 수 있음을 주장한다. 주요 결론은 HESS J1857+026이 PSR J1856+0245에 의해 구동되는 PWN일 가능성이 가장 높으며, 관측된 $\gamma$-선 방출이 상대론적 전자에 의한 역 컴프턴 산란에 의해 지배된다는 것이다.

연구는 저에너지 스펙트럼 성분($E < 10$ GeV)이 여전히 호스트 또는 인근 SNR으로부터의 하드론 기여를 포함하거나, 혹은 동일한 PWN 내의 별개 전자 집단을 포함할 수 있음을 강조한다. 그러나 고에너지 방출은 확고하게 렙토닉이다. 도출된 확산 억제 인자는 다른 TeV PWN(예: Geminga, Monogem)과 일치하며, 이는 일반적인 ISM에 비해 이들 소스 근처에서 입자 수송이 현저히 저해된다는 가설을 강화한다. 저자들은 중성미자 과잉이 낮은 유의성으로 보고되었으나, $\gamma$-선 영역에 비해 중성미자 방출 영역의 상한치가 매우 크다는 점은 잠재적인 중성미자 신호가 복잡하거나 관련 없는 기원을 가졌음을 시사한다고 언급하였다.