  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.
