Simulating AGN wind feedback with variable feedback efficiencies in idealised disc galaxies

以下は、提示された論文「Simulating AGN wind feedback with variable feedback efficiencies in idealised disc galaxies（理想的な円盤銀河における可変フィードバック効率を用いた AGN ウィンドフィードバックのシミュレーション）」の技術的な要約です。

1. 研究の背景と課題 (Problem)

• AGN フィードバックの重要性: 活動銀河核（AGN）からのフィードバックは、特に巨大銀河の成長を抑制し、星形成を停止させる（クエンチング）上で決定的な役割を果たすと考えられています。
• 既存モデルの限界: 現在の宇宙論的シミュレーション（EAGLE, IllustrisTNG など）では、AGN フィードバックは現象論的な prescriptions（経験則）で実装されており、フィードバック効率（結合効率 $\eta$）はブラックホール（BH）の質量や降着率に依存しない「定数」として扱われることが一般的です。
• 理論と観測のギャップ: 数値モデルや観測（特に超高速アウトフロー UFOs）は、AGN 降着円盤から放出されるウィンドの結合効率 $\eta$ が、エディントン比 $\dot{m}$（BH 降着率をエディントン率で規格化した値）に対してべき乗則（$\eta \propto \dot{m}^{\alpha}$）で依存することを示唆しています。しかし、この物理的に動機付けられた依存関係を宇宙論的・理想的な銀河形成シミュレーションに組み込んだ研究はこれまで行われていませんでした。
• 目的: 理想的な銀河シミュレーションにおいて、エディントン比に依存する可変結合効率モデルを実装し、それが BH の成長と銀河の進化にどのような影響を与えるかを検証すること。

2. 手法 (Methodology)

• シミュレーションコード: SWIFT（SPH 流体計算コード）を使用し、COLIBRE サブグリッド物理モデルを適用。
• 初期条件 (ICs):
  • 静止したダークマターハロー（Hernquist プロファイル）内に、Milky Way 質量（$M_{200} \approx 1.37 \times 10^{12} M_\odot$）の理想的な円盤銀河を構築。
  • 構成要素：恒星円盤、冷たいガス円盤、熱い銀河間ガス（CGM）、および中心の超大質量ブラックホール（SMBH）。
  • CGM は非球対称な重力ポテンシャル下で動的平衡に達するまで 3 Gyr 間緩和させる新しい手法を採用。
• 新しい AGN フィードバックモデル:
  • 熱的フィードバックのみを実装（AGN ウィンドの運動エネルギーがサブグリッドスケールで衝撃波により熱化すると仮定）。
  • 結合効率 $\eta$ の定義: 定数ではなく、エディントン比 $\dot{m}$ のべき乗則として定義：
    $$\eta = \min\{ N_\eta \dot{m}^{\alpha_\eta}, 1 \}$$
    ここで、$N_\eta$ は正規化係数、$\alpha_\eta$ は傾き。
  • パラメータ設定:
    • 基準モデル（Fiducial Variable）: 観測と Qwind モデルに基づき $\alpha_\eta = 2.6$、$\dot{m}=0.1$ で $\eta=0.1$ となるように $N_\eta$ を設定。
    • 比較対象：従来の定数結合効率モデル（$\eta=0.05$）。
    • 変数実験：BH 質量（$10^6$〜$10^9 M_\odot$）、$N_\eta$、$\alpha_\eta$ を変化させて実験。
• シミュレーション期間: 銀河形成物理をオンにして 3.3 Gyr 間進化させ、BH 導入後 3 Gyr 間のデータを分析。

3. 主要な貢献と結果 (Key Contributions & Results)

A. 結合効率パラメータの影響

