Spectra of high redshift QSOs show deep Gunn-Peterson absorptions on the blue sides of the \Lya emissions lines. They can be decomposed into components called \Lya leaks, defined to be emissive regions in complementary to otherwise zero-fluxed absorption gaps. Just like \Lya absorption forests at low redshifts, \Lya leaks are both easy to find in observations and containing rich sets of statistical properties that can be used to study the early evolution of the IGM. Among all properties of a leak profile, we investigate its equivalent width in this paper, since it is weakly affected by instrumental resolution and noise. Using 10 Keck QSO spectra at $z\sim6$, we have measured the number density distribution function $n(W,z)$, defined to be the number of leaks per equivalent width $W$ and per redshift $z$, in the redshift range $5.4 - 6.0$. These new observational statistics, in both the differential and cumulative forms, fit well to hydro numerical simulations of uniform ionizing background in the $\Lambda$CDM cosmology. In this model, Ly $\alpha$ leaks are mainly due to low density voids. It supports the early studies that the IGM at $z\simeq6$ would still be in a highly ionized state with neutral hydrogen fraction $\simeq 10^{-4}$. Measurements of $n(W,z)$ at $z>6$ would be effective to probe the reionization of the IGM.

Ly$\alpha$ absorption spectra of QSOs at redshifts $z\simeq6$ show complete Gunn-Peterson absorption troughs (dark gaps) separated by tiny leaks. The dark gaps are from the intergalactic medium (IGM) where the density of neutral hydrogen are high enough to produce almost saturated absorptions, however, where the transmitted leaks come from is still unclear so far. We demonstrate that leaking can originate from the lowest density voids in the IGM as well as the ionized patches around ionizing sources using semi-analytical simulations. If leaks were produced in lowest density voids, the IGM might already be highly ionized, and the ionizing background should be almost uniform; in contrast, if leaks come from ionized patches, the neutral fraction of IGM would be still high, and the ionizing background is significantly inhomogeneous. Therefore, the origin of leaking is crucial to determining the epoch of inhomogeneous-to-uniform transition of the the ionizing photon background. We show that the origin could be studied with the statistical features of leaks. Actually, Ly$\alpha$ leaks can be well defined and described by the equivalent width $W$ and the full width of half area $W_{\rm H}$, both of which are less contaminated by instrumental resolution and noise. It is found that the distribution of $W$ and $W_{\rm H}$ of Ly$\alpha$ leaks are sensitive to the modeling of the ionizing background. We consider four representative reionization models. It is concluded that the leak statistics provides an effective tool to probe the evolutionary history of reionization at $z\simeq5-6.5$. Similar statistics would also be applicable to the reionization of He II at $z \simeq 3$(Abridged)

We study the Gunn-Peterson effect of the photo-ionized intergalactic medium(IGM) in the redshift range 5< z <6.4 using semi-analytic simulations based on the lognormal model. Assuming a rapidly evolved and spatially uniform ionizing background, the simulation can produce all the observed abnormal statistical features near redshift z ~ 6. They include: 1) rapidly increase of absorption depths; 2) large scatter in the optical depths; 3) long-tailed distributions of transmitted flux and 4) long dark gaps in spectra. These abnormal features are mainly due to rare events, which correspond to the long-tailed probability distribution of the IGM density field, and therefore, they may not imply significantly spatial fluctuations in the UV ionizing background at z ~ 6.

We investigate the 21 cm emission from the epoch between reionization z_r and Gunn-Peterson transparency z_{gp}. According to the lognormal model of the thermal history around reionization, hydrogen clouds in z_r > z > z_{gp} are hot and a predominant part of baryonic gas is ionized, but still opaque to Ly-alpha photons. Therefore, 21 cm emission is a distinctive characteristic of this epoch. We show that the 21 cm emission comes from both uncollapsed and collapsing hydrogen clouds. The spatial distribution of the brightness temperature excess $\delta T_b$ is highly non-Gaussian. It consists of spikes with high $\delta T_b$ and a low $\delta T_b$ area between the spikes. The field has the following statistical features: (1) the one-point distributions of $\delta T_b$ are described approximately by power-law tailed probability distribution functions; (2) the nth-order moment of $\delta T_b$ is increasing much faster with n than that of a Gaussian field, but slower than that of a lognormal field; (3) the scale-scale correlation of the $\delta T_b$ field is significant for all scales larger than the Jeans length of the gas. These features would be useful for distinguishing the 21 cm emission of the early clustering from the noise of foreground contamination.

