We investigate spatially resolved specific star formation rate (SSFR) in the inner $\sim$ 40 pc for a nearby Seyfert 2 galaxy, M51 (NGC 5194) by analyzing spectra obtained with the \emph{Hubble Space Telescope (HST)} Space Telescope Imaging Spectrograph (STIS). We present 24 radial spectra measured along the STIS long slit in M51, extending $\sim 1\arcsec$ from the nucleus (i.e., -41.5 pc to 39.4 pc). By the simple stellar population synthesis, the stellar contributions in these radial optical spectra are modeled. Excluding some regions with zero young flux fraction near the center (from -6 pc to 2 pc), we find that the mean flux fraction of young stellar populations (younger than 24.5 Myr) is about 9 \%, the mean mass fraction is about 0.09\%. The young stellar populations are not required in the center inner $\sim$ 8 pc in M51, suggesting a possible SSFR suppression in the circumnuclear region ($\sim$ 10 pc) from the feedback of active galactic nuclei (AGNs). The radial distribution of SSFR in M51 is not symmetrical with respect to the long slit in STIS. This unsymmetrical SSFR distribution is possibly due to the unsymmetrical AGN feedback in M51, which is related to its jet.

We use a suite of cosmological hydrodynamic simulations, run by two fixed grid codes, to investigate the properties of solenoidal and dilatational motions of the intergalactic medium (IGM), and the impact of numerical viscosity on turbulence in a LCDM universe. The codes differ only in the spatial difference discretization. We find that (1) The vortical motion grows rapidly since $z=2$, and reaches $\sim 10 km/s -90 km/s$ at $z=0$. Meanwhile, the small-scale compressive ratio $r_{CS}$ drops from 0.84 to 0.47, indicating comparable vortical and compressive motions at present. (2) Power spectra of the solenoidal velocity possess two regimes, $\propto k^{-0.89}$ and $\propto k^{-2.02}$, while the total and dilatational velocity follow the scaling $k^{-1.88}$ and $k^{-2.20}$ respectively in the turbulent range. The IGM turbulence may contain two distinct phases, the supersonic and post-supersonic phases. (3) The non-thermal pressure support, measured by the vortical kinetic energy, is comparable with the thermal pressure for $\rho_b \simeq 10-100$, or $T <10^{5.5} K$ at $z=0.0$. The deviation of the baryon fraction from the cosmic mean shows a preliminary positive correlation with the turbulence pressure support. (4) A relatively higher numerical viscosity would dissipate both the compressive and vortical motions of the IGM into thermal energy more effectively, resulting in less developed vorticity, remarkably shortened inertial range, and leading to non-negligible uncertainty in the thermal history of gas accretion. Shocks in regions outside of clusters are significantly suppressed by numerical viscosity since $z=2$, which may directly cause the different levels of turbulence between two codes.

Using the properties of SDSS DR7 QSOs catalog from Shen et al., the Baldwin effect, its slope evolution, the underlying drive for a large sample of 35019 QSOs with reliable spectral analysis are investigated. We find that the Baldwin effect exists in this large QSOs sample, which is almost the same in 11 different redshift bins, up to $z\sim 5$. The slope is -0.238 by the BCES (\civ\ EW depends on the continuum), -0.787 by the BCES bisector. For 11 redshift-bins, there is an increasing of the Baldwin effect slope from $z\sim1.5$ to $z\sim2.0$. From $z\sim2.0$ to $z\sim5.0$, the slope change is not clear considering their uncertainties or larger redshift bins. There is a strong correlation between the rest-frame \civ\ EW and \civ-based \mbh while the relation between the \civ\ EW and \mgii-based \mbh is very weak. With the correction of \civ-based \mbh from the \civ\ blueshift relative to \mgii, we suggest that this strong correlation is due to the bias of the \civ-based \mbh, with respect to that from the \mgii\ line. Considering the \mgii-based \mbh, a medium strong correlation is found between the \civ\ EW and the Eddington ratio, which implies that the Eddington ratio seems to be a better underlying physical parameter than the central black hole mass.

With Monte Carlo simulation method, we investigate the time dependent behavior of Ly$\alpha$ photon transfer in optically thick medium of the concordance $\Lambda$CDM universe. At high redshift, the Ly$\alpha$ photon escaping from optically thick medium has a time scale as long as the age of the luminous object, or even comparable to the age of the universe. In this case, time-independent, or stationary solutions of the Ly$\alpha$ photon transfer with resonant scattering will overlook important features of the escaped Ly$\alpha$ photons in physical and frequency spaces. More seriously, the expansion of the universe leads to that the time-independent solutions of the Ly$\alpha$ photon transfer may not exist. We show that time-dependent solutions sometimes are essential for understanding the Ly$\alpha$ emission and absorption at high redshifts. For Ly\alpha photons from sources at redshift 1+z=10 and being surrounded by neutral hydrogen IGM of the $\Lambda$CDM universe, the escape coefficient is found to be always less, or much less than one, regardless of the age or life time of the sources. Under such environment, we also find that even when the Ly$\alpha$ photon luminosity of the sources is stable, the mean surface brightness is gradually increasing in the first 10^6 years, and then decreasing with a power law of time, but never approaches a stable, time-independent state. That is, all 1+z=10 sources in a neutral Hubble expanding IGM with Ly$\alpha$ luminosity L have their maximum of mean surface brightness ~ 10^{-21}(L/(10^{43}erg/s)) erg s^{-1} cm^{-2} arcsec^{-2} at the age of about 10^6 years. The time-dependent effects on the red damping wing profile are also addressed.

We investigate the effects of dust on Ly{\alpha} photons emergent from an optically thick medium by solving the integro-differential equation of the radiative transfer of resonant photons. To solve the differential equations numerically we use the Weighted Essentially Non-oscillatory method (WENO). Although the effects of dust on radia-tive transfer is well known, the resonant scattering of Ly{\alpha} photons makes the problem non-trivial. For instance, if the medium has the optical depth of dust absorption and scattering to be {\tau}a>>1, {\tau}>>1, and {\tau}>>{\tau}a, the effective absorption optical depth in a random walk scenario would be equal to \surd {\tau}a({\tau}a+{\tau}). We show, however, that for a resonant scattering at frequency {\nu}0, the effective absorption optical depth would be even larger than {\tau}({\nu}0). If the cross section of dust scattering and absorption is frequency-independent, the double-peaked structure of the frequency profile given by the resonant scattering is basically dust-independent. That is, dust causes neither narrowing nor widening of the width of the double peaked profile. One more result is that the time scales of the Ly{\alpha} photon transfer in the optically thick halo are also basically independent of the dust scattering, even when the scattering is anisotropic. This is because those time scales are mainly determined by the transfer in the frequency space, while dust scattering, either isotropic or anisotropic, does not affect the behavior of the transfer in the frequency space when the cross section of scattering is wavelength-independent. This result does not support the speculation that dust will lead to the smoothing of the brightness distribution of Ly{\alpha} photon source with optical thick halo.

We investigate the possibility of detecting the turbulent state of the IGM with the kinetic Sunyaev-Zel'dovich (kSZ) effect. Being sensitive to the divergence-free component of the momentum field of the IGM, the kSZ effect might be used to probe the vorticity of the turbulent IGM. With cosmological hydrodynamical simulation in the concordance $\Lambda$CDM universe, we find that the structure functions of 2D kSZ maps show strong intermittence, and the intermittent exponents follow a law similar to the She-Leveque scaling formula of fully developed turbulence. We also find that the intermittence is weak in the maps of thermal Sunyaev-Zel'dovich (tSZ) effect. Nevertheless, the superposition of the kSZ and tSZ effects still contain significant intermittence. We conclude that the turbulent behavior of the IGM may be revealed by the observation of SZ effect on angular scales equal to or less than 0.5 arcminute, corresponding to the multipole parameter $l\geq 2 \times10^4$.

The mean mass densities of cosmic dark matter is larger than that of baryonic matter by a factor of about 5 in the $\Lambda$CDM universe. Therefore, the gravity on large scales should be dominant by the distribution of dark matter in the universe. However, a series of observations incontrovertibly show that the velocity and density fields of baryonic matter are decoupling from underlying dark matter field. This paper shows our attemps to unveil the physics behind this puzzle. In linear approximation, the dynamics of the baryon fluid is completely governed by the gravity of the dark matter. Consequently, the mass density field of baryon matter $\rho_b({\bf r},t)$ will be proportional to that of dark matter $\rho_{\rm dm}({\bf r},t)$, even though they are different from each other initially. In weak and moderate nonlinear regime, the dynamics of the baryon fluid can be sketched by Burgers equation. A basic feature of the Burgers dynamics is to yield shocks. When the Reynolds number is large, the Burgers fluid will be in the state of Burgers turbulence, which consists of shocks and complex structures. On the other hand, the collisionless dark matter may not show such shock, but a multivalued velocity field. Therefore, the weak and moderate nonlinear evolution leads to the IGM-dark matter deviation. Yet, the velocity field of Burgers fluid is still irrotational, as gravity is curl-free. In fully nonlinear regime, the vorticity of velocity field developed, and the cosmic baryonic fluid will no longer be potential, as the dynamics of vorticity is independent of gravity and can be self maintained by the nonlinearity of hydrodynamics. In this case, the cosmic baryon fluid is in the state of fully developed turbulence, which is statistically and dynamically decoupling from dark matter. This scenario provides a mechanism of cohenent explanation of observations.

We investigate the non-Gaussian features of the IGM at redshift $z\sim 5 - 6$ using Ly$\alpha$ transmitted flux of quasar absorption spectra and cosmological hydrodynamic simulation of the concordance $\Lambda$CDM universe. We show that the neutral hydrogen mass density field and Ly$\alpha$ transmitted flux fluctuations possess all the non-Gaussian features predicted by the log-Poisson hierarchy, which depends only on two dimensionless parameters $\beta$ and $\gamma$, describing, respectively, the intermittence and singularity of the random fields. We find that the non-Gaussianity of the Ly$\alpha$ transmitted flux of quasars from $z=4.9$ to $z=6.3$ can be well reconstructed by the hydrodynamical simulation samples. Although the Gunn-Peterson optical depth and its variance underwent a significant evolution in the redshift range of $5 - 6$, the intermittency measured by $\beta$ is almost redshift-independent in this range. More interesting, the intermittency of quasar's absorption spectra on physical scales $0.1-1$ h$^{-1}$Mpc in redshift $5 - 6$ are found to be about the same as that on physical scales $1-10$ h$^{-1}$Mpc at redshifts $2 - 4$. Considering the Jeans length is less than 0.1 h$^{-1}$Mpc at $z\sim 5$, and $1$ h$^{-1}$Mpc at $z\sim 2$, these results imply that the nonlinear evolution in high and low redshifts will lead the cosmic baryon fluid to a state similar to fully developed turbulence. The log-Poisson high order behavior of current high redshift data of quasar's spectrum can be explained by uniform UV background in the redshift range considered. We also studied the log-Poisson non-Gaussianity by considering inhomogeneous background. With several simplified models of inhomogeneous background, we found the effect of the inhomogeneous background on the log-Poisson non-Gaussianity is not larger than 1-sigma.

In this paper we have analyzed the temporal and spectral behavior of 52 Fast Rise and Exponential Decay (FRED) pulses in 48 long-duration gamma-ray bursts (GRBs) observed by the CGRO/BATSE, using a pulse model with two shape parameters and the Band model with three shape parameters, respectively. It is found that these FRED pulses are distinguished both temporally and spectrally from those in long-lag pulses. Different from these long-lag pulses only one parameter pair indicates an evident correlation among the five parameters, which suggests that at least $\sim$4 parameters are needed to model burst temporal and spectral behavior. In addition, our studies reveal that these FRED pulses have correlated properties: (i) long-duration pulses have harder spectra and are less luminous than short-duration pulses; (ii) the more asymmetric the pulses are the steeper the evolutionary curves of the peak energy ($E_{p}$) in the $\nu f_{\nu}$ spectrum within pulse decay phase are. Our statistical results give some constrains on the current GRB models.

With a state-of-the-art numerical method for solving the integral-differential equation of radiative transfer, we investigate the flux of the Ly$\alpha$ photon $\nu_0$ emergent from an optically thick halo containing a central light source. Our focus is on the time-dependent effects of the resonant scattering. We first show that the frequency distribution of photons in the halo are quickly approaching to a locally thermalized state around the resonant frequency, even when the mean intensity of the radiation is highly time-dependent. Since initial conditions are forgotten during the thermalization, some features of the flux, such as the two peak structure of its profile, actually are independent of the intrinsic width and time behavior of the central source, if the emergent photons are mainly from photons in the thermalized state. In this case, the difference $|\nu_{\pm}-\nu_0|$, where $\nu_{\pm}$ are the frequencies of the two peaks of the flux, cannot be less than $2$ times of Doppler broadening. We then study the radiative transfer in the case where the light emitted from the central source is a flash. We calculate the light curves of the flux from the halo. It shows that the flux is still a flash. The time duration of the flash for the flux, however, is independent of the original time duration of the light source but depends on the optical depth of the halo. Therefore, the spatial transfer of resonant photons is a diffusion process, even though it is not a purely Brownian diffusion. This property enables an optically thick halo to trap and store thermalized photons around $\nu_0$ for a long time after the cease of the central source emission. The photons trapped in the halo can yield delayed emission, of which the profile also shows typical two peak structure as that from locally thermalized photons. Possible applications of these results are addressed.

