A new Monte Carlo code for star cluster simulations - I. Relaxation
A new Monte Carlo code for star cluster simulations - I. Relaxation
复制标题
用于星团模拟的新蒙特卡罗代码 - I. 松弛
DOI:
10.1051/0004-6361:20010706
复制
发表时间:
2001
影响因子:
6.5
通讯作者:
Bern University
中科院分区:
文献类型:
--
作者:
M. Freitag;Willy Benz Geneva Observatory;Bern University
We have developed a new simulation code aimed at studying the stellar dynamics of a galactic central star cluster surrounding a massive black hole. In order to include all the relevant physical ingredients (2-body relaxation, stellar mass spectrum, collisions, tidal disruption, ...) , we chose to revive an umerical scheme pioneered by H enon in the 70's (H enon 1971b,a; H enon 1973. It is basically a Monte Carlo resolution of the Fokker-Planck equation. It can cope with any stellar mass spectrum or velocity distribution. Being a particle-based method, it also allows one to take stellar collisions into account in a very realistic way. This rst paper covers the basic version of our code which treats the relaxation-driven evolution of a stellar cluster without a central BH. A technical description of the code is presented, as well as the results of test computations. Thanks to the use of a binary tree to store potential and rank information and of variable time steps, cluster models with up to 2 10 6 particles can be simulated on a standard personal computer and the CPU time required scales as Np ln(Np) with the particle number Np. Furthermore, the number of simulated stars needs not be equal to Np but can be arbitrarily larger. A companion paper will treat further physical elements, mostly relevant to galactic nuclei.