The nature of diffuse ionized gas in star-forming galaxies

Abstract

ABSTRACT We present an analysis of the diffuse ionized gas (DIG) in a high-resolution simulation of an isolated Milky Way-like galaxy, incorporating on-the-fly radiative transfer and non-equilibrium thermochemistry. We utilize the Monte-Carlo radiative transfer code colt to self-consistently obtain ionization states and line emission in post-processing. We find a clear bimodal distribution in the electron densities of ionized gas ($n_{\rm e}$), allowing us to define a threshold of $n_{\rm e}=10\, \mathrm{cm}^{-3}$ to differentiate DIG from ${\rm H\, {\small II}}$ regions. The DIG is primarily ionized by stars aged 5 – 25 Myr, which become exposed directly to low-density gas after ${\rm H\, {\small II}}$ regions have been cleared. Leakage from recently formed stars ($\lt 5$ Myr) is only moderately important for DIG ionization. We forward model local observations and validate our simulated DIG against observed line ratios in [${\rm S\, {\small II}}$]/H$\alpha$, [${\rm N\, {\small II}}$]/H$\alpha$, [${\rm O\, {\small I}}$]/H$\alpha$, and [${\rm O\, {\small III}}$]/H$\beta$ against $\Sigma _{\rm H\alpha }$. The mock observations not only reproduce observed correlations, but also demonstrate that such trends are related to an increasing temperature and hardening ionizing radiation field with decreasing $n_{\rm e}$. The hardening of radiation within the DIG is caused by the gradual transition of the dominant ionizing source with decreasing $n_{\rm e}$ from 0 to 25 Myr stars, which have progressively harder intrinsic ionizing spectra primarily due to the extended Wolf–Rayet phase caused by binary interactions. Consequently, the DIG line ratio trends can be attributed to ongoing star formation, rather than secondary ionization sources, and therefore present a potent test for stellar feedback and stellar population models.