Fact-checked by Grok 2 weeks ago

Reionization

Cosmic reionization is the in the early during which the neutral intergalactic medium (IGM), primarily composed of and , was ionized by (UV) and radiation from the first luminous sources, marking the last major in cosmic history. This process transformed the from an opaque, neutral state to a transparent, ionized one, allowing photons to propagate freely across cosmic distances. The Epoch of Reionization (EoR) spanned approximately from z ≈ 15 to z ≈ 6, corresponding to 200–900 million years after the , though recent observations suggest the onset may have begun as early as z ≈ 13, around 330 million years post-. It followed the cosmic "Dark Ages," a period after recombination (z ≈ 1100) when the universe was neutral and devoid of , and preceded the formation of the familiar large-scale observed today. The reionization process was extended and inhomogeneous, starting in underdense regions ("outside-in" mode) where ionizing photons could travel farther before being absorbed, forming expanding H II bubbles that eventually overlapped to ionize the entire IGM by z ≈ 6. The primary drivers of reionization were UV photons (with energies >13.6 eV) escaping from star-forming galaxies, particularly low-mass systems at high redshifts, which produced the necessary ionizing flux to balance recombinations in the clumpy IGM. Contributions from active galactic nuclei (AGNs), including quasars and accreting supermassive black holes, provided harder radiation that could penetrate denser regions, though their role remains secondary to stellar sources based on current luminosity function estimates. The escape fraction of ionizing photons (fesc) from these early galaxies is a key uncertain parameter, estimated at 5–20%, influencing the pace and topology of reionization. Observationally, the end of reionization is evidenced by the sharp increase in transmitted flux in the Lyman-α forest of spectra at z < 6, indicating a highly ionized IGM. The cosmic microwave background (CMB) provides an integrated measure through the optical depth to electron scattering, τ = 0.054 ± 0.007, from Planck 2018 measurements, implying reionization began before z ≈ 8. Recent James Webb Space Telescope (JWST) detections of galaxies at z > 10, such as those showing Lyman-α emission at z ≈ 13, are revealing the first ionizing sources and constraining the neutral hydrogen fraction. Future probes, including 21 cm intensity mapping with facilities like the (SKA), will map the spatial progression of reionization in detail.

Concept

Definition and Physical Process

Reionization refers to the in the early during which the neutral in the intergalactic medium (IGM), formed after recombination, becomes predominantly ionized. This process occurred roughly between redshifts z \approx 6 and z \approx 15, transforming the IGM from a neutral state into a highly ionized warm . The primary physical mechanism driving reionization is photoionization, where ultraviolet photons with energies exceeding the hydrogen ionization threshold of 13.6 interact with neutral atoms, ejecting electrons and producing free protons. This process can be described by the photoionization rate per neutral atom, C_{\rm HI} = 4\pi \int_{\nu_{\rm HI}}^\infty \frac{J_\nu}{h\nu} \sigma_{\rm HI}(\nu) \, d\nu, where J_\nu is the specific intensity of the radiation field, \sigma_{\rm HI}(\nu) is the cross-section, and \nu_{\rm HI} corresponds to 13.6 . The resulting free electrons and protons maintain a dynamic balance influenced by recombination, but the system is generally out of due to the influx of . In regions approaching ionization equilibrium, the Saha equation provides a framework for estimating the fraction, relating the densities of ionized and states through and . For , it takes the form \frac{n_{\rm HII} n_e}{n_{\rm HI}} = \left( \frac{2\pi m_e k_B T}{h^2} \right)^{3/2} \exp\left( -\frac{I_{\rm H}}{k_B T} \right), where n_{\rm HII}, n_{\rm HI}, and n_e are the number densities of protons, , and electrons, respectively; I_{\rm H} = 13.6 is the . However, this equation assumes local , which breaks down during reionization due to non-equilibrium conditions at propagating ionization fronts, where rapid and dynamic IGM evolution prevent full thermal balance. Locally, creates ionized bubbles analogous to Strömgren spheres around discrete sources, where the sphere's radius balances the production and recombination of ionizations within a uniform medium. The propagation of these bubbles is limited by the high \tau to ionizing photons, which arises from the large cross-section of and results in rapid over short mean free paths in the IGM. While primarily driven by , the process also sets the stage for reionization at lower redshifts (z ≈ 3), completing IGM ionization. Reionization concluded around z \approx 6, as evidenced by the emergence of the Gunn-Peterson trough in spectra, indicating a sharp decline in density.

Importance in Cosmic History

Reionization marked a transformative in the universe's thermal history, as the of in the intergalactic medium (IGM) by the first luminous sources heated the gas initially to ~20,000–30,000 K via the during the passage of ionization fronts, which cools to ~10^4 K by the end of the due to adiabatic expansion. The elevated temperatures increased the Jeans mass, providing thermal feedback that suppressed the collapse of gas into low-mass halos and inhibited in dwarf galaxies, thereby regulating small-scale and contributing to the scarcity of observed faint satellites in the local universe. The of the IGM during reionization fundamentally altered the of , transitioning the from an opaque, state to one where photons could travel freely beyond the wavelength. In a universe, the Gunn-Peterson effect would absorb nearly all shortward of , rendering high-redshift sources invisible; post-reionization, the ionized IGM enables the transmission of spectral features from distant galaxies, allowing astronomers to probe the early directly. This shift not only unveiled the high-redshift but also facilitated observations of emitters, whose visibility serves as a tracer of ionized bubbles amid lingering patches. On a broader scale, reionization bridged the cosmic dark ages—following recombination—to the era of modern , influencing the evolution of the large-scale structure through its effects on the IGM's and state and the observability of baryonic acoustic oscillations in post-reionization probes like the 21 cm line and Lyman-α forest. Additionally, the free electrons generated during reionization scattered (CMB) photons via , contributing the bulk of the measured \tau \approx 0.054 (as of Planck 2018), which generates a large-scale signal in the CMB and provides a key constraint on the timing of this epoch.

Timeline and Phases

Recombination and Dark Ages

