The study of an interaction between the jet and an interstellar medium around knot E and knot F of radio galaxy M87 by Chandra

기술 요약: M87 Knot E 및 F의 제트-ISM 상호작용에 대한 X선 분석

문제 정의
라디오 은하 M87의 상대론적 제트와 주변 성간 매질(ISM) 사이의 상호작용은 여전히 연구 대상이다. 2000년부터 2009년까지의 Chandra 아카이브 데이터를 병합하여 사용한 Dainotti 등(2012)의 선행 연구에서는, Knot E와 F 사이의 제트 외부 영역에서 "연엑스선 함몰(soft X-ray dip)" 현상을 보고했다. 이 특징은 제트가 ISM을 압축하여 연엑스선을 흡수하는 차갑고 밀도가 높은 구름을 생성한 결과로 해석되었다. 그러나 기존 분석에는 두 가지 잠재적 한계가 있었다: (1) 카운트율에서 "파일업(pile-up)" 이벤트(여러 광자가 CCD 픽셀에 동시에 충돌하여 에너지 및 이미지 데이터를 왜곡하는 현상)의 불완전한 제거, (2) 프리세션(precession) 및 속도로 인한 제트의 겉보기 운동을 고려하지 않아 장기 병합 이미지에서 특징이 흐려질 가능성이다. 또한, 압축된 ISM은 연엑스선을 흡수할 것으로 예상되지만, 제트 내부의 충격 가열(shock heating)은 이론적으로 열적 방출(thermal emission)을 생성할 것으로 기대된다. Knots E와 F의 스펙트럼에서 이러한 열적 방출의 존재 여부는 고통계수 데이터를 통해 확정적으로 확인되거나 배제되지 않은 상태였다.

방법론
저자 S. Osone은 이전의 발견을 검증하고 열적 방출을 탐색하기 위해 2000년부터 2018년까지의 Chandra 아카이브 데이터를 재분석하였다. 본 연구는 CIAO 4.15 및 CALDB 4.10.2를 사용하여 엄격한 데이터 선택 및 처리 파이프라인을 적용하였다:

