How to raise a supermassive black hole: interpreting early JWST AGN with the AESOPICA simulations

技術要約：超大質量ブラックホールをいかに育てるか：AESOPICAシミュレーションによる初期JWST AGNの解釈

問題提起
ジェイムズ・ウェッブ宇宙望遠鏡（JWST）は、高赤方偏移（$4 \lesssim z \lesssim 11$）において、局所的なブラックホール–銀河のスケーリング関係と比較して過大質量に見える活動銀河核（AGN）の集団を明らかにしている。これらの知見は、重いシード（$\sim 10^5 M_\odot$）とエディントン限界内の降着を通常想定する標準的なブラックホール形成・成長モデルに疑問を投げかけている。これらの天体が効率的な初期成長チャネルを辿っているのか、あるいは選択効果によって偏ったより広範な集団の観測可能な裾の部分を代表しているのか、あるいは質量推定における系統的な不確実性に起因するものなのかは依然として不明である。これらのシナリオを識別するには、観測的な選択を考慮しつつ、シード質量、降着効率、およびフィードバック機構を系統的に変化させることができる宇宙論的シミュレーションが必要である。

手法
著者らは、FABLE銀河形成モデル（Illustrisに由来）に基づく12個の中規模ボリューム（$L=60$ Mpc）の宇宙論的シミュレーション・スイートであるAESOPICAプロジェクトを紹介する。このスイートは、以下の3つの主要要素を変化させることで、ブラックホール進化のパラメータ空間を系統的に探索している。

1. シード質量： 理論的な形成チャネルを代表するために、4つの異なる初期質量がテストされている：極めて軽い（$10^2 M_\odot$、Pop III残骸）、軽い（$10^3 M_\odot$）、中間的（$10^4 M_\odot$、暴走恒星衝突）、および重い（$10^5 M_\odot$、直接崩壊）。
2. 降着物理： 3つの降着処方箋が用いられている：
   • Fiducial（標準）： ブースト係数 $\alpha=100$ を持つ標準的なボンディ・ホイル・リトロー・降着、およびエディントン限界（$f_{\rm Edd, max}=1$）。
   • SE-Boost（超エディントン・ブースト）： $f_{\rm Edd, max}=10$ までの超エディントン・バーストを許容する強化されたブースト係数（$\alpha=10^3$）。
   • SE-BoostMax（最適化された超エディントン・ブースト）： ボンディ降着の $M_{\rm BH}^2$ 依存性を補償するためにブースト係数がシード質量に反比例し（$10^2 M_\odot$ のシードに対して最大 $\alpha=10^5$）、かつ重いシードに対しては初期成長を促進するために超新星フィードバックが減少している（$\epsilon_{W,SN}=0.5$）「楽観的」なシナリオ。
3. 観測的選択： 著者らは、模擬JWSTサーベイ・パイプラインであるbalmersopicaを開発した。このツールは、様々な露出時間（1時間から100時間）に対してNIRSpecスペクトル（R1000およびR2700回折格子）をシミュレートし、信号対雑音比（SNR）と輝線幅の制約（FWHM $750-10,000$ km s$^{-1}$）に基づいて検出可能性を決定するためのベイズ推定を適用することで、広幅線AGNの検出をフォワードモデリングする。

主な貢献

• 系統的なパラメータ空間の探索： AESOPICAは、宇宙論的ボリューム全体にわたってシード質量、降着効率（超エディントン・バーストを含む）、およびフィードバック強度を同時に変化させる最初の流体計算シミュレーション・スイートの一つである。
• フォワードモデリングによる選択効果： balmersopicaを適用することで、本研究はシミュレーションの生の出力と観測を単に比較するのではなく、特定のJWSTサーベイ構成（CEERS, RUBIES, JADES）を模倣した合成カタログを生成している。
• シーディングと成長の分離： シミュレーションは、効率的な降着が初期の「シード質量の刻印」を急速に消し去り、特定の条件を満たさない限り、高赤方偏移の過大質量ブラックホールの起源がシード質量に対して退化することを示している。

結果

