Applying a star formation model calibrated on high-resolution interstellar medium simulations to cosmological simulations of galaxy formation

Abstract

Modern high-resolution simulations of the interstellar medium (ISM) have shown that key factors in governing star formation are the competing influences of radiative dissipation, pressure support driven by stellar feedback, and the relentless pull of gravity. Cosmological simulations of galaxy formation, such as IllustrisTNG or ASTRID, are however not able to resolve this physics in detail and therefore need to rely on approximate treatments. These have often taken the form of empirical subgrid models of the ISM expressed in terms of an effective equation of state (EOS) that relates the mean ISM pressure to the mean gas density. Here we seek to improve these heuristic models by directly fitting their key ingredients to results of the high-resolution TIGRESS simulations, which have shown that the dynamical equilibrium of the ISM can be understood in terms of a pressure-regulated, feedback modulated (PRFM) model for star formation. Here we explore a simple subgrid model that draws on the PRFM concept but uses only local quantities. It accurately reproduces PRFM for pure gas disks, while it predicts slightly less star formation than PRFM in the presence of an additional thin stellar disk. We compare the properties of this model with the older Springel and Hernquist and TNG prescriptions, and apply all three to isolated simulations of disk galaxies as well as to a set of high-resolution zoom-in simulations carried out with a novel 'multi-zoom' technique that we introduce in this study. The softer EOS implied by TIGRESS produces substantially thinner disk galaxies, which has important ramifications for disk stability and galaxy morphology. The total stellar mass of galaxies is however hardly modified at low redshift, reflecting the dominating influence of large-scale gaseous inflows and outflows to galaxies, which are not sensitive to the EOS itself