Following the and , the continued to expand and cool, reaching a of approximately 0.3 at a z \approx 1100, where the thermal energy became comparable to the 13.6 of neutral . At this point, known as the recombination epoch, free electrons and protons began to recombine efficiently to form neutral atoms, transitioning the from a fully ionized to a mostly neutral state. This process was governed by the Saha equation, which balances the ionization fraction based on , , and the , though detailed calculations account for non-equilibrium effects such as the finite and feedback from photons. Recombination occurred over a narrow range, \Delta z \approx 200, lasting about 200,000 years, during which the electron fraction dropped from nearly 1 to below 10^{-4}. The recombination process is characterized by the visibility function, g(z) = -\frac{d e^{-\tau}}{dz}, where \tau(z) is the due to Thomson scattering of photons by free electrons. This function peaks sharply at z \approx 1090, marking the moment when the universe became optically thin to photons, allowing them to from baryons and free-stream as the (CMB). Prior to decoupling, photons maintained tight coupling with electrons via , damping baryonic density perturbations on small scales, while after decoupling, the universe entered a matter-dominated era where baryons fell into potential wells, enabling further structure growth. The period immediately following recombination, spanning redshifts from z \approx 1100 to z \approx 30, is termed the cosmic dark ages, during which the universe was filled with neutral and gas in a homogeneous, opaque state devoid of luminous sources. This era lasted approximately 100 million years, with the universe expanding and cooling adiabatically, reaching gas temperatures around 80 K by z \approx 30, closely coupled to the temperature. Without stars or galaxies, the intergalactic medium remained dark to beyond the CMB, though primordial density perturbations—seeded by quantum fluctuations during —grew linearly under gravity in the matter-dominated regime, following the \delta \propto a scaling where a is the scale factor. During the dark ages, the first gravitationally bound structures, known as minihaloes, began to form at redshifts z \approx 20-30, with typical masses of $10^5 - 10^6 M_\odot arising from haloes that captured baryonic gas. These minihaloes cooled primarily via molecular emission but lacked sufficient mass and temperature to ignite massive stars capable of producing significant , thus preserving the neutrality of the intergalactic medium. The end of the dark ages was marked by the collapse of more massive haloes (> 10^6 M_\odot) around z \approx 20, where the first Population III stars formed, initiating the cosmic dawn and the onset of reionization.

Cosmic Dawn and Epoch of Reionization

The Cosmic Dawn marks the initial phase of reionization at redshifts z \approx 20-30, when the first Population III stars formed within the smallest halos, heralding the end of the neutral Dark Ages and the onset of cosmic structure illumination. Recent (JWST) observations suggest the onset of reionization may have begun as early as z \approx 13. These primordial stars, along with the earliest galaxies, emitted photons that began pockets of the intergalactic medium (IGM), forming small H II regions or "bubbles" around their host halos. This process initiated the transition from a fully neutral to one increasingly permeated by ionized , with the volume-averaged ionization fraction x_e starting from near zero and gradually rising as more sources ignited. The Epoch of Reionization (EoR) followed, spanning z \approx 6-15 over approximately 700 million years, during which the discrete nature of early sources led to patchy reionization characterized by inhomogeneous x_e . Ionization bubbles expanded and overlapped, percolating through the IGM and achieving roughly 50% global ionization by z \approx 8, as the ionized grew nonlinearly due to the clustering of sources in overdense regions. Two primary scenarios describe this progression: an inside-out mode, where denser regions ionize first around galaxies before expanding outward, versus an outside-in mode, where lower-density voids ionize earlier due to longer mean free paths for photons. Reionization culminated near z \approx 6, when the IGM became fully ionized, as indicated by the sharp decline in neutral hydrogen absorption in high-redshift spectra. Helium reionization occurred later and more abruptly, with He II ionizing primarily at z \approx 3 due to harder photons from quasars, distinct from the extended process. This phase left a legacy of patchy topology in the IGM, with x_e fluctuations persisting briefly after overlap, reflecting the finite number and distribution of discrete ionizing sources.

Ionizing Sources

Population III Stars

Population III stars represent the first generation of stars to form in the , emerging from pristine, metal-free gas clouds in small minihalos with virial masses of $10^5 to $10^6 masses at redshifts z > 20. These structures collapsed under their own , with cooling primarily facilitated by molecular (H₂) rovibrational line emission, enabling the fragmentation and subsequent formation of protostellar cores. The lack of heavy elements, which in later stellar generations promote efficient and smaller-scale fragmentation, resulted in an (IMF) heavily skewed toward massive stars, with characteristic masses spanning 10 to 1000 masses. This top-heavy IMF arises because the absence of metals limits fragmentation during collapse, favoring the accretion of larger gas reservoirs onto fewer, more massive protostars. These metal-free stars exhibited extreme physical properties, including effective temperatures exceeding $10^5 K, which produced a hard ultraviolet spectrum rich in high-energy photons capable of ionizing hydrogen and helium. For stars with initial masses above approximately 100 solar masses, this spectral hardness stemmed from their high surface temperatures and lack of line blanketing by metals, enhancing the production of photons with energies above 13.6 eV. Their lifetimes were exceptionally brief, typically lasting only a few million years—for instance, a 200 solar mass star evolves through its main sequence in about 2.4 million years—due to their high masses and rapid nuclear burning rates. At the end of their lives, many Population III stars in the mass range of 140 to 260 solar masses underwent pair-instability supernovae, explosive events triggered by electron-positron pair production in their cores, which completely disrupted the stars without leaving remnants. In the context of early cosmic reionization, Population III played a pivotal role during cosmic dawn by serving as the primary sources of , generating small, localized H II regions or bubbles around their minihalos. The escape fraction of these ionizing from the dense environments was relatively low, estimated at 10–20%, limited by the high gas densities and lack of pre-existing channels in the metal-free gas. Despite their efficiency in producing hard UV , the overall photon budget from Population III proves insufficient to drive the full reionization of the intergalactic medium, as their short-lived nature and limited numbers in minihalos contribute only modestly to the global ionization history. The explosions of these massive , however, injected the first metals into the intergalactic medium, enriching it to metallicities of $10^{-6} to $10^{-4} solar values and transitioning subsequent to metal-enriched Population II modes.

