Interaction between disk and extended corona in a general relativistic framework

技術要約：一般相対論的枠組みにおけるディスクと拡張コロナの相互作用（パートI）

問題提起
X線連星（XRB）における二相流降着モデルの現行モデルは、光学的に厚いディスクがブラックホール（BH）の近くまで存在し、硬X線を放射するコロナによって照らされているハード状態の観測結果を説明するのに苦慮していることが多い。先行研究（例：Haardt & Maraschi 1991; Poutanen et al. 2018）は、ディスク上に静的なスラブ状コロナが存在する場合、エネルギー平衡においてコンプトン化できるコンプトン化パワーに上限があることを確立した。もしコロナが光学的に厚くなりすぎたり、高温になりすぎたりすると、ディスクへの照射強化がシード光子を増加させ、それがさらにコロナの出力を押し上げるという暴走的なカスケードが発生し、自己整合的な平衡状態を妨げることになる。さらに、既存のモデルは多くの場合、局所的なエネルギーバランスを仮定しており、径方向の変化や一般相対論（GR）効果を無視しているか、あるいはコロナからディスクへのフィードバックを組み込んでいないコードに依存している。本研究は、完全なGRフレームワーク内で、ディスク・コロナ系のグローバルに自己整合的な平衡解を達成するというギャップに対処するものであり、特にハード状態における静的なスラブ幾何学の生存可能性を検証するものである。

手法
著者らは、モンテカルロ放射伝達コード MONK (Zhang et al. 2019) をアップグレードし、降着ディスクと拡張コロナの間の相互作用を組み込んだ。主な手法的な進展は以下の通りである：

1. 反復的なグローバル・エネルギー平衡： コードは、グローバルな平衡に達するまで、コロナによるコンプトン化とディスクへの照射を反復的に計算する。基本的な制約は、全観測輝度（$L_{tot}$）が全利用可能降着パワー（$L_{acc} = \eta \dot{M}c^2$）に等しくなければならないということである。
2. アルベドとフィードバック機構： 降着円盤によって反射される入射コロナフラックスの割合を定義するためにアルベドパラメータを導入する（Chandrasekharの定式化を用いた半無限電子大気としてモデル化）。残りの割合は吸収され、ディスクの実効温度と固有の散逸を増加させる。
3. 粘性散逸パラメータ ($\alpha$)： すべての降着エネルギーがコロナで散逸されると仮定する従来のモデルとは異なり、本モデルではディスク内で固有に散逸される粘性パワーの割合 $\alpha$ を調整する。コードは、$L_{tot} = L_{acc}$ となるように $\alpha$ を反復計算する。もし $L_{tot} > L_{acc}$ であれば、シード光子フラックスを下げるために $\alpha$ を減少させ、逆の場合も同様である。
4. 一般相対論的フレームワーク： フォトン（光子）の伝搬は、カー計量におけるヌル測地線を通じて追跡される。コードは、光の曲がり、重力赤方偏移、および4つの可能なフォトン軌道（無限遠へ到達する、コロナに衝突する、ディスクを自己照射する、またはBHへ落下する）を考慮に入れている。
5. スペクトル解析： コンプトン化成分のフォトン指数（$\Gamma$）は、複数の散乱次数が重なるエネルギー範囲内で、出現スペクトルをべき乗則（$N(E) \propto E^{-\Gamma}$）にフィットさせることで導出される。

主な結果
参照パラメータ（BH質量 $10 M_\odot$、スピン $a=0$、コロナ温度 $T_e = 120$ keV、アルベド $= 0.5$）を用いた結果、以下の知見が得られた：

• 平衡解： ハード状態および中間状態の両方において、グローバルな平衡解が存在する。ディスクを完全に覆う静的なスラブ状コロナの場合、典型的な条件下で達成可能な最低のフォトン指数は約 1.7–1.8 である。
• 光学的厚さの上限： コロナの光学的厚さには厳格な上限（参照ケースでは $\tau_{max} \approx 0.16$）が存在し、これを超えると物理的な平衡解が存在しなくなる（$\alpha < 0$ を必要とする）。$\tau$ が増加するにつれて、エネルギーバランスを維持するために固有のディスク散逸率 $\alpha$ は減少する必要があり、最終的にはディスク放射が完全に再処理されたコロナ照射によって支配される、ほぼゼロの値に達する。
• パラメータ依存性：
  • BHスピン： 高いスピン（$a=0.998$）は、光の曲がりとディスクへの照射が増加するため、スペクトルを硬くする（$\Delta \Gamma \sim -0.15$）。これは、より低い $\alpha$ を強制するためである。
  • コロナ温度： 高い $T_e$ は、より硬いスペクトルをもたらす。
  • アルベド： 高いアルベドは、スペクトルを著しく硬くする。アルベドが1.0（完全反射）の場合、システムはより高い光学的厚さ（$\tau \approx 0.6$）を維持でき、フォトン指数を $\Gamma \sim 1.0$ まで下げることが可能になる。これは、高いアルベドが再処理された軟らかいシード光子を抑制し、フィードバックループをデカップリング（分離）するためである。
• スペクトル形状： コンプトン化されたスペクトルはしば引数の単純なべき乗則から逸脱することが多く、高エネルギー側で「異方性ブレイク」を示す。これは、シード光子の片側からの注入とGR効果に起因する特徴である。
• 径方向のエネルギー流： 空間的に一様な拡張スラブコロナの仮定の下では、各径方向のアニュラス（環）において局所的な平衡を達成するために、グローバルな平衡を満たすことに加えて、内側から外側の領域へのエネルギー移動が必要となる。

意義と結論
アップグレードされたMONKコードを用いることで、ディスク・コロナ系のグローバルなエネルギー平衡を達成することが可能となり、システムパラメータを制約するための堅牢な枠組みを提供できる。しかし、本研究は、静的なスラブ状コロナは、XRBにおける最も硬い観測スペクトルを説明するための好ましい幾何学ではないという結論に至った。たとえグローバルなエネルギー平衡が強制されたとしても、静的なスラブがフォトン指数を1.7–1.8よりも有意に低く生成できないことは、代替の幾何学が必要であることを示唆している。著者らは、今後の課題として、これらの限界に対する潜在的な解決策として、**切断されたディスク幾何学（truncated disk geometries）および流出するコロナ（outflowing coronae）**を探索することを述べている。このアプローチは、フラックスに加えて偏光データを用いてコロナの幾何学を制約することにも応用できる可能性がある。