We investigate the vorticity of the IGM velocity field on large scales with cosmological hydrodynamic simulation of the concordance model of LCDM. We show that the vorticity field is significantly increasing with time as it can effectively be generated by shocks and complex structures in the IGM. Therefore, the vorticity field is an effective tool to reveal the nonlinear behavior of the IGM, especially the formation and evolution of turbulence in the IGM. We find that the vorticity field does not follow the filaments and sheets structures of underlying dark matter density field and shows highly non- Gaussian and intermittent features. The power spectrum of the vorticity field is used to measure the development of turbulence in Fourier space. We show that the relation between the power spectra of vorticity and velocity fields is perfectly in agreement with the prediction of a fully developed homogeneous and isotropic turbulence from 0.2 to 3 h^{-1} Mpc at z~0. This indicates that cosmic baryonic field is in the state of fully developed turbulence on scales less than about 3 h^{-1} Mpc. The random field of the turbulent fluid yields turbulent pressure to prevent the gravitational collapsing of the IGM. The vorticity and turbulent pressure are strong inside and even outside of high density regions. In IGM regions with 10 times mean overdensity, the turbulent pressure can be on an average equivalent to the thermal pressure of the baryonic gas with a temperature of 10^5 K. The fully developed turbulence would prevent the baryons in the IGM from falling into the gravitational well of dark matter halos. Moreover, turbulent pressure is dynamical and non-thermal, which makes it different from pre-heating mechanism as it does not affect the thermal state and ionizing process of hydrogen in the IGM.

We develop an algorithm of separating the $E$ and $B$ modes of the CMB polarization from the noisy and discretized maps of Stokes parameter $Q$ and $U$ in a finite area. A key step of the algorithm is to take a wavelet-Galerkin discretization of the differential relation between the $E$, $B$ and $Q$, $U$ fields. This discretization allows derivative operator to be represented by a matrix, which is exactly diagonal in scale space, and narrowly banded in spatial space. We show that the effect of boundary can be eliminated by dropping a few DWT modes located on or nearby the boundary. This method reveals that the derivative operators will cause large errors in the $E$ and $B$ power spectra on small scales if the $Q$ and $U$ maps contain Gaussian noise. It also reveals that if the $Q$ and $U$ maps are random, these fields lead to the mixing of the $E$ and $B$ modes. Consequently, the $B$ mode will be contaminated if the powers of $E$ modes are much larger than that of $B$ modes. Nevertheless, numerical tests show that the power spectra of both $E$ and $B$ on scales larger than the finest scale by a factor of 4 and higher can reasonably be recovered, even when the power ratio of $E$- to $B$-modes is as large as about 10$^2$, and the signal-to-noise ratio is equal to 10 and higher. This is because the Galerkin discretization is free of false correlations, and keeps the contamination under control. As wavelet variables contain information of both spatial and scale spaces, the developed method is also effective to recover the spatial structures of the $E$ and $B$ mode fields.

The 21 cm emission and absorption from gaseous halos around the first generation of star depend on the Wouthuysen-Field (W-F) coupling, which relates the spin temperature with the kinetic temperature of hydrogen gas via the resonant scattering between Lyman alpha photons and neutral hydrogen. Although the center object generally is a strong source of these photons, the transfer of these photons in the 21 cm region is inefficient, as the optical depth of the photons is large. Consequently, these photons from the source may not be able to transfer to the entire 21 cm region timely to provide the W-F coupling. This problem is important because the lifetime of first stars generally is short. The problem is investigated with numerical solution of the integro-differential equation, which describes the kinetics of these resonant photons in both physical and frequency spaces. We show that the photon transfer process in the physical space is actually coupled to that in the frequency space. Firstly diffusion in the frequency space provides a shortcut for the diffusion in the physical space. It makes the mean time for the escape of the resonant photon in optical depth \tau media roughly proportional to the optical depth \tau, not \tau^2. Secondly the resonant scattering is effective in bouncing photons with a frequency which is not equal to initial frequency back to the initial frequency. This process can restore initial frequency photons and establish the local Boltzmann distribution of the photon spectrum around the initial frequency. Therefore, the mechanism of 'escape via shortcut' plus 'bounce back' enables W-F coupling to be properly realized in the 21 cm region around first stars. This mechanism also works for photons injected into the 21 cm region by redshift.

Employing two samples containing of 56 and 59 well-separated FRED (fast rise and exponential decay) gamma-ray burst (GRB) pulses whose spectra are fitted by the Band spectrum and Compton model, respectively, we have investigated the evolutionary slope of $E_{p}$ (where $E_{p}$ is the peak energy in the $\nu F\nu$ spectrum) with time during the pulse decay phase. The bursts in the samples were observed by the Burst and Transient Source Experiment (BATSE) on the Compton Gamma-Ray Observatory. We first test the $E_{p}$ evolutionary slope during the pulse decay phase predicted by Lu et al. (2007) based on the model of highly symmetric expanding fireballs in which the curvature effect of the expanding fireball surface is the key factor concerned. It is found that the evolutionary slopes are normally distributed for both samples and concentrated around the values of 0.73 and 0.76 for Band and Compton model, respectively, which is in good agreement with the theoretical expectation of Lu et al. (2007). However, the inconsistence with their results is that the intrinsic spectra of most of bursts may bear the Comptonized or thermal synchrotron spectrum, rather than the Band spectrum. The relationships between the evolutionary slope and the spectral parameters are also checked. We show the slope is correlated with $E_{p}$ of time-integrated spectra as well as the photon flux but anticorrelated with the lower energy index $\alpha$. In addition, a correlation between the slope and the intrinsic $E_{p}$ derived by using the pseudo-redshift is also identified. The mechanisms of these correlations are unclear currently and the theoretical interpretations are required.

We develop a numerical solver for the integral-differential equations, which describes the radiative transfer of photon distribution in the frequency space with resonant scattering of Lyalpha photons by hydrogen gas in the early universe. The time-dependent solutions of this equation is crucial to the estimation of the effect of the Wouthuysen-Field (WF) coupling in relation to the 21 cm emission and absorption at the epoch of reionization. The resonant scattering leads to the photon distribution in the frequency space to be piecewise smooth containing sharp changes. The weighted essentially nonoscillatory (WENO) scheme is suitable to handle this problem, as this algorithm has been found to be highly stable and robust for solving Boltzmann equation. We test this numerical solver against analytic solutions of the evolution of the photon distribution in rest background, analytic solution in expanding background without resonant scattering and formation of local Boltzmann distribution around the resonant frequency with the temperature same as that of atom for recoil. We find that evolution of photon distribution undergoes three phases; profile is similar to the initial one, a flat plateau (without recoil) or local Boltzmann distribution (with recoil) forms around the resonant frequency, and finally the distribution around the resonant frequency is saturated when the photons from the source is balanced by the redshift of the expansion. This result indicates that the onset of the W-F coupling should not be determined by the third phase, but by the time scale of the second phase. We found that the time scale of the W-F coupling is equal to about a few hundreds of the mean free flight time of photons with resonant frequency, and is independent of the Sobolev parameter if this parameter is much less than 1.

We study the Wouthuysen-Field coupling at early universe with numerical solutions of the integrodifferential equation describing the kinetics of photons undergoing resonant scattering. The numerical solver is developed based on the weighted essentially non-oscillatory (WENO) scheme for the Boltzmann-like integrodifferential equation. We focus on the time evolution of the Wouthuysen-Field (W-F) coupling in relation to the 21 cm emission and absorption at the epoch of reionization. We show that a local Boltzmann distribution will be formed if photons with frequency \sim \nu_0 have undergone a ten thousand or more times of scattering, which corresponds to the order of 10^3 yrs for neutral hydrogen density of the concordance \Lambda CDM model. The time evolution of the shape and width of the local Boltzmann distribution actually doesn't dependent on the details of atomic recoil, photon sources, or initial conditions very much. However, the intensity of photon flux at the local Boltzmann distribution is substantially time-dependent. The time scale of approaching the saturated intensity can be as long as 10^5-10^6 yrs for typical parameters of the \Lambda CDM model. The intensity of the local Boltzmann distribution at time less than 10^5 yrs is significantly lower than that of the saturation state. Therefore, it may not be always reasonable to assume that the deviation of the spin temperature of 21 cm energy states from cosmic background temperature is mainly due to the W-F coupling if first stars or their emission/absorption regions evolved with a time scale equal to or less than Myrs.

Employing a sample presented by Kaneko et al. (2006) and Kocevski et al. (2003), we select 42 individual tracking pulses (here we defined tracking as the cases in which the hardness follows the same pattern as the flux or count rate time profile) within 36 Gamma-ray Bursts (GRBs) containing 527 time-resolved spectra and investigate the spectral hardness, $E_{peak}$ (where $E_{peak}$ is the maximum of the $\nu F_{\nu}$ spectrum), evolutionary characteristics. The evolution of these pulses follow soft-to-hard-to-soft (the phase of soft-to-hard and hard-to-soft are denoted by rise phase and decay phase, respectively) with time. It is found that the overall characteristics of $E_{peak}$ of our selected sample are: 1) the $E_{peak}$ evolution in the rise phase always start on the high state (the values of $E_{peak}$ are always higher than 50 keV); 2) the spectra of rise phase clearly start at higher energy (the median of $E_{peak}$ are about 300 keV), whereas the spectra of decay phase end at much lower energy (the median of $E_{peak}$ are about 200 keV); 3) the spectra of rise phase are harder than that of the decay phase and the duration of rise phase are much shorter than that of decay phase as well. In other words, for a complete pulse the initial $E_{peak}$ is higher than the final $E_{peak}$ and the duration of initial phase (rise phase) are much shorter than the final phase (decay phase). This results are in good agreement with the predictions of Lu et al. (2007) and current popular view on the production of GRBs. We argue that the spectral evolution of tracking pulses may be relate to both of kinematic and dynamic process even if we currently can not provide further evidences to distinguish which one is dominant. Moreover, our statistical results give some witnesses to constrain the current GRB model.

we study, in this paper, the non-Gaussian features of the mass density field of neutral hydrogen fluid and the Ly-alpha transmitted flux of QSO absorption spectrum from the point-of-view of self-similar log-Poisson hierarchy. It has been shown recently that, in the scale range from the onset of nonlinear evolution to dissipation, the velocity and mass density fields of cosmic baryon fluid are extremely well described by the She-Leveque's scaling formula, which is due to the log-Poisson hierarchical cascade. Since the mass density ratio between ionized hydrogen to total hydrogen is not uniform in space, the mass density field of neutral hydrogen component is not given by a similar mapping of total baryon fluid. Nevertheless, we show, with hydrodynamic simulation samples of the concordance $\Lambda$CDM universe, that the mass density field of neutral hydrogen, is also well described by the log-Poisson hierarchy. We then investigate the field of Ly$\alpha$ transmitted flux of QSO absorption spectrum. Due to redshift distortion, Ly$\alpha$ transmitted flux fluctuations are no longer to show all features of the log-Poisson hierarchy. However, some non-Gaussian features predicted by the log-Poisson hierarchy are not affected by the redshift distortion. We test these predictions with the high resolution and high S/N data of quasars Ly$\alpha$ absorption spectra. All results given by real data, including $\beta$-hierarchy, high order moments and scale-scale correlation, are found to be well consistent with the log-Poisson hierarchy. We compare the log-Poisson hierarchy with the popular log-normal model of the Ly$\alpha$ transmitted flux. The later is found to yield too strong non-Gaussianity at high orders, while the log-Poisson hierarchy is in agreement with observed data.

It has been revealed recently that, in the scale free range, i.e. from the scale of the onset of nonlinear evolution to the scale of dissipation, the velocity and mass density fields of cosmic baryon fluid are extremely well described by the self-similar log-Poisson hierarchy. As a consequence of this evolution, the relations among various physical quantities of cosmic baryon fluid should be scale invariant, if the physical quantities are measured in cells on scales larger than the dissipation scale, regardless the baryon fluid is in virialized dark halo, or in pre-virialized state. We examine this property with the relation between the Compton parameter of the thermal Sunyaev-Zel'dovich effect, $y(r)$, and X-ray luminosity, $L_{\rm x}(r)$, where $r$ being the scale of regions in which $y$ and $L_{\rm x}$ are measured. According to the self-similar hierarchical scenario of nonlinear evolution, one should expect that 1.) in the $y(r)$-$L_x(r)$ relation, $y(r)=10^{A(r)}[L_{\rm x}(r)]^{\alpha(r)}$, the coefficients $A(r)$ and $\alpha(r)$ are scale-invariant; 2.) The relation $y(r)=10^{A(r)}[L_{\rm x}(r)]^{\alpha(r)}$ given by cells containing collapsed objects is also available for cells without collapsed objects, only if $r$ is larger than the dissipation scale. These two predictions are well established with a scale decomposition analysis of observed data, and a comparison of observed $y(r)$-$L_x(r)$ relation with hydrodynamic simulation samples. The implication of this result on the characteristic scales of non-gravitational heating is also addressed.

This paper has been withdrawn

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)

In the nonlinear regime of cosmic clustering, the mass density field of the cosmic baryon fluid is highly non-Gaussian. It shows different dynamical behavior from collisionless dark matter. Nevertheless, the evolved field of baryon fluid is scale-covariant in the range from the Jeans length to a few ten h^{-1} Mpc, in which the dynamical equations and initial perturbations are scale free. We show that in the scale-free range, the non-Gaussian features of the cosmic baryon fluid, governed by the Navier-Stokes equation in an expanding universe, can be well described by a log-Poisson hierarchical cascade. The log-Poisson scheme is a random multiplicative process (RMP), which causes non-Gaussianity and intermittency even when the original field is Gaussian. The log-Poisson RMP contains two dimensionless parameters: $\beta$ for the intermittency and $\gamma$ for the most singular structure. All the predictions given by the log-Poisson RMP model, including the hierarchical relation, the order dependence of the intermittent exponent, the moments, and the scale-scale correlation, are in good agreement with the results given by hydrodynamic simulations of the standard cold dark matter model. The intermittent parameter $\beta$ decreases slightly at low redshift and indicates that the density field of baryon fluid contains more singular structures at lower redshifts. The applicability of the model is addressed.

