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.
Efficient and accurate simulations of the reionization epoch are crucial to exploring the vast uncharted parameter space that will soon be constrained by measurements of the 21-cm power spectrum. One of these parameters, Rmax, is meant to characterize the absorption of photons by residual neutral gas inside of ionized regions, but has historically been implemented in a very simplistic fashion acting only as a maximum distance that ionizing photons can travel. We leverage the correspondence between excursion set methods and the integrated flux from ionizing sources to define two physically motivated prescriptions of the mean free path (MFP) of ionizing photons that smoothly attenuate the contribution from distant sources. Implementation of these methods in seminumerical reionization codes requires only modest additional computational effort due to the fact that spatial filtering is still performed on scales larger than the characteristic absorption distance. We find that our smoothly defined MFP prescriptions more effectively suppress large-scale structures in the ionization field in seminumerical reionization simulations compared to the standard Rmax approach, and the magnitude of the MFP modulates the power spectrum in a much smoother manner. We show that this suppression of large-scale power is significant enough to be relevant for upcoming 21-cm power spectrum observations. Finally, we show that in our model, the MFP plays a larger role in regulating the reionization history than in models using Rmax.
more » « less- Award ID(s):
- 1812458
- NSF-PAR ID:
- 10367908
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 514
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 1302-1314
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
more » « less
Abstract Using the novel semi-numerical code for reionization AMBER, we model the patchy kinetic Sunyaev–Zel’dovich (kSZ) effect by directly specifying the reionization history with the redshift midpoint
z mid, duration Δz, and asymmetryA z. We further control the ionizing sources and radiation through the minimum halo massM hand the radiation mean free pathλ mfp. AMBER reproduces the free-electron number density and the patchy kSZ power spectrum of radiation–hydrodynamic simulations at the target resolution (1 Mpch −1) with matched reionization parameters. With a suite of (2 Gpc/h )3simulations using AMBER, we first constrain the redshift midpoint 6.0 <z mid< 8.9 using the Planck 2018 Thomson optical depth result (95% CL). Then, assumingz mid= 8, we find that the amplitude ofz and is consistent with the 1σ upper limit from South Pole Telescope (SPT) results up to Δz < 5.1 (Δz encloses 5%–95% ionization). Moreover, a shorterλ mfpcan lead to a ∼10% lowerz atℓ = 3000. Allowingz midandλ mfpto vary simultaneously, we get spectra consistent with the SPT result (95% CL) up to Δz = 12.8 (butA z > 8 is needed to ensure the end of reionization beforez = 5.5). We show that constraints on the asymmetry require ∼0.1μ k 2measurement accuracy at multipoles other thanℓ = 3000. Finally, we find that the amplitude and shape of the kSZ spectrum are only weakly sensitive toM h under a fixed reionization history and radiation mean free path. -
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
-
more » « less
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
ABSTRACT We measure the mean free path ($\lambda _{\rm mfp,H\, \small {I}}$), photoionization rate ($\langle \Gamma _{\rm H\, \small {I}} \rangle$), and neutral fraction ($\langle f_{\rm H\, \small {I}} \rangle$) of hydrogen in 12 redshift bins at 4.85 < z < 6.05 from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modelled by post-processing simulations from the Sherwood suite using our new code ‘EXtended reionization based on the Code for Ionization and Temperature Evolution’ (ex-cite). ex-cite uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution $\Gamma _{\rm H\, \small {I}}$ fluctuation models. Our simulation with ex-cite shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within 1σ uncertainty. We measure the three parameters by forward-modelling the Lyα forest and comparing the effective optical depth ($\tau _{\rm eff, H\, \small {I}}$) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modelling parameters, observational systematics, and cosmic variance. Our best-fitting parameters show significant evolution with redshift such that $\lambda _{\rm mfp,H\, \small {I}}$ and $\langle f_{\rm H\, \small {I}} \rangle$ decreases and increases by a factor ∼6 and ∼104, respectively from z ∼ 5 to z ∼ 6. By comparing our $\lambda _{\rm mfp,H\, \small {I}}$, $\langle \Gamma _{\rm H\, \small {I}} \rangle$ and $\langle f_{\rm H\, \small {I}} \rangle$ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best-fitting parameter evolution is consistent with a model in which reionization completes by z ∼ 5.2. Our best-fitting model that matches the $\tau _{\rm eff, H\, \small {I}}$ distribution also reproduces the dark gap length distribution and transmission spike height distribution suggesting robustness and accuracy of our measured parameters.