Early Galaxies

Early galaxies, which began forming at redshifts [z](/page/Z) \gtrsim 10, are compact systems with low metallicities typically below [Z](/page/Z) \approx 0.1 Z_\odot, fostering environments conducive to intense, bursty driven by rapid gas accretion onto halos of masses $10^8 - 10^{10} M_\odot. These galaxies exhibit irregular star formation histories characterized by short bursts interspersed with quiescence, influenced by the hierarchical merging of progenitors and limited metal enrichment from prior generations. The ultraviolet (UV) continuum emission from these galaxies traces their star formation and is described by the Schechter luminosity function, with faint-end slope \alpha \approx -2.0, characteristic magnitude M^* \approx -20.5 at [z](/page/Z) \sim 7, and normalization \phi^* \approx 10^{-3} Mpc^{-3} mag^{-1}, indicating a steep increase in the abundance of low-mass systems toward higher redshifts. This distribution highlights their role as numerous, small contributors to the cosmic UV background during reionization. The ionizing photon output from early galaxies primarily arises from Population II stars, which dominate the stellar mass in these metal-enriched systems and produce spectra softer than those of metal-free Population III stars, with production efficiencies \xi_{\rm ion} \approx 10^{25.5} Hz erg^{-1}. Low dust content in these young galaxies, due to inefficient enrichment and limited grain formation, enables high Lyman continuum escape fractions f_{\rm esc} > 20\%, allowing a substantial portion of ionizing photons to permeate the . Models integrating the UV function suggest that these galaxies supplied approximately 70-80% of the ionizing photons required for reionization, outpacing contributions from quasars and underscoring their dominance in sustaining the expanding ionized volume. Radiative feedback from young stars ionizes surrounding gas, creating H II regions that expand and regulate further accretion, while supernova explosions from massive stars inject momentum and energy, dispersing gas clouds and suppressing in low-mass halos through outflows reaching velocities of hundreds of km/s. These processes facilitate the overlap of ionized bubbles around neighboring galaxies, accelerating the transition to a fully ionized intergalactic medium by enhancing propagation and reducing neutral absorption. Recent (JWST) observations from 2024-2025 have identified thousands of such galaxies at z > 10, including catalogs exceeding 1,000 candidates up to z \sim 18, confirming their prevalence and bolstering models of galaxy-driven reionization.

Quasars and Active Galactic Nuclei

Quasars and active galactic nuclei (AGNs) are powered by accretion onto supermassive black holes (SMBHs), which form seeds at redshifts z > 10 through mechanisms such as the direct collapse of pristine, atomically cooled gas clouds or remnants of III stars, subsequently growing rapidly via gas accretion and mergers to reach masses of 10^8–10^9 solar masses by the end of reionization. These SMBHs achieve bolometric luminosities up to 10^47 erg/s, enabling them to emit intense radiation capable of influencing the intergalactic medium (IGM) during the late stages of cosmic reionization. The of quasars features a hard extending from (UV) to wavelengths, producing with exceeding 1 keV that can penetrate dense, neutral gas more effectively than softer stellar . This hard allows quasars to contribute approximately 10% to the total ionizing budget for reionization at z ≈ 6, with their role increasing to dominate (HeII) reionization due to the higher (54.4 ) required to ionize singly ionized . Overall, quasars provide less than 7–20% of the needed to sustain ionization in the IGM, underscoring their secondary but non-negligible influence compared to stellar sources. At redshifts z < 7, quasars play a key role in carving out large-scale ionized bubbles in the IGM, with their high luminosity ionizing volumes up to several megaparsecs across and facilitating the overlap of HII regions toward the end of reionization. Additionally, quasar-driven outflows and winds, reaching velocities of thousands of km/s, exert negative feedback by expelling gas from host galaxies, thereby suppressing subsequent star formation and regulating galaxy growth during this epoch. Quasars are rare at high redshifts due to the limited time available for SMBH growth, with over 300 confirmed examples known at z > 6 as of 2025, primarily selected from optical and near-infrared surveys. Recent (JWST) observations have identified candidate quasars or AGN at z ≈ 10, including obscured systems in luminous galaxies, providing new insights into the earliest SMBH activity and their potential contributions to pre-reionization .

Observational Probes

Quasar Spectra and Gunn-Peterson Trough

High-redshift quasars serve as backlights for probing the intergalactic medium (IGM) through absorption in their spectra, particularly by neutral hydrogen (H I) via the Lyman-alpha (Lyα) resonance line at 1216 Å. As quasar light travels through the IGM, neutral hydrogen atoms absorb photons at the redshifted Lyα wavelength, creating a series of absorption lines known as the Lyα forest at lower redshifts. At z > 6, when the ionized fraction x_e falls below unity during the final stages of reionization, this absorption becomes nearly complete, manifesting as the Gunn-Peterson trough—a broad, deep absorption feature in the quasar continuum spectrum blueward of the Lyα emission line. This trough arises from the damping wing of the Lyα line profile, where resonant scattering by neutral hydrogen suppresses transmission over a wide wavelength range, providing direct evidence of residual neutral gas. The depth of the Gunn-Peterson trough enables quantitative measurement of the in the IGM. The effective τ_eff is derived from the observed transmission T_obs, where T_obs ≈ e^{-τ_eff} in regions of complete , though actual spectra show sporadic transmission spikes due to ionized bubbles. For a uniform IGM, the Gunn-Peterson is approximated as \tau_{\rm GP} = \frac{n_{\rm HI} \sigma_\alpha (1+z)^2}{H(z)}, where n_HI is the proper of , σ_α ≈ 4.5 × 10^{-18} cm² is the Lyα cross-section at line center, and H(z) is the Hubble parameter at redshift z. In practice, during patchy reionization, τ_eff is related to the volume-averaged ⟨x_HI⟩ via models incorporating IGM inhomogeneities, with τ_eff ≈ 2.5–6 corresponding to ⟨x_HI⟩ ≈ 0.01–0.1 at z ≈ 6. Spectra from surveys like the (SDSS) are analyzed by fitting the trough profile, accounting for the quasar's intrinsic spectrum and foreground , to infer the evolution of x_HI. Key observations from SDSS quasars at z ≈ 5.7–6.4 reveal a sharp decline in Lyα transmission at z ≈ 6.2, marking the end of reionization, with the mean transmitted flux dropping from T ≈ 0.05 at z < 6 to T < 10^{-3} in complete troughs at z > 6.2, implying a volume-averaged ⟨x_HI⟩ ≈ 0.01–0.04 just after reionization completion. At higher redshifts (z > 6.5), the troughs exhibit patchy absorption with occasional transmission spikes (T ≈ 10^{-5}–10^{-4}), consistent with an incomplete, inhomogeneous reionization process where ionized regions coexist with pockets. Recent (JWST) spectra of z > 6 quasars confirm these findings, detecting Gunn-Peterson troughs with similarly low transmission levels (T ≈ 10^{-5}) and providing higher-resolution views of the damping wings, reinforcing that reionization concluded around z ≈ 6 while highlighting small-scale IGM variations.

