Ly$\alpha$ forest absorption spectra decomposed into multiple Voigt profile components (clouds) allow us to study clustering of intergalactic medium (IGM) as a function of HI column density ($N_{\rm HI}$). Here, we explore the transverse three-point correlation ($\zeta$) of these Ly$\alpha$ clouds using mock triplet spectra obtained from hydrodynamical simulations at $z \sim 2$ on scales of 1-5 $h^{-1}$cMpc. We find $\zeta$ to depend strongly on $N_{\rm HI}$ and scale and weakly on angle ($\theta$) of the triplet configuration. We show that the "hierarchical ansatz" is applicable for scales $\ge~ 3h^{-1}$cMpc, and obtain a median reduced three-point correlation (Q) in the range 0.2-0.7. We show, $\zeta$ is influenced strongly by the thermal and ionization state of the gas. As found in the case of galaxies, the influence of physical parameters on Q is weaker compared to that of $\zeta$. We show difference in $\zeta$ and Q between different simulations are minimized if we use appropriate $N_{\rm HI}$ cut-offs corresponding to a given baryon over-density ($\Delta$) using the measured $N_{\rm HI}~vs~\Delta$ relationship obtained from individual simulations. Additionally, we see the effect of pressure broadening on $\zeta$ in a model with artificially boosted heating rates. However, for models with realistic thermal and ionization histories the effect of pressure broadening on $\zeta$ is weak and sub-dominant compared to other local effects. We find strong redshift evolution shown by $\zeta$, mainly originating from the redshift evolution of thermal and ionization state of the IGM. We discuss the observational requirements for the detection of three-point correlation, specifically, in small intervals of configuration parameters and redshift.

We compare a sample of five high-resolution, high S/N Ly$\alpha$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.

Non-equilibrium ionization effects are important in cosmological hydrodynamical simulations but are computationally expensive. We study the effect of non-equilibrium ionization evolution and UV ionizing background (UVB) generated with different quasar spectral energy distribution (SED) on the derived physical conditions of the intergalactic medium (IGM) at $2\leq z \leq 6$ using our post-processing tool 'Code for Ionization and Temperature Evolution' (CITE). CITE produces results matching well with self-consistent simulations more efficiently. The HeII reionization progresses more rapidly in non-equilibrium model as compared to equilibrium models. The redshift of HeII reionization strongly depends on the quasar SED and occurs earlier for UVB models with flatter quasar SEDs. During this epoch the normalization of temperature-density relation, $T_0(z)$, has a maximum while the slope, $\gamma(z)$, has a minimum, but occurring at different redshifts. The $T_0$ is higher in non-equilibrium models using UVB obtained with flatter quasar SEDs. While our models produce the observed median HeII effective optical depth evolution and its scatter for equilibrium and non-equilibrium considerations, to explain the observed cumulative distributions we may need to consider fluctuating UVB. For a given UVB model, the redshift dependence of the HI photo-ionization rate derived from the observed HI effective optical depth ($\tau_{\rm eff,HI}$) for the equilibrium model is different from that for the non-equilibrium model. This may lead to different requirements on the evolution of ionizing emissivities of sources. We show that, in the absence of strong differential pressure smoothing effects, it is possible to recover the $T_0$ and $\gamma$ realised in non-equilibrium model from the equilibrium models generated by rescaling photo-heating rates while producing the same $\tau_{\rm eff,HI}$.

We present analysis of two- and three-point correlation functions of Ly$\alpha$ forest (at $2\le z\le 2.5$) using X-Shooter spectra of three background quasar triplets probing transverse separations of 0.5-1.6 pMpc. We present statistics based on transmitted flux and clouds identified using Voigt profile fitting. We show that the observed two-, three-point correlation functions and reduced three-point correlation (i.e Q) are well reproduced by our simulations. We assign probabilities for realising all the observed correlation properties simultaneously using our simulations. Our simulations suggest an increase in correlation amplitudes and Q with increasing $N_{\rm HI}$. We roughly see this trend in the observations too. We identify a concurrent gap of 17$\mathring{A}$ (i.e 14.2 $h^{-1}$cMpc, one of the longest reported) wide in one of the triplets. Such gap is realised only in 14.2% of our simulated sightlines and most of the time belongs to a void in the matter distribution. In the second triplet, we detect DLAs along all three sightlines (with spatial separations 0.64 to 1.6 pMpc) within a narrow redshift interval (i.e $\Delta z = 0.088$). Detection of a foreground quasar ($\sim$ 1 pMpc from the triplet sightlines) and excess partial Lyman Limit systems around these DLAs suggest that we may be probing a large over-dense region. We also report positive CIV- CIV correlations up to $\sim 500$ $km s^{-1}$ only in the longitudinal direction. Additionally, we conclude a positive CIV-Ly$\alpha$ correlations for higher $N_{\rm HI}$ thresholds up to a scale of $\sim 1000$ $km s^{-1}$ both in transverse and longitudinal directions.

