Spatially resolved spectral properties of M87* on event horizon scales

技術要約：イベントホライゾン・スケールにおけるM87*の空間分解スペクトル特性

問題と動機
M87の中心にある超大質量ブラックホール（M87*）は、一般相対性理論およびブラックホール天体物理学を検証するための主要な実験場となっている。これまでのイベントホライゾン・テレスコープ（EHT）およびGlobal Millimeter VLBI Array（GMVA）による観測では、1.3 mm（230 GHz）および3.5 mm（86 GHz）におけるM87*の形態と偏光のイメージングには成功しているが、事象の地平線付近における放射の空間分解されたスペクトル特性については、依然としてほとんど制約が付けられていない。地平線スケールにおけるスペクトル指数（$\alpha$、ここで $S_\nu \propto \nu^\alpha$）を特徴付けることは、プラズマの物理的条件（磁場、光学的厚さ、電子エネルギー分布）を診断し、光子リングのような周波数に依存しない成分を分離して、一般相対性理論に対するより厳格なテストを可能にするために極めて重要である。

手法
本研究では、2018年4月に取得されたM87*のほぼ同時期の高解像度画像（3.5 mm (86 GHz) のGM_VA画像（ALMAおよびグリーンランド・テレスコープを含む）および1.3 mm (230 GHz) のEHT画像）を解析している。

• データソース: 3.5 mm (86 GHz) 画像（ALMAおよびグリーンランド・テレスコープを含む）および1.3 mm (230 GHz) 画像（EHT）。
• 画像処理: ピクセルごとの比較可能性を確保するため、著者らは以下の3つの重要なステップを実行した。
  1. アライメント: 位相相互相関法を用いて、最も可能性の高いリング中心として定義されるブラックホール中心に画像を位置合わせした。両方の周波数における放射は降着流によって支配されているため、コアシフト効果は無視した。
  2. リサンプリング: 両方の画像を、共通の角解像度ピクセルスケールである0.2 $\mu$asにリサンプリングした。
  3. ブラーリング: 解像度に起因するスペクトルアーティファクトを軽減するため、両方の画像を共通の30 $\mu$asの円形ビームで畳み込み演算（コンボリューション）した。
• スペクトル指数の計算: スペクトル指数は、2つの周波数におけるフラックス密度を用いてピクセルごとに計算された。半径方向の傾向を特徴付け、方位角の非対称性や日単位の変動への感度を低減するため、著者らは以下の方法からもスペクトル指数を導出した。
  • 環状積分（Annular Integration）: ブラックホールを中心とする同心円状の環（幅2 $\mu$as）内に積分されたフラックス密度。
  • 方向性プロファイル（Directional Profiles）: ジェット軸（PA=288$^\circ$）およびそれに垂直な方向（PA=18$^\circ$）に沿ったストリップに沿って積分されたフラックス密度。
• 不確実性の推定: 系統的な不確実性（$\sim$10%）と熱雑音を考慮して誤差を推定し、典型的なスペクトル指数の不確実性を $\sigma_\alpha \sim 0.2$ とした。

主な結果
本解析は、M87*のコンパクトな領域（$\le$100 $\mu$as）内における初の空間分解スペクトル指数マップをもたらした。

• 形態: スペクトル指数のマップは、全強度画像と一致するリング状の形態を示している。内側領域（明るいリングの内側）は逆転した（正の）スペクトル指数を示し、リングの外側の領域は、次第に急峻な（負の）スペクトルを示す。
• 半径方向の勾配: $\sim$50 $\mu$as以内で堅牢な半径方向の勾配が検出された。
  • 内側領域 ($\lesssim 20 \mu$as): スペクトル指数は正（$\alpha \approx 1$）であり、これは光学的に厚いシンクロトロン放射を示している。$\sim$15 $\mu$as付近でわずかな増強が観察されている。
  • 遷移: スペクトル指数は半径とともに着実に減少し、$\sim$30 $\mu$as付近でゼロ（$\alpha \approx 0$）を通過する。これは3.5 mmのリングの特徴的な半径と一致している。
  • 外側領域: 3.5 mmのリングを超えると、スペクトルはさらに急峻になり、負の値となる。
• 方向の一致: ジェット軸および横方向のプロファイルは、全体的な半径方向の挙動を再現しており、この傾向が画像再構成のアーティファクトではなく、本質的なものであることを裏付けている。
• 積分スペクトル: コンパクトな領域（$<100 \mu$as）内の積分スペクトル指数は $\alpha \simeq -0.04$ であり、ほぼ平坦なスペクトルおよび従来の空間分解されていない多波長測定値と一致している。

意義と主張
本論文は、これらの結果がM87*における内側降着流およびジェット噴出の物理に関する新たな観測的制約を提供すると主張している。

• 物理的解釈: 観測された半径方向の進化（正から負への $\alpha$ の変化）は、内側領域を支配する光学的に厚いプラズマから、より大きな半径における次第に光学的に薄いシンクロトロン放射への遷移として解釈される。これは、周波数依存のシンクロトロン不透明度（高い周波数である1.3 mmの放射は、低い周波数である3.5 mmの放射よりもブラックホールに近い場所で光学的に薄くなる）と一致している。
• 理論的背景: 測定された半径方向の勾配（減少領域における $d\alpha/dr \sim -0.11 \mu$as$^{-1}$）は、一般相対論的磁気流体力学（GRMHD）シミュレーション（具体的には熱的電子SANEモデル）の予測に匹敵する。著者らは、理論モデルは光子リングにおいて鋭いスペクトルのスパイクを予測するが、現実的な観測解像度はこの特徴を平滑化するため、頑健な半径方向の減少が主要な検出可能なシグネチャとして残ると指摘している。
• 診断的有用性: スペクトル遷移の場所や半径方向の勾配を含む、空間分解されたスペクトル指数プロファイルは、降着シナリオ（例：MAD vs. SANE）を区別し、電子エネルギー分布を制約し、ブラックホールのスピンや磁場強度などの物理パラメータを推定するための強力な診断ツールとなる。
• 範囲: これらの知見は、M87*における地平線スケールでの系統的な半径方向のスペクトル遷移の最初の観測的証拠を提供し、地平線スケールの降着流の物理とパーセクスケールのジェット構造との間の溝を埋めるものである。