Cosmic Microwave Background Signals

The (CMB) provides key constraints on cosmic reionization through the effects of between CMB photons and free electrons produced during this epoch. This scattering process damps small-scale temperature anisotropies in the CMB by diffusing photon directions, while simultaneously generating large-scale E-mode from the quadrupole moment of the field. The resulting signal is particularly prominent at low multipoles (low-ℓ), where a characteristic "bump" in the E-mode power spectrum arises from scattering at redshifts corresponding to reionization. The primary observable from these interactions is the Thomson optical depth τ_e, which quantifies the integrated probability of scattering along the from the last scattering surface to . This is defined as \tau_e = \int n_e(z) \sigma_T \frac{dr}{dz} dz, where n_e(z) is the free electron density at z, σ_T is the Thomson cross-section, and dr/dz is the radial comoving distance element. Measurements from the Planck satellite yield τ_e ≈ 0.054, reflecting the cumulative fraction over reionization history. Updated analyses in 2024, incorporating refined foreground subtraction and polarization data, maintain this value within uncertainties of about 0.007. These CMB signals constrain the timing and duration of reionization, indicating a relatively late midpoint at z_reion ≈ 7–8 to match the position and amplitude of the low-ℓ polarization bump. However, this inference creates tensions with observations of high-redshift galaxies from the (JWST), which suggest more abundant ionizing sources at z > 10, implying an earlier or more extended reionization phase that would produce a higher τ_e. Recent 2025 analyses of data, combined with quasar absorption spectra, support an extended tail in helium reionization extending to z < 6, allowing partial reconciliation by distributing the optical depth over a broader redshift range.

Lyman-Alpha Emission and Forest

The Lyman-alpha forest consists of a dense series of absorption lines in the spectra of distant quasars and galaxies, primarily arising from neutral hydrogen (HI) in the (IGM), which traces the density fluctuations of the IGM along the line of sight. These absorption features appear blueward of the at 1216 Å, with the forest's opacity increasing at higher redshifts due to the rising neutral hydrogen fraction during the (EoR). Transmission spikes within the forest—regions of higher flux transmission—correspond to underdense, ionized voids in the IGM, where the neutral fraction is low, allowing to propagate more freely and providing a direct probe of the topology of reionized bubbles. In the extreme case of complete absorption, the forest transitions to the , marking regions of high neutrality. Lyman-alpha emitters (LAEs) are galaxies detected through their prominent Lyman-alpha emission line, which is resonant scattering of ultraviolet photons by neutral hydrogen in the galaxy's interstellar medium and surrounding IGM, making them key tracers of ionized regions during reionization. Narrow-band imaging surveys, such as those from the Subaru Intense Lyman-alpha Survey (SILVERRUSH) and its extensions, have identified thousands of LAEs at redshifts z ≈ 5–9 by isolating the redshifted Lyman-alpha line in narrow filters, revealing luminosity functions that decline toward higher redshifts as neutral gas increasingly attenuates the emission. The damping wing—a broad, redshifted absorption feature in LAE spectra—arises from resonant scattering by residual neutral hydrogen in the surrounding IGM, with its shape and extent probing the column density of neutral gas and the size of local ionized bubbles, typically on scales of 0.1–1 proper Mpc. For instance, models of damping wings in high-z spectra indicate that visible LAEs must reside within ionized bubbles to escape the damping absorption, constraining the neutral hydrogen fraction to below ~20% in their vicinity. Key metrics from LAE observations include the ionized coverage fraction f_{\rm cov}, which quantifies the volume fraction of the IGM that is ionized and allows LAE visibility, and the characteristic bubble size inferred from clustering and transmission statistics, often estimated at ~1–10 Mpc^3 comoving volume for z ≈ 6–8 sources. Recent James Webb Space Telescope (JWST) observations in the Cosmic Evolution Early Release Science (CEERS) survey reveal a sharp drop in the LAE fraction at z > 7, with only ~10% of UV-selected galaxies showing strong Lyman-alpha emission (equivalent widths > 50 Å) compared to ~40% at z ≈ 6, attributed to increasing IGM neutrality and smaller ionized bubbles around fainter galaxies. A 2025 study using JWST/NIRSpec spectroscopy of a galaxy at z ≈ 13 demonstrates the onset of reionization through Lyman-α emission, revealing a local ionized bubble of ~0.2 proper Mpc amid a mostly neutral IGM, with an inferred LyC escape fraction near unity. These findings highlight LAEs as sensitive indicators of bubble growth, with visibility requiring both high intrinsic production rates and local ionization to overcome damping.

21-cm Hydrogen Line

