A large amount of observations have constrained cosmological parameters and the initial density fluctuation spectrum to a very high accuracy. However, cosmological parameters change with time and the power index of the power spectrum varies with mass scale dramatically in the so-called concordance Lambda CDM cosmology. Thus, any successful model for its structural evolution should work well simultaneously for various cosmological models and different power spectra. We use a large set of high-resolution N-body simulations of a variety of structure formation models (scale-free, standard CDM, open CDM, and Lambda CDM) to study the mass accretion histories (MAHs), the mass and redshift dependence of concentrations and the concentration evolution histories of dark matter halos. We find that there is significant disagreement between the much-used empirical models in the literature and our simulations. According to two simple but tight correlations we find from the simulation results, we develop new empirical models for both the MAHs and the concentration evolution histories of dark matter halos, and the latter can also be used to predict the mass and redshift dependence of halo concentrations. These models are accurate and universal: the same set of model parameters works well for different cosmological models and for halos of different masses at different redshifts and the model predictions are highly accurate even when the histories are traced to very high redshift. These models are also simple and easy to implement. A web calculator and a user-friendly code to make the relevant calculations are available from http://www.shao.ac.cn/dhzhao/mandc.html . We explain why Lambda CDM halos on nearly all mass scales show two distinct phases in their evolution histories.

Using a combination of N-body simulations with different resolutions, we study in detail how the concentrations of cold dark matter (CDM) halos depend on halo mass at different redshifts. We confirm that halo concentrations at the present time depend strongly on halo mass, but our results also show marked differences from the predictions of some early empirical models. Our main result is that the mass dependence of the concentrations becomes weaker at higher redshifts, and at z >~ 3 halos of mass greater than 10^{11} Msun/h all have a similar median concentration, c ~ 3.5. While the median concentrations of low-mass halos grow significantly with time, those of massive halos change only little with redshifts. These results are quantitatively in good agreement with the empirical model proposed by Zhao et al. which shows that halos in the early fast accretion phase all have similar concentrations.

By comparing semi-analytic galaxy catalogues with data from the Sloan Digital Sky Survey (SDSS), we show that current galaxy formation models reproduce qualitatively the dependence of galaxy clustering and pairwise peculiar velocities on luminosity, but some subtle discrepancies with the data still remain. The comparisons are carried out by constructing a large set of mock galaxy redshift surveys that have the same selection function as the SDSS Data Release Four (DR4). The mock surveys are based on two sets of semi-analytic catalogues presented by Croton et al. and Kang et al. From the mock catalogues, we measure the redshift-space projected two-point correlation function, the power spectrum, and the pairwise velocity dispersion (PVD) in Fourier space and in configuration space, for galaxies in different luminosity intervals. We then compare these theoretical predictions with the measurements derived from the SDSS DR4. On large scales and for galaxies brighter than L*, both sets of mock catalogues agree well with the data. For fainter galaxies, however, both models predict stronger clustering and higher pairwise velocities than observed. We demonstrate that this problem can be resolved if the fraction of faint satellite galaxies in massive haloes is reduced by ~30% compared to the model predictions. A direct look into the model galaxy catalogues reveals that a signifcant fraction (15%) of faint galaxies ($-18<M_{^{0.1}r}<-17$) reside in haloes with $M_{vir}>10^{13}\msun$, and this population is predominantly red in colour. These faint red galaxies are responsible for the high PVD values of low-luminosity galaxies on small scales.

(abridged) We present measurements of the pairwise velocity dispersion (PVD) for different classes of galaxies in the Sloan Digital Sky Survey. For a sample of about 200,000 galaxies, we study the dependence of the PVD on galaxy properties such as luminosity, stellar mass (M_*), colour (g-r), 4000A break strength (D4000), concentration index (C), and stellar surface mass density (\mu_*). The luminosity dependence of the PVD is in good agreement with the results of Jing & B\"orner (2004) for the 2dFGRS catalog. The value of \sigma_{12} measured at k=1 h/Mpc decreases as a function of increasing galaxy luminosity for galaxies fainter than L*, before increasing again for the most luminous galaxies in our sample. Each of the galaxy subsamples selected according to luminosity or stellar mass is divided into two further subsamples according to colour, D4000, C and \mu_*. We find that galaxies with redder colours and higher D4000, C, and \mu_* values have larger PVDs on all scales and at all luminosities/stellar masses. The dependence of the PVD on parameters related to recent star formation(colour, D4000) is stronger than on parameters related to galaxy structure (C, \mu_*), especially on small scales and for faint galaxies. The reddest galaxies and galaxies with high surface mass densities and intermediate concentrations have the highest pairwise peculiar velocities, i.e. these move in the strongest gravitational fields. We conclude that the faint red population located in rich clusters is responsible for the high PVD values that are measured for low-luminosity galaxies on small scales.

