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 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.