---
title: 'Diffuse Ionized Gas (DIG): Astrophysical Insights'
url: https://www.emergentmind.com/topics/diffuse-ionized-gas-dig
type: topic
---

# Diffuse Ionized Gas (DIG): Astrophysical Insights

Diffuse ionized gas (DIG) is the faint, extended warm ionized component of the interstellar medium outside classical H II regions, often identified with the warm ionized medium in the Milky Way. In galaxies it occupies disks, interarm regions, bulges, extraplanar layers, and stripped tails, and it contributes a substantial fraction of the nebular emission budget. Relative to classical H II regions, DIG is characterized by lower density, lower surface brightness, and line ratios that are often enhanced in low-ionization forbidden lines such as \([\mathrm{N\,II}]\), \([\mathrm{S\,II}]\), \([\mathrm{O\,I}]\), and sometimes \([\mathrm{O\,II}]\) relative to Balmer recombination lines. Because integrated galaxy spectra mix DIG and H II-region emission, DIG is central to the interpretation of diagnostic diagrams, star-formation rates, and gas-phase metallicities [2005.06054][1711.07844].

## 1. Physical definition and astrophysical environments

DIG is observed in multiple galactic environments: early-type galaxies, bulges of late-type galaxies, interarm regions of disks, and gas layers above and below galactic planes. In the Milky Way context it is commonly referred to as the warm ionized medium or diffuse ionized medium [2005.06054]. In star-forming galaxies it is frequently described as warm gas with \(T_e \sim 10^4\) K, lower density than classical H II regions, lower ionization parameter, lower emission-line equivalent widths, and lower surface brightness, while also exhibiting a harder ionizing spectrum in many circumstances [2109.02684].

The physical contrast with classical H II regions is not solely geometric. DIG tends to show enhanced collisionally excited lines relative to recombination lines, especially \([\mathrm{N\,II}]/\mathrm{H}\alpha\), \([\mathrm{S\,II}]/\mathrm{H}\alpha\), and often \([\mathrm{O\,I}]/\mathrm{H}\alpha\), while \([\mathrm{O\,III}]/\mathrm{H}\beta\) can behave non-universally depending on environment and ionizing source [2005.06054][1612.02000]. In the Milky Way disk, the warm ionized medium is described as having density \(\sim 0.1\) and temperatures \(6000\)–\(10000\,\mathrm{K}\) [1709.09232]. In ram-pressure stripped systems, DIG is described as warm (\(\gtrsim 10^4\) K), low-density (\(\sim 10^{-1}\,\mathrm{cm^{-3}}\)) gas with vertical scale heights of \(1\)–\(2\) kpc [2011.08869].

Survey-based estimates show that DIG is not a minor residual. In spiral galaxies it has been cited as contributing roughly 30–60% of the total H\(\alpha\) emission [2005.06054]. In the CALIFA analysis of 391 galaxies, DIG defined as hDIG+mDIG contributes on average 56% of the total H\(\alpha\) flux, with strong morphology dependence: approximately \((100,0,0)\)% in E/S0, \((9,60,31)\)% in Sa–Sb, and \((0,13,87)\)% in later types for \((f_{\rm hDIG}, f_{\rm mDIG}, f_{\rm SFc})\) [1711.07844]. Other methodologies yield similarly broad but survey-dependent ranges, including 20%–90% of integrated disk flux in GASP galaxies and 40%–70% of total H\(\alpha\) in the BETIS showcase sample [2011.08869][2311.14254].

A simulation-based definition has also been proposed. In a high-resolution isolated Milky Way-like simulation, the ionized gas shows a bimodal electron-density distribution with peaks at \(n_{\rm e}\sim 0.1-1~\mathrm{cm^{-3}}\) and \(n_{\rm e}\sim 10^2-10^3~\mathrm{cm^{-3}}\), motivating a threshold \(n_{\rm e}=10~\mathrm{cm^{-3}}\) to separate DIG from H II-region gas [2403.03243]. This suggests that DIG may be physically distinct from classical H II regions rather than merely the low-surface-brightness tail of a continuous distribution.

## 2. Ionization sources and internal diversity

DIG is not a single ionization phenomenon. Proposed ionizing or heating mechanisms include leakage of ionizing photons from star-forming regions, hot low-mass evolved stars (HOLMES), shocks, turbulence dissipation, cosmic rays, dust-scattered light, old supernova remnants, magnetic reconnection, and turbulent mixing layers [2005.06054][2011.08869]. A useful tripartite distinction identifies “leaking-DIG,” ionized by photons escaping from H II regions; “HOLMES-DIG,” ionized by old post-AGB-like stars and commonly associated with LIER-like line ratios; and shocked DIG [2109.02684].