• 自己規制（Self-regulation）: 結合効率のパラメータ（$N_\eta$ や $\alpha_\eta$）を変化させても、銀河の星形成率（SFR）や BH の最終的な質量など、大域的な性質は僅かな変化しか示しませんでした。
• メカニズム: COLIBRE モデルにおける BH の自己規制が機能しており、結合効率 $\eta$ が変化すると、BH 降着率 $\dot{m}$ が逆方向に変化することで、注入される総エネルギーが調整されるためです。
• 例外: BH の降着率そのものや質量成長履歴、アウトフローの挙動には明確な違いが見られました。結合効率が低い場合、BH はより急速に成長し、アウトフローは弱くなります。

B. BH 質量の影響

• 支配的要因: 銀河と BH の進化を決定する最も重要な要因は、BH の初期質量であることが確認されました。
• フィードバック強度: 高質量 BH（$M_{BH} \ge 10^7 M_\odot$）では、フィードバックが強く働き、星形成が抑制され（クエンチング）、中心ガスが減少します。
• 可変モデルの特性: 可変結合効率モデルでは、定数モデルに比べて、高質量 BH であっても降着率がより高く維持され、BH の成長が速くなります。これは、フィードバックが弱まると $\eta$ も低下し、ガスが BH に再供給されやすくなるためです。

C. アウトフローと銀河構造

• アウトフロー特性: AGN フィードバックは主に極方向（ポラール）に強いアウトフローを駆動しますが、赤道面では流入（インフロー）が支配的です。これは「銀河噴水（Galactic Fountain）」現象として観測されます。
• モデル比較: 可変結合効率モデルは、定数モデルに比べて初期段階でのアウトフロー速度や質量流出率がやや低い傾向にありますが、長期的には中心部のガス密度をより高く維持し、持続的なフィードバック活動を示します。
• 熱的支配: 50 kpc でのアウトフローエネルギーの大部分は熱的エネルギー（約 70-80%）が占めており、運動エネルギーは 20-30% 程度です。

D. CGM の重要性

• CGM を含むシミュレーションと含まないシミュレーションを比較した結果、CGM は BH への燃料供給源として機能し、長期的な AGN 活動と星形成を維持する上で不可欠であることが示されました。CGM がなければ、中心ガスが枯渇し、AGN 活動が早期に停止します。

4. 意義と結論 (Significance & Conclusions)

• 高赤方偏移の超大質量ブラックホールへの示唆:

  • 可変結合効率モデル（特に結合効率が低い領域）は、BH の自己規制を強化しつつ、BH の成長を加速させる可能性があります。
  • これは、JWST 観測で発見されているような、高赤方偏移（初期宇宙）に存在する「過剰に質量の大きい SMBH」の存在を説明するメカニズムとなり得ます（低結合効率によりフィードバックが弱まり、BH が急速に成長するため）。
  • ただし、このモデルが低赤方偏移（現在）の観測（$M_{BH}-M_*$ 関係）とも整合するかどうかは、赤方偏移による結合効率パラメータの変化や、宇宙論的シミュレーションでの検証が必要です。

• モデルの妥当性:

  • 理想的な銀河シミュレーションにおいて、物理的に動機付けられた可変結合効率モデルを実装しても、定数モデルと比較して銀河の大域的な性質（SFR や恒星質量）には大きな違いをもたらさないことが示されました。
  • しかし、BH の成長履歴やアウトフローの詳細な挙動、特に「自己規制」のメカニズムにおいては、より物理的な記述が可能となり、高赤方偏移の BH 成長問題を解決する鍵となる可能性があります。

• 今後の展望:

  • 本研究は理想的な銀河を対象としているため、将来の宇宙論的シミュレーション（COLIBRE モデルを用いたもの）での検証が不可欠です。特に、高赤方偏移でのエディントン比の増加が結合効率にどう影響するか、そしてそれが BH 質量関数にどう影響するかを調べる必要があります。

総じて、この研究は AGN フィードバックの「定数」という単純化された仮定を見直し、観測・理論に基づく「可変」な効率を導入することで、BH の自己規制メカニズムをより精緻に捉え、初期宇宙の巨大ブラックホールの謎に迫る新たな道筋を示した点に意義があります。