Supermassive Black Hole Growth in Massive Galaxies at Cosmic Dawn

技術要約：宇宙の黎明期における大質量銀河内の超大質量ブラックホールの成長

問題提起
ジェイムズ・ウェッブ宇宙望遠鏡（JWST）による最近の観測は、宇宙の黎明期（$z \gtrsim 9$）において、高い星形成率（SFR）と高い恒星質量（$M_\star \gtrsim 10^8 M_\odot$）を持つ、大質量で明るい銀河の集団を明らかにしている。これらの銀河のいくつかは、活動銀河核（AGN）を示すスペクトル特徴を示しており、局所的なスケーリング関係に対して過大な質量を持つ超大質量ブラックホール（SMBH）の存在を示唆している。標準的な銀河形成モデルでは、これらSMBHの急速な成長や観測された高い星形成効率を説明することに苦慮している。理論的な提案では、重いシード（種）や超エディントン降着が呼び出されているが、初期の大質量銀河における、乱流的で多相的な星間物質（ISM）内でのSMBH成長、恒星フィードバック、およびAGNフィードバックの相互作用については、依然として理解が進んでいない。

手法
著者らは、適応格子細分化法（AMR）コードであるramsesを用いた、11組の宇宙論的流体ダイナミクス「ズームイン」シミュレーションを実施した。シミュレーションは、$z \sim 9$までに $\sim 10^{11} M_\odot$ のハローへと崩壊する特定の過密領域に焦点を当てており、そこには大質量銀河（Massive Dawn Galaxy: MDG、 $M_\star \gtrsim 10^8 M_\odot$）が存在する。

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

• 星形成モデル： 標準的なシミュレーションで使用されるケンカット・シュミット則とは異なり、本研究では物理的に動機付けられた、乱流ベースのサブグリッド・マルチフリーフォール・モデルを採用している。これにより、ISMの乱流運動エネルギーによって駆動される可変的な局所星形成効率（SFE）が可能となり、$\gtrsim 10\%$ に達することもある。
• SMBHの実装： SMBHは、$z \sim 12$において局所的なガスおよび恒星質量が $10^8 M_\odot$ のハロー内で $10^5 M_\odot$ を超えた際に配置される、シード粒子としてモデル化されている。シード質量は $10^4 M_\odot$ と $10^5 M_\odot$ の両方がテストされた。
• 降着物理： 降着は、未分解の密度ゆらぎを考慮するためのブースト係数を含む、修正されたボンディ・ホイル・リクトン公式に従う。降着率はエディントン限界によって制限されており、これはシミュレーション全体を通じて $\lambda_{\text{edd}}$（値：1, 3, および制限なしの $\infty$）を介して変化させる。
• フィードバック機構：
  • 恒星フィードバック： 熱的な光電離（初期型）およびタイプII超新星（後期型）を含む。「低恒星フィードバック」のバリアントでは、ISM構造への影響を孤立させるために光電離を除去している。
  • AGNフィードバック： シンク球への等方的な熱エネルギー注入としてモデル化され、結合効率 $\epsilon_c$（0%および15%でテスト）を持つ。
• シミュレーション・スイート： 11のランは、AGNフィードバックの強さ、シード質量、最大降着限界、および恒星フィードバック効率を系統的に変化させ、SMBHと銀河の共進化に対するそれらの影響を分離した。

主要な結果

1. ISM駆動の確率論的降着（饗宴と飢餓）：
   典型的なMDG環境において、恒星フィードバックは、冷たく高密度な（$n_H \sim 10^3 \text{ cm}^{-3}$）領域と、熱く希薄な（$n_H \sim 10 \text{ cm}^{-3}$）領域からなる、乱流的で多相的なISMを生成する。SMBHシンク球は、これらの相を確率的にサンプリングする。

   • 饗宴モード（Feast Mode）： シンクが冷たく高密度なパッチに入ると、音速が低下し、降着が（超）エディントン限界に達することが可能になる。
   • 飢餓モード（Starvation Mode）： シンクが熱く希薄なパッチに入ると、音速が上昇し、ボンディ降着率が著しく低下する。
   • 典型的なモデルでは、星形成の開始後（$t \sim 300$ Myr）、SMBHは約50%の時間、飢餓モード（$\alpha \approx 0.5$）にあり、滑らかなISMと比較して平均成長率を実質的に半分に抑制している。

