Monte-Carlo Simulations of Star Clusters I. First Results

Abstract

A revision of Stodolkiewicz's Monte-Carlo code is used to simulate evolution of star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The first calculations for isolated, equal-mass N-body systems with three-body energy generation according to Spitzer's formulae show good agreement with direct N-body calculations for N=2000, 4096 and 10000 particles. The density, velocity, mass distributions, energy generation, number of binaries etc. follow the N-body results. Only the number of escapers is slightly too high compared to N-body results 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 < 2000. For simulations with N > 10000 gravothermal oscillations are clearly visible. The calculations of N=2000, 4096, 10000, 32000 and 100000 models take about 2, 6, 20, 130 and 2500 hours, respectively. The Monte-Carlo code is at least 10^5 times faster than the N-body one for N=32768 with special-purpose hardware (Makino 1996ab). Thus it becomes possible to run several different models to improve statistical quality of the data and run individual models with N as large as 100000. The Monte-Carlo scheme can be regarded as a method which lies in the middle between direct N-body and Fokker-Planck models and combines most advantages of both methods.

11 pages, 8 PS-figures, To appear in MNRAS

Cited by in corpus (34)

• Formation of Massive Black Holes in Dense Star Clusters. I. Mass Segregation and Core Collapse
• The Evolution of Globular Clusters in the Galaxy
• Monte-Carlo Simulations of Globular Cluster Evolution - I. Method and Test Calculations
• Monte Carlo Simulations of Globular Cluster Evolution. V. Binary Stellar Evolution
• Population synthesis of black hole binary mergers from star clusters
• A Parallel Monte Carlo Code for Simulating Collisional N-body Systems
• A novel look at energy equipartition in globular clusters
• Monte Carlo Simulations of Star Clusters - VII. The globular cluster 47 Tuc
• Preparing the next gravitational million-body simulations: Evolution of single and binary stars in Nbody6++GPU, MOCCA and McLuster
• Monte Carlo cluster simulations to determine the rate of compact star inspiraling to a central galactic black hole
• Compact Binaries in Star Clusters II - Escapers and Detection Rates
• A prescription and fast code for the long-term evolution of star clusters
• MOCCA SURVEY Database I: Binary Black Hole Mergers from Globular Clusters with Intermediate Mass Black Holes
• The impact of stellar evolution on rotating star clusters: the gravothermal-gravogyro catastrophe and the formation of a bar of black holes
• MOCCA-SURVEY database I. Accreting white dwarf binary systems in globular clusters -- I. cataclysmic variables -- present-day population
• Modest-2: A Summary
• Rates of capture of stars by supermassive black holes in non-spherical galactic nuclei
• Monte Carlo simulations of multiple populations in globular clusters: constraints on the cooling flow vs. accretion scenario using million bodies simulations
• Performance analysis of direct N-body algorithms for astrophysical simulations on distributed systems
• Dynamical modelling of globular clusters: challenges for the robust determination of IMBH candidates
• MOCCA: Dynamics and evolution of binary stars of multiple stellar populations in tidally filling and underfilling globular star clusters
• Is there a Size Difference between Red and Blue Globular Clusters?
• Three-Body Binary Formation in Clusters: Analytical Theory
• Understanding the central kinematics of globular clusters with simulated integrated-light IFU observations
• MOCCA-SURVEY Database II -- Properties of Intermediate Mass Black Holes escaping from star clusters
• A Monte Carlo study of early gas expulsion and evolution of star clusters: new simulations with the MOCCA code in the AMUSE framework
• Properties of quasi-relaxed stellar systems in an external tidal field
• The initial conditions of observed star clusters - I. Method description and validation
• Monte Carlo simulations of multiple populations in globular clusters: constraints on the initial size of the second generation from binary stars
• Forming short period sub-stellar companions in 47 Tucanae -- II. Analytic expressions for the orbital evolution of planets in dense environments
• Gravitational Dynamics of Large Stellar Systems
• Studying Binary Systems in Omega Centauri with MUSE. I. Detection of Spectroscopic Binaries
• Forming short period sub-stellar companions in 47 Tuc: I. Dynamical model and brown dwarf tidal capture rates
• KRIOS: A new basis-expansion $N$-body code for collisional stellar dynamics