論文要約：埋め込まれた連星への自己制限的な降着（Self-Limited Accretion onto Embedded Binaries in a Uniform Medium）

著者: Marcus DuPont, Eliot Quataert
日付: 2026 年 3 月 19 日（ドラフト版）
対象: 天体物理学、降着流、連星ブラックホール

1. 研究の背景と問題設定

本研究は、静止した一様なガス中にある等質量の連星（Binary Bondi 問題）への降着効率を、断熱指数 $\gamma$ とコンパクトネス $\xi \equiv R_B/a$（$R_B$ は連星のボンディ半径、$a$ は連星間距離）の関数として解析したものである。

従来の期待では、コンパクトネスが大きい場合（$\xi \gg 1$）、連星は共有されたボンディ球内に深く埋め込まれ、単一の天体と同様に協調的に降着し、効率 $\eta \equiv \dot{M}_{\text{binary}}/\dot{M}_{\text{Bondi}}$ が 1 に近づくと考えられていた。しかし、この論文は、$\gamma > 1$（断熱的ガス）の場合、この期待が破綻することを示している。連星の軌道運動が衝撃波を駆動し、エントロピーを生成することで対流乱流が発生し、降着を抑制する「自己制限（self-limiting）」メカニズムが働くことを発見した。

2. 手法と数値シミュレーション

• 数値コード: 公開コード「SIMBI」を用いた、3 次元非粘性流体力学シミュレーション。
• 物理条件:
  • 断熱指数 $\gamma = 1$（等温）、$4/3$（相対論的）、$5/3$（単原子）。
  • コンパクトネス $\xi = \{0.1, 1, 10\}$。
  • 質量比 $q=1$（等質量）。
• 数値手法:
  • 広範囲のスケール（$r \sim a$ から $r \sim R_B$）を解像するため、固定メッシュリファインメント（FMR）を用いたテレスコーピンググリッドを採用。
  • 境界条件には、解析的なボンディ解へ緩和する「緩和境界条件」を使用し、人工的な反射や希薄化を防ぎ、無限の質量供給を模擬した。
  • 降着効率 $\eta$ の収束性を検証し、シンク（質量除去）のサイズ依存性が比値として相殺されることを確認。

3. 主要な理論的貢献：安定性基準の導出

著者は第一原理から、連星がボンディの鞍点（saddle point）を破壊し、降着流を乱すための条件を導出した。

1. 幾何学的閾値 ($\xi_m$):
   個々のコンポーネントの音速面が接触し、連星全体を囲む共有の音速面が形成される臨界コンパクトネスは、$\xi_m = 4/(5-3\gamma)$ で与えられる。

   • $\gamma=1$ の場合、$\xi_m = 4$。
   • $\gamma=5/3$ の場合、$\xi_m \to \infty$（分離した音速面は存在しない）。

2. 熱力学的不安定性:
   連星の軌道運動は、回転座標系において超音速でガスを掃き抜ける「プロペラ」として機能し、衝撃波を生成する。

   • 等温ガス ($\gamma=1$): 衝撃波によるエントロピー生成はゼロであり、降着は協調的に進行し、$\eta \to 1$ となる。
   • 断熱ガス ($\gamma > 1$): 衝撃波で生成されたエントロピーが連星を取り囲む領域に蓄積し、シュワルツシルト安定性基準（$ds/dr < 0$）を破る。これにより対流乱流が発生し、降着を抑制する。

3. 破壊までの軌道数 ($N$):
   連星が自身の燃料供給を遮断する（鞍点を消滅させる）までに必要な軌道数 $N$ は以下のように推定される：
   $$N \propto \frac{1}{\gamma - 1} \left( \sqrt{\frac{\xi}{\xi_m}} - 1 \right)$$

   • $\gamma \to 1$ で $N \to \infty$（等温では不安定化しない）。
   • $\gamma \to 5/3$ で $\xi_m \to \infty$ となり、この式は適用外となるが、$\gamma=5/3$ では中立成層のため、いかなる $\xi$ でも対流不安定となる。