The 21-cm hydrogen line arises from the hyperfine spin-flip transition between the parallel and antiparallel spin states of the electron and proton in neutral hydrogen atoms, with a rest-frame frequency of 1420.4057 MHz corresponding to a wavelength of 21 cm. This transition probes the neutral intergalactic medium (IGM) during the cosmic dark ages and the epoch of reionization, as the signal is emitted or absorbed by neutral hydrogen (HI) before it is ionized by ultraviolet photons from the first stars and galaxies. Due to cosmological redshift, observations at redshifts z ≈ 6–20 shift the frequency to the low-radio band of 70–200 MHz, enabling ground-based detection despite foreground contamination from galactic and extragalactic synchrotron emission. The primary observable is the differential brightness temperature δT_b, which quantifies the contrast between the 21-cm signal and the (CMB) in the Rayleigh-Jeans tail of the spectrum. In the limit of low , the is approximated as \delta T_b \approx 27 \, \mathrm{mK} \left(1 - \frac{T_{\mathrm{CMB}}}{T_s}\right) x_{\mathrm{HI}} (1 + \delta_b), where T_{\mathrm{CMB}} \approx 2.725 (1+z) K is the CMB temperature, T_s is the spin temperature of the atoms (which determines the excitation of the hyperfine levels), x_{\mathrm{HI}} is the neutral hydrogen fraction, and \delta_b is the baryonic overdensity. The \tau of the transition is small (\tau \ll 1), typically \tau \approx 0.01 x_{\mathrm{HI}} (1 + \delta_b) (1+z)^{3/2} / T_s at z ≈ 10, but it contributes a correction term -\tau in more precise expressions for δT_b when foregrounds or high densities are considered. The spin temperature T_s is set by couplings to the CMB, the kinetic temperature of the gas T_k (via collisions), and the color temperature of the Ly-α background (via Wouthuysen-Field effect), allowing δT_b to trace thermal evolution in the IGM. The 21-cm signal manifests in two regimes relative to the CMB: absorption when T_s < T_{\mathrm{CMB}} (yielding negative δT_b) and emission when T_s > T_{\mathrm{CMB}} (positive δT_b). During the cosmic dark ages (z ≳ 30), collisional decoupling keeps T_s \approx T_{\mathrm{CMB}}, producing negligible signal. At cosmic dawn (z ≈ 15–30), the first UV photons from Population III stars couple T_s to the cooling gas kinetic temperature (T_k < T_{\mathrm{CMB}} due to adiabatic expansion), resulting in global absorption signals with δT_b ≈ -100 to -500 . Subsequent X-ray heating from early sources raises T_k > T_{\mathrm{CMB}}, shifting to emission with δT_b > 0 . During reionization (z ≈ 6–15), the declining x_{\mathrm{HI}} suppresses the signal amplitude, but spatial fluctuations in density, ionization, and temperature enable 21-cm tomography—three-dimensional mapping of neutral hydrogen distribution to reveal ionized bubbles carved by UV sources. These fluctuations, quantified via the power spectrum P(k), show scale-dependent features: large-scale modes (k ≈ 0.1 h Mpc⁻¹) trace bubble morphology, while small scales probe minihalos and density peaks. Key experiments target both the global (sky-averaged) signal and statistical fluctuations. The Experiment to Detect the Global Epoch of Reionization (EDGES) reported in 2018 the first putative detection of an trough centered at 78 MHz ( ≈ 17), with amplitude ≈ -500 mK and duration ≈ 80 MHz in frequency (Δ ≈ 4), far deeper than standard predictions of -200 mK and suggesting enhanced cooling or radio backgrounds. This result remains debated due to potential foreground systematics and foreground removal challenges, with subsequent analyses questioning its astrophysical origin. For power spectrum measurements, the Epoch of Reionization Array () uses a 350-antenna interferometer in to target = 7–12, achieving upper limits such as Δ²_{21} ≤ 457 mK² (95% confidence level) at k ≈ 0.34 h Mpc^{-1} and ≈ 7.9 using 94 nights of Phase I data (as of 2023), with analyses of subsequent seasons, including Year 6 in 2025, providing deeper constraints and progressing toward detection of reionization-era fluctuations. The upcoming (SKA) Low will extend sensitivity to = 6–15, enabling high-fidelity imaging of ionized bubbles. Forecasts indicate that 21-cm with SKA-Low can detect signatures of reionization —regions of suppressed signal corresponding to ionized volumes of 10–100 Mpc —at = 6–12, distinguishing between inside-out (galaxy-driven) and outside-in (quasar-dominated) reionization scenarios with signal-to-noise ratios exceeding 10 after 1000 hours of observation.

Theoretical Models and Recent Insights

Semi-Analytic and Radiative Transfer Models

Semi-analytic models provide an efficient framework for simulating the large-scale evolution of reionization by combining halo functions from extended Press-Schechter theory with prescriptions for ionizing photon production and escape. These models track the volume-averaged ionized fraction x_e over time, incorporating parameters such as the escape fraction f_\mathrm{esc} of ionizing photons from early galaxies and Population III stars, which determine the overall photon budget. The extended Press-Schechter formalism is used to derive the distribution of ionized sizes, treating reionization as a where halos above a minimum host sources that expand H II regions until overlap occurs. A key application of these models is in semi-numeric simulations like 21cmFAST, which approximate the field by assigning to grid cells based on excursion and propagate them with a one-step filtering approach, enabling rapid exploration of parameter space. These methods evolve x_e by balancing production rates against recombinations, often assuming a uniform intergalactic medium density, and predict the transition from neutral to ionized phases around z \sim 6-10. By linking collapse fractions to source properties, such models highlight how variations in f_\mathrm{esc} influence the timing and patchiness of reionization. Radiative transfer (RT) simulations offer more detailed treatments by solving the equations of photon propagation, absorption, and scattering in three dimensions, typically post-processing N-body hydrodynamic simulations of structure formation. Codes such as Enzo, extended with ray-tracing modules like Moray, perform adaptive mesh refinement to resolve small-scale gas dynamics while capturing large ionized bubbles, coupling RT to hydrodynamics for self-consistent evolution of ionization fronts. Similarly, RAMSES-RT integrates multi-frequency RT into the adaptive mesh refinement framework of RAMSES, handling photon transport via a moment-based method that accounts for Doppler shifts and radiative feedback on galaxy formation. These RT approaches simulate the inhomogeneous propagation of ionizing radiation from discrete sources, revealing effects like shadowing and bubble merging that semi-analytic models approximate. Both semi-analytic and RT models generate predictions for key observables, including the topology of ionized regions, the 21-cm power spectrum from neutral hydrogen fluctuations, and the Thomson optical depth \tau_e to the . Reionization is parameterized by the efficiency \zeta, which quantifies the number of ionizing photons per collapsed atom in halos, and the R_\mathrm{mfp} of photons through the intergalactic medium, which sets the scale for recombination-limited bubble growth. For instance, typical values of \zeta \approx 10-100 and R_\mathrm{mfp} \approx 10-50 Mpc yield \tau_e \sim 0.05-0.09, consistent with Planck measurements, while the 21-cm power spectrum exhibits a characteristic "bump" from large-scale ionization modes during the mid-to-late stages. These models collectively predict a patchy reionization process, where ionized bubbles grow around overdensities and reach characteristic sizes of approximately 10–100 comoving Mpc by the end of the , reflecting the hierarchical buildup of and photon escape efficiencies.

