It is believed that massive galaxies have quenched their star formation because of active galactic nucleus feedback. However, recent studies have shown that some massive galaxies are still forming stars. We analyze the morphology of star formation regions for galaxies of stellar mass larger than 10$^{11.3}$ M$_{\odot}$ at around redshift $z_r=0.5$ using $u-z$ color images. We find that about $20\%$ of the massive galaxies are star-forming (SF) galaxies, and most of them ($\sim 85\%$) have asymmetric structures induced by recent mergers. Moreover, for these asymmetric galaxies, we find that the asymmetry of the SF regions becomes larger for bluer galaxies. Using the Illustris simulation, we can qualitatively reproduce the observed relation between asymmetry parameter and color. Furthermore, using the merger trees in the simulation, we find a correlation between the color of the main branch galaxies at $z_r=0.5$ and the sum of the Star Formation Rates (SFRs) of the recently accreted galaxies, which implies that star formation of the accreted galaxies has contributed to the observed star formation of the massive (host) galaxies (ex situ star formation). Furthermore, we find two blue and symmetric galaxies, candidates for massive blue disks, in our observed sample, which indicates that about $\sim 10\%$ of massive SF galaxies are forming stars in the normal mode of disk star formation (in situ star formation). With the simulation, we find that the disk galaxies at $z_r \approx 0.5$ should have experienced few major mergers during the last 4.3 Gyrs.

The 21~cm intensity mapping experiments promise to measure the baryon acoustic oscillations (BAO) through the mapping of H{\small I} gas at the post-reionization epoch. Its success depends on the assumption that H{\small I} gas traces the matter density distribution {\it linearly} on large scales. However, both nonlinear halo clustering and nonlinear effects modulating H{\small I} gas in halos may spoil this. We employ three approaches to generate the mock H{\small I} density from a large-scale N-body simulation at low redshifts, and demonstrate that the assumption of H{\small I} linearity is valid at the scale corresponding to the first BAO peak, but breaks down at $k \gtrsim 0.1\,h\, {\rm Mpc}^{-1}$. The nonlinear effects of halo clustering and H{\small I} content modulation counteract each other at small scales, and their competition results in a model-dependent "sweet-spot" redshift near $z$=1 where the H{\small I} bias is scale-independent down to small scales. We also find that the linear H{\small I} bias scales approximately linearly with redshift for $z\le 3$.

We estimate the Milky Way (MW) halo properties using satellite kinematic data including the latest measurements from Gaia DR2. With a simulation-based 6D phase-space distribution function (DF) of satellite kinematics, we can infer halo properties efficiently and without bias, and handle the selection function and measurement errors rigorously in the Bayesian framework. Applying our DF from the EAGLE simulation to 28 satellites, we obtain an MW halo mass of $M=1.23_{-0.18}^{+0.21}\times 10^{12} M_\odot$ and a concentration of $c=9.4_{ -2.1}^{ +2.8}$ with the prior based on the $M$-$c$ relation. The inferred mass profile is consistent with previous measurements but with better precision and reliability due to the improved methodology and data. Potential improvement is illustrated by combining satellite data and stellar rotation curves. Using our EAGLE DF and best-fit MW potential, we provide much more precise estimates of kinematics for those satellites with uncertain measurements. Compared to the EAGLE DF, which matches the observed satellite kinematics very well, the DF from the semi-analytical model based on the dark-matter-only simulation Millennium II (SAM-MII) over-represents satellites with small radii and velocities. We attribute this difference to less disruption of satellites with small pericenter distances in the SAM-MII simulation. By varying the disruption rate of such satellites in this simulation, we estimate a $\sim 5\%$ scatter in the inferred MW halo mass among hydrodynamics-based simulations.

Chiral symmetry is maximally violated in weak interactions, and such microscopic asymmetries in the early Universe might leave observable imprints on astrophysical scales without violating the cosmological principle. In this Letter, we propose a helicity measurement to detect primordial chiral violation. We point out that observations of halo-galaxy angular momentum directions (spins), which are frozen in during the galaxy formation process, provide a fossil chiral observable. From the clustering mode of large scale structure of the Universe, we construct a spin mode in Lagrangian space and show in simulations that it is a good probe of halo-galaxy spins. In standard model, a strong symmetric correlation between the left and right helical components of this spin mode and galaxy spins is expected. Measurements of these correlations will be sensitive to chiral breaking, providing a direct test of chiral symmetry breaking in the early Universe.

Galaxy clustering provides insightful clues to our understanding of galaxy formation and evolution, as well as the universe. The redshift assignment for the random sample is one of the key steps to measure the galaxy clustering accurately. In this paper, by virtue of the mock galaxy catalogs, we investigate the effect of two redshift assignment methods on the measurement of galaxy two-point correlation functions (hereafter 2PCFs), the Vmax method and the "shuffled" method. We found that the shuffled method significantly underestimates both of the projected 2PCFs and the two-dimensional 2PCFs in redshift space. While the Vmax method does not show any notable bias on the 2PCFs for volume-limited samples. For flux-limited samples, the bias produced by the Vmax method is less than half of the shuffled method on large scales. Therefore, we strongly recommend the Vmax method to assign redshifts to random samples in the future galaxy clustering analysis.

We propose a versatile and accurate method to estimate the halo mass and concentration from the kinematics of satellite galaxies. We construct the 6D phase-space distribution function of satellites from a cosmological simulation based on the similarity of internal dynamics for different halos. Within the Bayesian statistical framework, not only can we infer the halo mass and concentration efficiently, but also treat various observational effects, including the selection function, incomplete data, and measurement errors, in a rigorous and straightforward manner. Through tests with mock samples, we show that our method is valid and accurate, and more precise than pure steady-state methods. It can constrain the halo mass to within ~ 20% using only 20 tracers and has a small intrinsic uncertainty of ~ 10%. In addition to the clear application to the Milky Way and similar galaxies, our method can be extended to galaxy groups or clusters.

The CFHT Large Area U-band Deep Survey (CLAUDS) uses data taken with the MegaCam mosaic imager on CFHT to produce images of 18.60 deg2 with median seeing of FWHM=0.92 arcsec and to a median depth of U = 27.1 AB (5 sigma in 2 arcsec apertures), with selected areas that total 1.36 deg2 reaching a median depth of U=27.7 AB. These are the deepest U-band images assembled to date over this large an area. These data are located in four fields also imaged to comparably faint levels in grizy and several narrowband filters as part of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (HSC-SSP). These CFHT and Subaru datasets will remain unmatched in their combination of area and depth until the advent of the Large Synoptic Survey Telescope (LSST). This paper provides an overview of the scientific motivation for CLAUDS and gives details of the observing strategy, observations, data reduction, and data merging with the HSC-SSP. Three early applications of these deep data are used to illustrate the potential of the dataset: deep U-band galaxy number counts, z~3 Lyman break galaxy (LBG) selection, and photometric redshifts improved by adding CLAUDS U to the Subaru HSC grizy photometry.

Low mass galaxies are expected to be dark matter dominated even within their centrals. Recently two observations reported two dwarf galaxies in group environment with very little dark matter in their centrals. We explore the population and origins of dark-matter-deficient galaxies (DMDGs) in two state-of-the-art hydrodynamical simulations, the EAGLE and Illustris projects. For all satellite galaxies with $10^9<M_*<10^{10}$ M$_{\odot}$ in groups with $M_{200}>10^{13}$ M$_{\odot}$, we find that about $2.6\%$ of them in the EAGLE, and $1.5\%$ in the Illustris are DMDGs with dark matter fractions below $50\%$ inside two times half-stellar-mass radii. We demonstrate that DMDGs are highly tidal disrupted galaxies; and because dark matter has higher binding energy than stars, mass loss of the dark matter is much more rapid than stars in DMDGs during tidal interactions. If DMDGs were confirmed in observations, they are expected in current galaxy formation models.

We have obtained three-dimensional maps of the universe in $\sim200\times200\times80$ comoving Mpc$^3$ (cMpc$^3$) volumes each at $z=5.7$ and $6.6$ based on a spectroscopic sample of 179 galaxies that achieves $\gtrsim80$\% completeness down to the Ly$\alpha$ luminosity of $\log(L_{\rm Ly\alpha}/[\mathrm{erg\ s^{-1}}])=43.0$, based on our Keck and Gemini observations and the literature. The maps reveal filamentary large-scale structures and two remarkable overdensities made out of at least 44 and 12 galaxies at $z=5.692$ (z57OD) and $z=6.585$ (z66OD), respectively, making z66OD the most distant overdensity spectroscopically confirmed to date with $>10$ spectroscopically confirmed galaxies. We compare spatial distributions of submillimeter galaxies at $z\simeq 4-6$ with our $z=5.7$ galaxies forming the large-scale structures, and detect a $99.97\%$ signal of cross correlation, indicative of a clear coincidence of dusty star-forming galaxy and dust unobscured galaxy formation at this early epoch. The galaxies in z57OD and z66OD are actively forming stars with star formation rates (SFRs) $\gtrsim5$ times higher than the main sequence, and particularly the SFR density in z57OD is 10 times higher than the cosmic average at the redshift (a.k.a. the Madau-Lilly plot). Comparisons with numerical simulations suggest that z57OD and z66OD are protoclusters that are progenitors of the present-day clusters with halo masses of $\sim10^{14}\ \mathrm{M_\odot}$.

We develop a novel method to explore the galaxy-halo connection using the galaxy imaging surveys by modeling the projected two-point correlation function measured from the galaxies with reasonable photometric redshift measurements. By assuming a Gaussian form of the photometric redshift errors, we are able to simultaneously constrain the halo occupation distribution (HOD) models and the effective photometric redshift uncertainties. Tests with mock galaxy catalogs demonstrate that this method can successfully recover (within $\sim 1\sigma$) the intrinsic large-scale galaxy bias, as well as the HOD models and the effective photometric redshift uncertainty. This method also works well even for galaxy samples with 10 per cent catastrophic photometric redshift errors.

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles.

We present optical luminosity functions (LFs) of galaxies in the $^{0.1}g$, $^{0.1}r$, and $^{0.1}i$ bands, calculated using data in $\sim 40$ $deg^{2}$ sky area of LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in Southern Galactic Cap. Redshifts for galaxies brighter $r = 18.1$ were obtained mainly with LAMOST. In each band, LFs derived using both parametric and non-parametric maximum likelihood methods agree well with each other. In the $^{0.1}r$ band, our fitting parameters of the Schechter function are $\phi_{*}=(1.65\pm0.36)\times10^{-2}h^{3}Mpc^{-3}$, $M_{*}=-20.69\pm0.06$ mag, and $\alpha=-1.12\pm0.08$, in agreements with previous studies. Separate LFs are also derived for emission line galaxies and absorption line galaxies, respectively. The LFs of absorption line galaxies show a dip at $^{0.1}r \sim 18.5$ and can be well fitted by a double-Gaussian function, suggesting a bi-modality in passive galaxies.

Fiber collision is a persistent problem faced by modern spectroscopic galaxy surveys. In this work, we propose a new method to correct for this undesired effect, focusing on the clustering from the fiber-collision scale up to $\lesssim 10 \rm Mpc\,h^{-1}$. We assume that the fiber-collided galaxies are in association with their nearest three angular neighbors. Compared with the conventional nearest-neighbor method, we have properly accounted for the foreground (background) galaxies that are associated with the foreground (background) cosmic webs relative to the nearest neighbor. We have tested the new method with mock catalogs of the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). The test demonstrates that our new method can recover the projected two-point correlation functions at an accuracy better than 1\% on small (below the fiber-collision scale) to intermediate (i.e., $10 \rm Mpc\,h^{-1}$) scales, where the fiber collision takes effect and the SDSS main sample can probe. The new method also gives a better recovery of the redshift-space correlation functions almost on all scales that we are interested in.

We develop a novel approach in exploring the joint dependence of halo bias on multiple halo properties using Gaussian process regression. Using a $\Lambda$CDM $N$-body simulation, we carry out a comprehensive study of the joint bias dependence on halo structure, formation history and environment. We show that the bias is a multivariate function of halo properties that falls into three regimes. For massive haloes, halo mass explains the majority of bias variation. For early-forming haloes, bias depends sensitively on the recent mass accretion history. For low-mass and late-forming haloes, bias depends more on the structure of a halo such as its shape and spin. Our framework enables us to convincingly prove that $V_\mathrm{max}/V_\mathrm{vir}$ is a lossy proxy of formation time for bias modelling, whereas the mass, spin, shape and formation time variables are non-redundant with respect to each other. Combining mass and formation time largely accounts for the mass accretion history dependence of bias. Combining all the internal halo properties fully accounts for the density profile dependence inside haloes, and predicts the clustering variation of individual haloes to a $20\%$ level at $\sim 10\mathrm{Mpc}h^{-1}$. When an environmental density is measured outside $1\mathrm{Mpc}h^{-1}$ from the halo centre, it outperforms and largely accounts for the bias dependence on the internal halo structure, explaining the bias variation above a level of $30\%$.

I present a large set of high resolution simulations, called CosmicGrowth Simulations, which were generated with either 8.6 billion or 29 billion particles. As the nominal cosmological model that can match nearly all observations on cosmological scales, I have adopted a flat Cold Dark Matter (CDM) model with a cosmological constant $\Lambda$ (${\Lambda}$CDM). The model parameters have been taken either from the latest result of the WMAP satellite (WMAP ${\Lambda}$CDM) or from the first year's result of the Planck satellite (Planck ${\Lambda}$CDM). Six simulations are produced in the ${\Lambda}$CDM models with two in the Planck model and the others in the WMAP model. In order for studying the nonlinear evolution of the clustering, four simulations were also produced with $8.6$ billion particles for the scale-free models of an initial power spectrum $P(k)\propto k^n$ with $n=0$,$-1$,$-1.5$ or $-2.0$. Furthermore, two radical CDM models (XCDM) are simulated with 8.6 billion particles each. Since the XCDM have some of the model parameters distinct from those of the ${\Lambda}$CDM models, they must be unable to match the observations, but are very useful for studying how the clustering properties depend on the model parameters. The Friends-of-Friends (FoF) halos were identified for each snapshot and subhalos were produced by the Hierarchical Branch Tracing (HBT) algorithm. These simulations form a powerful database to study the growth and evolution of the cosmic structures both in theories and in observations.

The galaxy intrinsic alignment (IA) is a major challenge of weak lensing cosmology. To alleviate this problem, Zhang (2010, MNRAS, 406, L95) proposed a self-calibration method, independent of IA modeling. This proposal relies on several scaling relations between two-point clustering of IA and matter/galaxy fields, which were previously only tested with analytical IA models. In this paper, these relations are tested comprehensively with an $N$-body simulation of $3072^3$ simulation particles and boxsize 600 $h^{-1} \, \mathrm{Mpc}$. They are verified at the accuracy level of $\mathcal{O}(1)\%$ over angular scales and source redshifts of interest. We further confirm that these scaling relations are generic, insensitive to halo mass, weighting in defining halo ellipticities, photo-$z$ error, and misalignment between galaxy ellipticities and halo ellipticities. We also present and verify three new scaling relations on the B-mode IA. These results consolidate and complete the theory side of the proposed self-calibration technique.

A long-standing issue in peculiar velocity cosmology is whether the halo/galaxy velocity bias $b_v=1$ at large scale. The resolution of this important issue must resort to high precision cosmological simulations. However, this is hampered by another long-standing `sampling artifact' problem in volume weighted velocity measurement. We circumvent this problem with a hybrid approach. We first measure statistics free of sampling artifact, then link them to volume weighted statistics in theory, finally solve for the velocity bias. $b_v$ determined by our method is not only free of sampling artifact, but also free of cosmic variance. We apply this method to a $\Lambda$CDM N-body simulation of $3072^3$ particles and $1200 Mpc/{\rm h}$ box size. For the first time, we determine the halo velocity bias to $0.1\%$-$1\%$ accuracy. Our major findings are as follows: (1) $b_v\neq 1$ at $k>0.1 h/{\rm Mpc}$. The deviation from unity ($|b_v-1|$) increases with $k$. Depending on halo mass and redshift, it may reach $\mathcal{O}(0.01)$ at $k=0.2 h/{\rm Mpc}$ and $\mathcal{O}(0.05)$ at $k\sim 0.3 h/{\rm Mpc}$. The discovered $b_v\neq 1$ has statistically significant impact on structure growth rate measurement by spectroscopic redshift surveys, including DESI, Euclid and SKA. (2) Both the sign and the amplitude of $b_v-1$ depend on mass and redshift. These results disagree with the peak model prediction in that $b_v$ has much weaker deviation from unity, varies with redshift, and can be bigger than unity. (3) Most of the mass and redshift dependences can be compressed into a single dependence on the halo density bias. Based on this finding, we provide an approximate two-parameter fitting formula.

The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.

In a recent study, using the distribution of galaxies in the north galactic pole of SDSS DR7 region enclosed in a 500$\mpch$ box, we carried out our ELUCID simulation (Wang et al. 2016, ELUCID III). Here we {\it light} the dark matter halos and subhalos in the reconstructed region in the simulation with galaxies in the SDSS observations using a novel {\it neighborhood} abundance matching method. Before we make use of thus established galaxy-subhalo connections in the ELUCID simulation to evaluate galaxy formation models, we set out to explore the reliability of such a link. For this purpose, we focus on the following a few aspects of galaxies: (1) the central-subhalo luminosity and mass relations; (2) the satellite fraction of galaxies; (3) the conditional luminosity function (CLF) and conditional stellar mass function (CSMF) of galaxies; and (4) the cross correlation functions between galaxies and the dark matter particles, most of which are measured separately for all, red and blue galaxy populations. We find that our neighborhood abundance matching method accurately reproduces the central-subhalo relations, satellite fraction, the CLFs and CSMFs and the biases of galaxies. These features ensure that thus established galaxy-subhalo connections will be very useful in constraining galaxy formation processes. And we provide some suggestions on the three levels of using the galaxy-subhalo pairs for galaxy formation constraints. The galaxy-subhalo links and the subhalo merger trees in the SDSS DR7 region extracted from our ELUCID simulation are available upon request.

We study the radial acceleration relation (RAR) for early-type galaxies (ETGs) in the SDSS MaNGA MPL5 dataset. The complete ETG sample show a slightly offset RAR from the relation reported by McGaugh et al. (2016) at the low-acceleration end; we find that the deviation is due to the fact that the slow rotators show a systematically higher acceleration relation than the McGaugh's RAR, while the fast rotators show a consistent acceleration relation to McGaugh's RAR. There is a 1\sigma significant difference between the acceleration relations of the fast and slow rotators, suggesting that the acceleration relation correlates with the galactic spins, and that the slow rotators may have a different mass distribution compared with fast rotators and late-type galaxies. We suspect that the acceleration relation deviation of slow rotators may be attributed to more galaxy merger events, which would disrupt the original spins and correlated distributions of baryons and dark matter orbits in galaxies.

We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large scale environment. Matter densities are inferred from our ELUCID simulation, a constrained simulation of local Universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the `environmental quenching efficiency', which quantifies the quenched fraction as function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass-independence of density-based quenching efficiency, found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population and the difference between the two populations found previously mainly arises from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these halo-mass and stellar-mass effects, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites.

We show that the mass of a dark matter halo can be inferred from the dynamical status of its satellite galaxies. Using 9 dark-matter simulations of halos like the Milky Way (MW), we find that the present-day substructures in each halo follow a characteristic distribution in the phase space of orbital binding energy and angular momentum, and that this distribution is similar from halo to halo but has an intrinsic dependence on the halo formation history. We construct this distribution directly from the simulations for a specific halo and extend the result to halos of similar formation history but different masses by scaling. The mass of an observed halo can then be estimated by maximizing the likelihood in comparing the measured kinematic parameters of its satellite galaxies with these distributions. We test the validity and accuracy of this method with mock samples taken from the simulations. Using the positions, radial velocities, and proper motions of 9 tracers and assuming observational uncertainties comparable to those of MW satellite galaxies, we find that the halo mass can be recovered to within $\sim$40%. The accuracy can be improved to within $\sim$25% if 30 tracers are used. However, the dependence of the phase-space distribution on the halo formation history sets a minimum uncertainty of $\sim$20% that cannot be reduced by using more tracers. We believe that this minimum uncertainty also applies to any mass determination for a halo when the phase space information of other kinematic tracers is used.

The ground-breaking detections of gravitational waves from black hole mergers by LIGO have rekindled interest in primordial black holes (PBHs) and the possibility of dark matter being composed of PBHs. It has been suggested that PBHs of tens of solar masses could serve as dark matter candidates. Recent analytical studies demonstrated that compact ultra-faint dwarf galaxies can serve as a sensitive test for the PBH dark matter hypothesis, since stars in such a halo-dominated system would be heated by the more massive PBHs, their present-day distribution can provide strong constraints on PBH mass. In this study, we further explore this scenario with more detailed calculations, using a combination of dynamical simulations and Bayesian inference methods. The joint evolution of stars and PBH dark matter is followed with a Fokker-Planck code PhaseFlow. We run a large suite of such simulations for different dark matter parameters, then use a Markov Chain Monte Carlo approach to constrain the PBH properties with observations of ultra-faint galaxies. We find that two-body relaxation between the stars and PBH drives up the stellar core size, and increases the central stellar velocity dispersion. Using the observed half-light radius and velocity dispersion of stars in the compact ultra-faint dwarf galaxies as joint constraints, we infer that these dwarfs may have a cored dark matter halo with the central density in the range of 1-2 $\rm{M_{\odot}/pc^3}$, and that the PBHs may have a mass range of 2-14 $\rm{M_{\odot}}$ if they constitute all or a substantial fraction of the dark matter.

