2022MNRAS.510.5603K

ALMA observations of z > 6 galaxies have revealed abnormally high [O III]_88 µm_/[C II]_158 µm_ ratios and [C II]_158 µm_ deficits compared to local galaxies. The origin of this behaviour is unknown. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties, including ionization parameter, gas density, or photodissociation region (PDR) covering fraction. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX²⁰, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization which has been post-processed with CLOUDY to predict emission lines. We find that the observed z > 6 [O III]_88 µm_-SFR and [C II]_158 µm_-SFR relations can only be reproduced when the C/O abundance ratio is ∼8 x lower than Solar and the total metal production is ∼4 x higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [C II]_158 µm_-SFR and [C II]_158 µm_ luminosity functions are predicted to converge to observed values at z ∼ 4.5. While we demonstrate that ionization parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [O III]_88 µm_/[C II]_158 µm_, observed values at z > 6 can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [C II]_158 µm_ and [O III]_88 µm_ can be used to predict the values of ionization parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed z > 6 galaxies. Finally, we demonstrate that [O I]_63 µm_ can be used as a replacement for [C II]_158 µm_ in high-redshift galaxies where [C II]_158 µm_ is unobserved and argue that more observation time should be used to target [O I]_63 µm_ at z > 6. Future simulations will be needed to self-consistently address the numerous uncertainties surrounding a varying IMF at high redshift and the associated metal returns.