(abridged)We use a sample of ~200,000 galaxies drawn from the Sloan Digital Sky Survey to study how clustering depends on properties such as stellar mass (M*), colour (g-r), 4000A break strength (D4000), concentration index (C), and stellar surface mass density (\mu_*). We find that more massive galaxies cluster more strongly than less massive galaxies, with the difference increasing above the characteristic stellar mass of the Schechter mass function. When divided by physical quantities, galaxies with redder colours, larger D4000, higher C and larger \mu_* cluster more strongly. The clustering differences are largest on small scales and for low mass galaxies. At fixed stellar mass,the dependences of clustering on colour and 4000A break strength are similar. Different results are obtained when galaxies are split by concentration or surface density. The dependence of w(r_p) on g-r and D4000 extends out to physical scales that are significantly larger than those of individual dark matter haloes (> 5 Mpc/h). This large-scale clustering dependence is not seen for the parameters C or \mu_*. On small scales (< 1 Mpc/h), the amplitude of the correlation function is constant for ``young'' galaxies with 1.1 < D4000< 1.5 and a steeply rising function of age for ``older'' galaxies with D4000>1.5. In contrast, the dependence of the amplitude of w(r_p) on concentration on scales less than 1 Mpc/h is strongest for disk-dominated galaxies with C<2.6. This demonstrates that different processes are required to explain environmental trends in the structure and in star formation history of galaxies.

We model the galaxy formation in a series of high-resolution N-body simulations using the semi-analytical approach. Unlike many earlier investigations based on semi-analytical models, we make use of the subhalos resolved in the $N$-body simulations to follow the mergers of galaxies in dark halos, and we show that this is pivotal in modeling correctly the galaxy luminosity function at the bright end and the bimodal nature of galaxy color distribution. Mergers of galaxies based on subhalos also result in many more bright red galaxies at high $z$. The semi-analytical model we adopt is similar to those used in earlier semi-analytical studies, except that we consider the effect of a prolonged cooling in small halos and that we explicitly follow the chemical enrichment in the interstellar medium. We use our model to make predictions for the properties of the galaxy population at low redshift and compare them with various current observations. We find that our model predictions can match the luminosity functions of galaxies in various wavebands redder than the u-band. The shape of the luminosity function at bright end is well reproduced if galaxy mergers are modeled with the merger trees of subhalos and the steep faint-end slope can be moderated if the gas cooling time in low-mass halos is comparable to the age of the universe. The model with subhalos resolved can reproduce the main features in the observed color bimodal distribution, though it still predicts too many bright blue galaxies. The same model can also match the color-magnitude relation for elliptical galaxies in clusters, the metallicity-luminosity relation.

We present the first determination of the pairwise velocity dispersion (PVD) for galaxies in different luminosity intervals using the final release of the Two-Degree Field Galaxy Redshift Survey (2dFGRS). We have discovered quite surprisingly that the relative velocities of the faint galaxies at small separation are very high, around $700 \kms$, reaching similar values as the brightest galaxies. At intermediate luminosities $M^*-1$ ($M^*$ is the characteristic luminosity of the Schechter function), the relative velocities exhibit a well defined steep minimum near $400 \kms$. This result has been derived using a novel method to determine the real space power spectrum and the PVD from the redshift space power spectrum of the 2dFGRS. Both quantities can be determined quite reliably. We have taken the parameter $\beta$ equal to 0.45. But we have also checked that the results are changed very little, if we allow a variation of $\beta$ with luminosity. The result is a challenge to the current halo model of galaxies of Yang et al. that was obtained by matching the clustering and luminosity function of the 2dFGRS, but cannot reproduce the luminosity dependence of the PVD. It may also be an mportant constraint in general on theories of galaxy formation, such as semi-analytical models and hydro/N-body simulations of galaxy formation. Combined with the observed luminosity dependence of clustering, our result implies that a substantial fraction of faint galaxies, as well as the brightest ones, are in massive halos of galaxy cluster size, but most of the $M^*$ galaxies are in galactic halos.