We present a spectroscopic redshift catalog from the LAMOST Complete Spectroscopic Survey of Pointing Area (LaCoSSPAr) in the Southern Galactic Cap (SGC), which is designed to observe all sources (Galactic and extra-galactic) by using repeating observations with a limiting magnitude of $r=18.1~mag$ in two $20~deg^2$ fields. The project is mainly focusing on the completeness of LAMOST ExtraGAlactic Surveys (LEGAS) in the SGC, the deficiencies of source selection methods and the basic performance parameters of LAMOST telescope. In both fields, more than 95% of galaxies have been observed. A post-processing has been applied to LAMOST 1D spectrum to remove the majority of remaining sky background residuals. More than 10,000 spectra have been visually inspected to measure the redshift by using combinations of different emission/absorption features with uncertainty of $\sigma_{z}/(1+z)<0.001$. In total, there are 1528 redshifts (623 absorption and 905 emission line galaxies) in Field A and 1570 redshifts (569 absorption and 1001 emission line galaxies) in Field B have been measured. The results show that it is possible to derive redshift from low SNR galaxies with our post-processing and visual inspection. Our analysis also indicates that up to 1/4 of the input targets for a typical extra-galactic spectroscopic survey might be unreliable. The multi-wavelength data analysis shows that the majority of mid-infrared-detected absorption (91.3%) and emission line galaxies (93.3%) can be well separated by an empirical criterion of $W2-W3=2.4$. Meanwhile, a fainter sequence paralleled to the main population of galaxies has been witnessed both in $M_r$/$W2-W3$ and $M_*$/$W2-W3$ diagrams, which could be the population of luminous dwarf galaxies but contaminated by the edge-on/highly inclined galaxies ($\sim30\%$).

In this paper we use high-resolution cosmological simulations to study halo intrinsic alignment and its dependence on mass, formation time and large-scale environment. In agreement with previous studies using N-body simulations, it is found that massive halos have stronger alignment. For given mass, older halos have stronger alignment than younger ones. By identifying the cosmic environment of halo using Hessian matrix, we find that for given mass, halos in cluster regions also have stronger alignment than those in filament. The existing theory has not addressed these dependencies explicitly. In this work we extend the linear alignment model with inclusion of halo bias and find that the halo alignment with its mass and formation time dependence can be explained by halo bias. However, the model can not account for the environment dependence, as it is found that halo bias is lower in cluster and higher in filament. Our results suggest that halo bias and environment are independent factors in determining halo alignment. We also study the halo alignment correlation function and find that halos are strongly clustered along their major axes and less clustered along the minor axes. The correlated halo alignment can extend to scale as large as $100h^{-1}$Mpc where its feature is mainly driven by the baryon acoustic oscillation effect.

We study the dependence of galaxy clustering on atomic gas mass using a sample of $\sim$16,000 galaxies with redshift in the range of $0.0025<z<0.05$ and HI mass of $M_{\rm HI}>10^8M_{\odot}$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with $M_{\rm HI}$ above different thresholds, and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample selected from the Sloan Digital Sky Survey, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak HI-mass dependence, we find both the clustering amplitude on scales above a few Mpc and the bias factors to increase significantly with increasing HI mass for subsamples with HI mass thresholds above $10^9M_{\odot}$. For HI mass thresholds below $10^9M_{\odot}$, while the measurements have large uncertainties caused by the limited survey volume and sample size, the inferred galaxy bias factors are systematically lower than the minimum halo bias factor from mass-selected halo samples. The simple halo model, in which galaxy content is only determined by halo mass, has difficulties in interpreting the clustering measurements of the HI-selected samples. We extend the simple model by including the halo formation time as an additional parameter. A model that puts HI-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of HI mass with halo mass and formation time, as well as the halo occupation distributions and HI mass functions for central and satellite galaxies. These results are compared with the predictions from semi-analytic galaxy formation models and hydrodynamic galaxy formation simulations.

We take advantage of the statistical power of the large-volume dark-matter-only Millennium simulation, combined with a sophisticated semi-analytic galaxy formation model, to explore whether the recently reported $z=3.7$ quiescent galaxy ZF-COSMOS-20115 (ZF; Glazebrook et al. 2017) can be accommodated in current galaxy formation models. In our model, a population of quiescent galaxies (QGs) with stellar masses and star formation rates comparable to those of ZF naturally emerges at redshifts $z<4$. There are two and five ZF analogues at the redshift $3.86$ and $3.58$ in the Millennium simulation volume, respectively. We demonstrate that, while the $z>3.5$ massive QGs are rare (about 2\% of the galaxies with the similar stellar masses), the existing AGN feedback model implemented in the semi-analytic galaxy formation model can successfully explain the formation of the high-redshift QGs as it does on their lower redshift counterparts.

Massive spectroscopic redshift surveys open a promising window to accurately measure peculiar velocity at cosmological distances through redshift space distortion (RSD). In paper I of this series of work we proposed to decompose peculiar velocity into three eigen-modes (v_\delta, v_S and v_B) in order to facilitate the RSD modeling and peculiar velocity reconstruction. In the current paper we measure the dark matter RSD related statistics of the velocity eigen-modes through a set of N-body simulations, including the velocity power spectra, correlation functions, one-point probability distribution functions, cumulants and the damping functions describing the Finger of God effect. (1) The power spectrum measurement shows that these velocity components have distinctly different spatial distribution and redshift evolution. In particular, we measure the window function \tilde{W}(k,z), which describes the impact of nonlinear evolution on the v_\delta-density relation. We confirm that it can induce a significant systematic error of O(10%) in RSD cosmology. We demonstrate that \tilde{W} can be accurately described by a simple fitting formula with one or two free parameters. (2) The correlation function measurement shows that the correlation length is O(100), O(10) and O(1) Mpc for v_\delta, v_S and v_B respectively. These correlation lengths determine where we can treat the velocity fields as spatially uncorrelated. (3) The velocity PDFs and cumulants quantify non-Gaussianities of the velocity fields. We confirm speculation in paper I that v_\delta is largely Gaussian, nevertheless with non-negligible non-Gaussianity, v_B is significantly non-Gaussian. We also measure the damping functions. Despite the observed non-Gaussianities, the damping functions and hence the FOG effect are all well approximated as Gaussian ones at scales of interest.

The ELUCID project aims to build a series of realistic cosmological simulations that reproduce the spatial and mass distribution of the galaxies as observed in the Sloan Digital Sky Survey (SDSS). This requires powerful reconstruction techniques to create constrained initial conditions. We test the reconstruction method by applying it to several $N$-body simulations. We use 2 medium resolution simulations from each of which three additional constrained $N$-body simulations were produced. We compare the resulting friend of friend catalogs by using the particle indexes as tracers, and quantify the quality of the reconstruction by varying the main smoothing parameter. The cross identification method we use proves to be efficient, and the results suggest that the most massive reconstructed halos are effectively traced from the same Lagrangian regions in the initial conditions. Preliminary time dependence analysis indicates that high mass end halos converge only at a redshift close to the reconstruction redshift. This suggests that, for earlier snapshots, only collections of progenitors may be effectively cross-identified.

Measuring the volume-weighted peculiar velocity statistics from inhomogeneously and sparsely distributed galaxies/halos, by existing velocity assignment methods, suffers from a significant sampling artifact. As an alternative, the Kriging interpolation based on Gaussian processes was introduced and evaluated [Y. Yu, J. Zhang, Y. Jing, and P. Zhang, Phys. Rev. D 92, 083527 (2015)]. Unfortunately, the most straightforward application of Kriging does not perform better than the existing methods in the literature. In this work, we investigate two physically motivated extensions. The first takes into account of the anisotropic velocity correlations. The second introduces the nugget effect, on account of multi-streaming of the velocity field. We find that the performance is indeed improved. For sparsely sampled data [$n_P\lesssim 6\times 10^{-3}(h^{-1} {\rm Mpc})^{-3}$] where the sampling artifact is the most severe, the improvement is significant and is two-fold: 1) The scale of reliable measurement of the velocity power spectrum is extended by a factor $\sim 1.6$, and 2) the dependence on the velocity correlation prior is weakened by a factor of $\sim 2$. We conclude that such extensions are desirable for accurate velocity assignment by Kriging.

We study the integral Galactic extinction and reddening based on the galaxy catalog of the South Galactic Cap U-band Sky Survey (SCUSS), where $u$ band galaxy number counts and $u-r$ color distribution are used to derive the Galactic extinction and reddening respectively. We compare these independent statistical measurements with the reddening map of \citet{Schlegel1998}(SFD) and find that both the extinction and reddening from the number counts and color distribution are in good agreement with the SFD results at low extinction regions ($E(B-V)^{SFD}<0.12$ mag). However, for high extinction regions ($E(B-V)^{SFD}>0.12$ mag), the SFD map overestimates the Galactic reddening systematically, which can be approximated by a linear relation $\Delta E(B-V)= 0.43[E(B-V)^{SFD}-0.12$]. By combing the results of galaxy number counts and color distribution together, we find that the shape of the Galactic extinction curve is in good agreement with the standard $R_V=3.1$ extinction law of \cite{ODonnell1994}.

We present and analyze the possibility of using optical ${\it u}$-band luminosities to estimate star-formation rates (SFRs) of galaxies based on the data from the South Galactic Cap ${\it u }$ band Sky Survey (SCUSS), which provides a deep ${\it u}$-band photometric survey covering about 5000 $deg^2$ of the South Galactic Cap. Based on two samples of normal star-forming galaxies selected by the BPT diagram, we explore the correlations between ${\it u}$-band, H$\alpha$, and IR luminosities by combing SCUSS data with the Sloan Digital Sky Survey (SDSS) and ${\it Wide}$-${\it field\ Infrared\ Survey\ Explorer}$ (${\it WISE}$). The attenuation-corrected ${\it u}$-band luminosities are tightly correlated with the Balmer decrement-corrected H$\alpha$ luminosities with an rms scatter of $\sim$ 0.17 dex. The IR-corrected ${\it u }$ luminosities are derived based on the correlations between the attenuation of ${\it u}$-band luminosities and ${\it WISE}$ 12 (or 22) $\mu$m luminosities, and then calibrated with the Balmer-corrected H$\alpha$ luminosities. The systematic residuals of these calibrations are tested against the physical properties over the ranges covered by our sample objects. We find that the best-fitting nonlinear relations are better than the linear ones and recommended to be applied in the measurement of SFRs. The systematic deviations mainly come from the pollution of old stellar population and the effect of dust extinction; therefore, a more detailed analysis is needed in the future work.

To take full advantage of the unprecedented power of upcoming weak lensing surveys, understanding the noise, such as cosmic variance and geometry/mask effects, is as important as understanding the signal itself. Accurately quantifying the noise requires a large number of statistically independent mocks for a variety of cosmologies. This is impractical for weak lensing simulations, which are costly for simultaneous requirements of large box size (to cover a significant fraction of the past light cone) and high resolution (to robustly probe the small scale where most lensing signal resides). Therefore fast mock generation methods are desired and are under intensive investigation. We propose a new fast weak lensing map generation method, named the inverse-Gaussianization method, based on the finding that a lensing convergence field can be Gaussianized to excellent accuracy by a local transformation [Yu et al, Phys. Rev. D 84, 023523 (2011)]. Given a simulation, it enables us to produce as many as infinite statistically independent lensing maps as fast as producing the simulation initial conditions. The proposed method is tested against simulations for each tomography bin centered at lens redshift $z \sim 0.5$, 1, and 2, with various statistics. We find that the lensing maps generated by our method have reasonably accurate power spectra, bispectra, and power spectrum covariance matrix. Therefore, it will be useful for weak lensing surveys to generate realistic mocks. As an example of application, we measure the probability distribution function of the lensing power spectrum, from 16384 lensing maps produced by the inverse-Gaussianization method.

We present a measurement of two-dimensional (2D) redshift-space power spectrum for the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 CMASS galaxies in the North Galactic Cap (NGC) based on the method developed by Jing & Borner (2001). In this method, we first measure the 2D redshift-space correlation function for the CMASS galaxies, and obtain the 2D power spectrum based on Fourier Transform of the correlation function. The method is tested with an N-body mock galaxy catalog, which demonstrates that the method can yield an accurate and unbiased measurement of the redshift-space power spectrum given the input 2D correlation function is correct. Compared with previous measurements in literature that are usually based on direct Fourier Transform in redshift space, our method has the advantages that the window function and shot-noise are fully corrected. In fact, our 2D power spectrum, by its construction, can accurately reproduce the 2D correlation function, and in the meanwhile can reproduce, for example, the 2D power spectrum of Beutler et al. (2014) accurately if ours is convolved with the window function they provided. Thus, our measurement can facilitate a direct comparison with the theoretical predictions. With this accurate measurement of the 2D power spectrum, we then develop a method to measure the structure growth rate, by separating the anisotropic redshift-space power spectrum from the isotropic real-space power spectrum. We have also carefully corrected for the nonlinearities in the mapping from real space to redshift space, according to the theoretical model of Zhang et al. (2013). Finally, we obtain f(zeff)sigma_8(zeff)=0.438\pm0.037 at the effective redshift zeff=0.57, where f(zeff) is the linear growth rate at redshift zeff. The result is useful for constraining cosmological parameters. The measurements of 2D power spectrum will be released soon.

Based on the LAMOST survey and Sloan Digital Sky Survey (SDSS), we use low-resolution spectra of 130,043 F/G-type dwarf stars to study the kinematics and metallicity properties of the Galactic disk. Our study shows that the stars with poorer metallicity and larger vertical distance from Galactic plane tend to have larger eccentricity and velocity dispersion. After separating the sample stars into likely thin-disk and thick-disk sub-sample, we find that there exits a negative gradient of rotation velocity $V_{\phi}$ with metallicity [Fe/H] for the likely thin-disk sub-sample, and the thick-disk sub-sample exhibit a larger positive gradient of rotation velocity with metallicity. By comparing with model prediction, we consider the radial migration of stars appears to have influenced on the thin-disk formation. In addition, our results shows that the observed thick-disk stellar orbital eccentricity distribution peaks at low eccentricity ($e \sim 0.2$) and extends to a high eccentricity ($e \sim 0.8$). We compare this result with four thick-disk formation simulated models, and it imply that our result is consistent with gas-rich merger model.

A method we developed recently for the reconstruction of the initial density field in the nearby Universe is applied to the Sloan Digital Sky Survey Data Release 7. A high-resolution N-body constrained simulation (CS) of the reconstructed initial condition, with $3072^3$ particles evolved in a 500 Mpc/h box, is carried out and analyzed in terms of the statistical properties of the final density field and its relation with the distribution of SDSS galaxies. We find that the statistical properties of the cosmic web and the halo populations are accurately reproduced in the CS. The galaxy density field is strongly correlated with the CS density field, with a bias that depend on both galaxy luminosity and color. Our further investigations show that the CS provides robust quantities describing the environments within which the observed galaxies and galaxy systems reside. Cosmic variance is greatly reduced in the CS so that the statistical uncertainties can be controlled effectively even for samples of small volumes.

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

We develop a new Monte-Carlo-based method to convert the SDSS (Sloan Digital Sky Survey) $u$-band magnitude to the SCUSS (South Galactic Cap of $u$-band Sky Survey) $u$-band magnitude. Due to more accuracy of SCUSS $u$-band measurements, the converted $u$-band magnitude becomes more accurate comparing with the original SDSS $u$-band magnitude, in particular at the faint end. The average $u$ (both SDSS and SCUSS) magnitude error of numerous main-sequence stars with $0.2<g-r<0.8$ increase as $g$-band magnitude becomes fainter. When $g=19.5$, the average magnitude error of SDSS $u$ is 0.11. When $g=20.5$, the average SDSS $u$ error is up to 0.22. However, at this magnitude, the average magnitude error of SCUSS $u$ is just half as much as that of SDSS $u$. The SDSS $u$-band magnitudes of main-sequence stars with $0.2<g-r<0.8$ and $18.5<g<20.5$ are converted, therefore the maximum average error of converted $u$-band magnitudes is 0.11. The potential application of this conversion is to derive more accurate photometric metallicity calibration from SDSS observation, especially for those distant stars. Thus, we can explore stellar metallicity distributions either in the Galactic halo or some stream stars.

Based on the Sloan Digital Sky Survey (SDSS), we develop a new monte-carlo based method to estimate the photometric metallicity distribution function (MDF) for stars in the Milky Way. Compared with other photometric calibration methods, this method enables a more reliable determination of the MDF, in particular at the metal-poor and metal-rich ends. We present a comparison of our new method with a previous polynomial-based approach, and demonstrate its superiority. As an example, we apply this method to main-sequence stars with $0.2<g-r<0.6$, $6$ kpc$<R<9$ kpc, and in different intervals in height above the plane, $|Z|$. The MDFs for the selected stars within two relatively local intervals ($0.8$ kpc$<|Z|<1.2$ kpc, $1.5$ kpc$<|Z|<2.5$ kpc) can be well-fit by two Gaussians, with peaks at [Fe/H] $\approx-0.6$ and $-1.2$ respectively, one associated with the disk system, the other with the halo. The MDFs for the selected stars within two more distant intervals ($3$ kpc$<|Z|<5$ kpc, $6$ kpc$<|Z|<9$ kpc) can be decomposed into three Gaussians, with peaks at [Fe/H] $\approx-0.6$, $-1.4$ and $-1.9$ respectively, where the two lower peaks may provide evidence for a two-component model of the halo: the inner halo and the outer halo. The number ratio between the disk component and halo component(s) decreases with vertical distance from the Galactic plane, consistent with the previous literature.

Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.

We investigate the clustering and dark matter halo mass for a sample of $\sim$16,000 central galaxies selected from the SDSS/DR7 group catalog. We select subsamples of central galaxies on three two-dimensional planes, each formed by stellar mass ($M_{\star}$) and one of the other properties including optical color (g-r) , surface stellar mass density ($\mu_{\star}$) and central stellar velocity dispersion ($\sigma_{\star}$). For each subsample we measure both the projected cross-correlation function ($w_{p}(r_{p})$) relative to a reference galaxy sample, and an average mass of the host dark matter halos ($M_{halo}$) . For comparison we have also estimated the $w_{p}(r_{p})$ for the full galaxy population and the subset of satellite galaxies. We find that, for central galaxies, both $w_{p}(r_{p})$ and $M_{halo}$ show strongest dependence on $M_{\star}$, and there is no clear dependence on other properties when stellar mass is fixed. This result provides strong support to the previously-adopted assumption that, for central galaxies, stellar mass is the galaxy property that is best indicative of the host dark halo mass. The full galaxy population and the subset of satellites show similar clustering properties in all cases. However, they are similar to the centrals only at high masses ($M_{\star} \geq 10^{11} M_{\odot}$). At low-mass ($M_{\star} \leq 10^{11} M_{\odot}$), the results are different: the $w_{p}(r_{p})$ increases with both $\sigma$ and g-r when $M_{\star}$ is fixed, and depends very weakly on $M_{\star}$ when $\sigma_{\star}$ or g-r is fixed. At fixed $M_{star}$, the $\mu_{\star}$ shows weak correlations with clustering amplitude (and halo mass for central galaxies). Our results suggest that it is necessary to consider central and satellite galaxies separately when studying the link between galaxies and dark matter halos. (Abridged)

The SCUSS is a deep $u$-band imaging survey in the south Galactic cap using the 2.3m Bok telescope. The survey observations were completed in the end of 2013, covering an area of about 5000 square degrees. We release the data in the region with an area of about 4000 deg$^2$ that is mostly covered by the Sloan digital sky survey. The data products contain calibrated single-epoch images, stacked images, photometric catalogs, and a catalog of star proper motions derived by Peng et al, 2015. The median seeing and magnitude limit ($5\sigma$) are about 2".0 and 23.2 mag, respectively. There are about 8 million objects having measurements of absolute proper motions. All the data and related documentations can be accessed through the SCUSS data release website of \url{http://batc.bao.ac.cn/Uband/data.html}.

We present a new selection technique of producing spectroscopic target catalogues for massive spectroscopic surveys for cosmology. This work was conducted in the context of the extended Baryon Oscillation Spectroscopic Survey (eBOSS), which will use ~200 000 emission line galaxies (ELGs) at 0.6<zspec<1.0 to obtain a precise baryon acoustic oscillation measurement. Our proposed selection technique is based on optical and near-infrared broad-band filter photometry. We used a training sample to define a quantity, the Fisher discriminant (linear combination of colours), which correlates best with the desired properties of the target: redshift and [OII] flux. The proposed selections are simply done by applying a cut on magnitudes and this Fisher discriminant. We used public data and dedicated SDSS spectroscopy to quantify the redshift distribution and [OII] flux of our ELG target selections. We demonstrate that two of our selections fulfil the initial eBOSS/ELG redshift requirements: for a target density of 180 deg^2, ~70% of the selected objects have 0.6<zspec<1.0 and only ~1% of those galaxies in the range 0.6<zspec<1.0 are expected to have a catastrophic zspec estimate. Additionally, the stacked spectra and stacked deep images for those two selections show characteristic features of star-forming galaxies. The proposed approach using the Fisher discriminant could, however, be used to efficiently select other galaxy populations, based on multi-band photometry, providing that spectroscopic information is available. This technique could thus be useful for other future massive spectroscopic surveys such as PFS, DESI, and 4MOST.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

