Understanding Interstellar Metals during Reionization with Radiative SPH Simulation: Metallicity and Emission Lines from the ISM at $10 \geq z \geq 5$

Abstract: We compare the \texttt{Technicolor Dawn} cosmological simulations with recent observations of galactic nebular line emission during the Epoch of Reionization, providing stringent tests of the predicted ionization and metal enrichment levels. We validate the simulated population with the UVLF and $M_{\mathrm{UV}}-M_*$ relation and see that the simulated results are consistent with observations at lower masses. We extract local gas volumetric grids of density and mass-weighted metallicity, then we use \texttt{Cloudy} to produce synthetic emission spectra of \species{H}{ii} regions. The mass-metallicity relation does not evolve, which is also consistent with observations. The predicted oxygen abundance exceeds observational inferences by about 0.5 dex, suggesting either overly efficient enrichment or weak feedback. However, applying the O32 diagnostic directly to our synthetic spectra shows an offset of 1 dex from the correct outputted gas-phase metallicity. This suggests that O32 is biased high at a level that is more than sufficient to account for the simulation-observation offset. The simulated galaxies' line diagnostics show mostly weaker [\species{O}{iii}] lines and lower diagnostic values of O3 and Ne3O2 compared to observations. This suggests higher ionization parameters within the simulated galactic population in general.