We present the first determination of the pairwise velocity dispersion of galaxies at different luminosity with the final release of the Two-Degree Field Galaxy Redshift Survey (2dFGRS). Our result surprisingly shows that the random velocities of the faint galaxies are very high, around $ 700 \kms$, reaching similar values as the brightest galaxies. At intermediate luminosities slightly brighter than the characteristic luminosity $M_\star$, the velocities exhibit a well defined steep minimum near $ 400 \kms$. The result challenges the current halo model of galaxies of Yang et al. that was obtained by matching the clustering and luminosity function of 2dFGRS, and can be an important constraint in general on theories of galaxy formation, e.g., the semi-analytical model. Combining the observed luminosity dependence of clustering, our result implies that quite a fraction of faint galaxies are in massive halos of galaxy clusters as the brightest ones, but most of the $M_\star$ galaxies are in galactic halos.

In a detailed analysis of the three point correlation function (3PCF) for the 2dF Galaxy Redshift Survey we have accurately measured the 3PCF for galaxies of different luminosity. The 3PCF amplitudes [$\Qsu$ or $\Qrpu$] of the galaxies generally decrease with increasing triangle size and increase with the shape parameter $v$, in qualitative agreement with the predictions for the clustering of dark matter in popular hierarchical CDM models. The 2dFGRS results agree well with the results of Jing & B\"orner for the Las Camapanas Redshift Survey (LCRS), though the measurement accuracy is greatly improved in the present study. The dependence of the 3PCF on luminosity is not significant, but there seems to be a trend for the brightest galaxy sample to have a lower amplitude than the fainter ones. Comparing the measured 3PCF amplitudes [$\Qsu$ or $\Qrpu$] to the prediction of a WMAP concordance model, we find that the measured values are consistently lower than the predicted ones for dark matter. This is most pronounced for the brightest galaxies, for which about one-half of the predicted $Q$ value provides a good description of $\Qrpu$ for the 2dFGRS data. For a less luminous sample, the $Q$ values are also smaller than in the dark matter model on small scales, but on scales larger than $s=8 \mpc$ and $r_p=3.25 \mpc$ they reach the model values. We discuss implications for current theories of galaxy formation.

Using high--resolution N--body simulations, we test two theoretical models, based either on spherical or on ellipsoidal collapse model, for the higher--order moments of the dark matter halo distribution in CDM models. We find that a theoretical model based on spherical collapse describes accurately the simulated counts--in--cells moments for haloes of several mass ranges. It appears that the model using ellipsoidal collapse instead of spherical collapse in defining dark haloes is unable to improve the models for the higher--order moments of halo distribution, for haloes much smaller than $M^*$ (the mass scale on which the fluctuation of the density field has a rms about 1). Both models are particularly accurate for the descendants of haloes selected at high redshift, and so are quite useful in interpreting the high--order moments of galaxies. As an application we use the theoretical model to predict the higher--order moments of the Lyman break galaxies observed at $z\approx 3$ and their descendants at lower redshifts.

In this paper, we analyse in detail the mass-accretion histories and structural properties of dark haloes in high-resolution N-body simulations. Modeling the density distribution in individual haloes with the NFW profile, we find, for all main progenitors of a given halo, there is a tight correlation between its inner scale radius $r_s$ and the mass within it, $M_s$, which is the basic reason why halo structural properties are closely related to their mass-accretion histories. This correlation can be used to predict accurately the structural properties of a dark halo at any time from its mass-accretion history. We also test our model with a large sample of GIF haloes. The build-up of dark haloes in CDM models generally consists of an early phase of fast accretion and a late phase of slow accretion [where $M_h$ increases with time approximately as the expansion rate]. These two phases are separated at a time when the halo concentration parameter $c\sim 4$. Haloes in the two accretion phases show systematically different properties, for example, the circular velocity $v_h$ increases rapidly with time in the fast accretion phase but remain almost constant in the slow accretion phase,the inner properties of a halo, such as $r_s$ and $M_s$ increase rapidly with time in the fast accretion phase but change only slowly in the slow accretion phase. The potential well associated with a halo is built up mainly in the fast accretion phase, even though a large amount of mass (over 10 times) can be accreted in the slow accretion phase. We discuss our results in connection to the formation of dark haloes and galaxies in hierarchical models.