4. 冷却時間の役割:
   実際の天体物理系において、降着が協調的か乱流的かは、冷却時間 $t_{\text{cool}}$ と軌道周期 $T$ の積 $NT$ との比較で決まる。

   • $t_{\text{cool}} \ll NT$: 等温挙動（協調的降着）。
   • $t_{\text{cool}} \gg NT$: 断熱挙動（自己制限・乱流）。

4. シミュレーション結果

4.1 流れの形態

• $\gamma = 1$ (等温): 全ての $\xi$ で定常的かつ組織化された流れ。$\xi=10$ では、2 つの降着源を繋ぐ高密度の「橋（bridge）」が形成され、効率的な協調降着を実現する。
• $\gamma = 4/3$ (相対論的): $\xi < \xi_m$ では秩序だった螺旋構造が見られるが、$\xi > \xi_m$（特に $\xi=10$）では、衝撃波によるエントロピー蓄積が音速面を破壊し、完全な乱流状態となる。
• $\gamma = 5/3$ (単原子): 全ての $\xi$ で乱流化。分離した音速面が存在しないため、外側の流れも連星の重力攪拌の影響を直接受け、大規模な秩序構造は存在しない。

4.2 降着効率 ($\eta$)

• $\gamma = 1$: $\xi$ が増加するにつれて $\eta$ は単調増加し、$\xi=10$ で $\eta \approx 1$ に到達。
• $\gamma = 4/3$: 低 $\xi$ では独立降着に近いが、$\xi > \xi_m$ 以降、$\eta$ は急激に低下し、$\xi=10$ で $\eta \approx 0.3$ まで減少。乱流が降着を独立降着率以下に抑制する。
• $\gamma = 5/3$: 全ての $\xi$ で $\eta$ が抑制され、$\xi=10$ で $\eta \approx 0.1$ 程度まで低下。

4.3 時間変動と軌道信号

• $\gamma = 1$: 軌道周期の倍音（$4\Omega, 8\Omega$）に明確なピークが見られる（非線形応答による）。
• $\gamma = 4/3$: 軌道信号（$2\Omega$）は残存するが、乱流ノイズに埋もれ、検出が困難になる。
• $\gamma = 5/3$: 軌道信号は完全に広帯域ノイズに埋没し、検出不可能となる。

4.4 乱流特性

• エントロピーの角力率スペクトルは、$\gamma=4/3, 5/3$ ともに $k^{-3.5}$（$k^{-7/2}$）の傾きを示す。これはコルモゴロフ ($k^{-5/3}$) やバーガース ($k^{-2}$) 乱流よりも急峻であり、孤立した衝撃波面や対流境界によって支配されていることを示唆する。

5. 天体物理学的意義と結論

1. 降着効率の過大評価の是正:
   従来の連星ブラックホール（SMBH）の合体シナリオや、共通包絡星段階のモデルでは、層流的なボンディ降着率に基づいてガスによる軌道硬化（hardening）を推定してきた。しかし、断熱的な環境（冷却時間が長い場合）では、降着効率が 3〜10 倍も過大評価されていた可能性がある。これは「最終パーセク問題（Final Parsec Problem）」の解決策として、熱力学的なメカニズムが重要な役割を果たすことを示唆する。

2. 周期性 AGN の欠如の理由:
   連星 AGN の周期的な光度変動を検出する試みが失敗してきた理由の一つとして、断熱環境では軌道信号が乱流ノイズに埋もれる（$\gamma=5/3$）か、あるいは等温環境では信号が倍音（$4\Omega$）にシフトするため、従来の $2\Omega$ 探索では検出されなかった可能性が示された。

3. 対流支配降着流（CDAF）の新たな経路:
   粘性や磁場を介さず、連星の軌道運動による衝撃波エントロピー生成が、対流を駆動し降着を抑制する「第三の経路」として確立された。

結論:
埋め込まれた連星が降着するか、それとも自身の燃料供給を遮断して乱流化するかは、冷却時間 $t_{\text{cool}}$ と軌道周期 $T$ の積 $NT$ との比較によって決定される。冷却が速い系は協調的に成長するが、冷却が遅い断熱系は自己制限され、降着効率が大幅に低下し、軌道信号も消滅する。この単一の比較が、天体物理学的な連星の進化と観測可能性を支配する。