Numerical simulations of the Lyman α forest – a comparison of gadget-2 and enzo |
| |
| Authors: | John A Regan Martin G Haehnelt Matteo Viel |
| |
| Institution: | Institute of Astronomy, Madingley Road, Cambridge CB3 0HA;Kavli Institute for Theoretical Physics, Kohn Hall, UCSB, Santa Barbara CA 93106, USA;INAF –Osservatoria Astronomico di Trieste, Via G.B. Tiepolo 11, I-34131 Trieste, Italy |
| |
| Abstract: | We compare simulations of the Lyman α forest performed with two different hydrodynamical codes, gadget-2 and enzo . A comparison of the dark matter power spectrum for simulations run with identical initial conditions show differences of 1–3 per cent at the scales relevant for quantitative studies of the Lyman α forest. This allows a meaningful comparison of the effect of the different implementations of the hydrodynamic part of the two codes. Using the same cooling and heating algorithm in both codes, the differences in the temperature and the density probability distribution function are of the order of 10 per cent. The differences are comparable to the effects of box size and resolution on these statistics. When self-converged results for each code are taken into account, the differences in the flux power spectrum – the statistics most widely used for estimating the matter power spectrum and cosmological parameters from Lyman α forest data – are about 5 per cent. This is again comparable to the effects of box size and resolution. Numerical uncertainties due to a particular implementation of solving the hydrodynamic or gravitational equations appear therefore to contribute only moderately to the error budget in estimates of the flux power spectrum from numerical simulations. We further find that the differences in the flux power spectrum for enzo simulations run with and without adaptive mesh refinement are also of the order of 5 per cent or smaller. The latter require 10 times less CPU time making the CPU time requirement similar to that of a version of gadget-2 that is optimized for Lyman α forest simulations. |
| |
| Keywords: | methods: numerical cosmology: theory large-scale structure of Universe |
|
|