We investigate the feasibility of determining the pairwise velocity dispersion (PVD) for Lyman Break Galaxies (LBGs), and of using this quantity as a discriminator among theoretical models. We find that when the Two Point Correlation Function over the range of observable separations is fixed, different schemes of galaxy formation will lead to significant changes in PVD, while influence of different cosmogonies is invisible. Our test with mock samples shows furthermore that one can discriminate among such models already with currently available observational samples (if the measurement error of the redshift is negligible) which have a typical error of $80\kms$. The error will be reduced by a factor of 2 if the samples are increased four times. We also show that an erroneous assumption about the geometry of the universe and different infall models only slightly change the results. Therefore the PVD will become another promising statistic to test galaxy formation models with redshift samples of LBGs.

The stochasticity in the distribution of dark haloes in the cosmic density field is reflected in the distribution function $P_V(N_h|\delta_m)$ which gives the probability of finding $N_h$ haloes in a volume $V$ with mass density contrast $\delta_m$. We study the properties of this function using high-resolution $N$-body simulations, and find that $P_V(N_n|\delta_m)$ is significantly non-Poisson. The ratio between the variance and the mean goes from $\sim 1$ (Poisson) at $1+\delta_m\ll 1$ to $<1$ (sub-Poisson) at $1+\delta_m\sim 1$ to $>1$ (super-Poisson) at $1+\delta_m\gg 1$. The mean bias relation is found to be well described by halo bias models based on the Press-Schechter formalism. The sub-Poisson variance can be explained as a result of halo-exclusion while the super-Poisson variance at high $\delta_m$ may be explained as a result of halo clustering. A simple phenomenological model is proposed to describe the behavior of the variance as a function of $\delta_m$. Galaxy distribution in the cosmic density field predicted by semi-analytic models of galaxy formation shows similar stochastic behavior. We discuss the implications of the stochasticity in halo bias to the modelling of higher-order moments of dark haloes and of galaxies.

We report on the measurement of the two-point correlation function and the pairwise peculiar velocity of galaxies in the IRAS PSCz survey. The real space two-point correlation function can be fitted to a power law $\xi(r) = (r_0/r)^{\gamma}$ with $\gamma=1.69$ and $r_0=3.70 \mpc$. The pairwise peculiar velocity dispersion $\sigma_{12}(r_p)$ is close to $400 \kms$ at $r_p=3\mpc$ and decreases to about $150 \kms$ at $r_p \approx 0.2 \mpc$. These values are significantly lower than those obtained from the Las Campanas Redshift Survey, but agree very well with the results of blue galaxies reported by the SDSS team later on. We have constructed mock samples from N-body simulations with a cluster-weighted bias and from the theoretically constructed GIF catalog. We find that the two-point correlation function of the mock galaxies can be brought into agreemnt with the observed result, but the model does not reduce the velocity dispersions of galaxies to the level measured in the PSCz data. Thus we conclude that the peculiar velocity dispersions of the PSCz galaxies require a biasing model which substantially reduces the peculiar velocity dispersion on small scales relative to their spatial clustering. The results imply that either the cosmogony model needs to be revised or the velocity bias is important for the velocity dispersion of the IRAS galaxies.

