Gas-induced perturbations on the gravitational wave in-spiral of live post-Newtonian LISA massive black hole binaries: 0.1 disk aspect ratio

技術要約：厚い円盤を持つライブなポスト・ニュートン近似LISA質量ブラックホール連星の重力波インスパイラルにおけるガス誘起摂動

問題提起
レーザー干渉計宇宙アンテナ（LISA）は、特に合流付近のミリヘルツ帯における質量ブラックホール連星（MBHB）を検出することが期待されている。極めて重要な未解決の問いは、周囲の円盤（CBD）が重力波（GW）の波形に検出可能な痕跡を残すかどうかである。従来の研究では、半解析的に導出されたガス・トルクをGW駆動の進化に線形的に加算することで、ガスの影響をモデル化することが多かった。このアプローチは、ガス力学とGW放射の間の非線形な結合を見落とす可能性があり、両者のデカップリング（分離）を前提としている。さらに、薄い円盤（$H/R = 0.03$）の研究は進んでいるが、計算効率が高い（粘性時間が短い）厚い円盤（$H/R = 0.1$）が、結合したGW-ガス進化に与える影響については、まだ十分に探索されていない。

手法
著者らは、GIZMOコードを用いた3次元流体（HD）シミュレーションを実行し、前進回転の局所等温的なCBDに埋め込まれた、等質量（総質量 $10^6 M_\odot$）の準円軌道MBHBをモデル化した。

• 連星の進化: 連星は、ガストルクとGW放射の両方の下で進化する「ライブ」な系として扱われる。ダイナミクスは、2.5ポスト・ニュートン（PN）補正によって記述され、この2.5PN項がGW誘起のインスパイラルを駆動する。
• シミュレーションの設定: 以下の4つの構成を比較する：
  1. Gas + 2.5PN: 完全な結合進化（GW + ガス）。
  2. Gas + 2PN: ガストルクは存在するが、GW放射は無視される（2.5PN項を除去）。
  3. 2.5PN: 真空中でのGW駆動進化。
  4. 2PN: GWのない真空中の進化。
• 初期条件: シミュレーションは、赤方偏移 $z \sim 1$ における合流の約1年前に相当する、連星間隔 $54.5 r_s$（シュワルツシルト半径）から開始される。系は600軌道周期（$P_B$）の間進化し、$46.5 r_s$ まで縮小する。
• 解像度: 解像度スタディにより、中間解像度の実行（$3a$ において $\Delta x \approx 2.1 r_s$）が、トルクおよび位相シフトの測定に十分な収束性を備えていることが確認された。

主な貢献

1. $H/R=0.1$ における初の測定: 本研究は、$H/R = 0.1$ のCBD内にあるライブなMBHBにおける重力トルクおよび降着トルクの初の測定結果を報告しており、GW放射が同時に発生する場合としない場合を比較している。
2. 結合したGW-ガス力学: 従来の半解析的なアプローチとは異なり、本研究はGW駆動のインスパイラルとガス力学の間の非線形な結合を直接シミュレートしており、その結果生じる軌道位相シフトの精密な測定を可能にしている。
3. 降着形態: GW放射が支配的な進化メカニズムとなった際に、降着時系列の形態がどのように変化するかを分析している。

結果

• トルク測定: 無次元トルクパラメータ（$\xi$）を600軌道にわたって測定した。
  • 降着トルク ($\xi_{acc}$): $\sim 0.44$ (GW+gas) 対 $\sim 0.43$ (gas-only)。
  • 重力トルク ($\xi_{grav}$): $\sim 0.13$ (GW+gas) 対 $\sim 0.16$ (gas-only)。
  • 全トルク: 両方のケースにおいて純トルクは正であり、これはガスがインスパイラルを遅らせることを示している。著者らは、$H/R=0.1$ の場合、降着トルクの大きさは重力トルクの約2倍であり、重力トルクが支配的となる薄い円盤の結果とは対照的であると指摘している。
  • デカップリング: CBDは連星のインスパイラルに従うが、外側円盤（$R > a$）からの重力トルクは、連星の縮小に伴いわずかに減少しており、これはデカップリングの始まりと一致している。
• 降着率の変調: GWの存在は、降着率の形態をわずかに変化させる。スペクトログラムは、空洞の内縁変調（$\sim 3\%$）および連星軌道周期変調（時間の経過とともに$15\%$から$\sim 60\%$へ増加）におけるパワーの減少を示しており、これは円盤が連星から漸次デカップリングしていくことを示している。
• 位相シフト:
  • GW放射を伴うガス誘起の軌道位相シフト（$\delta\phi_{orb}^{(GW)}$）は、600サイクルで 0.12 rad と測定された。
  • GW放射がない場合の位相シフト（$\delta\phi_{orb}^{(NoGW)}$）は 0.08 rad である。
  • 結果としてのGW位相シフト（$\delta\phi_{GW} = 2\delta\phi_{orb}^{(GW)}$）は 0.24 rad である。
  • 解析的比較: 標準的な半解析的手法（GWとガスの率を線形に合算する手法）では、位相シフトは約 $\sim 0.002$ rad となる。数値計算の結果は、この解析的な予測よりも約60倍大きく、非線形なGW-ガス結合の重要性を浮き彫りにしている。
• 検出可能性: $z=1$ で信号対雑音比（SNR）が $\sim 1300$ の光源について、計算された位相シフト 0.24 rad は検出閾値（$8/\text{SNR} \approx 0.006$ rad）を超えており、LISAがこの環境的痕跡を検出できる可能性を示唆している。

意義および主張
著者らは、彼らの結果がマルチメッセンジャー天文学に直接的な影響を与えることを主張している。光学サーベイ（例：LSST/Roman）による降着率変調の検出と、LISAによる位相シフトの検出が組み合わさることで、MBHBを取り巻く複雑な環境に関する極めて重要な情報が得られる。

本論文は、数値的な位相シフトが半解析的な予測よりも大幅に大きいという事実は、単にガストルクをGW進化に線形的に加算する手法では不十分であることを強調している。本研究は、数百の軌道が観測可能な分離距離において、LISAソースの予測を正しく導き出すために、同時的なGW放射を流体シミュレーションに含めることが必要であることを示すベンチマークとして機能する。著者らは、本アプローチが定常状態の初期条件とシャクラ・スニヤ・ヴィスコーシティ（磁気流体力学的乱流は無視）を仮定しているものの、LISAソースの予測を推論するための最も実現可能な方法であると述べている。