JWST Observations and Model Tensions

Recent observations from the (JWST), particularly through the UNCOVER and CEERS surveys, have revealed an unexpectedly high abundance of faint, low-mass galaxies at redshifts z > 10, suggesting these systems played a pivotal role in driving cosmic reionization. The UNCOVER survey, targeting the cluster, has identified dozens of such galaxies with vigorous , exhibiting UV luminosity densities up to four times higher than pre-JWST predictions from semi-analytic models. Similarly, CEERS data confirm a surplus of UV-bright galaxies at z ≈ 9–11, with spectroscopic follow-up indicating these faint sources dominate the early ionizing photon budget. These findings create significant tensions with established reionization models, particularly regarding the timing and sources of . The abundance of early galaxies implies a reionization z_reion > 10, conflicting with (CMB) measurements of the τ_e ≈ 0.058 ± 0.006 from Planck PR4 (2025), which favor a later reionization ending around z ≈ 6–8. This discrepancy, dubbed the "JWST reionization crisis," suggests overproduction of ionizing photons in standard ΛCDM simulations unless adjusted for higher Lyman-continuum escape fractions (f_esc > 0.2–0.5) in these compact systems. Debates persist on f_esc , with JWST spectra showing UV continua consistent with near-unity escape in some z > 10 galaxies, challenging assumptions of low escape in dense early environments. As of 2025, Planck PR4 data refine τ_e constraints, while JWST surveys like UNCOVER and CEERS continue to reveal high galaxy abundances at z > 10, pushing simulations to incorporate bursty and higher efficiencies to address ongoing tensions. The implications extend to revising source efficiencies and simulation frameworks. Small, faint galaxies appear to contribute over 50% of the ionizing photons required for reionization, as inferred from their collective UV output and high f_esc, thereby diminishing the expected dominance of quasars and active galactic nuclei. This necessitates updated models incorporating bursty histories (SFH) to account for episodic high-efficiency ionization from these dwarfs, rather than steady-state assumptions. Ongoing studies as of November 2025 emphasize the need for deeper JWST surveys to resolve these tensions and refine the photon budget.

