skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Self-consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations

Abstract

The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. In this paper, we show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier ($$z\sim 15$$) than they should. This problem arises because at $$z\gt 6$$, where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyα forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. Finally, we make our new UVB photoionization and photoheating rates publicly available for use in future simulations.

Authors:
ORCiD logo [1];  [1];  [2]
  1. Max Planck Inst. for Astronomy, Heidelberg (Germany)
  2. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Max Planck Inst. for Astronomy, Heidelberg (Germany)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR) (SC-21); USDOE Office of Science (SC), High Energy Physics (HEP) (SC-25); Alexander von Humboldt Foundation (Germany); German Federal Ministry of Education and Research (BMBF)
OSTI Identifier:
1379770
Grant/Contract Number:
AC02-05CH11231
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
The Astrophysical Journal (Online)
Additional Journal Information:
Journal Name: The Astrophysical Journal (Online); Journal Volume: 837; Journal Issue: 2; Journal ID: ISSN 1538-4357
Publisher:
Institute of Physics (IOP)
Country of Publication:
United States
Language:
English
Subject:
79 ASTRONOMY AND ASTROPHYSICS; early universe; evolution of galaxies; formation of galaxies; intergalactic medium; large-scale structure of universe; numerical methods

Citation Formats

Oñorbe, Jose, Hennawi, Joseph F., and Lukić, Zarija. Self-consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations. United States: N. p., 2017. Web. doi:10.3847/1538-4357/aa6031.
Oñorbe, Jose, Hennawi, Joseph F., & Lukić, Zarija. Self-consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations. United States. doi:10.3847/1538-4357/aa6031.
Oñorbe, Jose, Hennawi, Joseph F., and Lukić, Zarija. Wed . "Self-consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations". United States. doi:10.3847/1538-4357/aa6031. https://www.osti.gov/servlets/purl/1379770.
@article{osti_1379770,
title = {Self-consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations},
author = {Oñorbe, Jose and Hennawi, Joseph F. and Lukić, Zarija},
abstractNote = {The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. In this paper, we show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier ($z\sim 15$) than they should. This problem arises because at $z\gt 6$, where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Lyα forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. Finally, we make our new UVB photoionization and photoheating rates publicly available for use in future simulations.},
doi = {10.3847/1538-4357/aa6031},
journal = {The Astrophysical Journal (Online)},
number = 2,
volume = 837,
place = {United States},
year = {Wed Mar 08 00:00:00 EST 2017},
month = {Wed Mar 08 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 8works
Citation information provided by
Web of Science

Save / Share:
  • We perform cosmological N-body simulations of the Dvali-Gabadadze-Porrati (DGP) braneworld model, by solving the full nonlinear equations of motion for the scalar degree of freedom in this model, the brane-bending mode. While coupling universally to matter, the brane-bending mode has self-interactions that become important as soon as the density field becomes nonlinear. These self-interactions lead to a suppression of the field in high-density environments, and restore gravity to general relativity. The code uses a multigrid relaxation scheme to solve the nonlinear field equation in the quasistatic approximation. We perform simulations of a flat self-accelerating DGP model without cosmological constant. However,more » the type of nonlinear interactions of the brane-bending mode, which are the focus of this study, are generic to a wide class of braneworld cosmologies. The results of the DGP simulations are compared with standard gravity simulations assuming the same expansion history, and with DGP simulations using the linearized equation for the brane-bending mode. This allows us to isolate the effects of the nonlinear self-couplings of the field which are noticeable already on quasilinear scales. We present results on the matter power spectrum and the halo mass function, and discuss the behavior of the brane-bending mode within cosmological structure formation. We find that, independently of cosmic microwave background constraints, the self-accelerating DGP model is strongly constrained by current weak lensing and cluster abundance measurements.« less
  • We examine the validity of the hydrostatic equilibrium (HSE) assumption for galaxy clusters using one of the highest-resolution cosmological hydrodynamical simulations. We define and evaluate several effective mass terms corresponding to the Euler equations of gas dynamics, and quantify the degree of the validity of HSE in terms of the mass estimate. We find that the mass estimated under the HSE assumption (the HSE mass) deviates from the true mass by up to {approx}30%. This level of departure from HSE is consistent with the previous claims, but our physical interpretation is rather different. We demonstrate that the inertial term inmore » the Euler equations makes a negligible contribution to the total mass, and the overall gravity of the cluster is balanced by the thermal gas pressure gradient and the gas acceleration term. Indeed, the deviation from the HSE mass is well explained by the acceleration term at almost all radii. We also clarify the confusion of previous work due to the inappropriate application of the Jeans equations in considering the validity of HSE from the gas dynamics extracted from cosmological hydrodynamical simulations.« less
  • Self-consistent models for elliptical galaxies have been constructed and applied to a sample of 16 observed E systems. The models are based on the hydrodynamical Jeans equations, and the E galaxies are assumed to be oblate ellipsoids of constant ellipticity. The majority of the galaxies in the sample are found to be anisotropic. Most of the anisotropic systems have greater velocity disperison in the equatorial plane than in the direction of the galaxy's symmetry axis, so that their flattenings are substantially determined by velocity anisotropy. It is found that half of the sample is well described by models with M/Lmore » ratios that increase with the distance from the galactic center. Such behavior is indicative of the existence of nonluminous matter in those early-type galaxies. The possibity that this result may be quite general for ellipticals is also discussed. 54 refs.« less
  • The problem of self-consistent solutions in higher-dimensional spacetimes is considered. The importance of quantum effects in this analysis is emphasized. It is shown how the one-loop effective action for conformally invariant matter fields on a time-dependent Kaluza-Klein background may be obtained from previously known results on static backgrounds. The field equations are obtained and solved for self-consistent solutions in cases for which the four-dimensional part of the space is a Friedmann-Robertson-Walker space. Stability of the solutions to small perturbations is studied. All solutions, including the product of flat four-dimensional space with a sphere, are shown to be unstable.