skip to main content


Title: The impact of inhomogeneous subgrid clumping on cosmic reionization
ABSTRACT

Cosmic reionization was driven by the imbalance between early sources and sinks of ionizing radiation, both of which were dominated by small-scale structure and are thus usually treated in cosmological reionization simulations by subgrid modelling. The recombination rate of intergalactic hydrogen is customarily boosted by a subgrid clumping factor, 〈n2〉/〈n〉2, which corrects for unresolved fluctuations in gas density n on scales below the grid-spacing of coarse-grained simulations. We investigate in detail the impact of this inhomogeneous subgrid clumping on reionization and its observables, as follows: (1) Previous attempts generally underestimated the clumping factor because of insufficient mass resolution. We perform a high-resolution N-body simulation that resolves haloes down to the pre-reionization Jeans mass to derive the time-dependent, spatially varying local clumping factor and a fitting formula for its correlation with local overdensity. (2) We then perform a large-scale N-body and radiative transfer simulation that accounts for this inhomogeneous subgrid clumping by applying this clumping factor-overdensity correlation. Boosting recombination significantly slows the expansion of ionized regions, which delays completion of reionization and suppresses 21 cm power spectra on large scales in the later stages of reionization. (3) We also consider a simplified prescription in which the globally averaged, time-evolving clumping factor from the same high-resolution N-body simulation is applied uniformly to all cells in the reionization simulation, instead. Observables computed with this model agree fairly well with those from the inhomogeneous clumping model, e.g. predicting 21 cm power spectra to within 20 per cent error, suggesting it may be a useful approximation.

 
more » « less
PAR ID:
10125954
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
491
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 1600-1621
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT A recent measurement of the Lyman-limit mean free path at z = 6 suggests it may have been very short, motivating a better understanding of the role that ionizing photon sinks played in reionization. Accurately modelling the sinks in reionization simulations is challenging because of the large dynamic range required if ∼104−108M⊙ gas structures contributed significant opacity. Thus, there is no consensus on how important the sinks were in shaping reionization’s morphology. We address this question with a recently developed radiative transfer code that includes a dynamical sub-grid model for the sinks based on radiative hydrodynamics simulations. Compared to assuming a fully pressure-smoothed intergalactic medium, our dynamical treatment reduces ionized bubble sizes by $10-20~{{\ \rm per\ cent}}$ under typical assumptions about reionization’s sources. Near reionization’s midpoint, the 21 cm power at k ∼ 0.1 hMpc−1 is similarly reduced. These effects are more modest than the $30-60~{{\ \rm per\ cent}}$ suppression resulting from the higher recombination rate if pressure smoothing is neglected entirely. Whether the sinks played a significant role in reionization’s morphology depends on the nature of its sources. For example, if reionization was driven by bright (MUV < −17) galaxies, the sinks reduce the large-scale 21 cm power by at most 20  per cent, even if pressure smoothing is neglected. Conveniently, when bright sources contribute significantly, the morphology in our dynamical treatment can be reproduced accurately with a uniform sub-grid clumping factor that yields the same ionizing photon budget. By contrast, if MUV ∼ −13 galaxies drove reionization, the uniform clumping model can err by up to 40  per cent. 
    more » « less
  2. ABSTRACT

    The formation of the first galaxies during cosmic dawn and reionization (at redshifts z = 5–30), triggered the last major phase transition of our universe, as hydrogen evolved from cold and neutral to hot and ionized. The 21-cm line of neutral hydrogen will soon allow us to map these cosmic milestones and study the galaxies that drove them. To aid in interpreting these observations, we upgrade the publicly available code 21cmFAST. We introduce a new, flexible parametrization of the additive feedback from: an inhomogeneous, H2-dissociating (Lyman–Werner; LW) background; and dark matter – baryon relative velocities; which recovers results from recent, small-scale hydrodynamical simulations with both effects. We perform a large, ‘best-guess’ simulation as the 2021 installment of the Evolution of 21-cm Structure (EOS) project. This improves the previous release with a galaxy model that reproduces the observed UV luminosity functions (UVLFs), and by including a population of molecular-cooling galaxies. The resulting 21-cm global signal and power spectrum are significantly weaker, primarily due to a more rapid evolution of the star formation rate density required to match the UVLFs. Nevertheless, we forecast high signal-to-noise detections for both HERA and the SKA. We demonstrate how the stellar-to-halo mass relation of the unseen, first galaxies can be inferred from the 21-cm evolution. Finally, we show that the spatial modulation of X-ray heating due to relative velocities provides a unique acoustic signature that is detectable at z ≈ 10–15 in our fiducial model. Ours are the first public simulations with joint inhomogeneous LW and relative-velocity feedback across the entire cosmic dawn and reionization, and we make them available at this link https://scholar.harvard.edu/julianbmunoz/eos-21.

     
    more » « less
  3. ABSTRACT

    We introduce the thesan project, a suite of large volume ($L_\mathrm{box} = 95.5 \, \mathrm{cMpc}$) radiation-magnetohydrodynamic simulations that simultaneously model the large-scale statistical properties of the intergalactic medium during reionization and the resolved characteristics of the galaxies responsible for it. The flagship simulation has dark matter and baryonic mass resolutions of $3.1 \times 10^6\, {\rm M_\odot }$ and $5.8 \times 10^5\, {\rm M_\odot }$, respectively. The gravitational forces are softened on scales of 2.2 ckpc with the smallest cell sizes reaching 10 pc at z = 5.5, enabling predictions down to the atomic cooling limit. The simulations use an efficient radiation hydrodynamics solver (arepo-rt) that precisely captures the interaction between ionizing photons and gas, coupled to well-tested galaxy formation (IllustrisTNG) and dust models to accurately predict the properties of galaxies. Through a complementary set of medium resolution simulations we investigate the changes to reionization introduced by different assumptions for ionizing escape fractions, varying dark matter models, and numerical convergence. The fiducial simulation and model variations are calibrated to produce realistic reionization histories that match the observed evolution of the global neutral hydrogen fraction and electron scattering optical depth to reionization. They also match a wealth of high-redshift observationally inferred data, including the stellar-to-halo-mass relation, galaxy stellar mass function, star formation rate density, and the mass–metallicity relation, despite the galaxy formation model being mainly calibrated at z = 0. We demonstrate that different reionization models give rise to varied bubble size distributions that imprint unique signatures on the 21 cm emission, especially on the slope of the power spectrum at large spatial scales, enabling current and upcoming 21 cm experiments to accurately characterize the sources that dominate the ionizing photon budget.

     
    more » « less
  4. Abstract

    The cross-correlation between the 21 cm field and the galaxy distribution is a potential probe of the Epoch of Reionization (EoR). The 21 cm signal traces neutral gas in the intergalactic medium and, on large spatial scales, this should be anticorrelated with the high-redshift galaxy distribution, which partly sources and tracks the ionized gas. In the near future, interferometers such as the Hydrogen Epoch of Reionization Array (HERA) are projected to provide extremely sensitive measurements of the 21 cm power spectrum. At the same time, the Nancy Grace Roman Space Telescope (Roman) will produce the most extensive catalog to date of bright galaxies from the EoR. Using seminumeric simulations of reionization, we explore the prospects for measuring the cross-power spectrum between the 21 cm and galaxy fields during the EoR. We forecast a 12σdetection between HERA and Roman, assuming an overlapping survey area of 500 deg2, redshift uncertainties ofσz= 0.01 (as expected for the high-latitude spectroscopic survey of Lyα-emitting galaxies), and an effective Lyαemitter duty cycle offLAE= 0.1. Thus the HERA–Roman cross-power spectrum may be used to help verify 21 cm detections from HERA. We find that the shot-noise in the galaxy distribution is a limiting factor for detection, and so supplemental observations using Roman should prioritize deeper observations, rather than covering a wider field of view. We have made a public GitHub repository containing key parts of the calculation, which accompanies this paper:https://github.com/plaplant/21cm_gal_cross_correlation.

     
    more » « less
  5. ABSTRACT

    Preparing for a first detection of the 21-cm signal during reionization by large-scale interferometer experiments requires rigorous testing of the data analysis and reduction pipelines. Validating that these do not erroneously add/remove features mimicking the signal (e.g. from side lobes or large-scale power leakage) requires simulations extending beyond the primary field of view. However, the Murchison Wide Field Array (MWA) with a field of view of ∼252 deg2 would require simulations spanning several Gpcs, which are currently infeasible. To address this, we developed a simplified version of the seminumerical reionization simulation code 21cmfast, sacrificing some physical accuracy (linear structure formation) in favour of extremely large volumes. We then constructed a 7.5 Gpc comoving volume specifically tailored to the binned spectral resolution of the MWA (∼1.17 cMpc), required for validating the pipeline used in the 2020 MWA 21-cm power spectrum (PS) upper limits. With this large-volume simulation, we then explored: (i) whether smaller volume simulations are biased by missing large-scale modes, (ii) non-Gaussianity in the cosmic variance uncertainty, (iii) biases in the recovered 21-cm PS following foreground wedge avoidance, and (iv) the impact of tiling smaller simulations to achieve large volumes. We found (i) no biases from missing large-scale power, (ii) significant contribution from non-Gaussianity, as expected, (iii) a 10–20 per cent overestimate of the 21-cm PS following wedge mode excision, and (iv) tiling smaller simulations underestimates the large-scale power and cosmic variance.

     
    more » « less