We report on the measurement of the two-point correlation function, and the pairwise peculiar velocity of galaxies in the IRAS PSCz survey. We compute these statistics first in redshift space, and then obtain the projected functions which have simple relations to the real-space correlation functions on the basis of the method developed earlier in analyzing the Las Campanas Redshift Survey (LCRS) by Jing, Mo, & B\"orner (1998). We find that the real space two-point correlation function can be fitted to a power law $\xi(r) = (r_0/r)^{\gamma}$ with $\gamma=1.69$ and $r_0=3.70 \mpc$. The pairwise peculiar velocity dispersion $\sigma_{12}(r_p)$ is close to $400 \kms$ at $r_p=3\mpc$ and decreases to about $150 \kms$ at $r_p \approx 0.2 \mpc$. These values are significantly lower than those obtained from the LCRS. In order to understand the implications of those measurements on the galaxy biasing, we construct mock samples for a low density spatially-flat cold dark matter model ($\Omega_0 = 0.3$, $\lambda_0=0.7$, $\Gamma=0.2$, $\sigma_8=1$) using a set of high-resolution N-body simulations. Applying a stronger cluster-underweight biasing ($\propto M^{-0.25}$) than for the LCRS ($\propto M^{-0.08}$), we are able to reproduce these observational data, except for the strong decrease of the pairwise peculiar velocity at small separations. This is qualitatively ascribed to the different morphological mixture of galaxies in the two catalogues. Disk-dominated galaxy samples drawn from the theoretically constructed GIF catalog yield results rather similar to our mock samples with the simple cluster-underweight biasing.

In a recent paper we have studied the redshift power spectrum $P^S(k,\mu)$ in three CDM models with the help of high resolution simulations. Here we apply the method to the largest available redshift survey, the Las Campanas Redshift Survey (LCRS). The basic model is to express $P^S(k,\mu)$ as a product of three factors P^S(k,\mu)=P^R(k)(1+\beta\mu^2)^2 D(k,\mu). Here $\mu$ is the cosine of the angle between the wave vector and the line of sight. The damping function $D$ for the range of scales accessible to an accurate analysis of the LCRS is well approximated by the Lorentz factor D=[1+{1\over 2}(k\mu\sigma_{12})^2]^{-1}. We have investigated different values for $\beta$ ($\beta=0.4$, 0.5, 0.6), and measured $P^R(k)$ and $\sigma_{12}(k)$ from $P^S(k,\mu)$ for different values of $\mu$. The velocity dispersion $\sigma_{12}(k)$ is nearly a constant from $k=0.5$ to 3 $\mpci$. The average value for this range is $510\pm 70 \kms$. The power spectrum $P^R(k)$ decreases with $k$ approximately with $k^{-1.7}$ for $k$ between 0.1 and 4 $\mpci$. The statistical significance of the results, and the error bars, are found with the help of mock samples constructed from a large set of high resolution simulations. A flat, low-density ($\Omega_0=0.2$) CDM model can give a good fit to the data, if a scale-dependent special bias scheme is used which we have called the cluster-under-weighted bias (Jing et al.).

We report the results of an analysis of the redshift power spectrum $P^S(k,\mu)$ in three typical Cold Dark Matter (CDM) cosmological models, where $\mu$ is the cosine of the angle between the wave vector and the line-of-sight. Two distinct biased tracers derived from the primordial density peaks of Bardeen et al. and the cluster-underweight model of Jing, Mo, & B\"orner are considered in addition to the pure dark matter models. Based on a large set of high resolution simulations, we have measured the redshift power spectrum for the three tracers from the linear to the nonlinear regime. We investigate the validity of the relation - guessed from linear theory - in the nonlinear regime $$ P^S(k,\mu)=P^R(k)[1+\beta\mu^2]^2D(k,\mu,\sigma_{12}(k)), $$ where $P^R(k)$ is the real space power spectrum, and $\beta$ equals $\Omega_0^{0.6}/b_l$. The damping function $D$ which should generally depend on $k$, $\mu$, and $\sigma_{12}(k)$, is found to be a function of only one variable $k\mu\sigma_{12}(k)$. This scaling behavior extends into the nonlinear regime, while $D$ can be accurately expressed as a Lorentz function - well known from linear theory - for values $D > 0.1$. The difference between $\sigma_{12}(k)$ and the pairwise velocity dispersion defined by the 3-D peculiar velocity of the simulations (taking $r=1/k$) is about 15%. Therefore $\sigma_{12}(k)$ is a good indicator of the pairwise velocity dispersion. The exact functional form of $D$ depends on the cosmological model and on the bias scheme. We have given an accurate fitting formula for the functional form of $D$ for the models studied.