In actively star-forming galaxies, several studies argue that leaked radiation from young stars is energetically dominant. PHANGS-MUSE data support a two-component picture in which an energetically dominant DIG component is powered by Lyman-continuum photons leaking from H II regions, while a harder but energetically weaker component from hot low-mass evolved stars is required to explain enhanced low-ionization ratios in central regions, flat or decreasing \([\mathrm{O\,III}]/\mathrm{H}\beta\) with \(\Sigma_{\mathrm{H}\alpha}\), and offsets into LI(N)ER-like regions of BPT space. In that analysis, HOLMES contribute about 2% of the galaxy-integrated H\(\alpha\) emission but are fundamental contributors to \([\mathrm{O\,III}]\) emission [2111.14876].

A different but partly convergent picture emerges from galaxy-scale simulation. In the 2024 isolated-galaxy calculation, DIG is primarily ionized by stars aged \(5\)–\(25\) Myr that become directly exposed to low-density gas after H II regions have been cleared; leakage from recently formed stars younger than 5 Myr is only moderately important. The model attributes DIG line-ratio trends to increasing temperature and a hardening radiation field at lower \(n_{\rm e}\), with the hardening driven by the age shift in the dominant ionizing population and by the harder intrinsic spectra of \(5\)–\(25\) Myr stars in BPASS-like binary-evolution models [2403.03243]. This suggests that recent star formation alone can explain much of the DIG phenomenology in normal star-forming disks, though it does not invalidate the observational evidence for HOLMES in bulges and low-EW environments.

Environmental dependence is substantial. In stripped tails, DIG with strong \([\mathrm{O\,I}]/\mathrm{H}\alpha\) excess, high LIER/LINER incidence, and projected distances up to 10 kpc from star-forming regions is interpreted as at least partly ionized by processes other than star formation, probably mixing, shocks, and accretion of intracluster and interstellar gas [2108.12433]. In two local analogues of high-redshift star-forming galaxies, DIG is interpreted as most likely dominated by photon leakage from H II regions with additional contributions from feedback-driven shocks [2601.19725]. The aggregate picture is therefore plural: leaked young-star radiation, intermediate-age exposed populations, HOLMES, and shocks all appear in the literature, but their relative roles vary strongly with galaxy type and location within galaxies.

## 3. Observational signatures and classification schemes

The most widely used practical DIG classifier in resolved optical spectroscopy is the H\(\alpha\) equivalent width, \(W_{\mathrm{H}\alpha}\), defined as
\[
W_{\mathrm{H}\alpha} \equiv \frac{F_{\mathrm{H}\alpha}}{C_{\lambda,\,\mathrm{cont}}}.
\]
Because it measures nebular emission relative to the local stellar continuum, it links line emission to the underlying stellar population and behaves as an intensive quantity under line-of-sight superposition [1711.07844].

In CALIFA-based work, three nebular regimes are defined by \(W_{\mathrm{H}\alpha}\) [1711.07844][2005.06054]:

| Regime | Criterion | Interpretation |
|---|---:|---|
| hDIG | \(W_{\mathrm{H}\alpha} < 3\,\AA\) | HOLMES-dominated diffuse gas |
| mDIG | \(3 < W_{\mathrm{H}\alpha} < 14\,\AA\) | Mixed diffuse regime |
| SFc | \(W_{\mathrm{H}\alpha} > 14\,\AA\) | Star-forming complexes |

The \(3\,\AA\) threshold is physically linked to the “retired galaxy” interpretation and to the low-\(W_{\mathrm{H}\alpha}\) peak near \(1\,\AA\), whereas the \(14\,\AA\) threshold is empirical and marks the mode of the high-\(W_{\mathrm{H}\alpha}\) population [1711.07844]. Low-\(W_{\mathrm{H}\alpha}\) zones have \(\xi \equiv L_{\mathrm{H}\alpha}^{\mathrm{obs}}/L_{\mathrm{H}\alpha}^{\mathrm{pred}(t>10^8\,\mathrm{yr})}\) of order unity, supporting the idea that HOLMES can energetically power them [1711.07844].

