Particle initialization effects on Lyman-α forest statistics in cosmological SPH simulations

Confronting measurements of the Lyman-α forest with cosmological hydrodynamical simulations has produced stringent constraints on models of particle dark matter and the thermal and ionization state of the intergalactic medium. We investigate the robustness of such models of the Lyman-α forest, focussing on the effect of particle initial conditions on the Lyman-α forest statistics in cosmological SPH simulations. We study multiple particle initialization algorithms in simulations that are designed to be identical in other respects. In agreement with the literature, we find that the correct linear theory evolution is obtained when a glass-like configuration is used for initial unperturbed gas particle positions alongside a regular grid configuration for dark matter particles and the use of non-identical initial density perturbations for gas and dark matter. However, we report that this introduces a large scale-dependent distortion in the one-dimensional Lyman-α transmission power spectrum at small scales (k > 0.05 s/km). The effect is close to 50 % at k ∼ 0.1 s/km, and persists at higher resolution. This can severely bias inferences in parameters such as the dark matter particle mass. By considering multiple initial conditions codes and their variations, we also study the impact of a variety of other assumptions and algorithmic choices, such as adaptive softening, background radiation density, particle staggering, and perturbation theory accuracy, on the matter power spectrum, the Lyman-α flux power spectrum, and the Lyman-α flux PDF. This work reveals possible pathways towards more accurate theoretical models of the Lyman-α forest to match the quality of upcoming measurements

Measuring the photo-ionization rate, neutral fraction and mean free path of HI ionizing photons at 4.9≤z≤6.04.9 \leq z \leq 6.0 from a large sample of XShooter and ESI spectra

We measure the mean free path ($\lambda_{\rm mfp,HI}$), photo-ionization rate ($\langle \Gamma_{\rm HI} \rangle$) and neutral fraction ($\langle f_{\rm HI} \rangle$) of hydrogen in 12 redshift bins at $4.85<z<6.05$ from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modeled by post-processing simulations from the Sherwood suite using our new code ''EXtended reionization based on the Code for Ionization and Temperature Evolution'' (EX-CITE). EX-CITE uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution $\Gamma_{\rm HI}$ fluctuation models. Our simulation with EX-CITE shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within $1\sigma$ uncertainty. We measure the three parameters by forward-modeling the Ly$\alpha$ forest and comparing the effective optical depth ($\tau_{\rm eff, HI}$) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modeling parameters, observational systematics and cosmic variance. Our best fit parameters show significant evolution with redshift such that $\lambda_{\rm mfp,HI}$ and $\langle f_{\rm HI} \rangle$ decreases and increases by a factor $\sim 6$ and $\sim 10^{4}$, respectively from $z \sim 5$ to $z \sim 6$. By comparing our $\lambda_{\rm mfp,HI}$, $\langle \Gamma_{\rm HI} \rangle$ and $\langle f_{\rm HI} \rangle$ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best fit parameter evolution is consistent with a model in which reionization completes by $z \sim 5.2$.Comment: 30 pages (+14 pages appendices), 14 figures (+18 figures appendices); submitted to MNRAS; Main results are summarized in Fig. 10, Fig. 11 and Table

Measuring the photo-ionization rate, neutral fraction and mean free path of HI ionizing photons at 4.9≤z≤6.0 from a large sample of XShooter and ESI spectra

We measure the mean free path (⁠λmfp,HI⁠), photo-ionization rate (⁠⟨ΓHI⟩⁠) and neutral fraction (⁠⟨fHI⟩⁠) of hydrogen in 12 redshift bins at 4.85 < z < 6.05 from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modeled by post-processing simulations from the Sherwood suite using our new code ‘EXtended reionization based on the Code for Ionization and Temperature Evolution’ (EX-CITE). EX-CITE uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution ΓHI fluctuation models. Our simulation with EX-CITE shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within 1σ uncertainty. We measure the three parameters by forward-modeling the Lyα forest and comparing the effective optical depth (⁠τeff,HI⁠) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modeling parameters, observational systematics and cosmic variance. Our best fit parameters show significant evolution with redshift such that λmfp,HI and ⟨fHI⟩ decreases and increases by a factor ∼6 and ∼104, respectively from z ∼ 5 to z ∼ 6. By comparing our λmfp,HI⁠, ⟨ΓHI⟩ and ⟨fHI⟩ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best fit parameter evolution is consistent with a model in which reionization completes by z ∼ 5.2. Our best fit model that matches the τeff,HI distribution also reproduces the dark gap length distribution and transmission spike height distribution suggesting robustness and accuracy of our measured parameters

Late-end reionization with ATON-HE: towards constraints from Lyman-α\alpha emitters observed with JWST

International audienceWe present a new suite of late-end reionization simulations performed with ATON-HE, a revised version of the GPU-based radiative transfer code ATON that includes helium. The simulations are able to reproduce the Ly$\alpha$ flux distribution of the E-XQR-30 sample of QSO absorption spectra at $5 \lesssim z \lesssim 6.2$, and show that a large variety of reionization models are consistent with these data. We explore a range of variations in source models and in the early-stage evolution of reionization. Our fiducial reionization history has a midpoint of reionization at $z = 6.5$, but we also explore an `Early' reionization history with a midpoint at $z = 7.5$ and an `Extremely Early' reionization history with a midpoint at $z = 9.5$. Haloes massive enough to host observed Ly$\alpha$ emitters are highly biased. The fraction of such haloes embedded in ionized bubbles that are large enough to allow high Ly$\alpha$ transmission becomes close to unity much before the volume filling factor of ionized regions. For our fiducial reionization history this happens at $z = 8$, probably too late to be consistent with the detection by JWST of abundant Ly$\alpha$ emission out to $z = 11$. A reionization history in our `Early' model or perhaps even our `Extremely Early' model may be required, suggesting a Thomson scattering optical depth in tension with that reported by Planck, but consistent with recent suggestions of a significantly higher value