We use observational measurements of thermal and kinetic Sunyaev-Zel'dovich effects, as well as soft X-ray emission associated with galaxy groups to constrain the gas density and temperature in the circumgalactic medium (CGM) for dark matter halos with masses above $10^{12.5}M_{\odot}$. A number of generic models are used together with a Bayesian scheme to make model inferences. We find that gas with a single temperature component cannot provide a consistent model to match the observational data. A simple two-phase model assuming a hot component and an ionized warm component can accommodate all the three observations. The total amount of the gas in individual halos is inferred to be comparable to the universal baryon fraction corresponding to the halo mass. The inferred temperature of the hot component is comparable to the halo virial temperature. The fraction of the hot component increases from $(15-40)\%$ for $10^{12.5}{M}_\odot$ halos to $(40-60)\%$ for $10^{14.5}{M}_\odot$ halos, where the ranges reflect uncertainties in the assumed gas density profile. Our results suggest that a significant fraction of the halo gas is in a non-thermalized component with temperature much lower than the virial temperature.

Identifying galaxy groups from redshift surveys of galaxies plays an important role in connecting galaxies with the underlying dark matter distribution. Current and future high-$z$ spectroscopic surveys, usually incomplete in redshift sampling, present both opportunities and challenges to identifying groups in the high-$z$ Universe. We develop a group finder that is based on incomplete redshift samples combined with photometric data, using a machine learning method to assign halo masses to identified groups. Test using realistic mock catalogs shows that $\gtrsim 90\%$ of true groups with halo masses $\rm M_h \gtrsim 10^{12} M_{\odot}/h$ are successfully identified, and that the fraction of contaminants is smaller than $10\%$. The standard deviation in the halo mass estimation is smaller than 0.25 dex at all masses. We apply our group finder to zCOSMOS-bright and describe basic properties of the group catalog obtained.

Wolf-Rayet (WR) galaxies are a rare population of galaxies that host living high-mass stars during their WR phase (i.e. WR stars) and are thus expected to provide interesting constraints on the stellar Initial Mass Function, massive star formation, stellar evolution models, etc. Spatially resolved spectroscopy should in principle provide a more efficient way of identifying WR galaxies than single-fiber surveys of galactic centers such as SDSS-I & II, as WR stars should be more preferentially found in discs. Using IFU data from the ongoing SDSS-IV MaNGA survey, we have performed a thorough search for WR galaxies. We first identify H II regions in each datacube and carry out full spectral fitting to the stacked spectra. We then visually inspect the residual spectrum of each H II region and identify WR regions that present a significant "blue bump" at 4600-4750 A. The resulting WR catalog includes 267 WR regions of ~500pc (radius) sizes, distributed in 90 galaxies from the current sample of MaNGA (MaNGA Product Launch 7). We find WR regions are exclusively found in galaxies that show bluest colors and highest star formation rates for their mass. Most WR galaxies have late-type morphologies and show relatively large asymmetry in their images, implying that WR regions are more preferentially found in interacting/merging galaxies. We estimate the stellar mass function of WR galaxies and the mass-dependent detection rate. The detection rate of WR galaxies is typically ~2%, with weak dependence on stellar mass. This detection rate is about 40 times higher than previous studies with SDSS single fiber data, and by a factor of 2 lower than the CALIFA-based WR catalog. We make comparisons with SDSS and CALIFA studies, and conclude that different detection rates can be explained mainly by three factors: spatial coverage, spectral signal-to-noise ratio, and redshift ranges of the parent sample.

We develop a method to estimate the dust attenuation curve of galaxies from full spectral fitting of their optical spectra. Motivated from previous studies, we separate the small-scale features from the large-scale spectral shape, by performing a moving average method to both the observed spectrum and the simple stellar population model spectra. The intrinsic dust-free model spectrum is then derived by fitting the observed ratio of the small-scale to large-scale (S/L) components with the S/L ratios of the SSP models. The selective dust attenuation curve is then determined by comparing the observed spectrum with the dust-free model spectrum. One important advantage of this method is that the estimated dust attenuation curve is independent of the shape of theoretical dust attenuation curves. We have done a series of tests on a set of mock spectra covering wide ranges of stellar age and metallicity. We show that our method is able to recover the input dust attenuation curve accurately, although the accuracy depends slightly on signal-to-noise ratio of the spectra. We have applied our method to a number of edge-on galaxies with obvious dust lanes from the ongoing MaNGA survey, deriving their dust attenuation curves and $E(B-V)$ maps, as well as dust-free images in $g$, $r$, and $i$ bands. These galaxies show obvious dust lane features in their original images, which largely disappear after we have corrected the effect of dust attenuation. The vertical brightness profiles of these galaxies become axis-symmetric and can well be fitted by a simple model proposed for the disk vertical structure. Comparing the estimated dust attenuation curve with the three commonly-adopted model curves, we find that the Calzetti curve provides the best description of the estimated curves for the inner region of galaxies, while the Milky Way and SMC curves work better for the outer region.

