Impact of initial mass function on the chemical evolution of high-redshift galaxies

Abstract: Star formation and metal enrichment in galaxies are regulated by supernova (SN) explosions and metal yields from massive stars, which are sensitive to the high-mass end of the initial mass function (IMF). Recent JWST observations found evidence for an invariant relation between stellar mass, metallicity, and star formation rate up to $z\sim 8$ and its breakdown at higher redshifts. It is crucial to understand the underlying physics, especially the role played by the IMF. We explore the impact of IMF on the chemical evolution of high-redshift galaxies and the interplay between IMF and galactic outflow parameters. The ultimate goal is to constrain the high-mass end of the IMF by the cosmic star formation history and stellar mass-metallicity-star formation rate relation (MZSFR) inferred from observations at $z\sim 4-10$. Using the semi-analytical galaxy evolution code A-SLOTH, we follow galactic baryon cycles along merger trees built from a high-resolution cosmological simulation. Stellar feedback is modeled with up-to-date stellar evolution tracks covering the full metallicity range ($Z \sim 10^{-11} - 0.03$) and a broad stellar mass range ($m_\star\sim2 - 600\ \rm M_\odot$) including the metal yields from stellar winds and all types of SNe. Assuming a Kroupa-like shape of the IMF with a varying upper mass limit $m_{\max}$, we find $m_{\max} \gtrsim 200\ \rm M_\odot$ is required to reproduce the observed MZSFR. Observational data at $z\gtrsim 6$ favor a galactic outflow model where the outflow rate is proportional to the supernova energy injection rate divided by the halo binding energy. We conclude that very massive ($\gtrsim 200\ \rm M_\odot$) stars can play important roles in the star formation and chemical enrichment histories of high-$z$ galaxies and discuss the implications of our results for reionization and transient sources of both electromagnetic waves and gravitational waves.