Introducing a new multi-particle collision method for the evolution of dense stellar systems II. Core collapse

Abstract

In a previous paper we introduced a new method for simulating collisional gravitational $N$-body systems with linear time scaling on $N$, based on the Multi-Particle Collision (MPC) approach. This allows us to simulate globular clusters with a realistic number of stellar particles in a matter of hours on a typical workstation. We evolve star clusters containing up to $10^6$ stars to core collapse and beyond. We quantify several aspects of core collapse over multiple realizations and different parameters, while always resolving the cluster core with a realistic number of particles. We run a large set of N-body simulations with our new code. The cluster mass function is a power-law with no stellar evolution, allowing us to clearly measure the effects of the mass spectrum. Leading up to core collapse, we find a power-law relation between the size of the core and the time left to core collapse. Our simulations thus confirm the theoretical self-similar contraction picture but with a dependence on the slope of the mass function. The time of core collapse has a non-monotonic dependence on the slope, which is well fit by a parabola. This holds also for the depth of core collapse and for the dynamical friction timescale of heavy particles. Cluster density profiles at core collapse show a broken power law structure, suggesting that central cusps are a genuine feature of collapsed cores. The core bounces back after collapse, and the inner density slope evolves to an asymptotic value. The presence of an intermediate-mass black hole inhibits core collapse. We confirm and expand on several predictions of star cluster evolution before, during, and after core collapse. Such predictions were based on theoretical calculations or small-size direct $N$-body simulations. Here we put them to the test on MPC simulations with a much larger number of particles, allowing us to resolve the collapsing core.

arXiv comment8 pages, 9 figs., 1 tab. Version accepted for publication in A&A

References in corpus (36)

• Dynamical evolution of star clusters in tidal fields
• Intermediate-Mass Black Holes
• Star Clusters Across Cosmic Time
• Internal Dynamics of Globular Clusters
• The DRAGON simulations: globular cluster evolution with a million stars
• NBODY6++GPU: Ready for the gravitational million-body problem
• Statistics of N-body Simulations I. Equal Masess Before Core Collapse
• Performance Analysis of Direct N-Body Algorithms on Special-Purpose Supercomputers
• MOCCA Code for Star Cluster Simulations - II. Comparison with N-body Simulations
• PeTar: a high-performance N-body code for modeling massive collisional stellar systems
• MOCCA code for star cluster simulations - I. Blue Stragglers, first results
• The Physics of Star Cluster Formation and Evolution
• Two body relaxation in CDM simulations
• Ratios of star cluster core and half-mass radii: a cautionary note on intermediate-mass black holes in star clusters
• On the mass-radius relation of hot stellar systems
• Post-Collapse evolution of globular clusters
• A comprehensive Nbody study of mass segregation in star clusters: Energy equipartition and escape
• Long-term evolution of isolated N-body sytems
• Tidal disruption, global mass function and structural parameters evolution in star clusters
• Statistics of N-Body Simulations. II. Equal Masses after Core Collapse
• Direct collisional simulation of 10 000 particles past core collapse
• A prescription and fast code for the long-term evolution of star clusters
• Monte Carlo Simulations of Globular Cluster Evolution. VI. The Influence of an Intermediate Mass Black Hole
• Two-Body Relaxation in Cosmological Simulations
• Towards an N-body Model for the Globular Cluster M4
• Different Fates of Young Star Clusters After Gas Expulsion
• A direct N-body model of core-collapse and core oscillations
• Parameters of core-collapse
• N-body models of globular clusters: metallicity, half-light radii and mass-to-light ratios
• Widespread presence of shallow cusps in the surface-brightness profile of globular clusters
• Multi-particle collision simulations of 2D one-component plasmas: anomalous transport and dimensional crossovers
• Surface Brightness Profiles for a sample of LMC, SMC and Fornax galaxy Globular Clusters
• Two Body Relaxation in Simulated Cosmological Haloes
• The hunt for self-similar core collapse
• Introducing a new multi-particle collision method for the evolution of dense stellar systems. Crash-test N-body simulations
• Dynamical Friction from field particles with a mass spectrum

Cited by in corpus (1)

• Dynamics of intermediate mass black holes in globular clusters. Wander radius and anisotropy profiles