A clear transition of galaxy quenching is identified in the multi-parameter space of stellar mass ($M_*$), bulge to total mass ratio ($B/T_{\rm m}$), halo mass ($M_{\rm h}$) and halo-centric distance ($r/r_{180}$). For given halo mass, the characteristic stellar mass ($M_{*, \rm ch}$) for the transition is about one-fifth of that of the corresponding central galaxy, and almost independent of $B/T_{\rm m}$. Once $B/T_{\rm m}$ is fixed, the quenched fraction of galaxies with $M_*<M_{*, \rm ch}$ increases with $M_{\rm h}$, but decreases with $M_*$ in the inner part of halos ($r/r_{180}<0.5$). In the outer part ($r/r_{180}>0.5$), the trend with $M_{\rm h}$ remains but the correlation with $M_*$ is absent or becomes positive. For galaxies above $M_{\rm *, ch}$ and with $B/T_{\rm m}$ fixed, the quenched fraction increases with $M_{\rm *}$, but depends only weakly on $M_{\rm h}$ in both the inner and outer regions. Moreover, at fixed $B/T_{\rm m}$ and $M_*$, the quenched fraction increases with the decreases of $r/r_{180}$ for galaxies with $M_*<M_{*, \rm ch}$, while remaining constant with $r/r_{180}$ for galaxies with $M_*>M_{*, \rm ch}$. Our finding provides a physically-motivated way to classify galaxies in halos into two classes based on their quenching properties: an `upper class' with $M_*>M_{\rm *,ch}$ and a `lower class' with $M_*<M_{\rm *,ch}$. Environmental quenching is important only for `lower class' galaxies, while internal quenching plays the dominating role for the `upper class'.

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

Satellites constitute an important fraction of the overall galaxy population and are believed to form in dark matter subhalos. Here we use the cosmological hydrodynamic simulation TNG100 to investigate how the formation histories of subhalos affect the properties and evolution of their host galaxies. We use a scaled formation time ($a_{\rm nf}$) to characterize the mass assembly histories of the subhalos before they are accreted by massive host halos. We find that satellite galaxies in young subhalos (low $a_{\rm nf}$) are less massive and more gas rich, and have stronger star formation and a higher fraction of ex situ stellar mass than satellites in old subhalos (high $a_{\rm nf}$). Furthermore, these low $a_{\rm nf}$ satellites require longer timescales to be quenched as a population than the high $a_{\rm nf}$ counterparts. We find very different merger histories between satellites in fast accretion (FA, $a_{\rm nf}<1.3$) and slow accretion (SA, $a_{\rm nf}>1.3$) subhalos. For FA satellites, the galaxy merger frequency dramatically increases just after accretion, which enhances the star formation at accretion. While, for SA satellites, the mergers occur smoothly and continuously across the accretion time. Moreover, mergers with FA satellites happen mainly after accretion, while a contrary trend is found for SA satellites. Our results provide insight into the evolution and star formation quenching of the satellite population.

We measure the star formation histories (SFH) of a sample of low-mass galaxies with $M_\ast<10^9M_\odot$ from the SDSS-IV MaNGA survey. The large number of IFU spectra for each galaxy are either combined to reach a high signal to noise ratio or used to investigate spatial variations. We use Bayesian inferences based on full spectrum fitting. Our analysis based on Bayesian evidence ratio indicates a strong preference for a model that allows the presence of an old stellar population, and that an improper model for the SFH can significantly underestimate the old population in these galaxies. The addition of NIR photometry to the constraining data can further distinguish between different SFH model families and significantly tighten the constraints on the mass fraction in the old population. On average more than half of the stellar mass in present-day low-mass galaxies formed 8 Gyrs ago, while about 30\% within the past 4 Gyrs. Satellite galaxies on average have formed their stellar mass earlier than central galaxies, and stars in the outer regions of galaxy are younger than those in the central part. Our results suggest that most of the low-mass galaxies have an early episode of active star formation that produces a large fraction of their present stellar mass.

We use the ELUCID N-Body simulation to study the relation of structural properties of dark matter halos to their assembly history and environment. We find that the complexity of halo assembly histories can be well described by a small number of principal components. Using decision trees built with the random ensemble method, we find that the correlation between halo concentration and assembly history is tight: more than 60% of the variance in halo concentration can be explained by assembly history alone. Halo shape and spin parameters are found to depend not only on assembly history, but also on environment and initial condition, and these dependencies are not orthogonal. The structures of more massive halos depend more strongly on higher-order components of the assembly history. The relations between halo intrinsic properties and environment are weak compared to their variances, with the anisotropy of local tidal fields having the strongest correlation with halo properties. Our method of dimension reduction and regression can help simplify the characterization of the halo population and clarify the degeneracy among halo properties.

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.

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

We present the detection of the kinetic Sunyaev-Zel'dovich effect (kSZE) signals from groups of galaxies as a function of halo mass down to $\log (M_{500}/{\rm M_\odot}) \sim 12.3$, using the {\it Planck} CMB maps and stacking about $40,000$ galaxy systems with known positions, halo masses, and peculiar velocities. The signals from groups of different mass are constrained simultaneously to take care of projection effects of nearby halos. The total kSZE flux within halos estimated implies that the gas fraction in halos is about the universal baryon fraction, even in low-mass halos, indicating that the `missing baryons' are found. Various tests performed show that our results are robust against systematic effects, such as contamination by infrared/radio sources and background variations, beam-size effects and contributions from halo exteriors. Combined with the thermal Sunyaev-Zel'dovich effect, our results indicate that the `missing baryons' associated with galaxy groups are contained in warm-hot media with temperatures between $10^5$ and $10^6\,{\rm K}$.

We study the relation between the properties of the bulge/disc components and the HI mass fraction of galaxies. We find that at fixed stellar mass, disc colours are correlated with the HI mass fraction, while bulge colours are not. The lack of a correlation between the bulge colour and the HI mass fraction is regardless whether the bulges are pseudo, or whether the galaxies host bars or are interacting with a neighbour. There is no strong correlation between the colours of the discs and bulges either. These results suggest that the current total amount of HI is closely related to the formation of discs, but does not necessarily fuel the formation of (pseudo) bulges in an efficient way. We do not find evidence for the star formation in the discs to be quenched by the bulges.

As demonstrated in Paper I, the quenching properties of central and satellite galaxies are quite similar as long as both stellar mass and halo mass are controlled. Here we extend the analysis to the size and bulge-to-total light ratio (B/T) of galaxies. In general central galaxies have size-stellar mass and B/T-stellar mass relations different from satellites. However, the differences are eliminated when halo mass is controlled. We also study the dependence of size and B/T on halo-centric distance and find a transitional stellar mass (M$_{*,t}$) at given halo mass (M$_h$), which is about one fifth of the mass of the central galaxies in halos of mass M$_h$. The transitional stellar masses for size, B/T and quenched fraction are similar over the whole halo mass range, suggesting a connection between the quenching of star formation and the structural evolution of galaxies. Our analysis further suggests that the classification based on the transitional stellar mass is more fundamental than the central-satellite dichotomy, and provide a more reliable way to understand the environmental effects on galaxy properties. We compare the observational results with the hydro-dynamical simulation, EAGLE and the semi-analytic model, L-GALAXIES. The EAGLE simulation successfully reproduces the similarities of size for centrals and satellites and even M$_{*,t}$, while L-GALAXIES fails to recover the observational results.

The three-point correlation function (3PCF) of the 21cm brightness temperature from the Epoch of Reionization (EoR) probes complementary information to the commonly studied two-point correlation function (2PCF) about the morphology of ionized regions. We investigate the 21cm 2PCF and 3PCF in configuration space using semi-numerical simulations and test whether they can be described by the local quadratic bias model. We find that fits of bias model predictions for the 2PCF and 3PCF deviate from our measurements by $\sim 20\%$ at scales above the typical size of ionized regions ($\simeq 30$ Mpc) and at early times with global neutral fractions of $\langle x_{\rm HI} \rangle \gtrsim 0.7$. At later times and smaller scales these deviations increase strongly, indicating a break down of the bias model. The 2PCF and 3PCF fits of the linear bias parameter agree at the $10\%$ level for different EoR model configurations. This agreement holds, when adding redshift space distortions to the simulations. The relation between spatial fluctuations in the matter density and the 21cm signal, as predicted by the bias model, is consistent with direct measurements of this relation in simulations for large smoothing scales ($\gtrsim 30$ Mpc). From this latter test we conclude that negative amplitudes of the 21cm 3PCF result from negative bias parameters, which describe the anti-correlation between the matter over-densities and the 21cm signal during the EoR. However, a more detailed interpretation of the bias parameters may require a description of non-local contributions to the bias model.

We develop a semi-analytic model to explore the physical properties of cool pressure-confined circumgalactic clouds with mass ranging from $10$ to $10^{8} \, \rm M_{\odot}$ in a hot diffuse halo. We consider physical effects that control the motion and mass loss of the clouds, and estimate the lifetime and the observed properties of individual cool gas clouds inferred from the CLOUDY simulation. Our results show that the cool pressure-confined gas clouds have physical properties consistent with absorption line systems with neutral hydrogen column densities $N_{\rm HI}\geq10^{18.5} \rm cm^{-2}$ such as strong metal absorbers, sub-DLAs, and DLAs. The cool circumgalactic clouds are transient due to evaporation and recycling and therefore a constant replenishment is needed to maintain the cool CGM. We further model the ensemble properties of the cool CGM with clouds originated from outflows, inflows, or/and in-situ formation with a range of initial cloud mass function and velocity distribution. We find that only with a certain combination of parameters, an outflow model can broadly reproduce three cool gas properties around star-forming galaxies simultaneously: the spatial distribution, down-the-barrel outflow absorption, and gas velocity dispersion. Both a constant insitu model and gas inflow model can reproduce the observed covering fractions of high $\rm N_{HI}$ gas around passive galaxies but they fail to reproduce sufficient number of low $\rm N_{HI}$ systems. The limitations and the failures of the current models are discussed. Our results illustrate that semi-analytic modeling is a promising tool to understand the physics of the cool CGM which is usually unresolved by state-of-the-art cosmological hydrodynamic simulations.

We investigate the clustering properties and close neighbour counts for galaxies with different types of bulges and stellar masses. We select samples of "classical" and "pseudo" bulges, as well as "bulge-less" disk galaxies, based on the bulge/disk decomposition catalog of SDSS galaxies provided by Simard et al. (2011). For a given galaxy sample we estimate: the projected two-point cross-correlation function with respect to a spectroscopic reference sample, w_p(r_p), and the average background-subtracted neighbour count within a projected separation using a photometric reference sample, N_neighbour(<r_p). We compare the results with the measurements of control samples matched in color, concentration and redshift. We find that, when limited to a certain stellar mass range and matched in color and concentration, all the samples present similar clustering amplitudes and neighbour counts on scales above ~0.1h^{-1}Mpc. This indicates that neither the presence of a central bulge, nor the bulge type is related to intermediate-to-large scale environments. On smaller scales, in contrast, pseudo-bulge and pure-disk galaxies similarly show strong excess in close neighbour count when compared to control galaxies, at all masses probed. For classical bulges, small-scale excess is also observed but only for M_stars < 10^{10} M_sun; at higher masses, their neighbour counts are similar to that of control galaxies at all scales. These results imply strong connections between galactic bulges and galaxy-galaxy interactions in the local Universe, although it is unclear how they are physically linked in the current theory of galaxy formation.

We analyze the stellar initial mass functions (IMF) of a large sample of early type galaxies (ETGs) provided by MaNGA. The large number of IFU spectra of individual galaxies provide high signal-to-noise composite spectra that are essential for constraining IMF and to investigate possible radial gradients of the IMF within individual galaxies. The large sample of ETGs also make it possible to study how the IMF shape depends on various properties of galaxies. We adopt a novel approach to study IMF variations in ETGs, use Bayesian inferences based on full spectrum fitting. The Bayesian method provides a statistically rigorous way to explore potential degeneracy in spectrum fitting, and to distinguish different IMF models with Bayesian evidence. We find that the IMF slope depends systematically on galaxy velocity dispersion, in that galaxies of higher velocity dispersion prefer a more bottom-heavy IMF, but the dependence is almost entirely due to the change of metallicity, $Z$, with velocity dispersion. The IMF shape also depends on stellar age, $A$, but the dependence is completely degenerate with that on metallicity through a combination $AZ^{-1.42}$. Using independent age and metallicity estimates we find that the IMF variation is produced by metallicity instead of age. The IMF near the centers of massive ETGs appears more bottom-heavy than that in the outer parts, while a weak opposite trend is seen for low-mass ETGs. Uncertainties produced by star formation history, dust extinction, $\alpha$-element abundance enhancement and noise in the spectra are tested.

We use a sample of ~3000 galaxies from the MaNGA MPL-7 internal data release to study the alpha abundance distribution within low-redshift galaxies. We use the Lick index ratio Mgb/<Fe> as an alpha abundance indicator to study relationships between the alpha abundance distribution and galaxy properties such as effective stellar velocity dispersion within 0.3 effective radii (sigma_*), galaxy environment, and dark matter halo formation time (z_f). We find that (1) all galaxies show a tight correlation between Mgb/<Fe> and sigma_*; (2) `old' (H_beta < 3) low-sigma_* galaxies in high local density environment and inner regions within galaxy groups are enhanced in Mgb/<Fe>, while `young' (H_beta>3) galaxies and high-mass galaxies show no or less environmental dependence; (3) `old' galaxies with high-z_f show enhanced Mgb/<Fe> over low- and medium-z_f; (4) Mgb/<Fe> gradients are close to zero and show dependence on sigma_* but no obvious dependence on the environment or z_f. Our study indicates that stellar velocity dispersion or galaxy mass is the main parameter driving the Mgb/<Fe> enhancement, although environments appear to have modest effects, particularly for low- and medium-mass galaxies.

We study the spatially resolved star formation history and molecular gas distribution of 58 nearby galaxies, using integral field spectroscopy from the CALIFA survey and CO $J=1\rightarrow 0$ intensity mapping from the CARMA EDGE survey. We use the 4000 \AA\ break (D$_n$4000), the equivalent width of the H$\delta$ absorption line (EW H$\delta_A$), and the equivalent width of the H$\alpha$ emission line (EW H$\alpha$) to measure the recent star formation history (SFH) of these galaxies. We measure radial profiles of the three SFH indicators and molecular gas mass surface density, from which we measure the level of centrally enhanced star formation and the molecular gas concentration. When we separate our galaxies into categories of barred (17 galaxies), unbarred (24 galaxies), and merging/paired (17 galaxies) we find that the galaxies which have centrally-enhanced star formation (19/58) are either barred (13/19) or in mergers/pairs (6/19) with relatively high molecular gas concentrations. A comparison between our barred galaxies and a snapshot of a hydrodynamic $N$-body simulation of a barred galaxy shows that the current theory of bar formation and evolution can qualitatively reproduce the main features of the observed galaxies in our sample, including both the sharp decrease of stellar age in the galactic center and the gradual decrease of age with increasing distance from center. These findings provide substantial evidence for a picture in which cold gas is transported inward by a bar or tidal interaction, which leads to the growth and rejuvenation of star formation in the central region.

We study the distribution and kinematics of the cool circumgalactic medium (CGM) of emission line galaxies (ELGs) traced by metal absorption lines. Using about 200,000 ELGs from SDSS-IV eBOSS and half a million background quasars from SDSS, we measure the median absorption strength of MgII and FeII lines in quasar spectra for impact parameters ranging from 10 kpc to 1 Mpc. For comparison we measure the same quantity around luminous red galaxies (LRGs). On scales greater than 100 kpc both ELGs and LRGs exhibit similar absorption profiles. However, metal absorption is 5-10 times stronger around ELGs on smaller scales. The metal absorption strength is anisotropic, with an excess along the minor-axis of the galaxies, indicating an outflow origin of the absorbing gas. The ratio between the velocity dispersion of the cool CGM and that of its host dark matter halo is about one for ELGs but about half for LRGs. These results show that the dichotomy of galaxy types is reflected in both the density distribution and kinematics of the CGM traced by metal absorption lines. Our results provide strong evidence that the CGM of ELGs is enriched by gas outflows generated by star formation.

Understanding how galaxies maintain the inefficiency of star formation with physically self-consistent models is a central problem for galaxy evolution. Although numerous theoretical models have been proposed in recent decades, the debate still exists. By means of high-resolution two-dimensional hydrodynamical simulations, we study the three feedback effects (the stellar wind heating, SNe feedback, and AGN feedback) in suppressing star formation activities on the evolution of early-type galaxies with different stellar masses. AGN feedback models are updated based on \citet{Yuan2018}. The gas sources comes exclusively from the mass losses of dying low-mass stars for most of our models. We find that SNe feedback can keep star formation at a significantly low level for low mass elliptical galaxies for a cosmological evolution time. For the high mass galaxies, AGN feedback can efficiently offset the radiative cooling and thus regulate the star formation activities. Such a suppression of star formation is extremely efficient in the inner region of the galaxies. AGB heating cannot account for this suppression for low and high mass galaxies. The X-ray temperature $T_{\rm X}$ and luminosity $L_{\rm X}$ of hot plasma can be in agreement with the observed data with the inclusion of effective feedback processes. These results thus suggest that we can use $T_{\rm X}$ and $L_{\rm X}$ to probe the role of different feedback processes. The inclusion of additional gas sources can make the mass scale between SNe and AGN feedback dominating in suppressing star formation decrease to an observationally inferred value of a few $10^{10}~M_{\odot}$.

Halo merger trees are constructed from ELUCID, a constrained $N$-body simulation in the Sloan Digital Sky Survey (SDSS) volume. These merger trees are used to populate dark matter halos with galaxies according to an empirical model of galaxy formation. Mock catalogs in the SDSS sky coverage are constructed, which can be used to study the spatial distribution of galaxies in the low-$z$ Universe. These mock catalogs are used to quantify the cosmic variance in the galaxy stellar mass function (GSMF) measured from the SDSS survey. The GSMF estimated from the SDSS magnitude-limited sample can be affected significantly by the presence of the under-dense region at $z<0.03$, so that the low-mass end of the function can be underestimated significantly. Several existing methods designed to deal with the effects of the cosmic variance in the estimate of GSMF are tested, and none is found to be able to fully account for the cosmic variance. We propose a method based on the conditional stellar mass functions in dark matter halos, which can provide an unbiased estimate of the global GSMF. The application of the method to the SDSS data shows that the GSMF has a significant upturn at $M_*< 10^{9.5} h^{-1}{\rm M}_\odot$, which has been missed in many earlier measurements of the local GSMF.

As we demonstrated in Paper I, the quenched fractions of central and satellite galaxies as function of halo mass are extremely similar, as long as one controls for stellar mass. The same holds for the quenched fractions as a function of central velocity dispersion, which is tightly correlated with black hole mass, as long as one controls for both stellar and halo mass. Here we use mock galaxy catalogs constructed from the latest semi-analytic model, L-GALAXIES, and the state-of-the-art hydrodynamical simulation, EAGLE, to investigate whether these models can reproduce the trends seen in the data. We also check how the group finder used to identify centrals and satellites impacts our results. We find that L-GALAXIES fails to reproduce the trends. The predicted quenched fraction of central galaxies increases sharply with halo mass around $10^{12.5}h^{-1}M_{\odot}$ and with black hole mass around $\sim10^{6.5}M_{\odot}$, while the predicted quenched fraction of satellites increases with both halo and black hole masses gradually. In contrast, centrals and satellites in EAGLE follow almost the same trend as seen in the data. We discuss the implications of our results for how feedback processes regulate galaxy quenching.

We extend the real-space mapping method developed in Shi et at. (2016) so that it can be applied to flux-limited galaxy samples. We use an ensemble of mock catalogs to demonstrate the reliability of this extension, showing that it allows for an accurate recovery of the real-space correlation functions and galaxy biases. We also demonstrate that, using an iterative method applied to intermediate-scale clustering data, we can obtain an unbiased estimate of the growth rate of structure $f\sigma_8$, which is related to the clustering amplitude of matter, to an accuracy of $\sim 10\%$. Applying this method to the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7), we construct a real-space galaxy catalog spanning the redshift range $0.01 \leq z \leq 0.2$, which contains 584,473 galaxies in the north Galactic cap (NGC). Using this data, we infer $\fss$ at a median redshift $z=0.1$, which is consistent with the WMAP9 cosmology at the $1\sigma$ level. By combining this measurement with the real-space clustering of galaxies and with galaxy-galaxy weak lensing measurements for the same sets of galaxies, we are able to break the degeneracy between $f$, $\sigma_8$, and $b$. From the SDSS DR7 data alone, we obtain the following cosmological constraints at redshift $z=0.1$: $f=$$0.464^{+0.040}_{-0.040}$, $\sigma_8=0.769^{+0.121}_{-0.089}$, and $b=1.910^{+0.234}_{-0.268}$, $1.449^{+0.194}_{-0.196}$, $1.301^{+0.170}_{-0.177}$, and $1.196^{+0.159}_{-0.161}~$ for galaxies within different absolute magnitude bins $^{0.1}{\rm M}_r-5\log h=[-23,0, -22.0], [-22,0, -21.0], [-21.0, -20.0]$ and $[-20.0, -19.0]$, respectively.

We investigate the origin of the type III (up-bending) discs based on a sample of 1808 galaxies from SDSS and a sub-sample of 286 galaxies with HI data from ALFALFA. We examine how the type III fraction f3, the fraction of disc galaxies which host up-bending disc breaks, depends on other galactic properties. We confirm that f3 strongly depends on the stellar concentration of galaxies. We find that HI-rich galaxies with low spins tend to have significantly more type III disc breaks than control galaxies, which are matched in concentration and stellar mass. This trend is independent of the existence of strong bars or environment of the galaxies. This result is broadly consistent with predictions from theoretical simulations, and indicates in-situ star formation fueled by gas accretion to be an important mechanism that builds the outer discs of type III galaxies.

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.

We investigate the quenching properties of central and satellite galaxies, utilizing the halo masses and central-satellite identifications from the SDSS galaxy group catalog of Yang et al. We find that the quenched fractions of centrals and satellites of similar stellar masses have similar dependence on host halo mass. The similarity of the two populations is also found in terms of specific star formation rate and 4000 \AA\ break. The quenched fractions of centrals and satellites of similar masses show similar dependencies on bulge-to-total light ratio, central velocity dispersion and halo-centric distance in halos of given halo masses. The prevalence of optical/radio-loud AGNs is found to be similar for centrals and satellites at given stellar masses. All these findings strongly suggest that centrals and satellites of similar masses experience similar quenching processes in their host halos. We discuss implications of our results for the understanding of galaxy quenching.

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 examine the thermal energy contents of the intergalactic medium (IGM) over three orders of magnitude in both mass density and gas temperature using thermal Sunyaev-Zel'dovich effect (tSZE). The analysis is based on {\it Planck} tSZE map and the cosmic density field, reconstructed for the SDSS DR7 volume and sampled on a grid of cubic cells of $(1h^{-1}{\rm Mpc})^3$, together with a matched filter technique employed to maximize the signal-to-noise. Our results show that the pressure - density relation of the IGM is roughly a power law given by an adiabatic equation of state, with an indication of steepening at densities higher than about $10$ times the mean density of the universe. The implied average gas temperature is $\sim 10^4\,{\rm K}$ in regions of mean density, $\rho_{\rm m} \sim {\overline\rho}_{\rm m}$, increasing to about $10^5\,{\rm K}$ for $\rho_{\rm m} \sim 10\,{\overline\rho}_{\rm m}$, and to $>10^{6}\,{\rm K}$ for $\rho_{\rm m} \sim 100\,{\overline\rho}_{\rm m}$. At a given density, the thermal energy content of the IGM is also found to be higher in regions of stronger tidal fields, likely due to shock heating by the formation of large scale structure and/or feedback from galaxies and AGNs. A comparison of the results with hydrodynamic simulations suggests that the current data can already provide interesting constraints on galaxy formation.

We use a 200 $h^{-1}Mpc$ a side N-body simulation to study the mass accretion history (MAH) of dark matter halos to be accreted by larger halos, which we call infall halos. We define a quantity $a_{\rm nf}\equiv (1+z_{\rm f})/(1+z_{\rm peak})$ to characterize the MAH of infall halos, where $z_{\rm peak}$ and $z_{\rm f}$ are the accretion and formation redshifts, respectively. We find that, at given $z_{\rm peak}$, their MAH is bimodal. Infall halos are dominated by a young population at high redshift and by an old population at low redshift. For the young population, the $a_{\rm nf}$ distribution is narrow and peaks at about $1.2$, independent of $z_{\rm peak}$, while for the old population, the peak position and width of the $a_{\rm nf}$ distribution both increases with decreasing $z_{\rm peak}$ and are both larger than those of the young population. This bimodal distribution is found to be closely connected to the two phases in the MAHs of halos. While members of the young population are still in the fast accretion phase at $z_{\rm peak}$, those of the old population have already entered the slow accretion phase at $z_{\rm peak}$. This bimodal distribution is not found for the whole halo population, nor is it seen in halo merger trees generated with the extended Press-Schechter formalism. The infall halo population at $z_{\rm peak}$ are, on average, younger than the whole halo population of similar masses identified at the same redshift. We discuss the implications of our findings in connection to the bimodal color distribution of observed galaxies and to the link between central and satellite galaxies.

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre<outer$), while the gradients for most massive galaxies are negative. The metallicity gradients show a clear peak around velocity dispersion $\log_{10} \sigma_{\rm e}\approx 2.0$, which corresponds to the critical mass $\sim 3\times 10^{10}M_{\odot}$ of the break in the mass-size relation. Spiral galaxies with large mass and size have the steepest gradients, while the most massive ETGs, especially above the critical mass $M_{\rm crit}\ge 2\times 10^{11} M_{\odot}$, where slow rotator ETGs start dominating, have much flatter gradients. This may be due to differences in their evolution histories, e.g. mergers.

As the second paper of a series on studying galaxy-galaxy lensing signals using the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our measurement and modelling of the lensing signals around groups of galaxies. We divide the groups into four halo mass bins, and measure the signals around four different halo-center tracers: brightest central galaxy (BCG), luminosity-weighted center, number-weighted center and X-ray peak position. For X-ray and SDSS DR7 cross identified groups, we further split the groups into low and high X-ray emission subsamples, both of which are assigned with two halo-center tracers, BCGs and X-ray peak positions. The galaxy-galaxy lensing signals show that BCGs, among the four candidates, are the best halo-center tracers. We model the lensing signals using a combination of four contributions: off-centered NFW host halo profile, sub-halo contribution, stellar contribution, and projected 2-halo term. We sample the posterior of 5 parameters i.e., halo mass, concentration, off-centering distance, sub halo mass, and fraction of subhalos via a MCMC package using the galaxy-galaxy lensing signals. After taking into account the sampling effects (e.g. Eddington bias), we found the best fit halo masses obtained from lensing signals are quite consistent with those obtained in the group catalog based on an abundance matching method, except in the lowest mass bin. Subject headings: (cosmology:) gravitational lensing, galaxies: clusters: general

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 apply a halo-based group finder to four large redshift surveys, the 2MRS, 6dFGS, SDSS and 2dFGRS, to construct group catalogs in the low-redshift Universe. The group finder is based on that of Yang et al. but with an improved halo mass assignment so that it can be applied uniformly to various redshift surveys of galaxies. Halo masses are assigned to groups according to proxies based on the stellar mass/luminosity of member galaxies. The performances of the group finder in grouping galaxies according to common halos and in halo mass assignments are tested using realistic mock samples constructed from hydrodynamical simulations and empirical models of galaxy occupation in dark matter halos. Our group finder finds $\sim 94\%$ of the correct true member galaxies for $90-95\%$ of the groups in the mock samples; the halo masses assigned by the group finder are un-biased with respect to the true halo masses, and have a typical uncertainty of $\sim0.2\,{\rm dex}$. The properties of group catalogs constructed from the observational samples are described and compared with other similar catalogs in the literature.

A matched filter technique is applied to the Planck all-sky Compton y-parameter map to measure the thermal Sunyaev-Zel'dovich (tSZ) effect produced by galaxy groups of different halo masses selected from large redshift surveys in the low-z Universe. Reliable halo mass estimates are available for all the groups, which allows us to bin groups of similar halo masses to investigate how the tSZ effect depends on halo mass over a large mass range. Filters are simultaneously matched for all groups to minimize projection effects. We find that the integrated y-parameter and the hot gas content it implies are consistent with the predictions of the universal pressure profile model only for massive groups above $10^{14}\,{\rm M}_\odot$, but much lower than the model prediction for low-mass groups. The halo mass dependence found is in good agreement with the predictions of a set of simulations that include strong AGN feedback, but simulations including only supernova feedback significantly over predict the hot gas contents in galaxy groups. Our results suggest that hot gas in galaxy groups is either effectively ejected or in phases much below the virial temperatures of the host halos.

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 show that the ratio between the stellar mass of central galaxy and the mass of its host halo, $f_c \equiv M_{*,c}/M_{\rm h}$, can be used as an observable proxy of halo assembly time, in that galaxy groups with higher $f_c$ assembled their masses earlier. Using SDSS groups of Yang et al., we study how $f_c$ correlates with galaxy properties such as color, star formation rate, metallicity, bulge to disk ratio, and size. Central galaxies of a given stellar mass in groups with $f_c>0.02$ tend to be redder in color, more quenched in star formation, smaller in size, and more bulge dominated, as $f_c$ increases. The trends in color and star formation appear to reverse at $f_c<0.02$, reflecting a down-sizing effect that galaxies in massive halos formed their stars earlier although the host halos themselves assembled later (lower $f_c$). No such reversal is seen in the size of elliptical galaxies, suggesting that their assembly follows halo growth more closely than their star formation. Satellite galaxies of a given stellar mass in groups of a given halo mass tend to be redder in color, more quenched in star formation and smaller in size as $f_c$ increases. For a given stellar mass, satellites also tend to be smaller than centrals. The trends are stronger for lower mass groups. For groups more massive than $\sim 10^{13}{\rm M}_\odot$, a weak reversed trend is seen in color and star formation. The observed trends in star formation are qualitatively reproduced by an empirical model based on halo age abundance matching, but not by a semi-analytical model tested here.

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.

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

Using a method to correct redshift space distortion (RSD) for individual galaxies, we mapped the real space distributions of galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7). We use an ensemble of mock catalogs to demonstrate the reliability of our method. Here as the first paper in a series, we mainly focus on the two point correlation function (2PCF) of galaxies. Overall the 2PCF measured in the reconstructed real space for galaxies brighter than $^{0.1}{\rm M}_r-5\log h=-19.0$ agrees with the direct measurement to an accuracy better than the measurement error due to cosmic variance, if the reconstruction uses the correct cosmology. Applying the method to the SDSS DR7, we construct a real space version of the main galaxy catalog, which contains 396,068 galaxies in the North Galactic Cap with redshifts in the range $0.01 \leq z \leq 0.12$. The Sloan Great Wall, the largest known structure in the nearby Universe, is not as dominant an over-dense structure as appears to be in redshift space. We measure the 2PCFs in reconstructed real space for galaxies of different luminosities and colors. All of them show clear deviations from single power-law forms, and reveal clear transitions from 1-halo to 2-halo terms. A comparison with the corresponding 2PCFs in redshift space nicely demonstrates how RSDs boost the clustering power on large scales (by about $40-50\%$ at scales $\sim 10 h^{-1}{\rm {Mpc}}$) and suppress it on small scales (by about $70-80\%$ at a scale of $0.3 h^{-1}{\rm {Mpc}}$).

We compare predictions of a number of empirical models and numerical simulations of galaxy formation to the conditional stellar mass functions (CSMF)of galaxies in groups of different masses obtained recently by Lan et al. to test how well different models accommodate the data. The observational data clearly prefer a model in which star formation in low-mass halos changes behavior at a characteristic redshift $z_c\sim 2$. There is also tentative evidence that this characteristic redshift depends on environment, becoming $z_c\sim 4$ in regions that eventually evolve into rich clusters of galaxies. The constrained model is used to understand how galaxies form and evolve in dark matter halos, and to make predictions for other statistical properties of the galaxy population, such as the stellar mass functions of galaxies at high $z$, the star formation and stellar mass assembly histories in dark matter halos. A comparison of our model predictions with those of other empirical models shows that different models can make vastly different predictions, even though all of them are tuned to match the observed stellar mass functions of galaxies.

Using low-redshift (z<0.09) samples of AGNs, normal galaxies and groups of galaxies selected from the Sloan Digital Sky Survey (SDSS), we study the environments of type 1 and type 2 AGNs both on small and large scales. Comparisons are made for galaxy samples matched in redshift, $r$-band luminosity, [OIII] luminosity, and also the position in groups (central or satellite). We find that type 2 AGNs and normal galaxies reside in similar environments. Type 1 and type 2 AGNs have similar clustering properties on large scales ($\gtrsim1$Mpc), but at scales smaller than 100 kpc, type 2s have significant more neighbors than type 1s ($3.09\pm0.69$ times more for central AGNs at $\lesssim30$kpc). These results suggest that type 1 and type 2 AGNs are hosted by halos of similar masses, as is also seen directly from the mass distributions of their host groups ($\sim10^{12}h^{-1} M_{\odot}$ for centrals and $\sim10^{13}h^{-1} M_{\odot}$ for satellites). Type~2s have significantly more satellites around them, and the distribution of their satellites is also more centrally concentrated. The host galaxies of both types of AGNs have similar optical properties, but their infrared colors are significantly different. Our results suggest that the simple unified model based solely on torus orientation is not sufficient, but that galaxy interactions in dark matter halos must have played an important role in the formation of the dust structure that obscures AGNs.

A galaxy group catalog is constructed from the 2MASS Redshift Survey (2MRS) with the use of a halo-based group finder. The halo mass associated with a group is estimated using a `GAP' method based on the luminosity of the central galaxy and its gap with other member galaxies. Tests using mock samples shows that this method is reliable, particularly for poor systems containing only a few members. On average 80% of all the groups have completeness >0.8, and about 65% of the groups have zero contamination. Halo masses are estimated with a typical uncertainty $\sim 0.35\,{\rm dex}$. The application of the group finder to the 2MRS gives 29,904 groups from a total of 43,246 galaxies at $z \leq 0.08$, with 5,286 groups having two or more members. Some basic properties of this group catalog is presented, and comparisons are made with other groups catalogs in overlap regions. With a depth to $z\sim 0.08$ and uniformly covering about 91% of the whole sky, this group catalog provides a useful data base to study galaxies in the local cosmic web, and to reconstruct the mass distribution in the local Universe.

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.

As the first paper in a series on the study of the galaxy-galaxy lensing from Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we present our image processing pipeline that corrects the systematics primarily introduced by the Point Spread Function (PSF). Using this pipeline, we processed SDSS DR7 imaging data in $r$ band and generated a background galaxy catalog containing the shape information of each galaxy. Based on our own shape measurements of the galaxy images from SDSS DR7, we extract the galaxy-galaxy (GG) lensing signals around foreground spectroscopic galaxies binned in different luminosity and stellar mass. The overall signals are in good agreement with those obtained by \citet{Mandelbaum2005, Mandelbaum2006} from the SDSS DR4. The results in this paper with higher signal to noise ratio is due to the larger survey area than SDSS DR4, confirm that more luminous/massive galaxies bear stronger GG lensing signal. We also divide the foreground galaxies into red/blue and star forming/quenched subsamples and measured their GG lensing signals, respectively. We find that, at a specific stellar mass/luminosity, the red/quenched galaxies have relatively stronger GG lensing signals than their counterparts especially at large radii. These GG lensing signals can be used to probe the galaxy-halo mass relations and their environmental dependences in the halo occupation or conditional luminosity function framework.

We characterize the luminosity functions of galaxies residing in $z\sim0$ groups and clusters over the broadest ranges of luminosity and mass reachable by the Sloan Digital Sky Survey. Our measurements cover four orders of magnitude in luminosity, down to about $M_r=-12$ mag or $L=10^7\,L_\odot$, and three orders of magnitude in halo mass, from $10^{12}$ to $10^{15} \, {\rm M}_\odot$. We find a characteristic scale, $M_r\sim-18$ mag or $L\sim10^9\, L_\odot$, below which the slope of the luminosity function becomes systematically steeper. This trend is present for all halo masses and originates mostly from red satellites. This ubiquitous faint-end upturn suggests that it is formation, rather than halo-specific environmental effect, that plays a major role in regulating the stellar masses of faint satellites. We show that the satellite luminosity functions can be described in a simple manner by a double Schechter function with amplitudes scaling with halo mass over the entire range of observables. Combining these conditional luminosity functions with the dark matter halo mass function, we accurately recover the entire field luminosity function over 10 visual magnitudes and reveal that satellite galaxies dominate the field luminosity function at magnitudes fainter than $-17$. We find that the luminosity functions of blue and red satellite galaxies show distinct shapes and we present estimates of the stellar mass fraction as a function of halo mass and galaxy type. Finally, using a simple model, we demonstrate that the abundances and the faint-end slopes of blue and red satellite galaxies can be interpreted in terms of their formation history, with two distinct modes separated by some characteristic time.

Large scale tidal field estimated directly from the distribution of dark matter halos is used to investigate how halo shapes and spin vectors are aligned with the cosmic web. The major, intermediate and minor axes of halos are aligned with the corresponding tidal axes, and halo spin axes tend to be parallel with the intermediate axes and perpendicular to the major axes of tidal field. The strengths of these alignments generally increase with halo mass and redshift, but the dependencies are only through the peak height, {\nu}. The scaling relations of the alignment strengths with the value of {\nu} indicate that the alignment strengths remain roughly constant when the structures within which the halos reside are still in quasi-linear regime, but decreases as nonlinear evolution becomes more important. We also calculate the alignments in projection so that our results can be compared directly with observations. Finally, we investigate the alignments of tidal tensors on large scales, and use the results to understand alignments of halo pairs separated at various distances. Our results suggest coherent structure of the tidal field is the underlying reason for the alignments of halos and galaxies seen in numerical simulations and in observations.