We study the evolution of baryonic gas before the reionization in the lognormal (LN) model of cosmic clustering. We show that the thermal history of the universe around the reionization can roughly be divided into three epochs: 1) cold dark age $z>z_r$, in which baryon gas is neutral, and opaque to Ly$\alpha$ photons; 2) hot dark age $z_r > z> z_{gp}$, in which a predominant part of baryon gas is ionized and hot, but it is still opaque to Ly$\alpha$ photons; 3) bright age $z<z_{gp}$, in which the universe is ionized highly enough to be transparent to Ly$\alpha$ photons. In the flat cold dark matter cosmological models given by WMAP and COBE, the difference of the two redshifts $z_r - z_{gp}$ is found to be as large as $\sim 10$ with $z_r\sim 17$ and $z_{gp}\sim 7$. This reionization history naturally yields a high optical depth to the CMB $\tau_e \simeq 0.12 - 0.19$ observed by the TE polarization of the WMAP, and a low redshift $z_{gp}$ of the appearance of the Ly$\alpha$ Gunn-Peterson trough $z_{gp} \simeq 6 - 8$ in QSO's absorption spectra. The reason why the universe stays long in an ionized, yet Ly$\alpha$ opaque, stage is because the first photo-ionization heats the intergalactic gas effectively and has balanced the gravitational clustering a long period of time. Therefore, the result of a high $\tau_e$ and low $z_{gp}$ is a common feature of all the models considered. Besides the cosmological parameters, the only free parameter we used in the calculation is $N_{ion}$, the mean ionization photons produced by each baryon in collapsed objects. We take it to be 40 - 80 in the calculation.

We study the baryonic gas clouds (the IGM) in the universe before the reionization with the lognormal model which is shown to be dynamcially legitimate in describing the fluctuation evolution in quasilinear as well as nonlinear regimes in recent years. The probability distribution function of the mass field in the LN model is long tailed and so plays an important role in rare events, such as the formation of the first generation of baryonic objects. We calculate density and velocity distributions of the IGM at very high spatial resolutions, and simulate the distributions at resolution of 0.15 kpc from z=7 to 15 in the LCDM cosmological model. We performed a statistics of the hydrogen clouds including column densities, clumping factors, sizes, masses, and spatial number density etc. One of our goals is to identify which hydrogen clouds are going to collapse. By inspecting the mass density profile and the velocity profile of clouds, we found that the velocity outflow significantly postpones the collapsing process in less massive clouds, in spite of their masses are larger than the Jeans mass. Consequently, only massive (> 10^5 M_sun) clouds can form objects at higher redshift, and less massive (10^4-10^5) collapsed objects are formed later. For example, although the mass fraction in clouds with sizes larger than the Jeans length is already larger than 1 at z=15, there is only a tiny fraction of mass (10^{-8}) in the clouds which are collapsed at that time. If all the ionizing photons, and the 10^{-2} metallicity observed at low redshift are produced by the first 1% mass of collapsed baryonic clouds, the majority of those first generation objects would not happen until z=10.

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.

Combining the observations of spatial abundance and two-point correlation function of QSOs, we can effectively set constraints on models of cosmic structure formation. Both the abundance of gravitationally confined halos and their two-point correlation functions can be calculated in the conventional Press-Schechter formalism. We apply this method to examine the properties of possible host halos of QSOs in three popular models: the standard cold dark matter (SCDM) model, the low density flat cold dark matter (LCDM) model and the cold-plus-hot dark matter (CHDM) model. The LCDM and CHDM models are normalized to the COBE-DMR observations, and the SCDM is normalized to $\sigma = 0.58$. We find that the SCDM and LCDM models can pass the abundance-plus-correlation test for QSOs. However, the CHDM are difficult to produce host halos to fit with the number density of high redshift QSOs and their clustering on large scales (10 h^{-1} Mpc) simultaneously. We studied various mechanisms, originated both gravitationally and non-gravitationally, which may lead to a biasing of the halo clustering. We conclude that these effects are too weak in order to release the trouble of the CHDM models.

