Characterizing the Scale Height and Filamentary Structure of Radiatively Cooled MADs

技術的サマリー：放射冷却を伴う MAD のスケールハイトとフィラメント構造の特性評価

問題提起
放射冷却は、現実的なブラックホール降着円盤の力学における基本的な構成要素であるが、それが磁気的に停止した円盤（MAD）の構造に及ぼす具体的な影響については未だ十分に探求されていない。以前の一般相対論的磁気流体力学（GR-MHD）シミュレーションでは MAD が広範に研究されてきたが、その多くは極めて低光度源に適した理想的な非放射流に焦点を当てていた。放射を取り入れた既存の研究は、計算の複雑さとパラメータ空間の次元性によって制限されることが多かった。さらに、MAD の垂直方向の厚さを特徴づける上で重要な問題が生じる。従来の円盤スケールハイトの定義は、赤道面に対して対称な密度分布を仮定している。しかし、MAD 状態では、強い磁束の蓄積が交換不安定を駆動し、赤道面から大きく逸脱することが多い離散的でフィラメント状の降着ストリームを形成する。その結果、標準的な指標は、特に放射冷却が熱的圧力支持を変化させる場合、円盤の真の幾何学的厚さや力学的状態について誤った結果をもたらす可能性がある。

手法
著者らは、質量 $M = 4 \times 10^6 M_\odot$ の非回転（シュワルツシルト）ブラックホール周りで、GPU 加速コード cuHARM を用いた一連の 3 次元 GR-MHD シミュレーションを実施した。シミュレーションは、質量降着率 ($\dot{M}$) が $10^{-7}$ から $10^{-4} \dot{M}_{\text{Edd}}$ の範囲にある MAD 状態をモデル化した。

主要な手法の構成要素は以下の通りである：

• 放射冷却： シミュレーションでは、シンクロトロン放射および制動放射による放射損失を一貫して含めた。コンプトン冷却は、これらの降着率では支配的ではなく、完全な放射輸送を必要とするため除外された。
• 熱力学： 陽子と電子の温度比 ($T_p/T_e = 3$) を固定すると仮定する単一温度 prescriptions を採用し、種族間の強い結合を仮定した。冷却率 $\Lambda$ は、局所的なガス密度、電子温度、および磁場エネルギー密度に基づいて計算された。
• 初期条件： シミュレーションは、静水圧平衡状態にある非磁化の軸対称トーラスから開始し、磁気回転不安定（MRI）を誘発するために摂動を与えた。MAD 状態の発達が保証されるよう、特定の磁場構成（単一ループのベクトルポテンシャル）が使用された。
• 解析： 時間平均された診断量は、$1.5 \times 10^4 t_g < t < 3.0 \times 10^4 t_g$ の区間内で計算された。著者らは、標準的な密度重み付き角スケールハイトと、フィラメント構造をよりよく捉えるために局所密度の極大の極方位置に基づく新たに提案された指標を比較した。

主要な貢献と結果

1. シンクロトロン冷却の支配性：
   この研究は、シミュレーションされたすべての降着率において、シンクロトロン放射が支配的な放射冷却メカニズムであることを確認した。シンクロトロン冷却は、制動放射による損失を数桁上回っている。シンクロトロン冷却は、高い磁場強度、密度、および温度に駆動され、ブラックホールに近い円盤内縁部および赤道面沿いでピークに達する。

2. 放射効率と円盤構造：
   放射効率 ($\eta_{\text{rad}}$) は、質量降着率とともに非線形的に増加する。低い率では ($\dot{M} \lesssim 10^{-6} \dot{M}_{\text{Edd}}$)、流れは放射非効率降着流（RIAF）として振る舞う。しかし、より高い率（$10^{-4} \dot{M}_{\text{Edd}}$ に近づく）では、冷却が効率的になり、熱的圧力を著しく減少させるほどになる。これにより、円盤が放射効率的となり（最高率では $r=5r_g$ で約 30% の効率に達する）、ガスがより薄く高密度なフィラメントに圧縮されるという遷移が生じる。

3. 従来のスケールハイト定義の失敗：
   著者らは、従来の密度重み付きスケールハイトの定義（式 14）が MAD に対して誤解を招くことを実証した。MAD 流れは、磁気交換不安定により赤道面から変位したフィラメント構造を示すため、赤道面全体で平均化すると、測定された円盤の高さが人工的に過大評価される。従来の指標は、降着ガスの物理的厚さとフィラメントの垂直方向の変位を区別することができない。

4. 新しいスケールハイト指標の提案：
   これに対処するため、本論文は、方位角スライス内の局所密度の極大の位置に基づく円盤高さの代替定義（式 15）を導入した。この手法は、高密度の coherent な「島」を同定し、低密度のコロナ領域を除外し、これらのフィラメントの厚さをその局所中心に対して計算する。

   • 結果： この新しい指標を用いると、円盤高さはブラックホールに向かって単調に減少し、事象の地平線ではゼロに近づくことが示された。これに対し、従来の指標は地平線近傍で厚さの人工的な増加を示していた。
   • 冷却の影響： 新しい指標は、降着率 $\dot{M} \gtrsim 10^{-4} \dot{M}_{\text{Edd}}$ において、放射冷却が内側の $\sim 30-40 r_g$ 内で円盤スケールハイトを著しく減少させることを明らかにした。一方、より低い降着率ではそのような薄化は見られなかった。

5. 臨界降着率：
   著者らは、放射冷却が力学的に有意（加熱よりも速い）となる臨界質量降着率に対する解析的条件を導出した。シンクロトロン冷却を乱流および磁気リコネクションを介した磁気加熱とバランスさせることで、この遷移が $r=30r_g$ で $\dot{M} \sim 2.17 \times 10^{-4} \dot{M}_{\text{Edd}}$ 付近で起こると推定された。この解析的推定値は、$10^{-4} \dot{M}_{\text{Edd}}$ の S-4 シミュレーションにおける構造変化の数値的観測と一致している。

意義と主張
本論文は、放射冷却が熱的圧力支持を減少させることで MAD の内部構造を根本的に変化させ、降着フィラメントを圧縮し、円盤のスケールハイトを減少させると主張している。この研究の主な意義は以下の点にある：

• 診断バイアスの修正： 従来の円盤厚さの測定値が、特に強い冷却が存在する MAD に対して不十分であることを確立し、フィラメントを考慮したより堅牢な定義を提案している。
• 遷移領域の特定： 冷却が minor な摂動から円盤の幾何学と力学を再形成する支配的な力へと変化する臨界降着率の閾値を特定している。
• 観測的含意： 著者らは、これらの冷却によって引き起こされる構造変化（より薄く高密度なフィラメント）が、スペクトルエネルギー分布や光変曲の揺らぎなどの観測可能な特性に刻印されるべきであると示唆しており、Sgr A* や X 線連星などの系における質量降着率を推定する新たな手法を提供する可能性がある。

本研究は範囲において控えめであり、非回転ブラックホールに解析を限定し、イオン - 電子温度比を固定と仮定している。また、ブラックホールのスピンや、より高い降着率におけるより複雑な電子 - 陽子結合モデルの影響を探るには、将来の研究が必要であると指摘している。