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

技术摘要：利用 Chandra 对 M87 喷流与星际介质相互作用的研究及光谱分析

问题陈述
位于室女座星系团中心的 M87 射电星系展现出一个延伸 20 角秒的显著喷流。虽然 M87 的 TeV $\gamma$ 射线起源通常归因于核心或特定的结节（如 HST-1），但驱动喷流与周围星际介质（ISM）相互作用的物理机制仍处于研究之中。以往使用 Chandra 数据对 M87 喷流进行的谱分析（例如 Wilson & Yang 2002; Perlman & Wilson 2005; Sun et al. 2018）主要通过拟合幂律模型来描述 X 射线能量谱，并假设其为同步辐射。然而，这些研究面临以下局限性：

1. 统计量不足： 许多分析依赖于较短的曝光时间（例如 13–38 ks），导致统计特性较差，使得复杂的模型可能仅仅因为数据质量低而显得符合预期。
2. 背景减除误差： 先前的研究通常对整个喷流使用统一的背景区域。这在物理上是不正确的，因为贡献于背景的室女座星系团热气体强度会随距离核心的距离而变化。
3. 忽略热成分： 存在一种可能性，即喷流通过冲击波与 ISM 相互作用，从而产生热成分（轫科效应/热制动辐射），除了非热同步辐射之外。这种成分在低统计量的研究中可能被忽略了。

方法论
为了解决这些问题，本研究利用了来自 Chandra X 射线天文台的显著更大的数据集，涵盖了 2000 年至 2018 年间的观测，总曝光时间约为 800 ks。

• 数据选择与堆积处理： 分析重点关注四个特定区域：核心、结节 HST-1、结节 D 和结节 A。由于核心和 HST-1 的亮度极高，数据经过仔细筛选以管理“堆积”（pile-up，即多个光子撞击单个像素，导致光谱失真）。数据经过过滤，以确保大多数分析中的堆积比例低于 5%，但在早期数据集中，为了最大化统计量，核心和 HST-1 的比例被接受到 10%。
• 背景建模： 与以往研究不同，每个结节都选择了单独的背景区域（位于源的北侧和南侧），以解释室女座星系团热气体的空间变化。对于结节 A，选择了一个特定的背景区域以避开已知的软 X 射线凹陷。
• 光谱分析： 使用 CIAO 软件（v4.15）提取能量谱，并使用 XSPEC 进行分析。能量范围设定在 0.5 keV 至 10 keV 之间，以避免低能端的定标问题。
• 模型拟合：
  1. 测试了带有光电吸收的简单幂律模型（代表同步辐射）。
  2. 如果幂律拟合在统计上被拒绝（表现为高 $\chi^2$ 和软过剩），则应用幂律与 APEC 模型（热制动辐射）的组合模型。
  3. 从最佳拟合模型中推导出物理参数，包括中性密度、等离子体密度和加速电子密度。
  4. 计算了由加速电子与 ISM 相互作用产生的非热制动辐射通量，并将其与 Fermi $\gamma$ 射线观测结果进行比较。

关键结果

1. 光谱特征：

   • 核心与结节 D： 核心和结节 D 的 X 射线能量谱均能被简单的幂律模型很好地描述，这与同步辐射一致。
   • 结节 HST-1 与结节 A： 对于这些区域，纯幂律模型在统计上被拒绝。
     • HST-1： 在组合数据集（包含 10% 堆积数据）中，幂律模型给出的概率为 $2.6 \times 10^{-4}$。添加热成分（APEC）可提供合理的拟合。
     • 结节 A： 幂律模型被强烈拒绝（对于 355 个自由度，$\chi^2 = 628.76$，概率为 $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}$。
   • 冲击加热： 存在一个相对于环境介质具有低金属丰度和高温度的热成分，这支持了喷流内的气体通过与 ISM 相互作用产生的冲击波被加热的假设。
   • 中性密度： 通过 X 射线推导的柱密度与 21 cm 射电观测值的差异计算得出，表明喷流周围存在冷物质。

4. 非热制动辐射与 Fermi 比较：

   • 研究计算了加速电子与 ISM 相互作用产生的预期非热制动辐射通量。
   • 计算出的来自核心和结节 HST-1 的通量比 Fermi 观测到的 GeV 通量低几个数量级。
   • 结论： 核心或 HST-1 中的非热制动辐射对观测到的 Fermi 通量没有贡献。论文指出，如果结节 A 的磁场显著低于假设值（<10 $\mu$G），它可能会有贡献，但这需要进一步研究。

意义与主张
本文声称，通过利用高统计量数据（800 ks）和正确的空间相关背景减除，它成功地在 M87 喷流结节 HST-1 和 A 的 X 射线光谱中识别出了热成分。这一发现挑战了以往仅依赖幂律模型的解释。

作者得出结论：

1. 热辐射很可能是喷流内的气体通过与 ISM 的冲击相互作用而被加热的结果，其低金属丰度和高温度证明了这一点。
2. 非热制动辐射机制（电子-ISM 相互作用）并不解释观测到的 M87 GeV $\gamma$ 射线通量，这强化了 TeV/GeV 发射起源可能在于其他地方（例如核心或不同的粒子加速机制）的观点。
3. 先前的研究可能由于统计量不足和不恰当的背景建模而遗漏了这些热特征。

该研究并未提出新的实验装置，而是强调了利用高质量存档数据进行分析对于解析相对论性喷流的多波段能量谱和物理条件的必要性。