We model dust formation in the core collapse supernova explosion SN 1987A by treating the gas-phase formation of dust grain nuclei as a chemical process. To compute the synthesis of fourteen species of grains we integrate a non-equilibrium network of nucleating and related chemical reactions and follow the growth of the nuclei into grains via accretion and coagulation. The effects of the radioactive cobalt, titanium, and sodium on the thermodynamics and chemistry of the ejecta are taken into account. The grain temperature, which we allow to differ from the gas temperature, affects the surface-tension-corrected evaporation rate. We also account for He$^+$, Ne$^+$, Ar$^+$, and O weathering. We combine our dust synthesis model with a crude prescription for anisotropic radioactive nickel dredge-up into the core ejecta, the so-called `nickel bubbles', to compute the total dust mass and molecular-species-specific grain size distribution. The total mass varies between $0.41\,M_\odot$ and $0.73\,M_\odot$, depending on the bubble shell density contrast. In the decreasing order of abundance, the grain species produced are: magnesia, silicon, forsterite, iron sulfide, carbon, silicon dioxide, alumina, and iron. The combined grain size distribution is a power law $dN/da\propto a^{-4.39}$. Early ejecta compaction by expanding radioactive nickel bubbles strongly enhances dust synthesis. This underscores the need for improved understanding of hydrodynamic transport and mixing over the entire pre-homologous expansion.

We present a simulation of the long-term evolution of a Population III supernova remnant in a cosmological minihalo. Employing passive Lagrangian tracer particles, we investigate how chemical stratification and anisotropy in the explosion can affect the abundances of the first low-mass, metal-enriched stars. We find that reverse shock heating can leave the inner mass shells at entropies too high to cool, leading to carbon-enhancement in the re-collapsing gas. This hydrodynamic selection effect could explain the observed incidence of carbon-enhanced metal-poor (CEMP) stars at low metallicity. We further explore how anisotropic ejecta distributions, recently seen in direct numerical simulations of core-collapse explosions, may translate to abundances in metal-poor stars. We find that some of the observed scatter in the Population II abundance ratios can be explained by an incomplete mixing of supernova ejecta, even in the case of only one contributing enrichment event. We demonstrate that the customary hypothesis of fully-mixed ejecta clearly fails if post-explosion hydrodynamics prefers the recycling of some nucleosynthetic products over others. Furthermore, to fully exploit the stellar-archaeological program of constraining the Pop III initial mass function from the observed Pop II abundances, considering these hydrodynamical transport effects is crucial. We discuss applications to the rich chemical structure of ultra-faint dwarf satellite galaxies, to be probed in unprecedented detail with upcoming spectroscopic surveys.

To constrain the properties of the first stars with the chemical abundance patterns observed in metal-poor stars, one must identify any non-trivial effects that the hydrodynamics of metal dispersal can imprint on the abundances. We use realistic cosmological hydrodynamic simulations to quantify the distribution of metals resulting from one Population III supernova and from a small number of such supernovae exploding in close succession. Overall, supernova ejecta are highly inhomogeneously dispersed throughout the simulations. When the supernova bubbles collapse, quasi-virialized metal-enriched clouds, fed by fallback from the bubbles and by streaming of metal-free gas from the cosmic web, grow in the centers of the dark matter halos. Partial turbulent homogenization on scales resolved in the simulation is observed only in the densest clouds where the vortical time scales are short enough to ensure true homogenization on subgrid scales. However, the abundances in the clouds differ from the gross yields of the supernovae. Continuing the simulations until the cloud have gone into gravitational collapse, we predict that the abundances in second-generation stars will be deficient in the innermost mass shells of the supernova (if only one has exploded) or in the ejecta of the latest supernovae (when multiple have exploded). This indicates that hydrodynamics gives rise to biases complicating the identification of nucleosynthetic sources in the chemical abundance spaces of the surviving stars.