We investigate the volume growth of ionized regions around UV photon sources with the WENO algorithm, which is an effective solver of photon kinetics in the phase space described by the radiative transfer equation. We show that the volume growth rate, either of isolated ionized regions or of clustered regions in merging, generally consists of three phases: fast or relativistic growth phase at the early stage, slow growth phase at the later stage, and a transition phase between the fast and slow phases. We also show that the volume growth of ionized regions around clustered sources with intensity $\dot{E}_i$ ($i=1, 2, ...$) would have the same behavior as a single source with intensity $\dot{E}=\sum_i\dot{E}_i$, if all the distances between nearest neighbor sources $i$ and $j$ are smaller than $c(t^i_c+t^j_c)$, $t^i_c$ being the time scale $t_c$ of source $i$. Therefore, a tightly clustered UV photon sources would lead to a slow growth of ionized volume. This effect would be important for studying the redshift-dependence of 21cm signals from the reionization epoch.

We investigate the 21 cm signals from the UV ionizing sources in the reionization epoch. The formation and evolution of 21 cm emission and absorption regions depend essentially on the kinetics of photons in the physical and frequency spaces. To solve the radiative transfer equation, we use the WENO algorithm, which is effective to capture the sharp ionization profile and the cut-off at the front of light $(r=ct)$ and to handle the small fraction of neutral hydrogen and helium in the ionized sphere. We show that a spherical shell of 21 cm emission and absorption will develop around a point source once the speed of the ionization front (I-front) is significantly lower than the speed of light. The 21 cm shell extends from the I-front to the front of light; its inner part is the emission region and its outer part is the absorption region. The 21 cm emission region depends strongly on the intensity, frequency-spectrum and life-time of the UV ionizing source. For a source of short life-time, no 21 cm emission region can be formed if the source dies out before the I-front speed is significantly lower than the speed of light. Yet, a 21 cm absorption region can form and develop even after the emission of the source ceases.

We estimate the power spectrum of SZ(Sunyaev-Zel'dovich)-effect-induced temperature fluctuations on sub-degree scales by using the cross correlation between the three-year WMAP maps and 2MASS galaxy distribution. We produced the SZ effect maps by hydrodynamic simulation samples of the $\Lambda$CDM model, and show that the SZ effect temperature fluctuations are highly non-Gaussian. The PDF of the temperature fluctuations has a long tail. More than 70% power of the SZ effect temperature fluctuations attributes to top $\sim 1%$ wavelet modes (long tail events). On the other hand, the CMB temperature fluctuations basically are Gaussian. Although the mean power of CMB temperature fluctuations on sub-degree scales is much higher than that of SZ effect map, the SZ effect temperature fluctuations associated with top 2MASS clusters is comparable to the power of CMB temperature fluctuations on the same scales. Thus, from noisy WMAP maps, one can have a proper estimation of the SZ effect power at the positions of the top 2MASS clusters. The power spectrum given by these top wavelet modes is useful to constrain the parameter of density fluctuations amplitude $\sigma_8$. We find that the power spectrum of these top wavelet modes of SZ effect on sub-degree scales basically is consistent with the simulation maps produced with $\sigma_8=0.84$. The simulation samples of $\sigma_8=0.74$ show, however, significant deviation from detected SZ power spectrum. It can be ruled out with confidence level 99% if all other cosmological parameters are the same as that given by the three-year WMAP results.

We develop a numerical solver for radiative transfer problems based on the weighted essentially nonoscillatory (WENO) scheme modified with anti-diffusive flux corrections, in order to solve the temperature and ionization profiles around a point source of photons in the reionization epoch. Algorithms for such simulation must be able to handle the following two features: 1. the sharp profiles of ionization and temperature at the ionizing front (I-front) and the heating front (T-front), and 2. the fraction of neutral hydrogen within the ionized sphere is extremely small due to the stiffness of the rate equations of atom processes. The WENO scheme can properly handle these two features, as it has been shown to have high order of accuracy and good convergence in capturing discontinuities and complicated structures in fluid as well as to be significantly superior over piecewise smooth solutions containing discontinuities. With this algorithm, we show the time-dependence of the preheated shell around a UV photon source. In the first stage the I-front and T-front are coincident, and propagate with almost the speed of light. In later stage, when the frequency spectrum of UV photons is hardened, the speeds of propagation of the ionizing and heating fronts are both significantly less than the speed of light, and the heating front is always beyond the ionizing front. In the spherical shell between the I- and T-fronts, the IGM is heated, while atoms keep almost neutral. The time scale of the preheated shell evolution is dependent on the intensity of the photon source. We also find that the details of the pre-heated shell and the distribution of neutral hydrogen remained in the ionized sphere are actually sensitive to the parameters used. The WENO algorithm can provide stable and robust solutions to study these details.

The power spectrum of the two-degree Field Galaxy Redshift Survey (2dFGRS) sample is estimated with the discrete wavelet transform (DWT) method. The DWT power spectra within $0.04 <k< 2.3 h$Mpc$^{-1}$ are measured for three volume-limited samples defined in connective absolute magnitude bins $-19 \sim -18$, $-20 \sim -19$ and $-21 \sim -20$. We show that the DWT power spectrum can effectively distinguish $\Lambda$CDM models of $\sigma_8=0.84$ and $\sigma_8=0.74$. We adopt maximum likelihood method to perform three-parameter fitting with bias parameter $b$, pairwise velocity dispersion $\sigma_{pv}$ and redshift distortion parameter $\beta=\Omega_m^{0.6}/b$ to the measured DWT power spectrum. Fitting results denotes that in a $\sigma_8=0.84$ universe the best fitted $\Omega_m$ given by the three samples are consistent in the range $0.28 \sim 0.36$, and the best fitted $\sigma_{pv}$ are $398^{+35}_{-27}$, $475^{+37}_{-29}$ and $550 \pm 20$km/s for the three samples, respectively. However in the model of $\sigma_8=0.74$, our three samples give very different values of $\Omega_m$. We repeat the fitting by using empirical formula of redshift distortion. The result of the model of low $\sigma_8$ is still poor, especially, one of the best value $\sigma_{pv}$ is as large as $10^3$km/s. The power spectrum of 2dFGRS seems in disfavor of models with low amplitude of density fluctuations.

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 study the X-ray emission of baryon fluid in the universe using the WIGEON cosmological hydrodynamic simulations. It has been revealed that cosmic baryon fluid in the nonlinear regime behaves like Burgers turbulence, i.e. the fluid field consists of shocks. Like turbulence in incompressible fluid, the Burgers turbulence plays an important role in converting the kinetic energy of the fluid to thermal energy and heats the gas. We show that the simulation sample of the $\Lambda$CDM model without adding extra heating sources can fit well the observed distributions of X-ray luminosity versus temperature ($L_{\rm x}$ vs. $T$) of galaxy groups and is also consistent with the distributions of X-ray luminosity versus velocity dispersion ($L_{\rm x}$ vs. $\sigma$). Because the baryonic gas is multiphase, the $L_{\rm x}-T$ and $L_{\rm x}-\sigma$ distributions are significantly scattered. If we describe the relationships by power laws $L_{\rm x}\propto T^{\alpha_{LT}}$ and $L_{\rm x}\propto \sigma^{\alpha_{LV}}$, we find $\alpha_{LT}>2.5$ and $\alpha_{LV}>2.1$. The X-ray background in the soft $0.5-2$ keV band emitted by the baryonic gas in the temperature range $10^5<T<10^7$ K has also been calculated. We show that of the total background, (1) no more than 2% comes from the region with temperature less than $10^{6.5}$ K, and (2) no more than 7% is from the region of dark matter with mass density $\rho_{\rm dm}<50 \bar{\rho}_{\rm dm}$. The region of $\rho_{\rm dm}>50\bar{\rho}_{\rm dm}$ is generally clustered and discretely distributed. Therefore, almost all of the soft X-ray background comes from clustered sources, and the contribution from truly diffuse gas is probably negligible. This point agrees with current X-ray observations.

We show that the algorithm based on the weighted essentially nonoscillatory (WENO) scheme with anti-diffusive flux corrections can be used as a solver of the radiative transfer equations. This algorithm is highly stable and robust for solving problems with both discontinuities and smooth solution structures. We test this code with the ionized sphere around point sources. It shows that the WENO scheme can reveal the discontinuity of the radiative or ionizing fronts as well as the evolution of photon frequency spectrum with high accuracy on coarse meshes and for a very wide parameter space. This method would be useful to study the details of the ionized patch given by individual source in the epoch of reionization. We demonstrate this method by calculating the evolution of the ionized sphere around point sources in physical and frequency spaces. It shows that the profile of the fraction of neutral hydrogen and the ionized radius are sensitively dependent on the intensity of the source.

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 study the SZ-effect-induced non-Gaussianity in the cosmic microwave background (CMB) fluctuation maps. If a CMB map is contaminated by the SZ effect of galaxies or galaxy clusters, the CMB maps should have similar non-Gaussian features as the galaxy and cluster fields. Using the WMAP data and 2MASS galaxy catalog we show that the non-Gaussianity of the 2MASS galaxies is imprinted on WMAP maps. The signature of non-Gaussianity can be seen with the 4^{th} order cross correlation between the wavelet variables of the WMAP maps and 2MASS clusters. The intensity of the 4^{th} order non-Gaussian features is found to be consistent with the contamination of the SZ effect of 2MASS galaxies. We also show that this non-Gaussianity can not be seen by the high order auto-correlation of the WMAP. This is because the SZ signals in the auto-correlations of the WMAP data generally is weaker than the WMAP-2MASS cross correlations by a factor f^2, which is the ratio between the powers of SZ effect map and the CMB fluctuations on the scale considered. Therefore, the ratio of high order auto-correlations of CMB maps to cross-correlations of the CMB maps and galaxy field would be effective to constrain the powers of SZ effect on various scales.

We investigate the statistical properties of cosmic baryon fluid in the nonlinear regime, which is crucial for understanding the large-scale structure formation of the universe. With the hydrodynamic simulation sample of the Universe in the cold dark matter model with a cosmological constant, we show that the intermittency of the velocity field of cosmic baryon fluid at redshift z=0 in the scale range from the Jeans length to about 16 Mpc/h can be extremely well described by She-Leveque's universal scaling formula. The baryon fluid also possesses the following features: (1) for volume weight statistics, the dissipative structures are dominated by sheets, and (2) the relation between the intensities of fluctuations is hierarchical. These results imply that the evolution of highly evolved cosmic baryon fluid is similar to a fully developed turbulence.

The dynamic evolution of the baryonic intergalactic medium (IGM) caused by the underlying dark matter gravity is governed by the Navier-Stokes equations in which many cooling and heating processes are involved. However, it has long been recognized that the growth mode dynamics of cosmic matter clustering can be sketched by a random force driven Burgers' equation if cooling and heating are ignored. Just how well the dynamics of the IGM can be described as a Burgers fluid has not been fully investigated probably because cooling and heating are essential for a detailed understanding of the IGM. Using IGM samples produced by a cosmological hydrodynamic simulation in which heating and cooling processes are properly accounted for, we show that the IGM velocity field in the nonlinear regime shows the features of a Burgers fluid, that is, when the Reynolds number is high, the velocity field consists of an ensemble of shocks. Consequently, (1) the IGM velocity $v$ is generally smaller than that of dark matter; (2) for the smoothed field, the IGM velocity shows tight correlation with dark matter given by $v \simeq s v_{dm}$, with $s<1$, such that the lower the redshift, the smaller $s$; (3) the velocity PDFs are asymmetric between acceleration and deceleration events; (4) the PDF of velocity difference $\Delta v=v(x+r)-v(x)$ satisfies the scaling relation for a Burgers fluid, i.e., $P(\Delta v)=(1 r^y)F(\Delta v/r^y)$. We find the scaling function and parameters for the IGM which are applicable to the entire scale range of the samples (0.26 - 8 h$^{-1}$ Mpc). These properties show that the similarity mapping between the IGM and dark matter is violated on scales much larger than the Jeans length of the IGM.

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.

We show that, in the halo model of large-scale structure formation, the difference between the Fourier and the DWT (discrete wavelet transform) power spectra provides a statistical measurement of the halos. This statistical quantity is free from parameters related to the shape of the mass profile and the identification scheme of halos. That is, the statistical measurement is invariant in the sense that models with reasonably defined and selected parameters of the halo models should yield the same difference of the Fourier and DWT spectra. This feature is useful to extract ensemble averaged properties of halos, which cannot be obtained with the identification of individual halo. To demonstrate this point, we show with WIGEON hydrodynamical simulation samples that the spectrum difference provides a quantitative measurement of the discrepancy of the distribution of baryonic gas from that of the underlying dark matter field within halos. We also show that the mass density profile of halos in physical space can be reconstructed with this statistical measurement. This profile essentially is the average over an ensemble of halos, including well virialized halos as well as halos with significant internal substructures. Moreover, this reconstruction is sensitive to the tail of the mass density profile. We showed that the profile with $1/r^3$ tail gives very different result from that of $1/r^2$. Other possible applications of this method are discussed as well.

(Abridged) The nonlinear evolution of a system consisting of baryons and dark matter is generally characterized by strong shocks and discontinuities. The baryons slow down significantly at postshock areas of gravitational strong shocks, which can occur in high overdense as well as low overdense regions. Consequently, the baryon fraction would be nonuniform on large scales. We studied these phenomena with simulation samples produced by the WENO hybrid cosmological hydrodynamic/N-body code. We find that the baryon fraction in high mass density regions is lower on average than the cosmic baryon fraction, and many baryons accumulate in the regions with moderate mass density to form a high baryon fraction phase (HBFP). In dense regions with rho>100, which are the possible hosts for galaxy clusters, the baryon fraction can be lower than the cosmic baryon fraction by about 10%--20% at z ~ 0. Our simulation samples show that about 3% of the cosmic baryon budget was hidden in the HBFP at redshift z=3, while this percentage increases to about 14% at the present day. The gas in the HBFP cannot be detected either by Ly-alpha forests of QSO absorption spectra or by soft X-ray background. That is, the HBFP would be missed in the baryon budget given by current observations.

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.

(Abridged)The temperature (T) and entropy (S) fields of baryonic gas/IGM are analyzed using simulation samples by a cosmological hydro/N-body code. We show that in the nonlinear regime the dynamical similarity between the IGM and dark matter will be broken by strong shocks. The heating and entropy production by the shocks breaks the IGM into multiple phases. The mean entropy, or the cosmological entropy floor, is found to be more than 100 h^(-1/3) keV cm^2 in all regions when z < 1. At redshift z ~ 2 - 3, high-entropy gas (S > 50 h^(-1/3) keV cm^2) mostly resides in areas on scales larger than 1 h^(-1) Mpc and with density $\rho_{dm}>10^2$. Therefore, gravitational shocks are an effective preheating mechanism of the IGM, and probably enough to provide the entropy excess of clusters and groups if the epoch of the gas falling in cluster cores is not earlier than z ~ 2 - 3. On the other hand, at redshifts z < 4, there is always a more than 90% volume of the low dark matter mass density regions filled by the IGM with T < 10^4.5 K. Therefore, the multiphased character and non-Gaussianity of the IGM field would explain the high-temperature and high-entropy gas observed in groups and clusters.

We investigate the relationship between the mass and velocity fields of the intergalactic medium (IGM) and dark matter. Although the evolution of the IGM is dynamically governed by the gravity of the underlying dark matter field, some statistical properties of the IGM inevitably decouple from those of the dark matter once the nonlinearity of the dynamical equations and the stochastic nature of the field is considered. With simulation samples produced by a hybrid cosmological hydrodynamic/N-body code, which is effective in capturing shocks and complicated structures with high precision, we find that the one-point distributions of the IGM field are systematically different from that of dark matter as follows: 1.) the one-point distribution of the IGM peculiar velocity field is exponential at least at redshifts less than 2, while the dark matter velocity field is close to a Gaussian field; 2.) although the one-point distributions of the IGM and dark matter are similar, the point-by-point correlation between the IGM and dark matter density fields significantly differs on all scales and redshifts analyzed; 3.) the one-point density distributions of the difference between IGM and dark matter fields are highly non-Gaussian and long tailed. These discrepancies violate the similarity between the IGM and dark matter and cannot be explained simply as Jeans smoothing of the IGM. However, these statistical discrepancies are consistent with the fluids described by stochastic-force driven nonlinear dynamics.

