Abstract
The chemical abundances of stars in galaxies are a fossil record of the star formation and stellar evolution processes that regulate galaxy formation, including the stellar initial mass function, the fraction and timing of type Ia supernovae (SNeIa), and nucleosynthesis inside massive stars. In this paper, we systematically explore uncertainties associated with modeling chemical enrichment in dwarf galaxies. We repeatedly simulate a single EDGE-INFERNO dwarf ( M ★ ≈ 10 5 M ⊙ ), varying the chemical yields of massive stars, the timing and yields of SNeIa, and the intrinsic stochasticity that arises from sampling individual stars and galaxy formation chaoticity. All simulations are high-resolution (3.6 pc), cosmological zoom-in hydrodynamical simulations that track the stellar evolution of all individual stars with masses of > 0.5 M ⊙ . We find that SNeIa make significant contributions to the iron content of low-mass, reionization-limited galaxies, with possible variations in mean abundance ratios and [Fe/H] related to minor changes in their evolutionary timescales. In contrast, different massive star yields, accounting (or not) for stellar rotation, result in mean abundance variations comparable to those arising from stochasticity, with the possible exception of extremely rapidly rotating stars. Nonetheless, massive stars significantly affect the shape of abundance trends with [Fe/H], for example, through the existence (or not) of a bimodality in the [X/Fe]–[Fe/H] planes, particularly in [Al/Fe]. Finally, we find that the variance arising from random sampling severely limits the interpretation of single galaxies. Our analysis showcases the power of star-by-star cosmological models to unpick how both systematic uncertainties (e.g., assumptions in low-metallicity chemical enrichment) and statistical uncertainties (e.g., averaging over enough galaxies and stars within a galaxy) affect the interpretation of chemical observables in ultra-faint dwarf galaxies.