Study of an interaction between the jet and an interstellar medium of M87 with a spectral analysis by using Chandra

기술 요약: Chandra를 이용한 M87의 제트와 성간 매질 간 상호작용 연구 및 스펙트럼 분석

문제 제기
처녀자리 은하단 중심에 위치한 M87 라디오 은하는 20초각(arcseconds)에 달하는 돌출된 제트를 보여준다. M87의 TeV 감마선 기원은 일반적으로 핵(nucleus) 또는 HST-1과 같은 특정 매듭(knot)에서 기인하는 것으로 여겨지지만, 제트와 주변 성간 매질(ISM) 사이의 상호작용을 구동하는 물리적 메커니즘은 여전히 조사 중이다. Chandra 데이터를 사용한 이전의 M87 제트 스펙트럼 분석(예: Wilson & Yang 2002; Perlman & Wilson 2005; Sun et al. 2018)은 주로 싱크로트론 방출을 가정하여 X선 에너지 스펙트럼을 파워 로우(power-law) 모델로 피팅하였다. 그러나 이러한 연구들은 다음과 같은 한계에 직면했다:

1. 불충분한 통계량: 많은 분석이 짧은 노출 시간(예: 13–38 ks)에 의존하였으며, 이는 데이터 품질이 낮아 복잡한 모델이 단순히 낮은 데이터 품질 때문에 적합해 보일 수 있는 불충분한 통계량 문제를 야기했다.
2. 배경 차감 오류: 이전 연구들은 종종 제트 전체에 공통된 배경 영역을 사용했다. 이는 물리적으로 옳지 않은데, 왜냐하면 배경에 기여하는 처녀자리 은하단의 뜨거운 가스는 핵으로부터의 거리에 따라 강도가 변하기 때문이다.
3. 열적 성분의 무시: 제트가 충격파(shock)를 통해 ISM과 상호작용하여, 비열적 싱크로트론 방출 외에도 열적 성분(브레름슈트랄룽, bremsstrahlung)을 생성하는 가스를 가열할 가능성이 있다. 이러한 성분은 저통계 연구에서 누락되었을 수 있다.

방법론
이러한 문제를 해결하기 위해, 본 연구는 2000년부터 2018년까지의 관측을 포함하여 총 노출 시간 약 800 ks에 달하는 상당히 큰 규모의 Chandra X선 천문대 데이터셋을 활용한다.

• 데이터 선택 및 파일 업(Pile-up) 처리: 분석은 핵(nucleus), 매듭 HST-1, 매듭 D, 매듭 A의 네 가지 특정 영역에 집중한다. 핵과 HST-1의 높은 밝기 때문에, "파일 업"(여러 광자가 단일 픽셀에 부딪혀 스펙트럼을 왜곡하는 현상)을 관리하기 위해 데이터를 신중하게 선택했다. 대부분의 분석에서 파일 업 비율이 5% 미만이 되도록 데이터를 필터링했으나, 통계량을 극대화하기 위해 초기 데이터셋의 핵과 HST-1에서는 10%까지 허용되었다.
• 배경 모델링: 이전 연구들과 달리, 본 연구에서는 처녀자리 은하단의 뜨거운 가스의 공간적 변화를 고려하기 위해 각 매듭(소스의 북쪽 및 남쪽)에 대해 개별적인 배경 영역을 선택했다. 매듭 A의 경우, 알려진 연 X선 저하(soft X-ray dips)를 피하기 위해 특정 배경 영역을 선택했다.
• 스펙트럼 분석: 에너지 스펙트럼은 CIAO 소프트웨어(v4.15)를 사용하여 추출되었으며 XSPEC으로 분석되었다. 저에너지에서의 보정 문제를 피하기 위해 에너지 범위는 0.5 keV에서 10 keV로 설정되었다.
• 모델 피팅:
  1. 싱크로트론 방출을 나타내는 광전 흡제가 포함된 단순 파워 로우 모델을 테스트했다.
  2. 만약 파워 로우 피팅이 통계적으로 거부될 경우(높은 $\chi^2$ 및 소프트 엑세스로 표시됨), 파워 로우와 APEC 모델(열적 브레름슈트랄룽)의 결합 모델을 적용했다.
  3. 최적 피팅 모델로부터 중성 밀도, 플라즈마 밀도, 가속 전자 밀도를 포함한 물리적 파라미터를 도출했다.
  4. 가속된 전자와 ISM의 상호작용으로부터 발생하는 비열적 브레름슈트랄룽의 플럭스를 계산하고 이를 Fermi 감마선 관측 결과와 비교했다.

