The Lyman α forest in optically thin hydrodynamical simulations
Abstract
We study the statistics of the Lyα forest in a flat Λ cold dark matter cosmology with the Nbody + Eulerian hydrodynamics code NYX. We produce a suite of simulations, covering the observationally relevant redshift range 2 ≤ z ≤ 4. We find that a grid resolution of 20 h^{1} kpc is required to produce 1 per cent convergence of Lyα forest flux statistics, up to k = 10 h^{1} Mpc. In addition to establishing resolution requirements, we study the effects of missing modes in these simulations, and find that box sizes of L > 40h^{1} Mpc are needed to suppress numerical errors to a subper cent level. Our optically thin simulations with the ionizing background prescription of Haardt & Madau reproduce an intergalactic medium densitytemperature relation with T_{0} ≈ 10^{4} K and γ ≈ 1.55 at z = 2, with a mean transmitted flux close to the observed values. When using the ionizing background prescription of FaucherGiguère et al., the mean flux is 1015 per cent below observed values at z = 2, and a factor of 2 too small at z = 4. We show the effects of the common practice of rescaling optical depths to the observed mean flux and how it affects convergence rates. We also investigate the practice of `splicing' results from a number of different simulations to estimate the 1D flux power spectrum and show it is accurate at the 10 per cent level. Finally, we find that collisional heating of the gas from dark matter particles is negligible in modern cosmological simulations.
 Publication:

Monthly Notices of the Royal Astronomical Society
 Pub Date:
 February 2015
 DOI:
 10.1093/mnras/stu2377
 arXiv:
 arXiv:1406.6361
 Bibcode:
 2015MNRAS.446.3697L
 Keywords:

 methods: numerical;
 intergalactic medium;
 quasars: absorption lines;
 largescale structure of universe;
 Astrophysics  Cosmology and Nongalactic Astrophysics
 EPrint:
 28 pages, 31 figures