AGN-driven BBH mergers: Black hole populations and hierarchical growth across the AGN parameter space

技術要約：AGN駆動型BBH合体：AGNパラメータ空間におけるブラックホール集団と階層的成長

問題と動機
活動銀河核（AGN）の円盤は、連星ブラックホール（BBH）の形成と合体における有望なダイナミカル・チャネルとして浮上している。これらの円盤内に埋め込まれた恒星質量ブラックホール（BH）は、ガストルクや重力波放射を介して移動、対形成、および合体を行う可能性があるが、予測されるこの集団の特性は依然として極めて不透明である。先行研究では、円盤構造、移動トルク、ガスの硬化効率、およびAGNの寿命といった物理的要素への強い依存性が指摘されているが、それぞれ異なる仮定を採用しているため、直接的な比較が困難となっている。近年の観測による制約は、AGNチャネルが重力波観測所によって検出されたBBH合体集団のすべてを説明することはできないことを示唆しているが、その具体的な寄与と特徴的なシグネチャーを特定するには、AGNパラメータ空間の系統的な探索を通じて、ロバストな観測的特徴と不確実性の源を特定する必要がある。

手法
著者らは、AGN円盤内におけるBBHの形成、進化、および階層的成長をモデル化するために、fastclusterソフトウェア環境内で実装された最新の半解析的フレームワークを提示する。このアプローチは、捕獲、半径方向の移動、連星形成、ガス駆動の硬化、ダイナミカルな摂動、そして最終的な合体に至るまでの、恒星起源（第1世代、1g）ブラックホールのライフサイクルを追跡する。

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

• 円盤モデリング: 吸収円盤は、定常状態の解析的モデルであるSirko & Goodman (2003) (SG) モデルを用いてモデル化される。著者らは、様々な超大質量ブラックホール質量（$M_\bullet$）、エディントン比（$\lambda_\bullet$）、および円盤粘性パラメータ（$\alpha$）に対して、1次元円盤方程式を自己整合的に解いている。粘性値については、$\alpha = 0.01$（低粘性）と$\alpha = 0.1$（高粘性）の2つの値がテストされている。
• 初期条件: シミュレーションは、核星団（NSC）のスケーリング関係に基づいたBH集団を初期化する。1g BHは、太陽金属量を仮定した集団合成カタログ（sevn）から抽出される。スピンはマクスウェル分布に基づいて割り当てられ、その向きは円盤の角運動量に整列している。
• 進化の物理学:
  • 移動: BHは、Grishin et al. (2024) および Kanagawa et al. (2018) の処方箋に従い、Type I および Type II レジームを通じて移動する。
  • ペアリング: 連星形成は、順行軌道の割合（$f_{progr}$）を考慮した、ガスによる補助基準に従う。
  • 硬化と遭遇: 形成後、連星はガス駆動の硬化（サーカムバイナリ円盤との相互作用）と単一BHとの三体遭遇の組み合わせの下で進化する。著者らは、ガス駆動の硬化（Ishibashi & Gröbner 2024を使用）とGW放射（Peters 1964）を組み合わせた軌道進化方程式を統合している。遭遇は、事前計算されたポスト・ニュートン散乱実験のグリッド（tsunami）を用いてモデル化される。
  • 階層的成長: 合体残骸は、整列スピンの処方箋と反跳キック速度を用いて計算される。円盤に拘束されたまま残る残骸は、AGNの寿命制限および放出基準に従って、後続の合体世代（$N$世代目）に参加する。
• 集団合成: 局所宇宙（$z \lesssim 0.3$）を代表する固有のBBH集団を構築するため、著者らは観測されたAGNの$M_\bullet$および$\lambda_\bullet$の分布（Greene & Ho 2007）に従ってシミュレーション出力を重み付けしている。

主な結果
本研究は、$M_\bullet$、$\lambda_\bullet$、および$\alpha$に対するBBH合体集団の依存性を系統的に探索し、以下の知見を得た：