References

  1. [1]
    [1907.06653] An Introductory Review on Cosmic Reionization - arXiv
    Jul 15, 2019 · This article gives an introduction to the theoretical and observational aspects of cosmic reionization and discusses their role in our understanding of early ...
  2. [2]
    Modeling cosmic reionization | Living Reviews in Computational ...
    Nov 28, 2022 · In this review, we survey the physics and key aspects of reionization-era modeling and describe the diverse range of computational techniques and tools ...
  3. [3]
    Witnessing the onset of reionization through Lyman-α emission at ...
    Mar 26, 2025 · Cosmic reionization began when ultraviolet (UV) radiation produced in the first galaxies began illuminating the cold, neutral gas that ...
  4. [4]
    Reionization and its sources
    ### Summary of Cosmic Reionization
  5. [5]
    None
    ### Summary of Cosmic Reionization from Cooray-Reionization.pdf
  6. [6]
    [1807.06209] Planck 2018 results. VI. Cosmological parameters - arXiv
    Jul 17, 2018 · Abstract:We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies.
  7. [7]
  8. [8]
    None
    Summary of each segment:
  9. [9]
    [PDF] Solution of the Saha equation from ionization fraction and ...
    In particular, we plan to apply this formalism to the precise computation of late CMB spectral distortions associated to the cosmological reionization process.
  10. [10]
    Heating of the Intergalactic Medium by Hydrogen Reionization - arXiv
    Jul 24, 2018 · During reionization, the intergalactic medium is heated impulsively by supersonic ionization fronts (I-fronts). The peak gas temperatures ...
  11. [11]
    Heating of the Intergalactic Medium by Hydrogen Reionization
    Apr 3, 2019 · During reionization, the intergalactic medium is heated impulsively by supersonic ionization fronts (I-fronts).
  12. [12]
    Suppression of star formation in low-mass galaxies caused by the ...
    ABSTRACT. Photoheating associated with reionization suppressed star formation in low-mass galaxies. Reionization was inhomogeneous, however, affecting diff.
  13. [13]
    Measurement of the baryonic acoustic oscillation scale in 21 cm ...
    However, the bias is expected to depend on scale during the process of reionization since galaxies preferentially form in overdense regions, resulting in early ...Missing: influence | Show results with:influence
  14. [14]
  15. [15]
    [1505.07568] The Cosmic Dawn and Epoch of Reionization ... - arXiv
    May 28, 2015 · Abstract page for arXiv paper 1505.07568: The Cosmic Dawn and Epoch of Reionization with the Square Kilometre Array.
  16. [16]
    Bridging the Gap between Cosmic Dawn and Reionization Favors ...
    Nov 28, 2023 · Bridging the gap between cosmic dawn and reionization favors models dominated by faint galaxies. Ankita Bera, Sultan Hassan, Aaron Smith, Renyue Cen, Enrico ...Introduction · Source Model · Source Model Calibration to... · Conclusions
  17. [17]
    Inside-out or outside-in: the topology of reionization in the photon ...
    Reionization proceeds then much more as an outside-in process. Low-density regions far from ionizing sources become ionized before regions of high gas density ...
  18. [18]
    [0711.1542] The History and Morphology of Helium Reionization
    Nov 9, 2007 · Abstract: A variety of observations now indicate that intergalactic helium was fully ionized by z~3. ... HeII reionization. Here we describe the ...
  19. [19]
    Building up the Population III initial mass function from cosmological ...
    Primordial star formation was instead driven by cooling through H2 line emission. Pop III stars are believed to have initially formed at z ∼ 20–30 in small ...Missing: H2 | Show results with:H2
  20. [20]
    The Birth Mass Function of Population III Stars - IOP Science
    Jan 21, 2022 · Population III stars are expected to form in 105–106 Me primordial minihalos at z = 20–30, as a result of gravitational collapse induced by H2.
  21. [21]
    [2010.02212] How the Population III Initial Mass Function Governs ...
    Oct 5, 2020 · We find that a top-heavy Pop III IMF results in earlier star formation but dimmer galaxies than a more conventional Salpeter-type IMF.
  22. [22]
    Pop III Stellar Masses and IMF - American Institute of Physics
    The mass scale for Pop III stars is fundamentally linked to the cooling properties of primordial gas that virializes to T ~ 1000 K in such halos. Trace amounts ...
  23. [23]
    [PDF] Generic spectrum and ionization efficiency of a heavy initial mass ...
    Since metal-free stars have systematically higher effective temperatures (Tef f ∼ 105 K for M > 100M⊙), they are very efficient at producing photons capable of ...
  24. [24]
    (PDF) Ultraviolet photon production rates of the first stars: Impact on ...
    Jul 30, 2025 · The very high effective temperatures ( ∼ 1 0 5 \sim 10^5 K) predicted for the most massive Population III stars could also give rise to ...
  25. [25]
    Evolution of massive Population III stars with rotation and magnetic ...
    For example, the main sequence lifetime (τH) and the helium core size (MHe − Core) of our non-rotating 200 M⊙ model are 2.40 Myr and 114 M⊙, respectively, (see ...
  26. [26]
    Pair Instability Supernovae of Very Massive Population III Stars
    Stellar models indicate that non-rotating Population III stars with initial masses of 140-260 Msun die as highly energetic pair-instability supernovae.
  27. [27]
    [PDF] Cosmic reionization by stellar sources: population III stars - NASA ADS
    In all cases, however, it appears that a large escape fraction is required during the population III epoch in order to contribute significantly to the ...
  28. [28]
    Cosmic reionization by stellar sources: population III stars
    Additionally, we generated synthetic spectra from a range of models with different escape fractions and showed that values in the range 10–20 per cent led to ...
  29. [29]
    Semi-analytical modelling of Pop. III star formation and metallicity ...
    We find Pop. III stars have little to no impact on the reionization history but significantly affect the thermal state of the intergalactic medium (IGM).
  30. [30]
    Baseline metal enrichment from Population III star formation in ...
    Metals left behind by the death of Pop III stars in supernova (SN) explosions enrich the early interstellar/intergalactic medium (ISM/IGM) (e.g. Yoshida, Bromm ...
  31. [31]
    THE SPECTRAL EVOLUTION OF THE FIRST GALAXIES. II ...
    The fraction of ionizing photons that escape (fesc) from z ≳ 6 galaxies is an important parameter for assessing the role of these objects in the reionization of ...Missing: output | Show results with:output
  32. [32]
    The Effects of Radiative Feedback and Supernova-induced ...
    At the same time, radiative transfer suppresses Population III star formation, via the development of ionized bubbles that slow gas accretion in these regions, ...
  33. [33]
    High-redshift Galaxy Formation with Self-consistently Modeled Stars ...
    We find that feedback from SFMCs and an accreting MBH suppresses runaway star formation locally in the galactic core region. Newly included radiation feedback ...
  34. [34]
    JWST Unveils Obscured Quasars in the Most Luminous Galaxies at ...
    Jul 14, 2025 · Bridging the Gap: JWST Unveils Obscured Quasars in the Most Luminous Galaxies at z > 6. Yoshiki Matsuoka, Masafusa Onoue, Kazushi Iwasawa, ...
  35. [35]
    Cosmic reionization after Planck II: contribution from quasars
    The fraction of photoionization rate contributed by the quasars, as shown in Fig. 3, is significantly small. It can take values at most ∼ 10 per cent ...Missing: percentage | Show results with:percentage
  36. [36]
    RELATIVE ROLE OF STARS AND QUASARS IN COSMIC ...
    ABSTRACT. We revisit the classical view that quasar contribution to the reionization of hydrogen is unimportant. When secondary ionizations are taken into ...
  37. [37]
    [PDF] Definitive upper bound on the negligible contribution of quasars to ...
    Quasar contribution to reionization is negligible, providing less than 7% of the total photons needed to keep the universe ionized at z=6.0-6.6.
  38. [38]
    On the Sizes of Ionized Bubbles Around Galaxies During the ...
    Aug 29, 2023 · The simulated bubble sizes then increase toward z ≃ 7, forming a fully ionized universe (effectively infinite bubble radii) by z = 6.5.
  39. [39]
    Local positive feedback in the overall negative: the impact of quasar ...
    In conclusion, our results suggest that positive and negative AGN feedback coexist in galaxies, where strong quasar winds can suppress global stellar mass ...
  40. [40]
    Discovery and characterization of 25 new quasars at 4.6 &lt; z &lt
    Jul 23, 2025 · We report the discovery of 25 new quasars at 4.6<z<6.9 (six at z≥6.5), with M1450 between −25.4 and −27.0. We also present new spectra of six z> ...
  41. [41]
  42. [42]
  43. [43]
    Constraining the Evolution of the Ionizing Background ... - IOP Science
    Constraining the evolution of the ionizing background and the epoch of reionization with z ∼ 6 quasars. II. A sample of 19 quasars.
  44. [44]
    IGM opacity constraints from the Lyman α forest of galaxies in legacy ...
    We present the first characterization of the Gunn–Peterson trough in high-redshift galaxies using public James Webb Space Telescope (JWST) NIRSpec spectroscopy.Missing: 2024 | Show results with:2024
  45. [45]
    XLVII. Planck constraints on reionization history
    The CMB is affected by the total column density of free electrons along each line of sight, parameterized by its Thomson scattering optical depth τ. This is one ...
  46. [46]
    Planck 2018 results. VI. Cosmological parameters - ADS
    Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision ...
  47. [47]
    Do JWST reionization (optical depth) puzzle, cosmological tensions ...
    Mar 10, 2025 · The James Webb Space Telescope (JWST) has observed massive galaxies at high redshifts, which implies an earlier epoch of reionization (EoR) ...
  48. [48]
  49. [49]
    SILVERRUSH. XIV. Lyα Luminosity Functions and Angular ...
    Mar 18, 2025 · Abstract. We present luminosity functions (LFs) and angular correlation functions (ACFs) derived from 18,960 Lyα emitters (LAEs) at z = 2.2− ...
  50. [50]
    The characteristic shape of damping wings during reionization
    ABSTRACT. Spectroscopic analysis of Ly α damping wings of bright sources at z > 6 is a promising way to measure the reionization history of the universe.
  51. [51]
    [2405.12273] Damping wings in the Lyman-α forest: a model ... - arXiv
    May 20, 2024 · Abstract page for arXiv paper 2405.12273: Damping wings in the Lyman-α forest: a model-independent measurement of the neutral fraction at 5.4.
  52. [52]
    [2210.06504] Lyman-alpha Emitters in Ionized Bubbles - arXiv
    Oct 12, 2022 · Lyman-alpha emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible.Missing: coverage f_cov
  53. [53]
    Insights into the Lyman-alpha line widths and the size of ionized ...
    Oct 23, 2024 · We find a positive correlation between bubble radius and Lya line width, which again suggests that large ionized bubbles are created around these LAEs.<|control11|><|separator|>
  54. [54]
    The 21 cm Transition and the High-Redshift Universe - astro-ph - arXiv
    Aug 1, 2006 · Here we review the physics of the 21 cm transition, focusing on processes relevant at high redshifts, and describe the insights to be gained ...
  55. [55]
    An absorption profile centred at 78 megahertz in the sky-averaged ...
    Mar 1, 2018 · The 21-cm absorption profile is detected in the sky-averaged radio spectrum, but is much stronger than predicted, suggesting that the ...
  56. [56]
    The EDGES measurement disfavors an excess radio background ...
    In 2018 the EDGES experiment claimed the first detection of the global cosmic 21 cm signal, which featured an absorption trough centered around z ∼ 17 with ...
  57. [57]
    Upper Limits on the Epoch of Reionization 21 cm Power Spectrum
    This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms ...
  58. [58]
    Radiative transfer simulations of cosmic reionization I - astro-ph - arXiv
    Dec 15, 2006 · We present a new hybrid code for large volume, high resolution simulations of cosmic reionization, which utilizes a N-body algorithm for dark matter.Missing: seminal Enzo RT
  59. [59]
    RAMSES-RT: Radiation hydrodynamics in the cosmological context
    Apr 26, 2013 · We present a new implementation of radiation hydrodynamics (RHD) in the adaptive mesh refinement (AMR) code RAMSES.Missing: simulations reionization Enzo
  60. [60]
  61. [61]
    First Release of Ultradeep JWST/NIRSpec PRISM spectra for ~700 ...
    Aug 7, 2024 · The UNCOVER Survey: First Release of Ultradeep JWST/NIRSpec PRISM spectra for ~700 galaxies from z~0.3-13 in Abell 2744. Authors:Sedona H. Price ...
  62. [62]
    NASA's Webb 'UNCOVERs' Galaxy Population Driving Cosmic ...
    Jun 11, 2025 · Astronomers using data from NASA's James Webb Space Telescope have identified dozens of small galaxies that played a starring role in a cosmic makeover.
  63. [63]
    JWST confirms a high abundance of galaxies and AGN at z ≃ 9–11 ...
    We present the JWST/NIRSpec PRISM follow-up of candidate galaxies at z ≃ 9–11 selected from deep JWST/NIRCam photometry in GLASS-JWST Early Release Science data ...
  64. [64]
    A Census from JWST of Extreme Emission-line Galaxies Spanning ...
    Oct 3, 2024 · We present a sample of 1165 extreme emission-line galaxies (EELGs) at 4 < z < 9 selected using James Webb Space Telescope (JWST) NIRCam photometry.
  65. [65]
    The cosmic timeline implied by the JWST reionization crisis
    Introduction. Observational constraints indicate that the epoch of reionization (EoR) probably began around redshift z ∼ 10 − 15 and ended by z ∼ 6 (McGreer ...
  66. [66]
    Redshift evolution of Lyman continuum escape fraction after JWST
    As the cosmic sSFR increases with redshift, so does f_{\rm esc}(z), which, when globally averaged, grows from 0.007 to 0.6 in 0 < z < 20.Missing: reionization 2024
  67. [67]
    Webb finds dwarf galaxies reionised the Universe
    Feb 28, 2024 · Webb found that dwarf galaxies are likely the main producers of radiation that reionized the early universe, with their radiation levels four ...
  68. [68]
    Inferring the ionizing photon contributions of high-redshift galaxies to ...
    ABSTRACT. JWST observations are providing unprecedented constraints on the history of reionization owing to the ability to detect faint galaxies at ...