This classifier maps cleanly onto excitation diagrams in many resolved data sets: SFc zones occupy the star-forming wing of the BPT plane, hDIG zones populate the LINER/LIER-like tip of the right wing, and mDIG forms the bridge between them [1711.07844]. The same framework underlies later MaNGA DIG-correction work, which is formally restricted to spectra with observed \(W_{\mathrm{H}\alpha}\gtrsim 10\,\AA\) because more weakly emitting spectra are already DIG-dominated [1907.08635].

Surface-brightness criteria remain common but are contested. A threshold such as \(\Sigma_{\mathrm{H}\alpha}>10^{39}\,\mathrm{erg\,s^{-1}\,kpc^{-2}}\) has been used to select H II-dominated emission [2109.02684], yet CALIFA-based analysis argues that \(\Sigma_{\mathrm{H}\alpha}\) is conceptually flawed because it behaves as an extensive projected quantity and can misclassify bright bulge DIG as star-forming while missing faint outer-disk star-forming regions [1711.07844]. The critique becomes even sharper in ram-pressure stripped systems: in GASP, neither a single H\(\alpha\) threshold, nor a single \([\mathrm{S\,II}]/\mathrm{H}\alpha\) threshold, nor an \(EW(\mathrm{H}\alpha)\) threshold was found to separate DIG-dominated from non-DIG emission reliably at kpc resolution [2011.08869].

More elaborate classifications have therefore been developed. GASP derives a spaxel-level DIG fraction \(C_{\mathrm{DIG}}\) by combining attenuation-corrected \(\Sigma H\alpha_{\rm corr}\) with metallicity-corrected \([\mathrm{S\,II}]/\mathrm{H}\alpha\), using
\[
C_{\mathrm{DIG}} = \left( \frac{f_0}{\Sigma H\alpha,corr} \right)^{\beta}
\qquad (\Sigma H\alpha_{\rm corr}>f_0),
\]
with galaxy-specific \(f_0\) and \(\beta\) [2011.08869]. BETIS instead combines adaptive binning based on \(S/N([\mathrm{S\,II}])=10\), morphological H II-region masks from pyHIIextractor, and an additional residual H\(\alpha\) surface-brightness threshold at \(3\sigma_{\Sigma(H\alpha)}\) to define lower and upper DIG limits over a wide range of spatial resolutions [2311.14254]. These developments underscore that DIG identification is operational and data-model dependent rather than fixed by a single universal observable.

## 4. Spatial structure, disk–halo connection, and kinematics

DIG is both a planar and an extraplanar phenomenon. In several edge-on spirals it forms vertically extended layers with kpc-scale structure, while in face-on systems it fills interarm and outer-disk regions and often dominates low-surface-brightness projected area [1304.7724][1711.07844]. CALIFA finds hDIG prevalent in ellipticals, S0s, bulges, and extraplanar regions, while the SF/mDIG proportion grows from early- to late-type spirals and from inner to outer radii [1711.07844].

Kinematic and structural modeling of NGC 4666 provides a detailed eDIG example. Its DIG is best fit by two ionized components: a thin disk with scale height \(0.2\) kpc and a thick disk with scale height \(0.8\) kpc, with a thin/thick flux ratio of 2.5, inclination \(i=78^\circ\), radial scale length \(3.0\) kpc, truncation radius \(11.5\) kpc, general linewidth \(35\,\mathrm{km\,s^{-1}}\), and maximum rotation velocity \(\sim195\,\mathrm{km\,s^{-1}}\) [1304.7724]. The same study identifies minor-axis line splitting with a component blueshifted by about \(100\,\mathrm{km\,s^{-1}}\), interpreted as outflowing ionized gas, and concludes that enhanced star formation both drives outflow and maintains a Reynolds-layer-like thick ionized disk [1304.7724].

The morphology of DIG is not universal. In NGC 4013 and NGC 4302 the extraplanar DIG is dominated by a smooth, diffuse component, while high-resolution dust images reveal strongly filamentary absorption structures with no counterpart in the H\(\alpha\) morphology. The conclusion is that the thick-disk DIG and the dusty extraplanar filaments trace physically distinct phases of the thick-disk ISM, with the latter representing denser neutral material [1301.0325].

