SCYON Abstract

Received on May 18 2001

A new Monte Carlo code for star cluster simulations: I.~Relaxation

Authors	Marc Freitag (^{1, 2}), Willy Benz (²)
Affiliation	(¹) Observatoire de Genève, CH-1290 Sauverny, Switzerland (²) Physikalisches Institut, Universität Bern, CH-3012 Bern, Switzerland
Accepted by	Astronomy & Astrophysics
Contact	marc.freitag@obs.unige.ch
URL	http://xxx.lanl.gov/abs/astro-ph/0102139v2
Links

Abstract

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 a numerical scheme pioneered by Hénon in the 70's. 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 first 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 x 10⁶ particles can be simulated on a standard personal computer and the CPU time required scales as N_p ln(N_p) with the particle number N_p. Furthermore, the number of simulated stars needs not be equal to N_p but can be arbitrarily larger. A companion paper will treat further physical elements, mostly relevant to galactic nuclei.