[论文解读] A new Monte Carlo code for star cluster simulations: II. Central black hole and stellar collisions

Marc Freitag, W. Benz|ArXiv.org|Apr 17, 2002

Astrophysical Phenomena and Observations参考文献 149被引用 50

一句话总结

本文提出了一种新的蒙特卡洛代码，用于模拟致密星系核，结合了不断增长的中心黑洞和基于10,000次SPH模拟的真实恒星碰撞。该代码可在普通个人电脑上实现百万粒子、宇宙学时标（Hubble timescale）的模拟，准确建模潮汐撕裂、碰撞驱动的质量转移以及恒星演化，以研究星系中心黑洞的增长和奇特恒星种群的形成。

ABSTRACT

We have recently written a new code to simulate the long term evolution of spherical clusters of stars. It is based on the pioneering Monte Carlo scheme proposed by Henon in the 70's. Our code has been devised in the specific goal to treat dense galactic nuclei. After having described how we treat relaxation in a first paper, we go on and include further physical ingredients that are mostly pertinent to galactic nuclei, namely the presence of a central (growing) black hole (BH) and collisions between MS stars. Stars that venture too close to the BH are destroyed by the tidal field. This process is a channel to feed the BH and a way to produce accretion flares. Collisions between stars have often been proposed as another mechanism to drive stellar matter into the central BH. To get the best handle on the role of this process in galactic nuclei, we include it with unpreceded realism through the use of a set of more than 10000 collision simulations carried out with a SPH (Smoothed Particle Hydrodynamics) code. Stellar evolution has also been introduced in a simple way, similar to what has been done in previous dynamical simulations of galactic nuclei. To ensure that this physics is correctly simulated, we realized a variety of tests whose results are reported here. This unique code, featuring most important physical processes, allows million particle simulations, spanning a Hubble time, in a few CPU days on standard personal computers and provides a wealth of data only rivalized by N-body simulations.

研究动机与目标

• 开发一种高保真蒙特卡洛代码，用于模拟致密星系核的长期演化，包含以往实现中缺失的关键物理过程。
• 利用10,000次SPH模拟，实现真实恒星碰撞的建模，以捕捉碰撞结果和向中心黑洞的质量转移。
• 建模中心黑洞增长过程中对恒星的潮汐撕裂，同时考虑安静星系核中吸积与耀斑产生的双重效应。
• 包含简化的恒星演化模型，并通过全面的测试模拟确保数值精度。
• 实现在普通个人电脑上支持高达百万粒子、长达宇宙时标（Hubble time）的大规模长期模拟。

提出的方法

• 该代码在Hénon于1970年代提出的蒙特卡洛框架基础上扩展，现包含双体松弛、恒星碰撞和潮汐撕裂过程。
• 恒星碰撞通过10,000次SPH模拟生成的数据库进行建模，以捕捉真实的流体动力学结果和质量转移。
• 潮汐撕裂通过识别轨道近日点位于中心黑洞潮汐半径内的恒星来实现，能量与动量转移通过轨道积分计算。
• 初始条件基于球对称自洽星系统分布函数（DF）的相空间分布函数生成，径向与速度分布通过拒绝采样方法获取。
• 通过重缩放速度以确保维里平衡，代码支持通过分段幂律初始质量函数（IMF）实现可变质量函数。
• 代码使用N体单位，并可通过星团质量与大小的缩放转换为物理单位，从而可与观测结果直接比较。

实验结果

研究问题

• RQ1 在致密星系核中，恒星碰撞在中心黑洞增长中起到何种作用？
• RQ2 潮汐撕裂在喂养中心黑洞及在原本安静的星系核中产生吸积耀斑的过程中扮演何种角色？
• RQ3 真实碰撞结果（包括非撕裂性并合）如何影响星系中心奇特恒星种群的形成？
• RQ4 基于SPH的碰撞物理的蒙特卡洛代码在多大程度上能再现百万粒子系统在宇宙时标内的动力学与演化？
• RQ5 引入不断增长的黑洞与真实恒星碰撞后，对星系核长期的结构与动力学演化产生何种影响？

主要发现

• 该代码可在普通个人电脑上仅用数个CPU天完成百万粒子系统在宇宙时标内的模拟，性能与N体模拟相当。
• 潮汐撕裂率主要由高度椭圆轨道的恒星主导，其速率符合相空间中由松弛驱动扩散的预期规律。
• 恒星碰撞显著促进向中心黑洞的质量转移，尤其在高密度环境中，且可能产生在安静星系核中可探测到的明亮耀斑。
• 非撕裂性碰撞可形成大质量、致密的恒星，可能解释银河系中心核星团中观测到的奇特恒星种群。
• 通过SPH模拟引入真实碰撞物理，相比解析近似，显著提升了质量转移与并合结果建模的准确性。
• 全面的测试模拟验证了代码的数值稳定性和正确性，证实其适用于研究星系核长期演化。