SPICE: the connection between cosmic reionisation and stellar feedback in the first galaxies

We present SPICE, a new suite of RHD cosmological simulations targeting the epoch of reionisation. The goal of these simulations is to systematically probe a variety of stellar feedback models, including "bursty" and "smooth" forms of supernova energy injection, as well as poorly-explored scenarios such as hypernova explosions and radiation pressure. Subtle differences in the behaviour of supernova feedback drive profound differences in reionisation histories, with burstier forms of feedback causing earlier reionisation. We also find that some global galaxy properties, such as the dust-attenuated luminosity functions and star formation main sequence, remain degenerate between models. Stellar feedback and its strength determine the morphological mix of galaxies emerging by z = 5 and that the reionisation history is inextricably connected to intrinsic properties such as galaxy kinematics and morphology. While star-forming, massive disks are prevalent if supernova feedback is "smooth", "bursty" feedback preferentially generates dispersion-dominated systems. Different modes of feedback produce different strengths of outflows, altering the ISM/CGM in different ways, and in turn strongly affecting the escape of LyC photons. We establish a correlation between galaxy morphology and LyC escape fraction, revealing that dispersion-dominated systems have escape fractions 10-50 times higher than their rotation-dominated counterparts at all redshifts. Dispersion-dominated systems should thus preferentially generate large HII regions as compared to their rotation-dominated counterparts. Since dispersion-dominated systems are more prevalent if stellar feedback is more explosive, reionisation occurs earlier in our simulation with burstier feedback. Statistical samples of post-reionisation galaxy morphologies probed with JWST, ALMA and MUSE can constrain stellar feedback and models of cosmic reionisation.