Emulating the global 21-cm signal from Cosmic Dawn and Reionization
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Reionization is an inhomogeneous process, thought to begin in small ionized bubbles of the intergalactic medium (IGM) around overdense regions of galaxies. Recent Lyα studies during the epoch of reionization show evidence that ionized bubbles formed earlier around brighter galaxies, suggesting higher IGM transmission of Lyα from these galaxies. We investigate this problem using IR slitless spectroscopy from the Hubble Space Telescope (HST) Wide-Field Camera 3 (WFC3) G102 grism observations of 148 galaxies selected via photometric redshifts at 6.0 < z < 8.2. These galaxies have spectra extracted from the CANDELS Lyα Emission at Reionization (CLEAR) survey. We combine the CLEAR data for 275 galaxies with the Keck Deep Imaging Multi-Object Spectrograph and MOSFIRE data set from the Texas Spectroscopic Search for Lyα Emission at the End of Reionization Survey. We constrain the Lyα equivalent width (EW) distribution at 6.0 < z < 8.2, which is described by an exponential form, dN / dEW ∝ exp(-EW) W(sub 0), with the characteristic e-folding scale width (W(sub 0)). We confirm a significant drop in the Lyα strength (i.e., W(sub 0)) at z > 6. Furthermore, we compare the redshift evolution of W(sub 0) between galaxies at different UV luminosities. UV-bright (M(sub UV) < −21 [i.e., L(sub UV) > L*]) galaxies show weaker evolution with a decrease of 0.4 ( ± 0.2) dex in W(sub 0) at z > 6, while UV-faint (M(sub UV) > −21 [LUV < L*]) galaxies exhibit a significant drop of 0.7–0.8 (±0.2) dex in W(sub 0) from z < 6 to z > 6. If the change in W(sub 0) is proportional to the change in the IGM transmission for Lyα photons, then this is evidence that the transmission is “boosted” around UV-brighter galaxies, suggesting that reionization proceeds faster in regions around such galaxies.
Reionization is an inhomogeneous process, thought to begin in small ionized bubbles of the intergalacticmedium (IGM) around overdense regions of galaxies. Recent Lyman-alpha (Lyα) studies during the epochof reionization show growing evidence that ionized bubbles formed earlier around brighter galaxies, suggest-ing higher IGM transmission of Lyαfrom these galaxies. We investigate this problem using IR slitless spec-troscopy from the Hubble Space Telescope (HST) Wide-Field Camera 3 (WFC3) G102 grism observations of148 galaxies selected via photometric redshifts at6.0< z <8.2. These galaxies have spectra extracted fromthe CANDELS LyαEmission at Reionization (CLEAR) survey. We combine the CLEAR data for 275 galaxieswith the Keck/DEIMOS+MOSFIRE dataset from the Texas Spectroscopic Search for LyαEmission at the Endof Reionization Survey. We then constrain the Lyαequivalent-width (EW) distribution at6.0< z <8.2, whichis described by an exponential form,dN/dEW∝exp(-EW)/W0, with the characteristice-folding scale width(W0). We confirm a significant drop of the Lyαstrength (orW0) atz >6. Furthermore, we compare the redshiftevolution ofW0between galaxies at different UV luminosities. The UV-bright (MUV<−21, orLUV> L∗)galaxies show weaker evolution with a decrease of 0.4 (±0.2) dex inW0atz >6while UV-faint (MUV>−21,orLUV< L∗) galaxies exhibit a significant drop by a factor of 0.7-0.8 (±0.2) dex inW0fromz <6toz >6.Our results add to the accumulating evidence that UV-bright galaxies exhibit boosted Lyαtransmission in theIGM, suggesting that reionization completes sooner in regions proximate to galaxies of higher UV luminosity.
Patchy reionization generates kinematic Sunyaev–Zel’dovich (kSZ) anisotropies in the cosmic microwave background (CMB). Large-scale velocity perturbations along the line of sight modulate the small-scale kSZ power spectrum, leading to a trispectrum (or four-point function) in the CMB that depends on the physics of reionization. We investigate the challenges in detecting this trispectrum and use tools developed for CMB lensing, such as realization-dependent bias subtraction and cross-correlation based estimators, to counter uncertainties in the instrumental noise and assumed CMB power spectrum. We also find that both lensing and extragalactic foregrounds can impart larger trispectrum contributions than the reionization kSZ signal. We present a range of mitigation methods for both of these sources of contamination, validated on microwave-sky simulations. We use ACT DR6 and Planck data to calculate an upper limit on the reionization kSZ trispectrum from a measurement dominated by foregrounds. The upper limit is about 50 times the signal predicted from recent simulations.
Abstract Dense gas in minihalos with masses of 10 6 −10 8 M ⊙ can shield themselves from reionization for ∼100 Myr after being exposed to the UV background. These self-shielded systems, often unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Ly α forest. We first derive the H i column density profile of photoevaporating minihalos by conducting 1D radiation–hydrodynamics simulations. We utilize these results to estimate the Ly α opacity from minihalos in a large-scale simulation that cannot resolve self-shielding. When the ionization rate of the background radiation is 0.03 × 10 −12 s −1 , as expected near the end of reionization at z ∼ 5.5, we find that the incidence rate of damped Ly α absorbers increases by a factor of ∼2−4 compared to at z = 4.5. The Ly α flux is, on average, suppressed by ∼3% of its mean due to minihalos. The absorption features enhance the 1D power spectrum up to ∼5% at k ∼ 0.1 h Mpc −1 (or 10 −3 km −1 s), which is comparable to the enhancement caused by inhomogeneous reionization. The flux is particularly suppressed in the vicinity of large halos along the line-of-sight direction at separations of up to 10 h −1 Mpc at r ⊥ ≲ 2 h −1 Mpc. However, these effects become much smaller for higher ionizing rates (≳0.3 × 10 −12 s −1 ) expected in the post-reionization Universe. Our findings highlight the need to consider minihalo absorption when interpreting the Ly α forest at z ≳ 5.5. Moreover, the sensitivity of these quantities to the ionizing background intensity can be exploited to constrain the intensity itself.
We investigate the properties of cosmological ionization fronts during the epoch of reionization using the Cosmic Reionization on Computers simulations. By analyzing reionization timing maps, we characterize ionization front velocities and curvatures and their dependence on the density structure of the intergalactic medium (IGM). The velocity distribution of ionization fronts in the simulations indicates that while the barrier-crossing analytical model captures the overall shape in high-velocity regions, it fails to reproduce the low-velocity tail, highlighting the non-Gaussian nature of the IGM’s density field. Ionization front velocities are inversely correlated with local density, propagating faster in underdense regions and more slowly in overdense environments. Faster ionization fronts also lead to higher postionization temperatures, reaching a plateau at ∼2 × 10$^{4}$ K for velocities exceeding 3000 km s$^{−1}$. Examining curvature statistics further establishes a connection between ionization front structure and the normalized density contrast ν, with trends in overdense regions aligning well with barrier-crossing model predictions, while deviations appear in underdense environments due to model limitations. These results provide a detailed characterization of ionization front dynamics and their interaction with the underlying density field, bridging small-scale reionization physics with large-scale observables such as the 21 cm signal and the IGM’s thermal history.
We present full sky microwave maps in five frequency bands (23 to 94 GHz) from the WMAP first year sky survey. Calibration errors are less than 0.5% and the low systematic error level is well specified. The cosmic microwave background (CMB) is separated from the foregrounds using multifrequency data. The sky maps are consistent with the 7 in. full-width at half-maximum (FWHM) Cosmic Background Explorer (COBE) maps. We report more precise, but consistent, dipole and quadrupole values. The CMB anisotropy obeys Gaussian statistics with -58 less than f(sub NL) less than 134 (95% CL). The 2 less than or = l less than or = 900 anisotropy power spectrum is cosmic variance limited for l less than 354 with a signal-to-noise ratio greater than 1 per mode to l = 658. The temperature-polarization cross-power spectrum reveals both acoustic features and a large angle correlation from reionization. The optical depth of reionization is tau = 0.17 +/- 0.04, which implies a reionization epoch of t(sub r) = 180(sup +220, sub -80) Myr (95% CL) after the Big Bang at a redshift of z(sub r) = 20(sup +10, sub -9) (95% CL) for a range of ionization scenarios. This early reionization is incompatible with the presence of a significant warm dark matter density. A best-fit cosmological model to the CMB and other measures of large scale structure works remarkably well with only a few parameters. The age of the best-fit universe is t(sub 0) = 13.7 +/- 0.2 Gyr old. Decoupling was t(sub dec) = 379(sup +8, sub -7)kyr after the Big Bang at a redshift of z(sub dec) = 1089 +/- 1. The thickness of the decoupling surface was Delta(sub z(sub dec)) = 195 +/- 2. The matter density of the universe is Omega(sub m)h(sup 2) = 0.135(sup +0.008, sub -0.009) the baryon density is Omega(sub b)h(sup 2) = 0.0224 +/- 0.0009, and the total mass-energy of the universe is Omega(sub tot) = 1.02 +/- 0.02. There is progressively less fluctuation power on smaller scales, from WMAP to fine scale CMB measurements to galaxies and finally to the Ly-alpha forest. This is accounted for with a running spectral index, significant at the approx. 2(sigma) level. The spectral index of scalar fluctuations is fit as n(sub s) = 0.93 +/-0.03 at wavenumber k(sub o) = 0.05/Mpc ((sub eff) approx. = 700), with a slope of dn(sub s)/d I(sub nk) = -0.031(sup + 0.016, sub -0.018) in the best-fit model.
The Square Kilometre Array (SKA) is expected to provide the first tomographic observations of the neutral intergalactic medium at redshifts z > 6 and pinpoint the locations of individual ionized bubbles during early stages of cosmic reionization. In scenarios where star-forming galaxies provide most of the ionizing photons required for cosmic reionization, one expects the first ionized bubbles to be centred on overdensities of such galaxies. Here, we model the properties of galaxy populations within isolated, ionized bubbles that SKA-1 should be able to resolve at z ≈ 7–10, and explore the prospects for galaxy counts within such structures with various upcoming near-infrared telescopes. We find that, for the bubbles that are within reach of SKA-1 tomography, the bubble volume is closely tied to the number of ionizing photons that have escaped from the galaxies within. In the case of galaxy-dominated reionization, galaxies are expected to turn up above the spectroscopic detection threshold of JWST and ELT in even the smallest resolvable bubbles at z ≤ 10. The prospects of detecting galaxies within these structures in purely photometric surveys with Euclid, WFIRST, JWST, or ELT are also discussed. While spectroscopy is preferable towards the end of reionization to provide a robust sample of bubble members, multiband imaging may be a competitive option for bubbles at z ≈ 10, due to the very small number of line-of-sight interlopers expected at these redshifts.
The Cosmology Large Angular Scale Surveyor (CLASS) is a four telescope array designed to characterize relic primordial gravitational waves from in ation and the optical depth to reionization through a measurement of the polarized cosmic microwave background (CMB) on the largest angular scales. The frequencies of the four CLASS telescopes, one at 38 GHz, two at 93 GHz, and one dichroic system at 145/217 GHz, are chosen to avoid spectral regions of high atmospheric emission and span the minimum of the polarized Galactic foregrounds: synchrotron emission at lower frequencies and dust emission at higher frequencies. Low-noise transition edge sensor detectors and a rapid front-end polarization modulator provide a unique combination of high sensitivity, stability, and control of systematics. The CLASS site, at 5200 m in the Chilean Atacama desert, allows for daily mapping of up to 70% of the sky and enables the characterization of CMB polarization at the largest angular scales. Using this combination of a broad frequency range, large sky coverage, control over systematics, and high sensitivity, CLASS will observe the reionization and recombination peaks of the CMB E- and B-mode power spectra. CLASS will make a cosmic variance limited measurement of the optical depth to reionization and will measure or place upper limits on the tensor-to-scalar ratio, r, down to a level of 0.01 (95% C.L.).
ABSTRACT The autocorrelation function of the Lyman-$\alpha$ (Ly $\alpha$) forest flux from high-z quasars probes the small-scale structure of the intergalactic medium (IGM). The thermal state of the IGM, determined by the physics of reionization, sets the small-scale power observed in the Ly $\alpha$ forest. To explore the sensitivity of the autocorrelation function to the IGM’s thermal state, we compute the autocorrelation function from a cosmological hydrodynamical simulation with an instantaneous reionization model and 135 post-processed thermal states. Using mock data sets of 20 quasars, we forecast constraints on $T_0$ and $\gamma$, which characterize the post-processed IGM thermal state, at $5.4 \le z \le 6$. While this model simplifies the IGM’s thermal state, it serves as a key first step in assessing future observational prospects. We also perform an inference test on mocks and re-weight out posterior distributions to guarantee that they exhibit statistically correct behaviour. At $z = 5.4$, we find that an idealized data set constrains $T_0$ to 59 per cent and $\gamma$ to 16 per cent at the 1$\sigma$ equivalent confidence level. To explore more realistic, non-instantaneous reionization scenarios, we analyse four models combining temperature and ultraviolet background (UVB) fluctuations at $z = 5.8$. We find that mock data generated from a model with both temperature and UVB fluctuations can rule out a model with only temperature fluctuations at the $> 1\sigma$ level 73.9 per cent of the time.
Understanding dwarf galaxy formation is crucial for testing dark matter models and reionization physics. However, constructing stellar-mass complete spectroscopic samples at low masses is increasingly difficult, and the potential existence of a local void complicates studies in an average environment. The Photometric object Around Cosmic webs (PAC) method, which combines deep photometric and spectroscopic data to measure the excess surface density $\bar{n}_2w_{\rm{p}}(r_{\rm{p}})$ of photometric objects around spectroscopic tracers, offers a promising path forward. We model 349 $\bar{n}_2w_{\rm{p}}(r_{\rm{p}})$ measurements from DESI Y1 BGS and DECaLS, reaching $M_*=10^{6.4}\,{\rm M}_{\odot}$, using a stellar mass-halo mass relation (SHMR)-based subhalo abundance matching framework applied to two high-resolution $N$-body simulations from the Jiutian suite. The resulting SHMR is constrained down to $M_{\rm h}\simeq10^{8.0}\,h^{-1}{\rm M}_{\odot}$, revealing a clear upturn at $\sim10^{10.0}\,h^{-1}{\rm M}_{\odot}$ toward lower masses, indicating rising star-formation efficiency (SFE) in small haloes. This feature persists under extensions of the model that allow mass-dependent scatter, reionization-induced suppression of the halo occupation fraction, galaxy assembly bias, and alternative cosmologies. Combining with the results from Paper I, we find that central red galaxies dominate the low-mass regime. Our results motivate a hypothesis in which SFE is significantly higher than previously thought prior to reionization, enabling relatively massive galaxies to form in small haloes. These systems are subsequently quenched by the UV background, producing the central red dwarf galaxies observed. Finally, we obtain $3σ$ and $5σ$ upper mass bounds of $10^{8.80}\,h^{-1}{\rm M}_{\odot}$ and $10^{10.24}\,h^{-1}{\rm M}_{\odot}$ on the smallest haloes required to exist.
The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are inconsistent with instrument noise and are significantly larger than the upper limits established for potential systematic errors. The correlations are present in all WAMP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations and fitting techniques. On small angular scales (theta less than 5 deg), the WMAP data show the temperature-polarization correlation expected from adiabatic perturbations in the temperature power spectrum. The data for l greater than 20 agree well with the signal predicted solely from the temperature power spectra, with no additional free parameters. We detect excess power on large angular scales (theta greater than 10 deg) compared to predictions based on the temperature power spectra alone. The excess power is well described by reionization at redshift 11 is less than z(sub r) is less than 30 at 95% confidence, depending on the ionization history. A model-independent fit to reionization optical depth yields results consistent with the best-fit ACDM model, with best fit value t = 0.17 +/- 0.04 at 68% confidence, including systematic and foreground uncertainties. This value is larger than expected given the detection of a Gunn-Peterson trough in the absorption spectra of distant quasars, and implies that the universe has a complex ionization history: WMAP has detected the signal from an early epoch of reionization.
The science objectives of the James Webb Space Telescope fall into four themes. The End of the Dark Ages: First Light and Reionization theme seeks to identify the first luminous sources to form and to determine the ionization history of the universe. The Assembly of Galaxies theme seeks to determine how galaxies and the dark matter, gas, stars, metals, morphological structures, and black holes within them evolved from the epoch of reionization to the present. The Birth of Stars and Protoplanetary Systems theme seeks to unravel the birth and early evolution of stars, from infall onto dust-enshrouded protostars, to the genesis of planetary systems. The Planetary Systems and the Origins of Life theme seeks to determine the physical and chemical properties of planetary systems around nearby stars and of our own, and to investigate the potential for life in those systems. These four science themes were used to establish the design requirements for the observatory and instrumentation. Since Webb's capabilities are unique, those science themes will remain relevant through launch and operations and goals contained within these themes will continue to guide the design and implementation choices for the mission. More recently, it has also become clear that Webb will make major contributions to other areas of research, including dark energy, dark matter, active galactic nuclei, stellar populations, exoplanet characterization and Solar System objects. In this paper, we review the original four science themes and discuss how the scientific output of Webb will extend to these new areas of research. The James Webb Space Telescope was designed to meet science objectives in four themes: The End of the Dark Ages: First Light and Reionization, The Assembly of Galaxies, The Birth of Stars and Protoplanetary Systems, and Planetary Systems and the Origins of Life. More recently, it has become clear that Webb will also make major contributions to studies of dark energy, dark matter, active galactic nuclei, stellar populations, exoplanet characterization and Solar System objects. We review the original four science themes and discuss how the scientific output of Webb will extend to these new areas of research.
The Cosmology Large Angular Scale Surveyor (CLASS) is a four telescope array designed to characterize relic primordial gravitational waves from inflation and the optical depth to reionization through a measurement of the polarized cosmic microwave background (CMB) on the largest angular scales. The frequencies of the four CLASS telescopes, one at 38 GHz, two at 93 GHz, and one dichroic system at 145217 GHz, are chosen to avoid spectral regions of high atmospheric emission and span the minimum of the polarized Galactic foregrounds: synchrotron emission at lower frequencies and dust emission at higher frequencies. Low-noise transition edge sensor detectors and a rapid front-end polarization modulator provide a unique combination of high sensitivity, stability, and control of systematics. The CLASS site, at 5200 m in the Chilean Atacama desert, allows for daily mapping of up to 70% of the sky and enables the characterization of CMB polarization at the largest angular scales. Using this combination of a broad frequency range, large sky coverage, control over systematics, and high sensitivity, CLASS will observe the reionization and recombination peaks of the CMB E- and B-mode power spectra. CLASS will make a cosmic variance limited measurement of the optical depth to reionization and will measure or place upper limits on the tensor-to-scalar ratio, r, down to a level of 0.01 (95% C.L.).
We present measurements of large-scale cosmic microwave background E-mode polarization from the Cosmology Large Angular Scale Surveyor 90 GHz data. Using 115 det-yr of observations collected through 2024 with a variable-delay polarization modulator, we achieved a polarization sensitivity of 82 μK arcimin, comparable to Planck at similar frequencies (100 and 143 GHz ). The analysis demonstrates effective mitigation of systematic errors and addresses challenges to large-angular-scale power recovery posed by time-domain filtering in maximum-likelihood map-making. A novel implementation of the pixel-space transfer matrix is introduced, which enables efficient filtering simulations and bias correction in the power spectrum using the quadratic cross-spectrum estimator. Overall, we achieved an unbiased time-domain filtering correction to recover the largest angular scale polarization, with the only power deficit, arising from map-making nonlinearity, being characterized as <3%. Through cross-correlation with Planck, we detected the cosmic reionization at 99.4% significance and measured the reionization optical depth τ = $0.053^{+0.018}_{-0.019}$, marking the first ground-based attempt at such a measurement. At intermediate angular scales (ℓ > 30), our results, both independently and in cross-correlation with Planck, remain fully consistent with Planck’s measurements.
X-rays from high-mass X-ray binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; z ∼10–15) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift observations, as well as theoretical models, imply that the integrated X-ray luminosity to star formation rate of HMXBs (L_X/SFR) should increase in metal-poor environments, typical of early galaxies. Here, we study the impact of the metallicity (Z) dependence of L_X/SFR during the EoH. For our fiducial models, galaxies with star formation rates of order 10^(−3)–10^(−1)M_⊙ yr^(−1) and metallicities of order 10^(−3) –10^(−2) Z_⊙ are the dominant contributors to the X-ray background (XRB) during this period. Different L_X/SFR–Z relations result in factors of ∼ 3 differences in these ranges, as well as in the mean IGM temperature and the large-scale 21-cm power, at a given redshift. We compute mock 21-cm observations adopting as a baseline a 1000 h integration with the upcoming Square Kilometer Array (SKA) for two different L_X/SFR–Z relations. We perform inference on these mock observations using the common simplification of a constant L_X/SFR, finding that constant L_X/SFR models can recover the IGM evolution of the more complicated L_X/SFR–Z simulations only during the EoH. At z < 10, where the typical galaxies are more polluted, constant L_X/SFR models overpredict the XRB and its relative contribution to the early stages of the reionization.
We discuss the model of astrophysical emission at millimeter wavelengths used to characterize foregrounds in the multi-frequency power spectra of the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6), expanding on Louis et al. (2025) (2503.14452). We detail several tests to validate the capability of the DR6 parametric foreground model to describe current observations and complex simulations, and show that cosmological parameter constraints are robust against model extensions and variations. We demonstrate consistency of the model with pre-DR6 ACT data and observations from Planck and the South Pole Telescope. We evaluate the implications of using different foreground templates and extending the model with new components and/or free parameters. In all scenarios, the DR6 ΛCDM and ΛCDM+N eff cosmological parameters shift by less than 0.5σ relative to the baseline constraints. Some foreground parameters shift more; we estimate their systematic uncertainties associated with modeling choices. From our constraint on the kinematic Sunyaev-Zel'dovich power, we obtain a conservative limit on the duration of reionization of Δz rei < 4.4, assuming a reionization midpoint consistent with optical depth measurements and a minimal low-redshift contribution, with varying assumptions for this component leading to tighter limits. Finally, we analyze realistic non-Gaussian, correlated microwave sky simulations containing Galactic and extragalactic foreground fields, built independently of the DR6 parametric foreground model. Processing these simulations through the DR6 power spectrum and likelihood pipeline, we recover the input cosmological parameters of the underlying cosmic microwave background field, a new demonstration for small-scale CMB analysis. These tests validate the robustness of the ACT DR6 foreground model and cosmological parameter constraints.