We use the shear catalog from the CFHT Stripe-82 Survey to measure the subhalo masses of satellite galaxies in redMaPPer clusters. Assuming a Chabrier Initial Mass Function (IMF) and a truncated NFW model for the subhalo mass distribution, we find that the sub-halo mass to galaxy stellar mass ratio increases as a function of projected halo-centric radius $r_p$, from $M_{\rm sub}/M_{\rm star}=4.43^{+ 6.63}_{- 2.23}$ at $r_p \in [0.1,0.3]$ $h^{-1}Mpc$ to $M_{\rm sub}/M_{\rm star}=75.40^{+ 19.73}_{- 19.09}$ at $r_p \in [0.6,0.9]$ $h^{-1}Mpc$. We also investigate the dependence of subhalo masses on stellar mass by splitting satellite galaxies into two stellar mass bins: $10<\log(M_{\rm star}/M_{\rm sun})<10.5$ and $11<\log(M_{\rm star}/M_{\rm sun})<12$. The best-fit subhalo mass of the more massive satellite galaxy bin is larger than that of the less massive satellites: $\log(M_{\rm sub}/M_{\rm sun})=11.14 ^{+ 0.66 }_{- 0.73}$ ($M_{\rm sub}/M_{\rm star}=19.5^{+19.8}_{-17.9}$) versus $\log(M_{\rm sub}/M_{\rm sun})=12.38 ^{+ 0.16 }_{- 0.16}$ ($M_{\rm sub}/M_{\rm star}=21.1^{+7.4}_{-7.7}$).

We study the quenching of satellite galaxies by gradual depletion of gas due to star formation, by ram-pressure striping and by tidally triggered starburst. Using progenitors constrained by the empirical model of Lu et al., in which outflow loading factor is low, we do not find an over-quenching problem in satellites even if there is no further cold gas supply from the cooling of the halo gas after a galaxy is accreted by its host. Gradual depletion alone predicts a unimodal distribution in specific star formation, in contrast to the bimodal distribution observed, and under-predicts the quenched fraction in low mass halos. Ram-pressure stripping nicely reproduces the bimodal distribution but under-predicts the quenched fraction in low-mass halos. Starbursts in gas-rich satellites triggered by tidal interactions with central galaxies can nicely reproduce the quenched satellite population in low-mass halos, but become unimportant for low-mass satellites in massive halos. The combined processes, together with the constrained progenitors, can reproduce the observed star formation properties of satellites in halos of different masses.

We investigate the ability of current implementations of galaxy group finders to recover colour-dependent halo occupation statistics. To test the fidelity of group catalogue inferred statistics, we run three different group finders used in the literature over a mock that includes galaxy colours in a realistic manner. Overall, the resulting mock group catalogues are remarkably similar, and most colour-dependent statistics are recovered with reasonable accuracy. However, it is also clear that certain systematic errors arise as a consequence of correlated errors in group membership determination, central/satellite designation, and halo mass assignment. We introduce a new statistic, the halo transition probability (HTP), which captures the combined impact of all these errors. As a rule of thumb, errors tend to equalize the properties of distinct galaxy populations (i.e. red vs. blue galaxies or centrals vs. satellites), and to result in inferred occupation statistics that are more accurate for red galaxies than for blue galaxies. A statistic that is particularly poorly recovered from the group catalogues is the red fraction of central galaxies as function of halo mass. Group finders do a good job in recovering galactic conformity, but also have a tendency to introduce weak conformity when none is present. We conclude that proper inference of colour-dependent statistics from group catalogues is best achieved using forward modelling (i.e., running group finders over mock data), or by implementing a correction scheme based on the HTP, as long as the latter is not too strongly model-dependent.

We study how halo intrinsic dynamical properties are linked to their formation processes for halos in two mass ranges, $10^{12}-10^{12.5}h^{-1}{\rm M_\odot}$ and $\ge 10^{13}h^{-1}{\rm M_\odot}$, and how both are correlated with the large scale tidal field within which the halos reside at present. Halo merger trees obtained from cosmological $N$-body simulations are used to identify infall halos that are about to merge with their hosts. We find that the tangential component of the infall velocity increases significantly with the strength of the local tidal field, but no strong correlation is found for the radial component. These results can be used to explain how the internal velocity anisotropy and spin of halos depend on environment. The position vectors and velocities of infall halos are aligned with the principal axes of the local tidal field, and the alignment depends on the strength of the tidal field. Opposite accretion patterns are found in weak and strong tidal fields, in the sense that in a weak field the accretion flow is dominated by radial motion within the local structure, while a large tangential component is present in a strong field. These findings can be used to understand the strong alignments we find between the principal axes of the internal velocity ellipsoids of halos and the local tidal field, and their dependence on the strength of tidal field. They also explain why halo spin increases with the strength of local tidal field, but only in weak tidal fields does the spin-tidal field alignment follow the prediction of the tidal torque theory. We discuss how our results may be used to understand the spins of disk galaxies and velocity structures of elliptical galaxies and their correlations with large-scale structure.

Using a semi-analytical approach we investigate the characteristics of predictions for the masses and metallicities of the baryonic matter in and around galaxies made by three galaxy formation models. These models represent three different feedback scenarios: one model with purely ejective feedback, one model with ejective feedback with reincorporation of ejected gas, and one preventative model. We find that, when the model parameters are adjusted to predict the correct stellar masses for a range of halo masses between 10^{10} to 10^{12}Msun, these three scenarios have very different predictions for the masses and metallicities of the interstellar and circum-galactic media. Compared with current observational data, the model implementing preventative feedback has a large freedom to match a broad range of observational data, while the ejective models have difficulties to match a number of observational constraints simultaneously, independent of how the ejection and reincorporation are implemented. Our results suggest that the feedback process which regulates the amounts of stars and cold gas in low-mass galaxies is preventative in nature.

The empirical model of Lu et al. 2014 is updated with recent data and used to study galaxy star formation and assembly histories. At $z > 2$, the predicted galaxy stellar mass functions are steep, and a significant amount of star formation is hosted by low-mass haloes that may be missed in current observations. Most of the stars in cluster centrals formed earlier than $z\approx 2$ but have been assembled much later. Milky Way mass galaxies have had on-going star formation without significant mergers since $z\approx 2$, and are thus free of significant (classic) bulges produced by major mergers. In massive clusters, stars bound in galaxies and scattered in the halo form a homogeneous population that is old and with solar metallicity. In contrast, in Milky Way mass systems the two components form two distinct populations, with halo stars being older and poorer in metals by a factor of $\approx 3$. Dwarf galaxies in haloes with $M_{\rm h} < 10^{11}h^{-1}M_{\odot}$ have experienced a star formation burst accompanied by major mergers at $z > 2$, followed by a nearly constant star formation rate after $z = 1$. The early burst leaves a significant old stellar population that is distributed in spheroids.

The empirical model of Lu et al. 2014a for the relation between star formation rate and halo mass growth is adopted to predict the classical bulge mass ($M_{\rm cb}$) - total stellar mass ($M_\star$) relation for central galaxies. The assumption that the supermassive black hole (SMBH) mass ($M_{\rm BH}$) is directly proportional to the classical bulge mass, with the proportionality given by that for massive galaxies, predicts a $M_{\rm BH}$ - $M_\star$ relation that matches well the observed relation for different types of galaxies. In particular, the model reproduces the strong transition at $M_\star=10^{10.5}$ - $10^{11}M_{\odot}$, below which $M_{\rm BH}$ drops rapidly with decreasing $M_\star$. Our model predicts a new sequence at $M_\star <10^{10.5}M_{\odot}$, where $M_{\rm BH} \propto M_\star$ but the amplitude is a factor of $\sim 50$ lower than the amplitude of the sequence at $M_\star>10^{11}M_{\odot}$. If all SMBH grow through similar quasar modes with a feedback efficiency of a few percent, the energy produced in low-mass galaxies at redshift $z\gtrsim 2$ can heat the circum-galactic medium up to a specific entropy level that is required to prevent excessive star formation in low-mass dark matter halos.

We use a set of observational data for galaxy cold gas mass fraction and gas phase metallicity to constrain the content, inflow and outflow of gas in central galaxies hosted by halos with masses between $10^{11} M_{\odot}$ to $10^{12} M_{\odot}$. The gas contents in high redshift galaxies are obtained by combining the empirical star formation histories of Lu et al. (2014) and star formation models that relate star formation rate with the cold gas mass in galaxies. We find that the total baryon mass in low-mass galaxies is always much less than the universal baryon mass fraction since $z = 2$, regardless of star formation model adopted. The data for the evolution of the gas phase metallicity require net metal outflow at $z\lesssim 2$, and the metal loading factor is constrained to be about $0.01$, or about $60\%$ of the metal yield. Based on the assumption that galactic outflow is more enriched in metal than both the interstellar medium and the material ejected at earlier epochs, we are able to put stringent constraints on the upper limits for both the net accretion rate and the net mass outflow rate. The upper limits strongly suggest that the evolution of the gas phase metallicity and gas mass fraction for low-mass galaxies at $z < 2$ is not compatible with strong outflow. We speculate that the low star formation efficiency of low-mass galaxies is owing to some preventative processes that prevent gas from accreting into galaxies in the first place.

Using a sample of spiral galaxies selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and Galaxy Zoo 2 (GZ2), we investigate the alignment of spin axes of spiral galaxies with their surrounding large scale structure, which is characterized by the large-scale tidal field reconstructed from the data using galaxy groups above a certain mass threshold. We find that the spin axes of only have weak tendency to be aligned with (or perpendicular to) the intermediate (or minor) axis of the local tidal tensor. The signal is the strongest in a \cluster environment where all the three eigenvalues of the local tidal tensor are positive. Compared to the alignments between halo spins and local tidal field obtained in N-body simulations, the above observational results are in best agreement with those for the spins of inner regions of halos, suggesting that the disk material traces the angular momentum of dark matter halos in the inner regions.

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 infer mechanisms of galaxy formation for a broad family of semi-analytic models (SAMs) constrained by the K-band luminosity function and HI mass function of local galaxies using tools of Bayesian analysis. Even with a broad search in parameter space the whole model family fails to match to constraining data. In the best fitting models, the star formation and feedback parameters in low-mass haloes are tightly constrained by the two data sets, and the analysis reveals several generic failures of models that similarly apply to other existing SAMs. First, based on the assumption that baryon accretion follows the dark matter accretion, large mass-loading factors are required for haloes with circular velocities lower than 200 km/s, and most of the wind mass must be expelled from the haloes. Second, assuming that the feedback is powered by Type-II supernovae with a Chabrier IMF, the outflow requires more than 25% of the available SN kinetic energy. Finally, the posterior predictive distributions for the star formation history are dramatically inconsistent with observations for masses similar to or smaller than the Milky-Way mass. The inferences suggest that the current model family is still missing some key physical processes that regulate the gas accretion and star formation in galaxies with masses below that of the Milky Way.

We introduce a semi-analytic galaxy formation model implementing a self-consistent treatment for the hot halo gas configuration and the assembly of central disks. Using the model, we explore a preventative feedback model, in which the circum-halo medium is assumed to be preheated up to a certain entropy level by early starbursts or other processes, and compare it with an ejective feedback model, in which baryons are first accreted into dark matter halos and subsequently ejected out by feedback. The model demonstrates that when the medium is preheated to an entropy comparable to the halo virial entropy the baryon accretion can be largely reduced and delayed. In addition, the preheated medium can establish an extended low density gaseous halo when it accretes into the dark matter halos, and result in a specific angular momentum of the cooling gas large enough to form central disks as extended as those observed. Combined with simulated halo assembly histories, the preventative feedback model can reproduce remarkably well a number of observational scaling relations. These include the cold baryon (stellar plus cold gas) mass fraction-halo mass relations, star formation histories, disk size-stellar mass relation and its evolution, and the number density of low-mass galaxies as a function of redshift. In contrast, the conventional ejective feedback model fails to reproduce these observational trends. Using the model, we demonstrate that the properties of disk galaxies are closely tied to the thermal state of hot halo gas and even possibly the circum-halo medium, which suggests that observational data for the disk properties and circum-galactic hot/warm medium may jointly provide interesting constraints for galaxy formation models.

We develop an empirical approach to infer the star formation rate in dark matter halos from the galaxy stellar mass function (SMF) at different redshifts and the local cluster galaxy luminosity function (CGLF), which has a steeper faint end relative to the SMF of local galaxies. As satellites are typically old galaxies which have been accreted earlier, this feature can cast important constraint on the formation of low-mass galaxies at high-redshift. The evolution of the SMFs suggests the star formation in high mass halos ($>10^{12}M_{\odot}h^{-1}$) has to be boosted at high redshift beyond what is expected from a simple scaling of the dynamical time. The faint end of the CGLF implies a characteristic redshift $z_c\approx2$ above which the star formation rate in low mass halos with masses $< 10^{11}M_{\odot}h^{-1}$ must be enhanced relative to that at lower z. This is not directly expected from the standard stellar feedback models. Also, this enhancement leads to some interesting predictions, for instance, a significant old stellar population in present-day dwarf galaxies with $M_* < 10^8 M_{\odot}h^{-2}$ and steep slopes of high redshift stellar mass and star formation rate functions.

We examine the statistical properties of the brightest group galaxies (BGGs) using a complete spectroscopic sample of groups/clusters of galaxies selected from the Data Release 7 of the Sloan Digital Sky Survey. We test whether BGGs and other bright members of groups are consistent with an ordered population among the total population of group galaxies. We find that the luminosity distributions of BGGs do not follow the predictions from the order statistics (OS). The average luminosities of BGGs are systematically brighter than OS predictions. On the other hand, by properly taking into account the brightening effect of the BGGs, the luminosity distributions of the second brightest galaxies are in excellent agreement with the expectations of OS. The brightening of BGGs relative to the OS expectation is consistent with a scenario that the BGGs on average have over-grown about 20 percent masses relative to the other member galaxies. The growth ($\Delta M$) is not stochastic but correlated with the magnitude gap ($G_{1,2}$) between the brightest and the second brightest galaxy. The growth ($\Delta M$) is larger for the groups having more prominent BGGs (larger $G_{1,2}$) and averagely contributes about 30 percent of the final $G_{1,2}$ of the groups of galaxies.

We use ROSAT All Sky Survey (RASS) broadband X-ray images and the optical clusters identified from SDSS DR7 to estimate the X-ray luminosities around $\sim 65,000$ candidate clusters with masses $\ga 10^{13}\msunh$ based on an Optical to X-ray (OTX) code we develop. We obtain a catalogue with X-ray luminosity for each cluster. This catalog contains 817 clusters (473 at redshift $z\le 0.12$) with $S/N> 3$ in X-ray detection. We find about $65\%$ of these X-ray clusters have their most massive member located near the X-ray flux peak; for the rest $35\%$, the most massive galaxy is separated from the X-ray peak, with the separation following a distribution expected from a NFW profile. We investigate a number of correlations between the optical and X-ray properties of these X-ray clusters, and find that: the cluster X-ray luminosity is correlated with the stellar mass (luminosity) of the clusters, as well as with the stellar mass (luminosity) of the central galaxy and the mass of the halo, but the scatter in these correlations is large. Comparing the properties of X-ray clusters of similar halo masses but having different X-ray luminosities, we find that massive halos with masses $\ga 10^{14}\msunh$ contain a larger fraction of red satellite galaxies when they are brighter in X-ray. ... A cluster catalog containing the optical properties of member galaxies and the X-ray luminosity is available at {\it http://gax.shao.ac.cn/data/Group.html}.

We select satellite galaxies from the galaxy group catalog constructed with the SDSS spectroscopic galaxies and measure the tangential shear around these galaxies with source catalog extracted from CFHT/MegaCam Stripe-82 Survey to constrain the mass of subhalos associated with them. The lensing signal is measured around satellites in groups with masses in the range [10^{13}, 5x10^{14}]h^{-1}M_{sun}, and is found to agree well with theoretical expectation. Fitting the data with a truncated NFW profile, we obtain an average subhalo mass of log M_{sub}= 11.68 \pm 0.67 for satellites whose projected distances to central galaxies are in the range [0.1, 0.3] h^{-1}Mpc, and log M_{sub}= 11.68 \pm 0.76 for satellites with projected halo-centric distance in [0.3, 0.5] h^{-1}Mpc. The best-fit subhalo masses are comparable to the truncated subhalo masses assigned to satellite galaxies using abundance matching and about 5 to 10 times higher than the average stellar mass of the lensing satellite galaxies.

Using a sample of galaxy groups selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7), we examine the alignment between the orientation of galaxies and their surrounding large scale structure in the context of the cosmic web. The latter is quantified using the large-scale tidal field, reconstructed from the data using galaxy groups above a certain mass threshold. We find that the major axes of galaxies in filaments tend to be preferentially aligned with the directions of the filaments, while galaxies in sheets have their major axes preferentially aligned parallel to the plane of the sheets. The strength of this alignment signal is strongest for red, central galaxies, and in good agreement with that of dark matter halos in N-body simulations. This suggests that red, central galaxies are well aligned with their host halos, in quantitative agreement with previous studies based on the spatial distribution of satellite galaxies. There is a luminosity and mass dependence that brighter and more massive galaxies in filaments and sheets have stronger alignment signals. We also find that the orientation of galaxies is aligned with the eigenvector associated with the smallest eigenvalue of the tidal tensor. These observational results indicate that galaxy formation is affected by large-scale environments, and strongly suggests that galaxies are aligned with each other over scales comparable to those of sheets and filaments in the cosmic web.

Our research objective in this paper is to reconstruct an initial linear density field, which follows the multivariate Gaussian distribution with variances given by the linear power spectrum of the current CDM model and evolves through gravitational instability to the present-day density field in the local Universe. For this purpose, we develop a Hamiltonian Markov Chain Monte Carlo method to obtain the linear density field from a posterior probability function that consists of two components: a prior of a Gaussian density field with a given linear spectrum, and a likelihood term that is given by the current density field. The present-day density field can be reconstructed from galaxy groups using the method developed in Wang et al. (2009a). Using a realistic mock SDSS DR7, obtained by populating dark matter haloes in the Millennium simulation with galaxies, we show that our method can effectively and accurately recover both the amplitudes and phases of the initial, linear density field. To examine the accuracy of our method, we use $N$-body simulations to evolve these reconstructed initial conditions to the present day. The resimulated density field thus obtained accurately matches the original density field of the Millennium simulation in the density range 0.3 <= rho/rho_mean <= 20 without any significant bias. Especially, the Fourier phases of the resimulated density fields are tightly correlated with those of the original simulation down to a scale corresponding to a wavenumber of ~ 1 h/Mpc, much smaller than the translinear scale, which corresponds to a wavenumber of ~ 0.15 h\Mpc.

Using the self-consistent modeling of the conditional stellar mass functions across cosmic time by Yang et al. (2012), we make model predictions for the star formation histories (SFHs) of {\it central} galaxies in halos of different masses. The model requires the following two key ingredients: (i) mass assembly histories of central and satellite galaxies, and (ii) local observational constraints of the star formation rates of central galaxies as function of halo mass. We obtain a universal fitting formula that describes the (median) SFH of central galaxies as function of halo mass, galaxy stellar mass and redshift. We use this model to make predictions for various aspects of the star formation rates of central galaxies across cosmic time. Our main findings are the following. (1) The specific star formation rate (SSFR) at high $z$ increases rapidly with increasing redshift [$\propto (1+z)^{2.5}$] for halos of a given mass and only slowly with halo mass ($\propto M_h^{0.12}$) at a given $z$, in almost perfect agreement with the specific mass accretion rate of dark matter halos. (2) The ratio between the star formation rate (SFR) in the main-branch progenitor and the final stellar mass of a galaxy peaks roughly at a constant value, $\sim 10^{-9.3} h^2 {\rm yr}^{-1}$, independent of halo mass or the final stellar mass of the galaxy. However, the redshift at which the SFR peaks increases rapidly with halo mass. (3) More than half of the stars in the present-day Universe were formed in halos with $10^{11.1}\msunh < M_h < 10^{12.3}\msunh$ in the redshift range $0.4 < z < 1.9$. (4) ... [abridged]

With galaxy groups constructed from the Sloan Digital Sky Survey (SDSS), we analyze the expected galaxy-galaxy lensing signals around satellite galaxies residing in different host haloes and located at different halo-centric distances. We use Markov Chain Monte Carlo (MCMC) method to explore the potential constraints on the mass and density profile of subhaloes associated with satellite galaxies from SDSS-like surveys and surveys similar to the Large Synoptic Survey Telescope (LSST). Our results show that for SDSS-like surveys, we can only set a loose constraint on the mean mass of subhaloes. With LSST-like surveys, however, both the mean mass and the density profile of subhaloes can be well constrained.

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.

We simultaneously constrain cosmology and galaxy bias using measurements of galaxy abundances, galaxy clustering and galaxy-galaxy lensing taken from the Sloan Digital Sky Survey. We use the conditional luminosity function (which describes the halo occupation statistics as function of galaxy luminosity) combined with the halo model (which describes the non-linear matter field in terms of its halo building blocks) to describe the galaxy-dark matter connection. We explicitly account for residual redshift space distortions in the projected galaxy-galaxy correlation functions, and marginalize over uncertainties in the scale dependence of the halo bias and the detailed structure of dark matter haloes. Under the assumption of a spatially flat, vanilla {\Lambda}CDM cosmology, we focus on constraining the matter density, {\Omega}m, and the normalization of the matter power spectrum, {\sigma}8, and we adopt WMAP7 priors for the spectral index, the Hubble parameter, and the baryon density. We obtain that \Omegam = 0.278_{-0.026}^{+0.023} and {\sigma}8 = 0.763_{-0.049}^{+0.064} (95% CL). These results are robust to uncertainties in the radial number density distribution of satellite galaxies, while allowing for non-Poisson satellite occupation distributions results in a slightly lower value for {\sigma}8 (0.744_{-0.047}^{+0.056}). These constraints are in excellent agreement (at the 1{\sigma} level) with the cosmic microwave background constraints from WMAP. This demonstrates that the use of a realistic and accurate model for galaxy bias, down to the smallest non-linear scales currently observed in galaxy surveys, leads to results perfectly consistent with the vanilla {\Lambda}CDM cosmology.