1. 이미징을 위한 데이터 선택:

   • 파일업 완화: 데이터는 다음 세 가지 기준에 따라 선택되었다: (a) 핵(nucleus) 또는 HST-1에서 발생하는 가시적인 파일업 라인을 제거하기 위한 시각적 검사, (b) 핵과 HST-1의 공간적 분리 보장, (c) 밝은 Knot들의 파일업 비율이 5%(일부 데이터셋에서는 10%) 미만임을 보장하기 위한 count/frame 비율 계산. 밝은 Knot의 경우 포화 현상 때문에 3.1–3.2초의 프레임 타임을 가진 관측은 제외되었으며, 0.4초 프레임 타임이 활용되었다.
   • 프리세션 보정: 저자는 제트의 겉보기 운동을 고려하였다. 약 11년의 프리세션 주기와 Knot A 위치에서의 약 0.6년 주기(0.5" 이동에 해당)를 고려하여, 흐림 현상을 방지하기 위해 0.6년 단위의 윈도우 내 데이터를 사용하여 병합 이미지를 구축하였다.
   • 정렬: 라디오 코어(radio core)의 위치를 고정 접점(fixed tangent point)으로 사용하였다. X선 핵이 라디오 코어로부터 0.5" 이상 벗어난 관측은 제외되었다.

2. 이미지 처리:

   • 병합 이미지는 두 가지 에너지 대역(0.5–2 keV(연X선), 2–7 keV(경X선))에서 생성되었다.
   • 저자는 이전 연구에서 사용된 csmooth 도구의 사용을 명시적으로 피했는데, 이는 Virgo 성단 내의 뜨거운 가스와 제트 방출과 같은 복잡한 환경에서 해당 도구의 배경 추정이 결함이 있다고 판단했기 때문이다. 대신, 보정된 이미지에 가우시안 평활화(Gaussian smoothing, sigma=1, radius=2)를 적용하였다.
   • 특징의 유의성은 의심되는 함몰(dip) 영역의 카운트를 동일한 면적과 핵으로부터의 거리를 가진 남쪽의 대조 영역(control region)과 비교하여 계산되었다.

3. 스펙트럼 분석:

   • 희미한 열적 방출을 탐색하기 위해, 총 480 ks에 달하는 더 긴 노출 시간(frame time 3.1 sec)을 가진 별도의 데이터셋을 사용하였다.
   • Knot E와 F에 대한 파일업을 점검하였으며, 선택된 모든 데이터는 5% 미만의 파일업 오염을 보였다.
   • 배경 영역은 Virgo 성단의 뜨거운 가스의 반경 방향 의존성을 고려하여 소스 영역에 인접하게 선택되었다.
   • 스펙트럼은 광전 흡수(phabs)로 수정된 파워 로(power-law) 모델(싱크로트론 방출)을 사용하여 XSPEC으로 피팅하였다. 칼럼 밀도(column density)는 자유 변수로 두거나 21 cm 라디오 관측값($1.27 \times 10^{20} \text{ cm}^{-2}$)으로 고정하였다.

주요 결과

1. 연엑스선 함몰의 확인:

   • 0.5–2 keV 대역에서, Knot E와 F 사이의 제트 외부 영역에서 연엑스선 카운트의 유의미한 함몰이 확인되었다.
   • 가장 품질이 높은 데이터셋(Data 1, 224 ks)에서 이 함몰의 유의성은 12$\sigma$였다. 통계치가 낮은 데이터셋(Data 2 및 3)에서도 각각 6$\sigma$ 및 7$\sigma$의 유의성을 유지하였다.
   • 2–7 keV 대역에서는 이러한 함몰이 관찰되지 않았으며(유의성 1$\sigma$), 이는 해당 특징이 연엑스선에 특이적임을 확인해 준다.
   • 이 결과는 제트가 ISM을 압축하여 연엑스선을 흡수하는 차가운 분자 구름을 생성한다는 가설을 뒷받침한다.

2. Knot E 및 F의 스펙트럼 분석:

   • 충격 가열된 가스로부터의 열적 방출을 탐지하기 위해 Knot E와 F의 합산된 X선 에너지 스펙트럼을 분석하였다.
   • 스펙트럼은 광전 흡수를 동반한 단순 파워 로 모델(싱크로트론 방출)로 잘 설명되었다.
   • Knot E의 경우, 적합 확률(chance probability)은 0.0657이었다. Knot F의 경우, 칼럼 밀도를 21 cm 값으로 고정했을 때 적합 확률은 0.141이었다.
   • 데이터를 설명하기 위해 열적 방출 성분을 추가하는 것은 통계적으로 요구되지 않았다. 결과는 이 특정 Knot들에서 방출이 싱크로트론 복사에 의해 지배됨을 나타내며, 충격 가열된 가스로부터의 열적 방출에 대한 유의미한 증거는 발견되지 않았다.

의의 및 주장
본 논문은 Dainotti 등(2012)이 보고한 대로, Knot E와 F 사이의 연엑스선 함몰을 견고하게 확인하였다고 주장한다. 본 연구의 의의는 특히 이전 연구에서 충분히 다뤄지지 않았던 완전한 파일업 제거와 제트 프리세션 고려를 포함한 엄격한 데이터 선택 과정에 있다. 이러한 아티팩트(artifact)를 제거함으로써, 저자는 이 함몰이 제트와 ISM 사이의 실제적인 물리적 상호작용임을 입증하였다.

그러나 본 연구는 압축된 ISM(아마도 분자 구름)에 의한 연엑스선 흡수는 확인되었으나, 제트 내부의 충격 가열로 인해 기대되는 열적 방출은 이 특정 Knot들의 스펙트럼에서 탐지되지 않았다고 겸허히 결론짓는다. 데이터는 싱크로트론 전용 모델과 일치하며, 이는 만약 열적 방출이 존재한다면 현재의 분석 탐지 임계값 미만이거나, 혹은 이 특정 영역에서 충격 가열 메커니즘이 예상보다 효율이 낮음을 시사한다. 본 논문은 새로운 응용 분야나 향대의 실험 제안을 제시하기보다는, 개선된 아카이브 데이터 분석을 통해 제트-ISM 상호작용의 물리적 상태를 명확히 하는 데 집중하였다.