1. 階層的成長と質量分布:

   • AGN円盤は繰り返しの合体を生成でき、対不安定型超新星（PISN）質量ギャップ（$50 \lesssim m_1/M_\odot \lesssim 130$）を超えて中間質量ブラックホール（IMBH）領域へと延びる高質量テールを形成できる。
   • 粘性依存性: 階層的成長は、ガスの表面密度が高く、ダイナミカルなタイムスケールが短縮される低粘性円盤（$\alpha = 0.01$）において、より効率的に進行する。これにより、より高い合体世代（$N_{gen} > 100$）が可能となる。高粘性円盤（$\alpha = 0.1$）では、高質量の残骸形成が抑制され、通常はより低い世代（$N_{gen} \approx 20$）で飽和する。
   • 質量比: 集団は、高質量においてますます不等な質量比へと進化する。一次質量は連続的な合体を通じて成長するが、二次質量はより低い値（通常はPISNギャップ以下）で飽和する。これは、高世代の一次質量が低世代の二次質量と合体するためである。

2. 合体効率とパラメータ空間:

   • 合体効率（$\eta$）は、一般にエディントン比（$\lambda_\bullet$）とともに増加する。
   • $M_\bullet$への依存性は非自明である。効率と最大残骸質量は中間から大質量領域のSMBH質量でピークに達する傾向があるが、単一の特性的な挙動が全パラメータ空間にわたって存在することはない。一部の高質量構成では合体が見られないことから、大きな$M_\bullet$だけでは効率的な集積を保証しないことが示されている。

3. スピンおよび離心率のシグネチャー:

   • 有効スピン（$\chi_{eff}$）: $\chi_{eff}$の分布は、初期の順行成分（$f_{progr}$）に強く依存する。順行構成は正の$\chi_{eff}$をもたらし、逆行構成は負の値をもたらす。中間的な割合では、二峰性の分布を示す。一次質量と$|\chi_{eff}|$の間には相関が存在し、繰り返しの合体がコヒーレントなスピンの向きを強化する。
   • 離心率: ほとんどの合体は、地上検出器のバンド（$f_{peak} = 10$ Hz）に入る際、低い離心率（$e \sim 10^{-3}–10^{-2}$）を持つ。しかし、有意な円形化が起こる前に合体したサブ集団は、高い離心率（$e \gtrsim 0.8$）を保持している。

4. GWイベントとの比較:

   • AGNチャネルは、GW190521やGW231123と一致する系を再現できるが、これらは予測分布の裾の部分を占める。
   • 低粘性モデルにおいて、GW190521に相当するシステムは$\sim 2 \times 10^{-3}$の割合で発生し、GW231123に似たシステムはより稀である（$\sim 3 \times 10^{-4}$）。
   • 本モデルは、これらの特定のイベントに対して二次質量（$m_2$）を過小評価する傾向があり、これはAGNチャネルが一次質量やスピンの特性を生成できる一方で、成分質量の完全な分布を一致させることは困難であることを示唆している。

5. モデルの堅牢性:

   • 集団の一般的な特性（高質量テール、スピン相関）は、移動トルクの形式（例：熱的補正）やガス硬化の処方箋（IG20 vs. IG24 vs. C24）の変化に対してロバストである。
   • 決定的な依存性: ガス硬化をオフにすると、高質量の残骸形成と高質量テールは完全に消失する。これは、ガス駆動の連星進化がAGN円盤における効率的な階層的成長のための鍵となる成分であることを裏付けている。

意義と主張
本論文は、AGNパラメータ空間を系統的にマッピングし、BBH合体予測の堅牢性を評価するための、包括的かつ最新のフレームワークを提供することを主張している。著者らは、AGNチャネルは単一の特性的な集団を生み出すのではなく、幅広い多様な結果をもたらすと断言している。本研究は、一次質量と$|\chi_{eff}|$の相関、高質量における不等質量比への傾向、およびスピン分布における初期軌道構成の記憶の保持といった、特筆すべき集団レベルのシグネチャーを特定している。

本研究は、AGNチャネルがGW190521やGW231123のような大規模なイベントを説明できる可能性がある一方で、これらのイベントはチャネル内における稀な事象であることを強調している。著者らは、ガス駆動の連星進化が効率的な階層的成長に不可ントであり、円盤構造や硬化に関する具体的な物理的仮定が、高質量残骸形成の効率に大きく影響するものの、高質量テールの存在自体はチャネルのロバストな特徴として残る、と結論づけている。