Simulating star clusters across cosmic time - I. Initial mass function, star formation rates, and efficiencies

Abstract

We present radiation-magneto-hydrodynamic simulations of star formation in self-gravitating, turbulent molecular clouds,modelling the formation of individual massive stars, including their UV radiation feedback. The set of simulations have cloud masses between mgas = 103Mθ and 3 × 105 Mθ and gas densities typical of clouds in the local Universe (ngas ∼ 1.8 × 102 cm-3) and 10× and 100× denser, expected to exist in high-redshift galaxies. The main results are as follows. (i) The observed Salpeter power-law slope and normalization of the stellar initial mass function at the high-mass end can be reproduced if we assume that each star-forming gas clump (sink particle) fragments into stars producing on average a maximum stellar mass about 40 per cent of the mass of the sink particle, while the remaining 60 per cent is distributed into smaller mass stars. Assuming that the sinks fragment according to a power-law mass function flatter than Salpeter, with log-slope 0.8, satisfy this empirical prescription. (ii) The star formation law that best describes our set of simulation is dρ/dt ∝ ρ1.5 gas if ngas < 0.1 Zθ, based on our limited sample, f Is reduced by a factor of ∼5. (iv) The most compact and massive clouds appear to form globular cluster progenitors, in the sense that star clusters remain gravitationally bound after the gas has been expelled.