We quantify the accuracy with which the cosmological parameters characterizing the energy density of matter (\Omega_m), the amplitude of the power spectrum of matter fluctuations (\sigma_8), the energy density of neutrinos (\Omega_{\nu}) and the dark energy equation of state (w_0) can be constrained using data from large galaxy redshift surveys. We advocate a joint analysis of the abundance of galaxies, galaxy clustering, and the galaxy-galaxy weak lensing signal in order to simultaneously constrain the halo occupation statistics (i.e., galaxy bias) and the cosmological parameters of interest. We parameterize the halo occupation distribution of galaxies in terms of the conditional luminosity function and use the analytical framework of the halo model described in our companion paper (van den Bosch et al. 2012), to predict the relevant observables. By performing a Fisher matrix analysis, we show that a joint analysis of these observables, even with the precision with which they are currently measured from the Sloan Digital Sky Survey, can be used to obtain tight constraints on the cosmological parameters, fully marginalized over uncertainties in galaxy bias. We demonstrate that the cosmological constraints from such an analysis are nearly uncorrelated with the halo occupation distribution constraints, thus, minimizing the systematic impact of any imperfections in modeling the halo occupation statistics on the cosmological constraints. In fact, we demonstrate that the constraints from such an analysis are both complementary to and competitive with existing constraints on these parameters from a number of other techniques, such as cluster abundances, cosmic shear and/or baryon acoustic oscillations, thus paving the way to test the concordance cosmological model.

We present a new method that simultaneously solves for cosmology and galaxy bias on non-linear scales. The method uses the halo model to analytically describe the (non-linear) matter distribution, and the conditional luminosity function (CLF) to specify the halo occupation statistics. For a given choice of cosmological parameters, this model can be used to predict the galaxy luminosity function, as well as the two-point correlation functions of galaxies, and the galaxy-galaxy lensing signal, both as function of scale and luminosity. In this paper, the first in a series, we present the detailed, analytical model, which we test against mock galaxy redshift surveys constructed from high-resolution numerical $N$-body simulations. We demonstrate that our model, which includes scale-dependence of the halo bias and a proper treatment of halo exclusion, reproduces the 3-dimensional galaxy-galaxy correlation and the galaxy-matter cross-correlation (which can be projected to predict the observables) with an accuracy better than 10 (in most cases 5) percent. Ignoring either of these effects, as is often done, results in systematic errors that easily exceed 40 percent on scales of $\sim 1 h^{-1}\Mpc$, where the data is typically most accurate. Finally, since the projected correlation functions of galaxies are never obtained by integrating the redshift space correlation function along the line-of-sight out to infinity, simply because the data only cover a finite volume, they are still affected by residual redshift space distortions (RRSDs). Ignoring these, as done in numerous studies in the past, results in systematic errors that easily exceed 20 perent on large scales ($r_\rmp \gta 10 h^{-1}\Mpc$). We show that it is fairly straightforward to correct for these RRSDs, to an accuracy better than $\sim 2$ percent, using a mildly modified version of the linear Kaiser formalism.

We present a method that uses observations of galaxies to simultaneously constrain cosmological parameters and the galaxy-dark matter connection (aka halo occupation statistics). The latter describes how galaxies are distributed over dark matter haloes, and is an imprint of the poorly understood physics of galaxy formation. A generic problem of using galaxies to constrain cosmology is that galaxies are a biased tracer of the mass distribution, and this bias is generally unknown. The great advantage of simultaneously constraining cosmology and halo occupation statistics is that this effectively allows cosmological constraints marginalized over the uncertainties regarding galaxy bias. Not only that, it also yields constraints on the galaxy-dark matter connection, this time properly marginalized over cosmology, which is of great value to inform theoretical models of galaxy formation. We use a combination of the analytical halo model and the conditional luminosity function to describe the galaxy-dark matter connection, which we use to model the abundance, clustering and galaxy-galaxy lensing properties of the galaxy population. We use a Fisher matrix analysis to gauge the complementarity of these different observables, and present some preliminary results from an analysis based on data from the Sloan Digital Sky Survey. Our results are complementary to and perfectly consistent with the results from the 7 year data release of the WMAP mission, strengthening the case for a true 'concordance' cosmology.

We present a new model to describe the galaxy-dark matter connection across cosmic time, which unlike the popular subhalo abundance matching technique is self-consistent in that it takes account of the facts that (i) subhalos are accreted at different times, and (ii) the properties of satellite galaxies may evolve after accretion. Using observations of galaxy stellar mass functions out to $z \sim 4$, the conditional stellar mass function at $z\sim 0.1$ obtained from SDSS galaxy group catalogues, and the two-point correlation function (2PCF) of galaxies at $z \sim 0.1$ as function of stellar mass, we constrain the relation between galaxies and dark matter halos over the entire cosmic history from $z \sim 4$ to the present. This relation is then used to predict the median assembly histories of different stellar mass components within dark matter halos (central galaxies, satellite galaxies, and halo stars). We also make predictions for the 2PCFs of high-$z$ galaxies as function of stellar mass. Our main findings are the following: (i) Our model reasonably fits all data within the observational uncertainties, indicating that the $\Lambda$CDM concordance cosmology is consistent with a wide variety of data regarding the galaxy population across cosmic time. (ii) ... [abridged]

We conduct Bayesian model inferences from the observed K-band luminosity function of galaxies in the local Universe, using the semi-analytic model (SAM) of galaxy formation introduced in Lu et al (2011). The prior distributions for the 14 free parameters include a large range of possible models. We find that some of the free parameters, e.g. the characteristic scales for quenching star formation in both high-mass and low-mass halos, are already tightly constrained by the single data set. The posterior distribution includes the model parameters adopted in other SAMs. By marginalising over the posterior distribution, we make predictions that include the full inferential uncertainties for the colour-magnitude relation, the Tully-Fisher relation, the conditional stellar mass function of galaxies in halos of different masses, the HI mass function, the redshift evolution of the stellar mass function of galaxies, and the global star formation history. Using posterior predictive checking with the available observational results, we find that the model family (i) predicts a Tully-Fisher relation that is curved; (ii) significantly over predicts the satellite fraction; (iii) vastly over predicts the HI mass function; (iv) predicts high-z stellar mass functions that have too many low mass galaxies and too few high mass ones. and (v) predicts a redshift evolution of the stellar mass density and the star formation history that are in moderate disagreement. These results suggest that some important processes are still missing in the current model family and we discuss a number of possible solutions to solve the discrepancies, such as interactions between galaxies and dark matter halos, tidal stripping, the bimodal accretion of gas, preheating, and a redshift-dependent initial mass function.

We use the ROSAT all sky survey X-ray cluster catalogs and the optical SDSS DR7 galaxy and group catalogs to cross-identify X-ray clusters with their optical counterparts, resulting in a sample of 201 X-ray clusters in the sky coverage of SDSS DR7. We investigate various correlations between the optical and X-ray properties of these X-ray clusters, and find that the following optical properties are correlated with the X-ray luminosity: the central galaxy luminosity, the central galaxy mass, the characteristic group luminosity ($\propto \Lx^{0.43}$), the group stellar mass ($\propto \Lx^{0.46}$), with typical 1-$\sigma$ scatter of $\sim 0.67$ in $\log \Lx$. Using the observed number distribution of X-ray clusters, we obtain an unbiased scaling relation between the X-ray luminosity, the central galaxy stellar mass and the characteristic satellite stellar mass as ${\log L_X} = -0.26 + 2.90 [\log (M_{\ast, c} + 0.26 M_{\rm sat}) -12.0]$ (and in terms of luminosities, as ${\log L_X} = -0.15 + 2.38 [\log (L_{c} + 0.72 L_{\rm sat}) -12.0]$). We find that the systematic difference between different halo mass estimations, e.g., using the ranking of characteristic group stellar mass or using the X-ray luminosity scaling relation can be used to constrain cosmology. Comparing the properties of groups of similar stellar mass (or optical luminosities) and redshift that are X-ray luminous or under-luminous, we find that X-ray luminous groups have more faint satellite galaxies and higher red fraction in their satellites. The cross-identified X-ray clusters together with their optical properties are provided in Appendix B.

[abridge]Cosmic velocity and tidal fields are important for the understanding of the cosmic web and the environments of galaxies, and can also be used to constrain cosmology. In this paper, we reconstruct these two fields in SDSS volume from dark matter halos represented by galaxy groups. Detailed mock catalogues are used to test the reliability of our method against uncertainties arising from redshift distortions, survey boundaries, and false identifications of groups by our group finder. We find that both the velocity and tidal fields, smoothed on a scale of ~2Mpc/h, can be reliably reconstructed in the inner region (~66%) of the survey volume. The reconstructed tidal field is used to split the cosmic web into clusters, filaments, sheets, and voids, depending on the sign of the eigenvalues of tidal tensor. The reconstructed velocity field nicely shows how the flows are diverging from the centers of voids, and converging onto clusters, while sheets and filaments have flows that are convergent along one and two directions, respectively. We use the reconstructed velocity field and the Zel'dovich approximation to predict the mass density field in the SDSS volume as function of redshift, and find that the mass distribution closely follows the galaxy distribution even on small scales. We find a large-scale bulk flow of about 117km/s in a very large volume, equivalent to a sphere with a radius of ~170Mpc/h, which seems to be produced by the massive structures associated with the SDSS Great Wall. Finally, we discuss potential applications of our reconstruction to study the environmental effects of galaxy formation, to generate initial conditions for simulations of the local Universe, and to constrain cosmological models. The velocity, tidal and density fields in the SDSS volume, specified on a Cartesian grid with a spatial resolution of ~700kpc/h, are available from the authors upon request.

We believe that a wide range of physical processes conspire to shape the observed galaxy population but we remain unsure of their detailed interactions. The semi-analytic model (SAM) of galaxy formation uses multi-dimensional parameterisations of the physical processes of galaxy formation and provides a tool to constrain these underlying physical interactions. Because of the high dimensionality, the parametric problem of galaxy formation may be profitably tackled with a Bayesian-inference based approach, which allows one to constrain theory with data in a statistically rigorous way. In this paper we develop a SAM in the framework of Bayesian inference. We show that, with a parallel implementation of an advanced Markov-Chain Monte-Carlo algorithm, it is now possible to rigorously sample the posterior distribution of the high-dimensional parameter space of typical SAMs. As an example, we characterise galaxy formation in the current $\Lambda$CDM cosmology using the stellar mass function of galaxies as an observational constraint. We find that the posterior probability distribution is both topologically complex and degenerate in some important model parameters, suggesting that thorough explorations of the parameter space are needed to understand the models. We also demonstrate that because of the model degeneracy, adopting a narrow prior strongly restricts the model. Therefore, the inferences based on SAMs are conditional to the model adopted. Using synthetic data to mimic systematic errors in the stellar mass function, we demonstrate that an accurate observational error model is essential to meaningful inference.

We study the behaviour of multiple radiative cooling algorithms implemented in seven Semi-Analytic Models (SAMs) of galaxy formation, including a new model we propose in this paper. We use versions of the models without feedback and apply them to dark matter haloes growing in a cosmological context, which have final masses that range from 10^{11}Msun to 10^{14}Msun. First, using simplified smoothly-growing halo models, we demonstrate that the different algorithms predict cooling rates and final cold gas masses that differ by a factor of ~5 for massive haloes (>10^{12}Msun). The algorithms are in better agreement for less massive haloes because they cool efficiently and, therefore, their cooling rates are largely limited by the halo accretion rate. However, for less massive haloes, all the SAMs predict less cooling than corresponding 1D hydrodynamic models. Second, we study the gas accretion history of the central galaxies of dark matter haloes using merger trees. The inclusion of mergers alters the cooling history of haloes by locking up gas in galaxies within small haloes at early times. For realistic halo models, the dispersion in the cold gas mass predicted by the algorithms is 0.5 dex for high mass haloes and 0.1 dex for low mass haloes, while the dispersion in the accretion rate is about two times larger. Comparing to cosmological SPH simulations, we find that most SAMs systematically under-predict the gas accretion rates for low-mass haloes but over-predict the gas accretion rates for massive haloes. Although the models all include both "rapid" and "slow" mode accretion, the transition between the two accretion modes varies between models and also differs from the simulations. Finally, we construct a new model that explicitly incorporates cold halo gas to illustrate that such a class of models can better match the results from cosmological hydrodynamic simulations.

An analytical model is developed for the mass function of cold dark matter subhalos at the time of accretion and for the distribution of their accretion times. Our model is based on the model of Zhao et al. (2009) for the median assembly histories of dark matter halos, combined with a simple log-normal distribution to describe the scatter in the main-branch mass at a given time for halos of the same final mass. Our model is simple, and can be used to predict the un-evolved subhalo mass function, the mass function of subhalos accreted at a given time, the accretion-time distribution of subhalos of a given initial mass, and the frequency of major mergers as a function of time. We test our model using high-resolution cosmological $N$-body simulations, and find that our model predictions match the simulation results remarkably well. Finally, we discuss the implications of our model for the evolution of subhalos in their hosts and for the construction of a self-consistent model to link galaxies and dark matter halos at different cosmic times.

It is generally assumed that the central galaxy in a dark matter halo, that is, the galaxy with the lowest specific potential energy, is also the brightest halo galaxy (BHG), and that it resides at rest at the centre of the dark matter potential well. This central galaxy paradigm (CGP) is an essential assumption made in various fields of astronomical research. In this paper we test the validity of the CGP using a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. For each group we compute two statistics, ${\cal R}$ and ${\cal S}$, which quantify the offsets of the line-of-sight velocities and projected positions of brightest group galaxies relative to the other group members. By comparing the cumulative distributions of $|{\cal R}|$ and $|{\cal S}|$ to those obtained from detailed mock group catalogues, we rule out the null-hypothesis that the CGP is correct. Rather, the data indicate that in a non-zero fraction $f_{\rm BNC}(M)$ of all haloes of mass $M$ the BHG is not the central galaxy, but instead, a satellite galaxy. In particular, we find that $f_{\rm BNC}$ increases from $\sim 0.25$ in low mass haloes ($10^{12} h^{-1} {\rm M_{\odot}} \leq M \lsim 2 \times 10^{13} h^{-1}{\rm M_{\odot}}$) to $\sim 0.4$ in massive haloes ($M \gsim 5 \times 10^{13} h^{-1} {\rm M_{\odot}}$). We show that these values of $f_{\rm BNC}$ are uncomfortably high compared to predictions from halo occupation statistics and from semi-analytical models of galaxy formation. We end by discussing various implications of a non-zero $f_{\rm BNC}(M)$, with an emphasis on the halo masses inferred from satellite kinematics.

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.

We investigate the potential of exploiting the Sunyaev-Zeldovich effect (SZE) to study the properties of hot gas in galaxy groups. It is shown that, with upcoming SZE surveys, one can stack SZE maps around galaxy groups of similar halo masses selected from large galaxy redshift surveys to study the hot gas in halos represented by galaxy groups. We use various models for the hot halo gas to study how the expected SZE signals are affected by gas fraction, equation of state, halo concentration, and cosmology. Comparing the model predictions with the sensitivities expected from the SPT, ACT and Planck surveys shows that a SPT-like survey can provide stringent constraints on the hot gas properties for halos with masses M ~> 10^{13} h^{-1}Msun. We also explore the idea of using the cross correlation between hot gas and galaxies of different luminosity to probe the hot gas in dark matter halos without identifying galaxy groups to represent dark halos. Our results show that, with a galaxy survey as large as the Sloan Digital Sky Survey and with the help of the conditional luminosity function (CLF) model, one can obtain stringent constraints on the hot gas properties in halos with masses down to 10^{13} h^{-1}Msun. Thus, the upcoming SZE surveys should provide a very promising avenue to probe the hot gas in relatively low-mass halos where the majority of L*-galaxies reside.

We use the kinematics of satellite galaxies that orbit around the central galaxy in a dark matter halo to infer the scaling relations between halo mass and central galaxy properties. Using galaxies from the Sloan Digital Sky Survey, we investigate the halo mass-luminosity relation (MLR) and the halo mass-stellar mass relation (MSR) of central galaxies. In particular, we focus on the dependence of these scaling relations on the colour of the central galaxy. We find that red central galaxies on average occupy more massive haloes than blue central galaxies of the same luminosity. However, at fixed stellar mass there is no appreciable difference in the average halo mass of red and blue centrals, especially for M* $\lsim$ 10^{10.5} h^{-2} Msun. This indicates that stellar mass is a better indicator of halo mass than luminosity. Nevertheless, we find that the scatter in halo masses at fixed stellar mass is non-negligible for both red and blue centrals. It increases as a function of stellar mass for red centrals but shows a fairly constant behaviour for blue centrals. We compare the scaling relations obtained in this paper with results from other independent studies of satellite kinematics, with results from a SDSS galaxy group catalog, from galaxy-galaxy weak lensing measurements, and from subhalo abundance matching studies. Overall, these different techniques yield MLRs and MSRs in fairly good agreement with each other (typically within a factor of two), indicating that we are converging on an accurate and reliable description of the galaxy-dark matter connection. We briefly discuss some of the remaining discrepancies among the various methods.

We compare the stellar masses of central and satellite galaxies predicted by three independent semianalytical models with observational results obtained from a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. In particular, we compare the stellar mass functions of centrals and satellites, the relation between total stellar mass and halo mass, and the conditional stellar mass functions, which specify the average number of galaxies of stellar mass M_* that reside in a halo of mass M_h. The semi-analytical models only predict the correct stellar masses of central galaxies within a limited mass range and all models fail to reproduce the sharp decline of stellar mass with decreasing halo mass observed at the low mass end. In addition, all models over-predict the number of satellite galaxies by roughly a factor of two. The predicted stellar mass in satellite galaxies can be made to match the data by assuming that a significant fraction of satellite galaxies are tidally stripped and disrupted, giving rise to a population of intra-cluster stars in their host halos. However, the amount of intra-cluster stars thus predicted is too large compared to observation. This suggests that current galaxy formation models still have serious problems in modeling star formation in low-mass halos.

Using a large SDSS galaxy group catalogue, we study how the stellar ages and metallicities of central and satellite galaxies depend on stellar mass and halo mass. We find that satellites are older and metal-richer than centrals of the same stellar mass. In addition, the slopes of the age-stellar mass and metallicity-stellar mass relations are found to become shallower in denser environments. This is due to the fact that the average age and metallicity of low mass satellite galaxies increase with the mass of the halo in which they reside. A comparison with the semi-analytical model of Wang et al. (2008) shows that it succesfully reproduces the fact that satellites are older than centrals of the same stellar mass and that the age difference increases with the halo mass of the satellite. This is a consequence of strangulation, which leaves the stellar populations of satellites to evolve passively, while the prolonged star formation activity of centrals keeps their average ages younger. The resulting age offset is larger in more massive environments because their satellites were accreted earlier. The model fails, however, in reproducing the halo mass dependence of the metallicities of low mass satellites, yields metallicity-stellar mass and age-stellar mass relations that are too shallow, and predicts that satellite galaxies have the same metallicities as centrals of the same stellar mass, in disagreement with the data. We argue that these discrepancies are likely to indicate the need to (i) modify the recipes of both supernova feedback and AGN feedback, (ii) use a more realistic description of strangulation, and (iii) include a proper treatment of the tidal stripping, heating and destruction of satellite galaxies. [Abridged]

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 use the "Millennium Simulation" to study the mass function of accreted sub-halos during merger events in the dark halo assembly history. Our study includes three kinds of sub-halo mergers: (1) mergers that happen to the main progenitor of dark halos; (2) mergers that happen on the entire merging history tree of dark halos; and (3) mergers that leave identifiable sub-halos in present-day dark halos. We estimate the unevolved sub-halo mass functions (USMFs), for which sub-halo masses are measured at the times of their accretion. For sub-halos that merge into the main branch of a present-day dark halo, their USMF can be well described by a universal functional form, in excellent agreement with previous results. The same conclusion can also be reached for the USMF of all progenitors that have merged to become sub-halos during the entire halo merging history. In both cases, the USMFs are also independent of the redshift of host halos. Due to tidal disruption, only a small fraction of the accreted halos survive as sub-halos identifiable in the present-day dark halos. In cluster-sized halos, about 30% of the survived sub-halos are sub-subhalos, and this fraction decreases with decreasing halo mass. For given halo and sub-halo masses, the accretion time has very broad distribution, but the survived sub-halos are all accreted quite recently.

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.

Using dark matter halos traced by galaxy groups selected from the Sloan Digital Sky Survey Data Release 4, we find that about 1/4 of the faint galaxies ($\rmag >-17.05$, hereafter dwarfs) that are the central galaxies in their own halo are not blue and star forming, as expected in standard models of galaxy formation, but are red. In contrast, this fraction is about 1/2 for dwarf satellite galaxies. Many red dwarf galaxies are physically associated with more massive halos. In total, about $\sim 45$% of red dwarf galaxies reside in massive halos as satellites, while another $\sim 25$% have a spatial distribution that is much more concentrated towards their nearest massive haloes than other dwarf galaxies. We use mock catalogs to show that the reddest population of non-satellite dwarf galaxies are distributed within about 3 times the virial radii of their nearest massive halos. We suggest that this population of dwarf galaxies are hosted by low-mass halos that have passed through their massive neighbors, and that the same environmental effects that cause satellite galaxies to become red are also responsible for the red colors of this population of galaxies. We do not find any significant radial dependence of the population of dwarf galaxies with the highest concentrations, suggesting that the mechanisms operating on these galaxies affect color more than structure. However, over 30% of dwarf galaxies are red and isolated and their origin remains unknown.

We use galaxy groups selected from the Sloan Digital Sky Survey (SDSS) together with mass models for individual groups to study the galaxy-galaxy lensing signals expected from galaxies of different luminosities and morphological types. We compare our model predictions with the observational results obtained from the SDSS by Mandelbaum et al. (2006) for the same samples of galaxies. The observational results are well reproduced in a $\Lambda$CDM model based on the WMAP 3-year data, but a $\Lambda$CDM model with higher $\sigma_8$, such as the one based on the WMAP 1-year data,significantly over-predicts the galaxy-galaxy lensing signal. We model, separately, the contributions to the galaxy-galaxy lensing signals from different galaxies: central versus satellite, early-type versus late-type, and galaxies in halos of different masses. We also examine how the predicted galaxy-galaxy lensing signal depends on the shape, density profile, and the location of the central galaxy with respect to its host halo.