We have performed a detailed statistical study of the evolution of structure in a photoionized intergalactic medium (IGM) using analytical simulations to extend the calculation into the mildly non-linear density regime found to prevail at z = 3. Our work is based on a simple fundamental conjecture: that the probability distribution function of the density of baryonic diffuse matter in the universe is described by a lognormal (LN) random field. The LN field has several attractive features and follows plausibly from the assumption of initial linear Gaussian density and velocity fluctuations at arbitrarily early times. Starting with a suitably normalized power spectrum of primordial fluc- tuations in a universe dominated by cold dark matter (CDM), we compute the behavior of the baryonic matter, which moves slowly toward minima in the dark matter potential on scales larger than the Jeans length. We have computed two models that succeed in matching observations. One is a non-standard CDM model with Omega=1, h=0.5 and \Gamma=0.3, and the other is a low density flat model with a cosmological constant(LCDM), with Omega=0.4, Omega_Lambda=0.6 and h=.65. In both models, the variance of the density distribution function grows with time, reaching unity at about z=4, where the simulation yields spectra that closely resemble the Ly-alpha forest absorption seen in the spectra of high z quasars. The calculations also successfully predict the observed properties of the Ly-alpha forest clouds and their evolution from z=4 down to at least z=2, assuming a constant intensity for the metagalactic UV background over this redshift range. However, in our model the forest is not due to discrete clouds, but rather to fluctuations in a continuous intergalactic medium. (This is an abreviated abstract; the complete abstract is included with the manuscript.)

We have developed a method for calculating the two-point correlation function of nonlinearly evolved mass and collapsed halos in the Press-Schechter formalism. The nonlinear gravitational interaction is treated as the sum of various individual spherical top-hat clustering. Because no collapsed halo of mass M can exist in initial regions (or top-hat spheres) of mass less than M, the bias that massive halos have stronger correlation than the background mass can be naturally introduced. We apply this method to derive constraints on popular dark-matter models from the spatial number density and the correlation function of C IV absorption systems in QSO spectra. Considering C IV systems should behosted by collapsed halos, one can obtain an upper limit to the threshold mass of the collapsed halos by requiring their number density to be larger than that of observed C IV systems. On the other hand, in order to explain the observed clustering of C IV systems, a lower limit to the threshold mass will be set for the hosting halos. The cold dark matter model, the Lamda CDM model and the cold-plus-hot dark matter model are tested using the C IV observations.

Our goal in this paper is to test some popular dark matter models by Ly-alpha forest in QSO spectra. Recent observations of the size and velocity of Ly-alpha forest clouds have indicated that the Ly-alpha absorption is probably not given by collapsed objects, but pre-collapsed regions in the baryonic density field. Therefore, a linear approximation description would be able to provide valuable information. We developed a technique to simulate Ly-alpha forest as the absorption of such pre-collapsed regions under linear approximation regime. The simulated Ly-alpha forests in models of the standard cold dark matter (SCDM), the cold plus hot dark matter (CHDM), and the low-density flat cold dark matter (LCDM) have been confronted with observational features, including 1) the number density of Ly-alpha lines and its dependencies on redshift and equivalent width; 2) the distribution of equivalent widths and its redshift dependence; 3) clustering; and 4) the Gunn-Peterson effect. The "standard" CHDM model, i.e. 60% cold and 30% hot dark matters and 10\% baryons, is found to be difficult to pass the Ly-alpha forest test, probably because it produces structures too late and favors to form structures on large scales instead of small scale objects like Ly-alpha clouds. Within a reasonable range of J_nu, the UV background radiation at high redshift, and delta_th, the threshold of the onset of gravitational collapse of the baryonic matter, the LCDM model is consistent with observational data in all above-mentioned four aspects. The model of SCDM can also fit with observation, but it requires a smaller J_nu and a higher delta_th. This suggests that whether a significant part of the Ly-alpha forest lines is located in the halos of collapsed objects would be crucial to the success of SCDM.

(\rl = Richardson-Lucy) We propose a simulation-based bootstrap method to access global significance levels of deconvolution models in the \rl and other iterative restoration algorithms that converge locally. These significance levels allow one to check at each iterative step how good the model is and when iterations can be stopped. Adding more iterations in the deconvolution improves the fitting but is very slow at later time; while too much entropy or smoothness will be lost in the models. A good deconvolution model should firstly have a significance level as high as possible ($\ge$ 20\%), and secondly, be as smooth as possible. We have used two examples to illustrate how such models can be derived in practice. We point out that maximizing the sum of the likelihood of fitting and {\em a priori} entropy does not guarantee an acceptable significance level for the resulting model. If one's {\em a priori} knowledge is too poor, the model may not be able to fit the data at a reasonable significance level. Instead, a maximum-entropy-like iterative restoration algorithm can be performed later by acquiring {\em a priori} knowledge from the \rl restoration. However, this is necessary only when it does increase the levels significantly.