In this paper we present Monte-Carlo simulations of Lyman alpha absorption systems which originate in galactic haloes, galaxy discs and dark matter (DM) satellites around big central haloes. It is found that for strong Lyman alpha absorption lines galactic haloes and satellites can explain ~20% and 40% of the line number density of QSO absorption line key project respectively. If big galaxies indeed possess such large numbers of DM satellites and they possess gas, these satellites may play an important role for strong Lyman alpha lines. However the predicted number density of Lyman-limit systems by satellites is \~0.1 (per unit redshift), which is four times smaller than that by halo clouds. Including galactic haloes, satellites and HI discs of spirals, the predicted number density of strong lines can be as much as 60% of the HST result. The models can also predict all of the observed Lyman-limit systems. The average covering factor within 250 kpc/h is estimated to be ~0.36. And the effective absorption radius of a galaxy is estimated to be ~150 kpc/h. The models predict W_r propto rho^{-0.5} L_B^{0.15} (1+z)^{-0.5}. We study the selection effects of selection criteria similar to the imaging and spectroscopic surveys. We simulate mock observations through known QSO lines-of-sight and find that selection effects can statistically tighten the dependence of line width on projected distance. (abridged)

We study the reliability of the reconstruction method which uses a modelling of the redshift distortions of the two-point correlation function to estimate the pairwise peculiar velocity dispersion of galaxies. In particular, the dependence of this quantity on different models for the infall velocity is examined for the Las Campanas Redshift Survey. We make extensive use of numerical simulations and of mock catalogs derived from them to discuss the effect of a self-similar infall model, of zero infall, and of the real infall taken from the simulation. The implications for two recent discrepant determinations of the pairwise velocity dispersion for this survey are discussed.

We report the measurement of the three-point correlation function (3PCF) of galaxies for the Las Campanas Redshift Survey (LCRS). We have not only measured the 3PCF in redshift space but also developed a method to measure the projected 3PCF which has simple relations to the real space 3PCF. Both quantities have been measured as a function of triangle size and shape with only a fractional uncertainty in each individual bin. Various tests derived from mock catalogs have been carried out to assure that the measurement is stable and that the errors are estimated reliably. Our results indicate that the 3PCFs both in redshift space and in real space have small but significant deviations from the well-known hierarchical form. The 3PCF in redshift space can be fitted by $\Qsu=0.5\cdot 10^{[0.2+0.1({s\over s+1})^2]v^2}$ for $0.8<s_{12}<8 \mpc$ and $s_{31}<16\mpc$, and the projected 3PCF by $\Qrpu=0.7r_p^{-0.3}$ for $0.2<r_{p12}<3\mpc$ and $r_{p31}<6\mpc$ ($s$ and $r_p$ are in unit of $\mpc$), though a systematic weak increase of $\Qrpu$ with $v$ at $r_p>1\mpc$ is noted. The real-space $\Qru$ for $0.2 \ls r_{12} \ls 3\mpc$ and $r_{31}\ls 6\mpc$ can be well described by half the mean 3PCF predicted by a CDM model with $\Omega_0 h=0.2$. Implications for the CDM cosmogonic models are discussed in detail with the help of mock simulation samples.

We show, with the help of large N-body simulations, that the real-space two-point correlation function and pairwise velocity dispersion of galaxies can both be measured reliably from the Las Campanas Redshift Survey. The real-space correlation function is well fitted by the power law $\xi(r)=(r_0/r)^\gamma$ with $r_0=(5.06\pm0.12)\mpc$ and $\gamma=1.862\pm 0.034$, and the pairwise velocity dispersion at $1\mpc$ is $(570\pm 80)\kms$. A detailed comparison between these observational results and the predictions of current CDM cosmogonies is carried out. We construct 60 mock samples for each theoretical model from a large set of high resolution N-body simulations, which allows us to include various observational selection effects in the analyses and to use exactly the same methods for both real and theoretical samples. We demonstrate that such a procedure is essential in the comparison between models and observations. The observed two-point correlation function is significantly flatter than the mass correlation function in current CDM models on scales $\la 1\mpc$. The observed pairwise velocity dispersion is also lower than that of dark matter particles in these models. We propose a simple antibias model to explain these discrepancies. This model assumes that the number of galaxies per unit dark matter mass, $N/M$, decreases with the mass of dark haloes. The predictions of CDM models with $\sigma_8\Omega_0^{0.6}\sim 0.4$-0.5 and $\Omega_0 h\sim 0.2$ are in agreement with the observational results, if the trend of $N/M$ with $M$ is at the level already observed for rich clusters of galaxies. Thus CDM models with cluster-abundance normalization are consistent with the observed correlation function and pairwise velocity dispersion of galaxies. A high level of velocity bias is not required in these models.