Fornax is the brightest Milky Way (MW) dwarf spheroidal galaxy and its star formation history (SFH) has been derived from observations. We estimate the time evolution of its gas mass and net inflow and outflow rates from the SFH using a simple star formation law that relates the star formation rate to the gas mass. We present a chemical evolution model on a 2D mass grid with supernovae (SNe) as sources of metal enrichment. We find that a key parameter controlling the enrichment is the mass M_x of the gas to mix with the ejecta from each SN. The choice of M_x depends on the evolution of SN remnants and on the global gas dynamics. It differs between the two types of SNe involved and between the periods before and after Fornax became an MW satellite at time t = t_sat . Our results indicate that due to the global gas outflow at t > t_sat , part of the ejecta from each SN may directly escape from Fornax. Sample results from our model are presented and compared with data.

N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halo-centric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped from individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the gravitational lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so called "anti-bias", as a simple selection effect. We provide a Python code SubGen for populating haloes with subhaloes at http://icc.dur.ac.uk/data/

A small fraction($<10\%$) of SDSS main sample galaxies(MGs) have not been targeted with spectroscopy due to the the fiber collision effect. These galaxies have been compiled into the input catalog of the LAMOST extra-galactic survey and named as the complementary galaxy sample. In this paper, we introduce the project and the status of the spectroscopies of the complementary galaxies in the first two years of the LAMOST spectral survey(till Sep. of 2014). Moreover, we present a sample of 1,102 galaxy pairs identified from the LAMOST complementary galaxies and SDSS MGs, which are defined as that the two members have a projected distance smaller than 100 kpc and the recessional velocity difference smaller than 500 $\rm kms^{-1}$. Compared with the SDSS only selected galaxy pairs, the LAMOST-SDSS pairs take the advantages of not being biased toward large separations and therefor play as a useful supplement to the statistical studies of galaxy interaction and galaxy merging.

We propose that detailed data on the star formation history of a dwarf spheroidal galaxy (dSph) may be used to estimate the evolution of the total mass M_g(t) for cold gas in its star-forming disk. Using Fornax as an illustrative example, we estimate its M_g(t) and the corresponding net gas flow rate Delta F(t) assuming a global star formation rate psi(t) = lambda_*(t)[M_g(t)/M_sun]^alpha consistent with observations of nearby galaxies. We identify the onset of the transition in Delta F(t) from a net inflow to a net outflow as the time t_sat at which the Fornax halo became a Milky Way satellite and estimate the evolution of its total mass M_h(t) at t < t_sat using the median halo growth history in the LambdaCDM cosmology and its present mass within the half-light radius derived from observations. We examine three different cases of alpha = 1, 1.5, and 2, and justify the corresponding lambda_*(t) by comparing the gas mass fraction f_g(t) = M_g(t)/M_h(t) at t < t_sat with results from simulations of gas accretion by halos in a reionized universe. We find that the Fornax halo grew to M_h(t_sat) ~ 2x10^9 M_sun at t_sat ~ 5 or 8 Gyr, in broad agreement with previous studies using data on its stellar kinematics and its orbital motion. We describe qualitatively the evolution of Fornax as a satellite and discuss potential extension of our approach to other dSphs.

[abridged] Volume-weighted statistics of large scale peculiar velocity is preferred by peculiar velocity cosmology, since it is free of uncertainties of galaxy density bias entangled in mass-weighted statistics. However, measuring the volume-weighted velocity statistics from galaxy (halo/simulation particle) velocity data is challenging. For the first time, we apply the Kriging interpolation to obtain the volume-weighted velocity field. Kriging is a minimum variance estimator. It predicts the most likely velocity for each place based on the velocity at other places. We test the performance of Kriging quantified by the E-mode velocity power spectrum from simulations. Dependences on the variogram prior used in Kriging, the number $n_k$ of the nearby particles to interpolate and the density $n_P$ of the observed sample are investigated. First, we find that Kriging induces $1\%$ and $3\%$ systematics at $k\sim 0.1h{\rm Mpc}^{-1}$ when $n_P\sim 6\times 10^{-2} ({\rm Mpc}/h)^{-3}$ and $n_P\sim 6\times 10^{-3} ({\rm Mpc}/h)^{-3}$, respectively. The deviation increases for decreasing $n_P$ and increasing $k$. When $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$, a smoothing effect dominates small scales, causing significant underestimation of the velocity power spectrum. Second, increasing $n_k$ helps to recover small scale power. However, for $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$ cases, the recovery is limited. Finally, Kriging is more sensitive to the variogram prior for lower sample density. The most straightforward application of Kriging on the cosmic velocity field does not show obvious advantages over the nearest-particle method (Zheng et al. 2013) and could not be directly applied to cosmology so far. However, whether potential improvements may be achieved by more delicate versions of Kriging is worth further investigation.

The South Galactic Cap u-band Sky Survey (SCUSS) is a deep u-band imaging survey in the Southern Galactic Cap, using the 90Prime wide-field imager on the 2.3m Bok telescope at Kitt Peak. The survey observations started in 2010 and ended in 2013. The final survey area is about 5000 deg2 with a median 5-sigma point source limiting magnitude of about 23.2. This paper describes the survey data reduction process, which includes basic imaging processing, astrometric and photometric calibrations, image stacking, and photometric measurements. Survey photometry is performed on objects detected both on SCUSS u-band images and in the SDSS database. Automatic, aperture, point-spread function (PSF), and model magnitudes are measured on stacked images. Co-added aperture, PSF, and model magnitudes are derived from measurements on single-epoch images. We also present comparisons of the SCUSS photometric catalog with those of the SDSS and CFHTLS.

SCUSS is a u-band photometric survey covering about 4000 square degree of the South Galactic Cap, reaching depths of up to 23 mag. By extending around 1.5 mag deeper than SDSS single-epoch u data, SCUSS is able to probe much a larger volume of the outer halo, i.e. with SCUSS data blue horizontal branch (BHB) stars can trace the outer halo of the Milky Way as far as 100-150 kpc. Utilizing this advantage we combine SCUSS u band with SDSS DR9 gri photometric bands to identify BHB stars and explore halo substructures. We confirm the existence of the Pisces overdensity, which is a structure in the outer halo (at around 80 kpc) that was discovered using RR Lyrae stars. For the first time we are able to determine its spatial extent, finding that it appears to be part of a stream with a clear distance gradient. The stream, which is ~5 degrees wide and stretches along ~25 degrees, consists of 20-30 BHBs with a total significance of around 6sigma over the background. Assuming we have detected the entire stream and that the progenitor has fully disrupted, then the number of BHBs suggests the original system was similar to smaller classical or a larger ultra-faint dwarf galaxy. On the other hand, if the progenitor still exists, it can be hunted for by reconstructing its orbit from the distance gradient of the stream. This new picture of the Pisces overdensity sheds new light on the origin of this intriguing system.

Based on SDSS g, r and SCUSS (South Galactic Cap of u-band Sky Survey) $u$ photometry, we develop a photometric calibration for estimating the stellar metallicity from $u-g$ and $g-r$ colors by using the SDSS spectra of 32,542 F- and G-type main sequence stars, which cover almost $3700$ deg$^{2}$ in the south Galactic cap. The rms scatter of the photometric metallicity residuals relative to spectrum-based metallicity is $0.14$ dex when $g-r<0.4$, and $0.16$ dex when $g-r>0.4$. Due to the deeper and more accurate magnitude of SCUSS $u$ band, the estimate can be used up to the faint magnitude of $g=21$. This application range of photometric metallicity calibration is wide enough so that it can be used to study metallicity distribution of distant stars. In this study, we select the Sagittarius (Sgr) stream and its neighboring field halo stars in south Galactic cap to study their metallicity distribution. We find that the Sgr stream at the cylindrical Galactocentric coordinate of $R\sim 19$ kpc, $\left| z\right| \sim 14$ kpc exhibits a relative rich metallicity distribution, and the neighboring field halo stars in our studied fields can be modeled by two-Gaussian model, with peaks respectively at [Fe/H]$=-1.9$ and [Fe/H]$=-1.5$.

Based on SDSS g, r and SCUSS (South Galactic Cap of u-band Sky Survey) $u$ photometry, we develop a photometric calibration for estimating the stellar metallicity from $u-g$ and $g-r$ colors by using the SDSS spectra of 32,542 F- and G-type main sequence stars, which cover almost $3700$ deg$^{2}$ in the south Galactic cap. The rms scatter of the photometric metallicity residuals relative to spectrum-based metallicity is $0.14$ dex when $g-r<0.4$, and $0.16$ dex when $g-r>0.4$. Due to the deeper and more accurate magnitude of SCUSS $u$ band, the estimate can be used up to the faint magnitude of $g=21$. This application range of photometric metallicity calibration is wide enough so that it can be used to study metallicity distribution of distant stars. In this study, we select the Sagittarius (Sgr) stream and its neighboring field halo stars in south Galactic cap to study their metallicity distribution. We find that the Sgr stream at the cylindrical Galactocentric coordinate of $R\sim 19$ kpc, $\left| z\right| \sim 14$ kpc exhibits a relative rich metallicity distribution, and the neighboring field halo stars in our studied fields can be modeled by two-Gaussian model, with peaks respectively at [Fe/H]$=-1.9$ and [Fe/H]$=-1.5$.

A profound assumption in peculiar velocity cosmology is $b_v=1$ at sufficiently large scales, where $b_v$ is the volume weighted halo(galaxy) velocity bias with respect to the matter velocity field. However, this fundamental assumption has not been robustly verified in numerical simulations. Furthermore, it is challenged by structure formation theory (BBKS, 1986, ApJ; Desjacques and Sheth, 2010, PRD), which predicts the existence of velocity bias (at least for proto-halos) due to the fact that halos reside in special regions (local density peaks). The major obstacle to measure the volume weighted velocity from N-body simulations is an unphysical sampling artifact. It is entangled in the measured velocity statistics and becomes significant for sparse populations. With recently improved understanding of the sampling artifact (Zhang, Zheng and Jing, 2015, PRD; Zheng, Zhang and Jing, 2015, PRD), for the first time we are able to {\it appropriately correct this sampling artifact and then robustly measure the volume weighted halo velocity bias}. (1) We verify $b_v=1$ within $2\%$ model uncertainty at $k\lesssim 0.1h/$Mpc and $z=0$-$2$ for halos of mass $\sim 10^{12}$-$10^{13} h^{-1} M_\odot$, and, therefore, consolidates a foundation of the peculiar velocity cosmology. (2) We also find statistically significant signs of $b_v\neq 1$ at $k\gtrsim 0.1h/$Mpc. Unfortunately, whether this is real or caused by residual sampling artifact requires further investigation. Nevertheless, cosmology based on $k\gtrsim 0.1h/$Mpc velocity data shall be careful this potential velocity bias.

The Large sky Area Multi-Object Spectroscopic Telescope (LAMOST) General Survey is a spectroscopic survey that will eventually cover approximately half of the celestial sphere and collect 10 million spectra of stars, galaxies and QSOs. Objects both in the pilot survey and the first year general survey are included in the LAMOST First Data Release (DR1). The pilot survey started in October 2011 and ended in June 2012, and the data have been released to the public as the LAMOST Pilot Data Release in August 2012. The general survey started in September 2012, and completed its first year of operation in June 2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336 spectra, of which 1,790,879 spectra have observed signal-to-noise S/N >10. All data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data policy. This data release contains a total of 2,204,696 spectra, of which 1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars. The DR1 includes not only spectra, but also three stellar catalogues with measured parameters: AFGK-type stars with high quality spectra (1,061,918 entries), A-type stars (100,073 entries), and M stars (121,522 entries). This paper introduces the survey design, the observational and instrumental limitations, data reduction and analysis, and some caveats. Description of the FITS structure of spectral files and parameter catalogues is also provided.

Cosmology based on large scale peculiar velocity preferes volume weighted velocity statistics. However, measuring the volume weighted velocity statistics from inhomogeneously distributed galaxies (simulation particles/halos) suffer from an inevitable and significant sampling artifact. We study this sampling artifact in the velocity power spectrum measured by the nearest-particle (NP) velocity assignment method(Zheng et al. 2013, PRD). We derive the analytical expression of leading and higher order terms. We find that the sampling artifact suppresses the $z=0$ E-mode velocity power spectrum by $\sim 10\%$ at $k=0.1h/$Mpc , for samples with number density $10^{-3}({\rm Mpc}/h)^{-3}$. This suppression becomes larger for larger $k$ and for sparser samples. We argue that, this source of systematic errors in peculiar velocity cosmology, albeit severe, can be self-calibrated in the framework of our theoretical modelling. We also work out the sampling artifact in the density-velocity cross power spectrum measurement. More robust evaluation of related statistics through simulations will be presented in a companion paper (Zheng, Zhang & Jing, 2015, PRD). We also argue that similar sampling artifact exists in other velocity assignment methods and hence must be carefully corrected to avoid systematic bias in peculiar velocity cosmology.

Measuring the volume weighted velocity power spectrum suffers from a severe systematic error, due to imperfect sampling of the velocity field from inhomogeneous distribution of dark matter particles/halos in simulations or galaxies with velocity measurement. This "sampling artifact" depends on both the mean particle number density $\bar{n}_P$ and the intrinsic large scale structure (LSS) fluctuation in the particle distribution. (1) We report robust detection of this sampling artifact in N-body simulations. It causes $\sim 12$% underestimation of the velocity power spectrum at $k=0.1$h/Mpc for samples with $\bar{n}_P=6\times10^{-3}$ (Mpc/h)$^{-3}$. This systematic underestimation increases with decreasing $\bar{n}_P$ and increasing $k$. Its dependence on the intrinsic LSS fluctuations is also robustly detected. (2) All these findings are expected by our theoretical modelling in paper I \cite{Zhang14}. In particular, the leading order theoretical approximation agrees quantitatively well with simulation result for $\bar{n}_P\gtrsim6\times 10^{-4}$(Mpc/h)$^{-3}$. Furthermore, we provide an ansatz to take high order terms into account. It improves the model accuracy to $\lesssim1$% at $k\lesssim0.1$h/Mpc over 3 orders of magnitude in $\bar{n}_P$ and over typical LSS clustering from $z=0$ to $z=2$. (3) The sampling artifact is determined by the deflection ${\bf D}$ field, which is straightforwardly available in both simulations and data of galaxy velocity. Hence the sampling artifact in the velocity power spectrum measurement can be self-calibrated within our framework. By applying such self-calibration in simulations, it becomes promising to determine the {\it real} large scale velocity bias of $10^{13}M_\odot$ halos with $\sim 1$% accuracy, and that of lower mass halos by better accuracy. ...[abridged]

We present the measurements of the redshift-space three-point correlation function (3PCF) for z~0.5 luminous red galaxies of the CMASS sample in the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Data Release 11. The 3PCF measurements are interpreted within the halo occupation distribution framework using high-resolution N-body simulations, and the model successfully reproduces the 3PCF on scales larger than 1Mpc/h. As with the case for the redshift-space two-point correlation functions, we find that the redshift-space 3PCF measurements also favour the inclusion of galaxy velocity bias in the model. In particular, the central galaxy in a halo is on average in motion with respect to the core of the halo. We discuss the potential of the small-scale 3PCF to tighten the constraints on the relation between galaxies and dark matter haloes and on the phase-space distribution of galaxies.

Absolute proper motions for $\sim$ 7.7 million objects were derived based on data from the South Galactic Cap u-band Sky Survey (SCUSS) and astrometric data derived from uncompressed Digitized Sky Surveys that the Space Telescope Science Institute (STScI) created from the Palomar and UK Schmidt survey plates. We put a great deal of effort into correcting the position-, magnitude-, and color-dependent systematic errors in the derived absolute proper motions. The spectroscopically confirmed quasars were used to test the internal systematic and random error of the proper motions. The systematic errors of the overall proper motions in the SCUSS catalog are estimated as -0.08 and -0.06 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}, respectively. The random errors of the proper motions in the SCUSS catalog are estimated independently as 4.2 and 4.4 mas/yr for {\mu}{\alpha} cos {\delta} and {\mu}{\delta}. There are no obvious position-, magnitude-, and color-dependent systematic errors of the SCUSS proper motions. The random error of the proper motions goes up with the magnitude from about 3 mas/yr at u < 18.0 mag to about 7 mas/yr at u = 22.0 mag. The proper motions of stars in SCUSS catalog are compared with those in the SDSS catalog, and they are highly consistent.

The South Galactic Cap $u$-band Sky Survey (SCUSS) provides a deep $u$-band imaging of about 5000 deg$^2$ in south Galactic cap. It is about 1.5 mag deeper than the SDSS $u$-band. In this paper we evaluate the capability of quasar selection using both SCUSS and SDSS data, based on considerations of the deep SCUSS $u$-band imaging and two-epoch $u$-band variability. We find that the combination of the SCUSS $u$-band and the SDSS $griz$ band allows us to select more faint quasars and more quasars at redshift around 2.2 than the selection only with the SDSS $ugriz$ data. Quasars have significant $u$-band variabilities. The fraction of quasars with large two-epoch variability is much higher than that of stars. The selection by variability can select both low-redshift quasars with ultraviolet excess and mid-redshift ($2 < z <3.5$) quasars where quasar selection by optical colors is inefficient. The above two selections are complementary and make full use of the SCUSS u-band advantages.

We derive the stellar-to-halo mass relations, SHMR, of local blue and red central galaxies separately, as well as the fraction of halos hosting blue/red central galaxies. We find that: 1) the SHMR of central galaxies is segregated by color, with blue centrals having a SHMR above the one of red centrals; at logMh~12, the Ms/Mh ratio of the blue centrals is ~0.05, which is ~1.7 times larger than the value of red centrals. 2) The intrinsic scatters of the SHMRs of red and blue centrals are ~0.14 and ~0.11dex, respectively. The intrinsic scatter of the average SHMR of all central galaxies changes from ~0.20dex to ~0.14dex in the 11.3<logMh<15 range. 3) The fraction of halos hosting blue centrals at Mh=1E11Msun is 87%, but at 2x1E12Msun decays to ~20%, approaching to a few per cents at higher masses. The characteristic mass at which this fraction is the same for blue and red galaxies is Mh~7x1E11Msun. Our results suggest that the SHMR of central galaxies at large masses is shaped by halo mass quenching (likely through shock virial heating and AGN feedback), but group richness also plays an important role: central galaxies living in less dense environments quenched their star formation later or did not quench it yet. At low masses, processes that delay star formation without invoking too strong supernova-driven outflows could explain the high Ms/Mh ratios of blue centrals as compared to those of the scarce red centrals.

We present the [OII] luminosity function measured in the redshift range 0.1<z<1.65 with unprecedented depth and accuracy. Our measurements are based on medium resolution flux-calibrated spectra of emission line galaxies with the FORS2 instrument at VLT and with the SDSS-III/BOSS spectrograph. The FORS2 spectra and the corresponding catalog containing redshifts and line fluxes are released along with this paper. In this work we use a novel method to combine the aforementioned surveys with GAMA, zCOSMOS and VVDS, which have different target selection, producing a consistent weighting scheme to derive the [OII] luminosity function. The measured luminosity function is in good agreement with previous independent estimates. The comparison with two state-of-the-art semi-analytical models is good, which is encouraging for the production of mock catalogs of [OII] flux limited surveys. We observe the bright end evolution over 8.5 Gyr: we measure the decrease of log L* from 42.4 erg/s at redshift 1.44 to 41.2 at redshift 0.165 and we find that the faint end slope flattens when redshift decreases. This measurement confirms the feasibility of the target selection of future baryonic acoustic oscillation surveys aiming at observing [OII] flux limited samples.

We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, $M_\star$, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium $\Lambda$CDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median $M_h|M_\star$ measured from weak lensing depends more sensitively on the lognormal dispersion in $M_\star$ at fixed $M_h$ than it does on the median $M_\star|M_h$. Our measurements suggest an intrinsic dispersion of $\sigma_{\log(M_\star)}\sim 0.15$; 4) Comparing our mass estimates with those in the catalogue, we find that the G$^3$Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass.

Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input non-linear density field, employing the Hamiltonian Markov Chain Monte Carlo (HMC) algorithm combined with Particle Mesh (PM) dynamics. The HMC+PM method is applied to cosmological simulations, and the reconstructed linear density fields are then evolved to the present day with N-body simulations. The constrained simulations so obtained accurately reproduce both the amplitudes and phases of the input simulations at various $z$. Using a PM model with a grid cell size of 0.75 Mpc/h and 40 time-steps in the HMC can recover more than half of the phase information down to a scale k~0.85 h/Mpc at high z and to k~3.4 h/Mpc at z=0, which represents a significant improvement over similar reconstruction models in the literature, and indicates that our model can reconstruct the formation histories of cosmic structures over a large dynamical range. Adopting PM models with higher spatial and temporal resolutions yields even better reconstructions, suggesting that our method is limited more by the availability of computer resource than by principle. Dynamic models of structure evolution adopted in many earlier investigations can induce non-Gaussianity in the reconstructed linear density field, which in turn can cause large systematic deviations in the predicted halo mass function. Such deviations are greatly reduced or absent in our reconstruction.