Using a large galaxy group catalogue constructed from the Sloan Digital Sky Survey Data Release 4 (SDSS DR4) with an adaptive halo-based group finder, we investigate the luminosity and stellar mass functions for different populations of galaxies (central versus satellite; red versus blue; and galaxies in groups of different masses) and for groups themselves. The conditional stellar mass function (CSMF), which describes the stellar distribution of galaxies in halos of a given mass for central and satellite galaxies can be well modeled with a log-normal distribution and a modified Schechter form, respectively. On average, there are about 3 times as many central galaxies as satellites. Among the satellite population, there are in general more red galaxies than blue ones. For the central population, the luminosity function is dominated by red galaxies at the massive end, and by blue galaxies at the low mass end. At the very low-mass end ($M_\ast \la 10^9 h^{-2}\Msun$), however, there is a marked increase in the number of red centrals. We speculate that these galaxies are located close to large halos so that their star formation is truncated by the large-scale environments. The stellar-mass function of galaxy groups is well described by a double power law, with a characteristic stellar mass at $\sim 4\times 10^{10}h^{-2}\Msun$. Finally, we use the observed stellar mass function of central galaxies to constrain the stellar mass - halo mass relation for low mass halos, and obtain $M_{\ast, c}\propto M_h^{4.9}$ for $M_h \ll 10^{11} \msunh$.

Using a representative sample of 911 central galaxies (CENs) from the SDSS DR4 group catalogue, we study how the structure of the most massive members in groups and clusters depend on (1) galaxy stellar mass (Mstar), (2) dark matter halo mass of the host group (Mhalo), and (3) their halo-centric position. We establish and thoroughly test a GALFIT-based pipeline to fit 2D Sersic models to SDSS data. We find that the fitting results are most sensitive to the background sky level determination and strongly recommend using the SDSS global value. We find that uncertainties in the background translate into a strong covariance between the total magnitude, half-light size (r50), and Sersic index (n), especially for bright/massive galaxies. We find that n depends strongly on Mstar for CENs, but only weakly or not at all on Mhalo. Less (more) massive CENs tend to be disk (spheroid)-like over the full Mhalo range. Likewise, there is a clear r50-Mstar relation for CENs, with separate slopes for disks and spheroids. When comparing CENs with satellite galaxies (SATs), we find that low mass (<10e10.75 Msun/h^2) SATs have larger median n than CENs of similar Mstar. Low mass, late-type SATs have moderately smaller r50 than late-type CENs of the same Mstar. However, we find no size differences between spheroid-like CENs and SATs, and no structural differences between CENs and SATs matched in both mass and colour. The similarity of massive SATs and CENs shows that this distinction has no significant impact on the structure of spheroids. We conclude that Mstar is the most fundamental property determining the basic structure of a galaxy. The lack of a clear n-Mhalo relation rules out a distinct group mass for producing spheroids, and the responsible morphological transformation processes must occur at the centres of groups spanning a wide range of masses. (abridged)

We investigate the correlation of star formation quenching with internal galaxy properties and large scale environment (halo mass) in empirical data and theoretical models. We make use of the halo-based Group Catalog of Yang and collaborators, which is based on the Sloan Digital Sky Survey. Data from the Galaxy Evolution Explorer (GALEX) are also used to extract the recent star formation rate. In order to investigate the environmental effects, we examine the properties of "central" and "satellite" galaxies separately. For central galaxies, we are unable to conclude whether star formation quenching is primarily connected with halo mass or stellar mass, because these two quantities are themselves strongly correlated. For satellite galaxies, a nearly equally strong dependence on halo mass and stellar mass is seen. We make the same comparison for five different semi-analytic models based on three independently developed codes. We find that the models with AGN feedback reproduce reasonably well the dependence of the fraction of central red and passive galaxies on halo mass and stellar mass. However, for satellite galaxies, the same models badly overproduce the fraction of red/passive galaxies and do not reproduce the empirical trends with stellar mass or halo mass. This {\em satellite overquenching problem} is caused by the too-rapid stripping of the satellites' hot gas halos, which leads to rapid strangulation of star formation.

Using the SDSS DR4 group catalogue of Yang et al. (2007), we investigate sizes, concentrations, colour gradients and surface brightness profiles of central and satellite galaxies. We compare central and satellite galaxies at fixed stellar mass, in order to disentangle environmental from stellar mass dependencies. Early and late type galaxies are defined according to concentration. We find that at fixed stellar mass, late type satellite galaxies have smaller radii and larger concentrations than late type central galaxies. No such differences are found for early-type galaxies. We have also constructed surface brightness and colour profiles for the central and satellite galaxies in our sample. We find that late-type satellite galaxies have a lower surface brightness and redder colours than late-type central galaxies. We show that all observed differences between satellite and central galaxies can be explained by a simple fading model, in which the star formation in the disk decreases over timescales of 2-3 Gyr after a galaxy becomes a satellite. Processes that induce strong morphological changes (e.g. harassment) and processes that strip the galaxy of its entire ISM need not to be invoked in order to explain the environmental dependencies we find.

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.

In the standard paradigm, satellite galaxies are believed to be associated with the population of dark matter subhalos. In this paper, we use the conditional stellar mass functions of {\it satellite galaxies} obtained from a large galaxy group catalogue together with models of the subhalo mass functions to explore the fraction and fate of stripped stars from satellites in galaxy groups and clusters of different masses. The majority of the stripped stars in massive halos are predicted to end up as intra-cluster stars, and the predicted amounts of the intra-cluster component as a function of the velocity dispersion of galaxy system match well the observational results obtained by Gonzalez et al. (2007). The fraction of the mass in the stripped stars to that remain bound in the central and satellite galaxies is the highest ($\sim 40%$ of the total stellar mass) in halos with masses $M_h\sim 10^{14}\msunh$. If all these stars end up in the intra-cluster component (Max), or maximum of them are accreted into the central galaxy (Min), then we can predict that a maximum $\sim 19%$ and a minimum $\sim 5%$ of the total stars in the whole universe are in terms of the diffused intra-cluster component. In the former case, in massive halos with $M_h \sim 10^{15} \msunh$, the stellar mass of the intra-cluster component is roughly 6 times as large as that of the central galaxy. This factor decreases to $\sim 2$, 1 and 0.1 in halos with $M_h \sim 10^{14}$, $10^{13}$, and $10^{12} \msunh$, respectively. The total amount of stars stripped from satellite galaxies is insufficient to build up the central galaxies in halos with masses $\la 10^{12.5}\msunh$, and so the quenching of star formation must occur in halos with higher masses. Abridged.

We demonstrate that the feedback from stellar bulges can play an essential role in shaping the halo gas of galaxies with substantial bulge components by conducting 1-D hydrodynamical simulations. The feedback model we consider consists of two distinct phases: 1) an early starburst during the bulge formation and 2) a subsequent long-lasting mass and energy injection from stellar winds of low-mass stars and Type Ia SNe. An energetic outward blastwave is initiated by the starburst and is maintained and enhanced by the long-lasting stellar feedback. For a MW-like galactic bulge, this blastwave sweeps up the halo gas in the proto-galaxy and heats up the surrounding medium to a scale much beyond the virial radius of the halo, thus the accretion of the halo hot gas can be completely stopped. In addition, the long-lasting feedback in the later phase powers a galactic bulge wind that is reverse-shocked at a large radius in the presence of surrounding intergalactic medium and hence maintains a hot gaseous halo. As the mass and energy injection decreases with time, the feedback evolves to a subsonic and quasi-stable outflow, which is enough to prevent halo gas from cooling. The two phases of the feedback thus re-enforce each-other's impact on the gas dynamics. The simulation results demonstrate that the stellar bulge feedback may provide a plausible solution to the long-standing problems in understanding the MW type galaxies, such as the "missing stellar feedback" problem and the "over-cooling" problem. The simulations also show that the properties of the hot gas in the subsonic outflow state depend sensitively on the environment and the formation history of the bulge. This dependence and variance may explain the large dispersion in the X-ray to B-band luminosity ratio of the low $L_X/L_B$ Es.

We use a modified version of the halo-based group finder developed by Yang et al. to select galaxy groups from the Sloan Digital Sky Survey (SDSS DR4). In the first step, a combination of two methods is used to identify the centers of potential groups and to estimate their characteristic luminosity. Using an iterative approach, the adaptive group finder then uses the average mass-to-light ratios of groups, obtained from the previous iteration, to assign a tentative mass to each group. This mass is then 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. Finally, each individual group is assigned two different halo masses: one based on its characteristic luminosity, and the other based on its characteristic stellar mass. Applying the group finder to the SDSS DR4, we obtain 301237 groups in a broad dynamic range, including systems of isolated galaxies. We use detailed mock galaxy catalogues constructed for the SDSS DR4 to test the performance of our group finder in terms of completeness of true members, contamination by interlopers, and accuracy of the assigned masses. This paper is the first in a series and focuses on the selection procedure, tests of the reliability of the group finder, and the basic properties of the group catalogue (e.g. the mass-to-light ratios, the halo mass to stellar mass ratios, etc.). The group catalogues including the membership of the groups are available at http://gax.shao.ac.cn/data/Group.html and http://www.astro.umass.edu/~xhyang/Group.html

Galaxy clustering and galaxy-galaxy lensing probe the connection between galaxies and their dark matter haloes in complementary ways. On one hand, the halo occupation statistics inferred from the observed clustering properties of galaxies are degenerate with the adopted cosmology. Consequently, different cosmologies imply different mass-to-light ratios for dark matter haloes. On the other hand, galaxy-galaxy lensing yields direct constraints on the actual mass-to-light ratios and it can be used to break this degeneracy, and thus to constrain cosmological parameters. In this paper we establish the link between galaxy luminosity and dark matter halo mass using the conditional luminosity function (CLF). We constrain the CLF parameters using the galaxy luminosity function and the luminosity dependence of the correlation lengths of galaxies. The resulting CLF models are used to predict the galaxy-galaxy lensing signal. For a cosmology with $(\Omega_{\rm m},\sigma_8)=(0.238,0.734)$, the model accurately fits the galaxy-galaxy lensing data obtained from the SDSS. For a comparison cosmology with $(\Omega_{\rm m},\sigma_8)=(0.3,0.9)$, however, we can accurately fit the luminosity function and clustering properties of the galaxy population, but the model predicts mass-to-light ratios that are too high, resulting in a strong overprediction of the galaxy-galaxy lensing signal. We conclude that the combination of galaxy clustering and galaxy-galaxy lensing is a powerful probe of the galaxy-dark matter connection, with the potential to yield tight constraints on cosmological parameters. Since this method mainly probes the mass distribution on non-linear scales, it is complementary to constraints obtained from the galaxy power-spectrum, which mainly probes the large-scale (linear) matter distribution.

The kinematics of satellite galaxies reflect the masses of the extended dark matter haloes in which they orbit, and thus shed light on the mass-luminosity relation (MLR) of their corresponding central galaxies. In this paper we select a large sample of centrals and satellites from the Sloan Digital Sky Survey (SDSS) and measure the kinematics (velocity dispersions) of the satellite galaxies as a function of the $r$-band luminosity of the central galaxies. Using the analytical framework presented in Paper I, we use these data to infer {\it both} the mean and the scatter of the MLR of central galaxies, carefully taking account of selection effects and biases introduced by the stacking procedure. As expected, brighter centrals on average reside in more massive haloes. In addition, we find that the scatter in halo masses for centrals of a given luminosity, $\sigma_{\log M}$, also increases with increasing luminosity. As we demonstrate, this is consistent with $\sigma_{\log L}$, which reflects the scatter in the conditional probability function $P(L_c|M)$, being independent of halo mass. Our analysis of the satellite kinematics yields $\sigma_{\log L}=0.16\pm0.04$, in excellent agreement with constraints from clustering and group catalogues, and with predictions from a semi-analytical model of galaxy formation. We thus conclude that the amount of stochasticity in galaxy formation, which is characterized by $\sigma_{\log L}$, is well constrained, is independent of halo mass, and is in good agreement with current models of galaxy formation.

We use a SDSS galaxy group catalogue to study the dependence of galaxy activity on stellar mass, halo mass, and group hierarchy (centrals vs. satellites). We split our galaxy sample in star-forming galaxies, galaxies with optical AGN activity and radio sources. We find a smooth transition in halo mass as the activity of central galaxies changes from star formation to optical AGN activity to radio emission. Star-forming centrals preferentially reside in haloes with M<10^{12} Msun, central galaxies with optical-AGN activity typically inhabit haloes with M \sim 10^{13} Msun, and centrals emitting in the radio mainly reside in haloes more massive than 10^{14} Msun. Although this seems to suggest that the environment (halo mass) determines the type of activity of its central galaxy, we find a similar trend with stellar mass: central star formers typically have stellar masses below 10^{10} Msun, while optical-AGN hosts and central radio sources have characteristic stellar masses of 10^{10.8} Msun and 10^{11.6} Msun, respectively. Since more massive haloes typically host more massive centrals, it is unclear whether the activity of a central galaxy is causally connected to its stellar mass or to its halo mass. In general, satellite galaxies have their activity suppressed wrt central galaxies of the same stellar mass. At fixed stellar mass, we find that the activity of satellite galaxies depends only weakly on halo mass. In fact, for satellite galaxies the dependence of galaxy activity on halo mass is more than four times weaker than the dependence on stellar mass. As we discuss, all these results are consistent with a picture in which low mass haloes accrete cold gas, while massive haloes have coronae of hot gas that promote radio activity of their central galaxies. [Abridged]

We use a sample of galaxy groups selected from the SDSS DR 4 with an adaptive halo-based group finder to probe how the clustering strength of groups depends on their masses and colors. In particular, we determine the relative biases of groups of different masses, as well as that of groups with the same mass but with different colors. In agreement with previous studies, we find that more massive groups are more strongly clustered, and the inferred mass dependence of the halo bias is in good agreement with predictions for the $\Lambda$CDM cosmology. Regarding the color dependence, we find that groups with red centrals are more strongly clustered than groups of the same mass but with blue centrals. Similar results are obtained when the color of a group is defined to be the total color of its member galaxies. The color dependence is more prominent in less massive groups and becomes insignificant in groups with masses $\gta 10^{14}\msunh$. We construct a mock galaxy redshift survey constructed from the large Millenium simulation that is populated with galaxies according to the semi-analytical model of Croton et al. Applying our group finder to this mock survey, and analyzing the mock data in exactly the same way as the true data, we are able to accurately recover the intrinsic mass and color dependencies of the halo bias in the model. This suggests that our group finding algorithm and our method of assigning group masses do not induce spurious mass and/or color dependencies in the group-galaxy correlation function. The semi-analytical model reveals the same color dependence of the halo bias as we find in our group catalogue. In halos with $M\sim 10^{12}\msunh$, though, the strength of the color dependence is much stronger in the model than in the data.

Using the Sloan Digital Sky Survey (SDSS) galaxy group catalogue of Yang et al. (2007), we study the average colour and average concentration of satellite galaxies as function of (i) their stellar mass, (ii) their group mass, and (iii) their group-centric radius. We find that the colours and concentrations of satellite galaxies are (almost) completely determined by their stellar mass. In particular, at fixed stellar mass, the average colours and concentrations of satellite galaxies are independent of either halo mass or halo-centric radius. We find clear evidence for mass segregation of satellite galaxies in haloes of all masses, and argue that this explains why satellites at smaller halo-centric radii are somewhat redder and somewhat more concentrated. In addition, the weak colour and concentration dependence of satellite galaxies on halo mass is simply a reflection of the fact that more massive haloes host, on average, more massive satellites. Combining these results with the fact that satellite galaxies are, on average, redder and somewhat more concentrated than central galaxies of the same stellar mass, the following picture emerges: galaxies become redder and somewhat more concentrated once they fall into a bigger halo (i.e., once they become a satellite galaxy). This is a clear manifestation of environment dependence. However, there is no indication that the magnitude of the transformation (or its timescale) depends on environment; a galaxy undergoes a transition when it becomes a satellite, but it does not matter whether it becomes a satellite of a small (Milky Way sized) halo, or of a massive cluster. We discuss the implication of this `dearth' of environment dependence for the physical processes responsible for transforming satellite galaxies.

In this paper we use the ``Millennium Simulation'' to re-examine the mass assembly history of dark matter halos and the age dependence of halo clustering. We use eight different definitions of halo formation times to characterize the different aspects of the assembly history of a dark matter halo. We find that these formation times have different dependence on halo mass. While some formation times characterize well the hierarchical nature of halo formation, in the sense that more massive halos have later formation, the trend is reversed for other definitions of the formation time. In particular, the formation times that are likely to be related to the formation of galaxies in dark halos show strong trends of ``down-sizing'', in that lower-mass halos form later. We also investigate how the correlation amplitude of dark matter halos depends on the different formation times. We find that this dependence is quite strong for some definitions of formation time but weak or absent for other definitions. In particular, the correlation amplitude of halos of a given mass is almost independent of their last major merger time. For the definitions that are expected to be more related to the formation of galaxies in dark halos, a significant assembly bias is found only for halos less massive than M_*. We discuss our results in connection to the hierarchical assembly of dark matter halos, the ``archaeological down-sizing'' observed in the galaxy population, and the observed color-dependence of the clustering strength of galaxy groups and clusters.

We study the three-dimensional and projected shapes of galaxy groups in the Sloan Digital Sky Survey Data Release 4, and examine the alignment between the orientation of the central galaxy and the spatial distribution of satellite galaxies. The projected ellipticity of a group is measured using the moments of the discrete distribution of its member galaxies. We infer the three-dimensional and projected axis ratios of their dark matter haloes by comparing the measured ellipticity distributions with those obtained from Monte Carlo simulations of projected, triaxial dark matter haloes with different axis ratios. We find that the halo shape has a strong dependence on the halo mass. While the haloes of low-mass groups are nearly spherical, those of massive groups tend to be prolate. For groups containing at least four members, the statistical distribution of their measured ellipticities does not have a strong dependence on the colors of their central galaxies. Our analysis further shows that the average three-dimensional axis ratio for haloes with $12<{\rm log}[M/(h^{-1}M_{\odot})]\leq15$ is about $1:0.46:0.46$, resulting in a projected axis ratio of $\sim 0.77$. Our results for the alignment between the orientation of the central galaxy of a group and the distribution of their satellite galaxies are in broad agreement with those obtained by Yang et al. The distribution of satellite galaxies preferentially aligns with the major axis of the central galaxy, with a clear dependence on both halo mass and galaxy colors. (abridged)

We investigate various halo occupation statistics using a large galaxy group catalogue constructed from the SDSS DR4 with an adaptive halo-based group finder. The conditional luminosity function (CLF) is measured separately for all, red and blue galaxies, as well as in terms of central and satellite galaxies. The CLFs for central and satellite galaxies can be well modelled with a log-normal distribution and a modified Schechter form, respectively. About 85% of the central galaxies and about 80% of the satellite galaxies in halos with masses $M_h\ga 10^{14}\msunh$ are red galaxies. These numbers decrease to 50% and 40%, respectively, in halos with $M_h \sim 10^{12}\msunh$. For halos of a given mass, the distribution of the luminosities of central galaxies, $L_c$, has a dispersion of about 0.15 dex. The mean luminosity (stellar mass) of the central galaxies scales with halo mass as $L_c\propto M_h^{0.17}$ ($M_{*,c}\propto M_h^{0.22}$) for halos with masses $M\gg 10^{12.5}\msunh$, and both relations are significantly steeper for less massive halos. We also measure the luminosity (stellar mass) gap between the first and second brightest (most massive) member galaxies, $\log L_1 - \log L_2$ ($\log M_{*,1}-\log M_{*,2}$). These gap statistics, especially in halos with $M_h \la 10^{14.0}\msunh$, indicate that the luminosities of central galaxies are clearly distinct from those of their satellites. The fraction of fossil groups, defined as those groups with $\log L_1 - \log L_2\ge 0.8$, ranges from $\sim 2.5%$ for groups with $M_h\sim 10^{14}\msunh$ to 18-60% for groups with $M_h\sim 10^{13}\msunh$. Finally, we measure the fraction of satellites, which changes from $\sim 5.0%$ for galaxies with $\rmag\sim -22.0$ to $\sim40%$ for galaxies with $\rmag\sim -17.0$. (abridged)

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)

In the current paradigm, red sequence galaxies are believed to have formed as blue disk galaxies that subsequently had their star formation quenched. Since red-sequence galaxies typically have an early-type morphology, the transition from the blue to the red sequence also involves a morphological transformation. In this paper we study the impact of transformation mechanisms that operate only on satellite galaxies, such as strangulation, ram-pressure stripping and galaxy harassment. Using a large galaxy group catalogue constructed from the SDSS, we compare the colors and concentrations of satellites galaxies to those of central galaxies of the same stellar mass, adopting the hypothesis that the latter are the progenitors of the former. On average, satellites are redder and more concentrated than central galaxies of the same stellar mass. Central-satellite pairs that are matched in both stellar mass and color, however, show no average concentration difference, indicating that the transformation mechanisms affect color more than morphology. The color and concentration differences of matched central-satellite pairs are completely independent of the halo mass of the satellite galaxy, indicating that satellite-specific transformation mechanisms are equally efficient in haloes of all masses. This strongly favors strangulation as the main quenching mechanism for satellite galaxies. Finally, we determine the relative importance of satellite quenching for the build-up of the red sequence. We find that roughly 70 percent of red sequence satellite galaxies with a stellar mass of 10^9 Msun had their star formation quenched as satellites. This drops rapidly to zero with increasing stellar mass, indicating that a significant fraction of red satellites were already quenched before they became a satellite.

We investigate the incidence of major mergers creating >10e11 Msun galaxies in present-day groups and clusters more massive than 2.5e13 Msun. We identify 38 pairs of massive galaxies with mutual tidal interaction signatures selected from >5000 galaxies with >5e10 Msun that reside in 845 such groups. We fit the images of each galaxy pair as the line-of-sight projection of symmetric models and identify mergers by the presence of residual asymmetries around each progenitor, such as off-center isophotes, broad tidal tails, and dynamical friction wakes. At the resolution and sensitivity of the SDSS, such mergers are found in 16% of high-mass, galaxy-galaxy pairs with magnitude differences of <1.5 and <30 kpc projected separations. We find that 90% of these mergers have nearly equal-mass progenitors with red-sequence colors and centrally-concentrated morphologies, the hallmarks of dissipationless merger simulations. Mergers at group centers are more common than between 2 satellites, but both are morphologically indistinguishable and we tentatively conclude that the latter are likely located at the dynamical centers of recently accreted subhalos. The frequency of central and satellite merging diminishes with group mass consistent with dynamical friction expectations. Based on reasonable assumptions, the centers of these massive halos are growing in stellar mass by 1-9% per Gyr, on average. Compared to all LRG-LRG mergers, we find a 2-9 times higher rate for their merging when restricted to these dense environments. Our results imply that the massive end of the galaxy population continues to evolve hierarchically at a measurable level, and that the centers of massive groups are the preferred environment for merger-driven galaxy assembly. (abridged)

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.

(abridged) We study the mass assembly history (MAH) of dark matter haloes. We compare MAHs obtained using (i) merger trees constructed with the extended Press-Schechter (EPS) formalism, (ii) numerical simulations, and (iii) the Lagrangian perturbation code PINOCCHIO. We show that the PINOCCHIO MAHs are in excellent agreement with those obtained using numerical simulations. Using a suite of 55 PINOCCHIO simulations, with 256^3 particles each, we study the MAHs of 12,924 cold dark matter haloes in a \LambdaCDM concordance cosmology. We show that haloes less massive than the characteristic non-linear mass scale establish their potential wells much before they acquire most of their mass. The time when a halo reaches its maximum virial velocity roughly divides its mass assembly into two phases, a fast accretion phase which is dominated by major mergers, and a slow accretion phase dominated by minor mergers. Each halo experiences about 3 \pm 2 major mergers since its main progenitor had a mass equal to one percent of the final halo mass. This major merger statistic is found to be virtually independent of halo mass. However, the average redshift at which these major mergers occur, is strongly mass dependent, with more massive haloes experiencing their major mergers later.

We study the cross-correlation between galaxies of different luminosities and colors, using a sample selected from the SDSS Dr 4. Galaxies are divided into 6 samples according to luminosity, and each of these samples is divided into red and blue subsamples. Projected auto-correlation and cross-correlation is estimated for these subsample. At projected separations r_p > 1\mpch, all correlation functions are roughly parallel, although the correlation amplitude depends systematically on luminosity and color. On r_p < 1\mpch, the auto- and cross-correlation functions of red galaxies are significantly enhanced relative to the corresponding power laws obtained on larger scales. Such enhancement is absent for blue galaxies and in the cross-correlation between red and blue galaxies. We esimate the relative bias factor on scales r > 1\mpch for each subsample using its auto-correlation function and cross-correlation functions. The relative bias factors obtained from different methods are similar. For blue galaxies the luminosity-dependence of the relative bias is strong over the luminosity range probed (-23.0<M_r < -18.0),but for red galaxies the dependence is weaker and becomes insignificant for luminosities below L^*. To examine whether a significant stochastic/nonlinear component exists in the bias relation, we study the ratio R_ij= W_{ii}W_{jj}/W_{ij}^2, where W_{ij} is the projected correlation between subsample i and j. We find that the values of R_ij are all consistent with 1 for all-all, red-red and blue-blue samples, however significantly larger than 1 for red-blue samples. For faint red - faint blue samples the values of R_{ij} are as high as ~ 2 on small scales r_p < 1 \mpch and decrease with increasing r_p.

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.

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.

Using a samples of 61506 spiral galaxies selected from the SDSS DR2, we examine the luminosity function (LF) of spiral galaxies with different inclination angles. We find that the characteristic luminosity of the LF, $L^*$, decreases with increasing inclination, while the faint-end slope, $\alpha$, depends only weakly on it. The inclination-dependence of the LF is consistent with that expected from a simple model where the optical depth is proportional to the cosine of the inclination angle, and we use a likelihood method to recover both the coefficient in front of the cosine, $\gamma$, and the LF for galaxies viewed face-on. The value of $\gamma$ is quite independent of galaxy luminosity in a given band, and the values of $\gamma$ obtained in this way for the 5 SDSS bands give an extinction curve which is a power law of wavelength ($\tau\propto\lambda^{-n}$), with a power index $n=0.96\pm0.04$. Using the dust extinction for galaxies obtained by Kauffmann et al. (2003), we derive an `extinction-corrected' luminosity function for spiral galaxies. Dust extinction makes $M^*$ dimmer by about 0.5 magnitudes in the $z$-band, and about 1.2 magnitudes in the $u$- band. Since our analysis is based on a sample where selection effects are well under control, the dimming of edge-on galaxies relative to face-on galaxies is best explained by assuming that galaxy disks are optically thick in dust absorptions.

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.

