Synthetic nebular emission from massive galaxies - I: Origin of the cosmic evolution of optical emission-line ratios

Galaxies occupy different regions of the [OIII]λ5007/H β-versus-[N II]λ6584/H α emissionline ratio diagram in the distant and local Universe.We investigate the origin of this intriguing result by modelling self-consistently, for the first time, nebular emission from young stars, accreting black holes (BHs) and older, post-asymptotic giant branch (post-AGB) stellar populations in galaxy formation simulations in a full cosmological context. In post-processing, we couple new-generation nebular-emission models with high-resolution, cosmological zoom-in simulations of massive galaxies to explore which galaxy physical properties drive the redshift evolution of the optical-line ratios [OIII]λ5007/H β, [NII]λ6584/H α, [SII]λλ6717, 6731/H α and [O I]λ6300/H α. The line ratios of simulated galaxies agree well with observations of both star-forming and active local Sloan Digital Sky Survey galaxies. Towards higher redshifts, at fixed galaxy stellar mass, the average [O III]/H β is predicted to increase and [N II]/H α, [S II]/H α and [O I]/H α to decrease - widely consistent with observations. At fixed stellar mass, we identify star formation history, which controls nebular emission from young stars via the ionization parameter, as the primary driver of the cosmic evolution of [OIII]/H β and [N II]/H α. For [SII]/H α and [O I]/H α, this applies only to redshifts greater than z = 1.5, the evolution at lower redshift being driven in roughly equal parts by nebular emission from active galactic nuclei and post-AGB stellar populations. Instead, changes in the hardness of ionizing radiation, ionized-gas density, the prevalence of BH αccretion relative to star formation and the dust-to-metal mass ratio (whose impact on the gas-phase N/O ratio we model at fixed O/H) play at most a minor role in the cosmic evolution of simulated galaxy line ratios.