We present a post-processing tool for GADGET-2 adiabatic simulations to model various observed properties of the Ly$\alpha$ forest at $2.5 \leq z \leq 4$ that enables an efficient parameter estimation. In particular, we model the thermal and ionization histories that are not computed self-consistently by default in GADGET-2. We capture the effect of pressure smoothing by running GADGET-2 at an elevated temperature floor and using an appropriate smoothing kernel. We validate our procedure by comparing different statistics derived from our method with those derived using self-consistent simulations with GADGET-3. These statistics are: line of sight density field power spectrum, flux probability distribution function, flux power spectrum, wavelet statistics, curvature statistics, HI column density (${\rm N_{HI}}$) distribution function, linewidth ($b$) distribution and $b$ versus $\log {\rm N_{HI}}$ scatter. For the temperature floor of $10^4$ K and typical signal-to-noise of 25, the results agree well within 20 percent of the self-consistent GADGET-3 simulation. However, this difference is smaller than the expected $1\sigma$ sample variance for an absorption path length of $\sim 5.35$ at $z=3$. Moreover for a given cosmology, we gain a factor of $\sim N$ in computing time for modelling the intergalactic medium under $N \gg 1$ different thermal histories. In addition, our method allows us to simulate the non-equilibrium evolution of thermal and ionization state of the gas and include heating due to non-standard sources like cosmic rays and high energy $\gamma$-rays from Blazars.

We have developed a parallel code called "VoIgt profile Parameter Estimation Routine (VIPER)" for automatically fitting the H I Ly-$\alpha$ forest seen in the spectra of QSOs. We obtained the H I column density distribution function (CDDF) and line width ($b$) parameter distribution for $z < 0.45$ using spectra of 82 QSOs obtained using Cosmic Origins Spectrograph and VIPER. Consistency of these with the existing measurements in the literature validate our code. By comparing this CDDF with those obtained from hydrodynamical simulation, we constrain the H I photoionization rate ($\Gamma_{\rm HI}$) at $z < 0.45$ in four redshift bins. The VIPER, together with the Code for Ionization and Temperature Evolution (CITE) we have developed for GADGET-2, allows us to explore parameter space and perform $\chi^2$ minimization to obtain $\Gamma_{\rm HI}$. We notice that the $b$ parameters from the simulations are smaller than what are derived from the observations. We show the observed $b$ parameter distribution and $b$ vs $\log {\rm N_{HI}}$ scatter can be reproduced in simulation by introducing sub-grid scale turbulence. However, it has very little influence on the derived $\Gamma_{\rm HI}$. The $\Gamma_{\rm HI}(z)$ obtained here, $(3.9 \pm 0.1) \times 10^{-14} \; (1+z)^{4.98 \pm 0.11} \;{\rm s^{-1}}$, is in good agreement with those derived by us using flux based statistics in the previous paper. These are consistent with the hydrogen ionizing ultra-violet (UV) background being dominated mainly by QSOs without needing any contribution from the non-standard sources of the UV photons.

We constrain the H I photoionization rate $(\Gamma_{\rm HI})$ at $z \lesssim 0.45$ by comparing the flux probability distribution function and power spectrum of the Ly-$\alpha$ forest data along 82 QSO sightlines obtained using Cosmic Origins Spectrograph with models generated from smoothed particle hydrodynamic simulations. We have developed a module named "Code for Ionization and Temperature Evolution (CITE)" for calculating the intergalactic medium (IGM) temperature evolution from high to low redshifts by post-processing the GADGET-2 simulation outputs. Our method, that produces results consistent with other simulations, is computationally less expensive thus allowing us to explore a large parameter space. It also allows rigorous estimation of the error covariance matrix for various statistical quantities of interest. We find that the best-fit $\Gamma_{\rm HI}(z)$ increases with $z$ and follows $(4 \pm 0.1) \times 10^{-14}\:(1+z)^{4.99 \pm 0.12}$ s$^{-1}$. At any given $z$ the typical uncertainties $\Delta \Gamma_{\rm HI} / \Gamma_{\rm HI}$ are $\sim 25$ per cent which contains not only the statistical errors but also those arising from possible degeneracy with the thermal history of the IGM and cosmological parameters and uncertainties in fitting the QSO continuum. These values of $\Gamma_{\rm HI}$ favour the scenario where only QSOs contribute to the ionizing background at $z<2$. Our derived $3\sigma$ upper limit on average escape fraction is $0.008$, consistent with measurements of low-$z$ galaxies.

Using our cosmological radiative transfer code, we study the implications of the updated QSO emissivity and star formation history for the escape fraction (f_esc) of hydrogen ionizing photons from galaxies. We estimate the f_esc that is required to reionize the Universe and to maintain the ionization state of the intergalactic medium in the post-reionization era. At z>5.5, we show that a constant f_esc of 0.14 to 0.22 is sufficient to reionize the Universe. At z<3.5, consistent with various observations, we find that f_esc can have values from 0 to 0.05. However, a steep rise in f_esc, of at least a factor of ~3, is required between z=3.5 to 5.5. It results from a rapidly decreasing QSO emissivity at z>3 together with a nearly constant measured H I photoionization rates at 3<z<5. We show that, this requirement of a steep rise in f_esc over a very short time can be relaxed if we consider the contribution from a recently found large number density of faint QSOs at z>4. In addition, a simple extrapolation of the contribution of such QSOs to high-z suggests that QSOs alone can reionize the Universe. This implies, at z>3.5, that either the properties of galaxies should evolve rapidly to increase the f_esc or most of the low-mass galaxies should host massive black holes and sustain accretion over a prolonged period. These results motivate a careful investigation of theoretical predictions of these alternate scenarios that can be distinguished using future observations. Moreover, it is also very important to revisit the measurements of H I photoionization rates that are crucial to the analysis presented here.