We study how to measure the galaxy merger rate from the observed close pair count. Using a high-resolution N-body/SPH cosmological simulation, we find an accurate scaling relation between galaxy pair counts and merger rates down to a stellar mass ratio of about 1:30. The relation explicitly accounts for the dependence on redshift (or time), on pair separation, and on mass of the two galaxies in a pair. With this relation, one can easily obtain the mean merger timescale for a close pair of galaxies. The use of virial masses, instead of stellar masses, is motivated by the fact that the dynamical friction time scale is mainly determined by the dark matter surrounding central and satellite galaxies. This fact can also minimize the error induced by uncertainties in modeling star formation in the simulation. Since the virial mass can be read from the well-established relation between the virial masses and the stellar masses in observation, our scaling relation can be easily applied to observations to obtain the merger rate and merger time scale. For major merger pairs (1:1-1:4) of galaxies above a stellar mass of 4*10^10 M_sun/h at z=0.1, it takes about 0.31 Gyr to merge for pairs within a projected distance of 20 kpc/h with stellar mass ratio of 1:1, while the time taken goes up to 1.6 Gyr for mergers with stellar mass ratio of 1:4. Our results indicate that a single timescale usually used in literature is not accurate to describe mergers with the stellar mass ratio spanning even a narrow range from 1:1 to 1:4.

Based on SDSS and South Galactic Cap of u-band Sky Survey (SCUSS) early data, we use star counts method to estimate the Galactic structure parameters in an intermediate latitude with 10,180 main-sequence (MS) stars in absolute magnitude interval of $4 \leq M_r \leq 13$. We divide the absolute magnitude into five intervals:$4 \leq M_r < 5$, $5 \leq M_r < 6$, $6 \leq M_r < 8$, $8 \leq M_r < 10$, $10 \leq M_r \leq 13$, and estimate the Galactic structure parameters in each absolute magnitude interval to explore their possible variation with the absolute magnitude. Our study shows the parameters depend on absolute magnitude. For the thin disk, the intrinsic faint MS stars have large local space density and they tend to stay close to the Galactic plane. A plausible explanation is that faint MS stars with long lifetime experience long gravitational interaction time result in a short scaleheight. However, for the thick disk, the parameters show a complex trend with absolute magnitude, which may imply the complicated original of the thick disk. For the halo, the intrinsic faint MS stars have large local density and small axial ratio, which indicate a flattened inner halo and a more spherical outer halo.

The three-point correlation function (3PCF) for galaxies provides an opportunity to measure the non-Gaussianity generated from nonlinear structure formation and also probes information about galaxy formation and evolution that is generally not available from the two-point correlation function (2PCF). We measure the 3PCF of the Sloan Digital Sky Survey DR7 main sample galaxies in both redshift and projected spaces on scales up to 40Mpc/h. We explore the dependence of the 3PCF on galaxy stellar mass and color in order to constrain the formation and evolution for galaxies of different properties. The study of the dependence on these properties also helps better constrain the relation between galaxy stellar mass and color and the properties of their hosting dark-matter halos. We focus on the study of the reduced 3PCF, Q, defined as the ratio between the 3PCF and the sum of the products of the 2PCFs. We find a very weak stellar mass dependence of Q in both redshift and projected spaces. On small scales, more massive galaxies tend to have slightly higher amplitudes of Q. The shape dependence of Q is also weak on these small scales, regardless of stellar mass and color. The reduced 3PCF has a strong color dependence for the low-mass galaxies, while no significant dependence on color is found for the high-mass galaxies. Low-mass red galaxies have higher amplitudes and stronger shape dependence of the reduced 3PCF than the blue galaxies, implying that these low-mass red galaxies tend to populate filamentary structures. The linear galaxy bias model fails to interpret the color dependence of Q, emphasizing the importance of a nonvanishing quadratic bias parameter in the correct modeling of the galaxy color distribution.

In this paper, we study the non-linear matter power spectrum in a specific family of $f(R)$ models that can reproduce the $\Lambda$CDM background expansion history, using high resolution $N$-body simulations based on the {\sc ecosmog} code. We measure the matter power spectrum in the range of $0.05h{\rm Mpc}^{-1}<k<10h{\rm Mpc}^{-1}$ from simulations for our $f(R)$ models and give theoretical explanations to their behaviour and evolution patterns. We also examine the chameleon mechanism for our models and find that it works throughout the cosmic history in dense regions, for our $f(R)$ models with $|f_{R0}|<10^{-4}$. On the other hand, for models with $|f_{R0}|>10^{-3}$, we find no chameleon screening in dense regions at late times ($z<3$), which means that those models could be ruled out due to the factor-of-$1/3$ enhancement to the strength of Newtonian gravity. We also give the best-fit parameters for a generalised PPF fitting formula which works well for the models studied here.

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.

We perform photometric analyses for a bright early-type galaxy (ETG) sample with 2949 galaxies ($M_{\rm r}<-22.5$ mag) in the redshift range of 0.05 to 0.15, drawn from the SDSS DR7 with morphological classification from Galaxy Zoo 1. We measure the Petrosian and isophotal magnitudes, as well as the corresponding half-light radius for each galaxy. We find that for brightest galaxies ($M_{\rm r}<-23$ mag), our Petrosian magnitudes, and isophotal magnitudes to 25 ${\rm mag/arcsec^2}$ and 1\% of the sky brightness are on average 0.16 mag, 0.20 mag, and 0.26 mag brighter than the SDSS Petrosian values, respectively. In the first case the underestimations are caused by overestimations in the sky background by the SDSS PHOTO algorithm, while the latter two are also due to deeper photometry. Similarly, the typical half-light radii ($r_{50}$) measured by the SDSS algorithm are smaller than our measurements. As a result, the bright-end of the $r$-band luminosity function is found to decline more slowly than previous works. Our measured luminosity densities at the bright end are more than one order of magnitude higher than those of Blanton et al. (2003), and the stellar mass densities at $M_{\ast}\sim 5\times10^{11} M_{\odot}$ and $M_{\ast}\sim 10^{12} M_{\odot}$ are a few tenths and a factor of few higher than those of Bernardi et al. (2010). These results may significantly alleviate the tension in the assembly of massive galaxies between observations and predictions of the hierarchical structure formation model.

We report on the detection of the alignment between galaxies and large-scale structure at z~0.6 based on the CMASS galaxy sample from the Baryon Oscillation Spectroscopy Survey data release 9. We use two statistics to quantify the alignment signal: 1) the alignment two-point correlation function which probes the dependence of galaxy clustering at a given separation in redshift space on the projected angle (theta_p) between the orientation of galaxies and the line connecting to other galaxies, and 2) the cos(2theta)-statistic which estimates the average of cos(2theta_p) for all correlated pairs at given separation. We find significant alignment signal out to about 70 Mpc/h in both statistics. Applications of the same statistics to dark matter halos of mass above 10^12 M_sun/h in a large cosmological simulation show similar scale-dependent alignment signals to the observation, but with higher amplitudes at all scales probed. We show that this discrepancy may be partially explained by a misalignment angle between central galaxies and their host halos, though detailed modeling is needed in order to better understand the link between the orientations of galaxies and host halos. In addition, we find systematic trends of the alignment statistics with the stellar mass of the CMASS galaxies, in the sense that more massive galaxies are more strongly aligned with the large-scale structure.

Analysis of the Pangu N-body simulation validates that the bulk flow of halos follows a Maxwellian distribution which variance is consistent with the prediction of the linear theory of structure formation. We propose that the consistency between the observed bulk velocity and theories should be examined at the effective scale of the radius of a spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the sample window function does. We compared some recently estimated bulk flows from observational samples with the prediction of the \Lambda CDM model we used; some results deviate from expectation at a level of ~ 3\sigma but the discrepancy is not as severe as previously claimed. We show that bulk flow is only weakly correlated with the dipole of the internal mass distribution, the alignment angle between the mass dipole and the bulk flow has a broad distribution peaked at ~ 30-50 deg., and also that the bulk flow shows little dependence on the mass of the halos used in the estimation. In a simulation of box size 1Gpc/h, for a cell of radius 100 Mpc/h the maximal bulk velocity is >500 km/s, dipoles of the environmental mass outside the cell are not tightly aligned with the bulk flow, but are rather located randomly around it with separation angles ~ 20-40 deg. In the fastest cell there is a slightly smaller number of low-mass halos; however halos inside are clustered more strongly at scales > ~ 20 Mpc/h, which might be a significant feature since the correlation between bulk flow and halo clustering actually increases in significance beyond such scales.

It was recently suggested that, compared to its stellar mass (M*), the central stellar velocity dispersion (sigma*) of a galaxy might be a better indicator for its host dark matter halo mass. Here we test this hypothesis by estimating the dark matter halo mass for central alaxies in groups as function of M* and sigma*. For this we have estimated the redshift-space cross-correlation function (CCF) between the central galaxies at given M* and sigma* and a reference galaxy sample, from which we determine both the projected CCF, w_p(r_p), and the velocity dispersion profile (VDP) of satellites around the centrals. A halo mass is then obtained from the average velocity dispersion within the virial radius. At fixed M*, we find very weak or no correlation between halo mass and sigma*. In contrast, strong mass dependence is clearly seen even when sigma* is limited to a narrow range. Our results thus firmly demonstrate that the stellar mass of central galaxies is still a good (if not the best) indicator for dark matter halo mass, better than the stellar velocity dispersion. The dependence of galaxy clustering on sigma* fixed M*, as recently discovered by Wake et al. (2012), may be attributed to satellite galaxies, for which the tidal stripping occurring within halos has stronger effect on stellar mass than on central stellar velocity dispersion.

We investigate the radial number density profile and the abundance distribution of faint satellites around central galaxies in the low redshift universe using the CFHT Legacy Survey. We consider three samples of central galaxies with magnitudes of M_r=-21, -22, and -23 selected from the Sloan Digital Sky Survey (SDSS) group catalog of Yang et al.. The satellite distribution around these central galaxies is obtained by cross-correlating these galaxies with the photometric catalogue of the CFHT Legacy Survey. The projected radial number density of the satellites obeys a power law form with the best-fit logarithmic slope of -1.05, independent of both the central galaxy luminosity and the satellite luminosity. The projected cross correlation function between central and satellite galaxies exhibits a non-monotonic trend with satellite luminosity. It is most pronounced for central galaxies with M_r=-21, where the decreasing trend of clustering amplitude with satellite luminosity is reversed when satellites are fainter than central galaxies by more than 2 magnitudes. A comparison with the satellite luminosity functions in the Milky Way and M31 shows that the Milky Way/M31 system has about twice as many satellites as around a typical central galaxy of similar luminosity. The implications for theoretical models are briefly discussed.

We develop a new code, the Hierarchical Bound-Tracing (HBT for short) code, to find and trace dark matter subhaloes in simulations based on the merger hierarchy of dark matter haloes. Application of this code to a recent benchmark test of finding subhaloes demonstrates that HBT stands as one of the best codes to trace the evolutionary history of subhaloes. The success of the code lies in its careful treatment of the complex physical processes associated with the evolution of subhaloes and in its robust unbinding algorithm with an adaptive source subhalo management. We keep a full record of the merger hierarchy of haloes and subhaloes, and allow growth of satellite subhaloes through accretion from its "satellite-of-satellites", hence allowing mergers among satellites. Local accretion of background mass is omitted, while rebinding of stripped mass is allowed. The justification of these treatments is provided by case studies of the lives of individual subhaloes and by the success in finding the complete subhalo catalogue. We compare our result to other popular subhalo finders and show that HBT is able to well resolve subhaloes in high density environment and keep strict physical track of subhaloes' merger history. This code is fully parallelized and freely available upon request to the authors.

We present measurements of the velocity dispersion profile (VDP) for galaxy groups in the final data release of the Sloan Digital Sky Survey (SDSS). For groups of given mass we estimate the redshift-space cross-correlation function (CCF) with respect to a reference galaxy sample, xi(r_p, pi), the projected CCF, w_p(r_p), and the real-space CCF, xi(r). The VDP is then extracted from the redshift distortion in xi(r_p, pi), by comparing xi(r_p, pi) with xi(r). We find that the velocity dispersion (VD) within virial radius (R_200) shows a roughly flat profile, with a slight increase at radii below ~0.3 R_200 for high mass systems. The average VD within the virial radius, sigma_v, is a strongly increasing function of central galaxy mass. We apply the same methodology to N-body simulations with the concordance Lambda cold dark matter cosmology but different values of the density fluctuation parameter sigma_8, and we compare the results to the SDSS results. We show that the sigma_v-M_* relation from the data provides stringent constraints on both sigma_8 and sigma_ms, the dispersion in log M_* of central galaxies at fixed halo mass. Our best-fitting model suggests sigma_8 = 0.86 +/- 0.03 and sigma_ms = 0.16 +/- 0.03. The slightly higher value of sigma_8 compared to the WMAP7 result might be due to a smaller matter density parameter assumed in our simulations. Our VD measurements also provide a direct measure of the dark matter halo mass for central galaxies of different luminosities and masses, in good agreement with the results obtained by Mandelbaum et al. (2006) from stacking the gravitational lensing signals of the SDSS galaxies.

LAMOST (Large sky Area Multi-Object fiber Spectroscopic Telescope) is a Chinese national scientific research facility operated by National Astronomical Observatories, Chinese Academy of Sciences (NAOC). After two years of commissioning beginning in 2009, the telescope, instruments, software systems and operations are nearly ready to begin the main science survey. Through a spectral survey of millions of objects in much of the northern sky, LAMOST will enable research in a number of contemporary cutting edge topics in astrophysics, such as: discovery of the first generation stars in the Galaxy, pinning down the formation and evolution history of galaxies especially theMilky Way and its central massive black hole, looking for signatures of dark matter distribution and possible sub-structures in the Milky Way halo. To maximize the scientific potential of the facility, wide national participation and international collaboration has been emphasized. The survey has two major components: the LAMOST ExtraGAlactic Survey (LEGAS), and the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE). Until LAMOST reaches its full capability, the LEGUE portion of the survey will use the available observing time, starting in 2012. An overview of the LAMOST project and the survey that will be carried out in next five to six years is presented in this paper. The science plan for the whole LEGUE survey, instrumental specifications, site conditions, the descriptions of the current on-going pilot survey, including its footprints and target selection algorithm, will be presented as separate papers in this volume.

We present measurements of projected autocorrelation functions w_p(r_p) for the stellar mass of galaxies and for their light in the U, B and V bands, using data from the third data release of the DEEP2 Galaxy Redshift Survey and the final data release of the Sloan Digital Sky Survey (SDSS). We investigate the clustering bias of stellar mass and light by comparing these to projected autocorrelations of dark matter estimated from the Millennium Simulations (MS) at z=1 and 0.07, the median redshifts of our galaxy samples. All of the autocorrelation and bias functions show systematic trends with spatial scale and waveband which are impressively similar at the two redshifts. This shows that the well-established environmental dependence of stellar populations in the local Universe is already in place at z=1. The recent MS-based galaxy formation simulation of Guo et al. (2011) reproduces the scale-dependent clustering of luminosity to an accuracy better than 30% in all bands and at both redshifts, but substantially overpredicts mass autocorrelations at separations below about 2 Mpc. Further comparison of the shapes of our stellar mass bias functions with those predicted by the model suggests that both the SDSS and DEEP2 data prefer a fluctuation amplitude of sigma_8 0.8 rather than the sigma_8=0.9 assumed by the MS.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We develop a new method which measures the projected density distribution w_p(r_p)n of photometric galaxies surrounding a set of spectroscopically-identified galaxies, and simultaneously the projected correlation function w_p(r_p) between the two populations. In this method we are able to divide the photometric galaxies into subsamples in luminosity intervals when redshift information is unavailable, enabling us to measure w_p(r_p)n and w_p(r_p) as a function of not only the luminosity of the spectroscopic galaxy, but also that of the photometric galaxy. Extensive tests show that our method can measure w_p(r_p) in a statistically unbiased way. The accuracy of the measurement depends on the validity of the assumption in the method that the foreground/background galaxies are randomly distributed and thus uncorrelated with those galaxies of interest. Therefore, our method can be applied to the cases where foreground/background galaxies are distributed in large volumes, which is usually valid in real observations. We applied our method to data from SDSS including a sample of 10^5 LRGs at z~0.4 and a sample of about half a million galaxies at z~0.1, both of which are cross-correlated with a deep photometric sample drawn from the SDSS. On large scales, the relative bias factor of galaxies measured from w_p(r_p) at z~0.4 depends on luminosity in a manner similar to what is found at z~0.1, which are usually probed by autocorrelations of spectroscopic samples. On scales smaller than a few Mpc and at both z~0.4 and z~0.1, the photometric galaxies of different luminosities exhibit similar density profiles around spectroscopic galaxies at fixed luminosity and redshift. This provides clear support for the assumption commonly-adopted in HOD models that satellite galaxies of different luminosities are distributed in a similar way, following the dark matter distribution within their host halos.

We use seven high-resolution $N$-body simulations to study the correlations among different halo properties (assembly time, spin, shape and substructure), and how these halo properties are correlated with the large-scale environment in which halos reside. The large-scale tidal field estimated from halos above a mass threshold is used as our primary quantity to characterize the large-scale environment, while other parameters, such as the local overdensity and the morphology of large-scale structure, are used for comparison. For halos at a fixed mass, all the halo properties depend significantly on environment, particularly the tidal field. The environmental dependence of halo assembly time is primarily driven by local tidal field. The mass of the unbound fraction in substructure is boosted in strong tidal force region, while the bound fraction is suppressed. Halos have a tendency to spin faster in stronger tidal field and the trend is stronger for more massive halos. The spin vectors show significant alignment with the intermediate axis of the tidal field, as expected from the tidal torque theory. Both the major and minor axes of halos are strongly aligned with the corresponding principal axes of the tidal field. In general, a halo that can accrete more material after the formation of its main halo on average is younger, is more elongated, spins faster, and contains a larger amount of substructure. Higher density environments not only provide more material for halo to accrete, but also are places of stronger tidal field that tends to suppress halo accretion. The environmental dependencies are the results of these two competing effects. The tidal field based on halos can be estimated from observation, and we discuss the implications of our results for the environmental dependence of galaxy properties.

The thermal Sunyaev-Zel'dovich (tSZ) effect directly measures the thermal pressure of free electrons integrated along the line of sight and thus contains valuable information on the thermal history of the universe. However, the redshift information is entangled in the projection along the line of sight. This projection effect severely degrades the power of the tSZ effect to reconstruct the thermal history. We investigate the tSZ tomography technique to recover this otherwise lost redshift information by cross correlating the tSZ effect with galaxies of known redshifts, or alternatively with matter distribution reconstructed from weak lensing tomography. We investigate in detail the 3D distribution of the gas thermal pressure and its relation with the matter distribution, through our adiabatic hydrodynamic simulation and the one with additional gastrophysics including radiative cooling, star formation and supernova feedback. (1) We find a strong correlation between the gas pressure and matter distribution, with a typical cross correlation coefficient r ~ 0.7 at k . 3h/Mpc and z < 2. This tight correlation will enable robust cross correlation measurement between SZ surveys such as Planck, ACT and SPT and lensing surveys such as DES and LSST, at ~20-100{\sigma} level. (2) We propose a tomography technique to convert the measured cross correlation into the contribution from gas in each redshift bin to the tSZ power spectrum. Uncertainties in gastrophysics may affect the reconstruction at ~ 2% level, due to the ~ 1% impact of gastrophysics on r, found in our simulations. However, we find that the same gastrophysics affects the tSZ power spectrum at ~ 40% level, so it is robust to infer the gastrophysics from the reconstructed redshift resolved contribution.

The kinetic Sunyaev Zel'dovich effect (kSZ) effect is a potentially powerful probe to the missing baryons. However, the kSZ signal is overwhelmed by various contaminations and the cosmological application is hampered by loss of redshift information due to the projection effect. We propose a kSZ tomography method to alleviate these problems, with the aid of galaxy spectroscopic redshift surveys. We propose to estimate the large scale peculiar velocity through the 3D galaxy distribution, weigh it by the 3D galaxy density and adopt the product projected along the line of sight with a proper weighting as an estimator of the true kSZ temperature fluctuation $\Theta$. We thus propose to measure the kSZ signal through the $\Hat{\Theta}$-$\Theta$ cross correlation. This approach has a number of advantages (see details in the abstract of the paper). We test the proposed kSZ tomography against non-adiabatic and adiabatic hydrodynamical simulations. We confirm that $\hat{\Theta}$ is indeed tightly correlated with $\Theta$ at $k\la 1h/$Mpc, although nonlinearities in the density and velocity fields and nonlinear redshift distortion do weaken the tightness of the $\hat{\Theta}$-$\Theta$ correlation. We further quantify the reconstruction noise in $\Hat{\Theta}$ from galaxy distribution shot noise. Based on these results, we quantify the applicability of the proposed kSZ tomography for future surveys. We find that, in combination with the BigBOSS-N spectroscopic redshift survey, the PLANCK CMB experiment will be able to detect the kSZ with an overall significance of $\sim 50\sigma$ and further measure its redshift distribution at many redshift bins over $0<z<2$.