We analyze the statistical properties and dynamical implications of galaxy distributions in phase space for samples selected from the 2MASS Extended Source Catalog. The galaxy distribution is decomposed into modes $\delta({\bf k, x})$ which describe the number density perturbations of galaxies in phase space cell given by scale band $\bf k$ to ${\bf k}+\Delta {\bf k}$ and spatial range $\bf x$ to ${\bf x}+\Delta {\bf x}$. In the nonlinear regime, $\delta({\bf k, x})$ is highly non-Gaussian. We find, however, that the correlations between $\delta({\bf k, x})$ and $\delta({\bf k', x'})$ are always very weak if the spatial ranges (${\bf x}$, ${\bf x}+\Delta {\bf x}$) and (${\bf x'}$, ${\bf x'}+\Delta {\bf x'}$) don't overlap. This feature is due to the fact that the spatial locality of the initial perturbations is memorized during hierarchical clustering. The highly spatial locality of the 2MASS galaxy correlations is a strong evidence for the initial perturbations of the cosmic mass field being spatially localized, and therefore, consistent with a Gaussian initial perturbations on scales as small as about 0.1 h$^{-1}$ Mpc. Moreover, the 2MASS galaxy spatial locality indicates that the relationship between density perturbations of galaxies and the underlying dark matter should be localized in phase space. That is, for a structure consisting of perturbations on scales from $k$ to $ k+\Delta {k}$, the nonlocal range in the relation between galaxies and dark matter should {\it not} be larger than $|{\Delta {\bf x}}|=2\pi/|\Delta {\bf k}|$. The stochasticity and nonlocality of the bias relation between galaxies and dark matter fields should be no more than the allowed range given by the uncertainty relation $|{\Delta {\bf x}|| \Delta{\bf k}}|=2\pi$.

Using the O(N) model as an example, we investigate the self-interaction effects of inflaton on the dynamics of the primordial perturbations. When taking interactions into account, it is essential to employ a self-consistent off-equilibrium formalism to study the evolution of the inflationary background field and its fluctuations with the back-reaction effects. Within the Hartree factorization scheme, we show that the O(N) model has at least two observable remains left behind the off-equilibrium processes: the running spectral index of primordial density perturbations and the correlations between perturbation modes in phase space. We find that the running of the spectral index is fully determined by the rate of the energy transfer from the inflationary background field to its fluctuations via particle creation processes as well as the dynamics of the background field itself. Furthermore, the amplitude of the field fluctuations turns out to be scale-dependent due to the off-equilibrium evolution. As a consequence, the scale-dependence of fluctuations yields a correlation between the phase space modes of energy density perturbations, while the one-point function of the fluctuations in each Hartree mode is still Gaussian. More importantly, the mode-mode correlation of the primordial perturbations depends upon the dynamics of the self-interaction {\it as well as} the initial conditions of the inflation. Hence, we propose that the running spectral index and the correlation between phase-space modes would be two observable fossils to probe the epoch of inflation, even beyond.

We study the off-equilibrium effects of inflaton on the dynamics of primordial perturbations in the O(N) model. A self-consistent off-equilibrium formalism is employed to investigate the evolution of the inflationary background field and its fluctuations with the back-reaction effects. We find two observable remains left behind the off-equilibrium processes: the running spectral index of primordial density perturbations and the correlations between perturbation modes in phase space, which would serve as the imprints to probe the epoch of inflation, even beyond.

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 investigate the primordial scalar perturbations in the thermal dissipative inflation where the radiation component (thermal bath) persists and the density fluctuations are thermally originated. The perturbation generated in this model is hybrid, i.e. it consists of both adiabatic and isocurvature components. We calculate the fractional power ratio ($S$) and the correlation coefficient ($\cos\Delta$) between the adiabatic and the isocurvature perturbations at the commencing of the radiation regime. Since the adiabatic/isocurvature decomposition of hybrid perturbations generally is gauge-dependent at super-horizon scales when there is substantial energy exchange between the inflaton and the thermal bath, we carefully perform a proper decomposition of the perturbations. We find that the adiabatic and the isocurvature perturbations are correlated, even though the fluctuations of the radiation component is considered uncorrelated with that of the inflaton. We also show that both $S$ and $\cos \Delta$ depend mainly on the ratio between the dissipation coefficient $\Gamma$ and the Hubble parameter $H$ during inflation. The correlation is positive ($\cos\Delta > 0$) for strong dissipation cases where $\Gamma/H >0.2$, and is negative for weak dissipation instances where $\Gamma/H <0.2$. Moreover, $S$ and $\cos \Delta$ in this model are not independent of each other. The predicted relation between $S$ and $\cos\Delta$ is consistent with the WMAP observation. Other testable predictions are also discussed.

We show that the nonlinear evolution of the cosmic gravitational clustering is approximately spatial local in the $x$-$k$ (position-scale) phase space if the initial perturbations are Gaussian. That is, if viewing the mass field with modes in the phase space, the nonlinear evolution will cause strong coupling among modes with different scale $k$, but at the same spatial area $x$, while the modes at different area $x$ remain uncorrelated, or very weakly correlated. We first study the quasi-local clustering behavior with the halo model, and demonstrate that the quasi-local evolution in the phase space is essentially due to the self-similar and hierarchical features of the cosmic gravitational clustering. The scaling of mass density profile of halos insures that the coupling between $(x-k)$ modes at different physical positions is substantially suppressed. Using high resolution N-body simulation samples in the LCDM model, we justify the quasi-locality with the correlation function between the DWT (discrete wavelet transform) variables of the cosmic mass field. Although the mass field underwent a highly non-linear evolution, and the DWT variables display significantly non-Gaussian features, there are almost no correlations among the DWT variables at different spatial positions. Possible applications of the quasi-locality have been discussed.

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 case is made that the intracluster medium (ICM) in spiral-rich galaxy groups today probably has undergone much slower evolution than that in elliptical-rich groups and clusters. The environments of proto-clusters and proto-groups at z > 2 are likely similar to spiral-rich group environments at lower redshift. Therefore, like the ICM in spiral-rich groups today, the ICM in proto-groups and proto-clusters at z > 2 is predicted to be significantly multiphased. The QSO Lyman alpha forest in the vicinity of galaxies is an effective probe of the ICM at a wide range of redshift. Two recent observations of Lyman alpha absorption around galaxies by Adelberger et al. and by Pascarelle et al are reconciled, and it is shown that observations support the multiphase ICM scenario. Galaxy redshifts must be very accurate for such studies to succeed. This scenario can also explain the lower metallicity and lower hot gas fraction in groups.

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.

It has been recently found that the local fluctuations of the QSO's Ly$\alpha$ absorption spectrum transmitted flux show spiky structures. This implies that the mass fields of the intergalactic medium (IGM) is intermittent. This feature cannot be explained by the clustering evolution of cosmic mass field in the linear regimes and is also difficult to incorporate into the hierarchical clustering scenario. We calculate the structure functions and intermittent exponent of the IGM and HI for full hydrodynamical simulation samples. The result shows the intermittent features of the Ly$\alpha$ transmitted flux fluctuations as well as the mass field of the IGM. We find that within the error bars of current data, all the intermittent behavior of the simulation samples are consistent with the observation. This result is different from our earlier result (Pando et al 2002), which shows that the intermittent behavior of samples generated by pseudo-hydro simulation cannot be fitted with observed data. One difference between the pseudo-hydro and full hydro simulations is in treating the dynamical relation between the IGM (or HI) and dark matter fields. The former assumes that the IGM density distribution traces the underlying dark matter point-by-point on scales larger than the Jeans length in either the linear or nonlinear regimes. However, hydrodynamic studies have found that a statistical discrepancy between the IGM field and underlying dark matter in nonlinear regime is possible. We find that the point-by-point correlation between the IGM density perturbations and dark matter become weaker on comoving scales less than 2 h$^{-1}$ Mpc (in LCDM model), which is larger than the IGM Jeans length.

In this paper, we develop a method of performing the one-point statistics of a perturbed density field with a multiresolutional decomposition based on the discrete wavelet transform (DWT). We establish the algorithm of the one-point variable and its moments in considering the effects of Poisson sampling and selection function. We also establish the mapping between the DWT one-point statistics in redshift space and real space, i.e. the algorithm for recovering the DWT one-point statistics from the redshift distortion of bulk velocity, velocity dispersion, and selection function. Numerical tests on N-body simulation samples show that this algorithm works well on scales from a few hundreds to a few Mpc/h for four popular cold dark matter models. Taking the advantage that the DWT one-point variable is dependent on both the scale and the shape (configuration) of decomposition modes, one can design estimators of the redshift distortion parameter (beta) from combinations of DWT modes. When the non-linear redshift distortion is not negligible, the beta estimator from quadrupole-to-monopole ratio is a function of scale. This estimator would not work without adding information about the scale-dependence, such as the power-spectrum index or the real-space correlation function of the random field. The DWT beta estimators, however, do not need such extra information. Numerical tests show that the proposed DWT estimators are able to determine beta robustly with less than 15% uncertainty in the redshift range 0 < z < 3.

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)

We investigate the constraint on the abundance of primordial black holes (PBHs) and the spectral index $n$ of primeval density perturbations given by the ionizing photon background at the epoch of reionization. Within the standard inflationary cosmogony, we show that the spectral index $n$ of the power-law power spectrum of primeval density perturbations should be $n<$1.27. Since the universe is still optical thick at the reionization redshift $z\sim 6$ - 8, this constraint is independent of the unknown parameter of reheating temperature of the inflation. The ionizing photon background from the PBHs can be well approximated by a power law spectrum $J(\nu)\propto{\nu}^3$, which is greatly different from those given by models of massive stars and quasars.

We show that if the primordial classical perturbations were generated by the gravitational particle creation during inflation, and followed by an evolution of quantum-to-classical transition,the time dependent correlation of these perturbations is long-tailed with a correlation time larger than the Hubble-time. Consequently, the inflationary perturbations are locally scale-scale correlated. Hence, the interaction of the fields during inflation can be explored via the detection of the local scale-scale correlation of the CMB fluctuations.

In this paper, we develop a theory of redshift distortion of the galaxy power spectrum in the discrete wavelet transform (DWT) representation. Because the DWT power spectrum is dependent of both the scale and shape (configuration) of the decomposition modes, it is sensitive to distortion of shape of the field. On the other hand, the redshift distortion causes a shape distortion of distributions in real space with respect to redshift space. Therefore, the shape-dependent DWT power spectrum is useful to detect the effect of redshift distortion. We first established the mapping between the DWT power spectra in redshift and real space. The mapping depends on the redshift distortion effects of (1) bulk velocity, (2) selection function and (3) pairwise peculiar velocity. We then proposed $\beta$-estimators using the DWT off-diagonal power spectra. These $\beta$-estimators are model-free even when the non-linear redshift distortion effect is not negligible. Moreover, these estimators do not rely on the assumption of whether the pairwise velocity dispersion being scale-dependent. The tests with N-body simulation samples show that the proposed $\beta$-estimators can yield reliable measurements of $\beta$ with about 20% uncertainty for all popular dark matter models. We also develop an algorithm for reconstruction of the power spectrum in real space from the redshift distorted power spectrum. The numerical test also shows that the real power spectrum can be well recovered from the redshift distorted power spectrum.

We propose a method for estimating the large-scale rms bulk velocity of the cosmic mass field from the transmitted fluxes of Ly$\alpha$ forests. It is based on two linear relationships on large scales: 1) the relation between the fluctuations of the transmission and the underlying density field, and 2) the relation between the density fluctuations and the peculiar velocity field. We show that, with a multiscale decomposition, the two relations can be effectively employed for predicting the rms bulk velocity. Since QSO's Ly$\alpha$ forest is due to the absorptions of diffusely distributed and photoionized IGM, this method provides an independent estimate of the rms bulk velocity at high redshifts, on large scales, and free from the bias of galaxies. Using the transmitted flux of 60 moderate-resolution QSO spectra, the rms bulk velocity is found to be 230$\pm$50 km s$^{-1}$ around redshift $z=2.25$ on scale 23 h$^{-1}$Mpc, and down to 110$\pm$45 km s$^{-1}$ around $z=3.25$ on scale 92 h$^{-1}$Mpc for an LCDM universe ($\Omega=0.3$ and $\Lambda=0.7$). The results are basically consistent with the linear evolution theory.

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.

