M. Giersz, MONTE-CARLO SIMULATIONS OF STAR-CLUSTERS - I - FIRST RESULTS, Monthly Notices of the Royal Astronomical Society, 298(4), 1998, pp. 1239-1248
A revision of Stodolkiewicz's Monte Carlo code is used to simulate evo
lution of star clusters. The new method treats each superstar as a sin
gle star and follows the evolution and motion of all individual stella
r objects. The first calculations for isolated, equal-mass N-body syst
ems with three-body energy generation according to Spitzer's formulae
show good agreement with direct N-body calculations for N = 2000, 4096
and 10 000 particles. The density, velocity, mass distributions, ener
gy generation, number of binaries, etc., follow the N-body results. On
ly the number of escapers is slightly too high compared with N-body re
sults, and there is no level-off anisotropy for advanced post-collapse
evolution of Monte Carlo models as is seen in N-body simulations for
N less than or equal to 2000. For simulations with N > 10 000 gravothe
rmal oscillations are clearly visible. The calculations of N = 2000, 4
096, 10 000, 32 000 and 100 000 models take about 2, 6, 20, 130 and 25
00 h, respectively. The Monte Carlo code is at least 10(5) times faste
r than the N-body one for N = 32768 with special-purpose hardware. Thu
s it becomes possible to run several different models to improve stati
stical quality of the data and run individual models with N as large a
s 100 000. The Monte Carlo scheme can be regarded as a method which li
es in the middle between direct N-body and Fokker-Planck models and co
mbines most advantages of both methods.