2. 自己調節の条件：
   AGNフィードバックがSMBHの成長を自己調節できるのは、ブラックホールがその加熱がシンク球内の冷却を支配するほど十分に重くなった場合（$M_{\text{sink}} > M_{\text{sink,crit}}^{\text{cool}}$）のみである。

   • この臨界質量は、ガスの密度と金属量に依存する。乱流的なMDG環境において、SMBHがハローからガスを解離させるほどの暴走的な加熱を実現するためには、熱く希薄なポケットに滞在する時間の $\gtrsim 50\%$ において、この閾値を超える必要がある。
   • エディントン制限付き降着（$\lambda_{\text{edd}}=1$）および標準的なシード質量を用いたシミュレーションでは、金属量の上昇に伴って臨界冷却質量閾値がSMBHの成長よりも速く上昇するか、あるいはSMBHが閾値を超えて滞在する時間が不十分であるため、自己調節に失敗することが多い。
   • 自己調節は、より高いシード質量（$10^5 M_\odot$）または超エディントン降着（$\lambda_{\text{edd}} \ge 3$）を持つシミュレーションにおいて、熱結合効率が十分に低い（$\epsilon_c \lesssim 0.15\%$）場合に達成される。

3. 銀河スケールのクエンチングの欠如：
   低赤方偏移のモデルからの予想に反して、これらのシミュレーションにおけるAGNフィードバックは、グローバルな星形成率（SFR）を抑制したり、ホスト銀河の恒星質量の成長を停止させたりすることはない。

   • ハローは、SMBHの最大成長率（$\dot{M}_{\text{SMBH}}$）を大幅に上回るレート（$\dot{M}_{\text{halo}}$）で、宇宙のフィラメントから冷たいガスを降着させている。
   • ハローの成長タイムスケール（$\tau_{\text{halo}}$）がSMBHの成長タイムスケールと同等かそれよりも短いため、新鮮な冷たいガスが継続的にISMを補充し、AGNフィードバックによる加熱を圧倒する。したがって、銀河はSMBHの活動に関わらず、指数関数的に成長を続ける。

4. 観測との比較：
   宇宙の黎明期における過大なSMBHの集団（例：UHZ-1, GN-z11）に最もよく一致するシミュレーション・スイートは、super_edd_no_AGN ラン（$M_{\text{seed}} \sim 10^4 M_\odot$, $\lambda_{\text{edd}} = 3$, AGNフィードバックなし）である。

   • 観測された質量に自己調節型のSMBHで一致させるためには、熱結合効率が極めて低い（$\epsilon_c \lesssim 0.15\%$）必要がある。
   • 結果は、SMBHが（$10^6 M_\odot$ ほど高くはない）シード質量を必要とせずに自己調節レジームに達するためには、超エディントン降着が必要である可能性が高いことを示唆している。

意義と結論
本論文は、宇宙の黎明期における大質量銀河内のSMBH成長は、AGNフィードバック単独ではなく、主に恒星フィードバックによって駆動される乱流的で多相的なISMの確率的な性質によって制御されていると主張している。著者らは以下のように結論付けている：

• これらの環境におけるSMBHは、頻繁に熱く希薄なガスのポケットによって「飢餓」状態に置かれ、それが超エディントン降着を維持できない限り、その成長を制限する。
• AGNフィードバックは、ハローの急速なガス降着が燃料供給を継続的に補充するため、この時代における銀河スケールの星形成を停止させるには効果的ではない。
• JWSTによって観測されたような大質量SMBHを実現するには、銀河がより安定した低赤方偏移のような状態へと進化する前に、上昇する臨界冷却質量閾値を克服するために、超エディントン降着と非効率なAGNフィードバック結合、あるいは非常に重いシードの組み合わせが必要である。

著者らは、自らの結果が解像度に依存していることを指摘しており、これらの知見を検証し、フィードバックの堆積の物理的スケールを精緻化するためには、非平衡熱化学、放射伝達、およびより高い空間解像度（$\sim 1$ pc）を取り入れた今後の研究が必要であるとしている。