In this paper, we develop the method of analyzing the velocity field of cosmic matter with a multiresolution decomposition. This is necessary in calculating the redshift distortion of power spectrum in the discrete wavelet transform (DWT) representation. We show that, in the DWT analysis, the velocity field can be described by discrete variables, which are given by assignment of the number density and velocity into the DWT modes. These DWT variables are complete and not redundant. In this scheme, the peculiar velocity and pairwise velocity of galaxies or particles are given by field variables. As a consequence, the velocity dispersion (VD) and pairwise velocity dispersion (PVD) are no longer measured by number-counting or pair-counting statistic, but with the ensemble of the field variables, and therefore, they are free from the bias due to the number-counting and pair-counting. We analyzed the VD and PVD of the velocity fields given by the N-body simulation for models of the SCDM, $\tau$CDM and $\Lambda$CDM. The spectrum (scale-dependence) of the VD and PVD show that the length scale of the two-point correlation of the velocity field is as large as few tens h$^{-1}$ Mpc. Although the VD and PVD show similar behavior in some aspects, they are substantially different from each other. The VD-to-PVD ratio shows the difference between the scale-dependencies of the VD and PVD. More prominent difference between the VD and PVD is shown by probability distribution function. The one-point distribution of peculiar velocity is approximately exponential, while the pairwise velocity's is lognormal, i.e. of long tail. This difference indicates that the cosmic velocity field is typically intermittent.

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 hierarchical clustering model of the cosmic mass field is examined in the context of intermittency. We show that the mass field satisfying the correlation hierarchy $\xi_n\simeq Q_n(\xi_2)^{n-1}$ is intermittent if $\kappa < d$, where $d$ is the dimension of the field, and $\kappa$ is the power-law index of the non-linear power spectrum in the discrete wavelet transform (DWT) representation. We also find that a field with singular clustering can be described by hierarchical clustering models with scale-dependent coefficients $Q_n$ and that this scale-dependence is completely determined by the intermittent exponent and $\kappa$. Moreover, the singular exponents of a field can be calculated by the asymptotic behavior of $Q_n$ when $n$ is large. Applying this result to the transmitted flux of HS1700 Ly$\alpha$ forests, we find that the underlying mass field of the Ly$\alpha$ forests is significantly intermittent. On physical scales less than about 2.0 h$^{-1}$ Mpc, the observed intermittent behavior is qualitatively different from the prediction of the hierarchical clustering with constant $Q_n$. The observations, however, do show the existence of an asymptotic value for the singular exponents. Therefore, the mass field can be described by the hierarchical clustering model with scale-dependent $Q_n$. The singular exponent indicates that the cosmic mass field at redshift $\sim 2$ is weakly singular at least on physical scales as small as 10 h$^{-1}$ kpc.

We investigate the weakly non-linear evolution of cosmic gravitational clustering in phase space by looking at the Zel'dovich solution in the discrete wavelet transform (DWT) representation. We show that if the initial perturbations are Gaussian, the relation between the evolved DWT mode and the initial perturbations in the weakly non-linear regime is quasi-local. That is, the evolved density perturbations are mainly determined by the initial perturbations localized in the same spatial range. Furthermore, we show that the evolved mode is monotonically related to the initial perturbed mode. Thus large (small) perturbed modes statistically correspond to the large (small) initial perturbed modes. We test this prediction by using QSO Ly$\alpha$ absorption samples. The results show that the weakly non-linear features for both the transmitted flux and identified forest lines are quasi-localized. The locality and monotonic properties provide a solid basis for a DWT scale-by-scale Gaussianization reconstruction algorithm proposed by Feng & Fang (Feng & Fang, 2000) for data in the weakly non-linear regime. With the Zel'dovich solution, we find also that the major non-Gaussianity caused by the weakly non-linear evolution is local scale-scale correlations. Therefore, to have a precise recovery of the initial Gaussian mass field, it is essential to remove the scale-scale correlations.

The power spectrum estimator based on the discrete wavelet transform (DWT) for 3-dimensional samples has been studied. The DWT estimator for multi-dimensional samples provides two types of spectra with respect to diagonal and off-diagonal modes, which are very flexible to deal with configuration-related problems in the power spectrum detection. With simulation samples and mock catalogues of the Las Campanas redshift survey (LCRS), we show (1) the slice-like geometry of the LCRS doesn't affect the off-diagonal power spectrum with ``slice-like'' mode; (2) the Poisson sampling with the LCRS selection function doesn't cause more than 1-$\sigma$ error in the DWT power spectrum; and (3) the powers of peculiar velocity fluctuations, which cause the redshift distortion, are approximately scale-independent. These results insure that the uncertainties of the power spectrum measurement are under control. The scatter of the DWT power spectra of the six strips of the LCRS survey is found to be rather small. It is less than 1-$\sigma$ of the cosmic variance of mock samples in the wavenumber range $0.1 < k < 2$ h Mpc$^{-1}$. To fit the detected LCRS diagonal DWT power spectrum with CDM models, we find that the best-fitting redshift distortion parameter $\beta$ is about the same as that obtained from the Fourier power spectrum. The velocity dispersions $\sigma_v$ for SCDM and $\Lambda$CDM models are also consistent with other $\sigma_v$ detections with the LCRS. A systematic difference between the best-fitting parameters of diagonal and off-diagonal power spectra has been significantly measured. This indicates that the off-diagonal power spectra are capable of providing information about the power spectrum of galaxy velocity field.

We investigated the evolution of the primordial density perturbations produced by inflation with thermal dissipation. A full relativistic analysis on the evolution of initial perturbations from the warm inflation era to a radiation-dominated universe has been developed. The emphasis is on tracking the ratio between the adiabatic and the isocurvature mode of the initial perturbations. This result is employed to calculate a testable factor: the super-Hubble suppression of the power spectrum of the primordial perturbations. We show that based on the warm inflation scenario, the super-Hubble suppression factor, $s$, for an inflation with thermal dissipation is at least 0.5. This prediction does not depend on the details of the model parameters. If $s$ is larger than 0.5, it implies that the friction parameter $\Gamma$ is larger than the Hubble expansion parameter $H$ during the inflation era.

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.

We present a method of measuring galaxy power spectrum based on the multiresolution analysis of the discrete wavelet transformation (DWT). Since the DWT representation has strong capability of suppressing the off-diagonal components of the covariance for selfsimilar clustering, the DWT covariance for popular models of the cold dark matter cosmogony generally is diagonal, or $j$(scale)-diagonal in the scale range, in which the second scale-scale correlations are weak. In this range, the DWT covariance gives a lossless estimation of the power spectrum, which is equal to the corresponding Fourier power spectrum banded with a logarithmical scaling. In the scale range, in which the scale-scale correlation is significant, the accuracy of a power spectrum detection depends on the scale-scale or band-band correlations. This is, for a precision measurements of the power spectrum, a measurement of the scale-scale or band-band correlations is needed. We show that the DWT covariance can be employed to measuring both the band-power spectrum and second order scale-scale correlation. We also present the DWT algorithm of the binning and Poisson sampling with real observational data. We show that the alias effect appeared in usual binning schemes can exactly be eliminated by the DWT binning. Since Poisson process possesses diagonal covariance in the DWT representation, the Poisson sampling and selection effects on the power spectrum and second order scale-scale correlation detection are suppressed into minimum. Moreover, the effect of the non-Gaussian features of the Poisson sampling can be calculated in this frame.

We investigate the effect of non-Gaussianity on the reconstruction of the initial mass field from the Ly$\alpha$ forest. We show that the transmitted flux of QSO absorption spectra are highly non-Gaussian in terms of the statistics, the kurtosis spectrum and scale-scale correlation. These non-Gaussianities can not be completely removed by the conventional algorithm of Gaussianization, and the scale-scale correlations are largely retained in the mass field recovered by the Gaussian mapping. Therefore, the mass power spectrum recovered by the conventional algorithm is systematically lower than the initial mass spectrum on scales at which the local scale-scale correlation is substantial. To reduce the non-Gaussian contamination, we present two methods. The first is to perform the Gaussianization scale-by-scale using the discrete wavelet transform (DWT) decomposition. We show that the non-Gaussian features of the Ly$\alpha$ forest basically will no longer exist in the scale-by-scale Gaussianized mass field. The second method is to choose a proper orthonormal basis (representation) to suppress the effect of the non-Gaussian correlations. In the quasilinear regime of cosmic structure formation, the DWT power spectrum is efficient for suppressing the non-Gaussian contamination. These two methods significantly improve the recovery of the mass power spectrum from the Ly$\alpha$ forest.

We show that scale-scale correlations are a generic feature of slow-roll inflation theories. These correlations result from the long-time tails characteristic of the time dependent correlations because the long wavelength density perturbation modes are diffusion-like. A relationship between the scale-scale correlations and time-correlations is established providing a way to reveal the time correlations of the perturbations during inflation. This mechanism provides for a testable prediction that the scale-scale correlations at two different spatial points will vanish.

Using samples of structures identified by a multi-scale decomposition from numerical simulation, we analyze the scale-dependence of the virialization of clusters. We find that beyond the scale of full virialization there exists a radius range over which clusters are quasi-virialized, i.e. while the internal structure of an {\it individual} cluster is at substantial departure from dynamical relaxation, some {\it statistical} properties of the multi-scale identified clusters are approximately the same as those for the virialized systems. The dynamical reason of the existence of quasi-virialization is that some of the scaling properties of dynamically relaxed systems of cosmic gravitational clustering approximately hold beyond the full virialization regime. The "individual-statistical" duality of the quasi-virialization provides an explanation of the observed puzzle that the total masses of clusters derived from virial theorem are statistically the same as the gravitational lensing determined masses, in spite of the presence of irregular configuration and substructures in individual clusters. It also explains the tight correlation between the velocity dispersion of optical galaxies and the temperature of X-ray emitting gas. Consequently, the virial mass estimators based on the assumptions of isothermal and hydrostatic model are statistically applicable to scales on which the clusters are quasi-virialized. In the quasi-virialization regime, the temperature functions of clusters also show scaling. This feature is a useful discriminator among cosmological models.

Recently, scale-scale correlations have been detected in the distributions of quasar's Ly$\alpha$ absorption lines and the maps of cosmic temperature fluctuations. In this paper, we investigate the scale-scale correlations in galaxy distributions. Using samples of mass field given by N-body simulation, we first show that the scale-scale correlation of 2-D and 3-D mass distributions at present day is capable of breaking the degeneracy between the SCDM (standard cold dark matter model) and OCDM (open CDM model) or LCDM (flat CDM model), and even show the difference between the OCDM and LCDM. Using biased galaxy samples produced in an appropriate bias model, we show that the scale-scale correlation of galaxy distribution at zero redshift is still a powerful tool to break the degeneracy in the parameter space, including both cosmological and biasing parameters. We analyze the scale-scale correlations of the APM bright galaxy catalog, and compare it with the mock catalog in the SCDM, OCDM and LCDM models. We find that all the spectra of local and non-local scale-scale correlations predicted by the OCDM model are in excellent agreement with those of the APM-BGC sample, while the SCDM and LCDM mock catalog appear to have somewhat weaker scale-scale correlations than the observation.

