Impacts of radiative cooling on the images of a black hole shadow and extended jets in two-temperature GRMHD simulations

技術的概要：二温度 GRMHD シミュレーションにおける放射冷却がブラックホールシャドウおよび拡張ジェット画像に与える影響

問題提起
最近のイベントホライズン望遠鏡（EHT）による M87* ブラックホールシャドウの観測は、高温降着流における放射効率を示しており、放射冷却が無視できないことを示唆している。一般相対性磁気流体力学（GRMHD）シミュレーションは降着物理をモデル化する標準的なツールであるが、多くの既存モデルは放射冷却を無視するか、パラメトリックな電子温度 prescriptions（例えば「R-β」モデル）を用いた単一流体近似に依存している。さらに、過去の研究ではしばしば熱的電子分布関数（eDF）が仮定されており、磁気リコネクションや乱流によって駆動される非熱的プロセスが見過ごされている可能性がある。本研究は、放射冷却が物理的に駆動された二温度モデルおよび非熱的 eDF と結合された場合、電子熱力学、ディスク構造、およびその結果生じるブラックホールシャドウと拡張ジェットの数値画像をどのように変化させるかという理解のギャップに対処する。

手法
著者らは、BHAC コードを用いて、回転ブラックホール（$a = 0.9375$）への磁化された降着流の一般相対性磁気流体力学（GRMHD）シミュレーションを実施した。シミュレーションでは、電子熱力学をイオンとは別に進化させる二温度枠組みを採用し、以下の要素を組み込んだ：

• 加熱メカニズム： 二つの異なる prescriptions がテストされた。乱流加熱（Kawazura et al. 2019）と磁気リコネクション加熱（Rowan et al. 2017）である。
• 冷却メカニズム： 放射冷却（制動放射、サイクロトロン・シンクロトロン、逆コンプトン散乱）および陽子と電子間のクーロン相互作用が、源/吸収項として含まれた。
• 降着率： シミュレーションは、エディントン率に規格化された質量降着率の範囲で実行された：$\dot{m} = (1.0 - 10) \times 10^{-6} \dot{M}_{Edd}$。
• 放射伝達： 一般相対性放射伝達（GRRT）計算が BHOSS コードを用いて行われ、$163^\circ$ の傾斜角で 230 GHz の画像を生成した。
• 電子分布： 本研究では、熱的マクスウェル・ユッテナー分布と可変の非熱的 $\kappa$ 分布を組み合わせたハイブリッド電子分布関数を用いた。$\kappa$ パラメータは、粒子インセル（PIC）シミュレーションから導出されたモデル（Ball et al. 2018; Meringolo et al. 2023）を用いて、磁化（$\sigma$）およびプラズマベータ（$\beta$）に基づいてパラメトリックに定義された。

主要な貢献と結果
放射冷却の有無を、異なる降着率および加熱 prescriptions に対して比較することにより、本研究は以下の知見をもたらす：

1. 電子温度と密度構造： 放射冷却は、密度が最も高い赤道面周辺の内部ディスク（$r \lesssim 10 r_g$）において電子温度を急激に低下させる。対照的に、ジェット鞘の温度はわずかにしか低下しない。冷却の組み込みは、冷却効率の増加により、より高い降着率においてより薄いディスク構造の形成をもたらす。
2. 画像形態とフラックス： 放射冷却は、より暗い降着リング（シャドウ）と、放射される総フラックスの減少をもたらす。しかし、より拡張され、相対的に明るいジェット構造を生み出す。冷却モデルでは、非冷却モデルと比較して、近側ジェットが画像全体に寄与する割合が増加する。
3. スペクトルエネルギー分布（SED）： 冷却の組み込みは、特定の降着率に対してターンオーバー周波数（SED のピーク）をより低い周波数へシフトさせ、シンクロトロン放射が自己吸収されにくくなることを示している。230 GHz において、非冷却モデルはより高い降着率で光学的に厚くなるのに対し、冷却モデルは光学的に薄い領域に留まる。
4. 時間変動： 本研究は、230 GHz の光度曲線における時間変動が、冷却および非冷却の両ケースにおいて主に中面領域に由来することを発見した。質量降着率が増加すると、両シナリオにおいて時間変動（変調指数）は減少する。
5. 加熱 prescriptions： 同じ降着率の下では、リコネクション加熱モデルは、冷却の有無にかかわらず、乱流加熱モデルよりも一般的に高い総フラックスを生み出す。しかし、放射冷却の定性的影響（ディスクの暗化、ジェットの明るさ）は、加熱 prescriptions にかかわらず一貫している。

意義と主張
著者らは、非熱的電子分布を伴う二温度 GRMHD シミュレーションに放射冷却を組み込むことが、M87* のブラックホールシャドウおよび大規模ジェット構造を正確にモデル化するために不可欠であると主張する。本研究は、冷却を無視することが、電子温度やフラックス変動に関する不完全なモデルや矛盾した予測につながる可能性があることを示している。

論文は、特定の特徴（拡張ジェットの可視性の向上や特定のフラックス減少など）は現在の EHT 観測のダイナミックレンジを下回っているが、ngEHT などの次世代アレイに対する検証可能な予測を表しているとの結論で締めくくられる。結果は、磁気的に停止したディスク（MAD）状態にある他の低光度活動銀河核（LLAGN）にも一般的に適用可能として提示されており、放射冷却はブラックホールの質量や距離に依存せず、普遍的にディスクを冷却し、ジェットを明るくし、総フラックスを減少させ、光度曲線を滑らかにすると示唆されている。著者らは、シミュレーションがまだ M87* のデータに完璧に適合するように調整されていないこと、および非熱的電子に対する冷却項をより自己整合的に統合して広帯域スペクトルを完全に再現するには将来の作業が必要であるという、控えめな立場を維持している。