Black Hole Feedback, Galaxy Quenching and Outflows at Cosmic Dawn: Analysis of the SEEDZ Simulations

技術要約：宇宙の夜明けにおけるブラックホール・フィードバック、銀河のクエンチング、およびアウトフロー

問題提起
ジェイムズ・ウェッブ宇宙望遠鏡（JWST）による最近の観測は、宇宙の最初の10億年以内（$z \gtrsim 6$、さらに $z \sim 12.3$ まで）に存在する、$10^6 M_\odot$ を超える質量を持つ大質量ブラックホール（MBH）の相当な集団を明らかにしている。これらの知見は、コンパクト天体の形成と成長に関する現在の理論モデルに疑問を投げかけている。具体的には、ブラックホールがいかにして初期のシード質量から、割り当てられた時間内にこれほど大きなスケールへと成長できたのかが不明である。さらに、高赤方偏移の観測は、局所宇宙の相関関係（$M_{BH}/M_{star} \simeq 10^{-4}$）と比較して、過大なブラックホール対恒星質量比（$M_{BH}/M_{star} \gtrsim 10^{-2}$）を示しており、これはブラックホールのシード形成の後に恒星成分が構築されるという、異なる進化経路を示唆している。SEEDZシミュレーション・スイートは、宇宙の夜明けにおける銀河およびMBHの進化経路を調査するために設計されており、特に、MBHの成長を制御または停止させる物理メカニズムに焦点を当てている。

手法
本研究は、移動メッシュコードArepo2を用いて実行された、3D流体計算による宇宙論的シミュレーションであるSEEDZシミュレーション・スイートを分析する。シミュレーションは「ズームイン」アプローチを採用しており、親ボックス（$40 \, \text{cMpc} \, h^{-1}$）内の3つの異なる領域（Normal1、Normal2、およびRarepeak）に焦点を当てている。高解像度領域のダークマター粒子質量は $\approx 6.9 \times 10^4 M_\odot h^{-1}$、最小セルサイズは $\approx 6 \, \text{pc}$（$z=15$ における物理量）である。

主な数値的および物理的構成要素は以下の通りである：

• シード形成（Seeding）： シミュレーションでは、「ライトシード」（ポピュレーションIII星の残骸）と「ヘビーシード」（高密度クラスター内での暴走衝突、超大質量星、または準星によって形成される）の両方をモデル化している。ヘビーシードは、べき乗則分布に基づき、$5 \times 10^3 M_\odot$ から $10^5 M_\odot$ の質量で形成される。
• 降着（Accretion）： ブラックホール降着は、カーネル半径を介してシミュレーションの解像度に合わせて調整されたボンディ・ホイル・リトルトン形式を用いてモデル化されており、高角運動量ガスによる過剰降着を防ぐためにガスの渦度によって修正されている。
• フィードバック（Feedback）： 熱エネルギーは降着半径内のガスセル内に等方的に注入される。放射効率（$\epsilon$）は動的に計算される。これは標準的なサブエディントン値（$\sim 0.1$）に従うが、超エディントン降着時には「スリム円盤」モデルへと移行し、光子トラッピングによって実効的な効率が低下し、フィードバックが弱まる。
• 力学（Dynamics）： ソフトニング長による重力相互作用の抑制を考慮するための力学的摩擦補正、および近接性と相対速度に基づく合体基準が含まれる。
• 範囲（Scope）： 解析は、スイート内の最も質量の大きいブラックホールの成長履歴を $z=12.5$ までに焦点を当てて行われる。

主要な結果

1. 成長と質量限界： シミュレーションにおける最も質量の大きいブラックホールは、$z=12.5$ までに $\sim 10^6 M_\odot$ の質量まで成長する。この成長は、5〜30 Myr間続く短時間の超エディントン降着の間に発生する。最終的な質量は初期のシード質量によって決まるのではなく、環境条件によって決定される。
2. フィードバックによるクエンチング： MBHの成長を停止させる主要なメカニズムは、超新星爆発やガスの枯渇ではなく、ブラックホール自身からのフィードバックである。最も活発に降着しているMBHにおいて、注入される熱エネルギーはホストハローの結合エネルギーを超える。これにより、ホストハローから銀河間物質へのガスの完全な排出が引き起こされ、事実上の銀河のクエンチング（活動停止）が生じる。
3. アウトフローの特性： 降着の遮断の瞬間、ガスのアウトフローは $\sim 10^3 \, \text{km s}^{-1}$ の速度と $10^7 \, \text{K}$ の温度に達する。これらのアウトフローは、冷却される前に $\sim 20 \, \text{kpc}$ のスケールまで拡大する。一部のハローは20〜100 Myrのタイムスケールでガスの流入率が回復する兆候を示すが、本シミュレーション内のブラックホールは、シミュレーション期間内に二度目の主要な成長フェーズを経験しない。
4. ハロー質量依存性： 最大ブラックホール質量は、フィードバックモデルとホストハローの結合ポテンシャルの組み合わせによって事実上決定される。シミュレートされたハロー（質量最大 $\sim 10^{10} M_\odot$）において、ガスを解離させるために必要なフィードバックエネルギーは、ブラックホールの成長を約 $10^6 M_\odot$ に制限する。
5. 星形成史： 最も質量の大きいブラックホールを宿すハローは、ヘビーシードの形成前に星形成がほとんど行われないため、効率的な降着が可能となる。一度降着が始まると、その結果としてのAGNフィードバックがさらなる星形成を抑制する。対照的に、恒星成分が支配的なハローは継続的な星形成を示すが、質量の大きいBHホストで見られるような急速かつ妨げのない降着は欠いている。

意義および結論
本論文は、実装されたフィードバックモデルの下では、初期の $\Lambda$CDM宇宙において、大質量ブラックホールが降着のみによって $z=12.5$ までに $\sim 10^6 M_\odot$ を超えることは困難であると結論付けている。急速な降着に伴うフィードバックは、ホストハローのガス貯蔵庫を解離させるのに十分なエネルギーを持ち、それによってブラックホールの成長を自己調節する。

著者らは、高赤方偏移で観測される質量超過したブラックホールの集団に対し、2つの解釈を提案している：

1. より弱いフィードバック： 宇宙初期のMBHフィードバックは、モデル化されているものよりも著しく効率が低い可能性がある（例えば、成長が降着主導ではなく合体主導である場合、あるいはフィードバックが狭い領域内に閉じ込められている場合など）。
2. 二段階形成： 高赤方偏移の銀河は、MBHがまずガス貯蔵庫を排出することでホスト銀河をクエンチし、その後、ガスに富む銀河との合体や宇宙網（コスミックウェブ）からの降着を通じてバリオン貯蔵庫を再構成し、その後の成長を可能にするという、二段階のプロセスを経て形成された可能性がある。

著者らは、現在の解像度が、後の成長段階において支配的なメカニズムとなり得るブラックホール合体を完全に解像する能力に制限されていることを指摘している。将来の高解像度シミュレーションは、これに対処し、特にJWSTによって観測されている「リトル・レッド・ドット（Little Red Dot）」系に関連して、フィードバックパラメータをより厳密に制約することを目的としている。結果は、フィードバックが極めて非効率であるか、あるいは限定的でない限り、赤方偏移12.5以前の大質量ブラックホールの最大質量が $10^6 M_\odot$ を大きく超えることはないことを示唆している。