It has been suggested that the mass-to-light ($M/L$) ratio of gravitationally clustering objects is scale-independent on scales beyond galaxy clusters, and may also be independent of the mass of the objects. In this paper, we show that the scale behavior of $M/L$ ratio is closely related to the scaling of cosmic structures larger than clusters. The scale dependence of the $M/L$ ratio can be determined by comparing the observed scaling of richness function (RF) of multi-scale identified objects with the model-predicted scaling of mass function (MF) of large scale structures. Using the multi-scale identified clusters from IRAS 1.2 Jy galaxy survey, we have made comparisons of the observed RF scaling of IRAS $r_{cl}$-clusters with the MF scalings given by simulations of three popular models SCDM, LCDM and OCDM. We find that, the M/L ratio basically is scale-independent from the Abell radius up to about 24 $h^{-1}$Mpc, while it seems to show a slight, but systematical, increase over this scale range. This result is weakly dependent on the cosmological parameters.

The recently observed bright optical transients(OT) of high redshift GRBs indicate that they are in a violent dynamical state. We think it is reasonable to assume that the GRBs form in the environment of gravitationally collapsed halos of the cosmic matter field, and we investigate the basic parameters of the halos which are favored to host GRBs. If the harboring coefficient $f$ of GRBs per halo is weakly dependent on the mass of the halo, the redshift data of GRB OTs can yield significant constraints on the massive halos hosting GRBs. We show that, in the framework of popular cold dark matter (CDM) models, the GRB-favored environments are newly collapsed halos (i.e. their ages less than about $2 \times 10^9$ yr) with masses around $10^9$ $h^{-1}$ M$_{\odot}$. In this scenario, low redshift GRBs, if they exist, could not have the same cosmic origin as the high redshift ones. To fit with the observed rate of GRBs, we conclude that each GRB halo can host probably no more than one GRB event on average. This result implies that GRBs may be related to the merging of the halos.