The linear growth factor of density perturbations is believed to be a powerful observable of future redshift surveys to probe physical properties of dark energy and to distinguish among gravity theories. We investigate systematic effects on determination of the growth factor f from a measurement of redshift-space distortions. Using N-body simulations we identify dark matter halos over a broad mass range. We compute the power spectra and correlation functions for the halos and then examine how well the redshift distortion parameter beta=f/b can be reconstructed as a function of halo mass. We find that beta measured for a fixed halo mass is generally a function of scale even on large scales, in contrast with the common expectation that beta approaches a constant described by Kaiser's formula on such scales. The scale dependence depends on the halo mass, being stronger for smaller halos. It also cannot be easily explained with the well-known distribution function of the halo peculiar velocities. We demonstrate that the biasing for smaller halos has larger nonlinearity and stochasticity, thus the linear bias assumption becomes worse for smaller halos. Only for massive halos with b>1.5, beta approaches the linear theory prediction on scales of r or pi/k>30Mpc/h. Luminous red galaxies (LRG), targeted by the SDSS-III's BOSS survey, tend to reside in very massive halos. Our results indicate that if the LRG is used for the measurement of redshift distortions, f can be measured unbiasedly. On the other hand, if one considers to use emission line galaxies, which are targeted by the BigBOSS survey and inhabited in halos of a broad mass range, the scale dependence of beta must be taken into account carefully; otherwise one might give incorrect constraints on dark energy or modified gravity theories. We also find that beta reconstructed in Fourier space behaves better than that in configuration space.

We study the dependence of the galaxy contents within halos on the halo formation time using two galaxy formation models, one being a semianalytic model utilizing the halo assembly history from a high resolution N-body simulation and the other being a smoothed particle hydrodynamics simulation including radiative cooling, star formation, and energy feedback from galactic winds. We confirm the finding by Gao et al. that at fixed mass, the clustering of halos depends on the halo formation time, especially for low-mass halos. This age dependence of halo clustering makes it desirable to study the correlation between the occupation of galaxies within halos and the halo age. We find that, in halos of fixed mass, the number of satellite galaxies has a strong dependence on halo age, with fewer satellites in older halos. The youngest one-third of the halos can have an order of magnitude more satellites than the oldest one-third. For central galaxies, in halos that form earlier, they tend to have more stars and thus appear to be more luminous, and the dependence of their luminosity on halo age is not as strong as that of stellar mass. The results can be understood through the star formation history in halos and the merging of satellites onto central galaxies. The age dependence of the galaxy contents within halos would constitute an important ingredient in a more accurate halo-based model of galaxy clustering.

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.

We apply the empirical method built for z=0 in the previous work of Wang et al. to a higher redshift, to link galaxy stellar mass directly with its hosting dark matter halo mass at z~0.8. The relation of the galaxy stellar mass and the host halo mass M_infall is constrained by fitting both the stellar mass function and the correlation functions at different stellar mass intervals of the VVDS observation, where M_infall is the mass of the hosting halo at the time when the galaxy was last the central galaxy. We find that for low mass haloes, their residing central galaxies are less massive at high redshift than those at low redshift. For high mass haloes, central galaxies in these haloes at high redshift are a bit more massive than the galaxies at low redshift. Satellite galaxies are less massive at earlier times, for any given mass of hosting haloes. Fitting both the SDSS and VVDS observations simultaneously, we also propose a unified model of the M_stars-M_infall relation, which describes the evolution of central galaxy mass as a function of time. The stellar mass of a satellite galaxy is determined by the same M_stars-M_infall relation of central galaxies at the time when the galaxy is accreted. With these models, we study the amount of galaxy stellar mass increased from z~0.8 to the present day through galaxy mergers and star formation. Low mass galaxies gain their stellar masses from z~0.8 to z=0 mainly through star formation. For galaxies of higher mass, the increase of stellar mass solely through mergers from z=0.8 can make the massive galaxies a factor ~2 larger than observed at z=0. We can also predict stellar mass functions of redshifts up to z~3, and the results are consistent with the latest observations.

We use a high-resolution $N$-body simulation to investigate the influence of background galaxy properties, including redshift, size, shape and clustering, on the efficiency of forming giant arcs by gravitational lensing of rich galaxy clusters. Two large sets of ray-tracing simulations are carried out for 10 massive clusters at two redshifts, i.e. $z_{\rm l} \sim 0.2$ and 0.3. The virial mass ($M_{\rm vir}$) of the simulated lens clusters at $z\sim0.2$ ranges from $6.8\times10^{14} h^{-1} {M_{\odot}}$ to $1.1\times 10^{15} h^{-1} M_{\odot}$. The information of background galaxies brighter than 25 magnitude in the $I$-band is taken from Cosmological Evolution Survey (COSMOS) imaging data. Around $1.7\times 10^5$ strong lensing realizations with these images as background galaxies have been performed for each set. We find that the efficiency for forming giant arcs for $z_{\rm l}=0.2$ clusters is broadly consistent with observations. The efficiency of producing giant arcs by rich clusters is weakly dependent on the source size and clustering. Our principal finding is that a small proportion ($\sim 1/3$) of galaxies with elongated shapes (e.g. ellipticity $\epsilon=1-b/a>0.5$) can boost the number of giant arcs substantially. Compared with recent studies where a uniform ellipticity distribution from 0 to 0.5 is used for the sources, the adoption of directly observed shape distribution increases the number of giant arcs by a factor of $\sim2$. Our results indicate that it is necessary to account for source information and survey parameters (such as point-spread-function, seeing) to make correct predictions of giant arcs and further to constrain the cosmological parameters.(abridged)

We perform Halo Occupation Distribution (HOD) modeling to interpret small-scale and intermediate-scale clustering of 35,000 luminous early-type galaxies and their cross-correlation with a reference imaging sample of normal L* galaxies in the Sloan Digital Sky Survey. The modeling results show that most of these luminous red galaxies (LRGs) are central galaxies residing in massive halos of typical mass M ~ a few times 10^13 to 10^14 Msun/h, while a few percent of them have to be satellites within halos in order to produce the strong auto-correlations exhibited on smaller scales. The mean luminosity Lc of central LRGs increases with the host halo mass, with a rough scaling relation of Lc \propto M^0.5. The halo mass required to host on average one satellite LRG above a luminosity threshold is found to be about 10 times higher than that required to host a central LRG above the same threshold. We find that in massive halos the distribution of L* galaxies roughly follows that of the dark matter and their mean occupation number scales with halo mass as M^1.5. The HOD modeling results also allows for an intuitive understanding of the scale-dependent luminosity dependence of the cross-correlation between LRGs and L_* galaxies. Constraints on the LRG HOD provide tests to models of formation and evolution of massive galaxies, and they are also useful for cosmological parameter investigations. In one of the appendices, we provide LRG HOD parameters with dependence on cosmology inferred from modeling the two-point auto-correlation functions of LRGs.

We present a general formalism to study the growth of dark matter perturbations when dark energy perturbations and interactions between dark sectors are present. We show that dynamical stability of the growth of structure depends on the type of coupling between dark sectors. By taking the appropriate coupling to ensure the stable growth of structure, we observe that the effect of the dark sectors' interaction overwhelms that of dark energy perturbation on the growth function of dark matter perturbation. Due to the influence of the interaction, the growth index can differ from the value without interaction by an amount within the observational sensibility, which provides a possibility to disclose the interaction between dark sectors through future observations on the growth of large structure.

We study whether the bias factors of galaxies can be unbiasedly recovered from their power spectra and bispectra. We use a set of numerical N-body simulations and construct large mock galaxy catalogs based upon the semi-analytical model of Croton et al. (2006). We measure the reduced bispectra for galaxies of different luminosity, and determine the linear and first nonlinear bias factors from their bispectra. We find that on large scales down to that of the wavenumber k=0.1h/Mpc, the bias factors b1 and b2 are nearly constant, and b1 obtained with the bispectrum method agrees very well with the expected value. The nonlinear bias factor b2 is negative, except for the most luminous galaxies with M<-23 which have a positive b2. The behavior of b2 of galaxies is consistent with the b2 mass dependence of their host halos. We show that it is essential to have an accurate estimation of the dark matter bispectrum in order to have an unbiased measurement of b1 and b2. We also test the analytical approach of incorporating halo occupation distribution to model the galaxy power spectrum and bispectrum. The halo model predictions do not fit the simulation results well on the precision requirement of current cosmological studies.

Using constrained simulations of the local Universe for generic cold dark matter and for 1keV warm dark matter, we investigate the difference in the abundance of dark matter halos in the local environment. We find that the mass function within 20 Mpc/h of the Local Group is ~2 times larger than the universal mass function in the 10^9-10^13 M_odot/h mass range. Imposing the field of view of the on-going HI blind survey ALFALFA in our simulations, we predict that the velocity function in the Virgo-direction region exceeds the universal velocity function by a factor of 3. Furthermore, employing a scheme to translate the halo velocity function into a galaxy velocity function, we compare the simulation results with a sample of galaxies from the early catalog release of ALFALFA. We find that our simulations are able to reproduce the velocity function in the 80-300 km/s velocity range, having a value ~10 times larger than the universal velocity function in the Virgo-direction region. In the low velocity regime, 35-80 km/s, the warm dark matter simulation reproduces the observed flattening of the velocity function. On the contrary, the simulation with cold dark matter predicts a steep rise in the velocity function towards lower velocities; for V_max=35 km/s, it forecasts ~10 times more sources than the ones observed. If confirmed by the complete ALFALFA survey, our results indicate a potential problem for the cold dark matter paradigm or for the conventional assumptions about energetic feedback in dwarf galaxies.

We use a set of numerical N-body simulations to study the large-scale behavior of the reduced bispectrum of dark matter and compare the results with the second-order perturbation theory and the halo models for different halo mass functions. We find that the second-order perturbation theory (PT2) agrees with the simulations fairly well on large scales of k<0.05 h/Mpc, but it shows a signature of deviation as the scale goes down. Even on the largest scale where the bispectrum can be measured reasonably well in our simulations, the inconsistency between PT2 and the simulations appears for the colinear triangle shapes. For the halo model, we find that it can only serve as a qualitative method to help study the behavior of Q on large scales and also on relatively small scales. The failure of second-order perturbation theory will also affect the precise determination of the halo models, since they are connected through the 3-halo term in the halo model. The 2-halo term has too much contribution on the large scales, which is the main reason for the halo model to overpredict the bispectrum on the large scales. Since neither of the models can provide a satisfying description for the bispectrum on scales of about 0.1h/Mpc for the requirement of precision cosmology, we release the reduced bispectrum of dark matter on a large range of scales for future analytical modeling of the bispectrum.

We investigate the orientation correlation of giant elliptical galaxies by measuring the intrinsic ellipticity correlation function of 83,773 luminous red galaxies (LRGs) at redshifts 0.16 -- 0.47 from the Sloan Digital Sky Survey. We have accurately determined the correlation up to 30 $h^{-1}$Mpc. Luminosity dependence of the ellipticity correlation is also detected although the error bars are large, while no evidence is found for its redshift evolution between $z=0.2$ and $z=0.4$. Then we use a cosmological $N$-body simulation to examine misalignment between the central LRGs and their parent dark matter halos. Central and satellite galaxies are assigned to simulated halos by employing a halo occupation distribution model for the LRGs. The ellipticity correlation is predicted to have the same shape as but an amplitude about 4 times higher than our observation if the central LRGs are perfectly aligned with their host halos. This indicates that the central LRG galaxies are preferentially but not perfectly aligned with their host halos. With the assumption that there is a misalignment angle between a central LRG and its host halo which follows a Gaussian distribution with a zero mean and a width $\sigma_\theta$, we obtain a tight constraint on the misalignment parameter, $\sigma_\theta={35.4}^{+4.0}_{-3.3}$ deg. This type of intrinsic ellipticity correlation, if not corrected, can lead to contamination at 5% level to the shear power spectrum in weak lensing surveys of limiting magnitude $R_{AB}=24.5$ if the source central galaxies follow the same misalignment distribution as the LRGs.

Using a high-resolution cosmological $N$-body simulation, we identify the ejected population of subhalos, which are halos at redshift $z=0$ but were once contained in more massive `host' halos at high redshifts. The fraction of the ejected subhalos in the total halo population of the same mass ranges from 9% to 4% for halo masses from $\sim 10^{11}$ to $\sim 10^{12}\msun$. Most of the ejected subhalos are distributed within 4 times the virial radius of their hosts. These ejected subhalos have distinct velocity distribution around their hosts in comparison to normal halos. The number of subhalos ejected from a host of given mass increases with the assembly redshift of the host. Ejected subhalos in general reside in high-density regions, and have a much higher bias parameter than normal halos of the same mass. They also have earlier assembly times, so that they contribute to the assembly bias of dark matter halos seen in cosmological simulations. However, the assembly bias is {\it not} dominated by the ejected population, indicating that large-scale environmental effects on normal halos are the main source for the assembly bias.

We examine whether the gravitational shear--intrinsic ellipticity (GI) correlation function of the luminous red galaxies (LRGs) can be modeled with the distribution function of a misalignment angle advocated recently by Okumura et al. For this purpose, we have accurately measured the GI correlation for the LRGs in the Data Release 6 (DR6) of the Sloan Digital Sky Survey (SDSS), which confirms the results of Hirata et al. who used the DR4 data. By comparing the GI correlation functions in the simulation and in the observation, we find that the GI correlation can be modeled in the current $\Lambda$CDM model if the misalignment follows a Gaussian distribution with a zero mean and a typical misalignment angle $\sigma_\theta=34.9^{+1.9}_{-2.1}$ degrees. We also find a correlation between the axis ratios and intrinsic alignments of LRGs. This effect should be taken into account in theoretical modeling of the GI and intrinsic ellipticity--ellipticity correlations for weak lensing surveys.

Following our previous work(Jiang et al.(2008)), in which we studied the merger time-scale of galaxies in a high-resolution cosmological hydro/N-body simulation, we investigate the potential influence of uncertainties in the numerical implementation of baryonic physics on the merger time-scale. The simulation used in the previous work suffers from the overcooling problem which causes the central galaxies of large halos too massive. This may result in a shortened merger time-scale compared to that in the real universe. We run a similar simulation, but the stellar mass is significantly reduced to model another extreme case of low stellar mass. Our result shows that the merger time-scale is little affected by the star formation recipes, except for the satellites in nearly radial orbits which show a 22 percent higher time-scale in the lower stellar mass case. Since the radial orbits only account for a small part of the satellites' orbits, the fitting formula in Jiang et al.(2008) is still applicable to a reasonable accuracy, nearly immune to the uncertainty in the baryonic physics.

Based on the Sloan Digital Sky Survey DR6 (SDSS) and Millennium Simulation (MS) we investigate the alignment between galaxies and large-scale structure. For this purpose we develop two new statistical tools, namely the alignment correlation function and the cos(2theta)-statistic. The former is a two-dimensional extension of the traditional two-point correlation function and the latter is related to the ellipticity correlation function used for cosmic shear measurements. Both are based on the cross correlation between a sample of galaxies with orientations and a reference sample which represents the large-scale structure. Applied to the SDSS galaxy catalog the alignment correlation function reveals an overabundance of reference galaxies along the major axes of red, luminous (L > L*) galaxies out to projected separations of 60 Mpc/h. No alignment signal is detected for blue galaxies. The cos(2\theta)-statistic yields very similar results. Starting from a MS semi-analytic galaxy catalog we assign an orientation to each red, luminous and central galaxy, based on the central region of the host halo. Alternatively, we use the orientation of the host halo itself. We find a mean projected misalignment between a halo and its central region of ~25 deg. Agreement with the SDSS results is good if the central orientations are used. Using the halo orientations overestimates the observed alignment by more than a factor of 2. The large volume of the MS allows to generate two-dimensional maps of the alignment correlation function which show the reference galaxy distribution to be flattened parallel to the orientations of red luminous galaxies with axis ratios of ~0.5 and ~0.75 for halo and central orientations,respectively. These ratios are almost independent of scale out to 60 Mpc/h.

We use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a \Lambda Cold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite simulation box size in the comparison. We determine the values of the maximum wavenumbers, k^{lim}_{1%} and k^{lim}_{3%}, below which the analytic models and the simulation results agree to within 1 and 3 percent, respectively. We then provide a simple empirical function which describes the convergence regime determined by comparison between our simulations and the analytical models. We find that if we use the Fourier modes within the convergence regime alone, the characteristic scale of baryon acoustic oscillations can be determined within 1% accuracy from future surveys with a volume of a few h^{-3}Gpc^3 at z\sim1 or z\sim3 in the absence of any systematic distortion of the power spectrum.

We develop a new method to reconstruct the cosmic density field from the distribution of dark matter haloes above a certain mass threshold. Our motivation is that well-defined samples of galaxy groups/clusters, which can be used to represent the dark halo population, can now be selected from large redshift surveys of galaxies, and our ultimate goal is to use such data to reconstruct the cosmic density field in the local universe. Our reconstruction method starts with a sample of dark matter haloes above a given mass threshold. Each volume element in space is assigned to the domain of the nearest halo according to a distance measure that is scaled by the virial radius of the halo. The distribution of the mass in and around dark matter haloes of a given mass is modelled using the cross-correlation function between dark matter haloes and the mass distribution within their domains. We use N-body cosmological simulations to show that the density profiles required in our reconstruction scheme can be determined reliably from large cosmological simulations, and that our method can reconstruct the density field accurately using haloes with masses down to $\sim 10^{12}\msun$ (above which samples of galaxy groups can be constructed from current large redshift surveys of galaxies). Working in redshift space, we demonstrate that the redshift distortions due to the peculiar velocities of haloes can be corrected in an iterative way. We also describe some applications of our method.

Galaxy formation inside dark matter halos, as well as the halo formation itself, can be affected by large-scale environments. Evaluating the imprints of environmental effects on galaxy clustering is crucial for precise cosmological constraints with data from galaxy redshift surveys. We investigate such an environmental impact on both real-space and redshift-space galaxy clustering statistics using a semi-analytic model derived from the Millennium Simulation. We compare clustering statistics from original SAM galaxy samples and shuffled ones with environmental influence on galaxy properties eliminated. Among the luminosity-threshold samples examined, the one with the lowest threshold luminosity (~0.2L_*) is affected by environmental effects the most, which has a ~10% decrease in the real-space two-point correlation function (2PCF) after shuffling. By decomposing the 2PCF into five different components based on the source of pairs, we show that the change in the 2PCF can be explained by the age and richness dependence of halo clustering. The 2PCFs in redshift space are found to change in a similar manner after shuffling. If the environmental effects are neglected, halo occupation distribution modeling of the real-space and redshift-space clustering may have a less than 6.5% systematic uncertainty in constraining beta from the most affected SAM sample and have substantially smaller uncertainties from the other, more luminous samples. We argue that the effect could be even smaller in reality. In the Appendix, we present a method to decompose the 2PCF, which can be applied to measure the two-point auto-correlation functions of galaxy sub-samples in a volume-limited galaxy sample and their two-point cross-correlation functions in a single run utilizing only one random catalog.

It has long been known that galaxy interactions are associated with enhanced star formation. In a companion paper, we explored this connection by applying a variety of statistics to SDSS data. In particular, we showed that specific star formation rates of galaxies are higher if they have close neighbours. Here we apply exactly the same techniques to AGN in the survey, showing that close neighbours are not associated with any similar enhancement of nuclear activity. Star formation is enhanced in AGN with close neighbours in exactly the same way as in inactive galaxies, but the accretion rate onto the black hole, as estimated from the extinction-corrected [O III] luminosity, is not influenced by the presence or absence of companions. Previous work has shown that galaxies with more strongly accreting black holes contain more young stars in their inner regions. This leads us to conclude that star formation induced by a close companion and star formation associated with black hole accretion are distinct events. These events may be part of the same physical process, for example a merger, provided they are separated in time. In this case, accretion onto the black hole and its associated star formation would occur only after the two interacting galaxies have merged. The major caveat in this work is our assumption that the extinction-corrected [O III] luminosity is a robust indicator of the bolometric luminosity of the central black hole. It is thus important to check our results using indicators of AGN activity at other wavelengths.

We investigate the clustering properties of a complete sample of 10^5 star-forming galaxies drawn from the SDSS DR4. On scales less than 100 kpc, the amplitude of the correlation function exhibits a strong dependence on the specific star formation rate of the galaxy. We interpret this as the signature of enhanced star formation induced by tidal interactions. We then explore how the average star formation rate in a galaxy is enhanced as the projected separation r_p between the galaxy and its companions decreases. We find that the enhancement depends strongly on r_p, but very weakly on the relative luminosity of the companions. The enhancement is also stronger in low mass galaxies than in high mass galaxies. In order to explore whether a tidal interaction is not only sufficient, but also necessary to trigger enhanced star formation in a galaxy, we compute background subtracted neighbour counts for the galaxies in our sample. The average number of close neighbours around galaxies with low to average values of SFR/M* is close to zero. At the highest specific star formation rates, however, more than 40% of the galaxies in our sample have a companion within a projected radius of 100 kpc. Visual inspection of the highest SFR/M* galaxies without companions reveals that more than 50% of these are clear interacting or merging systems. We conclude that tidal interactions are the dominant trigger of enhanced star formation in the most strongly star-forming systems. Finally, we find clear evidence that tidal interactions not only lead to enhanced star formation in galaxies, but also cause structural changes such as an increase in concentration.

Based on a cosmological N-body simulation we analyze spatial and kinematic alignments of satellite halos within six times the virial radius of group size host halos (Rvir). We measure three different types of spatial alignment: halo alignment between the orientation of the group central substructure (GCS) and the distribution of its satellites, radial alignment between the orientation of a satellite and the direction towards its GCS, and direct alignment between the orientation of the GCS and that of its satellites. In analogy we use the directions of satellite velocities and probe three further types of alignment: the radial velocity alignment between the satellite velocity and connecting line between satellite and GCS, the halo velocity alignment between the orientation of the GCS and satellite velocities and the auto velocity alignment between the satellites orientations and their velocities. We find that satellites are preferentially located along the major axis of the GCS within at least 6 Rvir (the range probed here). Furthermore, satellites preferentially point towards the GCS. The most pronounced signal is detected on small scales but a detectable signal extends out to 6 Rvir. The direct alignment signal is weaker, however a systematic trend is visible at distances < 2 Rvir. All velocity alignments are highly significant on small scales. Our results suggest that the halo alignment reflects the filamentary large scale structure which extends far beyond the virial radii of the groups. In contrast, the main contribution to the radial alignment arises from the adjustment of the satellite orientations in the group tidal field. The projected data reveal good agreement with recent results derived from large galaxy surveys. (abridged)

We study galaxy mergers using a high-resolution cosmological hydro/N-body simulation with star formation, and compare the measured merger timescales with theoretical predictions based on the Chandrasekhar formula. In contrast to Navarro et al., our numerical results indicate, that the commonly used equation for the merger timescale given by Lacey and Cole, systematically underestimates the merger timescales for minor mergers and overestimates those for major mergers. This behavior is partly explained by the poor performance of their expression for the Coulomb logarithm, \ln (m_pri/m_sat). The two alternative forms \ln (1+m_pri/m_sat) and 1/2\ln [1+(m_pri/m_sat)^2] for the Coulomb logarithm can account for the mass dependence of merger timescale successfully, but both of them underestimate the merger time scale by a factor 2. Since \ln (1+m_pri/m_sat) represents the mass dependence slightly better we adopt this expression for the Coulomb logarithm. Furthermore, we find that the dependence of the merger timescale on the circularity parameter \epsilon is much weaker than the widely adopted power-law \epsilon^{0.78}, whereas 0.94*{\epsilon}^{0.60}+0.60 provides a good match to the data. Based on these findings, we present an accurate and convenient fitting formula for the merger timescale of galaxies in cold dark matter models.

We use high-resolution N-body numerical simulations to study the number of predicted large-separation multiply-imaged systems produced by clusters of galaxies in the SDSS photometric and spectroscopic quasar samples. We incorporate the condensation of baryons at the centre of clusters by (artificially) adding a brightest central galaxy (BCG) as a truncated isothermal sphere. We make predictions in two flat cosmological models: a LCDM model with a matter density $\Omega_0=0.3$, and $\sigma_8=0.9$ (LCDM0), and a model favoured by the WMAP three-year data with $\Omega_0=0.238$, and $\sigma_8=0.74$ (WMAP3). We found that the predicted multiply-imaged quasars with separation >10" is about 6.2 and 2.6 for the SDSS photometric (with an effective area 8000 deg$^2$) and spectroscopic (with an effective area 5000 deg$^2$) quasar samples respectively in the LCDM0 model; the predicted numbers of large-separation lensed quasars agree well with the observations. These numbers are reduced by a factor of 7 or more in the WMAP3 model, and are consistent with data at < 8% level. The predicted cluster lens redshift peaks around redshift 0.5, and 90% are between 0.3 and 1. We find that the BCG creates a central circular region, comparable to the Einstein ring of the BCG, where the central image disappears in the usual three-image and five-image configurations. If we include four image systems as an extreme case of five-image systems (with an infinitely demagnified central image), we find that 68% of the central images are fainter by a factor of 100 than the brightest image, and about 80% are within 1.5"of the BCG.

Using a large galaxy group catalogue based on the Sloan Digital Sky Survey Data Release 4 we measure three different types of intrinsic galaxy alignment within groups: halo alignment between the orientation of the brightest group galaxies (BGG) and the distribution of its satellite galaxies, radial alignment between the orientation of a satellite galaxy and the direction towards its BGG, and direct alignment between the orientation of the BGG and that of its satellites. In agreement with previous studies we find that satellite galaxies are preferentially located along the major axis. In addition, on scales r < 0.7 Rvir we find that red satellites are preferentially aligned radially with the direction to the BGG. The orientations of blue satellites, however, are perfectly consistent with being isotropic. Finally, on scales r < 0.1 \Rvir, we find a weak but significant indication for direct alignment between satellites and BGGs. We briefly discuss the implications for weak lensing measurements.

Recent studies have shown that the distribution of satellite galaxies is preferentially aligned with the major axis of their central galaxy. The strength of this alignment has been found to depend strongly on the colours of the satellite and central galaxies, and only weakly on the mass of the halo in which the galaxies reside. In this paper we study whether these alignment signals, and their dependence on galaxy and halo properties, can be reproduced in a hierarchical structure formation model of a $\Lambda$CDM concordance cosmology. To that extent we use a large $N$-body simulation which we populate with galaxies following a semi-analytical model for galaxy formation. We find that if the orientation of the central galaxy is perfectly aligned with that of its dark matter halo, then the predicted central-satellite alignment signal is much stronger than observed. If, however, the minor axis of a central galaxy is perfectly aligned with the angular momentum vector of its dark matter halo, we can accurately reproduce the observed alignment strength as function of halo mass and galaxy color. Although this suggests that the orientation of central galaxies is governed by the angular momentum of their dark matter haloes, we emphasize that any other scenario in which the minor axes of central galaxy and halo are misaligned by $\sim 40^{\circ}$ (on average) will match the data equally well. Finally, we show that dependence of the alignment strength on the color of the central galaxy is most likely an artefact due to interlopers in the group catalogue. The dependence on the color of the satellite galaxies, on the other hand, is real and owes to the fact that red satellites are associated with subhaloes that were more massive at their time of accretion.

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.

We examine the age-dependence of dark matter halo clustering in an unprecedented accuracy using a set of 7 high-resolution cosmological simulations each with $N=1024^3$ particles. We measure the bias parameters for halos over a large mass range using the cross-power-spectrum method that can effectively suppress the random noise even in the sparse sampling of the most massive halos. This enables us to find, for the first time, that younger halos are more strongly clustered than older ones for halo masses $M>10M_{\ast}$, where $M_\ast$ is the characteristic nonlinear mass scale. For $M<M_{\ast}$, our results confirm the previous finding of Gao et al. that older halos are clustered more strongly than the younger ones. We also study the halo bias as a function of halo concentration, and find that the concentration dependence is weaker than the age dependence for $M<M_{\ast}$, but stronger for $M\ga 50 M_{\ast}$. The accurate and robust measurement of the age dependences of halo bias points to a limitation of the simple excursion set theory which predicts that the formation and structure of a halo of given mass is independent of its environment.

We consider nonlinear biasing models of galaxies with particular attention to a correlation between linear and quadratic biasing coefficients, b_1 and b_2. We first derive perturbative expressions for b_1 and b_2 in halo and peak biasing models. Then we compute power spectra and bispectra of dark matter particles and halos using N-body simulation data and of volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their b_1 and b_2. We find that the values of those coefficients at linear regimes (k<0.2h/Mpc) are fairly insensitive to the redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of bispectra, Q, for equilateral triangles, are insensitive to the values of b_1 implying that b_2 indeed correlates with b_1. The present results explain the previous finding of Kayo et al. (2004) for the hierarchical relation of three-point correlation functions of SDSS galaxies. While the relations between b_1 and b_2 are quantitatively different for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the biasing due to the gravity in primordial Gaussian density fields.

We use a high-resolution $N$-body simulation to study how the formation of cold dark matter (CDM) halos is affected by their environments, and how such environmental effects produce the age-dependence of halo clustering observed in recent $N$-body simulations. We estimate, for each halo selected at redshift $z=0$, an `initial' mass $M_{\rm i}$ defined to be the mass enclosed by the largest sphere which contains the initial barycenter of the halo particles and within which the mean linear density is equal to the critical value for spherical collapse at $z=0$. For halos of a given final mass, $M_{\rm h}$, the ratio $M_{\rm i}/M_{\rm h}$ has large scatter, and the scatter is larger for halos of lower final masses. Halos that form earlier on average have larger $M_{\rm i}/M_{\rm h}$, and so correspond to higher peaks in the initial density field than their final masses imply. Old halos are more strongly clustered than younger ones of the same mass because their initial masses are larger. The age-dependence of clustering for low-mass halos is entirely due to the difference in the initial/final mass ratio. Low-mass old halos are almost always located in the vicinity of big structures, and their old ages are largely due to the fact that their mass accretions are suppressed by the hot environments produced by the tidal fields of the larger structure. The age-dependence of clustering is weaker for more massive halos because the heating by large-scale tidal fields is less important.

