Implications of the nanoHertz Gravitational-Wave Background for Galactic Feedback and Massive Black Hole Growth

技術的概要：ナノヘルツ重力波背景が銀河フィードバックと超大質量ブラックホールの成長に及ぼす含意

問題提起
パルサータイミングアレイ（PTA）による最近の測定は、超大質量ブラックホール連星（SMBHB）の非干渉的な重ね合わせによって生成された可能性が高い、確率的なナノヘルツ重力波背景（GWB）の強力な証拠を提供した。この背景の振幅は、超大質量ブラックホール（SMBH）の質量関数に敏感である。しかし、緊張関係が存在する：宇宙論的流体力学シミュレーション（IllustrisTNG、MillenniumTNG、Simba など）は、PTA データから推定される GWB 振幅を一貫して過小評価している。この不一致は、SMBH の成長、種生成、およびそれらを調節するフィードバック機構（AGN フィードバックと恒星フィードバック）に関する現在のモデルが不完全であることを示唆している。具体的には、銀河フィードバックの処方箋の変化がブラックホール質量関数（BHMF）の高質量端をどのように変化させ、その結果として予測される GWB 信号にどのような影響を与えるかは不明である。

方法論
著者らは、Casey-Clyde ら（2022）に基づく半解析的枠組みを用いて、シミュレーションで予測された BHMF と GWB 振幅を結びつけた。方法論は以下の通りである：

1. シミュレーションスイート： 本研究では、IllustrisTNG、MillenniumTNG、Simba、および CAMELS（Cosmology and Astrophysics with MachinE Learning Simulations）プロジェクトの 4 つの主要なシミュレーションスイートを分析した。
2. ボックスサイズと収束性： GWB が BHMF の高質量端に敏感であるという点に対処するため、著者らは複数のボックスサイズ（25 Mpc/h から 500 Mpc/h の範囲）を用いて収束性研究を実施した。その結果、TNG の場合約 35 Mpc/h 未満、Simba の場合約 50 Mpc/h 未満のボックスは、ポアソンノイズの影響を受け、高質量ハロー（$M_{BH} > 10^9 M_\odot$）のサンプリングが不十分であることが判明した。
3. フィードバック変異： 本研究は、銀河フィードバックの影響を以下のように探求した：
   • 異なるシミュレーションコード間での基準モデル（fiducial models）の比較。
   • CAMELSスイートにおける「1P」（単一パラメータ）変異の分析。ここでは、恒星フィードバックとAGNフィードバックを制御する23の天体物理パラメータを個別に変化させた。
   • AGN ジェットモードを除去または変更した特定の Simba フィードバック変異（例：「SW のみ」、「SW + Rad」）の検討。
4. GWB 計算： シミュレーションの BHMF を、活動性クエーサーを同定するためのエディントン比分布関数（ERDF）を用いて連星集団に変換した。この集団は、連星とクエーサーの比率を考慮するために正規化係数（$N_{cal} \approx 4$）でスケーリングされ、二次質量は質量比分布を通じて割り当てられた。得られた連星集団を積分し、基準周波数 $f = 1 \text{ yr}^{-1}$ における特性ひずみスペクトルと GWB 振幅を計算した。

主要な貢献

• フィードバック影響の定量化： 本論文は、恒星フィードバックと AGN フィードバックのモデルの選択が BHMF の高質量端を著しく変化させることを実証した。非効率的または欠如した AGN フィードバックは、過剰な大量の SMBH を生み出し、予測される GWB 振幅を、基準変異では 2 倍、極端な CAMELS 変異では 10 倍増加させる。
• 収束性分析： 著者らは、頑健な GWB 予測に必要な最小ボックスサイズを確立し、標準的な 25 Mpc/h の CAMELS ボックスでは $M_{BH} > 10^9 M_\odot$ の統計が収束していないことを特定した。
• 体系的過小評価： 本研究は、IllustrisTNG、MillenniumTNG、および Simba の基準モデルが、PTA から推定される GWB 振幅を約 2 倍過小評価していることを確認した。観測と一致させるためには、これらのモデルは $M_{BH} > 10^8 M_\odot$ の範囲で 1.5 倍から 5.5 倍多くの SMBH を必要とする。

結果

• 基準モデル： 全ての標準的（基準）シミュレーションモデルは、NANOGrav 15 年データと比較して GWB 振幅を過小評価している。これは、現在のフィードバック処方箋が SMBH の成長を過度に効率的に抑制しているか、あるいは種生成が遅すぎることを意味する。
• フィードバック変異： AGN フィードバックが減少または欠如したシミュレーション（特にジェットモードなしの Simba 実行）は、PTA データと整合する GWB 振幅を生み出す。しかし、これらのモデルは現実的な銀河集団を再現できず（例えば、星形成の抑制が効果的でない）、物理的に非現実的である。
• パラメータ感度： CAMELS 1P セットにおいて、特定のパラメータが GWB 振幅を駆動する：
  • CAMELS-TNG の場合、低・高降着 AGN 状態間の切り替えを制御するパラメータ $\beta$ が最大の影響を持つ。$\beta=0$ の値は PTA データと一致するが、低質量の恒星質量関数を過小評価する。
  • CAMELS-Simba の場合、種質量比パラメータ（$A_{BHseed}$）が最も重要である。種閾値を低下させることで SMBH の成長を早期化させ、PTA データと一致させることができるが、初期の星形成を乱す可能性がある。
• 恒星フィードバック： しばしば低質量ブラックホールにのみ影響すると仮定されるが、本研究は、特定の恒星フィードバックパラメータ（例えば、風速度の正規化と熱的割合）が、ブラックホール成長に利用可能な燃料を調節することで、高質量 BHMF と GWB 振幅にも影響を与えることを発見した。

意義と主張
本論文は、PTA 測定が、従来の電磁気的制約を補完する、銀河フィードバック物理学の強力かつ独立したプローブを提供すると主張している。シミュレーション予測と PTA データの不一致は、初期宇宙における SMBH の成長が、現在のモデルで捉えられているよりも効率的であるか、あるいはより早期に起こっていることを示唆している。

著者らは、フィードバックパラメータや種生成処方箋を修正することでシミュレーションを PTA データと人工的に整合させることは可能だが、そのような修正はしばしば観測された銀河統計（例えば恒星質量関数）との整合性を破ると結論づけている。したがって、GWB 振幅とシミュレーション予測の間の緊張関係は、宇宙論的流体力学シミュレーションにおける SMBH の種生成と成長の処方箋を根本的に再評価する必要性を示している。本研究は、将来のモデルが現実的な銀河集団と観測された GWB 振幅の両方を同時に再現しなければならないことを強調しており、現在のフィードバックモデルが SMBH の成長を過度に抑制している可能性を示唆している。著者らは、自らの枠組みが環境硬化や合体タイミングの詳細を無視するなど、単純化された仮定に依存しており、将来の精度向上には連星進化の自己無撞着な処理が必要であると指摘している。