Perturber-Driven Dynamics of Supermassive Black Hole Binaries in Galaxy Merger

技術要約：銀河合体における摂動体駆動型の超大質量ブラックホール連星ダイナミクス

問題提起
超大質量ブラックホール連星（MBHB）の形成時における軌道離心率は、パルサー・タイミング・アレイ（PTA）によって検出可能な重力波背景放射（GWB）を決定する極めて重要なパラメータである。これまでの高解像度$N$体シミュレーションは、形成時の離心率が、特にほぼ径方向の合体軌道において顕著なラン・トゥ・ランのばらつきを示すことを明らかにしている。このばらつきは、解像度とともに$\sim 1/\sqrt{N}$に従って減少することが示されており、これはポアソン・サンプリング・ノイズと一致しているが、物理的なサブストラクチャ（具体的には、球状星団、剥ぎ取られた核、あるいは中間質量ブラックホール（IMBH）などの質量を持つ摂動体）が、数値的なノイズフロアを超える真の天体物理学的確率性を導入するかどうかは不明である。中心となる問いは、多数の低質量摂動体の累積効果（拡散レジーム）か、あるいは稀に発生する大規模な遭遇（衝撃レジーム）のどちらかが、大質量楕円銀河におけるMBHBの形成離心率を著しくランダム化できるのかという点である。

手法
著者らは、$Griffin$ $N$-bodyコードを用いて、IllustrisTNG100-1宇宙論的シミュレーションの主要な銀河合体（Merger 6）の高解像度再シミュレーションを実施した。シミュレーションの設定は以下の通りである：

• 初期条件： Hernquistプロファイルに従う衝突レスな系（星のバルジ、ダークマターハロー、および中心のMBH）としてモデル化された2つの銀河間の、ほぼ径方向の合体（$e_{orb} \approx 0.987$）。主銀河は小さなガス分率（$f_{gas} \approx 0.19$）を保持しているが、長寿命の恒星系に焦点を当てるため、本シミュレーションではシステムを衝突レスとして扱った。
• 摂動体シナリオ： 以下の3つのシミュレーション・スイートを比較した。各スイートには4つの独立した実現（realisations）が含まれる：
  1. コントロール： 大規模な摂動体なし。滑らかな星のバルジが背景として機能する。
  2. $10^7 M_\odot$ スイート： $3,162$個の等質量摂動体（$M_p = 10^7 M_\odot$）を用い、主バルジ質量の$10%$を再分配する。これは、摂動体対連星の質量比$\mu_p \approx 3.2 \times 10^{-3}$ に相当する。
  3. $10^8 M_\odot$ スイート： $316$個の等質量摂動体（$M_p = 10^8 M_\odot$）を用い、同じ総質量分率を再分配する。これは、$\mu_p \approx 3.2 \times 10^{-2}$ に相当する。
• 連星形成基準： 連星形成時間（$t_b$）は、永続的な重力結合後の離心率の最初の局所的最小値として定義された。形成離心率（$e_b$）およびその標準偏差（$\sigma_e$）が主要な指標となった。
• 動力学的診断： 本研究では、$t_b$に合わせて整列させた、軌道エネルギー（$E$）、角運動量（$J_z$）、トルク（$\tau_z$）の残差、および連星重心（COM）の揺らぎを分析した。

主な貢献と結果
本論文は、摂動体対連星の質量比（$\mu_p$）に基づく、拡散散乱レジームと近衝撃散乱レジームの間の遷移を確立している：

1. 数値ノイズフロア： 摂動体なしのコントロール・スイートは、採用された解像度（$N_{bulge} \approx 2 \times 10^6$）に対して期待されるポアソン・ノイズ・フロアと一致する、$\sigma_e \approx 0.11$ の形成離心率のばらつきを示した。
2. 拡散レジーム（$10^7 M_\odot$）： $10^7 M_\odot$ スイートは $\sigma_e \approx 0.115$ を生成し、これはコントロールと統計的に区別がつかない。エネルギーおよび角運動量の動力学的残差は滑らかであり、有意なトルクのスパイクも観察されなかった。これは、$\mu_p \approx 3.2 \times 10^{-3}$ において、システムが拡散レジームに留まり、個々の遭遇が形成離心率に識別可能な痕跡を残すには弱すぎることを裏付けている。
3. 近衝撃レジーム（$10^8 M_\odot$）： $10^8 M_\odot$ スイートは、有意に大きなばらつき $\sigma_e \approx 0.26$（コントロールに対し約2.4倍）を示した。この過剰なばらつきは、以下の現象を伴っていた：
   • 半長軸および軌道エネルギーの不連続なジャンプ。
   • $t_b$付近における、イベントに整列した角運動量の大きな残差およびトルクのスパイク。
   • 連星重心の揺らぎの増大（コントロールより$\sim 1.25$–$1.4$倍高い）。
     これらのシグネチャーは、このレジームにおける個別の近接遭遇が、硬化フェーズが始まる前に連星の角運動量と離心率を変化させるのに十分な衝撃的なキックを与えることを示している。
4. 理論的遷移： 結果は、$\mu_p$によって支配される連星・単体散乱理論と一致している。$10^8 M_\odot$ のケースは近衝撃限界（$\Delta e_{single} \sim \mu_p$）に近づいており、$10^7 M_\odot$ のケースは拡散レジームに留まっている。

意義と主張
本論文は、大規模な摂動体はMBHBの形成離心率に物理的な確率性を誘発し得るが、この効果は摂動体の質量スケールが連星に対してどの程度であるかに強く依存すると主張している。

• 天体物理学的文脈： （PTAに関連するMBHBの主要な宿主である）大質量楕円銀河において、支配的な摂動体集団は球状星団（$10^5$–$10^6 M_\odot$）や剥ぎ取られた核である。PTAに関連する連星質量（$10^8$–$10^{10} M_\odot$）において、これらの天体は $\mu_p \lesssim 10^{-4}$ であり、拡散レジームに明確に位置している。
• 結論： したがって、摂動体による離心率のランダム化は、GWB検出に関連する大多数のMBHB合体にとって、重要な経路となる可能性は低い。現在のシミュレーションで見られる確率性は、現実的な摂動体集団による物理的なサブストラクチャではなく、数値的なノイズによって支配されている。
• 注意点： 著者らは、彼らの結論が（単一の $10^8 M_\odot$ の増大の統計的有意性ではなく）質量比の余裕に基づいていることを注記している（わずか4回の実現によるため）。また、ガスが豊富な合体（ここでは完全にはモデル化されていない）は、巨大分子雲や円盤状の周囲円盤を介して異なるダイナミクスをもたらす可能性があること、および特定の軌道構成（ほぼ径方向）が正確な閾値に影響を与える可能性があることも認めている。しかし、主要な天体物理学的結論は、現実的な摂動体集団は、形成離心率を著しくランダム化するために必要な近衝撃閾値を下回っているということである。