We use the conditional luminosity function (CLF) and data from the 2dFGRS to constrain the average relation between light and mass in a LCDM cosmology with Omega_m=0.23 and sigma_8=0.74 (hereafter WMAP3 cosmology). Reproducing the observed luminosity dependence of the galaxy two-point correlation function results in average mass-to-light ratios that are about 35 percent lower than in a LCDM cosmology with Omega_m=0.3 and sigma_8=0.9 (hereafter WMAP1 cosmology). This removes an important problem with previous halo occupation models which had a tendency to predict cluster mass-to-light ratios that were too high. For the WMAP3 cosmology our model yields average mass-to-light ratios, central galaxy luminosities, halo occupation numbers, satellite fractions, and luminosity-gap statistics, that are all in excellent agreement with those obtained from a 2dFGRS group catalogue and from other independent studies. We also use our CLF model to compute the probability distribution P(M|L_cen), that a central galaxy of luminosity L_cen resides in a halo of mass M. We find this distribution to be much broader than what is typically assumed in HOD models, which has important implications for the interpretation of galaxy-galaxy lensing data. Finally, reproducing the luminosity dependence of the pairwise velocity dispersions in the 2dFGRS requires relatively low mass-to-light ratios for clusters and a satellite fraction that decreases strongly with increasing luminosity. This is only marginally consistent with our CLF constraints. We argue that a cosmology with parameters between those of the WMAP1 and WMAP3 cosmologies is likely to yield results with a higher level of consistency.

We consider the sensitivity of the circular-orbit adiabatic contraction approximation to the baryon condensation rate and the orbital structure of dark matter halos in the $\Lambda$CDM paradigm. Using one-dimensional hydrodynamic simulations including the dark matter halo mass accretion history and gas cooling, we demonstrate that the adiabatic approximation is approximately valid even though halos and disks may assemble simultaneously. We further demonstrate the validity of the simple approximation for $\Lambda$CDM halos with isotropic velocity distributions using three-dimensional N-body simulations. This result is easily understood: an isotropic velocity distribution in a cuspy halo requires more circular orbits than radial orbits. Conversely, the approximation is poor in the extreme case of a radial orbit halo. It overestimates the response a core dark matter halo, where radial orbit fraction is larger. Because no astronomically relevant models are dominated by low-angular momentum orbits in the vicinity of the disk and the growth time scale is never shorter than a dynamical time, we conclude that the adiabatic contraction approximation is useful in modeling the response of dark matter halos to the growth of a disk.

(abridged) We use a one-dimensional hydrodynamical code to investigate the effects of preheating on gas accretion and cooling in cold dark matter halos. In the absence of radiative cooling, preheating reduces the amount of gas that can be accreted into a halo, and the accreted gas fraction is determined by the ratio of the initial specific entropy of the gas to the virial entropy of the halo. In the presence of radiative cooling, preheating affects the gas fraction that can cool in two different ways. For small halos with masses <10^12Msun, preheating suppresses gas accretion, but most of the accreted gas can cool. For more massive halos, preheating not only reduces the amount of accreted gas, but also reduces the cooling efficiency. For both small and massive halos, gas cooling is delayed by preheating and in an inside-out fashion if the halo gas is assumed to be a single-phase medium. However, cooling can occur over a wider range of redshifts and radii, if we assume a multi-phase medium. As examples, two specific preheating cases are investigated. In the first case, the preheating entropy is assumed to be proportional to the virial entropy of the halo, as expected from AGN feedback. Such preheating effectively suppresses radiative cooling in halos with M>10^13Msun. We suggest that this may be the reason why the stellar mass function of galaxies breaks sharply at the massive end. Such preheating also helps create the hot diffused halos within which the "radio mode" feedback of AGNs can act effectively. In the second case, we assume the intergalactic medium is warm. Here the total amount of gas that can cool in a halo scales with halo mass as ~M^2, as would be required to match the observed stellar- and HI-mass functions in the current CDM model at the small mass end.

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 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 examine the relation between galaxy properties and environment in the SDSS DR2, quantifying environment in terms of the mass of the host halo, which is obtained with a new iterative group finder. We find that galaxy type fractions scale strongly and smoothly with halo mass, but, at fixed mass, not with luminosity. We compare these findings with the semi-analytical galaxy formation model of Croton et al. (2006). The discrepancies we find can be explained with an oversimplified implementation of strangulation, the neglect of tidal stripping, and shortcomings in the treatments of dust extinction and/or AGN feedback.

Successfully reproducing the galaxy luminosity function and the bimodality in the galaxy distribution requires a mechanism that can truncate star formation in massive haloes. Current models of galaxy formation consider two such truncation mechanisms: strangulation, which acts on satellite galaxies, and AGN feedback, which predominantly affects central galaxies. The efficiencies of these processes set the blue fraction of galaxies as function of galaxy luminosity and halo mass. In this paper we use a galaxy group catalogue extracted from the Sloan Digital Sky Survey (SDSS) to determine these fractions. To demonstrate the potential power of this data as a benchmark for galaxy formation models, we compare the results to the semi-analytical model for galaxy formation of Croton et al. (2006). Although this model accurately fits the global statistics of the galaxy population, as well as the shape of the conditional luminosity function, there are significant discrepancies when the blue fraction of galaxies as a function of mass and luminosity is compared between the observations and the model. In particular, the model predicts (i) too many faint satellite galaxies in massive haloes, (ii) a blue fraction of satellites that is much too low, and (iii) a blue fraction of centrals that is too high and with an inverted luminosity dependence. In the same order, we argue that these discrepancies owe to (i) the neglect of tidal stripping in the semi-analytical model, (ii) the oversimplified treatment of strangulation, and (iii) improper modeling of dust extinction and/or AGN feedback. The data presented here will prove useful to test and calibrate future models of galaxy formation and in particular to discriminate between various models for AGN feedback and other star formation truncation mechanisms.

We use the ROSAT All Sky Survey (RASS) to study the soft X-ray properties of a homogeneous sample of 46,420 quasars selected from the third data release of the Sloan Digital Sky Survey (SDSS). Optical luminosities, both at rest-frame 2500\AA ($L_{2500}$) and in [OIII] ($L_{[\rm{OIII}]}$) span more than three orders of magnitude, while redshifts range over $0.1<z<5.4$. We detect 3366 quasars directly in the observed 0.1--2.4 keV band. Sub-samples of radio-loud and radio-quiet objects (RLQs and RQQs) are obtained by cross-matching with the FIRST catalogue. We study the distribution of X-ray luminosity as a function of optical luminosity, redshift and radio power using both individual detections and stacks of complete sets of similar quasars. At every optical luminosity and redshift $\log L_{2\kev}$ is, to a good approximation, normally distributed with dispersion $\sim 0.40$, at least brightwards of the median X-ray luminosity. This median X-ray luminosity of quasars is a power law of optical luminosity with index $\sim 0.53$ for $L_{2500}$ and $\sim 0.30$ for $L_{[\rm{OIII}]}$. RLQs are systematically brighter than RQQs by about a factor of 2 at given optical luminosity. The zero-points of these relations increase systematically with redshift, possibly in different ways for RLQs and RQQs. Evolution is particularly strong at low redshift and if the optical luminosity is characterised by $L_{[\rm{OIII}]}$. At low redshift and at given $L_{[\rm{OIII}]}$ the soft X-ray emission from type II AGN is more than 100 times weaker than that from type I AGN.

We discuss an analytic approach for modeling structure formation in sheets, filaments and knots. This is accomplished by combining models of triaxial collapse with the excursion set approach: sheets are defined as objects which have collapsed along only one axis, filaments have collapsed along two axes, and halos are objects in which triaxial collapse is complete. In the simplest version of this approach, which we develop here, large scale structure shows a clear hierarchy of morphologies: the mass in large-scale sheets is partitioned up among lower mass filaments, which themselves are made-up of still lower mass halos. Our approach provides analytic estimates of the mass fraction in sheets, filaments and halos, and its evolution, for any background cosmological model and any initial fluctuation spectrum. In the currently popular $\Lambda$CDM model, our analysis suggests that more than 99% of the cosmic mass is in sheets, and 72% in filaments, with mass larger than $10^{10} M_{\odot}$ at the present time. For halos, this number is only 46%. Our approach also provides analytic estimates of how halo abundances at any given time correlate with the morphology of the surrounding large-scale structure, and how halo evolution correlates with the morphology of large scale structure.

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.

N-body simulations predict that CDM halo-assembly occurs in two phases: 1) a fast accretion phase with a rapidly deepening potential well; and 2) a slow accretion phase characterised by a gentle addition of mass to the outer halo with little change in the inner potential well. We demonstrate, using one-dimensional simulations, that this two-phase accretion leads to CDM halos of the NFW form and provides physical insight into the properties of the mass accretion history that influence the final profile. Assuming that the velocities of CDM particles are effectively isotropised by fluctuations in the gravitational potential during the fast accretion phase, we show that gravitational collapse in this phase leads to an inner profile rho(r) ~ r^{-1}. Slow accretion onto an established potential well leads to an outer profile with rho(r) ~ r^{-3}. The concentration of a halo is determined by the fraction of mass that is accreted during the fast accretion phase. Using an ensemble of realistic mass accretion histories, we show that the model predictions of the dependence of halo concentration on halo formation time, and hence the dependence of halo concentration on halo mass, and the distribution of halo concentrations all match those found in cosmological N-body simulations. Using a simple analytic model that captures much of the important physics we show that the inner r^{-1} profile of CDM halos is a natural result of hierarchical mass assembly with a initial phase of rapid accretion.

We study the dependence of the cross-correlation between galaxies and galaxy groups on group properties. Confirming previous results, we find that the correlation strength is stronger for more massive groups, in good agreement with the expected mass dependence of halo bias. We also find, however, that for groups of the same mass, the correlation strength depends on the star formation rate (SFR) of the central galaxy: at fixed mass, the bias of galaxy groups decreases as the SFR of the central galaxy increases. We discuss these findings in light of the recent findings by Gao et al (2005) that halo bias depends on halo formation time, in that halos that assemble earlier are more strongly biased. We also discuss the implication for galaxy formation, and address a possible link to galaxy conformity, the observed correlation between the properties of satellite galaxies and those of their central galaxy.

