[論文レビュー] Hydrodynamical Evolution of Black-Hole Binaries Embedded in AGN Discs

Rixin Li, Dong Lai|arXiv (Cornell University)|Feb 15, 2022

Astrophysical Phenomena and Observations被引用数 4

ひとこと要約

本研究では、活発な銀河中心核（AGN）円盤に埋め込まれた星質量ブラックホール二重星（BBH）の進化を、せん断ボックスフレームワークを用いた高解像度2次元流体ダイナミクスシミュレーションでモデル化する。円形で質量比が近いBBHは、質量倍加率の数倍のスケールで急速に収縮する（硬化する）が、離心率を有するBBHは離心率の減衰を経験する。すべての現実的なパラメータにおいて、軌道的減衰は移動より速く、降着は著しく変動的であり、離心率が増加するに従い、変動周波数は二重星の軌道周波数から半径方向の振動周波数にシフトする。

ABSTRACT

Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs are possible progenitors of black-hole mergers detected in gravitational waves by LIGO/VIRGO. To better understand the hydrodynamical evolution of BBHs interacting with the disc gas, we perform a suite of high-resolution 2D simulations of binaries in local disc (shearing-box) models, considering various binary mass ratios, eccentricities and background disc properties. We use the $γ$-law equation of state and adopt a robust post-processing treatment to evaluate the mass accretion rate, torque and energy transfer rate on the binary to determine its long-term orbital evolution. We find that circular comparable-mass binaries contract, with an orbital decay rate of a few times the mass doubling rate. Eccentric binaries always experience eccentricity damping. Prograde binaries with higher eccentricities or smaller mass ratios generally have slower orbital decay rates, with some extreme cases exhibiting orbital expansion. The averaged binary mass accretion rate depends on the physical size of the accretor. The accretion flows are highly variable, and the dominant variability frequency is the apparent binary orbital frequency (in the rotating frame around the central massive BH) for circular binaries but gradually shifts to the radial epicyclic frequency as the binary eccentricity increases. Our findings demonstrate that the dynamics of BBHs embedded in AGN discs is quite different from that of isolated binaries in their own circumbinary discs. Furthermore, our results suggest that the hardening timescales of the binaries are much shorter than their migration timescales in the disc, for all reasonable binary and disc parameters.

研究の動機と目的

• LIGO/Virgo重力波合体の可能性のある形成経路である、ガス状AGN円盤に埋め込まれた星質量ブラックホール二重星（BBH）の流体的進化を理解すること。
• 数値的解像度や降着モデルの違いに起因する、従来の研究で報告されたBBHがAGN円盤内で硬化するか拡張するかという矛盾する結果を解消すること。
• BBHの質量比、離心率、円盤特性が降着率、トルク、軌道的進化に与える影響を定量化すること。
• AGN円盤内BBHの硬化 timescale が移動 timescale より短いかどうかを特定すること。これは合体形成効率の鍵となる重要な問いである。

提案手法

• 中心の超大質量ブラックホールの周囲のAGN円盤の局所的領域をモデル化するため、せん断ボックス近似を用いた高解像度2次元流体ダイナミクスシミュレーションを実施する。
• γ法則状態方程式を採用し、後処理手法を用いてBBHに対する質量降着率、トルク、エネルギー移動率を正確に計算する。
• 重力ソフトニング長を二重星間隔のおよそ0.08に設定し、各二重星成分の周囲に形成される循環単一円盤（CSD）を適切に解像する高空間解像度を確保する。
• 質量比（q = 0.1 から 1.0）、離心率（e = 0.0 から 0.5）、円盤パラメータ（h = 0.01, λ = 2.5）を変化させたシミュレーションを実施し、最終300軌道周期にわたる時間平均量を算出する。
• 半長径、離心率、二重星の固有エネルギーの時間平均変化率を用いて軌道的進化を評価する。
• 一貫した後処理フレームワークを適用し、シミュレーション出力から降着およびトルク率を抽出することで、長期的進化の信頼性の高い推定を実現する。

実験結果

リサーチクエスチョン

• RQ1 AGN円盤内BBHの軌道的減衰率は、その移動 timescale を上回るか。どのような条件下でそうなるか？
• RQ2 BBHの質量比は、AGN円盤内での降着率および軌道的進化にどのように影響を与えるか？
• RQ3 二重星の離心率が、降着の主な変動周波数の特定および離心率減衰率に果たす役割は何か？
• RQ4 適切に解像されたCSDを有する高解像度シミュレーションは、低解像度研究と比較して、BBHがAGN円盤内で拡張するという従来の報告をどのように変えるか？
• RQ5 循環単一円盤（CSD）の熱力学的状態（例：温度）は、ネットトルクに影響を及ぼし、結果として二重星の硬化または拡張に寄与するか？

主な発見

• 円形で質量比が近いBBHは、質量倍加率の数倍のスケールで急速に収縮し、速やかな硬化を示す。
• 離心率を有する二重星は常に離心率の減衰を経験するが、より高い離心率またはより小さい質量比の前向き二重星では、軌道的減衰が遅くなり、一部のケースでは軌道の拡張が観測される。
• 時間平均降着率は降着体の物理的サイズに依存し、より大きな吸収半径ではより高い降着率を示す。
• 降着変動は、円形二重星では回転フレーム内での二重星の軌道周波数に支配的であるが、離心率が増加するに従い、半径方向振動周波数にシフトする。
• すべての現実的なパrameterの組み合わせにおいて、AGN円盤内BBHの硬化 timescale はその移動 timescale よりも著しく短く、急速な合体形成を促進する。
• 結果は、AGN円盤内BBHの力学的挙動が、循環二重星円盤内に存在する孤立した二重星とは本質的に異なることを示しており、より強いガストルクと特徴的な降着シグナルを有する。