Efficient Modelling of Lyman-{\alpha} opacity fluctuations during late reionization epoch

TL;DR

This paper introduces an efficient semi-numerical model for simulating Lyman-alpha opacity fluctuations during late reionization, matching observational data and constraining key reionization parameters.

Contribution

It develops a coarse simulation approach incorporating ionization and temperature fluctuations, calibrated with high-resolution hydrodynamic simulations, to analyze reionization epoch data.

Findings

Model matches observed opacity scatter at z=5.0-6.1

Constraints suggest higher photoionization rate and lower mean free path

Temperature fluctuations influence reionization parameter estimates

Abstract

The Lyman-$\alpha$ forest opacity fluctuations observed from high-redshift quasar spectra have been proven to be extremely successful in order to probe the late phase of the reionization epoch. For ideal modeling of these opacity fluctuations, one of the main challenges is to satisfy the extremely high dynamic range requirements of the simulation box, resolving the Lyman-$\alpha$ forest while probing the large cosmological scales. In this study, we adopt an efficient approach to model Lyman-$\alpha$ opacity fluctuations in a coarse simulation volume, utilizing the semi-numerical reionization model SCRIPT (including inhomogeneous recombination and radiative feedback) integrated with a realistic photoionization background fluctuation generating model. Our model crucially incorporates ionization and temperature fluctuations, which are consistent with the reionization model. After calibrating our method with respect to high-resolution full hydrodynamic simulation, Nyx, we compared the models with available observational data at the redshift range, $z=5.0-6.1$. With a fiducial reionization model (reionization end at $z=5.8$), we demonstrated that the observed scatter in the effective optical depth can be matched reasonably well by tuning the free parameters of our model, i.e., the effective ionizing photon mean free path and mean photoionization rate. We further pursued an MCMC-based parameter space exploration, utilizing the available data to put constraints on the above free parameters. Our estimation prefers a slightly higher photoionization rate and slightly lower mean free path than the previous studies, which is also a consequence of temperature fluctuations. This study holds significant promise for efficiently extracting important physical information about the Epoch of Reionization, utilizing the wealth of available and upcoming observational data.