P3M-SPH simulations of the Lyman-alpha Forest

(Abridged) We investigate the importance of several numerical artifacts such as lack of resolution on spectral properties of the Lyman alpha forest as computed from cosmological hydrodynamic simulations in a standard cold dark matter universe. We assume an ionising background produced by quasars as computed by Haardt & Madau. We use a new simulation code based on P3M and SPH, which we compare in detail with a modified version of HYDRA (Couchman et al.) and published results of TREESPH (Hernquist et al.). The agreement is very good between all three codes. We then use our new code to investigate several numerical effects such as resolution on spectral statistics deduced from Voigt profile fitting. Our highest resolution simulation has a mass resolution of 2.1x10^5 solar masses. The column density distribution is converged but the b-parameter distribution is only marginally converged. The simulation reproduces both the HI column density and b-parameter distribution when we assume a high baryon density, Omega_B h^2 > 0.028. In addition we need to impose a higher IGM temperature than predicted within our basic set of assumptions. The simulated HI optical depth is in good agreement with observations but the HeII optical depth is lower than observed. Fitting the latter requires a larger jump between the photon flux at the H and He edge than is present in the Haardt & Madau spectrum.