Using a set of high-resolution N-body/SPH cosmological simulations with identical initial conditions but run with different numerical setups, we investigate the influence of baryonic matter on the mass distribution of dark halos when radiative cooling is NOT included. We compare the concentration parameters of about 400 massive halos with virial mass from $10^{13}$ \Msun to $7.1 \times 10^{14}$ \Msun. We find that the concentration parameters for the total mass and dark matter distributions in non radiative simulations are on average larger by ~3% and 10% than those in a pure dark matter simulation. Our results indicate that the total mass density profile is little affected by a hot gas component in the simulations. After carefully excluding the effects of resolutions and spurious two-body heating between dark matter and gas particles, we conclude that the increase of the dark matter concentration parameters is due to interactions between baryons and dark matter. We demonstrate this with the aid of idealized simulations of two-body mergers. The results of individual halos simulated with different mass resolutions show that the gas profiles of densities, temperature and entropy are subjects of mass resolution of SPH particles. In particular, we find that in the inner parts of halos, as the SPH resolution increases the gas density becomes higher but both the entropy and temperature decrease.

Galaxy-galaxy lensing is rapidly becoming one of the most promising means to accurately measure the average relation between galaxy properties and halo mass. In order to obtain a signal of sufficient signal-to-noise, one needs to stack many lens galaxies according to their property of interest, such as luminosity or stellar mass. Since such a stack consists of both central and satellite galaxies, which contribute very different lensing signals, the resulting shear measurements can be difficult to interpret. In the past, galaxy-galaxy lensing studies have either completely ignored this problem, have applied rough isolation criteria in an attempt to preferentially select `central' galaxies, or have tried to model the contribution of satellites explicitely. However, if one is able to {\it a priori} split the galaxy population in central and satellite galaxies, one can measure their lensing signals separately. This not only allows a much cleaner measurement of the relation between halo mass and their galaxy populations, but also allows a direct measurement of the sub-halo masses around satellite galaxies. In this paper, we use a realistic mock galaxy redshift survey to show that galaxy groups, properly selected from large galaxy surveys, can be used to accurately split the galaxy population in centrals and satellites. Stacking the resulting centrals according to their group mass, estimated from the total group luminosity, allows a remarkably accurate recovery of the masses and density profiles of their host haloes. In addition, stacking the corresponding satellite galaxies according to their projected distance from the group center yields a lensing signal that can be used to accurate measure the masses of both sub-haloes and host haloes. (Abridged)

We have analyzed the clustering of ~ 90,000 narrow-line AGN drawn from the Data Release 4 (DR4) of the SDSS. We compute the cross-correlation between AGN and a reference sample of galaxies, and compare this to results for control samples of inactive galaxies matched simultaneously in redshift,stellar mass,concentration, velocity dispersion and the 4000A break strength. We also compare near-neighbour counts around AGN and around the control galaxies. On scales larger than a few Mpc, AGN have almost the same clustering amplitude as the control sample. This demonstrates that AGN host galaxies and inactive galaxies populate dark matter halos of similar mass.On scales between 100kpc and 1Mpc,AGN are clustered more weakly than the control galaxies. We use mock catalogues constructed from high-resolution N-body simulations to interpret this anti-bias, showing that the observed effect is easily understood if AGN are preferentially located at the centres of their dark matter halos. On scales less than 70 kpc, AGN cluster marginally more strongly than the control sample, but the effect is weak. When compared to the control sample, we find that only one in a hundred AGN has an extra neighbour within a radius of 70 kpc. This excess increases as a function of the accretion rate onto the black hole, but it does not rise above the few percent level. Although interactions between galaxies may be responsible for triggering nuclear activity in a minority of nearby AGN, some other mechanism is required to explain the activity seen in the majority of the objects in our sample. (abridged)

We propose a new smoothing method for obtaining surface densities from discrete particle positions from numerical simulations. This is an essential step for many applications in gravitational lensing. This method is based on the ``scatter'' interpretation of the discrete density field in the Smoothed Particle Hydrodynamics. We use Monte Carlo simulations of uniform density fields and one isothermal ellipsoid to empirically derive the noise properties, and best smoothing parameters (such as the number of nearest neighbors used). A cluster from high-resolution simulations is then used to assess the reality of high-order singularities such as swallowtails and butterflies in caustics, which are important for the interpretation of substructures in gravitational lenses. We also compare our method with the Delaunay tesselation field estimator using the galaxy studied by Bradac et al. (2004), and find good agreements. We show that higher order singularities are not only connected with bound subhaloes but also with the satellite streams. However, the presence of high-order singularities are sensitive to not only the fluctuation amplitude of the surface density, but also the detailed form of the underlying smooth lensing potential (such as ellipticity and external shear).

We use high-resolution $N$-body simulations to investigate the optical depth of giant arcs with length-to-width ratio larger than 7.5 and 10 in the `standard' $\LCDM$ model with $\sigma_8=0.9$ and $\Omega_{\rm m,0}=0.3$ and a model based on three-year Wilkinson Microwave Anisotropy Probe (WMAP) data. We find that, in dark-matter only simulations, the lensing probability in the three-year WMAP model (with $\sigma_8=0.74$ and $\Omega_{\rm m,0}=0.238)$ decreases by a factor of $\sim 6$ compared with that in the `standard' $\LCDM$ model. The effects of baryonic cooling, star formation and feedbacks are uncertain, but we argue that baryons will only increase the the lensing cross-section by a moderate factor, $\sim 2$. We conclude that the low central value of $\sigma_8$ and $\Omega_{\rm m,0}$ preferred by the WMAP three-year data may be too low to be compatible with observations if conventional assumptions of the background source population are correct.

We study whether hierarchical galaxy formation in a concordance $\Lambda$CDM universe can produce enough massive and red galaxies compared to the observations. We implement a semi-analytical model in which the central black holes gain their mass during major mergers of galaxies and the energy feedback from active galaxy nuclei (AGN) suppresses the gas cooling in their host halos. The energy feedback from AGN acts effectively only in massive galaxies when supermassive black holes have been formed in the central bulges. Compared with previous models without black hole formation, our model predicts more massive and luminous galaxies at high redshift, agreeing with the observations of K20 up to $z\sim 3$. Also the predicted stellar mass density from massive galaxies agrees with the observations of GDDS. Because of the energy feedback from AGN, the formation of new stars is stopped in massive galaxies with the termination of gas cooling and these galaxies soon become red with color $R-K>$5 (Vega magnitude), comparable to the Extremely Red Objects (EROs) observed at redshift $z\sim$1-2. Still the predicted number density of very EROs is lower than observed at $z\sim 2$, and it may be related to inadequate descriptions of dust extinction, star formation history and AGN feedback in those luminous galaxies.

