Modeling the Lyman-$\alpha$ forest with Eulerian and SPH hydrodynamical methods

Abstract

We compare two state-of-the-art numerical codes to study the overall accuracy in modeling the intergalactic medium and reproducing Lyman-$\alpha$ forest observables for DESI and high-resolution data sets. The codes employ different approaches to solving both gravity and modeling the gas hydrodynamics. The first code, Nyx, solves the Poisson equation using the Particle-Mesh (PM) method and the Euler equations using a finite volume method. The second code, \CRKHACC, uses a Tree-PM method to solve for gravity, and an improved Lagrangian smoothed particle hydrodynamics (SPH) technique, where fluid elements are modeled with particles, to treat the intergalactic gas. We compare the convergence behavior of the codes in flux statistics as well as the degree to which the codes agree in the converged limit. We find good agreement overall with differences being less than observational uncertainties, and a particularly notable $\lesssim$1\% agreement in the 1D flux power spectrum. This agreement was achieved by applying a tessellation methodology for reconstructing the density in \CRKHACC instead of using an SPH kernel as is standard practice. We show that use of the SPH kernel can lead to significant and unnecessary biases in flux statistics; this is especially prominent at high redshifts, $z\sim 5$, as the Lyman-$\alpha$ forest mostly comes from lower-density regions which are intrinsically poorly sampled by SPH particles.

submitted to MNRAS

References in corpus (19)

• The DESI Experiment Part I: Science,Targeting, and Survey Design
• Baryon acoustic oscillations from the complete SDSS-III Ly$\alpha$-quasar cross-correlation function at $z=2.4$
• How cold is cold dark matter? Small scales constraints from the flux power spectrum of the high-redshift Lyman-alpha forest
• Constraints from Ly-$\alpha$ forests on non-thermal dark matter including resonantly-produced sterile neutrinos
• Possible evidence for an inverted temperature-density relation in the intergalactic medium from the flux distribution of the Lyman-alpha forest
• An improved measurement of the flux distribution of the Ly-alpha forest in QSO absorption spectra: the effect of continuum fitting, metal contamination and noise properties
• Constraint on neutrino masses from SDSS-III/BOSS Ly$\alpha$ forest and other cosmological probes
• An Exact Integration Scheme for Radiative Cooling in Hydrodynamical Simulations
• Self-Consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations
• The Lyman-alpha forest power spectrum from the XQ-100 Legacy Survey
• IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function
• A New Precision Measurement of the Small-Scale Line-of-Sight Power Spectrum of the Ly$\alpha$ Forest
• How to estimate the 3D power spectrum of the Lyman-$\alpha$ forest
• Numerical Simulations of the Lyman-alpha forest - A comparison of Gadget-2 and Enzo
• Exploring the Thermal State of the Low-Density Intergalactic Medium at z=3 with an Ultra-High Signal-to-Noise QSO Spectrum
• The Probability Distribution of the Lya transmitted flux from a sample of SDSS quasars
• Simulating intergalactic gas for DESI-like small scale Lymanα forest observations
• WEAVE-QSO: A Massive Intergalactic Medium Survey for the William Herschel Telescope
• Extracting the astrophysics of reionization from the Ly$\alpha$ forest power spectrum: a first forecast