[論文レビュー] Merging of unequal mass binary black holes in non-axisymmetric galactic nuclei

Peter Berczik, Manuel Arca Sedda+4 more|arXiv (Cornell University)|Aug 10, 2020

Galaxies: Formation, Evolution, Phenomena参考文献 55被引用数 8

ひとこと要約

本研究は、ϕ-GPUコードを用いた高精度N体シミュレーションを用いて、非軸対称で回転する銀河核における非等質量の超大質量黒 holes二重星（SMBHBs）の力学的硬化と合体を調査している。著者らは、著しく停滞しない（「最終光年問題」は顕著でない）と判明し、非対称な二重星は等質量の二重星よりも最大10,000倍速く合体することが判明した。硬化は回転する系における星の散乱によって効率的に駆動されている。

ABSTRACT

In this work, we study the stellar-dynamical hardening of unequal mass supermassive black hole (SMBH) binaries in the central regions of merging galactic nuclei. We present a comprehensive set of direct $N$-body simulations of the problem, varying both the total mass and the mass ratio of the SMBH binary (SMBHB). Simulations were carried out with the $\varphi-$GPU $N$-body code, which enabled us to fully exploit supercomputers equipped with graphic processing units (GPUs). As a model for the galactic nuclei, we adopted initial axisymmetric, rotating models, aimed at reproducing the properties of a galactic nucleus emerging from a galaxy merger event, containing two SMBHs which were unbound initially. We found no 'final-parsec problem', as our SMBHs tend to pair and shrink without showing significant signs of stalling. This confirms earlier results and extends them to large particle numbers and rotating systems. We find that the SMBHB hardening depends on the binary-reduced mass ratio via a single parameter function. Our results suggest that, at a fixed value for the SMBHB primary mass, the merger time of highly asymmetric binaries is up to four order of magnitudes smaller than the equal-mass binaries. This can significantly affect the population of SMBHs potentially detectable as gravitational wave sources.

研究の動機と目的

• 銀河合体に起因する非軸対称で回転する銀河核における非等質量超大質量黒 holes二重星（SMBHBs）の星間力学的硬化を調査すること。
• 従来のシミュレーションを拡張し、質量比（q = 0.01～1）の広い範囲を探索するとともに、銀河核モデルに回転を組み込むこと。
• 大きな質量非対称性を示す回転する非軸対称系においても、「最終光年問題」——硬化の潜在的停滞——が継続するかどうかを検証すること。
• 合体時定数が二重星の質量比と全質量にどのように依存するかを定量化し、重力波検出への影響を評価すること。
• 全質量と換算質量を含む二パラメータスケーリング則の適用範囲を検証し、拡張すること。

提案手法

• ϕ-GPU N体コードを用いた直接N体シミュレーションを実施し、GPU加速を活用して最大10^6粒子の大きな粒子系をシミュレートした。
• 銀河核を初期的に軸対称で回転する系としてモデル化し、銀河合体の結果として生じる二つの非束縛なSMBHを再現した。
• 最終段階の進化において重力波放射効果を自己無撞着に含めるために、運動方程式に後ニュートン補正を組み込んだ。
• 全質量（10^4–10^9 M⊙）と質量比（q = 0.01–1）を複数のシミュレーションで変化させ、硬化時定数をマップした。
• 時間経過に伴う二重星の分離距離、束縛エネルギー、硬化速度を追跡し、系のパラメータに依存する合体時定数T_mergeを導出した。
• 換算質量µと全質量M12を用いた次元なしの硬化速度のスケーリングにより、普遍的なスケーリング関係を同定した。

実験結果

リサーチクエスチョン

• RQ1 非回転で非軸対称な銀河核に存在する非等質量SMBHBにおいて、「最終光年問題」は継続するか？
• RQ2 質量比（q = m2/m1）は、密度が高く回転する星間系におけるSMBHBの硬化速度と合体時定数にどのように影響するか？
• RQ3 銀河核内に回転が存在する場合、SMBHBの星間力学的硬化は促進されるか、あるいは抑制されるか？
• RQ4 合体時定数は、全質量M12と換算質量µを含む二パラメータスケーリング則で記述可能か？
• RQ5 極めて非対称なSMBHB（例：q = 0.01）の予想される合体時定数は、等質量二重星（q = 1）と比較してどの程度短いか？また、重力波源としての検出可能性にどのように影響するか？

主な発見

• シミュレーションでは著しい停滞（すなわち、「最終光年問題」）は観測されず、回転する非軸対称核内ではSMBHBが効率的に硬化し、合体を継続した。
• 主SMBH質量が一定の場合、二次質量比qが1から0.01に減少するに従い、合体時定数T_mergeは最大4桁短縮された。
• 硬化速度は換算質量µと全質量M12に依存し、異なる質量比にわたるデータを統合する二パラメータスケーリング則が支持された。
• 10^4–10^9 M⊙の成分質量と質量比q = 10^{-4}–1を有するSMBHBは、星間力学的硬化によりT_merge ≈ 10^3–10^7年という時定数で合体した。
• 極めて非対称な二重星（q ≈ 0.01）の合体時定数は、等質量二重星（q = 1）と比較して最大10,000倍短く、重力波源としての検出可能性が著しく向上した。
• 本結果は、より大きな粒子数と回転する系を用いた従来の発見を確認・拡張し、星間硬化が非対称で回転する核においても効果的であることを示した。