Using a large galaxy group catalogue constructed from the SDSS, we investigate the correlation between various galaxy properties and halo mass. We split the population of galaxies in early types, late types, and intermediate types, based on their colour and specific star formation rate. At fixed luminosity, the early type fraction increases with increasing halo mass. Most importantly, this mass dependence is smooth and persists over the entire mass range probed, without any break or feature at any mass scale. We argue that the previous claim of a characteristic feature on galaxy group scales is an artefact of the environment estimators used. At fixed halo mass, the luminosity dependence of the type fractions is surprisingly weak: galaxy type depends more strongly on halo mass than on luminosity. We also find that the early type fraction decreases with increasing halo-centric radius. Contrary to previous studies, we find that this radial dependence is also present in low mass haloes. The properties of satellite galaxies are strongly correlated with those of their central galaxy. In particular, the early type fraction of satellites is significantly higher in a halo with an early type central galaxy than in a halo of the same mass but with a late type central galaxy. This phenomenon, which we call `galactic conformity', is present in haloes of all masses and for satellites of all luminosities. Finally, the fraction of intermediate type galaxies is always ~20 percent, independent of luminosity, independent of halo mass, independent of halo-centric radius, and independent of whether the galaxy is a central galaxy or a satellite galaxy. We discuss the implications of all these findings for galaxy formation and evolution.

We measure the projected correlation function w_p(r_p) from the Sloan Digital Sky Survey for a flux-limited sample of 118,000 galaxies and for a volume limited subset of 22,000 galaxies with absolute magnitude M_r<-21. Both correlation functions show subtle but systematic departures from the best-fit power law, in particular a change in slope at r_p~1-2 Mpc/h. These departures are stronger for the volume-limited sample, which is restricted to relatively luminous galaxies. We show that the inflection point in w_p(r_p) can be naturally explained by contemporary models of galaxy clustering, according to which it marks the transition from a large scale regime dominated by galaxy pairs in separate dark matter halos to a small scale regime dominated by galaxy pairs in the same dark matter halo. For example, given the dark halo population predicted by an inflationary cold dark matter scenario, the projected correlation function of the volume-limited sample can be well reproduced by a model in which the mean number of M_r<-21 galaxies in a halo of mass M>M_1=4.74 X 10^{13}\msun/h is <N>_M=(M/M_1)^{0.89}, with 75% of the galaxies residing in less massive, single-galaxy halos. This physically motivated model has the same number of free parameters as a power law, and it fits the w_p(r_p) data better, with a \chi^2/d.o.f.=0.93 compared to 6.12 (for 10 degrees of freedom, incorporating the covariance of the correlation function errors). Departures from a power-law correlation function encode information about the relation between galaxies and dark matter halos. Higher precision measurements of these departures for multiple classes of galaxies will constrain galaxy bias and provide new tests of the theory of galaxy formation.

We use recent observations of the HI-mass function to constrain galaxy formation. The data conflicts with the standard model where most of the gas in a low-mass dark matter halo is assumed to settle into a disk of cold gas that is depleted by star formation and supernova-driven outflows until the disk becomes gravitationally stable. A consistent model can be found if low-mass haloes are embedded in a preheated medium, with a specific gas entropy ~ 10Kev cm^2. Such a model simultaneously matches the faint-end slope of the galaxy luminosity function. We propose a preheating model where the medium around low-mass haloes is preheated by gravitational pancaking. Since gravitational tidal fields suppress the formation of low-mass haloes while promoting that of pancakes, the formation of massive pancakes precedes that of the low-mass haloes within them. We demonstrate that the progenitors of present-day dark matter haloes with M<10^{12}h^{-1}\msun were embedded in pancakes of masses $~5x10^{12}h^{-1}\msun$ at z~2. The formation of such pancakes heats the gas to a temperature of 5x10^5K and compresses it to an overdensity of ~10. Such gas has a cooling time that exceeds the age of the Universe at z~2, and has a specific entropy of ~15Kev cm^2, almost exactly the amount required to explain the stellar and HI mass functions. (Abridged)

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.

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

What galaxy lives in what halo? The answer to this simple question holds important information regarding galaxy formation and evolution. We describe a new statistical technique to link galaxies to their dark matter haloes, or light to mass, using the clustering properties of galaxies as function of their luminosity. The galaxy-dark matter connection thus established, and parameterized through the conditional luminosity function, indicates the presence of two characteristic scales in galaxy formation: one at ~10^{11} Msun, where galaxy formation is most efficient, and another at ~10^{13} Msun, where a transition occurs from systems dominated by one brightest, central galaxy to systems with several dominant galaxies of comparable luminosity. We show that the relation between light and mass established from galaxy clustering alone is in excellent agreement with the dynamical masses inferred from satellite kinematics. We also present a new (halo-based) galaxy-group finder, and show that the multiplicity function of galaxy groups identified in the 2dFGRS suggests a relatively high mass-to-light ratio on the scales of galaxy clusters, or, alternatively, a relatively low value of sigma_8. These findings are also supported by our studies of pairwise peculiar velocities and satellite abundances.

We present results of age and metallicity gradient analysis inferred from both optical and near-infrared surface photometry. The analysis is based on a sample of 36 nearby early-type galaxies, obtained from the Early Data Release of the Sloan Digital Sky Survey and the Two Micron All Sky Survey. Surface brightness profiles were derived in each band, and used to study the color gradients of the galaxies. Using simple stellar population models with both optical and near infrared colors, we may interpret the color gradients in term of age and metallicity gradients of galaxies. Using $g_Z \equiv d \log Z_{\rm met} / d \log R $ and $g_A = d \log {\rm Age} / d \log R $ to represent the metallicity and age gradients, we found a median value of $g_Z=-0.25\pm 0.03$ for the metallicity gradient, with a dispersion $\sigma_{g_Z}=0.19\pm0.02$. The corresponding values for the age gradient were $g_A=0.02\pm 0.04$ and $\sigma_{g_A}=0.25\pm0.03$. These results are in good agreement with recent observational results, as well as with recent simulations that suggest both monolithic collapse and major merger have played important roles in the formation of early-type galaxies. Our results demonstrate the potential of using multi-waveband colors obtained from current and future optical and infrared surveys in constraining the age and metallicity gradients of early-type galaxies.

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.

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.

Recent N-body simulations show that the formation of a present-day, galaxy sized dark matter halo in the CDM cosmogony in general consists of an early fast collapse phase, during which the potential associated with a halo is established, followed by a slow accretion phase, during which mass is added rather gently in the outer region. Here we outline a scenario in which the fast collapse phase is accompanied with rapid formation of cold clouds and with starbursts that can eject a large amount of gas from the halo center.The loss of orbital energy of the cold clouds to the dark matter and the ejection of gas from halo center by starburst can significantly reduce the halo concentration. The outflow from the starburst can also heat the gas in the protogalaxy region. Subsequent formation of galaxies in the slow accretion regime is therefore in halos that have been pre-processed by these processes and may have properties different from that given by N-body simulations. This scenario can help to solve several outstanding problems in the standard Lambda CDM model of galaxy formation without compromising its success in allowing structure formation at high redshift. The predicted rotation curves are quite flat,alleviating the discrepancy of the Tully-Fisher relation predicted in the Lambda CDM model with observations. The flattened galaxy halos allow accreted minihalos to survive in their central regions longer, which may be helpful in producing the flux anomalies observed in some gravitational lensing systems. The preheating by the early starbursts effectively reduces the amount of gas that can be accreted into galaxy halos, which may explain why the baryon fraction in a spiral galaxy is in general much lower than the universal baryon fraction, f_B~0.16, in the standard Lambda CDM model.

Using detailed mock galaxy redshift surveys (MGRSs) we investigate the abundance and radial distribution of satellite galaxies. The mock surveys are constructed using the conditional luminosity function (CLF), and are compared against data from the 2dFGRS. We use Monte Carlo Markov Chains to explore the full posterior distribution of the CLF parameter space, and show that the average relation between light and mass is tightly constrained. The radial number density distribution of satellite galaxies in the 2dFGRS reveals a pronounced absence of satellites at small projected separations from their host galaxies. This is (at least partially) due to the overlap and merging of galaxy images in the APM catalogue. Because of the resulting close-pair incompleteness we are unfortunately unable to put meaningful constraints on the radial distribution of satellite galaxies. Marginalizing over the full CLF parameter space, we show that in a LCDM concordance cosmology the observed abundances of host and satellite galaxies in the 2dFGRS indicate a power spectrum normalization of sigma_8=0.7. The same cosmology but with sigma_8=0.9 is unable to simultanously match the abundances of host and satellite galaxies. This confirms our previous conclusions based on the pairwise peculiar velocity dispersions and the group multipliticy function.

We present models for the three-point correlation function (3PCF) of both dark matter and galaxies. We show that models based on the halo model can reasonably match the dark matter 3PCF obtained from high-resolution $N$-body simulations. On small scales ($r\la 0.5\mpch$) the 3PCF is sensitive to details regarding the density distributions of dark matter haloes. On larger scales ($r \ga 2.0\mpch$) the results are very sensitive to the abundance of the few most prominant haloes. Using the conditional luminosity function, we also construct models for the 3PCF of galaxies, which we test against large mock galaxy samples. The bias of the galaxy distrubution with respect to the dark matter, and the finite number of galaxies that can be hosted by individual haloes, significantly reduce the normalized three-point correlation function with respect to that of dark matter. Contrary to the 3PCF of the dark matter, the galaxy 3PCF is much less sensitive to details regarding the spatial number density distribution of galaxies in individual haloes or to the abundance of the few most massive systems. Finally, we show that our model based on the conditional luminosity function is in good agreement with results obtained from the 2-degree Field Galaxy Redshift Survey. In particular, the model nicely reproduces the observational finding that the 3PCF for early-type galaxies is slightly higher than that of late-type galaxies, and that there is no significant dependence of the 3PCF on galaxy luminosity.

Using detailed mock galaxy redshift surveys we investigate to what extent the kinematics of large samples of satellites galaxies extracted from flux-limited surveys can be used to constrain halo masses. Previous host-satellite selection criteria yield relatively large fractions of interlopers and with a velocity distribution that, contrary to what has typically been assumed, differs strongly from uniform. A new adaptive selection criterion is proposed which yields much larger host-satellite samples and with strongly reduced interloper fractions. We show that satellite weighting, which occurs naturally when stacking many host-satellite pairs, introduces a bias towards higher velocity dispersions compared to the true, host-averaged mean. A further bias, in the same direction, is introduced when using flux-limited, rather than volume-limited surveys. We apply our adaptive selection criterion to the 2dFGRS and obtain a sample of 12613 satellite galaxies and 8004 host galaxies. The satellite kinematics are in excellent agreement with predictions based on the conditional luminosity function (CLF), providing independent, dynamical confirmation of the average mass-to-light ratios predicted by the CLF formalism.

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 use the observed abundance and clustering of galaxies from the 2dFGRS to determine the matter density Omega_m and the linear amplitude of mass fluctuations sigma_8. We use a method based on the conditional luminosity function, which allows straightforward computation of the luminosity dependent bias of galaxies. In addition, it allows the inclusion of constraints on the redshift space distortion parameter beta=Omega_m^{0.6}/b, and yields average mass-to-light ratios as function of halo mass. Using only the luminosity function and the correlation lengths as function of luminosity we obtain constraints on Omega_m and sigma_8 that are in good agreement with COBE. Using additional constraints on the mass-to-light ratios of clusters and on the value of beta as determined from the 2dFGRS yields Omega_m = (0.27 +/- 0.06) and sigma_8 = (0.77 +/- 0.06). Adding further constraints from current CMB data, we obtain Omega_m = (0.25 +/- 0.04) and sigma_8 = (0.78 +/- 0.06). Thus, we find strong indications that both Omega_m and sigma_8 are significantly lower than their ``standard'' concordance values of 0.3 and 0.9, respectively. Finally, we show that the cosmologies favored here predict dark matter haloes that are significantly less centrally concentrated than for the standard LCDM concordance model. This may solve both the problem with the rotation curves of dwarf and LSB galaxies, as well as the problem of simultaneously matching the galaxy luminosity function and the Tully-Fisher zero-point. (Abridged)

(Abridged) Using data from the 2dFGRS we compute the conditional luminosity functions (CLFs) of early- and late-type galaxies. These functions give the average number of galaxies with luminosity L that reside in a halo of mass M, and are a powerful statistical tool to link the distribution of galaxies to that of dark matter haloes. We find that the average mass-to-light ratios of dark matter haloes have a minimum of about 100 Msun/Lsun around a halo mass of about 3x10^{11} Msun. In addition, haloes with M < 10^{10} Msun are virtually devoid of galaxies, while clusters have b_J band mass-to-light ratios in the range 500-1000 Msun/Lsun. Finally, the fact that early-type galaxies are more strongly clustered than late-type galaxies requires that the fraction of late-type galaxies is a strongly declining function of halo mass. We also compute two-point correlation functions as function of both luminosity and galaxy type, and find a weak (strong) luminosity dependence for the late (early) type galaxies. Finally we compare our CLFs with predictions from several semi-analytical models for galaxy formation. As long as these models accurately fit the 2dFGRS luminosity function the agreement with our predictions is remarkably good. This indicates that we have recovered a statistical description of how galaxies populate dark matter haloes which fits nicely within the standard framework for galaxy formation.

abridged: We use a complete sample of about 140,000 galaxies from the Sloan Digital Sky Survey (SDSS) to study the size distribution of galaxies and its dependence on their luminosity, stellar mass, and morphological type. The large SDSS database provides statistics of unprecedented accuracy. For each type of galaxy, the size distribution at given luminosity (or stellar mass) is well described by a log-normal function, characterized by its median $\bar{R}$ and dispersion $\sigma_{\ln R}$. For late-type galaxies, there is a characteristic luminosity at $M_{r,0}\sim -20.5$ (assuming $h=0.7$) corresponding to a stellar mass $M_0\sim 10^{10.6}\Msun$. Galaxies more massive than $M_0$ have $\bar{R}\propto M^{0.4}$ and $\sigma_{\ln R}\sim 0.3$, while less massive galaxies have $\bar{R}\propto M^{0.15}$ and $\sigma_{\ln R}\sim 0.5$. For early-type galaxies, the $\bar{R}$ - $M$ relation is significantly steeper, $\bar{R}\propto M^{0.55}$, but the $\sigma_{\ln R}$ - $M$ relation is similar to that of late-type galaxies. Faint red galaxies have sizes quite independent of their luminosities.

Galactic winds and mass outflows are observed both in nearby starburst galaxies and in high-redshift star-forming galaxies. In this paper we develop a simple analytic model to understand the observed superwind phenomenon. Our model is built upon the model of McKee & Ostriker (1977) for the interstellar medium. It allows one to predict how properties of a superwind, such as wind velocity and mass outflow rate, are related to properties of its star-forming host galaxy, such as size, gas density and star formation rate. The model predicts a threshold of star formation rate density for the generation of observable galactic winds. Galaxies with more concentrated star formation produce superwinds with higher velocities. The predicted mass outflow rates are comparable to (or slightly larger than) the corresponding star formation rates. We apply our model to both local starburst galaxies and high-redshift Lyman break galaxies, and find its predictions to be in good agreement with current observations. Our model is simple, and so can be easily incorporated into numerical simulations and semi-analytical models of galaxy formation.

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

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

We outline a scenario of galaxy formation in which the gas in galaxy-forming regions was preheated to high entropy by vigorous energy feedback associated with the formation of stars in old ellipticals and bulges and with AGN activity. Such preheating likely occurred at redshifts z ~ 2-3, and can produce the entropy excess observed today in low-mass clusters of galaxies without destroying the bulk of the Lyman alpha forest. Subsequent galaxy formation is affected by the preheating, because the gas no longer follows the dark matter on galaxy scales. The hot gas around galaxy haloes has very shallow profiles and emits only weakly in the X-ray. Cooling in a preheated halo is not inside-out, because the cooling efficiency does not change significantly with radius. Only part of the gas in a protogalaxy region can cool and be accreted into the final galaxy halo. The accreted gas is likely in diffuse clouds and so does not lose angular momentum to the dark matter. Cluster ellipticals are produced by mergers of stellar systems formed prior to the preheating, while large galaxy disks form in low-density environments where gas accretion can continue to the present time.

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 the conditional luminosity function (CLF), which gives the number of galaxies with luminosities in the range [L, L+dL] that reside in a halo of mass M, to link the distribution of galaxies to that of dark matter haloes. We seek the CLF that reproduces the galaxy luminosity function and the luminosity dependence of the galaxy clustering strength and test the models by comparing the resulting mass-to-light ratios (M/L) with constraints from the Tully-Fisher (TF) relation. We obtain a number of stringent constraints on both galaxy formation and cosmology. In particular, this method can break the degeneracy between Omega_0 and the power-spectrum normalization sigma_8, inherent in current weak-lensing and cluster-abundance studies. For flat LCDM cosmogonies with sigma_8 normalized by recent weak lensing observations, the best results are obtained for Omega_0~0.3; Omega_0 < 0.2 leads to too large galaxy correlation lengths, while Omega_0 > 0.4 gives too high M/L. The best-fit model for the LCDM concordance cosmology (Omega_0=0.3) predicts M/L that are slightly too high. We discuss a number of possible effects that might remedy this problem, including small modification of cosmological parameters, warm in stead of cold dark matter, systematic errors in current observational data, and the existence of dark galaxies. Finally we use the CLF to predict several statistics about the distribution of galaxy light in the local Universe.

A variety of approximate schemes for modelling the low-density Intergalactic Medium (IGM) in the high-redshift Universe is compared to the results of a large high-resolution hydro-dynamical simulation. These schemes use either an analytical description of the dark matter distribution and the IGM or numerical simulations of the DM distributions combined with different approximate relations between dark matter field and the gas distribution. Schemes based on a filtering of the dark matter distribution with a global Jeans scale result in a rather poor description of the gas distribution. An adaptive filtering which takes into account the density/temperature dependence of the Jeans scale is required. A reasonable description of the gas distribution can be achieved using a fit of the mean relation between the dark matter and gas densities in the hydro-dynamical simulation to relate dark matter and gas distribution. In the hydro-dynamical simulations deviations from this mean relation are correlated with gradients in the dark matter peculiar velocity field indicative of shocks in the gas component. A scheme which takes into account this correlation results in a further improved gas distribution. Such adaptive filtering schemes applied to dark matter simulations will be very well suited for studies of statistical properties of the Lyalpha forest which investigate the IGM and the underlying dark matter distribution and require a large dynamic range and/or an extensive parameter study.

McKay et al. (2002) have recently used measurements of weak galaxy-galaxy lensing in the Sloan Digital Sky Survey to estimate the cross correlation between galaxies and the projected dark matter density field. They derive a relation between aperture mass within a radius of $260\kpch$, $M_{260}$, and lens galaxy luminosity, that does not depend on galaxy luminosity, type or environment. In this paper, we study the cross-correlation between galaxies and dark matter using galaxy catalogs constructed from a high-resolution N-body simulation of a $\Lambda$CDM Universe. We show that our simulations reproduce the McKay et al. results reasonably well. In the simulation, $M_{260}$ is approximately equal to the halo virial mass for $L_*$ galaxies. $M_{260}$ overestimates the virial mass for fainter galaxies and underestimates it for brighter galaxies. We use the simulations to show that under certain circumstances the halo virial mass may be recovered by fitting an NFW profile to the projected galaxy-mass correlation function. If we apply our method to the observations, we find that $L_*$ galaxies typically reside in halos of $\sim 2\times 10^{12} h^{-1} \msun$, consistent with halo masses estimated from the observed Tully-Fisher relation. In the simulations, the halo virial mass scales with galaxy luminosity as $L^{1.5}$ for central galaxies in halos and for galaxies in low-density regions. For satellite galaxies, and for galaxies in high-density regions, there is no simple relation between galaxy luminosity and halo mass and care must be exercised when interpreting the lensing results.

Much evidence suggests that we live in a flat Cold Dark Matter universe with a cosmological constant. Accurate analytic formulae are now available for many properties of the dark halo population in such a Universe. Assuming current ``concordance'' values for the cosmological parameters, we plot halo abundance against redshift as a function of halo mass, of halo temperature, of the fraction of cosmic matter in haloes, of halo clustering strength, and of the clustering strength of the z=0 descendants of high redshift haloes. These plots are useful for understanding how nonlinear structure grows in the model. They demonstrate a number of properties which may seem surprising, for example: 10^9 solar mass haloes are as abundant at z=20 as L_* galaxies are today; 10^6K haloes are equally abundant at z=8 and at z=0; 10% of all matter is currently in haloes hotter than 1 keV, while more than half is in haloes too cool to trap photo-ionized gas; 1% of all matter at z=15 is in haloes hot enough to ionise hydrogen; haloes of given mass or temperature are more clustered at higher redshift; haloes with the abundance of present-day L_* galaxies are equally clustered at all z<20; the metals produced by star-formation at z>10 are more clustered at z=0 than are L_* galaxies.

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

We present the first measurements of clustering in the Sloan Digital Sky Survey (SDSS) galaxy redshift survey. Our sample consists of 29,300 galaxies with redshifts 5,700 km/s < cz < 39,000 km/s, distributed in several long but narrow (2.5-5 degree) segments, covering 690 square degrees. For the full, flux-limited sample, the redshift-space correlation length is approximately 8 Mpc/h. The two-dimensional correlation function \xi(r_p,\pi) shows clear signatures of both the small-scale, ``fingers-of-God'' distortion caused by velocity dispersions in collapsed objects and the large-scale compression caused by coherent flows, though the latter cannot be measured with high precision in the present sample. The inferred real-space correlation function is well described by a power law, \xi(r)=(r/6.1+/-0.2 Mpc/h)^{-1.75+/-0.03}, for 0.1 Mpc/h < r < 16 Mpc/h. The galaxy pairwise velocity dispersion is \sigma_{12} ~ 600+/-100 km/s for projected separations 0.15 Mpc/h < r_p < 5 Mpc/h. When we divide the sample by color, the red galaxies exhibit a stronger and steeper real-space correlation function and a higher pairwise velocity dispersion than do the blue galaxies. The relative behavior of subsamples defined by high/low profile concentration or high/low surface brightness is qualitatively similar to that of the red/blue subsamples. Our most striking result is a clear measurement of scale-independent luminosity bias at r < 10 Mpc/h: subsamples with absolute magnitude ranges centered on M_*-1.5, M_*, and M_*+1.5 have real-space correlation functions that are parallel power laws of slope ~ -1.8 with correlation lengths of approximately 7.4 Mpc/h, 6.3 Mpc/h, and 4.7 Mpc/h, respectively.

(abridged) We present an effective implementation of analytical calculations of the Lyalpha opacity distribution of the Intergalactic Medium (IGM) along multiple lines of sight (LOS) to distant quasars in a cosmological setting. This method fully accounts for the expected correlations between LOS and the cosmic variance in the large-scale modes of the dark matter distribution. Strong correlations extending up to ~300 kpc (proper) and more are found at redshift z~2->3, in agreement with observations. These correlations are investigated using the cross-correlation coefficient and the cross-power spectrum of the flux distribution along different LOS and by identifying coincident absorption features as fitted with a Voigt profile fitting routine. The cross-correlation coefficient between the LOS can be used to constrain the shape-parameter Gamma of the power spectrum if the temperature and the temperature density relation of the IGM can be determined indepedently. We also propose a new technique to recover the 3D linear dark matter power spectrum by integrating over 1D flux cross-spectra which is complementary to the usual `differentiation' of 1D auto-spectra. The cross-power spectrum suffers much less from errors uncorrelated in different LOS, like those introduced by continuum fitting. Investigations of the flux correlations in adjacent LOS should thus allow to extend studies of the dark matter power spectrum with the Lyalpha forest to significantly larger scales than is possible with flux auto-power spectra. 30 pairs with separation of 1-2 arcmin should be sufficient to determine the 1D cross-spectrum at scales of 60 Mpc/h with an accuracy of about 30% if the error is dominated by cosmic variance.

(abridged) We use a theoretical model to predict the clustering properties of galaxy clusters. Our technique accounts for past light-cone effects on the observed clustering and follows the non-linear evolution of the dark matter correlation function and cluster bias factor. A linear treatment of redshift-space distortions is also included. We perform a maximum-likelihood analysis by comparing the theoretical predictions to a set of observational data, both in the optical (APM and EDCC catalogues) and X-ray band (RASS1 Bright Sample, BCS, XBACs, REFLEX). In the framework of CDM models, we compute the constraints on cosmological parameters. Our results show that X-ray data are more powerful than optical ones, allowing smaller regions in the parameter space. We find that the clustering properties of clusters are almost independent of the matter density parameter and of the presence of a cosmological constant, while they appear to be strongly dependent on the shape parameter. Using the X-ray data only, we obtain Gamma~0.1 and 0.4<sigma_8<1.1 for the Einstein-de Sitter model, while 0.14<Gamma<0.22 and 0.6<sigma_8<1.3 for open and flat models with Omega_0m=0.3. Finally, we use our model to make predictions on the correlation length of galaxy clusters expected in future surveys. In particular, we show the results for an optical catalogue with characteristics similar to the EIS project and for a very deep X-ray catalogue with the characteristics of the XMM/LSS survey. We find that clusters at high redshifts are expected to have larger a correlation length than local ones.

Current theory of disk galaxy formation is used to study fundamental-plane (FP) type of relations for disk galaxies. We examine how the changes in model parameters affect these relations and explore the possibility of using such relations to constrain theoretical models. The distribution of galaxy disks in the space of their fundamental properties are predicted to be concentrated in a plane, with the Tully-Fisher (TF) relation (a relation between luminosity $L$ and maximum rotation velocity $V_m$) being an almost edge-on view. Using rotation velocities at larger radii generally leads to larger TF scatter. In searching for a third parameter, we find that both the disk scale-length $R_d$ (or surface brightness) and the rotation-curve shape are correlated with the TF scatter. The FP relation in the $(\Log L, \Log V_m, \Log R_d)$-space obtained from the theory is $L\propto R_d^{\alpha'} V_m^{\beta'}$, with ${\alpha'}\sim 0.50$ and ${\beta'}\sim 2.60$, consistent with the preliminary result we obtain from observational data. Among the model parameters we probe, variation in any of them can generate significant scatter in the TF relation, but the effects of the spin parameter and halo concentration can be reduced significantly by introducing $R_d$ while the scatter caused by varying $m_d$ (the ratio between disk mass and halo mass) is most effectively reduced by introducing the parameters which describes the rotation-curve shape. The TF and FP relations combined should therefore provide useful constraints on models of galaxy formation.

We use the observed sizes and star formation rates of Lyman-break galaxies (LBGs) to understand their properties expected in the hierarchical scenario of galaxy formation. The observed size distribution constrains the masses of the host haloes of LBGs from below, while the observed star-formation distribution constrains the masses from above. Assuming a flat CDM model with a cosmological constant ($\Omega_\Lambda=0.7$), we find that consistency with observations requires the circular velocities of LBG host haloes to be in the range 100--$300 {\rm km s^{-1}}$. The predicted comoving correlation length of these objects is $\sim 3 h^{-1}{\rm Mpc}$, and the predicted velocity dispersion of their stellar contents is typically $70 {\rm km s^{-1}}$. The same prescription applied to larger haloes in the CDM cosmogony predicts the existence of galaxies with star formation rates $\sim 1000 {\rm M_\odot yr^{-1}}$ at redshift $z\sim 3$. We explore the possibility of identifying these galaxies to be the bright sub-millimeter (sub-mm) sources detected by SCUBA. The model predicts that the host haloes of these sub-mm sources are massive. These galaxies are predicted to be strongly correlated, with a comoving correlation length of $\sim 7 h^{-1}{\rm Mpc}$. The descendants of the bright sub-mm sources should reside in clusters of galaxies at the present time, and it is likely that these objects are the progenitors of giant ellipticals. The model predicts the existence of a relatively bright and red population of galaxies at $z\sim 3$, which may be observed in the K-band.

