Star-Formation-Driven Ionization Scenario

• Star-Formation-Driven Ionization Scenario is a framework that links massive star formation to the generation of ionizing photons, shaping nebular emission and cosmic reionization.
• The approach employs distributions of ionizing flux and hydrogen density through models like the LOC formalism to interpret emission lines and feedback effects.
• Observational diagnostics and simulations validate the scenario by linking strong-line ratios, cosmic-ray ionization, and feedback-regulated starbursts to evolving intergalactic conditions.

The star-formation-driven ionization scenario encompasses a suite of astrophysical processes and models in which the ionization state of the interstellar and circumgalactic medium is regulated, and in many cases dominated, by ionizing photons, mechanical feedback, and cosmic rays originating from sites of active star formation. This framework provides a physical basis for interpreting nebular emission lines, thermal and non-thermal feedback signatures, the evolution of the ionized volume filling factor within and around galaxies, and the reionization of the intergalactic medium. Star-formation-driven ionization plays a critical role from small-scale triggering in molecular clouds and bright-rimmed clumps up to the statistical properties of emission-line galaxies during cosmic dawn.

1. Theoretical Foundations and Formulation

At its core, the star-formation-driven ionization scenario links the production of ionizing photons (primarily Lyman continuum, $\lambda < 912$ Å) to the formation and evolution of massive stars and the resulting stellar populations. The canonical ionization parameter is defined as

$U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$

where $Q(\mathrm{H^0})$ is the rate of hydrogen-ionizing photons, $R$ the characteristic radius of the H II region, $n_H$ the hydrogen density, and $c$ the speed of light. For a volume-integrated ensemble, $U$ tracks the ratio of ionizing photon density to hydrogen-atom density, capturing the interplay between stellar birthrate and gas reservoir (Richardson et al., 2016, Kaasinen et al., 2018, Shen et al., 2024).

Modern models replace the simplistic assumption of a single $U$ per galaxy with distributions over both $n_H$ and incident ionizing photon flux $\Phi_H$, as in the Locally Optimally Emitting Cloud (LOC) formalism: $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$0 yielding total nebular emission

$U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$1

where $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$2 is the Cloudy-predicted emissivity. Along emission-line galaxy sequences (e.g. the BPT diagram), $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$3 varies systematically with galaxy ionization level, with metallicity as a secondary parameter (Richardson et al., 2016).

For global ionization of the intergalactic medium (IGM), the mass-weighted ionized fraction $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$4 evolves per: $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$5 with $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$6 integrated $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$7 over the population, and $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$8 the recombination timescale (Hartley et al., 2016, Chakraborty et al., 2024, Stanway et al., 2015).

2. Star Formation, Stellar Populations, and the Ionizing Photon Budget

The ionizing photon production rate per unit star-formation rate, $U \equiv \frac{Q(\mathrm{H^0})}{4\pi R^2 n_H c}$9, depends sensitively on metallicity, stellar population age, initial mass function (IMF), and binary evolution physics:

• For continuous star formation, $Q(\mathrm{H^0})$0 at $Q(\mathrm{H^0})$1, but rises to $Q(\mathrm{H^0})$2 at $Q(\mathrm{H^0})$3.
• Rapidly rotating and interacting binaries further boost $Q(\mathrm{H^0})$4 by $Q(\mathrm{H^0})$560% at $Q(\mathrm{H^0})$6 (Stanway et al., 2015, 1803.02340).
• In dwarf galaxies and metal-poor starbursts, short ($Q(\mathrm{H^0})$7 Myr), efficient ($Q(\mathrm{H^0})$8–$Q(\mathrm{H^0})$9\% mass converted) bursts generate $R$0–$R$1 s$R$2, $R$3–$R$4, and high $R$5; these conditions produce strong forbidden and recombination lines (Martín-Manjón et al., 2011, 1803.02340).
• The observed UV continuum and line emission allow calibration of $R$6, escape fractions ($R$7), and photon budgets required for reionization and feedback (Chakraborty et al., 2024, Stanway et al., 2015).

3. Ionization in Galactic, Circumgalactic, and Extragalactic Contexts

Diffuse Ionized Gas and Galaxy Disks

Ionizing photons from massive stars maintain the warm ionized medium (WIM/DIG) in spiral disks. MHD/radiative transfer simulations confirm that realistic star formation rates ($R$8) can sustain a WIM with volume filling factor $R$9 at $n_H$0 kpc, matching observed H$n_H$1 emission and electron density profiles. The photon escape fraction to halo scales is low ($n_H$2) (Kado-Fong et al., 2020). A clumping factor $n_H$3 corrects for inhomogeneity in emission measure analyses.

Global Reionization

High-redshift galaxies dominate the ionizing emissivity of the IGM. Semi-analytical models constrained by the UV luminosity function, electron-scattering optical depth $n_H$4, and the evolving neutral fraction require:

• Rapidly rising star-formation efficiencies at $n_H$5
• Halo mass–dependent escape fractions ($n_H$6 with $n_H$7)
• Early reionization completed by $n_H$8
• Cumulative contributions dominated by galaxies fainter than $n_H$9 (Hartley et al., 2016, Chakraborty et al., 2024).

