[论文解读] The erratic path to coalescence of LISA massive black hole binaries in sub-pc resolution simulations of smooth circumnuclear gas disks

Rafael Souza Lima, Lucio Mayer+3 more|arXiv (Cornell University)|Mar 30, 2020

Astrophysical Phenomena and Observations参考文献 1被引用 4

一句话总结

本研究利用高分辨率流体动力学模拟，探究了在平滑的核区气体盘（CND）中大质量黑洞（MBH）双星的轨道衰减过程，揭示了由盘面力矩驱动的不规则衰减行为。关键发现是，当分离尺度小于1 pc时，衰减过程变得低效且脆弱，仅在高分辨率模拟中出现最终的缓慢阶段，从而促成周围双星盘（CBD）的形成，暗示在富含气体的环境中，恒星三体相互作用可能是双星合并所必需的。

ABSTRACT

We report on high-resolution simulations that explore the orbital decay of massive black hole (MBH) pairs with masses between $10^5$ and $10^7 M_{\odot}$ embedded in a circumnuclear gas disk (CND). An adiabatic equation of state is adopted, with a range of adiabatic indices, which maintains a smooth flow. Mergers between MBHs in this mass range would be detectable by the upcoming Laser Inteferometer Space Antenna (LISA). The orbital evolution is followed from the CND scale ($100$~pc) down to separations of $0.1$--$0.01$~pc at which a circumbinary disk (CBD) could form. The decay is erratic and strongly dependent on the gas flow within the disk, that ultimately determines the net torques experienced by the sinking MBH. Overall, we can identify three different evolutionary stages: (i) an initially slow decay that leads to no significant change in the orbital angular momentum, resulting in some circularization; (ii) a fast migration phase in which the orbital angular momentum decreases rapidly; and (iii) a final, very slow decay phase, in which orbital angular momentum can even increase, and a CBD can form. The fast migration phase owes to disk-driven torques originating primarily from the co-orbital region of the secondary MBH, at a distance of 1--3 Hill radii. We find strong analogies with fast Type III migration for massive planets in protoplanetary disks. The CBD forms only when the decay rate becomes small enough to allow it enough time to carve a cavity around the primary MBH, at scales $\lesssim 1$~pc; when this happens, the MBH separation nearly stalls in our higher-resolution run. We suggest an empirically modified gap-opening criterion that takes into account such timescale effects as well as other deviations from standard assumptions made in the literature. [Abriged]

研究动机与目标

• 研究从约100 pc至0.01 pc的核区气体盘（CND）中大质量黑洞（MBH）双星的轨道衰减过程。
• 确定气体盘结构与分辨率如何影响亚pc分离尺度下MBH双星迁移的效率与机制。
• 评估形成周围双星盘（CBD）的条件，这是引力波辐射的关键前兆。
• 评估标准间隙开启准则在高分辨率、平滑CND模拟中的局限性，并提出经验性修正准则。
• 探讨对‘最终光年问题’的启示，以及三体恒星相遇等额外硬化机制的作用。

提出的方法

• 在质量介于10⁵至10⁷ M☉之间的平滑、压力支撑的CND中，对MBH双星进行高分辨率光滑粒子流体动力学（SPH）模拟。
• 采用可变绝热指数的绝热状态方程，以维持平滑的气流并控制盘面温度与压力。
• 追踪从约100 pc至0.01 pc的轨道演化，重点关注力矩驱动迁移与角动量演化。
• 引入一个现象学上的修正参数g，以经验方式调整标准间隙开启准则，考虑迁移 timescales 与非理想条件。
• 通过比较高分辨率与低分辨率模拟结果，分离出分辨率对空腔与CBD形成的影响。
• 分析不同盘区（特别是1–3倍希尔半径范围内的共轨区域）的力矩贡献，以识别快速迁移的起源。

实验结果

研究问题

• RQ1 在平滑CND中，MBH双星的轨道衰减从约100 pc至亚pc尺度如何演化？不同阶段的衰减由何种机制驱动？
• RQ2 在何种条件下会形成周围双星盘（CBD）？分辨率如何影响该过程？
• RQ3 为何在亚pc分离尺度下衰减速率变得低效且不规则？导致最终减速或停滞的物理机制是什么？
• RQ4 由次级MBH共轨区域产生的盘面力矩如何促成快速迁移阶段？
• RQ5 标准间隙开启准则能否可靠地应用于高分辨率、平滑CND模拟中？还是必须进行经验性调整？

主要发现

• 轨道衰减经历三个不同阶段：初始缓慢阶段（伴随轨道圆化）、由1–3倍希尔半径范围内的共轨区域力矩驱动的快速迁移阶段，以及角动量甚至可能增加的最终缓慢阶段。
• 快速迁移阶段类似于III型行星迁移，主要由次级MBH共轨区域的气体力矩驱动。
• 仅在高分辨率模拟中，当衰减速率足够慢以允许在≤1 pc尺度上形成主MBH周围的空腔时，才会形成周围双星盘（CBD）。
• 在低分辨率模拟中，快速衰减抑制了空腔与CBD的形成，导致轨道衰减更快且无停滞。
• 由于timescale效应，标准间隙开启准则在本模拟中失效；引入修正参数g的经验性修正准则能成功预测当迁移变慢时CBD的形成。
• 即使存在CBD与空腔，MBH双星仍保持非零偏心率，且在亚pc尺度上，恒星三体相互作用仍是可行的硬化机制。