Models for the formation of galaxies and clusters of galaxies require strong feedback in order to explain the observed properties of these systems. We investigate whether such feedback has observational consequences for the intergalactic medium, as probed in absorption towards background quasars. A typical quasar sight-line intersects one proto-cluster per unit redshift, and significant feedback from forming galaxies or AGN, heating the proto-cluster gas, will result in a large clearing of reduced absorption in the Lyman alpha-forest. Such a gap could be detected at redshift > 3 when the mean opacity is high. Feedback from Lyman-break galaxies in proto-clusters can be probed by the absorption lines produced in their winds. Strong feedback from galaxies has a major impact on the number and properties of absorption lines with neutral hydrogen column densities of 10^16/ cm^2. This feedback can be probed with multiple sight-lines and by studying the unsaturated higher-order lines of the Lyman series. Galactic winds from dwarf galaxies should break-up into clouds, in order not to over produce the number of absorption lines. These clouds can then coast to large distances.

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

We show that the Tully-Fisher relation observed for spiral galaxies can be explained in the current scenario of galaxy formation without invoking subtle assumptions, provided that galactic-sized dark haloes have shallow, core-like central profiles, with a core radius proportional to halo circular velocity. In such a system, both disk and halo contribute significantly to the maximum rotation of the disk, and the interaction between the disk and halo components acts to reduce the scatter in the Tully-Fisher relation. With model parameters chosen in plausible ranges, the model can well accommodate the zero-point, slope, scatter of the observed Tully-Fisher relation, as well as the large range of disk sizes. The model predicts that LSB disks obey a Tully-Fisher relation similar to that of normal disks, if disk mass-to-light ratio is properly taken into account. The halo profile required by the Tully-Fisher relation is as shallow as that required by the rotation curves of faint disks, but much shallower than that predicted by conventional CDM models. Our results cannot be explained by some of the recent proposals for resolving the conflict between conventional CDM models and the rotation-curve shapes of faint galaxies. If dark matter self-interaction (either scattering or annihilation) is responsible for the shallow profile, the observed Tully-Fisher relation requires the interaction cross section \sigma_X to satisfy <\sigma_{X}|v|>/m_{X}~10^{-16} cm^3/s/GeV.

The Press--Schechter, excursion set approach allows one to make predictions about the shape and evolution of the mass function of bound objects. It combines the assumption that objects collapse spherically with the assumption that the initial density fluctuations were Gaussian and small. While the predicted mass function is reasonably accurate at the high mass end, it has more low mass objects than are seen in simulations of hierarchical clustering. This discrepancy can be reduced substantially if bound structures are assumed to form from an ellipsoidal, rather than a spherical collapse. In the spherical model, a region collapses if the initial density within it exceeds a threshold value, delta_sc. This value is independent of the initial size of the region, and since the mass of the collapsed object is related to its initial size, delta_sc is independent of final mass. In the ellipsoidal model, the collapse of a region depends on the surrounding shear field, as well as on its initial overdensity. Therefore, there is a relation between the density threshold value required for collapse, and the mass of the final object. We provide a fitting function to this delta_ec(m) relation for initially Gaussian fields which simplifies the inclusion of ellipsoidal dynamics in the excursion set approach. We discuss the relation between the excursion set predictions and the halo distribution in high resolution N-body simulations, and show that our simple parametrization of the ellipsoidal collapse model represents a significant improvement on the spherical model on an object-by-object basis. Finally, we show that the associated statistical predictions, the mass function and the large scale halo-to-mass bias relation, are also more accurate than the standard predictions.

We review some of our recent progress in modelling the formation of disk galaxies in the framework of hierarchical structure formation. Our model is not only consistent with the local observations, but also provides a good description of disk evolution out to redshift z~1. We use this model to interpret recent observational results on Lyman-break galaxies (LBGs) at z~3. Assuming that these galaxies are associated with the most massive halos and adopting an empirical law for their star formation rates, we find that many properties of the LBG population, including their correlations, sizes and kinematics can be accommodated in the model.

It has been suggested recently that the rate of gamma-ray bursts (GRBs) is proportional to the star formation rate in the universe. In this paper, we study the nature of GRB hosts expected in this scenario. We improve upon previous studies by incorporating a luminosity function for the GRBs, as required by observations. This model provides a good match to the observed number counts of GRBs as a function of peak-count rate. The model predicts that the host galaxies have their redshift distribution peaked around z~1, and about 15 percent have z>2.5. This high-redshift fraction have the same properties as the star-forming galaxies recently discovered by the Lyman-break technique. At z<~1, many of the GRBs may be hosted by faint blue galaxies. Using a photometric redshift sample of galaxies from the Hubble Deep Field, we find that the host galaxies have magnitudes in the range from 21.5 to 28 in the I-band, and about 90 percent of them have semi-major axis smaller than 1.3 arcsec. Assuming isotropic emission, the typical peak-luminosity and total energy of GRBs are in an Einstein-de Sitter universe with H_0=100km/s/Mpc. We also discuss further observational tests of this scenario.

The number density and clustering properties of Lyman-break galaxies (LBGs) are consistent with them being the central galaxies of the most massive dark halos present at z~3. This conclusion holds in all currently popular hierarchical models for structure formation, and is almost independent of the global cosmological parameters. We examine whether the sizes, luminosities, kinematics and star-formation rates of LBGs are also consistent with this identification. Simple formation models tuned to give good fits to low redshift galaxies can predict the distribution of these quantities in the LBG population. The LBGs should be small (with typical half-light radii of 0.6-2 kpc/h), should inhabit haloes of moderately high circular velocity (180-290 km/s) but have low stellar velocity dispersions (70-120 km/s) and should have substantial star formation rates (15-100 Msun/yr). The numbers here refer to the predicted median values in the LBG sample of Adelberger et al. (1998); the first assumes an Omega=1 universe and the second a flat universe with Omega=0.3. For either cosmology these predictions are consistent with the current (rather limited) observational data. Following the work of Kennicutt (1998) we assume stars to form more rapidly in gas of higher surface density. This predicts that LBG samples should preferentially contain objects with low angular momentum, and so small size, for their mass. In contrast, samples of damped Lyman alpha systems (DLSs), should be biased towards objects with large angular momentum. Bright LBGs and DLSs may therefore form distinct populations, with very different sizes and star formation rates, LBGs being smaller and more metal-rich than DLSs of similar mass and redshift.

We use data from the literature to constrain theoretical models of galaxy formation. We show how to calculate the dimensionless spin parameter $\lambda$ of the halos of disk galaxies and we compare the distribution of $\lambda$ with that observed in cosmological $N$-body simulations. The agreement is excellent, which provides strong support for the hierarchical picture of galaxy formation. Assuming only that the radial surface density distribution of disks is exponential, we estimate crudely the maximum-disk mass-to-light ratio in the $I$-band and obtain $<M/L> \la 3.56h$, for a Hubble constant of $100h~\kms \mpc^{-1}$. We discuss this result and its limitations in relation to other independent determinations of $M/L$. We also define a dimensionless form of the Tully-Fisher relation, and use it to derive a value of the baryon fraction in disk galaxies. For galaxies with circular velocity $\vmax>100\kms$, the median value is $md = 0.086 (M/L/3.56h)$. Assuming that the gas fraction in galactic halos is at most as large as that in clusters, we also conclude that $<M/L> \la 2.48h^{-1/2}$.

We suggest that the bulge-to-disc ratios of spiral galaxies are primarily determined by the angular momenta of their host haloes predicted in current hierarchical clustering models for structure formation. Gas with low specific angular momentum becomes self-gravitating and presumably forms stars before it can settle into a rotationally supported disc. We assume this part of the gas in a dark halo to form a bulge, while the rest is assumed to settle into a rotationally supported disc. With these assumptions the predicted bulge-to-disc ratios in mass and in size, and other correlations between the bulge and disc components can match current observational results. This model predicts the existence of a population of low-surface-brightness galaxies which are bulgeless. The model also predicts that the bulge component has many properties in common with the disc component, because both form through similar processes. In particular, many bulges should be supported (at least partially) by rotation.

We use recent observations of high-redshift galaxies to study the evolution of galactic disks over the redshift range 0<z<1. The data are inconsistent with models in which disks were already assembled at z=1 and have evolved only in luminosity since that time. Assuming that disk properties change with redshift as powers of 1+z and analysing the observations assuming an Einstein-de Sitter universe, we find that for given rotation speed, disk scalelength decreases with z as ~(1+z)^{-1}, total B-band mass-to-light ratio decreases with z as ~(1+z)^{-1}, and disk luminosity (again in B) depends only weakly on z. These scalings are consistent with current data on the evolution of disk galaxy abundance as a function of size and luminosity. Both the scalings and the abundance evolution are close to the predictions of hierarchical models for galaxy formation. If different cosmogonies are compared, the observed evolution in disk-size and disk abundance favours a flat low-density universe over an Einstein-de Sitter 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.

We propose that a large fraction of the QSO Lyman limit absorption systems (LLS) observed at high redshift (z > 3) originate from gas trapped in small objects, such as minihalos, that form prior to reionization. In the absence of a strong UV flux, the gas is predominantly neutral and may form clouds with HI column density NHI > 10^18 cm^-2. Due to their high densities and high HI column densities, these clouds are not destroyed by the onset of the UV background at a later time. Thus, if not disrupted by other processes, such as mergers into larger systems or `blow away' by supernovae, they will produce LLS. We show that the observed number density of LLS at high redshifts can well be reproduced by the survived `minihalos' in hierarchical clustering models such as the standard cold dark matter model. The number density of LLS in such a population increases with z even beyond the redshifts accessible to current observations and dies off quickly at z < 2. This population is distinct from other populations because the absorbing systems have small velocity widths and a close to primordial chemical composition. The existence of such a population requires that the reionization of the universe occurs late, at z < 20.

We examine the global stability of an exponential stellar disk embedded in a dark matter halo and constrain the mass-to-light ratio of observed disks in the I-band, Y_I. Assuming only that the radial surface density distribution of disks is exponential, we derive an analytic upper limit: Y_I <~ 2.4h, for a Hubble constant of 100h~\kms\mpc^{-1}. Using N-body simulations we derive a stability criterion significantly different from that of previous authors. Using this criterion we argue that, almost independent of the concentration of the halo, Y_I <~ 1.9h. We discuss this result in relation to other independent determinations of Y_I, its limitations and its implications for theories of disk formation and barred galaxies.

Recent observations have shown that (some) metal line absorption systems in QSO spectra arise from gas in galactic halos with radii much larger than the optical radii of galaxies. We show that the observed galaxy-absorber connections are natural results of galaxy formation in hierarchical cosmogonies, where galaxies form from gas cooling and condensation in dark matter halos.

Galaxies can be classified in two broad sequences which are likely to reflect their formation mechanism. The `main sequence', consisting of spirals, irregulars and all dwarf galaxies, is probably produced by gas settling within dark matter haloes. We show that the sizes and surface densities along this sequence are primarily determined by the distributions of the angular momentum and formation time of dark haloes. They are well reproduced by current cosmogonies provided that galaxies form late, at $z \la 2$. In this scenario, dwarf ellipticals were small `disks' at $z\sim 1$ and become `ellipticals' after they fall into cluster environments. The strong clustering of dwarf ellipticals is then a natural by-product of the merging and transformation process. The number of dwarf galaxies predicted in a cluster such as Virgo is in good agreement with the observed number. On the other hand, the `giant branch', consisting of giant ellipticals and bulges, is probably produced by the merging of disk galaxies. Based on the observed phase-space densities of galaxies, we show that the main bodies of all giant ellipticals can be produced by dissipationless mergers of high-redshift disks. However, high-redshift disks, although denser than present-day ones, are still not compact enough to produce the high central phase space density of some low-luminosity ellipticals. Dissipation must have occurred in the central parts of these galaxies during the merger which formed them.

We study the population of galactic disks expected in current hierarchical clustering models for structure formation. A rotationally supported disk with exponential surface density profile is assumed to form with a mass and angular momentum which are fixed fractions of those of its surrounding dark halo. We assume that haloes respond adiabatically to disk formation, and that only stable disks can correspond to real systems. With these assumptions the predicted population can match both present-day disks and the damped Lyman alpha absorbers in QSO spectra. Good agreement is found provided: (i) the masses of disks are a few percent of those of their haloes; (ii) the specific angular momenta of disks are similar to those of their haloes; (iii) present-day disks were assembled recently (at z<1). In particular, the observed scatter in the size-rotation velocity plane is reproduced, as is the slope and scatter of the Tully-Fisher relation. The zero-point of the TF relation is matched for a stellar mass-to-light ratio of 1 to 2 h in the I-band, consistent with observational values derived from disk dynamics. High redshift disks are predicted to be small and dense, and could plausibly merge together to form the observed population of elliptical galaxies. In many (but not all) currently popular cosmogonies, disks with rotation velocities exceeding 200 km/s can account for a third or more of the observed damped Lyman alpha systems at z=2.5. Half of the lines-of-sight to such systems are predicted to intersect the absorber at r>3kpc/h and about 10% at r>10kpc/h. The cross-section for absorption is strongly weighted towards disks with large angular momentum and so large size for their mass. The galaxy population associated with damped absorbers should thus be biased towards low surface brightness systems.

We discuss some recent progress in constructing analytic approximations to the galaxy clustering. We show that successful models can be constructed for the clustering of both dark matter and dark matter haloes. Our understanding of galaxy clustering and galaxy biasing can be greatly enhanced by these models.

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.

Hamilton et al. (1991) proposed a simple formula relating the nonlinear autocorrelation function of the mass distribution to the primordial spectrum of density fluctuations for gravitational clustering in an $\Omega=1$ universe. High resolution N-body simulations show this formula to work well for scale-free spectra $P(k) \propto k^n$ when the spectral index $n\sim 0$, but not when $n \lsim -1$. We show that a modified version of the formula can work well provided its form depends on $n$. This dependence can be derived from a simple physical model for collapse from Gaussian initial conditions. Our modified formula is easy to apply and is an excellent fit to N-body simulations with $0 \leq n \leq -2$. It can also be applied to non-power law initial spectra such as that of the standard Cold Dark Matter model by using the local spectral index at the current nonlinear scale as the effective value of $n$ at any given redshift. We give analytic expressions both for the nonlinear correlation function and for the nonlinear power spectrum.

A recent observation of Steidel et al. indicates that a substantial fraction of giant galaxies were formed at an epoch as early as redshift $z>3-3.5$. We show that this early formation of giant galaxies gives strong constraints on models of cosmic structure formation. Adopting the COBE normalization for the density perturbation spectrum, we argue that the following models do not have large enough power on galactic scales to yield the observed abundance: (i) standard cold dark matter (CDM) models (where mass density $\Omega_0=1$ and power index $n=1$) with the Hubble constant $h\lsim 0.35$; (ii) tilted CDM models with $h=0.5$ and $n\lsim 0.75$; (iii) open CDM models with $h\lsim 0.8$ and $\Omega_0 \lsim 0.3$, and (iv) mixed dark matter models with $h=0.5$ and $\Omega _\nu \gsim 0.2$. Flat CDM models with a cosmological constant $\lambda_0 \sim 0.7$ are consistent with the observation, provided that $h\gsim 0.6$. Combined with constraints from large-scale structure formation, these results imply that the flat CDM model with a low $\Omega_0$ is the only one that is fully consistent with observations. We predict that these high-redshift galaxies are more strongly clustered than normal galaxies observed today.

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.

A model of Lyman limit QSO absorption systems is investigated where they are produced in gaseous galactic halos with a two-phase structure: a hot phase at the halo virial temperature in approximate pressure equilibrium, and a cold, photoionized phase in the form of clouds confined by the pressure of the hot medium, and falling through it to the halo center. We show that the masses of these clouds must be over a relatively narrow range so that they are stable against gravitational collapse and evaporation. These masses, near $10^6 \msun$, also lead to a covering factor near unity. The hot phase is required to have a core radius such that its cooling time in the core is equal to the age of the halo, and the mass in the cold phase is determined by the rate at which the hot gas cools. We calculate the number of Lyman limit systems arising in halos of different masses in a CDM model and their impact parameters implied by these assumptions, and find them to be in reasonable agreement with observations. The observed low ionization systems such as MgII are well reproduced as arising from the photoionized phase, and are mostly in halos around massive galaxies, while CIV selected systems are predicted to be more commonly associated with lower mass galaxies and larger impact parameters. The hot phase may give rise to detectable absorption lines in OVI, while the column densities predicted for other highly ionized species are low and difficult to observe.

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.

We develop a simple analytic model for the gravitational clustering of dark matter haloes to understand how their spatial distribution is biased relative to that of the mass. The statistical distribution of dark haloes within the initial density field (assumed Gaussian) is determined by an extension of the Press-Schechter formalism. Modifications of this distribution caused by gravitationally induced motions are treated using a spherical collapse approximation. We test this model against results from a variety of N-body simulations, and find that it gives an accurate description of a bias function. This bias function is sufficient to calculate the cross-correlation between dark haloes and mass, and again we find excellent agreement between simulation results and analytic predictions. Because haloes are spatially exclusive, the variance in the count of objects within spheres of fixed radius and overdensity is significantly smaller than the Poisson value. This seriously complicates any analytic calculation of the autocorrelation function of dark halos. Our simulation results show that this autocorrelation function is proportional to that of the mass over a wide range in $R$, even including scales where both functions are significantly greater than unity. The constant of proportionality is very close to that predicted on large scales by the analytic model. This result permits an entirely analytic estimate of the autocorrelation function of dark haloes. We use our model to study how the distribution of galaxies may be biased with respect to that of the mass. In conjunction with other data these techniques should make it possible to measure the amplitude of cosmic mass fluctuations and the density of the Universe.

We develop a simple analytic model for the gravitational clustering of dark haloes. The statistical properties of dark haloes are determined from the initial density field (assumed to be Gaussian) through an extension of the Press-Schechter formalism. Gravitational clustering is treated by a spherical model which describes the concentration of dark haloes in collapsing regions. We test this model against results from a variety of N-body simulations. The autocorrelation function of dark haloes in such simulations depends significantly on how haloes are identified. Our predictions agree well with results based on algorithms which break clusters into subgroups more efficiently than the standard friends-of-friends algorithm. The agreement is better than that found by assuming haloes to lie at the present positions of peaks of the linear density field. We use these techniques to study how the distribution of haloes is biased with respect to that of the mass. The initial (Lagrangian) positions of haloes identified at a given redshift and having circular velocities $v_c=v_c^*(z)$ (i.e. mass equal to the characteristic nonlinear mass $M^*$ at that redshift) are very weakly correlated with the linear density field or among themselves. As a result of dynamical evolution, however, the present-day correlations of these haloes are similar to those of the mass. Haloes with lower $v_c$ are biased toward regions with negative overdensity, while those with higher $v_c$ are biased toward regions with positive overdensity. Among the haloes identified at any given epoch, those with higher circular velocities are more strongly correlated today. Among the haloes of given circular velocity, those at higher redshifts are also more strongly clustered today. In the ``standard CDM'' model, haloes with

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.

Motivated by recent observation of Lanzetta et al. that most luminous galaxies at low redshifts produce $\lya$ absorption at impact parameter $l\lsim 160 \kpch$, we propose that these absorbers are clouds confined by the pressure of ambient hot gas in galactic haloes. We determine the properties of this hot gas and of the absorption systems on the basis of observational and theoretical constraints. The absorbing clouds need to be replenished on about one orbital time ($\sim 10^9$ yrs) in the galactic halo. The pressure and temperature of the gas at radius $r\sim 100\kpc$ are $P=(10-100){\rm cm^{-3} K}$, $T=10^{(5.5-6.5)}{\rm K}$. The model requires that about 10 per cent of the gas in low-redshift galactic haloes is in the hot phase. Such gas in galactic haloes emits x-ray with bolometric luminosity of the order $10^{37-40}{\rm erg\,s^{-1}}$. The plausibility for such gas to exist in current models of galaxy formation is discussed.

We examine the constraints on theories of galaxy formation that are obtained from observations of damped $\lya$ (DL) systems, assuming they are gaseous protodisks in dark matter halos. Using the Press-Schechter formalism, we find that the mixed dark matter model, with $\ohdm = 0.3$, $\ocdm = 0.65$, $\obaryon = 0.05$, and $h=0.5$, is ruled out because the number of galactic halos at $z\simeq 3$ is too small to account for the total gaseous mass in DL systems, even under the assumption that all the gas in collapsed halos has settled into disks of neutral gas. The standard CDM model can account for the gas in DL systems if the bias is $b\lsim 2$; the same is true for the CDM model with a cosmological constant, if $b\lsim 1.5$ for $\Lambda = 0.8$. However, one still needs to assume that a fraction $\gsim 0.4$ of the baryons in collapsed halos at $z\simeq 3$ is in the form of neutral gas in disks. We also calculate the column density distribution $f(\nhi)$ of the DL systems, in terms of the surface density profiles of disks and the distribution of their central column densities. It is shown that the form of $f(\nhi)$ at the high end of column density is a diagnostic for the nature of DL systems.

Using redshift samples, we calculate the cross-correlation functions of LSB galaxies with normal galaxies in complete samples (i.e. CfA and IRAS), which enables us to compare directly the amplitudes and shapes of the correlation functions. For pair separations $r\gs 2\mpc$, we find that the shape is in agreement with that of the correlation functions for other galaxies. The amplitudes ($A$) of $\xiab (r)$ are lower than those of the autocorrelation functions for the CfA and IRAS samples, with $A_{\rm LSB-CfA}:A_{\rm CfA-CfA}\approx 0.4$ and $A_{\rm LSB-IRAS}:A_{\rm IRAS-IRAS}\approx 0.6$. These results suggest that LSB galaxies are imbedded in the same large scale structure as other galaxies, but are less strongly clustered. This offers the hope that LSB galaxies may be unbiased tracers of the mass density on large scales. For $r\ls 2\mpc$, the cross-correlation functions are significantly lower than that expected from the extrapolation of $\xiab$ on larger scales, showing that the formation and survival of LSB galaxies may be inhibited by interaction with neighboring galaxies. The model which matches these observations suggests that strong luminosity segregation in galaxy clustering is not a necessary consequence of biased galaxy formation, unless the effect of surface brightness is taken in to account. It is also implies a significant mass density in LSB galaxies.

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.