DAWN. I. Simulating the formation and early evolution of stellar clusters with Phantom N-Body

Abstract

Context. Simulating stellar dynamics in a molecular cloud environment is numerically challenging due to the strong coupling between young stars and their surrounding gas, and the large range of length and time scales. Aims. This paper is the first of a suite aimed at investigating the complex early stellar dynamics in star-forming regions. We present a new simulation framework which is the key to generating a larger set of simulations, enabling statistical analysis. Methods. Methods originating from the stellar dynamics community, including regularisation and slowdown methods (SDAR), have been added to the hydrodynamical code Phantom to produce simulations of embedded cluster early dynamics. This is completed by a novel prescription of star formation to initialise stars with a low numerical cost, but in a way that is consistent with the gas distribution. Finally, a prescription for H ii region expansion has been added to model the gas removal. Results. We have run testcase simulations following the dynamical evolution of stellar clusters from the cloud collapse to a few Myr. Our new numerical methods fulfil their function by speeding up the calculation. The N-body dynamics with our novel implementation never appear as a bottleneck. Our first simulations show that massive stars largely impact the star formation process and shape the dynamics of the resulting cluster. Depending on the position of these massive stars and the strength of their feedback, they can prematurely dismantle part of the cloud or trigger a second event of cloud collapse, preferentially forming low-mass stars. This stochastic behaviour confirms the need for statistical studies. Conclusions. Our new Phantom N-Body framework enables efficient simulation of the formation and evolution of star clusters. It enables the statistical analysis needed to build models of the dynamical evolution of embedded star clusters.

References in corpus (16)

• SPLASH: An interactive visualisation tool for Smoothed Particle Hydrodynamics simulations
• A Spitzer Survey of Young Stellar Clusters within One Kiloparsec of the Sun: Cluster Core Extraction and Basic Structural Analysis
• Smoothed Particle Hydrodynamics and Magnetohydrodynamics
• The Importance of Radiative Feedback for the Stellar Initial Mass Function
• Evidence for the Strong Effect of Gas Removal on the Internal Dynamics of Young Stellar Clusters
• PeTar: a high-performance N-body code for modeling massive collisional stellar systems
• Towards a more realistic sink particle algorithm for the RAMSES code
• N-body simulations of gravitational dynamics
• The origin of the Milky Way globular clusters
• Smoothed Particle Hydrodynamics simulations of expanding HII regions. I. Numerical methods and tests
• The W43-MM1 mini-starburst ridge, a test for star formation efficiency models
• A slow-down time-transformed symplectic integrator for solving the few-body problem
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• Stellar mass spectrum within massive collapsing clumps III. Effects of temperature and magnetic field
• Does turbulence determine the initial mass function?
• Gas and star kinematics in cloud-cloud collisions