주요 결과

1. 스펙트럼 특성:

   • 핵 및 매듭 D: 핵과 매듭 D의 X선 에너지 스펙트럼은 모두 싱크로트론 방출과 일치하는 단순 파워 로우 모델로 잘 설명된다.
   • 매듭 HST-1 및 매듭 A: 이 영역들의 경우, 순수 파워 로우 모델은 통계적으로 거부된다.
     • HST-1: 결합 데이터셋(10% 파일 업 데이터 포함)에서 파워 로우 모델은 $2.6 \times 10^{-4}$의 우연 확률(chance probability)을 보인다. 열적 성분(APEC)을 추가하면 합리적인 피팅이 가능하다.
     • 매듭 A: 파워 로우 모델이 강력하게 거부된다 ($\chi^2 = 628.76$, 자유도 355, 확률 $1.6 \times 10^{-17}$). 파워 로우와 열적 브레름슈트랄룽의 결합 모델은 데이터를 잘 설명한다.

2. 열적 성분의 특성:

   • 매듭 HST-1과 A에서 검출된 열적 성분은 0.2 keV의 온도와 0.00의 금속 함량(사실상 제로)을 특징으로 한다.
   • 금속 함량이 주변 처녀자리 은하단의 뜨거운 가스(약 1 태양 함량)보다 현저히 낮은 것은, 가열된 가스가 주변 매질이 아니라 제트 자체 내부에서 기원함을 시사한다.
   • 매듭 A에서의 열적 플럭스 비율은 총 플럭스의 약 22%이다.

3. 물리적 파라미터:

   • 플라즈마 밀도: APEC 정규화로부터 도출된 매듭 A의 전자 밀도는 약 $11.0 \, \text{cm}^{-3}$이다.
   • 충격 가열(Shock Heating): 낮은 금속 함량과 주변 매질에 비해 높은 온도를 가진 열적 성분의 존재는, 제트 내부의 가스가 ISM과의 상호작용으로 인한 충격파에 의해 가열된다는 가설을 뒷받침한다.
   • 중성 밀도: X선으로 도출된 칼럼 밀도(column density)와 21cm 라디오 관측값의 차이로부터 계산되었으며, 이는 제트 주변에 차가운 물질이 존재함을 나타낸다.

4. 비열적 브레름슈트랄룽 및 Fermi 비교:

   • 가속된 전자가 ISM과 상호작용하며 발생하는 비열적 브레름슈트랄룽의 예상 플럭스를 계산했다.
   • 핵과 매듭 HST-1에서 계산된 플럭스는 Fermi에서 관측된 GeV 플럭스보다 수십 배 이상 낮다.
   • 결론: 핵이나 HST-1에서 관측된 GeV 감마선 플럭스에 대한 비열적 브레름슈트랄룽의 기여는 없다. 논문은 매듭 A의 자기장이 가정된 것(<10 $\mu$G)보다 현저히 낮다면 기여할 수도 있다고 언급했으나, 이는 추가 연구가 필요하다.

의의 및 주장
본 논문은 고통계 데이터(800 ks)와 올바른 공간 의존적 배경 차감을 활용함으로써, M87 제트의 매듭인 HST-1과 A에서 X선 스펙트럼의 열적 성분을 성공적으로 식별했다고 주장한다. 이 발견은 파워 로우 모델에만 의존했던 기존의 해석에 도전한다.

저자들은 다음과 같이 결론짓는다:

1. 열적 방출은 낮은 금속 함량과 높은 온도로 볼 때, 제트 내부의 가스가 ISM과의 충격 상호작용에 의해 가열된 것으로 보인다.
2. 비열적 브레름슈트랄룽 메커니즘(전자-ISM 상호작용)은 M87의 관측된 GeV 감마선 플럭스를 설명하지 못하며, 이는 TeV/GeV 방출의 기원이 다른 곳(예: 핵 또는 다른 입자 가속 메커니즘)에 있을 가능성을 강화한다.
3. 이전 연구들은 불충분한 통계량과 부적절한 배경 모델링으로 인해 이러한 열적 특징을 놓쳤을 수 있다.

이 연구는 새로운 실험 설정을 제안하기보다는, 상대론적 제트의 다파장 에너지 스펙트럼과 물리적 조건을 해결하기 위해 고품질의 아카이브 데이터 분석이 필수적임을 강조한다.