Search for: All records

Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics codeFLASHto theN-body codeph4and the stellar evolution codeSeBausing the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisionalN-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds, and supernovae inFLASH. We here describe the algorithms used for each of these processes. We denote this integrated package Torch. We then use this novel numerical method to simulate the formation and early evolution of several examples of open clusters of ∼1000 stars formed from clouds with a mass range of 10³M_⊙to 10⁵M_⊙. Analyzing the effects of stellar feedback on the gas and stars of the natal clusters, we find that in these examples, the stellar clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H iiregion shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations.

 

Context. Supersonic disordered flows accompany the formation and evolution of molecular clouds (MCs). It has been argued that this is turbulence that can support against gravitational collapse and form hierarchical sub-structures. Aims. We examine the time evolution of simulated MCs to investigate: What physical process dominates the driving of turbulent flows? How can these flows be characterised? Are they consistent with uniform turbulence or gravitational collapse? Do the simulated flows agree with observations? Methods. We analysed three MCs that have formed self-consistently within kiloparsec-scale numerical simulations of the interstellar medium (ISM). The simulated ISM evolves under the influence of physical processes including self-gravity, stratification, magnetic fields, supernova-driven turbulence, and radiative heating and cooling. We characterise the flows using velocity structure functions (VSFs) with and without density weighting or a density cutoff, and computed in one or three dimensions. However, we do not include optical depth effects that can hide motions in the densest gas, limiting comparison of our results with observations. Results. In regions with sufficient resolution, the density-weighted VSFs initially appear to follow the expectations for uniform turbulence, with a first-order power-law exponent consistent with Larson’s size-velocity relationship. Supernova blast wave impacts on MCs produce short-lived coherent motions at large scales, increasing the scaling exponents for a crossing time. Gravitational contraction drives small-scale motions, producing scaling coefficients that drop or even turn negative as small scales become dominant. Removing the density weighting eliminates this effect as it emphasises the diffuse ISM. Conclusions. We conclude that two different effects coincidentally reproduce Larson’s size velocity relationship. Initially, uniform turbulence dominates, so the energy cascade produces VSFs that are consistent with Larson’s relationship. Later, contraction dominates and the density-weighted VSFs become much shallower or even inverted, but the relationship of the global average velocity dispersion of the MCs to their radius follows Larson’s relationship, reflecting virial equilibrium or free-fall collapse. The injection of energy by shocks is visible in the VSFs, but decays within a crossing time.