• JWST AGNの再現： 降着効率が十分に高い場合、どのようなシード質量（$10^2$–$10^5 M_\odot$）からでも、JWSTで検出されるAGN集団の大部分を再現できる。ただし、非常に軽いシード（$10^2 M_\odot$）は、観測された質量スケールに達するために最も好ましい条件（SE-BoostMax）を必要とする。
• スケーリング関係：
  • $M_{\rm BH}$–$M_{\rm stellar}$： Fiducial降着では、中間的または重いシードのみが $z \gtrsim 4$ までに過大質量の領域に到達する。SE-BoostMaxを用いると、軽いシードもこの領域に入り、シードの痕跡は $z \sim 6$ までに消失する。$z=0$ において、すべてのモデルは局所的な関係へと収束する。
  • $M_{\rm BH}$–$\sigma_{\rm stellar}$： 選択されたAGNは、局所的で弱く進化する $M_{\rm BH}$–$\sigma_{\rm stellar}$ 関係上に位置しており、これはブラックホールが恒星成分が完全に確立される前に組み立てられる「ブラックホール・ファースト」のシナリオを支持している。
• 選択バイアス： JWST AGNの明らかな過大質量性は、主に観測的選択の結果である。広幅線選択は、自然に最も質量が大きく、最も活発に降着しているブラックホールを回収するが、これらは非効率な降着モデルにおいて $M_{\rm BH}$–$M_{\rm stellar}$ 関係から最も遠く離れた位置にある。
• 光度および質量関数： 推定された高赤方偏移AGNの存在量は、すべての降着モデルにおいて中間的および重いシードと一致している。非常に軽いシードは、高度に効率的な超エディントン成長を経ない限り、否定的である。ブラックホール質量関数は、低質量端（$M_{\rm BH} < 10^7 M_\odot$）においてシード質量への感受性を保持しているが、高質量端では収束する。
• 金属量制約： 気相金属量は、潜在的な識別子となり得る。重いシードモデルのみが、極めて低い金属量（$Z < 0.01 Z_\odot$）を持つホストにおいて $M_{\rm BH}/M_{\rm stellar} > 0.1$ となる系を生み出す。しかし、最も楽観的な重いシードモデルであっても、最も極端な「原始的」なAGN（例：QSO1）を再現するには至らず、これは現在のISMモデリングの限界、あるいは代替的な起源（例：原始ブラックホール）の必要性を示唆している。
• クエンチングと休眠AGN： 高赤方偏移における強くクエンチされた大質量銀河は、超エディントン・バーストを許容するシミュレーションにおいてのみ出現する。休眠状態の過大質量ブラックホール（例：GN-1001830）は、急速な降着バーストの減少期に観測され、その後の強力なフィードバックが星形成と降着の両方を停止させるシステムとして自然に説明される。

意義と主張
本論文は、効率的な降着チャネルが許容されるならば、初期JWST AGNの「過大質量」な性質は必ずしも特定のシーディング機構（例：排他的な重いシード）を必要としないと主張している。むしろ、観測された集団は、明るい広幅線天体を好む選択効果によって大きく偏った、より広い分布の「氷山の一角」である可能性が高い。

著者らは、今後の制約は以下の点に依存すると強調している：

1. 低質量端の統計： ブラックホール質量関数の低質量端を調査する深いサーベイ。ここではシード質量の痕跡が残存している。
2. 金属量測定： $M_{\rm BH}$–金属量関係をテストするための高分解能分光。これは現在、最も極端な天体に対してさえ、最も楽観的な重いシードモデルに課題を突きつけている。
3. 改良された物理学： 現象論的なブースト係数を、物理的根拠に基づいた処方箋に置き換えるための、次世代の多相ISMおよび降着円盤物理を解像するシミュレーションの必要性。
4. 重力波： $10^4$–$10^7 M_\odot$ の範囲での合体を捉える将来のLISA観測は、電磁波的な選択バイアスから自由な、初期の組み立て履歴に対する独立した制約を提供する。

本研究は、最初の超大質量ブラックホールの起源を解明するには、より深い観測的センサス、より洗練されたサブグリッド物理、およびマルチメッセンジャーの制約を組み合わせることが必要であると結論付けている。