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

Title: Elucid—exploring the local universe with the reconstructed initial density field. I. Hamiltonian Markov chain Monte Carlo method with particle mesh dynamics

Abstract

Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input nonlinear density field, employing the Hamiltonian Markov Chain Monte Carlo (HMC) algorithm combined with Particle-mesh (PM) dynamics. The HMC+PM method is applied to cosmological simulations, and the reconstructed linear density fields are then evolved to the present day with N-body simulations. These constrained simulations accurately reproduce both the amplitudes and phases of the input simulations at various z. Using a PM model with a grid cell size of 0.75 h {sup –1} Mpc and 40 time steps in the HMC can recover more than half of the phase information down to a scale k ∼ 0.85 h Mpc{sup –1} at high z and to k ∼ 3.4 h Mpc{sup –1} at z = 0, which represents a significant improvement over similar reconstruction models in the literature, and indicates that our model can reconstruct the formation histories of cosmic structures over a large dynamical range. Adopting PM models with higher spatial and temporal resolutions yields even better reconstructions, suggesting that our method ismore » limited more by the availability of computer resource than by principle. Dynamic models of structure evolution adopted in many earlier investigations can induce non-Gaussianity in the reconstructed linear density field, which in turn can cause large systematic deviations in the predicted halo mass function. Such deviations are greatly reduced or absent in our reconstruction.« less

Authors:
;  [1]; ;  [2]
  1. Key Laboratory for Research in Galaxies and Cosmology, Department of Astronomy, University of Science and Technology of China, Hefei, Anhui 230026 (China)
  2. Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Shanghai 200030 (China)
Publication Date:
OSTI Identifier:
22370451
Resource Type:
Journal Article
Journal Name:
Astrophysical Journal
Additional Journal Information:
Journal Volume: 794; Journal Issue: 1; Other Information: Country of input: International Atomic Energy Agency (IAEA); Journal ID: ISSN 0004-637X
Country of Publication:
United States
Language:
English
Subject:
79 ASTROPHYSICS, COSMOLOGY AND ASTRONOMY; ALGORITHMS; AMPLITUDES; DENSITY; EVOLUTION; GALAXIES; HAMILTONIANS; MARKOV PROCESS; MASS; MONTE CARLO METHOD; NONLINEAR PROBLEMS; NONLUMINOUS MATTER; PARTICLES; RESOLUTION; SIMULATION; UNIVERSE

Citation Formats

Wang, Huiyuan, Mo, H. J., Yang, Xiaohu, Lin, W. P., and Jing, Y. P., E-mail: whywang@mail.ustc.edu.cn. Elucid—exploring the local universe with the reconstructed initial density field. I. Hamiltonian Markov chain Monte Carlo method with particle mesh dynamics. United States: N. p., 2014. Web. doi:10.1088/0004-637X/794/1/94.
Wang, Huiyuan, Mo, H. J., Yang, Xiaohu, Lin, W. P., & Jing, Y. P., E-mail: whywang@mail.ustc.edu.cn. Elucid—exploring the local universe with the reconstructed initial density field. I. Hamiltonian Markov chain Monte Carlo method with particle mesh dynamics. United States. https://doi.org/10.1088/0004-637X/794/1/94
Wang, Huiyuan, Mo, H. J., Yang, Xiaohu, Lin, W. P., and Jing, Y. P., E-mail: whywang@mail.ustc.edu.cn. 2014. "Elucid—exploring the local universe with the reconstructed initial density field. I. Hamiltonian Markov chain Monte Carlo method with particle mesh dynamics". United States. https://doi.org/10.1088/0004-637X/794/1/94.
@article{osti_22370451,
title = {Elucid—exploring the local universe with the reconstructed initial density field. I. Hamiltonian Markov chain Monte Carlo method with particle mesh dynamics},
author = {Wang, Huiyuan and Mo, H. J. and Yang, Xiaohu and Lin, W. P. and Jing, Y. P., E-mail: whywang@mail.ustc.edu.cn},
abstractNote = {Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input nonlinear density field, employing the Hamiltonian Markov Chain Monte Carlo (HMC) algorithm combined with Particle-mesh (PM) dynamics. The HMC+PM method is applied to cosmological simulations, and the reconstructed linear density fields are then evolved to the present day with N-body simulations. These constrained simulations accurately reproduce both the amplitudes and phases of the input simulations at various z. Using a PM model with a grid cell size of 0.75 h {sup –1} Mpc and 40 time steps in the HMC can recover more than half of the phase information down to a scale k ∼ 0.85 h Mpc{sup –1} at high z and to k ∼ 3.4 h Mpc{sup –1} at z = 0, which represents a significant improvement over similar reconstruction models in the literature, and indicates that our model can reconstruct the formation histories of cosmic structures over a large dynamical range. Adopting PM models with higher spatial and temporal resolutions yields even better reconstructions, suggesting that our method is limited more by the availability of computer resource than by principle. Dynamic models of structure evolution adopted in many earlier investigations can induce non-Gaussianity in the reconstructed linear density field, which in turn can cause large systematic deviations in the predicted halo mass function. Such deviations are greatly reduced or absent in our reconstruction.},
doi = {10.1088/0004-637X/794/1/94},
url = {https://www.osti.gov/biblio/22370451}, journal = {Astrophysical Journal},
issn = {0004-637X},
number = 1,
volume = 794,
place = {United States},
year = {Fri Oct 10 00:00:00 EDT 2014},
month = {Fri Oct 10 00:00:00 EDT 2014}
}