(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 use galaxy groups selected from the Sloan Digital Sky Survey to examine the alignment between the orientation of the central galaxy (defined as the brightest group member) and the distribution of satellite galaxies. By construction, we therefore only address the alignment on scales smaller than the halo virial radius. We find a highly significant alignment of satellites with the major axis of their central galaxy. This is in qualitative agreement with the recent study of Brainerd (2005), but inconsistent with several previous studies who detected a preferential minor axis alignment. The alignment strength in our sample is strongest between red central galaxies and red satellites. On the contrary, the satellite distribution in systems with a blue central galaxy is consistent with isotropic. We also find that the alignment strength is stronger in more massive haloes and at smaller projected radii from the central galaxy. In addition, there is a weak indication that fainter (relative to the central galaxy) satellites are more strongly aligned. We present a detailed comparison with previous studies, and discuss the implications of our findings for galaxy formation.

Future weak lensing measurements of cosmic shear will reach such high accuracy that second order effects in weak lensing modeling, like the influence of baryons on structure formation, become important. We use a controlled set of state of the art cosmological simulations to quantify this effect by comparing pure N-body dark matter runs with corresponding hydrodynamical simulations, carried out both in non-radiative, and in dissipative form with cooling and star formation. In both hydrodynamical simulations, the clustering of the gas is suppressed while that of dark matter is boosted at scales k>1 h/Mpc. Despite this counterbalance between dark matter and gas, the clustering of the total matter is suppressed by up to 1 percent at 1<k<10 h/Mpc, while for k ~ 20 h/Mpc it is boosted, up to 2 percent in the non-radiative run and 10 percent in the run with star formation. The stellar mass formed in the latter is highly biased relative to the dark matter in the pure N-body simulation. Using our power spectrum measurements to predict the effect of baryons on the weak lensing signal at 100<l<10000, we find that baryons may change the lensing power spectrum by less than 0.5 percent at l<1000, but by 1 to 10 percent at 1000<l<10000. The size of the effect exceeds the predicted accuracy of future lensing power spectrum measurements and will likely be detected. Precise determinations of cosmological parameters with weak lensing, and studies of small-scale fluctuations and clustering, therefore rely on properly including baryonic physics.

We use high-resolution numerical simulations to study the physical properties of subhalos when they merge into their host halos. An improved algorithm is used to identify the subhalos. We then examine their spatial and velocity distributions in spherical and triaxial halo models. We find that the accretion of satellites preferentially occurs along the major axis and perpendicular to the spin axis of the host halo. Furthermore, the massive subhalos show a stronger preference to be accreted along the major axis of the host halo than the low-mass ones. Approximate fitting formulae are provided for the physical properties of subhalos. Combined with analytical and semi-analytic techniques, these empirical formulae provide a useful basis for studying the subsequent evolution of subhalos and satellite galaxies in their hosts. Future studies should however account for satellites that may not be undergoing the first infall in their evolution.

We use high-resolution N-body simulations to study the galaxy-cluster cross-sections and the abundance of giant arcs in the $\Lambda$CDM model. Clusters are selected from the simulations using the friends-of-friends method, and their cross-sections for forming giant arcs are analyzed. The background sources are assumed to follow a uniform ellipticity distribution from 0 to 0.5 and to have an area identical to a circular source with diameter $1\arcsec$. We find that the optical depth scales as the source redshift approximately as $\tau_{1''} = 2.25 \times 10^{-6}/[1+(\zs/3.14)^{-3.42}]$ ($0.6<\zs<7$). The amplitude is about 50% higher for an effective source diameter of $0.5\arcsec$. The optimal lens redshift for giant arcs with the length-to-width ratio ($L/W$) larger than 10 increases from 0.3 for $\zs=1$, to 0.5 for $\zs=2$, and to 0.7-0.8 for $\zs>3$. The optical depth is sensitive to the source redshift, in qualitative agreement with Wambsganss et al. (2004). However, our overall optical depth appears to be only $\sim$ 10% to 70% of those from previous studies. The differences can be mostly explained by different power spectrum normalizations ($\sigma_8$) used and different ways of determining the $L/W$ ratio. Finite source size and ellipticity have modest effects on the optical depth. We also found that the number of highly magnified (with magnification $|\mu|>10$) and ``undistorted'' images (with $L/W<3$) is comparable to the number of giant arcs with $|\mu|>10$ and $L/W>10$. We conclude that our predicted rate of giant arcs may be lower than the observed rate, although the precise `discrepancy' is still unclear due to uncertainties both in theory and observations.

We investigate the intrinsic alignments of dark halo substructures with their host halo major-axis orientations both analytically and numerically. Analytically, we derive the probability density distribution of the angles between the minor axes of the substructures and the major axes of their host halos from the physical principles, under the assumption that the substructure alignment on galaxy scale is a consequence of the tidal fields of the host halo gravitational potential. Numerically, we use a sample of four cluster-scale halos and their galaxy-scale substructures from recent high-resolution N-body simulations to measure the probability density distribution. We compare the numerical distribution with the analytic prediction, and find that the two results agree with each other very well. We conclude that our analytic model provides a quantitative physical explanation for the intrinsic alignment of dark halo substructures. We also discuss the possibility of discriminating our model from the anisotropic infall scenario by testing it against very large N-body simulations in the future.

We present an analytic expression for the axis ratio distribution of triaxial dark matter halos driven from physical principles. Adopting the picture of triaxial collapse based on the Zel'dovich approximation, we derive analytically both the minor-to-major and the conditional intermediate-to-major axis ratio distributions, and examine how they depend on the halo mass, redshift, and cosmology. Our analytic model is tested against the simulation data given by Jing & Suto in 2002, and found to reproduce the conditional intermediate-to-major axis ratio distribution successfully and the minor-to-major axis ratio distribution approximately. However, the trends of our analytic axis-ratio distributions with mass and redshift are opposite to what is found in N-body simulations. This failure of our analytic model puts a limitation on analytic approaches based on the Lagrangian theory to the halo ellipticity. Given the overall agreement with the simulation results, our model provides a new theoretical step toward using the axis-ratio distribution of dark halos as a cosmological probe. We also discuss several possibilities to improve the model.

We determine the cross-correlation function between galaxies and galaxy groups, using both the Two-Degree Field Galaxy Redshift Survey (2dFGRS) and the Sloan Digital Sky Survey (SDSS). We study the cross-correlation as a function of group mass, and as a function of the luminosity, stellar mass, colour, spectral type and specific star formation rate of the galaxies. All these cross-correlation functions show a clear transition from the `1-halo' to the `2-halo' regimes on a scale comparable to the virial radius of the groups in consideration. On scales larger than the virial radius, all cross-correlation functions are roughly parallel, consistent with the linear bias model. In particular, the large scale correlation amplitudes are higher for more massive groups, and for brighter and redder galaxies. In the `1-halo' regime, the cross-correlation function depends strongly on the definition of the group center. We consider both a luminosity-weighted center (LWC) and a center defined by the location of the brightest group galaxy (BGC). With the first definition, the bright early-type galaxies in massive groups are found to be more centrally concentrated than the fainter, late-type galaxies. Using the BGC, and excluding the brightest galaxy from the cross correlation analysis, we only find significant segregation in massive groups ($M \gta 10^{13}h^{-1}\msun$) for galaxies of different spectral types (or colours or specific star formation rates). In haloes with masses $\la 10^{13}h^{-1}\msun$, there is a significant deficit of bright satellite galaxies. Comparing the results from the 2dFGRS with those obtained from realistic mock samples, we find that the distribution of galaxies in groups is much less concentrated than dark matter haloes predicted by the current $\Lambda$CDM model. (Abridged)

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 study the spatial distribution of satellite galaxies by assuming that they follow the dark matter distribution. This assumption is supported by semi-analytical studies based on high-resolution numerical simulations. We find that for a Milky-Way type halo, if only a dozen satellite galaxies are observed, then they can lie on a ``great'' disk with an rms height of about 40 kpc. The normal to the plane is roughly isotropic on the sky. These results are consistent with the observed properties of the satellite galaxies in the Milky Way. If, however, the satellite galaxies follow the distribution of substructure selected by present mass, then great disks similar to the one in the Milky Way are rare and difficult to reproduce, in agreement with the conclusion reached by Kroupa et al. (2004).

The brightest galaxy in a dark matter halo is expected to reside at rest at the center of the halo. In this paper we test this `Central Galaxy Paradigm' using group catalogues extracted from the 2dFGRS and the SDSS. For each group we compute a parameter, R, defined as the difference between the velocity of the brightest group galaxy and the average velocity of the other group members (hereafter satellites), normalized by the unbiased estimator of the velocity dispersion of the satellite galaxies. Using detailed mock galaxy redshift surveys, we show that the central galaxy paradigm is inconsistent with the data at high confidence, and that instead the brightest halo galaxies have a specific kinetic energy that is about 25 percent of that of the satellites. This indicates that either central galaxies reside at the minimum of the dark matter potential, but that the halo itself is not yet fully relaxed, or, that the halo is relaxed, but that the central galaxy oscillates in its potential well. The former is consistent with the fact that we find a weak hint that the velocity bias of brightest halo galaxies is larger in more massive haloes, while the latter may be indicative of cored, rather than cusped, dark matter haloes. We discuss several implications of these findings, including mass estimates based on satellite kinematics, strong gravitational lensing, halo occupation models, and the frequency and longevity of lopsidedness in disk galaxies.

We have accurately measured the bispectrum for four scale-free models of structure formation with the spectral index $n=1$, 0, -1, and -2. The measurement is based on a new method that can effectively eliminate the alias and numerical artifacts, and reliably extend the analysis into the strongly non-linear regime. The work makes use of a set of state-of-the art N-body simulations that have significantly increased the resolution range compared with the previous studies on the subject. With these measured results, we demonstrated that the measured bispectrum depends on the shape and size of $k$-triangle even in the strongly nonlinear regime. It increases with wavenumber and decreases with the spectral index. These results are in contrast with the hypothesis that the reduced bispectrum is a constant in the strongly non-linear regime. We also show that the fitting formula of Scoccimarro & Frieman (1999) does not describe our simulation results well (with a typical error about 40 percent). In the end, we present a new fitting formula for the reduced bispectrum that is valid for $-2 \leq n \leq 0$ with a typical error of 10 percent only.

We study the occupation statistics of galaxies in dark matter haloes using galaxy groups identified from the 2-degree Field Galaxy Redshift Survey with the halo-based group finder of Yang et al. The occupation distribution is considered separately for early and late type galaxies, as well as in terms of central and satellite galaxies. The mean luminosity of the central galaxies scales with halo mass approximately as $L_c\propto M^{2/3}$ for haloes with masses $M<10^{13}h^{-1}\msun$, and as $L_c\propto M^{1/4}$ for more massive haloes. The characteristic mass of $10^{13} h^{-1} \Msun$ is consistent with the mass scale where galaxy formation models suggest a transition from efficient to inefficient cooling. Another characteristic halo mass scale, $M\sim 10^{11} h^{-1}\msun$, which cannot be probed directly by our groups, is inferred from the conditional luminosity function (CLF) that matches the observed galaxy luminosity function and clustering. For a halo of given mass, the distribution of $L_c$ is rather narrow. The satellite galaxies are found to follow a Poissonian number distribution. The central galaxies in low-mass haloes are mostly late type galaxies, while those in massive haloes are almost all early types. We also measure the CLF of galaxies in haloes of given mass. Over the mass range that can be reliably probed with the present data ($13.3 \lta {\rm log}[M/(h^{-1}\Msun)] \lta 14.7$), the CLF is reasonably well fit by a Schechter function. Contrary to recent claims based on semi-analytical models of galaxy formation, the presence of central galaxies does not show up as a strong peak at the bright end of the CLF. (Abridged)

We present a simple model that accurately describes various statistical properties of peculiar velocities of dark matter halos. We pay particular attention to the following two effects; first, the evolution of the halo peculiar velocity depends on the local matter density, instead of the global density. Second, dark matter halos are biased tracers of the underlying mass distribution, thus halos tend to be located preferentially at high density regions. For the former, we develop an empirical model calibrated with N-body simulations, while for the latter, we use a conventional halo bias models based on the extended Press-Schechter model combined with an empirical log-normal probability distribution function of the mass density distribution. We find that compared with linear theory, the present model significantly improves the accuracy of predictions of statistical properties of the halo peculiar velocity field including the velocity dispersion, the probability distribution function, and the pairwise velocity dispersion at large separations. Thus our model predictions may be useful in analyzing future observations of the peculiar velocities of galaxy clusters.

We analyze the 2-point correlation function (2PCF) of galaxy groups identified from the 2-degree Field Galaxy Redshift Survey with the halo based group finder recently developed by Yang et al (2004b). With this group catalogue we are able to estimate the 2PCFs for systems ranging from isolated galaxies to rich clusters of galaxies. The real-space correlation length obtained for these systems ranges from $\sim 4\mpch$ to $\sim 15\mpch$, respectively. The observed correlation amplitude (and the corresponding bias factor) as a function of group abundance is well reproduced by associating galaxy groups with dark matter haloes in the standard $\Lambda$CDM model. Redshift distortions are clearly detected in the redshift-space correlation function, the degree of which is consistent with the assumption of gravitational clustering and halo bias in the cosmic density field. In agreement with previous studies we find a strong increase of the correlation length with the mean inter-group separation. Although well determined observationally, we show that current theoretical predictions are not yet accurate enough to allow for stringent constraints on cosmological parameters. Finally, we use our results to explore the power-law nature of the 2PCF of galaxies. We split the 2PCF in 1-group and 2-group terms, equivalent to the 1-halo and 2-halo terms in halo occupation models, and show that the power-law form of the 2PCF is broken, when only including galaxies in the more massive systems.

Because of mass assignment onto grid points in the measurement of the power spectrum using the Fast Fourier Transform (FFT), the raw power spectrum $\la |\delta^f(k)|^2\ra$ estimated with FFT is not the same as the true power spectrum $P(k)$. In this paper, we derive the formula which relates $\la |\delta^f(k)|^2\ra$ to $P(k)$. For a sample of $N$ discrete objects, the formula reads: $\la |\delta^f(k)|^2\ra=\sum_{\vec n} [|W(\kalias)|^2P(\kalias)+1/N|W(\kalias)|^2]$, where $W(\vec k)$ is the Fourier transform of the mass assignment function $W(\vec r)$, $k_N$ is the Nyquist wavenumber, and $\vec n$ is an integer vector. The formula is different from that in some of previous works where the summation over $\vec n$ is neglected. For the NGP, CIC and TSC assignment functions, we show that the shot noise term $\sum_{\vec n} 1/N|W(\kalias)|^2]$ can be expressed by simple analytical functions. To reconstruct $P(k)$ from the alias sum $\sum_{\vec n}|W(\kalias)|^2 P(\kalias)$, we propose an iterative method. We test the method by applying it to an N-body simulation sample, and show that the method can successfully recover $P(k)$. The discussion is further generalized to samples with observational selection effects.

We use the halo occupation model to calibrate galaxy group finders in magnitude-limited redshift surveys. Since, according to the current scenario of structure formation, galaxy groups are associated with cold dark matter halos, we make use of the properties of the halo population in the design of our group finder. The method starts with an assumed mass-to-light ratio to assign a tentative mass to each group. This mass is used to estimate the size and velocity dispersion of the underlying halo that hosts the group, which in turn is used to determine group membership (in redshift space). This procedure is repeated until no further changes occur in group memberships. We find that the final groups selected this way are insensitive to the mass-to-light ratio assumed. We use mock catalogues, constructed using the conditional luminosity function (CLF), to test the performance of our group finder in terms of completeness of true members and contamination by interlopers. Our group finder is more successful than the conventional Friends-of-Friends group finder in assigning galaxies in common dark matter halos to a single group. We apply our group finder to the 2-degree Field Galaxy Redshift Survey and compare the resulting group properties with model predictions based on the CLF. For the $\Lambda$CDM `concordance' cosmology we find a clear discrepancy between the model and data in the sense that the model predicts too many rich groups. In order to match the observational results we have to either increase the mass-to-light ratios of rich clusters to a level significantly higher than current observational estimates, or to assume $\sigma_8 \simeq 0.7$, compared to the `concordance' value of 0.9.

We analyze the Minkowski functionals with a large $N$-body simulation of a standard $\Lambda$CDM model, focusing on transition scales between linear and non-linear gravitational evolution. We numerically calculate the Minkowski functionals with sufficient accuracies to investigate the transition scales, 10--$50\himpc$. The results are compared with analytic formulae of linear and second-order perturbation theories. We first show that the skewness parameters of the density fluctuations, which are important in second-order analytic formulae of the Minkowski functionals, are in good agreement with the perturbation theory. Considering relative differences between the Minkowski functionals of the analytic formulae and that of the simulation data, we evaluate accuracy levels of the predictions of the perturbation theory. When the straightforward threshold $\nu$ by density value is used in Minkowski functionals, the accuracy of the second-order perturbation theory is within 10% for smoothing length $R > 15\himpc$, and within a several % for $R > 20\himpc$. The accuracies of the linear theory are 2--5 times worse than that. When the rescaled threshold by volume fraction, $\nu_{\rm f}$ is used, accuracies of both linear and second-order theories are within a few % on all scales of $10\himpc < R < 50\himpc$.

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.

A basic assumption in current halo occupation model is that the properties of a galaxy depend only on the mass of its dark matter halo. An important consequence of this is that the segregation of the galaxy population by large-scale environment is entirely due to the environmental dependence of the halo population. In this paper we use such a model to predict how the galaxy luminosity function depends on large-scale environment. The latter is represented by the density contrast (delta) averaged over a spherical volume of radius R=8Mpc/h. The model predicts that the Schechter function is a good approximation to the luminosity functions of galaxies brighter than ~10^9 h^{-2}L_sun (b_j-band) in virtually all environments. The characteristic luminosity, L^{\star}, increases moderately with \delta. The faint-end slope, \alpha, on the other hand, is quite independent of $\delta$. However, when splitting the galaxy population into early and late types, it is found that for late-types \alpha is virtually constant, whereas for early-types \alpha increases from ~-0.3 in underdense regions (\delta~-0.5) to ~-0.8 in highly overdense regions with \delta ~ 10. The luminosity function at L_{b_j}<10^9 h^{-2} L_sun is significantly steeper than the extrapolation of the Schechter function that fits the brighter galaxies. This steepening is more significant for early-types and in low-density environments. The model also predicts that the luminosity density and mass density are closely correlated. The relation between the two is monotonic but highly non-linear. This suggests that one can use the luminosity density, averaged over a large volume, to rank the mass density. This, in turn, allows the environmental effects predicted here to be tested by observations.

In two recent papers, we developed a powerful technique to link the distribution of galaxies to that of dark matter haloes by considering halo occupation numbers as function of galaxy luminosity and type. In this paper we use these distribution functions to populate dark matter haloes in high-resolution N-body simulations of the standard LCDM cosmogony with Omega_m=0.3 and sigma_8=0.9. Stacking simulation boxes of 100 Mpc/h and 300 Mpc/h with 512^3 CDM particles each we construct Mock Galaxy Redshift Surveys out to a redshift of z=0.2. We use these mock surveys to investigate various clustering statistics. The projected correlation functions for galaxies with different luminosities and types match the observations well on scales larger than about 3 Mpc/h. On smaller scales, however, the model overpredicts the clustering power by about a factor two. Furthermore, the model predicts pairwise velocity dispersions (PVD) that are about 400 km/s too high at projected pair separations of ~1 Mpc/h. A strong velocity bias in massive haloes can reduce the predicted PVD to the observed level, but does not help to resolve the over-prediction of clustering power on small scales. Consistent results can be obtained within the standard LCDM model only when the mass-to-light ratio of clusters is almost twice the observed value. Alternatively, a LCDM model with sigma_8=0.75 can also reproduce the observational results.

We review the evidence for substructures from the anomalous flux ratios in gravitational lenses. Using high-resolution numerical simulations, we show that at typical image positions, the fraction of surface mass density in substructures is <~ 0.5% with mass above 10^{-4} virial masses in the ``concordance'' LCDM cosmology. Substructures outside the virial radius (but projected at typical lens image positions) only increase the fraction moderately. Several effects, in particular baryonic settling and the requirement of compactness, may further decrease the predictions by a factor of few. The predicted fraction with appropriate properties thus appears to be lower than that required by lensing, although both are still uncertain. More speculative substructures such as massive black holes (M ~ 10^5-10^6 solar masses) in the halo may offer viable alternatives.

We present a detailed non-spherical modeling of dark matter halos on the basis of a careful analysis of state-of-the art N-body simulations. The fitting formulae presented here form a complete and accurate description of the triaxial density profiles of halos in Cold Dark Matter(CDM) models. This modeling allows us to quantitatively discuss implications for shape observations of galactic and cluster halos. The predictions of the concordance $\lambda$CDM model are confronted with the shape observations from the Milky Way to X-ray clusters.

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.

We study the angular momentum profiles both for dark matter and for gas within virialized halos, using a statistical sample of halos drawn from cosmological hydrodynamics simulations. Three simulations have been analyzed, one is the ``non-radiative'' simulation, and the other two have radiative cooling. We find that the gas component on average has a larger spin and contains a smaller fraction of mass with negative angular momentum than its dark matter counterpart in the non-radiative model. As to the cooling models, the gas component shares approximately the same spin parameter as its dark matter counterpart, but the hot gas has a higher spin and is more aligned in angular momentum than dark matter, while the opposite holds for the cold gas. After the mass of negative angular momentum is excluded, the angular momentum profile of the hot gas component approximately follows the universal function originally proposed by Bullock et al. for dark matter, though the shape parameter $\mu$ is much larger for hot gas and is comfortably in the range required by observations of disk galaxies. Since disk formation is related to the distribution of hot gas that will cool, our study may explain the fact that the disk component of observed galaxies contains a smaller fraction of low angular momentum material than dark matter in halos.

This paper uses numerical simulations to test the formation time distribution of dark matter haloes predicted by the analytic excursion set approaches. The formation time distribution is closely linked to the conditional mass function and this test is therefore an indirect probe of this distribution. The excursion set models tested are the extended Press-Schechter (EPS) model, the ellipsoidal collapse (EC) model, and the non-spherical collapse boundary (NCB) model. Three sets of simulations (6 realizations) have been used to investigate the halo formation time distribution for halo masses ranging from dwarf-galaxy like haloes ($M=10^{-3} M_*$, where $M_*$ is the characteristic non-linear mass scale) to massive haloes of $M=8.7 M_*$. None of the models can match the simulation results at both high and low redshift. In particular, dark matter haloes formed generally earlier in our simulations than predicted by the EPS model. This discrepancy might help explain why semi-analytic models of galaxy formation, based on EPS merger trees, under-predict the number of high redshift galaxies compared with recent observations.

We present velocity statistics of galaxies and their biases inferred from the statistics of the underlying dark matter using a cosmological hydrodynamic simulation of galaxy formation in low-density and spatially flat cold dark matter cosmogony. We find that the pairwise velocity dispersion (PVD) of all galaxies is significantly lower than that of the dark matter particles, and that the PVD of the young galaxies is lower than that of the old types, and even of all galaxies together, especially at small separations. These results are in reasonable agreement with the recent measurements of PVDs in the Las Campanas redshift survey, the PSCz catalogue and the SDSS data. We also find that the low PVD of young galaxies is due to the effects of dynamical friction as well as the different spatial distribution. We also consider the mean infall velocity and the POTENT density reconstruction that are often used to measure the cosmological parameters, and investigate the effects of spatial bias and dynamical friction. In our simulation, the mean infall velocity of young galaxies is significantly lower than that of all the galaxies or of the old galaxies, and the dynamical bias becomes important on scales less than 3Mpc/h. The mass density field reconstructed from the velocity field of young galaxies using the POTENT-style method suffers in accuracy both from the spatial bias and the dynamical friction on the smoothing scale of R_s=8Mpc/h. On the other hand, in the case of R_s=12Mpc/h, which is typically adopted in the actual POTENT analysis, the density reconstruction based on various tracers of galaxies is reasonably accurate.

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 present a detailed non-spherical modeling of dark matter halos on the basis of a combined analysis of the high-resolution halo simulations (12 halos with $N\sim 10^6$ particles within their virial radius) and the large cosmological simulations (5 realizations with $N=512^3$ particles in a $100h^{-1}$Mpc boxsize). The density profiles of those simulated halos are well approximated by a sequence of the concentric triaxial distribution with their axis directions being fairly aligned. We characterize the triaxial model quantitatively by generalizing the universal density profile which has previously been discussed only in the framework of the spherical model. We obtain a series of practically useful fitting formulae in applying the triaxial model; the mass and redshift dependence of the axis ratio, the mean of the concentration parameter, and the probability distribution functions of the the axis ratio and the concentration parameter. These accurate fitting formulae form a complete description of the triaxial density profiles of halos in Cold Dark Matter models. Our current description of the dark halos will be particularly useful in predicting a variety of nonsphericity effects, to a reasonably reliable degree, including the weak and strong lens statistics, the orbital evolution of galactic satellites and triaxiality of galactic halos, and the non-linear clustering of dark matter. In addition, this provides a useful framework for the non-spherical modeling of the intra-cluster gas, which is crucial in discussing the gas and temperature profiles of X-ray clusters and the Hubble constant estimated via the Sunyaev -- Zel'dovich effect.

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.

We use N-body simulations to test the predictions of the redshift distortion in the power spectrum given by the halo model in which the clustering of dark matter particles is considered as a result both of the clustering of dark halos in space and of the distribution of dark matter particles in individual dark halo. The predicted redshift distortion depends sensitively on several model parameters in a way different from the real-space power spectrum. An accurate model of the redshift distortion can be constructed if the following properties of the halo population are modelled accurately: the mass function of dark halos, the velocity dispersion among dark halos, and the non-linear nature of halo bias on small scales. The model can be readily applied to interpreting the clustering properties and velocity dispersion of different populations of galaxies once a cluster-weighted bias (or equivalently an halo occupation number model) is specified for the galaxies. Some non-trivial bias features observed from redshift surveys of optical galaxies and of IRAS galaxies relative to the standard low-density cold dark matter model can be easily explained in the cluster weighted bias model. The halo model further indicates that a linear bias can be a good approximation only on for k <= 0.1 hMpc^{-1}.

We use a large set of state-of-the-art cosmological N-body simulations [512^3 particles] to study the intrinsic ellipticity correlation functions of halos. With the simulations of different resolutions, we find that the ellipticity correlations converge once the halos have more than 160 members. For halos with fewer members, the correlations are underestimated, and the underestimation amounts to a factor of 2 when the halos have only 20 particles. After correcting for the resolution effects, we show that the ellipticity correlations of halos in the bigger box (L=300 mpc) agree very well with those obtained in the smaller box (L=100 mpc). Combining these results from the different simulation boxes, we present accurate fitting formulae for the ellipticity correlation function c_{11}(r) and for the projected correlation functions Sigma_{11}(r_p) and Sigma_{22}(r_p) over three orders of magnitude in halo mass. The latter two functions are useful for predicting the contribution of the intrinsic correlations to deep lensing surveys. With reasonable assumptions for the redshift distribution of galaxies and for the mass of galaxies, we find that the intrinsic ellipticity correlation can contribute significantly not only to shallow surveys but also to deep surveys. Our results indicate that previous similar studies significantly underestimated this contribution for their limited simulation resolutions.

We study the angular momentum profile of dark matter halos for a statistical sample drawn from a set of high-resolution cosmological simulations of $256^3$ particles. Two typical Cold Dark Matter (CDM) models have been analyzed, and the halos are selected to have at least $3\times 10^4$ particles in order to reliably measure the angular momentum profile. In contrast with the recent claims of Bullock et al., we find that the degree of misalignment of angular momentum within a halo is very high. About 50 percent of halos have more than 10 percent of halo mass in the mass of negative angular momentum $j$. After the mass of negative $j$ is excluded, the cumulative mass function $M(<j)$ follows approximately the universal function proposed by Bullock et al., though we still find a significant fraction of halos ($\sim 50%$) which exhibit systematic deviations from the universal function. Our results, however, are broadly in good agreement with a recent work of van den Bosch et al.. We also study the angular momentum profile of halos in a Warm Dark Matter (WDM) model and a Self-Interacting Dark Matter (SIDM) model. We find that the angular momentum profile of halos in the WDM is statistically indistinguishable from that in the CDM model, but the angular momentum of halos in the SIDM is reduced by the self-interaction of dark matter.

