The first supernovae enrich the previously pristine gas with metals, out of which the next generation of stars form. Based on hydrodynamical simulations, we develop a new stochastic model to predict the metallicity of star-forming gas in the first galaxies. On average, in internally enriched galaxies, the metals are well mixed with the pristine gas. However, in externally enriched galaxies, the metals can not easily penetrate into the dense gas, which yields a significant metallicity difference between the star-forming and average gas inside a halo. To study the consequences of this effect, we apply a semi-analytical model to Milky Way-like dark matter merger trees and follow stellar fossils from high redshift until the present day with a novel realistic metal mixing recipe. We calibrate the model to reproduce the metallicity distribution function (MDF) at low metallicities and find that a primordial IMF with slope of $\mathrm{d}N/\mathrm{d}M \propto M^{-0.5}$ from $2 Msun$ to $180 Msun$ best reproduces the MDF. Our improved model for inhomogeneous mixing can have a large impact for individual minihalos, but does not significantly influence the modelled MDF at [Fe/H]$\gtrsim -4$ or the best-fitting Pop~III IMF.

We study the enrichment and mixing of r-process elements in ultrafaint dwarf galaxies (UFDs). We assume that r-process elements are produced by neutron-star mergers (NSMs), and examine multiple models with different natal kick velocities and explosion energies. To this end, we perform cosmological simulations of galaxy formation to follow mixing of the dispersed r-process elements driven by star formation and the associated stellar feedback in progenitors of UFDs. We show that the observed europium abundance in Reticulum II is reproduced by our inner explosion model where a NSM is triggered at the centre of the galaxy, whereas the relatively low abundance in Tucana III is reproduced if a NSM occurs near the virial radius of the progenitor galaxy. The latter case is realised only if the neutron-star binary has a large natal kick velocity and travels over a long distance of a kilo-parsec before merger. In both the inner and outer explosion cases, it is necessary for the progenitor galaxy to sustain prolonged star formation over a few hundred million years after the NSM, so that the dispersed r-process elements are well mixed within the inter-stellar medium. Short-duration star formation results in inefficient mixing, and then a large variation is imprinted in the stellar europium abundances, which is inconsistent with the observations of Reticulum II and Tucana III.