Using cosmological hydrodynamic simulations of the LCDM model, we present a comparison between the simulation sample and real data sample of HI and HeII Ly\alpha transmitted flux in the absorption spectra of the QSO HE2347-4342. The LCDM model is successful in simultaneously explaining the statistical features of both HI and HeII Ly\alpha transmitted flux. It includes: 1.) the power spectra of the transmitted flux of HI and HeII can be well fitted on all scales > 0.28h^{-1} Mpc for H, and > 1.1h^{-1} Mpc for He; 2.) the Doppler parameters of absorption features of HeII and HI are found to be turbulent-broadening; 3.) the ratio of HeII to HI optical depths are substantially scattered, due to the significant effect of noise. A large part of the \eta-scatter is due to the noise in the HeII flux. However, the real data contain more low-\eta events than simulation sample. This discrepancy may indicate that the mechanism leading extra fluctuations upon the simulation data, such as a fluctuating UV radiation background, is needed. Yet, models of these extra fluctuations should satisfy the constraints: 1.) if the fluctuations are Gaussian, they should be limited by the power spectra of observed HI and HeII flux; 2.) if the fluctuations are non-Gaussian, they should be limited by the observed non-Gaussian features of the HI and HeII flux.

We have studied the power spectrum and the intermittent behavior of the fluctuations in the transmitted flux of HE2347-4342 ${\rm Ly}{\alpha}$ absorption in order to investigate if there is any discrepancy between the LCDM model with parameters given by the WMAP and observations on small scales. If the non-Gaussianity of cosmic mass field is assumed to come only from halos with an universal mass profile of the LCDM model, the non-Gaussian behavior of mass field would be effectively measured by its intermittency, because intermittency is a basic statistical feature of the cuspy structures. We have shown that the Ly$\alpha$ transmitted flux field of HE2347-4342 is significantly intermittent on small scales. With the hydrodynamic simulation, we demonstrate that the LCDM model is successful in explaining the power spectrum and intermittency of ${\rm Ly}{\alpha}$ transmitted flux. Using statistics ranging from the second to eighth order, we find no discrepancy between the LCDM model and the observed transmitted flux field, and no evidence to support the necessity of reducing the power of density perturbations relative to the standard LCDM model up to comoving scales as small as about $0.08 {\rm h}^{-1} {\rm Mpc}$. Moreover, our simulation samples show that the intermittent exponent of the Ly$\alpha$ transmitted flux field is probably scale-dependent. This result is different from the prediction of universal mass profile with a constant index of the central cusp. The scale-dependence of the intermittent exponent indicates that the distribution of baryonic gas is decoupled from the underlying dark matter.

Using a set of 28 high resolution, high signal to noise ratio (S/N) QSO Ly-alpha absorption spectra, we investigate the non-Gaussian features of the transmitted flux fluctuations, and their effect upon the power spectrum of this field. We find that the spatial distribution of the local power of the transmitted flux on scales k >= 0.05 s/km is highly spiky or intermittent. The probability distribution functions (PDFs) of the local power are long-tailed. The power on small scales is dominated by small probability events, and consequently, the uncertainty in the power spectrum of the transmitted flux field is generally large. This uncertainty arises due to the slow convergence of an intermittent field to a Gaussian limit required by the central limit theorem (CLT). To reduce this uncertainty, it is common to estimate the error of the power spectrum by selecting subsamples with an "optimal" size. We show that this conventional method actually does not calculate the variance of the original intermittent field but of a Gaussian field. Based on the analysis of intermittency, we propose an algorithm to calculate the error. It is based on a bootstrap re-sampling among all independent local power modes. This estimation doesn't require any extra parameter like the size of the subsamples, and is sensitive to the intermittency of the fields. This method effectively reduces the uncertainty in the power spectrum when the number of independent modes matches the condition of the CLT convergence.

