A panchromatic view of star cluster formation in a simulated dwarf galaxy starburst

We present a photometric analysis of star and star cluster (SC) formation in a high-resolution simulation of a dwarf galaxy starburst that allows the formation of individual stars to be followed. Previous work demonstrated that the properties of the SCs formed in the simulation are in good agreement with observations. In this paper, we create mock spectral energy distributions and broad-band photometric images using the radiative transfer code skirt 9. We test several observational star formation rate (SFR) tracers and find that 24 μm, total infrared and Hα trace the underlying SFR during the (post)starburst phase, while UV tracers yield a more accurate picture of star formation during quiescent phases prior to and after the merger. We then place the simulated galaxy at distances of 10 and 50 Mpc and use aperture photometry at Hubble Space Telescope resolution to analyse the simulated SC population. During the starburst phase, a hierarchically forming set of SCs leads inaccurate source separation because of crowding. This results in estimated SC mass function slopes that are up to ∼0.3 shallower than the true slope of ∼−1.9 to −2 found for the bound clusters identified from the particle data in the simulation. The masses of the largest clusters are overestimated by a factor of up to 2.9 due to unresolved clusters within the apertures. The aperture-based analysis also produces a relation between cluster formation efficiency and SFR surface density that is slightly flatter than that recovered from bound clusters. The differences are strongest in quiescent SF environments

The challenge of simulating the star cluster population of dwarf galaxies with resolved interstellar medium

We present results on the star cluster properties from a series of high resolution smoothed particles hydrodynamics (SPH) simulations of isolated dwarf galaxies as part of the griffin project. The simulations at sub-parsec spatial resolution and a minimum particle mass of 4 M⊙ incorporate non-equilibrium heating, cooling, and chemistry processes, and realize individual massive stars. The simulations follow feedback channels of massive stars that include the interstellar-radiation field variable in space and time, the radiation input by photo-ionization and supernova explosions. Varying the star formation efficiency per free-fall time in the range ϵff = 0.2–50 per cent neither changes the star formation rates nor the outflow rates. While the environmental densities at star formation change significantly with ϵff, the ambient densities of supernovae are independent of ϵff indicating a decoupling of the two processes. At low ϵff, gas is allowed to collapse more before star formation, resulting in more massive, and increasingly more bound star clusters are formed, which are typically not destroyed. With increasing ϵff, there is a trend for shallower cluster mass functions and the cluster formation efficiency Γ for young bound clusters decreases from 50 per cent to ∼1 per cent showing evidence for cluster disruption. However, none of our simulations form low mass (3 M⊙) clusters with structural properties in perfect agreement with observations. Traditional star formation models used in galaxy formation simulations based on local free-fall times might therefore be unable to capture star cluster properties without significant fine tuning