Understanding of the structure evolution in the Universe has been greatly advanced with a rapid progress in high-resolution cosmological simulations. In this contribution, we report a new set of cosmological simulations with $512^3$ ($\sim 1.4\times 10^8$) particles in a simulation box of $100\mpc$. With this unprecedented resolution, we have successfully overcomed the over-merging problem while retaining a cosmological volume size that is crucial for the statistical analysis. The simulations are being applied to studying important astrophysical problems that are sensitive to the simulation resolution.

Cold Dark Matter (CDM) models of galaxy formation had been remarkably successful to explain a number of observations in the past decade. However, with both the theoretical modeling and the observations being improved, CDM models have been very recently shown to have excessive clustering on the sub-galactic scale. Here I discuss a solution, based on our high-resolution numerical simulations, to this outstanding problem by considering Warm Dark Matter (WDM). Our results show that the over-clustering problem on sub-galactic scales can be overcome by WDM models, and all the advantages of CDM models are preserved by WDM models. Therefore, the WDM model will become an interesting alternative to the well-studied CDM models.

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.

The fundamental plane of spiral galaxies was searched from observational data, which can be approximately represented in observational parameters by $ L \propto V^2 R $, where $L$, $R$, and $V$ are the luminosity, the linear size of galactic disk, and the rotation velocity. This plane exists at all optical bands in our samples of more than 500 spiral galaxies in total. It is more fundamental than the relationship between any two of the parameters, and can reduce the residual of the Tully--Fisher $L$--$V$ relations by about 50\%. Noticed that the power index of $V$ is doubled from that of $R$, which implies that the total mass and the mass distribution of galaxies plays an important role in forming the fundamental plane. Involving the third parameter of galactic size has a strong physical implication on galaxy formation and the dark--mass distribution.

We perform an extensive analysis of nonlinear and stochastic biasing of galaxies and dark halos in spatially flat low-density CDM universe using cosmological hydrodynamic simulations. We compare their biasing properties with the predictions of an analytic halo biasing model. Dark halos in our simulations exhibit reasonable agreement with the predictions only on scales larger than 10h^{-1}Mpc, and on smaller scales the volume exclusion effect of halos due to their finite size becomes substantial. Interestingly the biasing properties of galaxies are better described by extrapolating the halo biasing model predictions. We also find the clear dependence of galaxy biasing on their formation epoch; the distribution of old populations of galaxies tightly correlates with the underlying mass density field, while that of young populations is slightly more stochastic and anti-biased relative to dark matter. The amplitude of two-point correlation function of old populations is about 3 times larger than that of the young populations. Furthermore, the old population of galaxies reside within massive dark halos while the young galaxies are preferentially formed in smaller dark halos. Assuming that the observed early and late-type galaxies correspond to the simulated old and young populations of galaxies, respectively, all of these segregations of galaxies are consistent with observational ones for the early and late-type of galaxies such as the morphology--density relation of galaxies.

We apply the morphological descriptions of two-dimensional contour map, the so-called Minkowski functionals (the area fraction, circumference, and Euler characteristics), to the convergence field $\kappa(\bm{\theta})$ of the large-scale structure reconstructed from the shear map produced by the ray-tracing simulations. The perturbation theory of structure formation has suggested that the non-Gaussian features on the Minkowski functionals with respect to the threshold in the weakly nonlinear regime are induced by the three skewness parameters of $\kappa$ that are sensitive to the density parameter of matter, $\Omega_{\rm m}$. We show that, in the absence of noise due to the intrinsic ellipticities of source galaxies with which the perturbation theory results can be recovered, the accuracy of $\Omega_{\rm m}$ determination is improved by $\sim 20%$ using the Minkowski functionals compared to the conventional method of using the direct measure of skewness.

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

The stable clustering hypothesis is a fundamental assumption about the nonlinear clustering of matter in cosmology. It states that the mean physical separation of particles is a constant on sufficiently small scales. While many authors have attempted to test the hypothesis with cosmological N-body simulations, no consensus has been reached on whether and where the hypothesis is valid, because of the limited dynamical range this type of simulations can achieve. In this Letter, we propose to test the hypothesis with high resolution halo simulations, since the individual halo simulations can resolve much better the fine structures of the halos and since almost all pairs of particles with small separations are presumed to be inside virialized halos. We calculated the mean pair velocity for 14 high resolution halos of $\sim 1$ million particles in a low-density flat cold dark matter model. The result agrees very well with the stable clustering prediction within the measurement uncertainty $\sim 30%$ over a large range of scales where the overdensity is $10^3$ to $10^6$. The accuracy of the test can be improved to $\sim 10%$ if some 100 halos with a similar resolution are analyzed.

We quantify the degree of nonlinearity and stochasticity of the clustering of biased objects, using cosmological N-body simulations. Adopting the peaks and the halos as representative biasing models, we focus on the two-point correlation of the biased objects, dark matter and their cross-correlation. Especially, we take account of the effect of redshift-space distortion and attempt to clarify the scale-dependence and the time-dependence by analyzing the biasing factor and the cross-correlation factor. On small scales, stochasticity and nonlinearity become appreciable and strongly object-dependent, especially in redshift space due to the pair-wise velocity dispersion of the biased objects. Nevertheless, an approximation of deterministic linear biasing $\delta_{\rm obj}\simeq b_{\rm obj} \delta_{\rm mass}$ works reasonably well even in the quasi-linear regime $r > 10 h^{-1}$ Mpc, and linear redshift-space distortion explains the clustering amplitudes in redshift space in this regime.

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.

We report a series of high-resolution N-body simulations designed to examine the internal physical properties of dark matter halos. A total of fifteen halos, each represented by $\sim 1$ million particles within the virial radius, have been simulated covering the mass range of $2\times10^{12}\sim 5\times10^{14} h^{-1}{\rm M_\odot}$. As the first application of these simulations, we have examined the density profiles of the halos. We found a clear systematic correlation between the halo mass and the slope of the density profile at one percent of the virial radius, in addition to the variations of the slope among halos of the similar mass. More specifically, the slope is $\sim - 1.5$, -1.3, and -1.1 for galaxy, group, and cluster mass halos, respectively. Thus we conclude that the dark matter density profiles, especially in the inner region, are not universal.

We study the diversity of the density profiles of dark matter halos based on a large set of high-resolution cosmological simulations of 256^3 particles. The cosmological models include four scale-free models and three representative cold dark matter models. The simulations have good force resolution, and there are about 400 massive halos with more than 10^4 particles within the virial radius in each cosmological model. Our unbiased selection of all massive halos enables to quantify how well the bulk of dark matter halos can be described by the Navarro, Frenk & White (NFW) profile which was established for equilibrium halos. We find that about seventy percent of the halos can be fitted by the NFW profile with a fitting residual dvi_{max} less than 30% in Omega_0=1 universes. This percentage is higher in lower density cosmological models. The rest of the halos exhibits larger deviations from the NFW profile for more significant internal substructures. There is a considerable amount of variation in the density profile even for the halos which can be fitted by the NFW profile (i.e. dvi_{max}<0.30). The distribution of the profile parameter, the concentration $c$, can be well described by a lognormal function with the mean value \bar c slightly smaller (15%) than the NFW result and the dispersion \sigma_c in \ln c about 0.25. The more virialized halos with dvi_{max}<0.15 have the mean value \bar c in good agreement with the NFW result and a slightly smaller dispersion \sigma_c (about 0.2). Our results can alleviate some of the conflicts found recently between the theoretical NFW profile and observational results. Implications for theoretical and observational studies of galaxy formation are discussed.

We present results from a series of cosmological SPH (smoothed particle hydrodynamics) simulations coupled with the P3M (Particle-Particle-Particle-Mesh) solver for the gravitational force. The simulations are designed to predict the statistical properties of X-ray clusters of galaxies as well as to study the formation of galaxies. We have seven simulation runs with different assumptions on the thermal state of the intracluster gas. Following the recent work by Pearce et al., we modify our SPH algorithm so as to phenomenologically incorporate the galaxy formation by decoupling the cooled gas particles from the hot gas particles. All the simulations employ 128^3 particles both for dark matter and for gas components, and thus constitute the largest systematic catalogues of simulated clusters in the SPH method performed so far. These enable us to compare the analytical predictions on statistical properties of X-ray clusters against our direct simulation results in an unbiased manner. We find that the luminosities of the simulated clusters are quite sensitive to the thermal history and also to the numerical resolution of the simulations, and thus are not reliable. On the other hand, the mass-temperature relation for the simulated clusters is fairly insensitive to the assumptions of the thermal state of the intracluster gas, robust against the numerical resolution, and in fact agrees well with the analytic prediction. Therefore the prediction for the X-ray temperature function of clusters on the basis of the Press-Schechter mass function and the virial equilibrium is fairly reliable.

We perform a series of high-resolution N-body simulations designed to examine the density profiles of dark matter halos. From 12 simulated halos ranging the mass of $2\times10^{12}\sim 5\times10^{14} h^{-1}{\rm M_\odot}$ (represented by $\sim 1$ million particles within the virial radius), we find a clear systematic correlation between the halo mass and the slope of the density profile at 1% of the virial radius, in addition to the variations of the slope among halos of the similar mass. More specifically, the slope is $\sim -1.5$, -1.3, and -1.1 for galaxy, group, and cluster mass halos, respectively. While we confirm the earlier simulation results that the inner slope is steeper than the {\it universal} profile originally proposed by Navarro, Frenk & White, this mass dependence is inconsistent with the several analytical arguments attempting to link the inner slope with the primordial index of the fluctuation spectrum. Thus we conclude that the dark matter density profiles, especially in the inner region, are not universal.

We have developed a model to describe two-point correlation functions of clusters of galaxies in X-ray flux-limited surveys. Our model properly takes account of nonlinear gravitational evolution of mass fluctuations, redshift-space distortion due to linear peculiar velocity field and to finger-of-god, cluster abundance and bias evolution on the basis of the Press -- Schechter theory, the light-cone effect, and the selection function due to the X-ray flux, temperature and luminosity limits. Applying this model in representative cosmological models, we have presented quantitative predictions for X-ray selected samples feasible from the future surveys with the X-ray satellites including Astro-E, Chandra, and XMM. The comparison of these predictions and the observed cluster clustering will place important cosmological constraints which are complementary to the cluster abundance and the cosmic microwave background.

We examine the cosmological redshift-space distortion effect on the power spectrum of the objects at high-redshifts, which is an unavoidable observational contamination in general relativistic cosmology. In particular, we consider the nonlinear effects of density and velocity evolution using high-resolution N-body simulations in cold dark matter models. We find that the theoretical modeling on the basis of the fitting formulae of nonlinear density and velocity fields accurately describes the numerical results, especially in quasi-nonlinear regimes. These corrections for nonlinear effects are essential in order to use the the cosmological redshift-space distortion as a cosmological test. We perform a feasibility test to derive constraints from the future catalogues of high-redshift quasars using our theoretical modeling and results of N-body simulations. Applying the present methodology to the future data from on-going surveys of high-redshift galaxies and quasars will provide a useful tool to constrain a geometry of the universe.

The three-dimensional distribution of astronomical objects observed in redshift space significantly differs from the true distribution since the distance to each object cannot be determined by its redshift $z$ only; for $z \ll 1$ the peculiar velocity field contaminates the true recession velocity of the Hubble flow, while the true distance for objects at $z < 1$ sensitively depends on the (unknown and thus assumed) cosmological parameters. This hampers the effort to understand the true distribution of large-scale structure of the universe. In addition, all cosmological observations are carried out on a light-cone, the null hypersurface of an observer at $z=0$. This implies that their intrinsic properties and clustering statistics should change even within the survey volume. Therefore a proper comparison taking account of the light-cone effect is important to extract any cosmological information from redshift catalogues, especially for $z < 1$. We present recent theoretical development on the two effects -- the cosmological light-cone effect and the cosmological redshift-space distortion -- which should play key roles in observational cosmology in the 21st century.

We use a new method, the cross power spectrum between the linear density field and the halo number density field, to measure the Lagrangian bias for dark matter halos. The method has several important advantages over the conventional correlation function analysis. By applying this method to a set of high-resolution simulations of 256^3 particles, we have accurately determined the Lagrangian bias, over 4 magnitudes in halo mass, for four scale-free models with the index n=-0.5, -1.0, -1.5 and -2.0 and three typical CDM models. Our result for massive halos with $M \ge M_*$ ($M_*$ is a characteristic non-linear mass) is in very good agreement with the analytical formula of Mo & White for the Lagrangian bias, but the analytical formula significantly underestimates the Lagrangian clustering for the less massive halos $M < M_*. Our simulation result however can be satisfactorily described, with an accuracy better than 15%, by the fitting formula of Jing for Eulerian bias under the assumption that the Lagrangian clustering and the Eulerian clustering are related with a linear mapping. It implies that it is the failure of the Press-Schechter theories for describing the formation of small halos that leads to the inaccuracy of the Mo & White formula for the Eulerian bias. The non-linear mapping between the Lagrangian clustering and the Eulerian clustering, which was speculated as another possible cause for the inaccuracy of the Mo & White formula, must at most have a second-order effect. Our result indicates that the halo formation model adopted by the Press-Schechter theories must be improved.

The density profile and the clustering of dark matter (DM) halos are two very important ingredients for understanding many observations of galaxies. In this talk, I present our new results for these two quantities from a large set of high quality N-body simulations. We show that in an Einstein de Sitter Universe about 35 percent of the halos could not be fitted by the Navarro, Frenk & White (NFW) profile because of too significant substructures, another 50 percent with less substructures can be reasonably described by the NFW profile, and the rest with the least substructures can be fitted by the NFW profile very nicely. For the halos which can be reasonably or better fitted by the NFW profile, the probability distribution function of the concentration parameter $c$ could be well described by a lognormal function with a dispersion $\sigma=0.27$. The mean value of $c$ is quite close to that predicted by NFW, though it decreases with the increase of the substructures. Our clustering analysis shows a significantly stronger clustering for low mass halos than Mo & White predict. A formula is therefore presented for the two-point correlation function of halos which can fit our simulation results (both CDM models and scale-free clustering models) very accurately.

An accurate fitting formula is reported for the two-point correlation function of the dark matter halos in hierarchical clustering models, It is valid for the the linearly clustering regime, and its accuracy is about 10% in the correlation amplitude for the halos with mass M greater than 1/100 -- 1/1000 of the characteristic non-linear mass M_c. The result is found on the basis of a careful analysis for a large set of scale-free simulations with 17 million particles. The fitting formula has a weak explicit dependence on the index n of the initial power spectrum, but can be equally well applied to the cold dark matter (CDM) cosmological models if the effective index n_{eff} of the CDM power spectrum at the scale of the halo mass replaces the index n. The formula agrees with the analytical formula of Mo & White (MW96) for massive halos with M>M_c, but the MW96 formula significantly underpredicts the correlation for the less massive halos. The difference between the fitting and the analytical formulae amounts to a factor => 2 in the correlation amplitude for M=0.01 M_c. One of the most interesting applications of this fitting formula would be the clustering of galaxies since the majority of halos hosting galaxies satisfies M<< M_c.

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.

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 perform a detailed analysis of the statistical significance of a concentration of Lyman break galaxies at $z \sim 3$ recently discovered by Steidel et al. (1997), using a series of N-body simulations with $N=256^3$ particles in a $(100\himpc)^3$ comoving box. While the observed number density of Lyman break galaxies at $z\sim3$ implies that they correspond to systems with dark matter halos of $\simlt 10^{12}M_\odot$, the resulting clustering of such objects on average is not strong enough to be reconciled with the concentration if it is fairly common; we predict one similar concentration approximately per ($6\sim 10$) fields in three representative cold dark matter models. Considering the current observational uncertainty of the frequency of such clustering at $z\sim3$, it would be premature to rule out the models, but the future spectroscopic surveys in a dozen fields could definitely challenge all the existing cosmological models a posteriori fitted to the $z=0$ universe.

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.

In this talk, I will show how to determine the biasing factor $b$ from the high-order moments of galaxies. The determination is based on the analytical modeling of primordial peaks and virialized halos and is independent of the currently unknown density parameter $\Omega_0$ and other cosmological parameters. The observed high-oder moments of the APM galaxies require that the biasing factor $b$ be very close to 1, i.e. the optical galaxies are an unbiased tracer of the underlying mass distribution (on quasilinear scale). The theoretical argument can be easily generalized to the three-point correlation function and the bispectrum both of which can used as further observational tests to the important conclusion of $b\approx 1$ drawn from the high-order moments. Finally I present our preliminary results of the three-point correlation functions for the Las Campanas Redshift Survey.

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.

We show that in a hierarchical clustering model the low-order statistics of the density and the peculiar velocity fields can all be modelled semianalytically for a given cosmology and an initial density perturbation power spectrum $P(k)$. We present such models for the two-point correlation function $\xi(r)$, the amplitude $Q$ of the three-point correlation function, the mean pairwise peculiar velocity $< v_{12}(r)> $, the pairwise peculiar velocity dispersion $< v_{12}^2(r)>$, and the one-point peculiar velocity dispersion $< v_1^2 >$. We test our models against results derived from N-body simulations. These models allow us to understand in detail how these statistics depend on $P(k)$ and cosmological parameters. They can also help to interpret, and maybe correct for, sampling effects when these statistics are estimated from observations. The dependence of the small-scale pairwise peculiar velocity dispersion on rich clusters in the sample, for instance, can be studied quantitatively. There are also significant implications for the reconstruction of the cosmic density field from measurements in redshift space.

We discuss the effect of the finite size of galaxies on estimating small-scale relative pairwise peculiar velocity dispersions from the cosmic virial theorem (CVT). Specifically we evaluate the effect by incorporating the finite core radius $r_c$ in the two-point correlation function of mass, i.e. $\xi_\rho(r) \propto (r+r_c)^{-\gamma}$ and the effective gravitational force softening $r_s$ on small scales. We analytically obtain the lowest-order correction term for $\gamma <2$ which is in quantitative agreement with the full numerical evaluation. With a nonzero $r_s$ and/or $r_c$ the cosmic virial theorem is no longer limited to the case of $\gamma<2$. We present accurate fitting formulae for the CVT predicted pairwise velocity dispersion for the case of $\gamma>2$. Compared with the idealistic point-mass approximation ($r_s=r_c=0$), the finite size effect can significantly reduce the small-scale velocity dispersions of galaxies at scales much larger than $r_s$ and $r_c$. Even without considering the finite size of galaxies, nonzero values for $r_c$ are generally expected, for instance, for cold dark matter (CDM) models with a scale-invariant primordial spectrum. For these CDM models, a reasonable force softening $r_s\le 100 \hikpc$ would have rather tiny effect. We present the CVT predictions for the small-scale pairwise velocity dispersion in the CDM models normalized by the COBE observation. The implication of our results for confrontation of observations of galaxy pair-wise velocity dispersions and theoretical predictions of the CVT is also discussed.

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

We develop an analytic model for the hierarchical correlation amplitudes S_j(R)= \bxi_j(R)/\bxi_2^{j-1}(R) of density peaks and dark matter halos in the quasi-linear regime. The statistical distribution of density peaks and dark matter halos within the initial density field are determined by the peak formalism and by an extension of the Press-Schechter formalism, respectively. Modifications of these distributions caused by gravitationally induced motions are treated using a spherical collapse model. We test our model against results from a variety of N-body simulations. The model works well for peaks and for halos that are identified earlier than the time when the moments are calculated. Because halos are spatially exclusive at the time of their identification, our model is only qualitatively correct for halos identified at the same time as the moments are calculated. The S_j depend only weakly on the bias parameter b for large b but increase rapidly with decreasing b at b\sim 1. Thus if galaxies are associated with peaks in the initial density field, or with dark halos formed at high redshifts, a measurement of S_j in the quasilinear regime should determine whether galaxies are significantly biased relative to the mass. We use our model to interpret the observed high order correlation functions of galaxies and clusters. We find that if the values of S_j for galaxies are as high as those given by the APM survey, then APM galaxies should not be significantly biased.

We show that if a sample of galaxy clusters is complete above some mass threshold, then hierarchical clustering theories for structure formation predict its autocorrelation function to be determined purely by the cluster abundance and by the spectrum of linear density fluctuations. Thus if the shape of the initial fluctuation spectrum is known, its amplitude $\sigma_8$ can be estimated directly from the correlation length of a cluster sample in a way which is independent of the value of $\Omega_0$. If the cluster mass corresponding to the sample threshold is also known, it provides an independent estimate of the quantity $\sigma_8\Omega_0^{0.6}$. Thus cluster data should allow both $\sigma_8$ and $\Omega_0$ to be determined observationally. We explore these questions using N-body simulations together with a simple but accurate analytical model based on extensions of Press-Schechter theory. Applying our results to currently available data we find that if the linear fluctuation spectrum has a shape similar to that suggested by the APM galaxy survey, then a correlation length $r_0$ in excess of $20\mpch$ for Abell clusters would require $\sigma_8>1$, while $r_0<15\mpch$ would require $\sigma_8<0.5$. With conventional estimates of the relevant mass threshold these imply $\Omega_0\la 0.3$ and $\Omega_0\ga 1$ respectively.

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.

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

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.