A formation scenario of black hole-envelope systems --viscous hydrodynamics simulation in general relativity--

技術的サマリー：ブラックホール・エンベロープ系の形成シナリオ

問題提起
最近のジェイムズ・ウェッブ宇宙望遠鏡（JWST）の観測により、「リトル・レッド・ドット（LRD）」として知られるコンパクトで高光度の赤方偏移天体の集団が特定された。これらの天体は、独特の V 字型スペクトル、紫外線/可視光バンドでの輝き、X 線/電波バンドでの暗さという特徴を示し、高密度ガスに埋め込まれた超大質量ブラックホール（SMBH）の存在を示唆している。解釈は低質量のココナッツから巨大な静水圧エンベロープ（準星）まで多様であるが、超エディントン降着下におけるこれらのブラックホール・エンベロープ系の形成と進化を支配する物理メカニズムは依然として不明である。具体的には、質量 $M \sim 10^5 - 10^7 M_\odot$ の SMBH への準球状降着流における粘性加熱、光子閉じ込め、およびアウトフロー生成の相互作用は、一般相対性理論において厳密な調査を必要とする。

手法
著者らは、SACRA-2D コードを用いて、一般相対性理論における長期的な軸対称粘性流体力学シミュレーションを実行した。シミュレーションは、質量が $10^5 M_\odot$ から $10^7 M_\odot$ の範囲にある非回転ブラックホールへの降着流をモデル化した。

• 初期条件: 設定には、音速点における固定された静止質量注入率 $\dot{M} \approx 1.5 \times 10^{25} \text{ g/s}$（$10^7 M_\odot$ ブラックホールに対するエディントン率に相当）を持つ準球状のボンディ型流を用い、落下ガス温度を $T \lesssim 10^4 \text{ K}$ と仮定した。
• 角運動量: 特定の角運動量プロファイル $j = j_0 \sin^{2n}\theta$（$n=1$）を課し、トーラスの形成を探るために $j_0$ を変えた。
• 物理モデル:
  • 状態方程式（EOS）: 冷却を現象論的に取り込むためにハイブリッド EOS を採用した。光子閉じ込め領域（$\tau > 1$）では、断熱指数を $\Gamma_{\text{th}} = 4/3$（放射優勢）に設定した。光学的に薄い領域（$\tau < 1$）では、ほぼ等温の EOS（$\Gamma_{\text{th}} \approx 1.01$）を使用した。
  • 粘性: せん断粘性係数 $\nu_{\text{vis}} = \alpha_{\text{vis}} c_s^2 \Omega^{-1}$ を実装し、$\alpha_{\text{vis}}$ の値を 0.01 から 0.05 の範囲で変えた。
  • 放射: 光学深度 $\tau$ をアイコナル方程式により計算し、自由落下時間スケールが拡散時間スケールよりも短い閉じ込め半径 $r_{\text{trap}}$ を定義することで、光子の閉じ込めと拡散をモデル化した。
• シミュレーションパラメータ: シミュレーションは $5 \times 10^6 t_g$（ここで $t_g = GM/c^3$）を超える時間スケールで実行され、ブラックホールの質量を $10^5, 10^6, 10^{6.5},$ および $10^7 M_\odot$ とし、様々な角運動量および粘性パラメータを適用した。

主要な結果
シミュレーションは、ブラックホールの質量に依存する明確な進化経路を明らかにした。

1. 低質量ブラックホール（$M \lesssim 10^6 M_\odot$）:

   • 内部領域に光子閉じ込め領域が形成される。
   • 赤道面近傍での粘性散逸と衝撃加熱が、極軸（$z$ 軸）に沿って顕著なアウトフローを駆動する。
   • このアウトフローは落下物質のラム圧力を克服し、明確な流入・流出構造を形成する。$M=10^5 M_\odot$ の場合、アウトフローは計算領域の外部境界に到達する。$M=10^6 M_\odot$ の場合、$r \sim 3\text{--}4 \times 10^4 r_g$ で停滞する。
   • ブラックホールへの質量降着率（$\dot{M}_b$）はエディントン率（$\dot{M}_{\text{Edd}}$）の約 10% に抑制され、一方エンベロープ質量は急速に成長する。

2. 高質量ブラックホール（$M \gtrsim 3 \times 10^6 M_\odot$）:

   • 光子閉じ込め領域は狭く、極方向のアウトフローを起動させるのに十分なエネルギー生成効率が得られない。
   • アウトフローの代わりに、ブラックホールの周りに高温の対流性エンベロープ（またはトーラス）が徐々に発達する。
   • 熱エネルギーが放出されるのではなく対流によって輸送される、準定常状態へと系が進化する。
   • 光度はエディントン限界（$L_{\text{Edd}}$）に近い値に緩和し、光球（$\tau=1$）における平均温度は $\sim 10^4 \text{ K}$ となる。

3. 一般的な知見:

   • 降着効率: ブラックホールの質量や特定の角運動量（$j_0 \gtrsim 4M$ の場合）に関わらず、妥当な粘性値に対してブラックホールへの質量降着率は一貫して $\dot{M}_{\text{Edd}}$ のオーダー 10% である。
   • エンベロープの成長: $M \lesssim 10^6 M_\odot$ の場合、エンベロープ質量の成長率はブラックホールの成長率を大幅に上回る。もし高降着率が $\gtrsim 10^8$ 年間持続すれば、エンベロープ質量は最終的にブラックホール質量を上回り、準星の構成に至る可能性がある。
   • 対流優勢降着流（CDAF）: $M=10^7 M_\odot$ モデルにおいて、流入率と流出率の半径依存性（$\dot{M} \propto r^\beta$、$\beta \sim 0.6$）は、内部領域に CDAF 状態が確立されていることを示している。

意義と主張
本論文は、LRD の性質を説明するブラックホール・エンベロープ系の形成シナリオを提案する。著者らは以下を主張する。

• 質量依存進化: 超エディントン降着流の運命はブラックホールの質量によって決定される。低質量系は降着を調節する極方向のアウトフローを駆動するのに対し、高質量系は顕著なアウトフローなしに対流性エンベロープを発達させる。
• 準星の形成: 低質量ブラックホール（$M \lesssim 10^6 M_\odot$）において、遅いブラックホールの成長（エディントン率の約 10%）と急速なエンベロープ蓄積との格差は、降着が $\sim 10^8$ 年間継続すれば、これらの系が自然に準星（巨大な自己重力を持つエンベロープに囲まれたブラックホール）へと進化することを示唆している。
• 観測との整合性: 模擬された系はエディントン限界に近い光度と $\sim 10^4 \text{ K}$ の光球温度を示し、LRD の性質と整合する。ただし、模擬された温度は一部の LRD で観測される $\sim 5000 \text{ K}$ よりもわずかに高い。これは、LRD が高度に活発な降着段階にないことを示唆する可能性がある。
• 限界: 著者らは、結果が簡略化された状態方程式と現象論的粘性モデルに依存していることに言及している。また、磁場（例えば、ブランフォード・ズナジェック機構）を考慮するためには将来の磁気流体力学（MHD）シミュレーションが必要であり、エンベロープの自己重力がブラックホール質量と同等になった場合、本研究では完全に含まれていなかったことを認めている。