Introducing a new multi-particle collision method for the evolution of dense stellar systems. Crash-test N-body simulations

Abstract

Stellar systems are broadly divided into collisional and non-collisional. The latter are large-N systems with long relaxation timescales and can be simulated disregarding two-body interactions, while either computationally expensive direct N-body simulations or approximate schemes are required to properly model the former. Large globular clusters and nuclear star clusters, with relaxation timescales of the order of a Hubble time, are small enough to display some collisional behaviour and big enough to be impossible to simulate with direct $N$-body codes and current hardware. We introduce a new method to simulate collisional stellar systems, and validate it by comparison with direct $N$-body codes on small-$N$ simulations. The Multi-Particle collision for Dense stellar systems Code (MPCDSS) is a new code for evolving stellar systems with the Multi-Particle Collision method. Such method amounts to a stochastic collision rule that allows to conserve exactly the energy and momentum over a cluster of particles experiencing the collision. The code complexity scales with $N\log N$ in the number of particles. Unlike Monte-Carlo codes, MPCDSS can easily model asymmetric, non-homogeneous, unrelaxed and rotating systems, while allowing us to follow the orbits of individual stars. We evolve small ($N=3.2\times 10^4$) star clusters with MPCDSS and with the direct-summation code NBODY6, finding a similar evolution of key indicators. We then simulate different initial conditions in the $10^4-10^6$ star range. MPCDSS bridges the gap between small, collisional systems that can be simulated with direct $N$-body codes and large noncollisional systems. MPCDSS in principle allows us to simulate globular clusters such as Omega Cen and M54 and even the nuclear star cluster, beyond the limits of current direct N-body codes in terms of the number of particles.

12 pages, 13 figures. Matching the version accepted for publication in A&A

References in corpus (29)

• Evidence for Two Early Accretion Events That Built the Milky Way Stellar Halo
• MOCCA-SURVEY Database I: Coalescing Binary Black Holes Originating From Globular Clusters
• Very massive stars, pair-instability supernovae and intermediate-mass black holes with the SEVN code
• Merging black hole binaries: the effects of progenitor's metallicity, mass-loss rate and Eddington factor
• Merging black holes in young star clusters
• Dynamics of stellar black holes in young star clusters with different metallicities - II. Black hole-black hole binaries
• Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids
• Sausage & Mash: The dual origin of the Galactic thick disc and halo from the gas-rich Gaia-Enceladus-Sausage merger
• The Globular Cluster Migratory Origin of Nuclear Star Clusters
• Newton vs the machine: solving the chaotic three-body problem using deep neural networks
• Dynamical Formation of Low-Mass Merging Black Hole Binaries like GW151226
• Dynamics of black hole - neutron star binaries in young star clusters
• Self-consistent simulations of Nuclear Cluster formation through Globular Cluster orbital decay and merging
• The possible role of stellar mergers for the formation of multiple stellar populations in globular clusters
• Effects of Intermediate Mass Black Holes on Nuclear Star Clusters
• A new Monte Carlo method for dynamical evolution of non-spherical stellar systems
• The Main Sequence of Star Clusters
• Formation of super-massive black holes in galactic nuclei I: delivering seed intermediate-mass black holes in massive stellar clusters
• A deep view into the nucleus of the Sagittarius Dwarf Spheroidal Galaxy with MUSE. II. Kinematic characterization of the stellar populations
• A dynamical gravitational wave source in a dense cluster
• Probing the role of dynamical friction in shaping the BSS radial distribution. I - Semi-analytical models and preliminary N-body simulations
• Treatment of realistic tidal field in Monte Carlo simulations of star clusters
• Dynamical origin of non-thermal states in galactic filaments
• Discreteness effects, $N-$body chaos and the onset of radial-orbit instability
• Stellar-Encounter Driven Red-Giant Star Mass-Loss in Globular Clusters
• X-ray sources in Galactic globular clusters and old open clusters
• Brownian motion of supermassive black holes in galaxy cores
• Radial Dependence of the Proto-Globular Cluster Contribution to the Milky Way Formation
• Taking apart the dynamical clock. Fat-tailed dynamical kicks shape the blue-straggler star bimodality

Cited by in corpus (3)

• Symplectic coarse graining approach to the dynamics of spherical self-gravitating systems
• Dynamics of intermediate mass black holes in globular clusters. Wander radius and anisotropy profiles
• Introducing a new multi-particle collision method for the evolution of dense stellar systems II. Core collapse