We investigated the structures on scales beyond the typical clusters of galaxies. These structures are crucial to understand the cosmic gravitational clustering in the pre-virialized stage, or quasilinear r\`egime. Based on the multi-resolution analysis of the discrete wavelet transformation, we got statistical available ensembles of $r_{cl}$-clusters, i.e. the structures on scale $r_{cl}$, in the range $1 \leq r_{cl} \leq 24 h^{-1}$ Mpc for both N-body simulation and the IRAS 1.2 Jy galaxy survey samples. If the mass-to-light ratio on scales larger than clusters asymptotically reaches a constant, we found that the abundances and correlations of these IRAS $r_{cl}$-clusters to be basically consistent with the predictions of the flat low-density CDM model (LCDM) and the open CDM model (OCDM), except the model-predicted abundance of $r_{cl} = 24 h^{-1}$ Mpc clusters seems to be higher than IRAS data. The standard CDM (SCDM) gives too much power, and too weak correlations on all the scales. For a given $r_{cl}$, the amplitude of two-point correlation function of $r_{cl}$-clusters is increasing with their richness. However, for a given richness (defined by the mean separation of neighbor objects), the clustering strengths of both simulation and observation sample are found to be declining with $r_{cl}$ when $r_{cl}$ is larger than 3 - 4 $h^{-1}$ Mpc. Therefore, the ``universal'' increase of the correlation amplitude with the scale of objects from galaxies, groups, to poor, rich clusters is broken down for structures of $r_{cl} > 3 - 4 h^{-1}$ Mpc. Supercluster should not be a member of the ``universal'' increase family.

The relationships between the X-ray determined bolometric luminosity Lx, the temperature T of the intracluster gas, and the optical measured velocity dispersion sigma of the cluster galaxies are updated for galaxy clusters using the largest sample of 256 clusters drawn from literature. The newly established relationships, based on the doubly weighted orthogonal distance regression (ODR) method, are justified by both their self-consistency and co-consistency, which can then be used to test the theoretical models of cluster formation and evolution. The observationally determined Lx-T and Lx-sigma relationships, $L_x\propto T^{2.72\pm0.05}\propto \sigma^{5.24\pm0.29}$, are marginally consistent with those predicted in the scenario that both intracluster gas and galaxies are in isothermal and hydrostatic equilibrium with the underlying gravitational potential of clusters. A comparison between these observed and predicted Lx-T relationships also suggests that the mean cluster baryon fraction fb remains approximately constant among different clusters, $f_b\approx0.17$, which gives rise to a low-mass density universe of $\Omega_m\approx0.3$.

We study the power spectrum of the mass density perturbations in an inflation scenario that includes thermal dissipation. We show that the condition on which the thermal fluctuations dominate the primordial density perturbations can easily be realized even for weak dissipation, i.e., the rate of dissipation is less than the Hubble expansion. We find that our spectrum of primordial density perturbations follows a power law behavior, and exhibits a ``thermodynamical'' feature -- the amplitude and power index of the spectrum depend mainly on the thermodynamical variable $M$, the inflation energy scale. Comparing this result with the observed temperature fluctuations of the cosmic microwave background, we find that both amplitude and index of the power spectrum can be fairly well fitted if $M \sim 10^{15}-10^{16}$ GeV.

We perform a discrete wavelet analysis of the COBE-DMR 4yr sky maps and find a significant scale-scale correlation on angular scales from about 11 to 22 degrees, only in the DMR face centered on the North Galactic Pole. This non-Gaussian signature does not arise either from the known foregrounds or the correlated noise maps, nor is it consistent with upper limits on the residual systematic errors in the DMR maps. Either the scale-scale correlations are caused by an unknown foreground contaminate or systematic errors on angular scales as large as 22 degrees, or the standard inflation plus cold dark matter paradigm is ruled out at the $> 99%$ confidence level.

We present an analysis of the non-Gaussianity in the distribution of Ly$\alpha$ forest lines in the QSO absorption spectra. Statistical tests performed on this data indicate that there may be large scale structure even though the power spectrum of the Ly$\alpha$ line distribution on large scales is found to be flat. It is apparent that higher (than two) order statistics are crucial in quantifying the clustering behavior of Ly$\alpha$ clouds. Using the discrete wavelet on three independent data sets of Ly$\alpha$ forests, we find that the distribution of Ly$\alpha$ forests does show non-Gaussian behavior on scales from 5 to 10 h$^{-1}$ Mpc with confidence level larger than 95%. Two data sets available on large scales are found to be non-Gaussian on even larger scales. These techniques are effective in discriminating among models of the Ly$\alpha$ forest formation, which are degenerate at second and lower order statistics (abridged).

The relationship between the X-ray determined temperature $T$ of the intracluster gas and the optical measured velocity dispersion $\sigma$ of the cluster galaxies is often believed to be not only a straightforward but also robust test for the dynamical properties of galaxy clusters. Here, we present the $\sigma$-$T$ relationship using the 94 clusters drawn from the largest sample of 149 clusters in literature, for which both $\sigma$ and $T$ are observationally determined. Employment of the doubly weighted orthogonal distance regression to our sample yields $\sigma=10^{2.47\pm0.06}T^{0.67\pm0.09}$, indicating an apparent deviation of dynamical state from that predicted by the isothermal and hydrostatic equilibrium model for galaxy clusters, though the average ratio $\beta_{spec}$ of specific energy in galaxies to that in gas is found to be in excellent agreement with unity. It shows that a nonisothermal gas distribution with a mean polytropic index of $\gamma=1.3$ can account for the reported $\sigma$-$T$ relationship, while overall clusters can still be regarded as dynamically-relaxed systems.

Rich and massive clusters of galaxies at intermediate redshift are capable of magnifying and distorting the images of background galaxies. A comparison of different mass estimators among these clusters can provide useful information about the distribution and composition of cluster matter and their dynamical evolution. Using a hitherto largest sample of lensing clusters drawn from literature, we compare the gravitating masses of clusters derived from the strong/weak gravitational lensing phenomena, from the X-ray measurements based on the assumption of hydrostatic equilibrium, and from the conventional isothermal sphere model for the dark matter profile characterized by the velocity dispersion and core radius of galaxy distributions in clusters. While there is an excellent agreement between the weak lensing, X-ray and isothermal sphere model determined cluster masses, these methods are likely to underestimate the gravitating masses enclosed within the central cores of clusters by a factor of 2--4 as compared with the strong lensing results. Such a mass discrepancy has probably arisen from the inappropriate applications of the weak lensing technique and the hydrostatic equilibrium hypothesis to the central regions of clusters as well as an unreasonably large core radius for both luminous and dark matter profiles. Nevertheless, it is pointed out that these cluster mass estimators may be safely applied on scales greater than the core sizes. Namely, the overall clusters of galaxies at intermediate redshift can still be regarded as the dynamically relaxed systems, in which the velocity dispersion of galaxies and the temperature of X-ray emitting gas are good indicators of the underlying gravitational potentials of clusters.

We investigate the dependence of cluster abundance $n(>M,r_{cl})$, i.e., the number density of clusters with mass larger than $M$ within radius $r_{cl}$, on scale parameter $r_{cl}$. Using numerical simulations of clusters in the CDM cosmogonic theories, we notice that the abundance of rich clusters shows a simple scale invariance such that $n[>(r_{cl}/r_0)^{\alpha}M, r_{cl}]= n(>M,r_0)$, in which the scaling index $\alpha$ remains constant in a scale range where halo clustering is fully developed. The abundances of scale $r_{cl}$ clusters identified from IRAS are found basically to follow this scaling, and yield $\alpha \sim 0.5$ in the range $1.5 < r_{cl} < 4 h^{-1}$Mpc. The scaling gains further supports from independent measurements of the index $\alpha$ using samples of X-ray and gravitational lensing mass estimates. We find that all the results agree within error limit as: $\alpha \sim 0.5 - 0.7$ in the range of $1.5 < r_{cl} < 4 h^{-1}$Mpc. These numbers are in good consistency with the predictions of OCDM ($\Omega_M=0.3$) and LCDM ($\Omega_M+\Omega_{\Lambda} =1$), while the standard CDM model has different behavior. The current result seems to favor models with a low mass density.

We present an investigation of the scale-dependence of bias described by the linear model: $(\delta \rho({\bf x})/\bar{\rho})_{g} = b (\delta \rho(x)/\bar{\rho})_{m}$, $b$ being the bias parameter, and $\rho({\bf x})_{g}$ and $\rho({\bf x})_{m}$ are the galaxy number density and mass density, respectively. Using a discrete wavelet decomposition, we show that the behavior of bias scale-dependence cannot be described by one parameter $b$. In the linear bias model the scale-dependence should be measured by the $j$-spectra of wavelet-coefficient-represented bias parameters $\tilde{b}^{(n)}_j$ and $b_j^{(n)}$, $n$ being positive integers. Because $\tilde{b}^{(n)}_j$ with different $n$ are independent from each other, a systematic analysis of the $j$-spectra of $\tilde{b}^{(n)}_j$ and $b_j^{(n)}$ is necessary. We performed a $j$-spectrum analysis for samples of elliptical and lenticular (EL), and spiral (SP) galaxies listed in the APM bright galaxy catalog. We found that, for statistics of two-point correlation functions or DWT power spectrum, the scale-independence holds within 1 $\sigma$. However, the bias scale-dependence becomes substantial when phase-sensitive statistics (e.g. $\tilde{b}^{(n)}_j$ with $n>2$ or $b_j^{(n)}$) are applied. These results indicate that the bias scale-dependence has the same origin as the non-Gaussianity of galaxy distributions. This is generally consistent with the explanation that the bias scale-dependence originated from non-linear and non-local relationship between galaxy formation and their environment.

He II Ly-alpha 304/H-alpha 1640 emission lines are mainly produced by recombination, and their canonical ratio of about 10 may be a sensitive reddening indicator. We obtain the high S/N optical spectra of two quasars and combine them with the far-UV spectra that show the He II 304 emission. For HS 1700+64, the He II 1640 emission is not detected, and an upper limit to it sets the ratio greater than 20. For Q0302-003, the ratio is very low, on the order of unity. The most plausible cause for such a low ratio is extinction in the EUV band by very fine grains of dust. Q0302-003 has a prominent narrow component of FWHM ~ 2000 km/s in its major emission lines, and it appears that reddening is associated only with the line-emitting region. We suggest that the geometry of the line-emitting region in high-z quasars resembles that in the low-luminosity active galaxies, with the presence of dust mostly in the outer part.

We propose that the majority of quasars at redshift $z\sim 1 - 5$ formed in the environment of new born collapsed halos with 1-D velocity dispersion $\sigma_v^{1d} \sim 400 \kms$. The harboring coefficient $f$ of quasars per halo and the lifetime of quasars depend only on local process, not modulated by the density inhomogeneities on scales larger than the size of the halos. Thus, the bias of quasars on scale larger than the size of these halos is mainly determined by the parameter $\sigma_v$ used for quasar environment identification. With this model, the popular structure formation models, like SCDM and LCDM, can be fairly well reconciled with the data of quasars, including a. observed feature of the environment for quasars; b. redshift evolution of quasar abundance; c. the two-point correlation functions of quasars. This bias model predicts that the correlation function of quasars doesn't significantly evolve, or only slightly increases with redshift.

The energy source has remained to be the great mystery in understanding of the gamma-ray bursts (GRBs) if the events are placed at cosmological distances as indicated by a number of recent observations. The currently popular models include (1)the merger of two neutron stars or a neutron star and a black hole binary and (2)the hypernova scenario of the collapse of a massive member in a close binary. Since a neutron star will inevitably collapse into a black hole if its mass exceeds the limit $M_{max}\approx3M_{\odot}$, releasing a total binding gravitational energy of $\sim10^{54}$ erg, we explore semi-empirically the possibility of attributing the energy source of GRB to the accretion- induced collapse of a neutron star (AICNS) in a massive X-ray binary system consisting of a neutron star and a type O/B companion. This happens because a significant mass flow of $\sim10^{-3}$--$10^{-4}M_{\odot}$ yr$^{-1}$ may be transferred onto the neutron star through the Roche-lobe overflow and primarily during the spiral-in phase when it plunges into the envelope of the companion, which may eventually lead to the AICNS before the neutron star merges with the core of the companion. In this scenario, a ``dirty'' fireball with a moderate amount of beaming is naturally expected because of the nonuniformity of the stellar matter surrounding the explosion inside the companion, and a small fraction ($\sim0.1%$) of the energy is sufficient to create the observed GRBs. In addition, the bulk of the ejecting matter of the companion star with a relatively slow expansion rate may act as the afterglow. Assuming a non-evolutionary model for galaxies, we estimate that the birthrate of the AICNS events is about 2 per day within a volume to redshift $z=1$ for an $\Omega_0=1$ universe, consistent with the reported GRB rate.

It is of fundamental importance to determine if and how hierarchical clustering is involved in large-scale structure formation of the universe. Hierarchical evolution is characterized by rules which specify how dark matter halos are formed by the merging of halos at smaller scales. We show that scale-scale correlations of the matter density field are direct and sensitive measures to quantify this merging tree. Such correlations are most conveniently determined from discrete wavelet transforms. Analyzing two samples of Ly-alpha forests of QSO's absorption spectra, we find significant scale-scale correlations whose dependence is typical for a branching process. Therefore, models which predict a "history" independent evolution are ruled out and the halos hosting the Ly-alpha clouds must have gone through a "history" dependent merging process during their formation.

The theoretically expected amplitude of the associations of background quasars with foreground galaxies as a result of gravitational lensing has been updated in this paper. Since the galactic matter alone yields an amplitude of quasar overdensity smaller than that observed, a special attention has been paid to the examination or re-examination of the uncertainties in the estimate of the quasar enhancement factor arising from the cosmic evolution of galaxies, the core radius and velocity bias of galactic matter distributions, the clusters of galaxies, the obstruction effect by galactic disks, the non-zero cosmological constant, etc. Unfortunately, none of these factors has been shown to be able to significantly improve the situation, although a combination of some effects may provide a result that marginally agrees with observations. It is concluded that the quasar-galaxy association still remains to be an unsolved puzzle in today's astronomy, if the reported quasar-galaxy associations are not due to the statistical variations and/or the observed quasar number counts as a whole have not been seriously contaminated by gravitational lensing.

Inflationary cosmology predicts that the particle horizon should be generically much bigger than the present-day Hubble radius, $1/H_0$. This implies a special regime of super-Hubble scale energy density fluctuations imprinted on the cosmic microwave background radiation (CMBR), which from present theory could only be explained by inflation Causality constraints are used to determine models for the power spectrum that accommodate a suppression scale. A three parameter likelihood analysis is performed of the COBE-DMR 4-year data with respect to the amplitude, spectral index, and suppression scale. It is found that all suppression length scales larger than $1/H_0$ are consistent with the data, but that scales of order $1/H_0$ are slightly preferred, at roughly the one-sigma level. Many non-inflation models would be consistent with a small suppression length scale, whereas for standard inflation models, the duration of the inflation epoch would have to be bounded by a fairly small upper limit. Suppression scales smaller than $1/H_0$ are strongly excluded by the anisotrophy data.

We study thermally induced density perturbations during inflation. This scenario is characterized by two thermodynamical conditions: (1) The primordial perturbations originate in the epoch when the inflationary universe contains a thermalized heat bath. (2) The perturbations of the inflationary scalar field are given by the fluctuation-dissipation relation. We show that the spectrum of the primordial perturbations is of power law, but tilted, and there is a relation between the amplitude and the index of the power spectrum. Aside from the mass scale of the inflation, the amplitude-index relation does not depend on other parameters like $g$-factor. These results are found to be well consistent with observations of the temperature fluctuations of cosmic microwave background if the mass scale of the inflation is about $10^{15}$ GeV. Instead of the purely adiabatic case, the consequent density perturbation is an admixture of adiabatic and isocurvature one. Therefore, the detection of super-Hubble suppression of the spectrum would be effective for further discrimination between the thermally originated models and others.

Using numerical simulations of cluster formation in the standard CDM model (SCDM) and in a low-density, flat CDM model with a cosmological constant (LCDM), we investigate the gravitational lensing explanation for the reported associations between background quasars and foreground clusters. Under the thin-lens approximation and the unaffected background hypothesis , we show that the recently detected quasar overdensity around clusters of galaxies on scales of $\sim10$ arcminutes cannot be interpreted as a result of the gravitational lensing by cluster matter and/or by their environmental and projected matter along the line of sight, which is consistent with the analytical result based on the observed cluster and galaxy correlations (Wu, et al. 1996). It appears very unlikely that uncertainties in the modeling of the gravitational lensing can account for the disagreement between the theoretical predictions and the observations. We conclude that either the detected signal of the quasar-cluster associations is a statistical fluke or the associations are are generated by mechanisms other than the magnification bias.

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 present a detailed review of large-scale structure (LSS) study using the discrete wavelet transform (DWT). After describing how one constructs a wavelet decomposition we show how this bases can be used as a complete statistical discription of LSS. Among the topics studied are the the DWT estimation of the probability distribution function; the reconstruction of the power spectrum; the regularization of complex geometry in observational samples; cluster identification; extraction and identification of coherent structures; scale-decomposition of non-Gaussianity, such as spectra of skewnes and kurtosis and scale-scale correlations. These methods are applied to both observational and simulated samples of the QSO Lyman-alpha forests. It is clearly demonstrated that the statistical measures developed using the DWT are needed to distinguish between competing models of structure formation. The DWT also reveals physical features in these distributions not detected before. We conclude with a look towards the future of the use of the DWT in LSS.

We present a statistical comparison of three different estimates of cluster mass, namely, the dynamical masses obtained from the velocity dispersion of optical galaxies, the X-ray masses measured from the temperature of X-ray emitting gas under the assumption of isothermal hydrostatic equilibrium, and the gravitational lensing masses derived from the strong/weak distortions of background galaxy images. Using a sample of 29 lensing clusters available in literature, we have shown that the dynamical masses are in agreement with the gravitational lensing masses, while the X-ray method has systematically underestimated cluster masses by a factor 2-3 as compared with the others. These results imply that galaxies indeed trace the gravitational potential of their clusters, and there is no bias between the velocities of the dark matter particles and the galaxies in clusters. The X-ray cluster mass discrepancy is probably from the simplification in the models for the X-ray gas distribution and dynamical evolution.

Gravitational lensing by clusters of galaxies has been detected on scales ranging from $\sim10^{-1}$ Mpc to $\sim10$ Mpc, namely, arcs/arclets, weak lensing and quasar-cluster associations. This allows us to derive an overall radius matter distribution of clusters of galaxies. While the dynamical analysis of the X-ray observations has yielded a great number of data for the virial cluster masses, it becomes possible to statistically compare the cluster mass determinations by these two independent methods. In this letter we show that as compared with gravitational lensing, the dynamical analysis under the assumption of isothermal and hydrostatic equilibrium has systematically underestimated the cluster masses inside the Abell radius by a factor of $\sim2$ with scatter between $0.7$ and $5$. Because the same correction factor should be applicable to the gas baryon fraction of clusters of galaxies obtained from the X-ray data, it is probably too premature to claim a baryon crisis in today's cosmology.

Using a discrete wavelet based space-scale decomposition (SSD), the spectrum of the skewness and kurtosis is developed to describe the non-Gaussian signatures in cosmologically interesting samples. Because the basis of the discrete wavelet is compactly supported, the one-point distribution of the father function coefficients (FFCs) taken from one realization is a good estimate of the probability distribution function of the density if the ``fair sample hypothesis" holds. These FFC one-point distributions can also avoid the constraints of the central limit theorem on the detection of non-Gaussianity. Thus the FFC one-point distributions are effective in detecting non-Gaussian behavior in samples such as non-Gaussian clumps embedded in a Gaussain background, regardless of the number or density of the clumps. We demonstrate that the non-Gaussianity can reveal not only the magnitudes but also the scales of non-Gaussianity. Also calculated are the FFC one-point distributions, skewness and kurtosis spectra for real data and linearly simulated samples of QSO Ly$\alpha$ forests. When considering only second and lower order of statistics, such as the number density and two-point correlation functions, the simulated data show the same features as the real data. However, the the kurtosis spectra of samples given by different models are found to be different. On the other hand, the spectra of skewness and kurtosis for independent observational data sets are found to be the same. Moreover, the real data are significantly different from the non-Gaussianity spectrum of various posssible random samples. Therefore the non-Gaussain spectrum is necessary and valuable for model discrimination.

The associations of the angular positions of background quasars with foreground galaxies, clusters of galaxies and quasars are often attributed to the statistical lensing by gravitational potentials of the matter along the lines of sight, although it has been known that none of the individual objects (galaxies, clusters or quasars) are able to fully explain the reported amplitudes of the quasar number enhancements. This probably arises from the fact that the gravitational lensing effect by the environmental matter surrounding these objects has been ignored. In this paper we conduct an extensive study of the influence of the environmental matter on the prediction of quasar enhancement factor by employing the spatial two-point correlation function. Assuming a singular isothermal sphere for mass density profile in galaxy and cluster of galaxies, we estimate the average surface mass density $\overline{\Sigma}$ around galaxies, clusters and quasars from the galaxy-galaxy, cluster-cluster, cluster-galaxy and quasar-galaxy correlations. Our results show that the $\overline{\Sigma}$ induced quasar number enhancement in the scenario of gravitational magnification depends critically on the mass density parameters of galaxies ($\Omega_g$) and clusters of galaxies ($\Omega_c$) in the universe. For a flat cosmological model of $\Omega_0=1$ the environmental matter can indeed play an important role in the lensing origin of the quasar-quasar and quasar-galaxy associations if $\Omega_g\sim\Omega_c\sim\Omega_0$, while it is unlikely that $\overline{\Sigma}$ is sufficient to account for the reported quasar overdensity behind quasars/galaxies if galaxies and clusters of galaxies contribute no more than $25\%$ to the matter of the universe. Nonetheless, the

We present nine color CCD intermediate-band spectrophotometry of a two square degree field centered on the old open cluster M67, from 3890$\rm \AA$ to nearly 1$\mu$. These observations are taken as a part of the BATC (Beijing-Arizona-Taipei-Connecticut) Color Survey of the Sky, for both scientific and calibration reasons. With these data we show that the BATC survey can reach its goal of obtaining spectrophotometry to a zero point accuracy of 0.01 mag, and down to V = 21 with 0.3 mag random error. We fit the color-magnitude diagrams (CMDs) with Worthey's theoretical models. The net result is the excellent fit of the 4.0 Gyr, [Fe/H] = $-0.10$ model to our data, including a good fit to the main sequence (MS) turn-off. Our data are consistent with a toy model with 50\% of the stars in M67 being binaries and a random distribution of binary mass-ratios, although other models with different mass-ratio distributions cannot be ruled out. The spatial distribution and mass function (MF) of stars in M67 show marked effects of dynamical evolution and evaporation of stars from the cluster. Blue stragglers and binary stars are the most condensed within the cluster, with degree of condensation depending on mass.We find M67 to have an elongated shape, oriented at an angle of $15^{\circ}$ relative to the galactic plane. Within its tidal radius, the observed MF of M67 between 1.2 $\rm M_\odot$ and $\rm 0.8 M_\odot$ has a Salpeter slope $\rm \eta = -1.93 \pm 0.66$. For stars of mass below 0.8 $\rm M_\odot$, $\rm \eta \sim 0$. It is plausible that the leveling-off of the MF at lower masses is a result of evaporation of lower mass stars in this mass range at a rate of one every $\sim 10^7$ years. If so, it is plausible that the IMF of M67 has the canonical field value of $\rm \eta = -2.0$.

Current measurements of the Hubble constant $H_0$ on scale less than $\sim100$ Mpc appear to be controversial, while the observations made at high redshift seem to provide a relatively low value. On the other hand, the Hubble expansion is driven by the matter content of the universe. The dynamical analysis on scale of a few $\sim10$ Mpc indicates that the matter density $\Omega_0$ is only $\sim0.2$--$0.3$, which is significantly smaller than $\Omega_0=1$ predicted in the standard inflation model. This might support the tendency of a decreasing Hubble constant towards distance. In this paper, we discuss the influence of a possible variant Hubble constant on two fundamental relations in astronomy: the magnitude-redshift ($m$--$z$) and the number-magnitude relations. Using a distant type Ia supernova at $z=0.458$, we show that the deceleration parameter $q_0$ or $\Omega_0$ cannot be determined from the $m$--$z$ relation at moderate/high redshift unless the variation of the Hubble constant is {\it a priori} measured. It is further demonstrated that the number density of distant sources would be underestimated when their local calibration is employed, which may partially account for the number excess of the faint blue galaxies observed at moderate/high redshift.

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.

Motivated by the significant overdensity of background bright quasars recently detected behind the foreground clusters of galaxies on scale of $10$ arcminutes, we have investigated the possibility of attributing the quasar-cluster associations to gravitational lensing by large-scale matter inhomogeneities. Based on the conventional lensing models, we have shown that the reported quasar overdensity is unlikely to be generated by cluster matter alone. The situation does not change even if all the clusters of galaxies which follow their spatial two-point correlation function are taken into account, while matter clumps on scale of $>20$ Mpc are also found to be unable to provide the required mass surface density since their density contrast is strictly limited by the anisotropy measurements of the cosmic background radiation. Moreover, we have pointed out that the influence of a nonzero cosmological constant on the quasar-cluster associations is very minor. We conclude that either the observed quasar number counts have been seriously contaminated by the magnification bias of matter inhomogeneities of the universe or there should exists some intercluster matter on scale of less than $\sim20$ Mpc, e.g. from cluster-galaxy correlation, whose mean cosmic density is about an order of magnitude higher than that of clusters of galaxies.

Using the Boltzmann equation with a Langevin-like term describing the stochastic force in a baryon-photon plasma, we investigate the influence of the incoherent electron-photon scattering on the subhorizon evolution of the cosmic microwave radiation. The stochastic fluctuation caused by each collision on average is found to be small. Nevertheless, it leads to a significant Brownian drifting of the phase in the acoustic oscillation, and the coherent oscillations cannot be maintained during their dynamical evolution. As a consequence, the proposed Doppler peaks probably do not exist.

A method for measuring the spectrum of a density field by a discrete wavelet space-scale decomposition (SSD) has been studied. We show how the power spectrum can effectively be described by the father function coefficients (FFC) of the wavelet SSD. We demonstrate that the features of the spectrum, such as the magnitude, the index of a power law, and the typical scales, can be determined with high precision by the FFC reconstructed spectrum. This method does not require the mean density, which normally is poorly determined. The problem of the complex geometry of observed samples can also be easily solved because the basis are always orthogonal, regardless the geometry of the samples. Using this method, we examine the spectra inferred from Ly$\alpha$ forests of both simulated and real samples. We find that 1.) the magnitude of the 1-D spectra is significantly dfferent froma Poisson process; 2) the 1-d spectra are flat on scales less than 5 h$^{-1}$ Mpc, and increase slowly at larger ranges; 3.) the reconstructed 3-D spectra have about the same power as the COBE normalized linear spectrum of the SCDM model on scales less than 40 h$^{-1}$ Mpc, but is larger than the SCDM model on scales larger than 40 h$^{-1}$ Mpc; 4) the magnitudes of high redshift ($z>2.51$) spectra generally are larger than those of low redshift ($z<2.51$).

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.

A wavelet space-scale decomposition (SSD) analysis of large scale structures in the universe has been developed. The SSD method of identifying and measuring structures in the spatial distribution of objects has been demonstrated. The position and strength (richness) of the identified clusters can be described by the corresponding coefficient of the wavelet transform. Using this technique, we systematically detected the clustering and its evolution of QSO's Ly$\alpha$ forest lines in real data and simulated samples. We showed that the clusters of Ly$\alpha$ absorbers do exist on scales as large as at least 20 h$^{-1}$ Mpc at significance levels of 2-4 $\sigma$. Independent data sets show about the same strength distribution of the decomposed clusters. The number densities of the clusters on scales of 10 - 20 h$^{-1}$ Mpc are found to evolve in an opposite sense as that of the lines themselves, i.e. they decrease with redshift. We also showed that the number density and the strength distribution of clusters can play an important role in testing or discriminating models, i.e. it can distinguish real data and simulated samples, which cannot be discriminated by traditional ways. We used Daubechies 4 and Mallat wavelets as the bases of the SSD. All above-mentioned conclusions do not depend on either wavelet basis.

The possibility of thermally induced initial density perturbations in inflationary cosmology is examined. The fluctuation dynamics of a scalar field plus thermal bath system during slow roll is described by a Langevin-like equation. Fluctuation-dissipation arguments show that for a wide parameter range within the standard inflation model, the thermal fluctuations of the scalar field can dominate its quantum fluctuations. The initial amplitude of density perturbations is found to lie in a range which is consistent with the recent observations of cosmic temperature fluctuations.

Astrophysical observations indicate that the ``Local Universe" has a relatively lower matter density ($\Omega_0$) than the predictions of the standard inflation cosmology and the large-scale motions of galaxies which provide a mean mass density to be very close to unity. In such a local underdense region the Hubble expansion may not be representative of the global behaviour. Utilizing an underdense sphere embedded in a flat universe as the model of our ``Local Universe", we show that the local Hubble constant would be 1.2 -- 1.4 times larger than the global value on scale of $\sim80$ Mpc, depending on the variation of $\Omega_0$. This may account for the recent measurements of the unpleasantly large Hubble constant of $\sim$80 km/s/Mpc using the Cepheid variables in the Virgo cluster and the relative distance between Virgo and Coma cluster and removes the resulted apparent paradox of the age of our universe.

In this paper we investigate, using high resolution N-body simulations, the density profiles and the morphologies of galaxy clusters in seven models of structure formation. We show that these properties of clusters are closely related to the occurrence of a significant merging event in the recent past. The seven models are: (1) the standard CDM model (SCDM) with $\Omega_0 = 1$, $\Lambda_0=0$ and $h=0.5$; (2) a low-density flat model (FL03) with $\Omega_0=0.3$, $\Lambda_0=0.7$ and $h=0.75$; (3) an open model (OP03) with $\Omega_0=0.3$, $\Lambda_0=0$ and $h=0.75$; (4) a low-density flat model (FL02) with $\Omega_0=0.2$, $\Lambda_0=0.8$ and $h=1$; (5) an open model (OP02) with $\Omega_0=0.2$, $\Lambda_0=0$ and $h=1$; (6) a low-density flat model (FL01) with $\Omega_0=0.1$ and $\Lambda_0=0.9$; (7) an open model (OP01) with $\Omega_0=0.1$ and $\Lambda_0=0$. We find that the density profiles and morphologies of clusters depend both on $\Omega_0$ and on $\Lambda_0$. For $\Lambda_0=0$, these properties are a monotonic function of $\Omega_0$. Clusters in OP01 have the steepest density profiles, their density contours are the roundest and show the smallest center shifts. The other extreme case is SCDM, where clusters show the least steep density profiles and the most elongated contours. For a given $\Omega_0$ ($<1$), clusters in the flat model (i.e. with $\Lambda_0=1-\Omega_0$) have flatter density profiles and less substructures than in the corresponding open model. In particular, our results show that low-density flat models with $\Omega_0\sim 0.3$, which are currently considered as a successful alternative to SCDM, can produce a substantial fraction of clusters with substructures. This is in contrast to the conception that this kind of models may have serious problem in this aspect.

We investigate, using high resolution N-body simulations, the density profiles and the morphologies of galaxy clusters in three models of structure formation. We show that these properties are closely related to the occurrence of a significant merging event to a cluster in the recent past. The three models are: (1) the standard CDM model (SCDM) with $\Omega_0 = 1$, $\Lambda_0=0$ and $h=0.5$; (2) a low-density flat model (LCDM) with $\Omega_0=0.3$, $\Lambda_0=0.7$ and $h=0.75$; (3) a CDM dominated open model (OCDM) with $\Omega_0=0.1$, $\Lambda_0=0$. We find that the clusters in the OCDM model have the steepest density profiles, their density contours are the roundest and show the smallest center shifts, among the three models. The clusters in the SCDM model have contours that are most elongated and show the largest center shifts, indicating the existence of frequent mergers in the recent past. The clusters in the LCDM model have density profiles and center shifts (in their density contours) similar to those in the SCDM model, although the density contours tend to be rounder. Our results show that, although the density profiles and morphologies of clusters depend on models of structure formation, the LCDM model, which is currently considered as a successful alternative to the SCDM model, can do as well in producing a substantial fraction of clusters with substructures. This is in contrast to the conception that this model may have serious problem in this aspect.

One of the puzzles in understanding the spectra of active galactic nuclei (AGN) is the origin of the FeII emission. FeI emission, if present, will help reveal the physical conditions of the emitting gas. In an attempt to verify the presence of FeI lines, high S/N spectra of two FeII-strong quasars, IRAS 07598-6508 and PHL 1092, were obtained at the Multiple Mirror Telescope and the Steward 2.3 m Telescope. We have identified emission lines of FeI and TiII. The source of energy for FeII, FeI and TiII emission is probably not from ionization by the photon continuum, but heat. The high rate of energy generation and the presence of both high and low velocity gas indicate that the heat is generated not over a large area, but a narrow band in accretion disk, in which the rotational speed decreases rapidly.

We show that a good fitting to the first two years of {\it COBE}-DMR observations of the two-point angular correlation function of CBR temperature is given by models with a non-zero infrared cutoff $k_{min}$ in the spectrum of the primordial density perturbations. If this cutoff comes from the finiteness of the universe, say, a topological T3 model, we find $k_{min} \sim (0.3 - 1.1)\pi$H$_0/c$ with confidence level 95\%. Such a non-zero $k_{min}$ universe would also give a better match to the observations both of the RMS quadrupole anisotropy of CBR and of galaxy clustering.

Using a set of $\pppm$ simulation which accurately treats the density evolution of two components of dark matter, we study the evolution of clusters in the Cold + Hot dark matter (CHDM) model. The mass function, the velocity dispersion function and the temperature function of clusters are calculated for four different epochs of $z\le 0.5$. We also use the simulation data to test the Press-Schechter expression of the halo abundance as a function of the velocity dispersion $\sigma_v$. The model predictions are in good agreement with the observational data of local cluster abundances ($z=0$). We also tentatively compare the model with the Gunn and his collaborators' observation of rich clusters at $z\approx 0.8$ and with the x-ray luminous clusters at $z\approx 0.5$ of the {\it Einstein} Extended Medium Sensitivity Survey. The important feature of the model is the rapid formation of clusters in the near past: the abundances of clusters of $\sigma_v\ge 700\kms$ and of $\sigma_v\ge 1200 \kms$ at $z=0.5$ are only 1/4 and 1/10 respectively of the present values ($z=0$). Ongoing ROSAT and AXAF surveys of distant clusters will provide sensitive tests to the model. The abundance of clusters at $z\approx 0.5$ would also be a good discriminator between the CHDM model and a low-density flat CDM model both of which show very similar clustering properties at $z=0$.

We present results of searching for the possible typical scales in the spatial distribution of QSOs. Our method is based on the second derivative of the two-point correlation function. This statistic is sensitive to the scale of the maximum in the spectrum $P(k)$ of the density perturbation in the universe. This maximum or bend scale can be detected as the wavelengths of the periodic component in the second derivative of the integral correlation function. For various QSO samples compiled from surveys of pencil-beam and bright QSOs, a typical scale of about 93 $\pm$ 10 h$^{-1}$Mpc for $q_0=0.5$ has been detected. This typical scale is in good agreement with that found in the spatial distributions of galaxies, clusters of galaxies, and CIV absorption systems of QSOs if $q_0$ is taken to be $\sim 0.2$. Therefore, it is likely a common or universal scale in the large scale structure traced by these objects. This result is consistent with the assumption that the typical scale comes from a characteristic scale in the spectrum of the density perturbation in the universe.

It is shown that the evolution of the density perturbations during certain eras of substantial entropy generation in the universe can be described in the scheme of the KPZ equation. Therefore, the influence on cosmological structure formation by stochastic forces arising from various dissipations can be studied through the universal characteristics of surface growth in $d=3+1$ dimensions. We identify eras of strong stochastic fluctuations and describe dynamically how these other dissipative sources of noise, besides initial (inflationary) quantum fluctuations, generate seeds of density perturbation with power law spectrum, including the Harrison-Zeldovich spectrum.

We present an exact solution of the anisotropies of cosmic background radiation (CBR) from a local collapse described by a spherical over-dense region embedded in a flat universe, with the emphasis on the relationship between the dipole $(\Delta {\sf T}/{\sf T})_d$ and the quadrupole $(\Delta {\sf T}/{\sf T})_q$ anisotropy. This result has been used to examine the kinematic quadrupole correction $(\Delta {\sf T}/{\sf T})_q=(\Delta {\sf T}/{\sf T})_d^2/2$, which is usually applied to remove the contamination of the quadrupole produced by local density inhomogeneities when finding the cosmic amplitude of the quadrupole at the surface of last scattering. We have found that the quadrupole of local collapse origin cannot always be approximately described by the kinematic quadrupole. Our numerical result shows that the difference between the kinematic and local quadrupoles depends on the size and matter density in the peculiar field, and the position of the observer. For a given dipole, the local quadrupole can be different from the kinematic quadrupole by a factor as large as 3. Therefore, the kinematic quadrupole correction remains an uncertain factor in the determination of the amplitude of a cosmic quadrupole. Nevertheless, a preliminary analysis shows that this uncertainty might not dramatically change the cosmological origin of the COBE-DMR's quadrupole, unless a huge peculiar gravitational field is assumed.

Using numerical simulations, we investigate the large-scale gravitational clustering in a flat universe dominated by cold plus hot dark matter (i.e., $\Omega_0=\ocdm+\ohdm+\obaryon=1$). Primordial density fluctuation spectrum is taken to have the Zel'dovich-Harrison form. Three models are studied, with Model I having $\ocdm=0.69$, $\obaryon=0.01$, and $\ohdm=0.30$ in one flavor of neutrinos; Model II having $\ocdm=0.60$, $\obaryon=0.10$, and $\ohdm=0.30$ in one flavor of neutrinos; Model III having $\ocdm=0.69$, $\obaryon=0.01$, and $\ohdm=0.30$ in three flavors of neutrinos. The initial density spectra are normalized by the COBE quadrupole measurement, and galaxies are identified from the peaks of initial density fields above a certain threshold chosen, to match the observed two-point correlation on scales $\ls 10\mpc$. Thus the clustering properties of both the mass and the galaxies are completely specified. The biasing parameter (for the `galaxies') determined in this way is $b_g\approx 1.2$ for Model I, 1.5 for Model II and 1.6 for Model III. The clustering and motions of the simulated `galaxies' are compared with recent observations. The spatial distributions of galaxies in the hybrid models are very frothy; filaments, sheets, voids etc. of sizes 10 -- 50$\mpc$ are frequently seen in the simulations.