Using the second-order Eulerian perturbation theory (SEPT), we study the three-point correlation function $\zeta$ in the quasilinear regime for the SCDM, LCDM and MDM models, with the interesting result that these three models have distinctive three-point correlation functions. We test this SEPT prediction using a large set of high-resolution N-body simulations. The N-body results show that the SEPT prediction for $\zeta$ is not accurate even in the quasilinear regime ($r\la 10 \mpc$), in contrast to previous N-body tests on the skewness. However, similar to the perturbation theory, our N-body results still predict a strong dependence of the three-point correlation on the triangle shape which is observable in the distribution of galaxies if the galaxies trace the distribution of the underlying mass

The gravitational instability of inhomogeneities in the expanding universe is studied by the relativistic second-order approximation. Using the tetrad formalism we consider irrotational dust universes and get equations very similar to those given in the Lagrangian perturbation theory in Newtonian cosmology. Neglecting the cosmological constant and assuming a flat background model we give the solutions of the nonlinear dynamics of cosmological perturbations. We present the complete second-order solutions, which extend and improve earlier works.

In this paper we investigate the velocity dispersion profiles of clusters of galaxies for seven cosmological models. One model is the SCDM model, and the others are six low-density models with the density parameter $\Omega=0.1$, 0.2 or 0.3 and with or without a cosmological constant $\Lambda=1-\Omega$. We find that the velocity dispersion profiles depend both on $\Omega$ and on $\Lambda$. For $\Lambda=0$, the profiles are steeper in a lower-$\Omega$ model than in a higher-$\Omega$ one. The cosmological constant significantly weakens the dependence on $\Omega$: the difference in the profile distributions between two flat models is much smaller than that between the two corresponding open models with the same $\Omega$. These results in principle can be used to constrain the cosmological parameters when a large sample of the velocity dispersion profiles is available. Motivated by the practical situation that a sample of $\sim 100$ clusters with $\sim 100$ measured redshifts per cluster is still the best sample available in the foreseen future, we examine carefully to what degree the cosmological parameters can be constrained with the velocity dispersion profiles of such a sample of clusters. The limited sampling around clusters and the limited number of clusters seriously degrade the discriminative power of the velocity dispersion profiles among cosmological models. We find that the five models of $\Omega\ge 0.2$ cannot be distinguished by this type of observation. Due to the limited sampling, one should be very cautious in extracting information about the density profile and/or the dynamics around single clusters from the diluted velocity dispersion profiles

In this paper, we use large $\pppm$ N-body simulations to study the three-point correlation function $\zeta$ of clusters in two theoretical models. The first model (LCDM) is a low-density flat model of $\Omega_0=0.3$, $\Lambda_0=0.7$ and $h=0.75$, and the second model (PIM) is an Einstein-de-Sitter model with its linear power spectrum obtained from observations. We found that the scaled function $Q(r,u,v)$, which is defined as the ratio of $\zeta (r, ru, ru+rv)$ to the hierarchical sum $\xi (r)\xi (ru)+ \xi (ru) \xi (ru+rv) +\xi (ru+rv)\xi (r)$ (where $\xi$ is the two-point correlation function of clusters), depends weakly on $r$ and $u$, but very strongly on $v$. $Q(r,u,v)$ is about 0.2 at $v=0.1$ and 1.8 at $v=0.9$. A model of $Q(r,u,v)=\Theta 10^{1.3v^2}$ can fit the data of $\zeta$ very nicely with $\Theta\approx 0.14$. This model is found to be universal for the LCDM clusters and for the PIM clusters. Furthermore, $Q(r,u,v)$ is found to be insensitive to the cluster richness. We have compared our N-body results with simple analytical theories of cluster formation, like the peak theories or the local maxima theories. We found that these theories do not provide an adequate description for the three-point function of clusters. We have also examined the observational data of $\zeta$ presently available, and have not found any contradiction between the observations and our model predictions. The $v$-dependence of $q$ in a projected catalogue of clusters is shown to be much weaker than the $v$-dependence of $Q$ found in the three-dimensional case. It would be important to search for the $v$-dependence of $Q$ in redshift samples of rich clusters.

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.

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

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.