Gas Accretion versus BH Merger driven Growth Modes of Supermassive Black Holes and Implications for the Little Red Dots

技術的サマリー：超大質量ブラックホールのガス降着対ブラックホール合体駆動成長モードおよび「リトル・レッド・ドット」への示唆

問題提起
超大質量ブラックホール（SMBH）の形成および成長メカニズムは、特に宇宙の最初の 10 億年以内（$z \sim 6$）にどのようにして $10^9 M_\odot \ge$ の質量に達したかという点において、未解決の問いとして残されています。恒星質量の種に対する急速な超エディントン降着から、直接崩壊による巨大な種（$10^5 M_\odot$）の形成、あるいは中間質量ブラックホール（IMBH）の合体に至るまで、様々なシナリオが存在しますが、具体的な経路は不明のままです。最近のジェームズ・ウェッブ宇宙望遠鏡（JWST）による観測は、宿主銀河の恒星質量に対して「過剰質量」であるか、あるいは局所的なスケーリング関係と一致する「正常な」質量を持つ SMBH を含む高赤方偏移（$z \sim 4-12$）の「リトル・レッド・ドット」の集団を明らかにし、状況をさらに複雑にしました。ガス降着とブラックホール（BH）合体駆動の成長モードを区別することは、SMBH とその宿主銀河の共進化を理解し、これらの多様な高赤方偏移集団を説明するために不可欠です。

手法
本研究では、Tree-PM および平滑化粒子流体（SPH）法を採用した改変版 GADGET-3 コードを用いて、$(50 \text{ Mpc})^3$ の共動体積の宇宙論的流体力学シミュレーションを実施しました。シミュレーションは、平坦な $\Lambda$CDM 宇宙論の下で、赤方偏移 $z=100$ から $z=0$ まで進化させました。

• 解像度: 標準的な解像度は、$256^3$ のダークマター粒子と $256^3$ のガス粒子を使用し、ガス粒子の質量は $4.61 \times 10^7 M_\odot$、重力軟化長は $3.3 \text{ kpc}$（共動）です。
• 物理: モデルには金属線冷却、時間依存の光電離背景、多相サブ解像度星形成（SF）モデル、および 11 種の原子種を追跡する化学的富化が含まれています。超新星（SN）フィードバックは、質量負荷係数 $\eta=2$ および風速 $v_{SN}=350 \text{ km/s}$ で運動学的に実装されています。
• ブラックホールの種付けと成長: BH は、$M_{halo} > 10^{10} M_\odot$、$M_* > 5 \times 10^7 M_\odot$ の巨大ハローに、初期質量 $10^4$、$10^5$、または $10^6 M_\odot$ で種付けされます。成長は、ボンド・ホーレイ・リットルトン降着（エディントン率に制限される）および、粒子が軟化長以内かつ相対速度が音速未満である場合の BH 合体によって行われます。
• フィードバック: 運動学的 AGN フィードバックは、二円錐体積内のガスに運動量を付与することでモデル化されます。フィードバック効率（$\epsilon_f$）および風キック速度（$v_w$）は、$M_{BH} - M_*$ 関係への影響をテストするために、シミュレーション実行間で変化させます。
• シミュレーション実行: 種付け質量（$10^4-10^6 M_\odot$）、フィードバック効率（$0, 0.1, 1$）、およびキック速度（$1000, 5000, 10000 \text{ km/s}$）を変化させる一連の実行（表 1）に加え、BH がない、あるいは BH フィードバックがない対照実行も行われました。

主要な結果

1. 二相成長メカニズム: シミュレーションは、SMBH にとって明確な二段階の成長履歴を明らかにしました。
   • 合体支配相: $10^5 M_\odot$ の種から $\sim 10^7 M_\odot$ への初期成長は、主に BH 合体によって起こります。この相の間、ガス降着は低く、極めてエディントン以下（$\dot{M}_{BH}/\dot{M}_{Edd} < 0.001$）です。
   • 降着支配相: BH が $\sim 10^7 M_\odot$ に達すると、効率的なガス降着によって駆動される急速な成長バーストを経験します。この相は約 $600-700 \text{ Myr}$ 続き、その間エディントン比は 1 に達します。この窓において、BH 質量は $10^2 - 10^3$ 倍増加し、$10^7 M_\odot$ から $(10^9 - 10^{10}) M_\odot$ へと成長します。
2. 成長の停止: この急速な降着相に続いて、SMBH の成長は停止します。シミュレーションは、$z=0$ までに SMBH が $\sim 10^{10} M_\odot$ を超えないことを示しています。この停止は、AGN フィードバックによって銀河中心のガスが枯渇し、さらなる成長を自己規制することに起因しています。
3. $M_{BH} - M_*$ 共進化: シミュレーションされた銀河は当初、局所的な $M_{BH} - M_*$ 相関より下に位置しています。宇宙時間を通じて、それらは関係に向かって上方へ進化します。相関の傾きは比較的一定のままですが、切片は進化します。標準的な実行（$10^5 M_\odot$ の種、$\epsilon_f=0.1$、$v_w=5000 \text{ km/s}$）は、局所的な $z=0$ の観測関係（Marconi & Hunt 2003）を最もよく再現します。
4. パラメータ感度: BH 種付け質量は、最終的な $M_{BH} - M_*$ 関係に最も大きな影響を与えます。より大きな種（$10^6 M_\odot$）は過剰質量の BH につながり、より小さな種（$10^4 M_\odot$）は過小質量の BH をもたらします。フィードバック効率とキック速度の影響は比較的小さいです。
5. 星形成: AGN フィードバックは、$z \le 4$ において星形成率（SFR）を著しく抑制し、BH フィードバックのないシミュレーションと比較して、全球的な星形成率密度（SFRD）を $1.3 - 15$ 倍減少させます。直接降着によるガスの除去は抑制において minor な役割を果たすのみであり、主要なメカニズムはフィードバック誘起の加熱と放出です。

「リトル・レッド・ドット」への示唆
著者らは、これらの成長モードを JWST による高赤方偏移の「リトル・レッド・ドット」の観測解釈に応用しました。彼らは、形成チャネルの二項性を提案します。

• 正常質量 SMBH: これらの天体は、シミュレーションで観測された合体支配相と一致する、主に BH 合体駆動の進化を経た可能性が高いです。
• 過剰質量 SMBH: これらの天体は、シミュレーションで特定された急速成長相と同様の、エディントン制限または超エディントンの効率的ガス降着のバーストを経験した可能性が高いです。初期宇宙の銀河における高い核ガス密度は、この効率的な降着を促進し、宿主恒星質量に対して過剰質量に見える BH へと導きました。

重要性
本論文は、合体駆動と降着駆動の成長領域を区別する物理的枠組みを提供し、JWST によって観測された高赤方偏移 SMBH の多様な集団（正常対過剰質量）が、異なる支配的な成長経路の結果である可能性を示唆しています。本研究は、AGN フィードバックが中央ガス貯留層を枯渇させることで SMBH の成長を $\sim 10^{10} M_\odot$ 以下に自己規制し、宇宙時間を通じて観測された $M_{BH} - M_*$ スケーリング関係および銀河の星形成履歴を成功裡に再現することを示しています。