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

技術要約：Chandraを用いたM87のジェットと星間物質（ISM）の相互作用に関する研究：スペクトル解析による検討

問題提起
処女座銀河団の中心に位置する電波銀河M87は、20秒角にわたって延びる顕著なジェットを有している。M87におけるTeVガンマ線の起源は、一般に核またはHST-1のような特定の結節点（knot）に由来するとされているが、ジェットと周囲の星間物質（ISM）との間の相互作用を駆動する物理的メカニズムについては、依然として調査が進められている段階である。Chandraのデータを用いたM87ジェットの過去のスペクトル解析（例：Wilson & Yang 2002; Perlman & Wilson 2005; Sun et al. 2018）では、主にシンクロトロン放射を仮定した冪関数（パワーロー）モデルを用いてX線エネルギー・スペクトルをフィッティングしてきた。しかし、これらの研究には以下の限界があった：

1. 統計的不足： 多くの解析は短い露光時間（例：13〜38 ks）に依存しており、データの質が低いために、複雑なモデルが単にデータの統計量の乏しさゆえに適合しているように見えてしまうという問題があった。
2. 背景放射減算の誤り： 従来の解析では、ジェット全体に対して共通の背景領域を使用していた。これは、M87のジェットと周囲の高温ガス（Virgo銀河団のガス）の輝度は核からの距離に応じて変化するため、物理的に不正確である。
3. 熱的成分の無視： ジェットが衝撃波を介してISMと相互作用し、非熱的なシンクロトロン放射に加えて、熱的成分（制動放射／熱的ブレムストラールング）を生じさせる可能性がある。低統計の研究では、この成分が見逃されていた可能性がある。

手法
これらの問題に対処するため、本研究では2000年から2018年までの観測を網羅する、総露光時間約800 ksという大幅に増大したChandra X線天文台のデータセットを利用する。

• データ選択およびパイルアップ処理： 解析は、核（nucleus）、結節点HST-1、結節点D、および結節点Aの4つの特定の領域に焦点を当てる。核およびHST-1の明るさが非常に高いため、パイルアップ（複数の光子が単一のピクセルに衝突し、スペクトルを歪ませる現象）を管理するためにデータを慎重に選択した。ほとんどの解析においてパイルアップの割合が5%未満になるようデータをフィルタリングしたが、統計量を最大化するために、初期のデータセットでは核とHST-1に対して10%まで許容した。
• 背景モデル化： 以前の研究とは異なり、Virgo銀河団の高温ガスの空間的な変化を考慮するため、各結節点のそれぞれに対して個別の背景領域を選択した。結節点Aについては、既知の軟X線ディップを避けるために特定の背景領域を選択した。
• スペクトル解析： エネルギー・スペクトルは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}$の偶然確率を示す。熱的成分（APEC）を加えることで、妥当なフィッティングが得られる。
     • 結節点A： 冪関数モデルは強く棄却される（自由度355に対し$\chi^2 = 628.76$、確率$1.6 \times 10^{-17}$）。冪関数と熱的制動放射の組み合わせモデルは、データを良好に説明する。

2. 熱的成分の特性：

   • 結節点HST-1およびAで検出された熱的成分は、温度0.2 keV、金属量0.00（実質的にゼロ）という特徴を持つ。
   • 金属量は周囲のVirgo銀河団の高温ガス（約1太陽組成）よりも著しく低く、これは加熱されたガスが周囲の媒体ではなく、ジェット自体に由来することを示唆している。
   • 結節点Aにおける熱的フラックス比は、全フラックスの約22%である。

3. 物理パラメータ：

   • プラズマ密度： APECの正規化量から導出された結節点Aにおける電子密度は約$11.0 \, \text{cm}^{-3}$である。
   • 衝撃加熱： 低金属量かつ周囲の媒体に対して高温である熱的成分の存在は、ジェット内のガスがISMとの相互作用による衝撃波によって加熱されているという仮説を支持している。
   • 中性密度： X線から導出されたカラム密度と21cm電波観測との差から算出され、ジェットの周囲に冷たい物質が存在することを示している。

4. 非熱的制動放射およびFermiとの比較：

   • 加速された電子とISMの相互作用による非熱的制動放射の期待フラックスを計算した。
   • 核および結節点HST-1から計算されたフラックスは、Fermiで観測されたGeVフラックスよりも数桁低い。
   • 結論： 核またはHST-1における非熱的制動放射は、観測されたFermiのフラックスへの寄与はない。論文では、もし結節点Aの磁場が想定よりも著しく低い（<10 $\mu$G）ならば、結節点Aが寄与している可能性があると述べているが、これにはさらなる研究が必要である。

意義と主張
本論文は、高統計データ（800 ks）を利用し、空間依存的な正しい背景減算を行うことで、M87ジェットの結節点HST-1およびAにおいて熱的成分を特定することに成功したと主張している。この発見は、冪関数モデルのみに依存していた従来の解釈に疑問を投げかけるものである。

著者らは以下のように結論付けている：

1. 熱的放射は、低金属量と高温によって裏付けられるように、ジェット内のガスがISMとの衝撃相互作用によって加熱されたものである可能性が高い。
2. 非熱的制動放射メカニズム（電子とISMの相互作用）は、M87から観測されるGeVガンマ線フラックスを説明できず、TeV/GeV放射の起源が他の場所（例：核または異なる粒子加速メカニズム）にあることを補強している。
3. 従来の研究は、統計不足および不適切な背景モデリングのために、これらの熱的特徴を見逃していた可能性がある。

本研究は新しい実験装置を提案するものではなく、マルチ波長エネルギー・スペクトルと相対論的ジェットの物理的条件を解明するために、高品質なアーカイブ・データの解析がいかに不可欠であるかを強調している。