We calculate the structure function and intermittent exponent of the 1.) Keck data, which consists of 29 high resolution, high signal to noise ratio (S/N) QSO Ly$\alpha$ absorption spectra, and 2.)the Ly$\alpha$ forest simulation samples produced via the pseudo hydro scheme for the low density cold dark matter (LCDM) model and warm dark matter (WDM) model with particle mass $m_W=300, 600, 800$ and 1000 eV. These two measures detect not only non-gaussianities, but also the type of non-gaussianty in the the field. We find that, 1.) the structure functions of the simulation samples are significantly larger than that of Keck data on scales less than about 100 h$^{-1}$ kpc, 2.) the intermittent exponent of the simulation samples is more negative than that of Keck data on all redshifts considered, 3.) the order-dependence of the structure functions of simulation samples are closer to the intermittency of hierarchical clustering on all scales, while the Keck data are closer to a lognormal field on small scales. These differences are independent of noise and show that the intermittent evolution modeled by the pseudo-hydro simulation is substantially different from observations, even though they are in good agreement in terms of second and lower order statistics. (Abridged)

The calculation of the transmission power spectrum of QSO's Lyman alpha absorption requires two parameters for the normalization: the continuum Fc and mean transmission, i.e. average of e^{-tau}. Traditionally, the continuum is obtained by a polynomial fitting truncating it at a lower order, and the mean transmission is calculated over the entire wavelength range considered. The flux F is then normalized by the average of Fc e^{-tau}. However, the fluctuations in the transmitted flux are significantly correlated with the local background flux on scales for which the field is intermittent. In this paper, we develop a self-normalization algorithm of the transmission power spectrum based on a multiresolution analysis. This self-normalized power spectrum estimator needs neither a continuum fitting, nor pre-determining the mean transmission. With simulated samples, we show that the self-normalization algorithm can perfectly recover the transmission power spectrum from the flux regardless of how the continuum varies with wavelength. We also show that the self-normalized power spectrum is also properly normalized by the mean transmission. Moreover, this power spectrum estimator is sensitive to the non-linear behavior of the field. That is, the self-normalized power spectrum estimator can distinguish between fields with or without the fluctuation-background correlation. This cannot be accomplished by the power spectrum with the normalization by an overall mean transmission. Therefore, the self-normalized power spectrum would be useful for the discrimination among models without the uncertainties caused by free (or fitting) parameters.

We analyze the power spectrum of a QSO's Ly_alpha transmitted flux in the discrete wavelet transform (DWT) representation. Although the mean DWT power spectrum is consistent with its counterpart in Fourier representation, the spatial distribution of the local power varies greatly, i.e. the local DWT power spectra show remarkably spiky structures on small scales. To measure these spiky features, we introduce the quantities, roughness, of the local power spectrum, and the correlation between spikes on different scales. We then test the predictions made by the correlation hierarchy model on the roughness and the scale-scale correlations of the local power spectrum. Using the Ly_alpha transmitted flux of the QSO HS1700, we find that the underlying cosmic mass field of the transmitted flux at redshift around z ~ 2.2 can be described by the hierarchical clustering model on physical scales from 2.5 h^{-1} Mpc to few tens h^{-1} kpc in an Einstein-de Sitter universe. However, the non-linear features of the clustering show differences on different scale ranges; 1. On physical scales larger than ~ 1.3 h^{-1} Mpc, the field is almost Gaussian. 2. On scales 1.3 h^{-1} Mpc - 0.3 h^{-1} Mpc, the field is consistent with the correlation hierarchy with a constant value for the coefficient Q_4. 3. On scales less than 300 h^{-1} kpc, the field is no longer Gaussian, but essentially intermittent. In this case, the field can still be fitted by the correlation hierarchy, but the coefficient, Q_4, should be scale-dependent. These three points are strongly supported by the following result: the scale dependencies of Q_4 given by two statistically independent measures, i.e. Q_4^R by the roughness and Q_4^C by scale-scale correlation, are the same in the entire scale range considered.

The intermittent behavior of the space-scale distribution of Ly$\alpha$ transmitted flux of QSO HS1700+64 has been analyzed via a discrete wavelet transform. We found that there are strong indications of intermittency on scales down to about 10 $h^{-1}$ kpc. These are: 1.) the probability distribution function of the local fluctuations of the flux is significantly long-tailed on small scales, and 2.) the local power spectrum of the flux shows prominent spiky structures on small scales. Moreover, the local power spectrum averaged on regions with different sizes shows similar spiky structures. Therefore, the random mass density field traced by the Ly$\alpha$ forests is rougher on smaller scales, consistent with singular clustering.