Numerical Simulations of the Lyman-alpha forest - A comparison of Gadget-2 and Enzo

Abstract

We compare simulations of the Lyman-alpha 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% at the scales relevant for quantitative studies of the Lyman-alpha 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%. These differences can be mainly attributed to a slight mismatch in the resolution. The differences are smaller than or equal to the effects of boxsize and resolution on these statistics. Once resolution effects 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-alpha forest data - are about 5%. This is again smaller than or equal to the effects of boxsize 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 order 5% 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 optimised for Lyman-alpha forest simulations.