Milky Way studies likewise show that DIG encodes both galactic structure and disk–halo transport. In the first Galactic quadrant, Green Bank Telescope radio recombination-line data reveal two dominant DIG velocity components centered around \(100\,\mathrm{km\,s^{-1}}\) and \(45\,\mathrm{km\,s^{-1}}\), with the higher-velocity component associated with W43 and the lower-velocity component interpreted either as gas at a different distance or as bar-driven streaming near W43 [1709.09232]. In the anti-center, a LAMOST sample of 17,821 DIG spectra shows that \([\mathrm{N\,II}]/\mathrm{H}\alpha\) and \([\mathrm{S\,II}]/\mathrm{H}\alpha\) are enhanced in the interarm region between the Local and Perseus arms near \(R_{\rm gal}\sim 9.1\) kpc, while \([\mathrm{N\,II}]/\mathrm{H}\alpha\), \([\mathrm{S\,II}]/\mathrm{H}\alpha\), and \([\mathrm{S\,II}]/[\mathrm{N\,II}]\) all increase with \(|z|\) [2412.05692].

Cluster environments add a further DIG regime. In stripped tails, DIG-dominated regions can lie at projected distances \(>2\) kpc and up to 10 kpc from star-forming regions, with high \(\Delta \log([\mathrm{O\,I}]/\mathrm{H}\alpha)\) and a high fraction of LIER/LINER-like emission [2108.12433]. This suggests that stripped tails probe ionization channels—mixing, shocks, and ICM/ISM interaction—that are harder to isolate in normal disks.

## 5. Consequences for diagnostic diagrams, metallicity, and star-formation rates

DIG directly affects optical diagnostic ratios. In MaNGA, DIG-dominated low-\(\Sigma_{\mathrm{H}\alpha}\) regions display enhanced \([\mathrm{S\,II}]/\mathrm{H}\alpha\), \([\mathrm{N\,II}]/\mathrm{H}\alpha\), \([\mathrm{O\,II}]/\mathrm{H}\beta\), and \([\mathrm{O\,I}]/\mathrm{H}\alpha\) at fixed metallicity, and DIG contamination moves H II regions toward composite or LI(N)ER-like areas of BPT space [1612.02000]. Standard H II-region photoionization grids fail to reproduce these DIG line ratios, and leaky H II-region models shift them only slightly, favoring a harder ionizing spectrum for DIG with LI(N)ER-like emission [1612.02000].

The impact on metallicity diagnostics is strongly calibrator-dependent. In the MaNGA-based DIG-correction framework built from 1,409 star-forming galaxies, the DIG effect is negligible for O3N2 but reaches \(\sim 0.1\) dex at the high-metallicity end for N2, and the real effect is expected to be larger than the measured one because MaNGA’s \(\sim\)kpc resolution still mixes H II regions and DIG [1907.08635]. Applying line-by-line DIG corrections to SDSS galaxies modifies the inferred \(M_\star\)–\(Z\)–SFR relation, especially at high stellar mass, and with the N2 indicator the corrected relation yields oxygen abundance increasing with SFR at high mass, contrary to previous claims [1907.08635].

Earlier MaNGA work reached a closely related but more categorical conclusion about calibrators. Metallicities derived using N2O2 are described as optimal because they exhibit the smallest bias and error, whereas O3N2, \(R_{23}\), N2, and the Dopita et al. (2016) N2S2H\(\alpha\) calibration can introduce biases in derived metallicity gradients as large as the gradient itself. The strong-line method IZI cannot be applied to DIG accurately because it contains only H II-region models [1612.02000]. In the same study, DIG-enhanced N2 can overestimate metallicity by roughly \(0.24\) dex in an illustrative case [1612.02000].

DIG also biases star-formation rates when total H\(\alpha\) is interpreted as massive-star emission. The 2020 review emphasizes that DIG can contribute a large fraction of total H\(\alpha\), so treating all H\(\alpha\) as star-formation emission distorts or overestimates the inferred SFR, especially when old stars ionize part of the gas [2005.06054]. In the GASP analysis, removing DIG lowers SFR estimates by about 0.2 dex, although the shift is similar in stripped and control samples [2011.08869]. DIG fractions in GASP disks span \(0.2 \lesssim C_{\mathrm{DIG,disk}} \lesssim 0.9\), and \(C_{\mathrm{DIG,disk}}\) anti-correlates with both sSFR and \(\Sigma\mathrm{SFR}\) [2011.08869].