Bursty star-formation histories, with short ($c$0 Myr) intense episodes, produce larger H II regions and long-lived relics that maintain partial ionization, raising $c$1 for a given $c$2 (Hartley et al., 2016).

4. Feedback, Outflows, and the Regulation of Ionization

Stellar feedback—via radiation pressure, mechanical energy, and photoionization—drives complex structures (expanding bubbles, pillars, outflows) within star-forming environments:

• In massive compact systems at high $c$3, outflows of ionized gas are directly powered by stellar feedback, with no evidence for AGN contribution (e.g., EGSY8p7/CEERS-1019), and mass-loading factors $c$4, electron densities $c$5 cm$c$6, and $c$7 (Zamora et al., 9 Dec 2025).
• In local giant H II regions (e.g., SMC-N66), photoionization by known O stars balances the observed recombination rate, with $c$8 leakage through low-density channels (Geist et al., 2022).
• High-resolution ALMA studies of bright-rimmed clouds reveal triggering of sequential massive star formation via radiation-driven implosion (RDI), with overpressured ionized boundary layers and coeval filamentary inflow (“clump-fed” scenarios) mediating mass growth in cores (Zhang et al., 2023, Bisbas et al., 2010).
• Numerical simulations of ionization-induced star formation in bound and unbound clusters establish that while feedback lowers global star-formation efficiencies, it induces local triggering (10–40% of the final stellar mass) at bubble rims and pillar tips, and imprints lower-mass stellar initial mass functions (Dale et al., 2012, Dale et al., 2013).
• In high-mass star forming regions, photoionization alone can launch molecular outflows, but with much lower momentum than observed; observed outflows likely require collective low/intermediate-mass jet activity (Klaassen et al., 2013, Peters et al., 2012).

5. Observational Diagnostics and Empirical Trends

Empirical studies across a range of redshifts leverage strong-line ratios and broadband indicators to diagnose ionization conditions:

• The [O III]/[O II] ratio ($c$9), a robust $U$0-tracer, correlates tightly with sSFR, EW(H$U$1), and SFR surface density. At fixed mass or SFR, galaxies at $U$2–$U$3 exhibit higher $U$4 than local counterparts. High $U$5 and hard radiation fields are generic in early star formation (Shen et al., 2024, Kaasinen et al., 2018).
• The correlation $U$6 holds, even at constant metallicity, indicating that increased sSFR—not simply lower $U$7—drives elevated ionization in high-$U$8 main sequence galaxies (Kaasinen et al., 2018).
• Extreme high-ionization rest-frame UV line emission in metal-poor galaxies (C IV, O III], He II) is explained by short ($U$9 Myr) binary-rich starburst models, ruling out significant contributions from AGN and fast shocks (1803.02340).

6. Ionization by Cosmic Rays in Star Formation Environments

In starburst galaxy nuclei (e.g., NGC 253), cosmic rays—created predominantly by supernovae accompanying vigorous star formation—establish extreme ionization rates ($U$0–$U$1, $U$2–$U$3 Milky Way values). The [H$U$4O$U$5]/[SO] abundance ratio is a direct probe: chemical modeling and ALMA observations exclude pure PDR, XDR, or shock-only mechanisms (Holdship et al., 2022). Such high CRIR modifies the chemistry, thermal balance, and magnetic coupling of dense gas, with implications for star-formation efficiency and feedback regulation.

7. Extensions, Limitations, and Model Comparisons

LOC models and bursty SF prescriptions unify emission line phenomenology throughout the BPT locus and explain the majority of observed variation with continuous variations in the distribution of ionizing flux at fixed (or mildly variable) metallicity (Richardson et al., 2016). However, at the extreme ends of the galaxy ionization sequence, additional physical parameters (dust, cloud geometry, and SED) must be invoked. The cosmic reionization history can only be reconciled with CMB and quasar absorption constraints by invoking time-variable, efficiency-enhanced, or halo-mass-sensitive escape fractions, and a population-averaged photon budget shaped by both sustained and bursty star formation (Hartley et al., 2016, Chakraborty et al., 2024). In molecular cloud outflow contexts, pure ionization-driven models cannot account for the dynamical range seen in massive star-forming regions; the collective feedback of cluster members is required (Peters et al., 2012).

---

In summary, the star-formation-driven ionization scenario quantitatively and qualitatively explains the excitation and dynamical state of H II regions, diffuse ionized parameters, cosmic reionization, and a variety of feedback-driven structures in star-forming environments. This framework is robustly supported by emission-line, continuum, and dynamical observations; self-consistent theoretical/numerical models; and detailed spectroscopic diagnostics across a broad range of cosmic environments (Richardson et al., 2016, Hartley et al., 2016, 1803.02340, Stanway et al., 2015, Geist et al., 2022, Martín-Manjón et al., 2011, Bisbas et al., 2010, Dale et al., 2012, Dale et al., 2013, Zamora et al., 9 Dec 2025, Chakraborty et al., 2024, Kado-Fong et al., 2020, Zhang et al., 2023, Shen et al., 2024, Klaassen et al., 2013, Kaasinen et al., 2018, Holdship et al., 2022, Peters et al., 2012).