The extent of DIG contamination in integrated spectra of bright star-forming galaxies remains debated. One line of argument emphasizes that DIG is a serious contaminant for metallicity, SFR, and AGN classification [2005.06054][1612.02000]. Another argues that in normal bright star-forming galaxies the effect on integrated spectra is smaller than previously claimed because much of the apparent discrepancy between local H II-region slit spectra and galaxy-integrated spectra is actually an aperture effect. Using MUSE data for 11 bright H II regions in three nearby galaxies, Mannucci et al. show that when large apertures are used while still masking DIG-dominated spaxels, \([\mathrm{N\,II}]/\mathrm{H}\alpha\) and \([\mathrm{O\,III}]/\mathrm{H}\beta\) change little, whereas sulfur ratios shift strongly and the DIG-free large-aperture H II-region spectra move into agreement with SDSS and MaNGA galaxy loci [2109.02684]. The strongest practical implication is that sulfur-based discrepancies do not by themselves prove DIG contamination; aperture mismatch and internal H II-region stratification must be considered first [2109.02684].

## 6. Numerical modeling, synthesis, and current debates

Radiation-hydrodynamic and post-processed radiative-transfer models now reproduce several canonical DIG properties, but they do so with differing physical emphases. Idealized radiation-hydrodynamic simulations of disk patches show that photoionization feedback can drive low levels of turbulence in the dense disk and supply thermal pressure support for an extended diffuse ionized layer, creating a natural coupling between the ionizing photon budget and the mass in different ionization phases [1907.02067]. In that framework, photoionization raises the gas temperature from \(500\) K to \(8000\) K and lowers the mean molecular weight, thereby increasing pressure support enough to build a warm quasi-hydrostatic layer [1907.02067].

Cosmic rays address a related but distinct problem: maintaining the large vertical scale heights of DIG. Post-processing of SILCC simulations shows that cosmic-ray feedback produces more extended gaseous disks and H\(\alpha\) DIG scale heights of \(0.643\)–\(0.723\) kpc, far larger than thermal-feedback-only models and closer to observed values. However, adding a fiducial cosmic-ray heating term increases temperature with height but fails to reproduce observed nitrogen and sulfur forbidden-line intensities; the required heating must affect gas over a broader density range, or the total ionizing luminosity must be fine-tuned so that the ionizing spectrum hardens appropriately with height [1802.07749].

A newer isolated-galaxy simulation with on-the-fly radiative transfer and non-equilibrium thermochemistry reaches a different synthesis. It reproduces observed correlations of \([\mathrm{S\,II}]/\mathrm{H}\alpha\), \([\mathrm{N\,II}]/\mathrm{H}\alpha\), \([\mathrm{O\,I}]/\mathrm{H}\alpha\), and \([\mathrm{O\,III}]/\mathrm{H}\beta\) with \(\Sigma_{\mathrm{H}\alpha}\), and it attributes these trends to increasing temperature and a hardening radiation field with decreasing \(n_{\rm e}\). In that model the DIG contributes only 26% of the intrinsic H\(\alpha\) but 60% of the observed H\(\alpha\), because about \(\sim 80\%\) of H\(\alpha\) from H II regions is absorbed by dust versus only \(\sim 10\%\) for DIG emission [2403.03243]. This suggests that observational DIG prominence can be amplified by differential attenuation even when its intrinsic recombination budget is smaller.

Current observational synthesis remains two-sided rather than uniform. PHANGS-MUSE argues for “two DIGs,” with leaked H II-region radiation dominating the DIG H\(\alpha\) energy budget and HOLMES required for harder line-ratio phenomena [2111.14876]. The 2024 simulation instead argues that ongoing star formation, specifically stars aged \(5\)–\(25\) Myr, can account for both the low-ionization and the \([\mathrm{O\,III}]/\mathrm{H}\beta\) trends without invoking secondary ionization sources in normal star-forming disks [2403.03243]. A plausible synthesis is that the dominant ionization channel depends on galaxy type and environment: young-star leakage and exposed intermediate-age populations in actively star-forming disks, HOLMES in bulges and low-EW systems, and shocks or mixing in outflows and stripped tails [2109.02684][2108.12433].

The main conceptual convergence is that DIG is neither negligible nor homogeneous. It is a structurally extended, spectroscopically distinct, and environmentally sensitive component of the ionized interstellar medium. Its observational definition remains operational, its inferred ionization source depends on spatial scale and host environment, and its impact on nebular diagnostics ranges from secondary to dominant depending on whether the target is a bright star-forming disk, a bulge-dominated system, or a low-surface-brightness extraplanar or stripped-gas structure [2005.06054][2109.02684].

Source: https://www.emergentmind.com/topics/diffuse-ionized-gas-dig