X-ray and microwave cluster scaling relations are immensely valuable for cosmological analysis. However, their power is limited by astrophysical systematics that bias mass estimates and introduce additional scatter. Turbulence injected into the intracluster medium via mass assembly contributes substantially to cluster non-thermal pressure support, a significant source of such uncertainties. We use an analytical model to compute the assembly-driven non-thermal pressure profiles of haloes based on Monte Carlo-generated accretion histories. We introduce a fitting function for the average non-thermal pressure fraction profile, which exhibits minimal dependence on redshift at fixed peak height. Using the model, we predict deviations from self-similarity and the intrinsic scatter in the Sunyaev--Zel'dovich effect observable-mass scaling relation ($Y_\mathrm{SZ}-M$) due solely to inter-cluster variation in mass accretion histories. We study the dependence of $Y_\mathrm{SZ}-M$ on aperture radius, cosmology, redshift, and mass limit. The model predicts $5-9\%$ scatter in $Y_\mathrm{SZ}-M$ at $z=0$, increasing as the aperture used to compute $Y_\mathrm{SZ}$ increases from $R_\mathrm{500c}$ to $5R_\mathrm{500c}$. The predicted scatter lies slightly below that of studies based on non-radiative hydro-simulations, illustrating that assembly history variance is likely responsible for a substantial fraction of scatter in $Y_\mathrm{SZ}-M$. This should be regarded as a lower bound, which will likely increase with the use of an updated gas density model that incorporates a more realistic response to halo assembly. As redshift increases, $Y_\mathrm{SZ}-M$ deviates more from self-similarity and scatter increases. We show that the $Y_\mathrm{SZ}-M$ residuals correlate strongly with the recent halo mass accretion rate, potentially providing an opportunity to infer the latter.

We present a semi-analytic model of satellite galaxies, SatGen, which can generate large samples of satellite populations for a host halo of desired mass, redshift, and assembly history. The model combines dark-matter halo merger trees, empirical relations for the galaxy-halo connection, and analytic prescriptions for tidal effects, dynamical friction, and ram pressure stripping. SatGen emulates cosmological zoom-in hydro-simulations in certain aspects. Satellites can reside in cored or cuspy DM subhaloes, depending on the halo response to baryonic physics that can be formulated from hydro-simulations and physical modeling. The subhalo profile and the stellar mass and size of a satellite evolves depending on its tidal mass loss and initial structure. The host galaxy can include a baryonic disc and a stellar bulge, each described by a density profile that allows analytic orbit integration. SatGen complements simulations by propagating the effect of halo response found in simulated field galaxies to satellites (not properly resolved in simulations) and outperforms simulations by sampling the halo-to-halo variance of satellite statistics and overcoming artificial disruption due to insufficient resolution. As a first application, we use the model to study satellites of Milky Way sized hosts, making it emulate simulations of bursty star formation and of smooth star formation, respectively, and to experiment with a disc potential in the host halo. Our model reproduces the observed satellite statistics reasonably well. Different physical recipes make a difference in satellite abundance and spatial distribution at the 25% level, not large enough to be distinguished by current observations given the halo-to-halo variance. The MW disc depletes satellites by 20% and has a subtle effect of diversifying the internal structure of satellites, important for alleviating certain small-scale problems.

The concentration parameter is a key characteristic of a dark matter halo that conveniently connects the halo's present-day structure with its assembly history. Using 'Dark Sky', a suite of cosmological $N$-body simulations, we investigate how halo concentration evolves with time and emerges from the mass assembly history. We also explore the origin of the scatter in the relation between concentration and assembly history. We show that the evolution of halo concentration has two primary modes: (1) smooth increase due to pseudo-evolution; and (2) intense responses to physical merger events. Merger events induce lasting and substantial changes in halo structures, and we observe a universal response in the concentration parameter. We argue that merger events are a major contributor to the uncertainty in halo concentration at fixed halo mass and formation time. In fact, even haloes that are typically classified as having quiescent formation histories experience multiple minor mergers. These minor mergers drive small deviations from pseudo-evolution, which cause fluctuations in the concentration parameters and result in effectively irreducible scatter in the relation between concentration and assembly history. Hence, caution should be taken when using present-day halo concentration parameter as a proxy for the halo assembly history, especially if the recent merger history is unknown.

We study the effects of Kelvin Helmholtz Instability (KHI) on the cold streams that feed massive halos at high redshift, generalizing our earlier results to include the effects of radiative cooling and heating from a UV background, using analytic models and high resolution idealized simulations. We currently do not consider self-shielding, thermal conduction or gravity. A key parameter in determining the fate of the streams is the ratio of the cooling time in the turbulent mixing layer which forms between the stream and the background following the onset of the instability, $t_{\rm cool,mix}$, to the time in which the mixing layer expands to the width of the stream in the non-radiative case, $t_{\rm shear}$. This can be converted into a critical stream radius, $R_{\rm s,crit}$, such that $R_{\rm s}/R_{\rm s,crit}=t_{\rm shear}/t_{\rm cool,mix}$. If $R_{\rm s}<R_{\rm s,crit}$, the non-linear evolution proceeds similarly to the non-radiative case studied by Mandelker et al. 2019a. If $R_{\rm s}>R_{\rm s,crit}$, which we find to almost always be the case for astrophysical cold streams, the stream is not disrupted by KHI. Rather, background mass cools and condenses onto the stream, and can increase the mass of cold gas by a factor of $\sim 3$ within 10 stream sound crossing times. The mass entrainment induces thermal energy losses from the background and kinetic energy losses from the stream, which we model analytically. Roughly half of the dissipated energy is radiated away from gas with $T<5\times 10^4{\rm K}$, suggesting much of it will be emitted in Ly$\alpha$.

We present an analytic toy model for the radiation produced by the interaction between the cold streams thought to feed massive halos at high redshift and their hot CGM. We begin by deriving cosmologically motivated parameters for the streams as they enter the halo virial radius, $R_{\rm v}$, as a function of halo mass and redshift. For $10^{12}M_{\odot}$ halos at $z=2$, we find the Hydrogen number density in streams to be $n_{\rm H,s}\sim (0.1-5)\times 10^{-2}{\rm cm}^{-3}$, a factor of $\delta \sim (30-300)$ times denser than the hot CGM density, while the stream radii are in the range $R_{\rm s}\sim (0.03-0.50)R_{\rm v}$. As the streams accelerate towards the halo centre, they become denser and narrower. The stream-CGM interaction induces Kelvin-Helmholtz Instability (KHI), which leads to entrainment of CGM mass by the stream and therefore to stream deceleration by momentum conservation. Assuming that the entrainment rates derived by Mandelker et al. 2019 in the absence of gravity can be applied locally at each halocentric radius, we derive equations of motion for the stream in the halo. Using these, we derive the net acceleration, mass growth, and energy dissipation induced by the stream-CGM interaction, as a function of halo mass and redshift, for different CGM density profiles. For the range of model parameters considered, we find that the interaction can induce dissipation luminosities $L_{\rm diss}>10^{42}~{\rm erg~s^{-1}}$ within $\le 0.6 R_{\rm v}$ of halos with $M_{\rm v}>10^{12}M_{\odot}$ at $z=2$, with the emission scaling with halo mass and redshift approximately as $\propto M_{\rm v}\,(1+z)^2$. The magnitude and spatial extent of the emission produced in massive halos at high redshift is consistent with observed Ly$\alpha$ blobs, though better treatment of the UV background and self-shielding is needed to solidify this conclusion.

We investigate dynamical self-friction, the process by which material that is stripped from a subhalo torques its remaining bound remnant, which causes it to lose orbital angular momentum. By running idealized simulations of a subhalo orbiting within an analytical host halo potential, we isolate the effect of self-friction from traditional dynamical friction due to the host halo. While at some points in a subhalo's orbit the torque of the stripped material can boost the orbital angular momentum of the remnant, the net effect over the long term is orbital decay regardless of the initial orbital parameters or subhalo mass. In order to quantify the strength of self-friction, we run a suite of simulations spanning typical host-to-subhalo mass ratios and orbital parameters. We find that the time-scale for self-friction, defined as the exponential decay time of the subhalo's orbital angular momentum, scales with mass ratio and orbital circularity similar to standard dynamical friction. The decay time due to self-friction is roughly an order of magnitude longer, suggesting that self-friction only contributes at the 10 percent level. However, along more radial orbits, self-friction can occasionally dominate over dynamical friction close to pericentric passage, where mass stripping is intense. This is also the epoch at which the self-friction torque undergoes large and rapid changes in both magnitude and direction, indicating that self-friction is an important process to consider when modeling pericentric passages of subhaloes and their associated satellite galaxies.

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.

Accurately predicting the abundance and structural evolution of dark matter subhaloes is crucial for understanding galaxy formation, modeling galaxy clustering, and constraining the nature of dark matter. Due to the nonlinear nature of subhalo evolution, cosmological $N$-body simulations remain its primary method of investigation. However, it has recently been demonstrated that such simulations are still heavily impacted by artificial disruption, diminishing the information content on small scales and reducing the reliability of all simulation-calibrated semi-analytical models. In this paper, we utilize the recently released DASH library of high-resolution, idealized simulations of the tidal evolution of subhaloes, which are unhindered by numerical overmerging due to discreteness noise or force softening, to calibrate an improved, more-accurate model of the evolution of the density profiles of subhaloes that undergo tidal heating and stripping within their host halo. By testing previous findings that the structural evolution of a tidally truncated subhalo depends solely on the fraction of mass stripped, independent of the details of the stripping, we identify an additional dependence on the initial subhalo concentration. We provide significantly improved fitting functions for the subhalo density profiles and structural parameters ($V_\mathrm{max}$ and $r_\mathrm{max}$) that are unimpeded by numerical systematics and applicable to a wide range of parameter space. This model will be an integral component of a future semi-analytical treatment of substructure evolution, which can be used to predict key quantities, such as the evolved subhalo mass function and annihilation boost factors, and validate such calculations performed with cosmological simulations.

Extracting accurate cosmological information from galaxy-galaxy and galaxy-matter correlation functions on non-linear scales ($\lesssim 10 h^{-1} \mathrm{Mpc}$) requires cosmological simulations. Additionally, one has to marginalise over several nuisance parameters of the galaxy-halo connection. However, the computational cost of such simulations prohibits naive implementations of stochastic posterior sampling methods like Markov chain Monte Carlo (MCMC) that would require of order $\mathcal{O}(10^6)$ samples in cosmological parameter space. Several groups have proposed surrogate models as a solution: a so-called emulator is trained to reproduce observables for a limited number of realisations in parameter space. Afterwards, this emulator is used as a surrogate model in an MCMC analysis. Here, we demonstrate a different method called Cosmological Evidence Modelling (CEM). First, for each simulation, we calculate the Bayesian evidence marginalised over the galaxy-halo connection by repeatedly populating the simulation with galaxies. We show that this Bayesian evidence is directly related to the posterior probability of cosmological parameters. Finally, we build a physically motivated model for how the evidence depends on cosmological parameters as sampled by the simulations. We demonstrate the feasibility of CEM by using simulations from the Aemulus simulation suite and forecasting cosmological constraints from BOSS CMASS measurements of redshift-space distortions. Our analysis includes an exploration of how galaxy assembly bias affects cosmological inference. Overall, CEM has several potential advantages over the more common approach of emulating summary statistics, including the ability to easily marginalise over highly complex models of the galaxy-halo connection and greater accuracy, thereby reducing the number of simulations required.

Most models for the connection between galaxies and their haloes ignore the possibility that galaxy properties may be correlated with halo properties other than mass, a phenomenon known as galaxy assembly bias. Yet, it is known that such correlations can lead to systematic errors in the interpretation of survey data. At present, the degree to which galaxy assembly bias may be present in the real Universe, and the best strategies for constraining it remain uncertain. We study the ability of several observables to constrain galaxy assembly bias from redshift survey data using the decorated halo occupation distribution (dHOD), an empirical model of the galaxy--halo connection that incorporates assembly bias. We cover an expansive set of observables, including the projected two-point correlation function $w_{\mathrm{p}}(r_{\mathrm{p}})$, the galaxy--galaxy lensing signal $\Delta \Sigma(r_{\mathrm{p}})$, the void probability function $\mathrm{VPF}(r)$, the distributions of counts-in-cylinders $P(N_{\mathrm{CIC}})$, and counts-in-annuli $P(N_{\mathrm{CIA}})$, and the distribution of the ratio of counts in cylinders of different sizes $P(N_2/N_5)$. We find that despite the frequent use of the combination $w_{\mathrm{p}}(r_{\mathrm{p}})+\Delta \Sigma(r_{\mathrm{p}})$ in interpreting galaxy data, the count statistics, $P(N_{\mathrm{CIC}})$ and $P(N_{\mathrm{CIA}})$, are generally more efficient in constraining galaxy assembly bias when combined with $w_{\mathrm{p}}(r_{\mathrm{p}})$. Constraints based upon $w_{\mathrm{p}}(r_{\mathrm{p}})$ and $\Delta \Sigma(r_{\mathrm{p}})$ share common degeneracy directions in the parameter space, while combinations of $w_{\mathrm{p}}(r_{\mathrm{p}})$ with the count statistics are more complementary. Therefore, we strongly suggest that count statistics should be used to complement the canonical observables in future studies of the galaxy--halo connection.

We present a Bayesian hierarchical inference formalism (Basilisk) to constrain the galaxy-halo connection using satellite kinematics. Unlike traditional methods, Basilisk does not resort to stacking the kinematics of satellite galaxies in bins of central luminosity, and does not make use of summary statistics, such as satellite velocity dispersion. Rather, Basilisk leaves the data in its raw form and computes the corresponding likelihood. In addition, Basilisk can be applied to flux-limited, rather than volume-limited samples, greatly enhancing the quantity and dynamic range of the data. And finally, Basilisk is the only available method that simultaneously solves for halo mass and orbital anisotropy of the satellite galaxies, while properly accounting for scatter in the galaxy-halo connection. Basilisk uses the conditional luminosity function to model halo occupation statistics, and assumes that satellite galaxies are a relaxed tracer population of the host halo's potential with kinematics that obey the spherical Jeans equation. We test and validate Basilisk using mocks of varying complexity, and demonstrate that it yields unbiased constraints on the galaxy-halo connection and at a precision that rivals galaxy-galaxy lensing. In particular, Basilisk accurately recovers the full PDF of the relation between halo mass and central galaxy luminosity, and simultaneously constrains the orbital anisotropy of the satellite galaxies. Basilisk's inference is not affected by potential velocity bias of the central galaxies, or by slight errors in the inferred, radial profile of satellite galaxies that arise as a consequence of interlopers and sample impurity.

We present a new cosmological zoom-in simulation, where the zoom region consists of two halos with virial mass M_v~5x10^{12}M_{sun} and a ~Mpc long cosmic filament connecting them at z~2. Using this simulation, we study the evolution of the intergalactic medium in between these two halos at unprecedented resolution. At 5>z>3, the two halos are found to lie in a large intergalactic sheet, or "pancake", consisting of multiple co-planar dense filaments along which nearly all halos with M_v>10^9M_{sun} are located. This sheet collapses at z~5 from the merger of two smaller sheets. The strong shock generated by this merger leads to thermal instabilities in the post-shock region, and to a shattering of the sheet resulting in <~kpc scale clouds with temperatures of T>~2x10^4K and densities of n>~10^{-3}cm^{-3}, which are pressure confined in a hot medium with T~10^6K and n>~10^{-5}cm^{-3}. When the sheet is viewed face on, these cold clouds have neutral hydrogen column densities of N_{HI}>10^{17.2}cm^{-2}, making them detectable as Lyman limit systems, though they lie well outside the virial radius of any halo and even well outside the dense filaments. Their chemical composition is pristine, having zero metalicity, similar to several recently observed systems. Since these systems form far from any galaxies, these results are robust to galaxy formation physics, resulting purely from the collapse of large scale structure and radiative cooling, provided sufficient spatial resolution is available.

We investigate the abundance, small-scale clustering and galaxy-galaxy lensing signal of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). To this end, we present new measurements of the redshift and stellar mass dependence of the lensing properties of the galaxy sample. We analyse to what extent models assuming the Planck18 cosmology fit to the number density and clustering can accurately predict the small-scale lensing signal. In qualitative agreement with previous BOSS studies at redshift $z \sim 0.5$ and with results from the Sloan Digital Sky Survey, we find that the expected signal at small scales ($0.1 < r_{\rm p} < 3 \, h^{-1} \mathrm{Mpc}$) is higher by $\sim 25\%$ than what is measured. Here, we show that this result is persistent over the redshift range $0.1 < z < 0.7$ and for galaxies of different stellar masses. If interpreted as evidence for cosmological parameters different from the Planck CMB findings, our results imply $S_8 = \sigma_8 \sqrt{\Omega_{\rm m} / 0.3} = 0.744 \pm 0.015$, whereas $S_8 = 0.832 \pm 0.013$ for Planck18. However, in addition to being in tension with CMB results, such a change in cosmology alone does not accurately predict the lensing amplitude at larger scales. Instead, other often neglected systematics like baryonic feedback or assembly bias are likely contributing to the small-scale lensing discrepancy. We show that either effect alone, though, is unlikely to completely resolve the tension. Ultimately, a combination of the two effects in combination with a moderate change in cosmological parameters might be needed.

The dark matter content of the ultra diffuse galaxy NGC1052-DF2, as inferred from globular cluster (GC) and stellar kinematics, carries a considerable amount of uncertainty, with current constraints also allowing for the complete absence of dark matter. We test the viability of such a scenario by examining whether in a `baryon-only' mass model, the observed GC population experiences rapid orbital decay due to dynamical friction. Using a suite of 50 multi-GC N-body simulations that match observational constraints on both the stellar component of NGC1052-DF2 and its GC population but differ in the initial line-of-sight positions and the tangential velocities of the GCs, we show that there is a substantial amount of realization-to-realization variance in the evolution of the GCs. Nevertheless, over 10 Gyr, some of the GCs experience significant orbital evolution. Others evolve less. A combination of reduced dynamical friction in the galaxy core and GC-GC scattering keeps the GCs afloat, preventing them from sinking all the way to the galaxy center. While the current phase-space coordinates of the GCs are not unlikely for a baryon-only mass model, the GC system does evolve over time. Therefore, if NGC1052-DF2 has no dark matter, some of its GCs must have formed further out, and the GC system must have been somewhat more extended in the past. The presence of a low mass cuspy halo, while allowed by the kinematics, seems improbable as significantly shorter inspiral timescales in the central region would quickly lead to the formation of a nuclear star cluster.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological $N$-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that accurately describe subhalo evolution exist either. We therefore present a library of idealized, high resolution $N$-body simulations of the tidal evolution of individual subhaloes that can be used to calibrate semi-analytical models and to complement cosmological simulations. The simulations focus on minor mergers, i.e., the mass of the subhalo is much smaller than that of the host halo, such that the impact of dynamical friction is negligible. This setup allows the adoption of a fixed analytical potential for modelling the host halo. The dynamical evolution of subhaloes is followed with $N$-body computations. In the library, four parameters, two of which characterize the subhalo orbit with respect to the host halo, and the two concentrations of the host- and subhalo, are varied over the ranges encountered in cosmological simulations. We show several representative examples from the library that illustrate the evolution of the subhalo mass and velocity dispersion profiles. Additionally, we make publicly available a pre-trained non-parametric model of the subhalo mass evolution based on random forest regression. This model is able to interpolate the simulation data at the 0.1\,dex level and provides efficient access to the data for further use in modelling.

We propose that cold filamentary accretion in massive galaxies at high redshift can lead to the formation of star-forming clumps in the halos of these galaxies without the presence of dark matter sub-structure. In certain cases, these clumps can be the birth places of metal poor globular-clusters (MP GCs). Using cosmological simulations, we show that narrow streams of dense gas feeding massive galaxies from the cosmic web can fragment, producing star-forming clumps. We then derive an analytical model for the properties of streams as a function of halo mass and redshift, and assess when these are gravitationally unstable, when this can lead to collapse and star-formation in the halo, and when it may result in the formation of MP GCs. For stream metalicities >~1% solar, this is likely to occur at z>4.5. At z~6, the collapsing clouds have masses of (5-10)x10^7 M_{sun} and the average stream pressure is ~10^6 K cm^{-3}. The conditions for GC formation are met in the extremely turbulent "eyewall" at ~0.3 R_{vir}, where counter-rotating streams can collide, driving very large densities. Our scenario can account for the observed kinematics and spatial distribution of MP GCs, the correlation between their mass and metalicity, and the mass ratio between the GC system and the host halo. For MW mass halos, we infer that ~30% of MP GCs could have formed in this way, the remainder likely accreted in mergers. Our predictions for GC formation along circumgalactic filaments at high-redshift are testable with JWST.

Strong lensing of active galactic nuclei in the radio can result in razor-thin arcs, with a thickness of less than a milli-arcsecond, if observed at the resolution achievable with very long baseline interferometry (VLBI). Such razor-thin arcs provide a unique window on the coarseness of the matter distribution between source and observer. In this paper, we investigate to what extent such razor-thin arcs can constrain the number density and mass function of `free-floating' black holes, defined as black holes that do not, or no longer, reside at the centre of a galaxy. These can be either primordial in origin or arise as by-products of the evolution of super-massive black holes in galactic nuclei. When sufficiently close to the line of sight, free-floating black holes cause kink-like distortions in the arcs, which are detectable by eye in the VLBI images as long as the black hole mass exceeds $\sim 1000$ Solar masses. Using a crude estimate for the detectability of such distortions, we analytically compute constraints on the matter density of free-floating black holes resulting from null-detections of distortions along a realistic, fiducial arc, and find them to be comparable to those from quasar milli-lensing. We also use predictions from a large hydrodynamical simulation for the demographics of free-floating black holes that are not primordial in origin, and show that their predicted mass density is roughly four orders of magnitude below the constraints achievable with a single razor-thin arc.

We present new results on the relationship between central galaxies and dark matter haloes inferred from observations of satellite kinematics in the Sloan Digital Sky Survey (SDSS) DR7. We employ an updated analysis framework that includes detailed mock catalogues to model observational effects in SDSS. Our results constrain the colour-dependent conditional luminosity function (CLF) of dark matter haloes, as well as the radial profile of satellite galaxies. Confirming previous results, we find that red central galaxies live in more massive haloes than blue galaxies at fixed luminosity. Additionally, our results suggest that satellite galaxies have a radial profile less centrally concentrated than dark matter but not as cored as resolved subhaloes in dark matter-only simulations. Compared to previous works using satellite kinematics by More et al., we find much more competitive constraints on the galaxy-halo connection, on par with those derived from a combination of galaxy clustering and galaxy-galaxy lensing. We compare our results on the galaxy-halo connection to other studies using galaxy clustering and group catalogues, showing very good agreement between these different techniques. We discuss future applications of satellite kinematics in the context of constraining cosmology and the relationship between galaxies and dark matter haloes.

The kinematics of satellite galaxies moving in a dark matter halo are a direct probe of the underlying gravitational potential. Thus, the phase-space distributions of satellites represent a powerful tool to determine the galaxy-halo connection from observations. By stacking the signal of a large number of satellite galaxies this potential can be unlocked even for haloes hosting a few satellites on average. In this work, we test the impact of various modelling assumptions on constraints derived from analysing satellite phase-space distributions in the non-linear, 1-halo regime. We discuss their potential to explain the discrepancy between average halo masses derived from satellite kinematics and gravitational lensing previously reported. Furthermore, we develop an updated, more robust analysis to extract constraints on the galaxy-halo relation from satellite properties in spectroscopic galaxy surveys such as the SDSS. We test the accuracy of this approach using a large number of realistic mock catalogues. Furthermore, we find that constraints derived from such an analysis are complementary and competitive with respect to the commonly used galaxy clustering and galaxy-galaxy lensing observables.

The circumgalactic medium (CGM), i.e. the gaseous haloes around galaxies, is both the reservoir of gas that fuels galaxy growth and the repository of gas expelled by galactic winds. Most cosmological, hydrodynamical simulations focus their computational effort on the galaxies themselves and treat the CGM more coarsely, which means small-scale structure cannot be resolved. We get around this issue by running zoom-in simulations of a Milky Way-mass galaxy with standard mass refinement and additional uniform spatial refinement within the virial radius. This results in a detailed view of its gaseous halo at unprecedented (1 kpc) uniform resolution with only a moderate increase in computational time. The improved spatial resolution does not impact the central galaxy or the average density of the CGM. However, it drastically changes the radial profile of the neutral hydrogen column density, which is enhanced at galactocentric radii larger than 40 kpc. The covering fraction of Lyman-Limit Systems within 150 kpc is almost doubled. We therefore conclude that some of the observational properties of the CGM are strongly resolution dependent. Increasing the resolution in the CGM, without increasing the resolution of the galaxies, is a promising and computationally efficient method to push the boundaries of state-of-the-art simulations.

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 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 investigate how Einstein rings and magnified arcs are affected by small-mass dark-matter haloes placed along the line-of-sight to gravitational lens systems. By comparing the gravitational signature of line-of-sight haloes with that of substructures within the lensing galaxy, we derive a mass-redshift relation that allows us to rescale the detection threshold (i.e. lowest detectable mass) for substructures to a detection threshold for line-of-sight haloes at any redshift. We then quantify the line-of-sight contribution to the total number density of low-mass objects that can be detected through strong gravitational lensing. Finally, we assess the degeneracy between substructures and line-of-sight haloes of different mass and redshift to provide a statistical interpretation of current and future detections, with the aim of distinguishing between CDM and WDM. We find that line-of-sight haloes statistically dominate with respect to substructures, by an amount that strongly depends on the source and lens redshifts, and on the chosen dark matter model. Substructures represent about 30 percent of the total number of perturbers for low lens and source redshifts (as for the SLACS lenses), but less than 10 per cent for high redshift systems. We also find that for data with high enough signal-to-noise ratio and angular resolution, the non-linear effects arising from a double-lens-plane configuration are such that one is able to observationally recover the line-of-sight halo redshift with an absolute error precision of 0.15 at the 68 per cent confidence level.

To gain understanding of the complicated, non-linear and numerical processes associated with the tidal evolution of dark matter subhaloes in numerical simulation, we perform a large suite of idealized simulations that follow individual N-body subhaloes in a fixed, analytical host halo potential. By varying both physical and numerical parameters, we investigate under what conditions the subhaloes undergo disruption. We confirm the conclusions from our more analytical assessment in van den Bosch et al. that most disruption is numerical in origin; as long as a subhalo is resolved with sufficient mass and force resolution, a bound remnant survives. This implies that state-of-the-art cosmological simulations still suffer from significant overmerging. We demonstrate that this is mainly due to inadequate force-softening, which causes excessive mass loss and artificial tidal disruption. In addition, we show that subhaloes in N-body simulations are susceptible to a runaway instability triggered by the amplification of discreteness noise in the presence of a tidal field. These two processes conspire to put serious limitations on the reliability of dark matter substructure in state-of-the-art cosmological simulations. We present two criteria that can be used to assess whether individual subhaloes in cosmological simulations are reliable or not, and advocate that subhaloes that satisfy either of these two criteria be discarded from further analysis. We discuss the potential implications of this work for several areas in astrophysics.

Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight velocity and spatial offsets between brightest galaxies and their neighbours. We compare those to detailed mock catalogues, constructed from high-resolution, dark-matter-only $N$-body simulations, in which it is assumed that satellite galaxies trace dark matter subhaloes. This allows us to place constraints on the fraction $f_{\rm BNC}$ of haloes in which the brightest galaxy is not the central. Compared to previous studies we explicitly take into account the unrelaxed state of the host haloes, velocity offsets of halo cores and correlations between $f_{\rm BNC}$ and the satellite occupation. We find that $f_{\rm BNC}$ strongly decreases with the luminosity of the brightest galaxy and increases with the mass of the host halo. Overall, in the halo mass range $10^{13} - 10^{14.5} h^{-1} M_\odot$ we find $f_{\rm BNC} \sim 30\%$, in good agreement with a previous study by Skibba et al. We discuss the implications of these findings for studies inferring the galaxy--halo connection from satellite kinematics, models of the conditional luminosity function and galaxy formation in general.

We study 379 central and 159 satellite early-type galaxies with two-dimensional kinematics from the integral-field survey Mapping Nearby Galaxies at APO (MaNGA) to determine how their angular momentum content depends on stellar and halo mass. Using the Yang et. al. (2007) group catalog, we identify central and satellite galaxies in groups with halo masses in the range 10^12.5 h^-1 M_sun < M_200b < 10^15 h^-1 M_sun. As in previous work, we see a sharp dependence on stellar mass, in the sense that ~ 70% of galaxies with stellar mass M_* > 10^11 h^-2 M_sun tend to have very little rotation, while nearly all galaxies at lower mass show some net rotation. The ~ 30% of high-mass galaxies that have significant rotation do not stand out in other galaxy properties except for a higher incidence of ionized gas emission. Our data are consistent with recent simulation results suggesting that major merging and gas accretion have more impact on the rotational support of lower-mass galaxies. When carefully matching the stellar mass distributions, we find no residual differences in angular momentum content between satellite and central galaxies at the 20\% level. Similarly, at fixed mass, galaxies have consistent rotation properties across a wide range of halo mass. However, we find that errors in classification of centrals and satellites with group finders systematically lowers differences between satellite and central galaxies at a level that is comparable to current measurement uncertainties. To improve constraints, the impact of group finding methods will have to be forward modeled via mock catalogs.

Accurately predicting the demographics of dark matter (DM) substructure is of paramount importance for many fields of astrophysics, including gravitational lensing, galaxy evolution, halo occupation modeling, and constraining the nature of dark matter. Because of its strongly non-linear nature, DM substructure is typically modeled using N-body simulations, which reveal that large fractions of DM subhaloes undergo complete disruption. In this paper we use both analytical estimates and idealized numerical simulations to investigate whether this disruption is mainly physical, due to tidal heating and stripping, or numerical (i.e., artificial). We show that, contrary to naive expectation, subhaloes that experience a tidal shock $\Delta E$ that exceeds the subhalo's binding energy, $|E_{\rm b}|$, do not undergo disruption, even when $\Delta E/|E_{\rm b}|$ is as large as 100. Along the same line, and contrary to existing claims in the literature, instantaneously stripping matter from the outskirts of a DM subhalo also does not result in its complete disruption, even when the instantaneous remnant has positive binding energy. In addition, we show that tidal heating due to high-speed (impulsive) encounters with other subhaloes (`harassment'), is negligible compared to the tidal effects due to the host halo. Hence, we conclude that, in the absence of baryonic processes, the complete, physical disruption of CDM substructure is extremely rare, and that most disruption in numerical simulations therefore must be artificial. We discuss various processes that have been associated with numerical overmerging, and conclude that inadequate force-softening is the most likely culprit.

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.

Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $\Delta$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $\Delta$. For conventional halo definitions ($\Delta \sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $\Delta \sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $\Delta$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $\Delta \gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass.

We estimate the stellar mass that satellite galaxies lose once they enter groups (and clusters) by identifying groups in a high-resolution cosmological N-body simulation, assigning entry masses to satellite galaxies with halo abundance matching at the entry time, and comparing the predicted conditional stellar mass function of satellite galaxies at $z\simeq0$ with observations. Our results depend on the mass of the stars that form in satellite galaxies after the entry time. A model in which star formation shuts down completely as soon a galaxy enters a group environment is ruled out because it underpredicts the stellar masses of satellite galaxies even in the absence of tidal stripping. The greater is the stellar mass that is allowed to form, the greater the fraction that needs to be tidally stripped. The stellar mass fraction lost by satellite galaxies after entering a group or cluster environment is consistent with any value in the range $0-40\%$. To place stronger constraints, we consider a more refined model of tidal stripping of galaxies on elongated orbits (where stripping occurs at orbit pericentres). Our model predicts less tidal stripping: satellite galaxies lose $\sim 20-25\%$ of their stellar mass since their entry into the group. This finding is consistent with a slow-starvation delayed-quenching picture, in which galaxies that enter a group or cluster environment keep forming stars until at least the first pericentric passage.

Galaxy clustering on small scales is significantly under-predicted by sub-halo abundance matching (SHAM) models that populate (sub-)haloes with galaxies based on peak halo mass, $M_{\rm peak}$. SHAM models based on the peak maximum circular velocity, $V_{\rm peak}$, have had much better success. The primary reason $M_{\rm peak}$ based models fail is the relatively low abundance of satellite galaxies produced in these models compared to those based on $V_{\rm peak}$. Despite success in predicting clustering, a simple $V_{\rm peak}$ based SHAM model results in predictions for galaxy growth that are at odds with observations. We evaluate three possible remedies that could "save" mass-based SHAM: (1) SHAM models require a significant population of "orphan" galaxies as a result of artificial disruption/merging of sub-haloes in modern high resolution dark matter simulations; (2) satellites must grow significantly after their accretion; and (3) stellar mass is significantly affected by halo assembly history. No solution is entirely satisfactory. However, regardless of the particulars, we show that popular SHAM models based on $M_{\rm peak}$ cannot be complete physical models as presented. Either $V_{\rm peak}$ truly is a better predictor of stellar mass at $z\sim 0$ and it remains to be seen how the correlation between stellar mass and $V_{\rm peak}$ comes about, or SHAM models are missing vital component(s) that significantly affect galaxy clustering.

We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations.

We present a comprehensive analysis of the evolution of dark matter subhaloes in the cosmological Bolshoi simulation. We identify a complete set of 12 unique evolution channels by which subhaloes evolve in between simulation outputs, and study their relative importance and demographics. We show that instantaneous masses and maximum circular velocities of individual subhaloes are extremely noisy, despite the use of a sophisticated, phase-space-based halo finder. We also show that subhaloes experience frequent penetrating encounters with other subhaloes (on average about one per dynamical time), and that subhaloes whose apo-center lies outside the virial radius of their host (the 'ejected' or 'backsplash' haloes) experience tidal forces that modify their orbits. This results in an average fractional subhalo exchange rate among host haloes of roughly 0.01 per Gyr (at the present time). In addition, we show that there are three distinct disruption channels; one in which subhaloes drop below the mass resolution limit of the simulation, one in which subhaloes merge with their host halo largely driven by dynamical friction, and one in which subhaloes abruptly disintegrate. We estimate that roughly 80 percent of all subhalo disruption in the Bolshoi simulation is numerical, rather than physical. This over-merging is a serious road-block for the use of numerical simulations to interpret small scale clustering, or for any other study that is sensitive to the detailed demographics of dark matter substructure.

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

The baryonic Tully-Fisher relation (BTFR) is a valuable observational tool and a critical test of galaxy formation theory. We explore the systematic uncertainty in the slope and the scatter of the observed line width BTFR utilizing homogeneously measured, unresolved HI observations for 930 isolated galaxies. We measure a fiducial relation of $\log_{10}{M_{\rm baryon}} = 3.24\log_{10}{V_{\rm rot}}+3.21$ with observed scatter of 0.25 dex with $10^{7.4} < M_{\rm baryon} < 10^{11.3}M_{\odot}$ where $V_{\rm rot}$ is measured from 20\% HI line widths. We vary the definitions of $M_{\rm baryon}$ and $V_{\rm rot}$, the sample selection, and the linear fitting algorithm. We obtain slopes ranging from 2.64 to 3.53 and scatter measurements ranging from 0.14 to 0.41~dex, indicating a systematic uncertainty of 0.25 in the BTFR slope derived from unresolved HI line widths. We compare our fiducial slope to literature measurements, where reported slopes range from 3.0 to 4.3 and scatter is either unmeasured, immeasurable or as large as 0.4~dex. Measurements derived from unresolved HI line widths tend to produce slopes of 3.3, while measurements derived from asymptotic rotation curves tend to produce slopes of 3.9. The largest factor affecting the BTFR slope is the definition of $V_{\rm rot}$. The sample definition, the mass range, and the linear fitting algorithm also significantly affect the measured BTFR. We find that galaxies in our sample with $V_{\rm rot} < 100~\rm{km~s^{-1}}$ are consistent with the line width BTFR of more massive galaxies, but these galaxies drive most of the observed scatter. It is critical when comparing predictions to an observed BTFR that the $V_{\rm rot}$ definition, the sample selection and the fitting algorithm are similarly defined. We recommend direct statistical comparisons between data sets with commensurable properties as opposed to comparing BTFR power-law fits.

We present a study of unprecedented statistical power regarding the halo-to-halo variance of dark matter substructure. Using a combination of N-body simulations and a semi-analytical model, we investigate the variance in subhalo mass fractions and subhalo occupation numbers, with an emphasis on how these statistics scale with halo formation time. We demonstrate that the subhalo mass fraction, f_sub, is mainly a function of halo formation time, with earlier forming haloes having less substructure. At fixed formation redshift, the average f_sub is virtually independent of halo mass, and the mass dependence of f_sub is therefore mainly a manifestation of more massive haloes assembling later. We compare observational constraints on f_sub from gravitational lensing to our model predictions and simulation results. Although the inferred f_sub are substantially higher than the median LCDM predictions, they fall within the 95th percentile due to halo-to-halo variance. We show that while the halo occupation distribution of subhaloes, P(N|M), is super-Poissonian for large <N>, a well established result, it becomes sub-Poissonian for <N> < 2. Ignoring the non-Poissonity results in systematic errors of the clustering of galaxies of a few percent, and with a complicated scale- and luminosity-dependence. Earlier-formed haloes have P(N|M) closer to a Poisson distribution, suggesting that the dynamical evolution of subhaloes drives the statistics towards Poissonian. Contrary to a recent claim, the non-Poissonity of subhalo occupation statistics does not vanish by selecting haloes with fixed mass and fixed formation redshift. Finally, we use subhalo occupation statistics to put loose constraints on the mass and formation redshift of the Milky Way halo. Using observational constraints on the V_max of the most massive satellites, we infer that 0.25<M_vir/10^12M_sun/h<1.4 and 0.1<z_f<1.4 at 90% confidence.

We propose that one of the sources in the recently detected system CR7 by Sobral et al. (2015) through spectro-photometric measurements at $z = 6.6$ harbors a direct collapse blackhole (DCBH). We argue that the LW radiation field required for direct collapse in source A is provided by sources B and C. By tracing the LW production history and star formation rate over cosmic time for the halo hosting CR7 in a $\Lambda$CDM universe, we demonstrate that a DCBH could have formed at $z\sim 20$. The spectrum of source A is well fit by nebular emission from primordial gas around a BH with MBH $\sim 4.4 \times 10^6 \ \rm M_{\odot}$ accreting at a 40 % of the Eddington rate, which strongly supports our interpretation of the data. Combining these lines of evidence, we argue that CR7 might well be the first DCBH candidate.

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 constrain the newly-introduced decorated Halo Occupation Distribution (HOD) model using SDSS DR7 measurements of projected galaxy clustering or r-band luminosity threshold samples. The decorated HOD is a model for the galaxy-halo connection that augments the HOD by allowing for the possibility of galaxy assembly bias: galaxy luminosity may be correlated with dark matter halo properties besides mass, Mvir. We demonstrate that it is not possible to rule out galaxy assembly bias using DR7 measurements of galaxy clustering alone. Moreover, galaxy samples with Mr < -20 and Mr < -20.5 favor strong central galaxy assembly bias. These samples prefer scenarios in which high-concentration are more likely to host a central galaxy relative to low-concentration halos of the same mass. We exclude zero assembly bias with high significance for these samples. Satellite galaxy assembly bias is significant for the faintest sample, Mr < -19. We find no evidence for assembly bias in the Mr < -21 sample. Assembly bias should be accounted for in galaxy clustering analyses or attempts to exploit galaxy clustering to constrain cosmology. In addition to presenting the first constraints on HOD models that accommodate assembly bias, our analysis includes several improvements over previous analyses of these data. Therefore, our inferences supersede previously-published results even in the case of a standard HOD analysis.

The connection between galaxies and dark matter halos is often inferred from data using probabilistic models, such as the Halo Occupation Distribution (HOD). Conventional HOD formulations assume that only halo mass governs the galaxy-halo connection. Violations of this assumption, known as galaxy assembly bias, threaten the HOD program. We introduce decorated HODs, a new, flexible class of models designed to account for assembly bias. Decorated HODs minimally expand the parameter space and maximize the independence between traditional and novel HOD parameters. We use decorated HODs to quantify the influence of assembly bias on clustering and lensing statistics. For SDSS-like samples, the impact of assembly bias on galaxy clustering can be as large as a factor of two on r ~ 200 kpc scales and ~15% in the linear regime. Assembly bias can either enhance or diminish clustering on large scales, but generally increases clustering on scales r <~ 1 Mpc. We performed our calculations with Halotools, an open-source, community-driven python package for studying the galaxy-halo connection (http://halotools.readthedocs.org). We conclude by describing the use of decorated HODs to treat assembly bias in otherwise conventional likelihood analyses.

We present the first comprehensive analysis of the segregation of dark matter subhaloes in their host haloes. Using numerical simulations, we examine the segregation of twelve different subhalo properties with respect to both orbital energy and halo-centric radius (in real space as well as in projection). Subhaloes are strongly segregated by accretion redshift, which is an outcome of the inside-out assembly of their host haloes. Since subhaloes that were accreted earlier have experienced more tidal stripping, subhaloes that have lost a larger fraction of their mass at infall are on more bound orbits. Subhaloes are also strongly segregated in their masses and maximum circular velocities at accretion. We demonstrate that part of this segregation is already imprinted in the infall conditions. For massive subhaloes it is subsequently boosted by dynamical friction, but only during their first radial orbit. The impact of these two effects is counterbalanced, though, by the fact that subhaloes with larger accretion masses are accreted later. Because of tidal stripping, subhaloes reveal little to no segregation by present-day mass or maximum circular velocity, while the corresponding torques cause subhaloes on more bound orbits to have smaller spin. There is a weak tendency for subhaloes that formed earlier to be segregated towards the center of their host halo, which is an indirect consequence of the fact that (sub)halo formation time is correlated with other, strongly segregated properties. We discuss the implications of our results for the segregation of satellite galaxies in galaxy groups and clusters.

Large surveys have shown that red galaxies are preferentially aligned with their halos while blue galaxies have a more isotropic distribution. Since halos generally align with their filaments this introduces a bias in the measurement of the cosmic shear from weak lensing. It is therefore vitally important to understand why this difference arises. We explore the stability of different disc orientations within triaxial halos. We show that, in the absence of gas, the disc orientation is most stable when its spin is along the minor axis of the halo. Instead when gas cools onto a disc it is able to form in almost arbitrary orientation, including off the main planes of the halo (but avoiding an orientation perpendicular to the halo's intermediate axis). Substructure helps gasless galaxies reach alignment with the halo faster, but have less effect on galaxies when gas is cooling onto the disc. Our results provide a novel and natural interpretation for why red, gas poor galaxies are preferentially aligned with their halo, while blue, star-forming, galaxies have nearly random orientations, without requiring a connection between galaxies' current star formation rate and their merger history.

We use a semi-analytical model for the substructure of dark matter haloes to assess the too-big-to-fail (TBTF) problem. The model accurately reproduces the average subhalo mass and velocity functions, as well as their halo-to-halo variance, in N-body simulations. We construct thousands of realizations of Milky Way (MW) size host haloes, allowing us to investigate the TBTF problem with unprecedented statistical power. We examine the dependence on host halo mass and cosmology, and explicitly demonstrate that a reliable assessment of TBTF requires large samples of hundreds of host haloes. We argue that previous statistics used to address TBTF suffer from the look-elsewhere effect and/or disregard certain aspects of the data on the MW satellite population. We devise a new statistic that is not hampered by these shortcomings, and, using only data on the 9 known MW satellite galaxies with $V_{\rm max}>15{\rm kms}^{-1}$, demonstrate that $1.4^{+3.3}_{-1.1}\%$ of MW-size host haloes have a subhalo population in statistical agreement with that of the MW. However, when using data on the MW satellite galaxies down to $V_{\rm max}=8{\rm kms}^{-1}$, this MW consistent fraction plummets to $<5\times10^{-4}$ (at 68% CL). Hence, if it turns out that the inventory of MW satellite galaxies is complete down to 8km/s, then the maximum circular velocities of MW satellites are utterly inconsistent with $\Lambda$CDM predictions, unless baryonic effects can drastically increase the spread in $V_{\rm max}$ values of satellite galaxies compared to that of their subhaloes.

A $^{238}$U projectile beam was used to create cadmium isotopes via abrasion-fission at 410 MeV/u in a beryllium target at the entrance of the in-flight separator FRS at GSI. The fission fragments were separated with the FRS and injected into the isochronous storage ring ESR for mass measurements. The Isochronous Mass Spectrometry (IMS) was performed under two different experimental conditions, with and without B$\rho$-tagging at the dispersive central focal plane of the FRS. In the experiment with B$\rho$-tagging the magnetic rigidity of the injected fragments was determined by an accuracy of $2\times 10^{-4}$. A new method of data analysis, using a correlation matrix for the combined data set from both experiments, has provided mass values for 25 different isotopes for the first time. The high selectivity and sensitivity of the experiment and analysis has given access even to rare isotopes detected with a few atoms per week. In this letter we present for the $^{129,130,131}$Cd isotopes mass values directly measured for the first time. The Cd results clearly show a very pronounced shell effect at $N=82$ which is in agreement with the conclusion from $\gamma$-ray spectroscopy of $^{130}$Cd and confirms the assumptions of modern shell-model calculations.

We present Ha maps at 1kpc spatial resolution for star-forming galaxies at z~1, made possible by the WFC3 grism on HST. Employing this capability over all five 3D-HST/CANDELS fields provides a sample of 2676 galaxies. By creating deep stacked Halpha (Ha) images, we reach surface brightness limits of 1x10^-18\erg\s\cm^2\arcsec^2, allowing us to map the distribution of ionized gas out to >10kpc for typical L* galaxies at this epoch. We find that the spatial extent of the Ha distribution increases with stellar mass as r(Ha)[kpc]=1.5(Mstars/10^10Msun)^0.23. Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out assembly of galactic disks. This effect, however, is mass dependent with r(Ha)/r(stars)=1.1(M/10^10Msun)^0.054, such that at low masses r(Ha)~r(stars). We map the Ha distribution as a function of SFR(IR+UV) and find evidence for `coherent star formation' across the SFR-M plane: above the main sequence, Ha is enhanced at all radii; below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at fixed mass. At high masses (10^10.5<M/Msun<10^11), above the main sequence, Ha is particularly enhanced in the center, plausibly building bulges and/or supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across the entirety of the SFR-M plane, the absolute star formation rate as traced by Ha is always centrally peaked, even in galaxies below the main sequence.

This work presents a direct measurement of the $^{96}$Ru($p, \gamma$)$^{97}$Rh cross section via a novel technique using a storage ring, which opens opportunities for reaction measurements on unstable nuclei. A proof-of-principle experiment was performed at the storage ring ESR at GSI in Darmstadt, where circulating $^{96}$Ru ions interacted repeatedly with a hydrogen target. The $^{96}$Ru($p, \gamma$)$^{97}$Rh cross section between 9 and 11 MeV has been determined using two independent normalization methods. As key ingredients in Hauser-Feshbach calculations, the $\gamma$-ray strength function as well as the level density model can be pinned down with the measured ($p, \gamma$) cross section. Furthermore, the proton optical potential can be optimized after the uncertainties from the $\gamma$-ray strength function and the level density have been removed. As a result, a constrained $^{96}$Ru($p, \gamma$)$^{97}$Rh reaction rate over a wide temperature range is recommended for $p$-process network calculations.

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.

Correlations between the star formation rates (SFRs) of nearby galaxies (so-called galactic conformity) have been observed for projected separations up to 4 Mpc, an effect not predicted by current semi-analytic models. We investigate correlations between the mass accretion rates (dMvir/dt) of nearby halos as a potential physical origin for this effect. We find that pairs of host halos "know about" each others' assembly histories even when their present-day separation is greater than thirty times the virial radius of either halo. These distances are far too large for direct interaction between the halos to explain the correlation in their dMvir/dt. Instead, halo pairs at these distances reside in the same large-scale tidal environment, which regulates dMvir/dt for both halos. Larger halos are less affected by external forces, which naturally gives rise to a mass dependence of the halo conformity signal. SDSS measurements of galactic conformity exhibit a qualitatively similar dependence on stellar mass, including how the signal varies with distance. Based on the expectation that halo accretion and galaxy SFR are correlated, we predict the scale-, mass- and redshift-dependence of large-scale galactic conformity, finding that the signal should drop to undetectable levels by z > 1. These predictions are testable with current surveys to z ~ 1; confirmation would establish a strong correlation between dark matter halo accretion rate and central galaxy SFR.

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.

It is common practice for methods that use galaxy clustering to constrain the galaxy-halo relationship, such as the halo occupation distribution (HOD) and/or conditional luminosity function (CLF), to assume that halo mass alone suffices to determine a halo's resident galaxy population. Yet, the clustering strength of cold dark matter halos depends upon halo properties in addition to mass, such as formation time, an effect referred to as assembly bias. If galaxy characteristics are correlated with any of these auxiliary halo properties, the basic assumption of HOD/CLF methods is violated. We estimate the potential for assembly bias to induce systematic errors in inferred halo occupation statistics. We use halo abundance matching and age matching to construct fiducial mock galaxy catalogs that exhibit assembly bias as well as additional mock catalogs with identical HODs, but with assembly bias removed. We fit a parameterized HOD to the projected two-point clustering of mock galaxies in each catalog to assess the systematic errors induced by reasonable levels of assembly bias. In the absence of assembly bias, the inferred HODs generally describe the true underlying HODs well, validating the basic methodology. However, in all of the cases with assembly bias, the inferred HODs have systematic errors that are statistically significant. In most cases, these systematic errors cannot be identified using void statistics as auxiliary observables. We conclude that the galaxy-halo relationship inferred from galaxy clustering should be subject to a non-negligible systematic error induced by assembly bias. Our work suggests that efforts to model and/or constrain assembly bias should be high priorities as it is a threatening source of systematic error in galaxy evolution studies as well as the precision cosmology program.

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.

In this paper, we test the age matching hypothesis that the star formation rate (SFR) of a galaxy of fixed stellar mass is determined by its dark matter halo formation history, and as such, that more quiescent galaxies reside in older halos. This simple model has been remarkably successful at predicting color-based galaxy statistics at low redshift as measured in the Sloan Digital Sky Survey (SDSS). To further test this method with observations, we present new SDSS measurements of the galaxy two-point correlation function and galaxy-galaxy lensing as a function of stellar mass and SFR, separated into quenched and star-forming galaxy samples. We find that our age matching model is in excellent agreement with these new measurements. We also employ a galaxy group finder and show that our model is able to predict: (1) the relative SFRs of central and satellite galaxies, (2) the SFR-dependence of the radial distribution of satellite galaxy populations within galaxy groups, rich groups, and clusters and their surrounding larger scale environments, and (3) the interesting feature that the satellite quenched fraction as a function of projected radial distance from the central galaxy exhibits an ~ r^-.15 slope, independent of environment. The accurate prediction for the spatial distribution of satellites is intriguing given the fact that we do not explicitly model satellite-specific processes after infall, and that in our model the virial radius does not mark a special transition region in the evolution of a satellite, contrary to most galaxy evolution models. The success of the model suggests that present-day galaxy SFR is strongly correlated with halo mass assembly history.

We present a detailed study of how dark matter haloes assemble their mass and grow their (central) potential well. We characterize these via their mass accretion histories (MAHs) and potential well growth histories (PWGHs), which we extract from the Bolshoi simulation and from semi-analytical merger trees supplemented with a method to compute the maximum circular velocity, Vmax, of progenitor haloes. The results of both methods are in excellent agreement, both in terms of the average and the scatter. We show that the MAH and PWGH are tightly correlated, and that growth of the central potential precedes the assembly of mass; the maximum circular velocity is already half the present day value by the time the halo has accreted only 2 percent of its final mass. Finally, we demonstrate that MAHs have a universal form, which we use to develop a new and improved universal model that can be used to compute the average or median MAH and PWGH for a halo of any mass in any LCDM cosmology, without having to run a numerical simulation or a set of halo merger trees.

We use the Sloan Digital Sky Survey to investigate the properties of massive elliptical galaxies in the local Universe (z\leq0.08) that have unusually blue optical colors. Through careful inspection, we distinguish elliptical from non-elliptical morphologies among a large sample of similarly blue galaxies with high central light concentrations (c_r\geq2.6). These blue ellipticals comprise 3.7 per cent of all c_r\geq2.6 galaxies with stellar masses between 10^10 and 10^11 h^{-2} {\rm M}_{\sun}. Using published fiber spectra diagnostics, we identify a unique subset of 172 non-star-forming ellipticals with distinctly blue urz colors and young (< 3 Gyr) light-weighted stellar ages. These recently quenched ellipticals (RQEs) have a number density of 2.7-4.7\times 10^{-5}\,h^3\,{\rm Mpc}^{-3} and sufficient numbers above 2.5\times10^{10} h^{-2} {\rm M}_{\sun} to account for more than half of the expected quiescent growth at late cosmic time assuming this phase lasts 0.5 Gyr. RQEs have properties that are consistent with a recent merger origin (i.e., they are strong `first-generation' elliptical candidates), yet few involved a starburst strong enough to produce an E+A signature. The preferred environment of RQEs (90 per cent reside at the centers of < 3\times 10^{12}\,h^{-1}{\rm M}_{\sun} groups) agrees well with the `small group scale' predicted for maximally efficient spiral merging onto their halo center and rules out satellite-specific quenching processes. The high incidence of Seyfert and LINER activity in RQEs and their plausible descendents may heat the atmospheres of small host halos sufficiently to maintain quenching.

Quenched central galaxies tend to reside in a preferentially quenched large-scale environment, a phenomenon that has been dubbed galactic conformity. Remarkably, this tendency persists out to scales far larger than the virial radius of the halo hosting the central. Therefore, conformity manifestly violates the widely adopted assumption that the dark matter halo mass Mvir exclusively governs galaxy occupation statistics. This paper is the first in a series studying the implications of the observed conformity signal for the galaxy-dark matter connection. We show that recent measurements of conformity on scales R ~1-5 Mpc imply that central galaxy quenching statistics cannot be correctly predicted with the knowledge of Mvir alone. We also demonstrate that ejected (or `backsplash') satellites cannot give rise to the signal. We then invoke the age matching model, which is predicated on the co-evolution of galaxies and halos. We find that this model produces a strong signal, and that central galaxies are solely responsible. We conclude that large-scale `2-halo' conformity represents a smoking gun of central galaxy assembly bias, and indicates that contemporary models of satellite quenching have systematically over-estimated the influence of post-infall processes.

We compare subhalo mass and velocity functions obtained from different simulations with different subhalo finders among each other, and with predictions from the new semi-analytical model of Jiang & van den Bosch (2014). We find that subhalo mass functions (SHMFs) obtained using different subhalo finders agree with each other at the level of ~ 20 percent, but only at the low mass end. At the massive end, subhalo finders that identify subhaloes based purely on density in configuration space dramatically underpredict the subhalo abundances by more than an order of magnitude. These problems are much less severe for subhalo velocity functions (SHVFs), indicating that they arise from issues related to assigning masses to the subhaloes, rather than from detecting them. Overall the predictions from the semi-analytical model are in excellent agreement with simulation results obtained using the more advanced subhalo finders that use information in six dimensional phase-space. In particular, the model accurately reproduces the slope and host-mass-dependent normalization of both the subhalo mass and velocity functions. We find that the SHMFs and SHVFs have power-law slopes of 0.82 and 2.6, respectively, significantly shallower than what has been claimed in several studies in the literature.

We present a new, semi-analytical model describing the evolution of dark matter subhaloes. The model uses merger trees constructed using the method of Parkinson et al. (2008) to describe the masses and redshifts of subhaloes at accretion, which are subsequently evolved using a simple model for the orbit-averaged mass loss rates. The model is extremely fast, treats subhaloes of all orders, accounts for scatter in orbital properties and halo concentrations, and uses a simple recipe to convert subhalo mass to maximum circular velocity. The model accurately reproduces the average subhalo mass and velocity functions in numerical simulations. The inferred subhalo mass loss rates imply that an average dark matter subhalo loses in excess of 80 percent of its infall mass during its first radial orbit within the host halo. We demonstrate that the total mass fraction in subhaloes is tightly correlated with the `dynamical age' of the host halo, defined as the number of halo dynamical times that have elapsed since its formation. Using this relation, we present universal fitting functions for the evolved and unevolved subhalo mass and velocity functions that are valid for any host halo mass, at any redshift, and for any {\Lambda}CDM cosmology.

Galaxies that are several virial radii beyond groups/clusters show preferentially quiescent star formation rates. Using a galaxy group/cluster catalog from the Sloan Digital Sky Survey, together with a cosmological N-body simulation, we examine the origin of this environmental quenching beyond the virial radius. Accounting for the clustering of groups/clusters, we show that central galaxies show enhanced SFR quenching out to 2.5 virial radii beyond groups/clusters, and we demonstrate that this extended environmental enhancement can be explained simply by 'ejected' satellite galaxies that orbit beyond their host halo's virial radius. We show that ejected satellites typically orbit for several Gyr beyond the virial radius before falling back in, and thus they compose up to 40% of all central galaxies near groups/clusters. We show that a model in which ejected satellites experience the same SFR quenching as satellites within a host halo can explain essentially all environmental dependence of galaxy quenching. Furthermore, ejected satellites (continue to) lose significant halo mass, an effect that is potentially observable via gravitational lensing. The SFRs/colors and stellar-to-halo masses of ejected satellites highlight the importance of environmental history and present challenges to models of galaxy occupation that ignore such history.

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.

Halo merger trees describe the hierarchical mass assembly of dark matter haloes, and are the backbone for modeling galaxy formation and evolution. Merger trees constructed using Monte Carlo algorithms based on the extended Press-Schechter (EPS) formalism are complementary to those extracted from N-body simulations, and have the advantage that they are not trammeled by limited numerical resolution and uncertainties in identifying (sub)haloes and linking them between snapshots. This paper compares multiple EPS-based merger tree algorithms to simulation results using four diagnostics: progenitor mass function (PMF), mass assembly history (MAH), merger rate per descendant halo, and the unevolved subhalo mass function (USMF). In general, algorithms based on spherical collapse yield major-merger rates that are too high by a factor of two, resulting in MAHs that are systematically offset. Assuming ellipsoidal collapse solves most of these issues, but the particular algorithm investigated here that incorporates ellipsoidal collapse dramatically overpredicts the minor-merger rate for massive haloes. The only algorithm in our comparison that yields MAHs, merger rates, and USMFs in good agreement with simulations, is that by Parkinson et al. (2008). However, this is not a true EPS-based algorithm as it draws its progenitor masses from a PMF calibrated against simulations, rather than `predicted' by EPS. Finally we emphasize that the benchmarks used to test the EPS algorithms are obtained from simulations and are hampered by significant uncertainties themselves. In particular, MAHs and halo merger rates obtained from simulations by different authors reveal discrepancies that easily exceed 50 percent, even when based on the same simulation. Given this status quo, merger trees constructed using the Parkinson et al. algorithm are as accurate as those extracted from N-body simulations.

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]

Satellite galaxies in groups and clusters are more likely to have low star formation rates (SFR) and lie on the red-sequence than central (field) galaxies. Using galaxy group/cluster catalogs from SDSS DR7, together with a cosmological N-body simulation to track satellite orbits, we examine the star formation histories and quenching timescales of satellites of M_star > 5 x 10^9 M_sun at z=0. We first explore satellite infall histories: group preprocessing and ejected orbits are critical aspects of satellite evolution, and properly accounting for these, satellite infall typically occurred at z~0.5, or ~5 Gyr ago. To obtain accurate initial conditions for the SFRs of satellites at their time of first infall, we construct an empirical parametrization for the evolution of central galaxy SFRs and quiescent fractions. With this, we constrain the importance and efficiency of satellite quenching as a function of satellite and host halo mass, finding that satellite quenching is the dominant process for building up all quiescent galaxies at M_star < 10^10 M_sun. We then constrain satellite star formation histories, finding a 'delayed-then-rapid' quenching scenario: satellite SFRs evolve unaffected for 2-4 Gyr after infall, after which star formation quenches rapidly, with an e-folding time of < 0.8 Gyr. These quenching timescales are shorter for more massive satellites but do not depend on host halo mass: the observed increase in satellite quiescent fraction with halo mass arises simply because of satellites quenching in a lower mass group prior to infall (group preprocessing), which is responsible for up to half of quenched satellites in massive clusters. Because of the long time delay before quenching starts, satellites experience significant stellar mass growth after infall, nearly identical to central galaxies. This fact provides key physical insight into the subhalo abundance matching method.

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.

Recent Hubble Space Telescope (HST) observations have revealed that a majority of active galactic nuclei (AGN) at z ~ 1-3 are resident in isolated disk galaxies, contrary to the usual expectation that AGN are triggered by mergers. Here we develop a new test of the cosmic evolution of supermassive black holes (SMBHs) in disk galaxies by considering the local population of SMBHs. We show that substantial SMBH growth in spiral galaxies is required as disks assemble. SMBHs exhibit a tight relation between their mass and the velocity dispersion of the spheroid within which they reside, the M_BH-sigma relation. In disk galaxies the bulge is the spheroid of interest. We explore the evolution of the M_BH-sigma relation when bulges form together with SMBHs on the M_BH-sigma relation and then slowly reform a disk around them. The formation of the disk compresses the bulge raising its sigma. We present evidence for such compression in the form of larger velocity dispersion of classical bulges compared with elliptical galaxies at the same mass. This compression leads to an offset in the M_BH-sigma relation if it is not accompanied by an increased M_BH. We quantify the expected offset based on photometric data and show that, on average, SMBHs must grow by ~ 50-65% just to remain on the M_BH-sigma relation. We find no significant offset in the M_BH-sigma relations of classical bulges and of ellipticals, implying that SMBHs have been growing along with disks. Our simulations demonstrate that SMBH growth is necessary for the local population of disk galaxies to have remained on the M_BH-sigma relation.

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.

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 exploit the SDSS galaxy groups catalogue of Yang et al. to study how the gas-phase metallicities of star-forming galaxies depend on environment. We find that satellite and central galaxies follow a qualitatively similar stellar mass (M_*) - gas-phase metallicity relation. Satellites, though, have higher gas-phase metallicities than equally massive centrals, and this difference increases with decreasing M_*. We also find that the gas-phase metallicity of satellites increases with halo mass at fixed stellar mass. This increment is more pronounced for less massive galaxies. We also show that low mass satellite galaxies have higher gas-phase metallicities than central galaxies of the same stellar metallicity. This difference becomes negligible for more massive galaxies of roughly solar metallicity. We demonstrate that the observed differences in gas-phase metallicity between centrals and satellites at fixed M_* are not a consequence of stellar mass stripping (advocated by Pasquali et al. in order to explain similar differences but in stellar metallicity), nor to the past star formation history of these galaxies as quantified by their surface mass density or gas mass fraction. Rather, we argue that these trends probably originate from a combination of three environmental effects: (i) strangulation, which prevents satellite galaxies from accreting new, low metallicity gas which would otherwise dilute their ISM, (ii) ram-pressure stripping of the outer gas disk, thereby inhibiting radial inflows of low-metallicity gas, and (iii) external pressure provided by the hot gas of the host halo which prevents metal-enriched outflows from escaping the galaxies. [abridged]

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 formulate the concept of non-linear and stochastic galaxy biasing in the framework of halo occupation statistics. Using two-point statistics in projection, we define the galaxy bias function, b_g(r_p), and the galaxy-dark matter cross-correlation function, R_{gm}(r_p), where r_p is the projected distance. We use the analytical halo model to predict how the scale dependence of b_g and R_{gm}, over the range 0.1 Mpc/h < r_p < 30 Mpc/h, depends on the non-linearity and stochasticity in halo occupation models. In particular we quantify the effect due to the presence of central galaxies, the assumption for the radial distribution of satellite galaxies, the richness of the halo, and the Poisson character of the probability to have a certain number of satellite galaxies in a halo of a certain mass. Overall, brighter galaxies reveal a stronger scale dependence, and out to a larger radius. In real-space, we find that galaxy bias becomes scale independent, with R_{gm} = 1, for radii r > 1 - 5 Mpc/h, depending on luminosity. However, galaxy bias is scale-dependent out to much larger radii when one uses the projected quantities defined in this paper. These projected bias functions have the advantage that they are more easily accessible observationally and that their scale dependence carries a wealth of information regarding the properties of galaxy biasing. To observationally constrain the parameters of the halo occupation statistics and to unveil the origin of galaxy biasing we propose the use of the bias function Gamma_{gm}(r_p)=b_g(r_p)/R_{gm}(r_p). This function is obtained via a combination of weak gravitational lensing and galaxy clustering, and it can be measured using existing and forthcoming imaging and spectroscopic galaxy surveys.

We present a statistical study of the luminosity functions of galaxies surrounding luminous red galaxies (LRGs) at average redshifts <z>=0.34 and <z>=0.65. The luminosity functions are derived by extracting source photometry around more than 40,000 LRGs and subtracting foreground and background contamination using randomly selected control fields. We show that at both studied redshifts the average luminosity functions of the LRGs and their satellite galaxies are poorly fitted by a Schechter function due to a luminosity gap between the centrals and their most luminous satellites. We utilize a two-component fit of a Schechter function plus a log-normal distribution to demonstrate that LRGs are typically brighter than their most luminous satellite by roughly 1.3 magnitudes. This luminosity gap implies that interactions within LRG environments are typically restricted to minor mergers with mass ratios of 1:4 or lower. The luminosity functions further imply that roughly 35% of the mass in the environment is locked in the LRG itself, supporting the idea that mass growth through major mergers within the environment is unlikely. Lastly, we show that the luminosity gap may be at least partially explained by the selection of LRGs as the gap can be reproduced by sparsely sampling a Schechter function. In that case LRGs may represent only a small fraction of central galaxies in similar mass halos.

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.

We combine constraints on the galaxy-dark matter connection with structural and dynamical scaling relations to investigate the angular momentum content of disc galaxies. For haloes with masses in the interval 10^{11.3} < M_vir/M_sun < 10^{12.7} we find that the galaxy spin parameters are independent of halo mass with <\lambda'_gal> = (J_gal/M_gal) / (\sqrt{2} R_vir V_vir) = 0.019^{+0.004}_{-0.003} (1sigma). This is significantly lower than for relaxed LCDM haloes, which have an average spin parameter <\lambda'_halo> = 0.031. The average ratio between the specific angular momentum of disk galaxies and their host dark matter haloes is therefore R_j = \lambda'_gal/\lambda'_halo = 0.61^{+0.13}_{-0.10}. This calls into question a standard assumption made in the majority of all (semi-analytical) models for (disc) galaxy formation, namely that R_j=1. Using simple disc formation models we show that it is particularly challenging to understand why R_j is independent of halo mass, while the galaxy formation efficiency (\epsilon_GF, proportional to the ratio of galaxy mass to halo mass) reveals a strong halo mass dependence. We argue that the empirical scaling relations between \epsilon_GF, R_j and halo mass require both feedback (i.e., galactic outflows) and angular momentum transfer from the baryons to the dark matter (i.e., dynamical friction). The efficiency of angular momentum loss need to decrease with increasing halo mass. Such a mass dependence may reflect a bias against forming stable discs in high mass, low spin haloes or a transition from cold-mode accretion in low mass haloes to hot-mode accretion at the massive end. However, current hydrodynamical simulations of galaxy formation, which should include these processes, seem unable to reproduce the empirical relation between \epsilon_GF and R_j. We conclude that the angular momentum build-up of galactic discs remains poorly understood.

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

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 investigate the origin of the relations between stellar mass and optical circular velocity for early-type (ETG) and late-type (LTG) galaxies --- the Faber-Jackson (FJ) and Tully-Fisher (TF) relations. We combine measurements of dark halo masses (from satellite kinematics and weak lensing), and the distribution of baryons in galaxies (from a new compilation of galaxy scaling relations), with constraints on dark halo structure from cosmological simulations. The principle unknowns are the halo response to galaxy formation and the stellar initial mass function (IMF). The slopes of the TF and FJ relations are naturally reproduced for a wide range of halo response and IMFs. However, models with a universal IMF and universal halo response cannot simultaneously reproduce the zero points of both the TF and FJ relations. For a model with a universal Chabrier IMF, LTGs require halo expansion, while ETGs require halo contraction. A Salpeter IMF is permitted for high mass (sigma > 180 km/s) ETGs, but is inconsistent for intermediate masses, unless V_circ(R_e)/sigma_e > 1.6. If the IMF is universal and close to Chabrier, we speculate that the presence of a major merger may be responsible for the contraction in ETGs while clumpy accreting streams and/or feedback leads to expansion in LTGs. Alternatively, a recently proposed variation in the IMF disfavors halo contraction in both types of galaxies. Finally we show that our models naturally reproduce flat and featureless circular velocity profiles within the optical regions of galaxies without fine-tuning.

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 review recent results on the dependence of various galaxy properties on environment at low redshift. As environmental indicators, we use group mass, group-centric radius, and the distinction between centrals and satellites; examined galaxy properties include star formation rate, colour, AGN fraction, age, metallicity and concentration. In general, satellite galaxies diverge more markedly from their central counterparts if they reside in more massive haloes. We show that these results are consistent with starvation being the main environmental effect, if one takes into account that satellites that reside in more massive haloes and at smaller halo-centric radii on average have been accreted a longer time ago. Nevertheless, environmental effects are not fully understood yet. In particular, it is puzzling that the impact of environment on a galaxy seems independent of its stellar mass. This may indicate that the stripping of the extended gas reservoir of satellite galaxies predominantly occurs via tidal forces rather than ram-pressure.

AGN feedback is believed to play an important role in shaping a variety of observed galaxy properties, as well as the evolution of their stellar masses and star formation rates. In particular, in the current theoretical paradigm of galaxy formation, AGN feedback is believed to play a crucial role in regulating the levels of activity in galaxies, in relatively massive halos at low redshift. Only in recent years, however, has detailed statistical information on the dependence of galaxy activity on stellar mass, parent halo mass and hierarchy has become available. In this paper, we compare the fractions of galaxies belonging to different activity classes (star-forming, AGN and radio active) with predictions from four different and independently developed semi-analytical models. We adopt empirical relations to convert physical properties into observables (H_alpha emission lines, OIII line strength and radio power). We demonstrate that all models used in this study reproduce the overall distributions of galaxies belonging to different activity classes as a function of stellar mass and halo mass: star forming galaxies and the strongest radio sources are preferentially associated with low-mass and high-mass galaxies/halos respectively. However, model predictions differ from observational measurements in a number of ways. All models used in our study predict that almost every >1.e12 Msun dark matter halo and/or >1.e11 Msun galaxy should host a bright radio source, while only a small fraction of galaxies belong to this class in the data. In addition, radio brightness is expected to depend strongly on the mass of the parent halo mass in the models, while strong and weak radio galaxies are found in similar environments in data. Our results highlight that the distribution of AGN as a function of stellar mass provides one of the most promising discriminants between different gas accretion schemes.

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

We study the evolution of the scaling relations between maximum circular velocity, stellar mass and optical half-light radius of star-forming disk-dominated galaxies in the context of LCDM-based galaxy formation models. Using data from the literature combined with new data from the DEEP2 and AEGIS surveys we show that there is a consistent observational and theoretical picture for the evolution of these scaling relations from z\sim 2 to z=0. The evolution of the observed stellar scaling relations is weaker than that of the virial scaling relations of dark matter haloes, which can be reproduced, both qualitatively and quantitatively, with a simple, cosmologically-motivated model for disk evolution inside growing NFW dark matter haloes. In this model optical half-light radii are smaller, both at fixed stellar mass and maximum circular velocity, at higher redshifts. This model also predicts that the scaling relations between baryonic quantities evolve even more weakly than the corresponding stellar relations. We emphasize, though, that this weak evolution does not imply that individual galaxies evolve weakly. On the contrary, individual galaxies grow strongly in mass, size and velocity, but in such a way that they move largely along the scaling relations. Finally, recent observations have claimed surprisingly large sizes for a number of star-forming disk galaxies at z \sim 2, which has caused some authors to suggest that high redshift disk galaxies have abnormally high spin parameters. However, we argue that the disk scale lengths in question have been systematically overestimated by a factor \sim 2, and that there is an offset of a factor \sim 1.4 between H\alpha sizes and optical sizes. Taking these effects into account, there is no indication that star forming galaxies at high redshifts (z\sim 2) have abnormally high spin parameters.

Using an analytical model, we study the evolution of subhalo, including its mass, angular momentum and merging time-scale. This model considers the dominant processes governing subhalo evolution, such as dynamical friction, tidal stripping and tidal heating. We find that in order to best match the evolution of angular momentum measured from N-body simulation, mass stripping by tidal force should become inefficient after subhalo has experienced a few passages of pericenter. It is also found that the often used Coulomb logarithm $\ln M/m$ has to be revised to best fit the merging time-scales from simulation. Combining the analytical model with the Extended Press-Schechter (EPS) based merger trees, we study the subhalo mass function, and their spatial distribution in a Milky-Way (MW) type halo. By tuning the tidal stripping efficiency, we can gain a better match to the subhalo mass function from simulation. The predicted distribution of subhaloes is found to agree with the distribution of MW satellites, but is more concentrated than the simulation results. The radial distribution of subhaloes depends weakly on subhaloes mass at both present day and the time of accretion, but strongly on the accretion time. Using the improved model, we measure the second moment of the subhalo occupation distribution, and it agrees well with the results of Kravtsov et al. (2004a) and Zheng et al. (2005).

We present a new systematic way of setting up galactic gas disks based on the assumption of detailed hydrodynamic equilibrium. To do this, we need to specify the density distribution and the velocity field which supports the disk. We first show that the required circular velocity has no dependence on the height above or below the midplane so long as the gas pressure is a function of density only. The assumption of disks being very thin enables us to decouple the vertical structure from the radial direction. Based on that, the equation of hydrostatic equilibrium together with the reduced Poisson equation leads to two sets of second-order non-linear differential equation, which are easily integrated to set-up a stable disk. We call one approach `density method' and the other one `potential method'. Gas disks in detailed balance are especially suitable for investigating the onset of the gravitational instability. We revisit the question of global, axisymmetric instability using fully three-dimensional disk simulations. The impact of disk thickness on the disk instability and the formation of spontaneously induced spirals is studied systematically with or without the presence of the stellar potential. In our models, the numerical results show that the threshold value for disk instability is shifted from unity to 0.69 for self-gravitating thick disks and to 0.75 for combined stellar and gas thick disks. The simulations also show that self-induced spirals occur in the correct regions and with the right numbers as predicted by the analytic theory.

Using estimates of dark halo masses from satellite kinematics, weak gravitational lensing, and halo abundance matching, combined with the Tully-Fisher and Faber-Jackson relations, we derive the mean relation between the optical, V_opt, and virial, V_200, circular velocities of early- and late-type galaxies at redshift z~0. For late-type galaxies V_opt ~ V_200 over the velocity range V_opt=90-260 km/s, and is consistent with V_opt = V_maxh (the maximum circular velocity of NFW dark matter haloes in the concordance LCDM cosmology). However, for early-type galaxies V_opt \ne V_200, with the exception of early-type galaxies with V_opt simeq 350 km/s. This is inconsistent with early-type galaxies being, in general, globally isothermal. For low mass (V_opt < 250 km/s) early-types V_opt > V_maxh, indicating that baryons have modified the potential well, while high mass (V_opt > 400 km/s) early-types have V_opt < V_maxh. Folding in measurements of the black hole mass - velocity dispersion relation, our results imply that the supermassive black hole - halo mass relation has a logarithmic slope which varies from ~1.4 at halo masses of ~10^{12} Msun/h to ~0.65 at halo masses of 10^{13.5} Msun/h. The values of V_opt/V_200 we infer for the Milky Way and M31 are lower than the values currently favored by direct observations and dynamical models. This offset is due to the fact that the Milky Way and M31 have higher V_opt and lower V_200 compared to typical late-type galaxies of the same stellar masses. We show that current high resolution cosmological hydrodynamical simulations are unable to form galaxies which simultaneously reproduce both the V_opt/V_200 ratio and the V_opt-M_star (Tully-Fisher/Faber-Jackson) relation.

We use a semi-analytic model for disk galaxies to explore the origin of the time evolution and small scatter of the galaxy SFR sequence -- the tight correlation between star-formation rate (SFR) and stellar mass (M_star). The steep decline of SFR from z~2 to the present, at fixed M_star, is a consequence of the following: First, disk galaxies are in a steady state with the SFR following the net (i.e., inflow minus outflow) gas accretion rate. The evolution of the SFR sequence is determined by evolution in the cosmological specific accretion rates, \propto (1+z)^{2.25}, but is found to be independent of feedback. Although feedback determines the outflow rates, it shifts galaxies along the SFR sequence, leaving its zero point invariant. Second, the conversion of accretion rate to SFR is materialized through gas density, not gas mass. Although the model SFR is an increasing function of both gas mass fraction and gas density, only the gas densities are predicted to evolve significantly with redshift. Third, star formation is fueled by molecular gas. Since the molecular gas fraction increases monotonically with increasing gas density, the model predicts strong evolution in the molecular gas fractions, increasing by an order of magnitude from z=0 to z~2. On the other hand, the model predicts that the effective surface density of atomic gas is ~10 M_sun pc^{-2}, independent of redshift, stellar mass or feedback. Our model suggests that the scatter in the SFR sequence reflects variations in the gas accretion history, and thus is insensitive to stellar mass, redshift or feedback. The large scatter in halo spin contributes negligibly, because it scatters galaxies along the SFR sequence. An observational consequence of this is that the scatter in the SFR sequence is independent of the size (both stellar and gaseous) of galaxy disks.

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.

We use a disk formation model to study the effects of galactic outflows (a.k.a. feedback) on the rotation velocity - stellar mass - disk size, gas fraction - stellar mass, and gas phase metalicity - stellar mass scaling relations of disk galaxies. We show that models without outflows are unable to explain these scaling relations, having both the wrong slopes and normalization. The problem can be traced to the model galaxies having too many baryons. Models with outflows can solve this "over-cooling" problem by removing gas before it has time to turn into stars. Models with both momentum and energy driven winds can reproduce the observed scaling relations. However, these models predict different slopes which, with better observations, may be used to discriminate between these models.

We present a new algorithm for identifying the substructure within simulated dark matter haloes. The method is an extension of that proposed by Tormen et al. (2004) and Giocoli et al. (2008a), which identifies a subhalo as a group of self-bound particles that prior to being accreted by the main progenitor of the host halo belonged to one and the same progenitor halo (hereafter satellite). However, this definition does not account for the fact that these satellite haloes themselves may also have substructure, which thus gives rise to sub-subhaloes, etc. Our new algorithm identifies substructures at all levels of this hierarchy, and we use it to determine the mass function of all substructure (counting sub-haloes, sub-subhaloes, etc.). On average, haloes which formed more recently tend to have a larger mass fraction in substructure and to be less concentrated than average haloes of the same mass. We provide quantitative fits to these correlations. Even though our algorithm is very different from that of Gao et al. (2004), we too find that the subhalo mass function per unit mass at redshift z = 0 is universal. This universality extends to any redshift only if one accounts for the fact that host haloes of a given mass are less concentrated at higher redshifts, and concentration and substructure abundance are anti-correlated. This universality allows a simple parametrization of the subhalo mass function integrated over all host halo masses, at any given time. We provide analytic fits to this function which should be useful in halo model analyses which equate galaxies with halo substructure when interpreting clustering in large sky surveys. Finally, we discuss systematic differences in the subhalo mass function that arise from different definitions of (host) halo mass.

We use a statistical approach to determine the relationship between the stellar masses of galaxies and the masses of the dark matter halos in which they reside. We obtain a parameterized stellar-to-halo mass (SHM) relation by populating halos and subhalos in an N-body simulation with galaxies and requiring that the observed stellar mass function be reproduced. We find good agreement with constraints from galaxy-galaxy lensing and predictions of semi-analytic models. Using this mapping, and the positions of the halos and subhalos obtained from the simulation, we find that our model predictions for the galaxy two-point correlation function (CF) as a function of stellar mass are in excellent agreement with the observed clustering properties in the SDSS at z=0. We show that the clustering data do not provide additional strong constraints on the SHM function and conclude that our model can therefore predict clustering as a function of stellar mass. We compute the conditional mass function, which yields the average number of galaxies with stellar masses in the range [m, m+dm] that reside in a halo of mass M. We study the redshift dependence of the SHM relation and show that, for low mass halos, the SHM ratio is lower at higher redshift. The derived SHM relation is used to predict the stellar mass dependent galaxy CF and bias at high redshift. Our model predicts that not only are massive galaxies more biased than low mass ones at all redshifts, but the bias increases more rapidly with increasing redshift for massive galaxies than for low mass ones. We present convenient fitting functions for the SHM relation as a function of redshift, the conditional mass function, and the bias as a function of stellar mass and redshift.

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]

We use a disk galaxy evolution model to investigate the impact of mass outflows (a.k.a. feedback) on disk galaxy scaling relations. Our model follows the accretion, cooling, star formation and ejection of baryonic mass inside growing dark matter haloes, with cosmologically motivated specific angular momentum distributions. Models without feedback produce disks that are too small and rotate too fast. Feedback reduces the baryonic masses of galaxies, resulting in larger disks with lower rotation velocities. Models with feedback can reproduce the zero points of the scaling relations between rotation velocity, stellar mass and disk size, but only in the absence of adiabatic contraction. Our feedback mechanism is maximally efficient in expelling mass, but our successful models require 25% of the SN energy, or 100% of the SN momentum, to drive the outflows. It remains to be seen whether such high efficiencies are realistic or not. Our energy and momentum driven wind models result in different slopes of various scaling relations, such as size - stellar mass, stellar mass - halo mass, and metallicity - stellar mass. Observations favor the energy driven wind at stellar masses below Mstar = 10^{10.5} Msun, but the momentum driven wind model at high masses. The ratio between the specific angular momentum of the baryons to that of the halo, (j_gal/m_gal), is not unity in our models. Yet this is the standard assumption in models of disk galaxy formation. Feedback preferentially ejects low angular momentum material because star formation is more efficient at smaller galactic radii. This results in (j_gal/m_gal) increasing with decreasing halo mass. This effect helps to resolve the discrepancy between the high spin parameters observed for dwarf galaxies with the low spin parameters predicted from LCDM. [Abridged]

We present a simple, empirically motivated model that simultaneously predicts the evolution of the mean size and the comoving mass density of massive early-type galaxies from z=2 to the present. First we demonstrate that some size evolution of the population can be expected simply due to the continuous emergence of early-type galaxies. SDSS data reveal that in the present-day universe more compact early-type galaxies with a given dynamical mass have older stellar populations. In contrast, at a given stellar velocity dispersion, SDSS data show that there is no relation between size and age, which implies that the velocity dispersion can be used to estimate the epoch at which galaxies stopped forming stars, turning into early-type galaxies. Applying such a 'formation' criterion to a large sample of nearby early-type galaxies, we predict the redshift evolution in the size distribution and the comoving mass density. The resulting evolution in the mean size is roughly half of the observed evolution. Then we include a prescription for the merger histories of galaxies between the 'formation' redshift and the present, based on cosmological simulations of the assembly of dark matter halos. Such mergers after the transformation into an early-type galaxy are presumably dissipationless ('dry'), where the increase in size is expected to be approximately proportional to the increase in mass. This model successfully reproduces the observed evolution since z~2 in the mean size and in the comoving mass density of massive early-type galaxies. We conclude that the recently measured, substantial size evolution of early-type galaxies can be explained by the combined effect of the continuous emergence of galaxies as early types and their subsequent growth through dry merging.

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 investigate the effects of changes in the cosmological parameters between the WMAP 1st, 3rd, and 5th year results on the structure of dark matter haloes. We use a set of simulations that cover 5 decades in halo mass ranging from the scales of dwarf galaxies (V_c ~30 km/s) to clusters of galaxies (V_c ~ 1000 km/s). We find that the concentration mass relation is a power law in all three cosmologies. However the slope is shallower and the zero point is lower moving from WMAP1 to WMAP5 to WMAP3. For haloes of mass log(M_200/Msun) = 10, 12, and 14 the differences in the concentration parameter between WMAP1 and WMAP3 are a factor of 1.55, 1.41, and 1.29, respectively. As we show, this brings the central densities of dark matter haloes in good agreement with the central densities of dwarf and low surface brightness galaxies inferred from their rotation curves, for both the WMAP3 and WMAP5 cosmologies. We also show that none of the existing toy models for the concentration-mass relation can reproduce our simulation results over the entire range of masses probed. In particular, the model of Bullock et al (B01) fails at the higher mass end (M > 1e13 Msun), while the NFW model of Navarro, Frenk & White (1997) fails dramatically at the low mass end (M < 1e12 Msun). We present a new model, based on a simple modification of that of B01, which reproduces the concentration-mass relations in our simulations over the entire range of masses probed (1e10 Msun < M < 1e15 Msun). Haloes in the WMAP3 cosmology (at a fixed mass) are more flatted compared to the WMAP1 cosmology, with a medium to long axis ration reduced by ~10 %. Finally, we show that the distribution of halo spin parameters is the same for all three cosmologies.

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.

Satellite kinematics can be used to probe the masses of dark matter haloes of central galaxies. In order to measure the kinematics with sufficient signal-to-noise, one uses the satellite galaxies of a large number of central galaxies stacked according to similar properties (e.g., luminosity). However, in general the relation between the luminosity of a central galaxy and the mass of its host halo will have non-zero scatter. Consequently, this stacking results in combining the kinematics of satellite galaxies in haloes of different masses, which complicates the interpretation of the data. In this paper we present an analytical framework to model satellite kinematics, properly accounting for this scatter and for various selection effects. We show that in the presence of scatter in the halo mass-luminosity relation, the commonly used velocity dispersion of satellite galaxies can not be used to infer a unique halo mass-luminosity relation. In particular, we demonstrate that there is a degeneracy between the mean and the scatter of the halo mass-luminosity relation. We present a new technique that can break this degeneracy, and which involves measuring the velocity dispersions using two different weighting schemes: host-weighting (each central galaxy gets the same weight) and satellite-weighting (each central galaxy gets a weight proportional to its number of satellites). The ratio between the velocity dispersions obtained using these two weighting schemes is a strong function of the scatter in the halo mass-luminosity relation, and can thus be used to infer a unique relation between light and mass from the kinematics of satellite galaxies.

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.

We derive new constraints on the mass of the Milky Way's dark matter halo, based on a set of halo stars from SDSS as kinematic tracers. Our sample comprises 2401 rigorously selected Blue Horizontal-Branch (BHB) halo stars drawn from SDSS DR-6. To interpret these distributions, we compare them to matched mock observations drawn from two different cosmological galaxy formation simulations designed to resemble the Milky Way, which we presume to have an appropriate orbital distribution of halo stars. We then determine which value of $\rm V_{cir}(r)$ brings the observed distribution into agreement with the corresponding distributions from the simulations. This procedure results in an estimate of the Milky Way's circular velocity curve to $\sim 60$ kpc, which is found to be slightly falling from the adopted value of $\rm 220 km s^{-1}$ at the Sun's location, and implies M$(<60 \rm kpc) = 4.0\pm 0.7\times 10^{11}$M$_\odot$. The radial dependence of $\rm V_{cir}(r)$, derived in statistically independent bins, is found to be consistent with the expectations from an NFW dark matter halo with the established stellar mass components at its center. If we assume an NFW halo profile of characteristic concentration holds, we can use the observations to estimate the virial mass of the Milky Way's dark matter halo, M$_{\rm vir}=1.0^{+0.3}_{-0.2} \times 10^{12}$M$_\odot$, which is lower than many previous estimates. This estimate implies that nearly 40% of the baryons within the virial radius of the Milky Way's dark matter halo reside in the stellar components of our Galaxy. A value for M$_{\rm vir}$ of only $\sim 1\times10^{12}$M$_\odot$ also (re-)opens the question of whether all of the Milky Way's satellite galaxies are on bound orbits.

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.

Using data from the Sloan Digital Sky Survey (SDSS) it has recently been shown that the red fraction of satellite galaxies increases with stellar mass. Semi-analytical models, however, predict red satellite fractions that are independent of stellar mass, and much higher than observed. It has been argued that this discrepancy owes to the fact that the models assume that satellite galaxies are instantaneously stripped of their hot gas reservoirs at the moment they are accreted into a bigger halo. In this letter we show that the fraction of red satellites can be brought in better agreement with the data by simply decreasing this stripping efficiency. However, this also results in a red fraction of massive centrals that is much too low. This owes to the fact that the massive centrals now accrete satellite galaxies that are bluer and more gas-rich. However, if a significant fraction of low mass satellite galaxies is tidally disrupted before being accreted by their central host galaxy, as suggested by recent studies, the red fractions of both centrals and satellites can be reproduced reasonably well. A problem remains with the red fraction of centrals of intermediate mass, which is likely to reflect an oversimplified treatment of AGN feedback.

A long standing problem for hierarchical disk galaxy formation models has been the simultaneous matching of the zero point of the Tully-Fisher relation and the galaxy luminosity function (LF). We illustrate this problem for a typical disk galaxy and discuss three solutions: low stellar mass-to-light ratios, low initial dark halo concentrations, and no halo contraction. We speculate that halo contraction may be reversed through a combination of mass ejection through feedback and angular momentum exchange brought about by dynamical friction between baryons and dark matter during the disk formation process.

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)

Using a cosmological N-Body simulation and a sample of re-simulated cluster-like haloes, we study the mass loss rates of dark matter subhaloes, and interpret the mass function of subhaloes at redshift zero in terms of the evolution of the mass function of systems accreted by the main halo progenitor. When expressed in terms of the ratio between the mass of the subhalo at the time of accretion and the present day host mass the unevolved subhalo mass function is found to be universal. However, the subhalo mass function at redshift zero clearly depends on $M_0$, in that more massive host haloes host more subhaloes. To relate the unevolved and evolved subhalo mass functions, we measure the subhalo mass loss rate as a function of host mass and redshift. We find that the average, specific mass loss rate of dark matter subhaloes depends mainly on redshift. These results suggest a pleasingly simple picture for the evolution and mass dependence of the evolved subhalo mass function. Less massive host haloes accrete their subhaloes earlier, which are thus subjected to mass loss for a longer time. In addition, their subhaloes are typically accreted by denser hosts, which causes an additional boost of the mass loss rate. To test the self-consistency of this picture, we use a merger trees constructed using the extended Press-Schechter formalism, and evolve the subhalo populations using the average mass loss rates obtained from our simulations, finding the subhalo mass functions to be in good agreement with the simulations. [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)

We construct a large data set of global structural parameters for 1300 field and cluster spiral galaxies and explore the joint distribution of luminosity L, optical rotation velocity V, and disk size R at I- and 2MASS K-bands. The I- and K-band velocity-luminosity (VL) relations have log-slopes of 0.29 and 0.27, respectively with sigma_ln(VL)~0.13, and show a small dependence on color and morphological type in the sense that redder, early-type disk galaxies rotate faster than bluer, later-type disk galaxies for most luminosities. The VL relation at I- and K-bands is independent of surface brightness, size and light concentration. The log-slope of the I- and K-band RL relations is a strong function of morphology and varies from 0.25 to 0.5. The average dispersion sigma_ln(RL) decreases from 0.33 at I-band to 0.29 at K, likely due to the 2MASS selection bias against lower surface brightness galaxies. Measurement uncertainties are sigma_ln(V)~0.09, sigma_ln(L)~0.14 and somewhat larger and harder to estimate for ln(R). The color dependence of the VL relation is consistent with expectations from stellar population synthesis models. The VL and RL residuals are largely uncorrelated with each other; the RV-RL residuals show only a weak positive correlation. These correlations suggest that scatter in luminosity is not a significant source of the scatter in the VL and RL relations. The observed scaling relations can be understood in the context of a model of disk galaxies embedded in dark matter halos that invokes low mean spin parameters and dark halo expansion, as we describe in our companion paper (Dutton et al. 2007). We discuss in two appendices various pitfalls of standard analytical derivations of galaxy scaling relations, including the Tully-Fisher relation with different slopes. (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.

Using a semi-analytical model for galaxy formation, combined with a large N-body simulation, we investigate the origin of the dichotomy among early-type galaxies. We find that boxy galaxies originate from mergers with a progenitor mass ratio $n < 2$ and with a combined cold gas mass fraction $F_{\rm cold} < 0.1$. Our model accurately reproduces the observed fraction of boxy systems as a function of luminosity and halo mass, for both central galaxies and satellites. After correcting for the stellar mass dependence, the properties of the last major merger of early-type galaxies are independent of their halo mass. This provides theoretical support for the conjecture of Pasquali et al (2007) that the stellar mass of an early-type galaxy is the main parameter that governs its isophotal shape. We argue that the observed dichotomy of early-type galaxies has a natural explanation within hierarchical structure formation, and does not require AGN feedback. Rather, we argue that it owes to the fact that more massive systems (i) have more massive progenitors, (ii) assemble later, and (iii) have a larger fraction of early-type progenitors. Each of these three trends causes the cold gas mass fraction of the progenitors of more massive early-types to be lower, so that their last major merger was dryer. Finally, our model predicts that (i) less than 10 percent of all early-type galaxies form in major mergers that involve two early-type progenitors, (ii) more than 95 percent of all boxy early-type galaxies with $M_* < 2 \times 10^{10} h^{-1} \Msun$ are satellite galaxies, and (iii) about 70 percent of all low mass early-types do not form a supermassive black hole binary at their last major merger. The latter may help to explain why low mass early-types have central cusps, while their massive counterparts have cores.

We study the dependence of the isophotal shape of early-type galaxies on their absolute B-band magnitude, their dynamical mass, and their nuclear activity and environment, using an unprecedented large sample of 847 early-type galaxies identified in the SDSS by Hao et al (2006). We find that the fraction of disky galaxies smoothly decreases with increasing luminosity. The large sample allows us to describe these trends accurately with tight linear relations that are statistically robust against the uncertainty in the isophotal shape measurements. There is also a host of significant correlations between the disky fraction and indicators of nuclear activity (both in the optical and in the radio) and environment (soft X-rays, group mass, group hierarchy). Our analysis shows however that these correlations can be accurately matched by assuming that the disky fraction depends only on galaxy luminosity or mass. We therefore conclude that neither the level of activity, nor group mass or group hierarchy help in better predicting the isophotal shape of early-type galaxies.

We use a series of cosmological N-body simulations for a flat LCDM cosmology to investigate the properties of dark matter haloes in the mass range 3.0e9-3.0e13 Msun. These properties include the concentration parameter (c), the spin parameter (lambda) and the mean axis ratio (q). For the concentration-mass relation we find c~M^(-0.11) in agreement with the model proposed by Bullock et al. even if we find a lower normalization (15%). The results for lambda and q are in good agreement with previous studies, while c and lambda are anti-correlated. In an attempt to remove unrelaxed haloes, we use the offset parameter (xoff), defined as the distance between the most bound particle and the center of mass. Removing haloes with large xoff increases the c by ~10%, lowers the lambda by ~15%, and removes the most prolate haloes. In addition, it largely removes the anti-correlation between c and lambda though not entirely. We also investigate the effects of the large-scale environment. We find that more concentrated haloes live in denser environments. Note, however, that the trend is weak compared to the scatter. For the spin parameters we find no environment dependence, while there is a weak indication that the most spherical haloes reside in denser region. Finally, using a simple model for disk galaxy formation we show that haloes that host low surface brightness galaxies are expected to be hosted by a biased sub-set of haloes. Not only do these haloes have spin parameters that are larger than average, they also have c that are 15% percent lower than the average at a given halo mass. We discuss the implications of all these findings for the claimed disagreement between halo concentrations inferred from LSB rotation curves, and those expected for a LCDM cosmology. (abridged)

We use observed rotation velocity-luminosity (VL) and size-luminosity (RL) relations to single out a specific scenario for disk galaxy formation in the LCDM cosmology. Our model involves four independent log-normal random variables: dark-halo concentration c, disk spin lam_gal, disk mass fraction m_gal, and stellar mass-to-light ratio M/L_I. A simultaneous match of the VL and RL zero points with adiabatic contraction requires low-c halos, but this model has V_2.2~1.8 V_vir (where V_2.2 and V_vir are the circular velocity at 2.2 disk scale lengths and the virial radius, respectively) which will be unable to match the luminosity function (LF). Similarly models without adiabatic contraction but standard c also predict high values of V_2.2/V_vir. Models in which disk formation induces an expansion rather than the commonly assumed contraction of the dark-matter halos have V_2.2~1.2 V_vir which allows a simultaneous fit of the LF. This may result from non-spherical, clumpy gas accretion, where dynamical friction transfers energy from the gas to the dark matter. This model requires low lam_gal and m_gal values, contrary to naive expectations. However, the low lam_gal is consistent with the notion that disk galaxies predominantly survive in halos with a quiet merger history, while a low m_gal is also indicated by galaxy-galaxy lensing. The smaller than expected scatter in the RL relation, and the lack of correlation between the residuals of the VL and RL relations, respectively, imply that the scatter in lam_gal and in c need to be smaller than predicted for LCDM halos, again consistent with the idea that disk galaxies preferentially reside in halos with a quiet merger history.

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.

Stellar-population analyses of today's galaxies show "downsizing", where the stars in more massive galaxies tend to have formed earlier and over a shorter time span. We show that this phenomenon is not necessarily "anti-hierarchical" but rather has its natural roots in the bottom-up clustering process of dark-matter haloes. While the main progenitor does indeed show an opposite effect, the integrated mass in all the progenitors down to a given minimum mass shows a robust downsizing that is qualitatively similar to what has been observed. These results are derived analytically from the standard extended Press Schechter (EPS) theory, and are confirmed by merger trees based on EPS or drawn from N-body simulations. The downsizing is valid for any minimum mass, as long as it is the same for all haloes at any given time, but the effect is weaker for smaller minimum mass. If efficient star formation is triggered by atomic cooling, then a minimum halo mass arises naturally from the minimum virial temperature for cooling, T~10^4 K, though for such a small minimum mass the effect is weaker than observed. Baryonic feedback effects, which are expected to stretch the duration of star formation in small galaxies and shut it down in massive haloes at late epochs, are likely to play a subsequent role in shaping up the final downsizing behaviour. Other appearances of downsizing, such as the decline with time of the typical mass of star-forming galaxies, may not be attributed to the gravitational clustering process but rather arise from the gas processes.

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

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

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 present a simple, semi-analytical model to compute the mass functions of dark matter subhaloes. The masses of subhaloes at their time of accretion are obtained from a standard merger tree. During the subsequent evolution, the subhaloes experience mass loss due to the combined effect of dynamical friction, tidal stripping, and tidal heating. Rather than integrating these effects along individual subhalo orbits, we consider the average mass loss rate, where the average is taken over all possible orbital configurations. This allows us to write the average mass loss rate as a simple function that depends only on redshift and on the instantaneous mass ratio of subhalo and parent halo. After calibrating the model by matching the subhalo mass function (SHMF) of cluster-sized dark matter haloes obtained from numerical simulations, we investigate the predicted mass and redshift dependence of the SHMF.We find that, contrary to previous claims, the subhalo mass function is not universal. Instead, both the slope and the normalization depend on the ratio of the parent halo mass, M, and the characteristic non-linear mass M*. This simply reflects a halo formation time dependence; more massive parent haloes form later, thus allowing less time for mass loss to operate. We analyze the halo-to-halo scatter, and show that the subhalo mass fraction of individual haloes depends most strongly on their accretion history in the last Gyr. Finally we provide a simple fitting function for the average SHMF of a parent halo of any mass at any redshift and for any cosmology, and briefly discuss several implications of our findings.

We use a simple model of disc-galaxy formation in the LCDM cosmology to simultaneously reproduce the slopes, zero-points, scatter and uncorrelated residuals in the observed velocity-luminosity (VL) and radius-luminosity (RL) relations. Observed I-band luminosities are converted to stellar masses using the IMF-dependent relation between stellar mass-to-light ratio, M/L, and color. The model treats halo concentration, spin parameter, and disc mass fraction as independent log-normal random variables. Our main conclusion is that the VL and RL zero-points and the uncorrelated residuals can only be reproduced simultaneously if adiabatic contraction is avoided. One or the other could be fixed by appealing to unrealistic values for M/L, halo concentration c, and spin parameter lambda, but not all together. The small VL scatter is naturally determined by the predicted scatter in c and in M/L, quite independent of the large scatter in lambda. However, the RL scatter, driven by the scatter in lambda, can be as low as observed only if sigma_lambda simeq 0.25, about half the value predicted for CDM haloes. This may indicate that discs form in a special subset of haloes. The VL slope is reproduced once star formation occurs only above a threshold surface density, and the RL slope implies that the disc mass fraction is increasing with halo mass, consistent with feedback effects. A model that incorporates the above ingredients provides a simultaneous fit to all the observed features. In particular, the elimination of adiabatic contraction allows 80% disc contribution to the observed rotation velocity of bright discs at 2.2 disc scale lengths. The lack of halo contraction may indicate that disc formation is not as smooth as typically envisioned, but instead may involve clumpy, cold streams.

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

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) This paper presents an exploration of two fundamental scaling relations of spiral galaxies, the luminosity-rotation speed (or Tully-Fisher; TF) relation, and the size-luminosity (SL) relation, and the dependences of their scatter, at red and IR bands. We verify that the observed virial relations of disk galaxies are given by Vobs \prop L^0.31 and Rd \prop L^0.33 using distance-redshift surveys of high surface brightness (HSB) and low surface brightness (LSB) non-interacting galaxies. Thus the galaxy surface brightness (SB) is Sigma \prop Rd \prop Vobs \prop Mvir^0.33. We derive the first extensive J-band TF and SL relations using 2MASS data products. The TF residuals are fully independent of SB and other tested observables. TF relations (TFRs) for various morphological types have different zero points, but a common slope, such that early-type (redder) disks rotate faster than later-types (bluer) of the same luminosity. This dependence of the TFR is a direct consequence of the more fundamental dependence on color, which itself is related to the galaxy star formation history. The SL relation scatter is mostly dominated by SB, but color also plays a small role. The TF and SL residuals for HSB and LSB galaxies are weakly correlated, in agreement with the earlier claim by Courteau & Rix (1999). This suggests that spiral disks of all color, brightness, and barredness are, on average, dominated by dark matter even in their inner parts. The observed scaling relations of disk galaxies are likely the result of the simplest scenario of galaxy formation, but their color dependence remains a challenge to hierarchical structure formation models. The scaling relations and their derivatives are contrasted against galaxy formation models in our companion paper.

We examine the effect of preheating of the intergalactic medium on galaxy formation using cosmological hydrodynamical simulations. By performing simulations both with and without a simple model for preheating, we analyse and compare the angular momentum distributions of the dark matter and the baryons. Preheating unbinds baryons from their dark matter haloes, yielding a baryonic mass fraction that declines with decreasing halo mass. In addition, the spin parameter of the gas is reduced with respect to the case without preheating, while the misalignment between the angular momentum directions of the gas and dark matter increases strongly. The angular momentum distributions of individual haloes reveal that preheating decreases (increases) the mass fraction with low (negative) specific angular momentum. We discuss the implications of these findings for the formation of disk galaxies in a preheated intergalactic medium, and compare our results to the predictions of Maller & Dekel (2002), who propose an alternative interpretation for the origin of the angular momentum of (proto)-galaxies.

We present Halpha rotation curves for a sample of 15 dwarf and LSB galaxies. From these, we derive limits on the slopes of the central mass distributions. Assuming the density distributions of dark matter halos follow a power-law at small radii, rho(r)~r^(-alpha), we find inner slopes in the range 0<alpha<1 for most galaxies. In general, halos with constant density cores (\alpha=0) provide somewhat better fits, but the majority of our galaxies (~75%) are also consistent with alpha=1, provided that the R-band mass-to-light ratios are smaller than about 2. Halos with alpha=1.5, however, are ruled out in virtually every case. To investigate the robustness of these results we discuss and model several possible causes of systematic errors including non-circular motions, slit width, seeing, and slit alignment errors. Taking the associated uncertainties into account, we conclude that even for the 25% of the cases where alpha=1 seems inconsistent with the rotation curves, we cannot rule out cusp slopes this steep. Inclusion of literature samples similar to the one presented here leads to the same conclusion when possible systematic errors are taken into account. In the ongoing debate on whether the rotation curves of dwarf and LSB galaxies are consistent with predictions for a CDM universe, we argue that our sample and the literature samples discussed in this paper provide insufficient evidence to rule out halos with alpha=1. At the same time, we note that none of the galaxies in these samples require halos with steep cusps, as most are equally well or better explained by constant density cores. (abridged)

The standard picture for the formation of disk galaxies has come to face a number of problems related to angular momentum. Numerical simulations produce disks that are too small, and analytical models predict too much low angular momentum compared to observations. In addition, angular momentum distributions of proto-galaxies in numerical simulations reveal large mass fractions with negative specific angular momentum, which implies that the gas cannot conserve its detailed distribution of specific angular momentum when cooling to form the disk. This crumples one of the main pillars of our standard picture, and indicates that a new spin on the angular momentum acquisition of disk galaxies may be required.

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.

We use non-radiative N-body/SPH simulations of structure formation in a LCDM cosmology to compare the angular momentum distributions of dark matter and gas in a large sample of halos. We show that the two components have identical spin parameter distributions and that their angular momentum distributions within individual halos are very similar, all in excellent agreement with standard assumptions. Despite these similarities, however, we find that the angular momentum vectors of the gas and dark matter are poorly aligned, with a median misalignment angle of about 30 degrees, which might have important implications for spin correlation statistics used in weak lensing studies. We present distributions for the component of the angular momentum that is aligned with the total angular momentum of each halo, and find that for between 5 and 50 percent of the mass this component is negative. This disagrees with the `Universal' angular momentum distribution suggested by Bullock et al. (2001), for which the mass fraction with negative specific angular momentum is zero. We discuss the implications of our results for the formation of disk galaxies. Since galactic disks generally do not contain counter-rotating stars or gas, disk formation cannot occur under detailed conservation of specific angular momentum. We suggest that the material with negative specific angular momentum combines with positive angular momentum material to build a bulge component, and show that in such a scenario the remaining material can form a disk with a density distribution that is very close to exponential.

We investigate the hypothesis that the cores of elliptical galaxies and bulges are created from the binding energy liberated by the coalescence of supermassive binary black holes during galaxy mergers. Assuming that the central density profiles of galaxies were initially steep power laws, we define the ``mass deficit'' as the mass in stars that had to be removed from the nucleus in order to produce the observed core. We use nonparametric deprojection to compute the mass deficit in a sample of 35 early-type galaxies with high-resolution imaging data. We find that the mass deficit correlates well with the mass of the nuclear black hole, consistent with the predictions of merger models. We argue that core sizes in halos of non-interacting dark matter particles should be comparable to those observed in the stars.

The puzzling correlation between the spin parameter lambda of galactic disks and the disk-to-halo mass fraction fdisk is investigated. We show that such a correlation arises naturally from uncertainties in determining the virial masses of dark matter halos. This result leads to the conclusion that the halo properties derived from fits to observed rotation curves are still very uncertain which might explain part of the disagreements between cosmological models and observations. We analyse lambda and fdisk as function of the adopted halo virial mass. Reasonable halo concentrations require fdisk=0.01-0.07 which is significantly smaller than the universal baryon fraction. Most of the available gas either never settled into the galactic disks or was ejected subsequently. In both cases it is not very surprising that the specific angular momentum distribution of galactic disks does not agree with the cosmological predictions which neglect these effects.

I present detailed models for the formation of disk galaxies, and investigate which observables are best suited as virial mass estimators. Contrary to naive expectations, the luminosities and circular velocities of disk galaxies are extremely poor indicators of total virial mass. Instead, I show that the product of disk scale length and rotation velocity squared yields a much more robust estimate. Finally, I show how this estimator may be used to put limits on the efficiencies of cooling and feedback during the process of galaxy formation.

We present detailed, analytical models for the formation of disc galaxies to investigate the impact that cooling and feedback have on their structural properties. In particular, we investigate which observables extracted directly from the models are best suited as virial mass estimators, and to what extent they allow the recovery of the model input parameters regarding the feedback and cooling efficiencies. Contrary to naive expectations, the luminosities and circular velocities of disc galaxies are extremely poor indicators of total virial mass. Instead, we show that the product of disc scale length and rotation velocity squared yields a much more robust estimate. We show that feedback can cause a narrow correlation between galaxy mass fraction and halo spin parameter, similar to that found recently by van den Bosch, Burkert and Swaters from an analysis of dwarf galaxy rotation curves. Finally we investigate the impact that cooling and feedback have on the colors, metallicities, star formation histories and Tully-Fisher relation of disc galaxies.

We use the extended Press-Schechter formalism to investigate the rate at which cold dark matter haloes accrete mass. We discuss the shortcomings of previous methods that have been used to compute the mass accretion histories of dark matter haloes, and present an improved method based on the N-branch merger tree algorithm of Somerville & Kolatt. We show that this method no longer suffers from inconsistencies in halo formation times, and compare its predictions with high resolution N-body simulations. Although the overall agreement is reasonable, there are slight inconsistencies which are most easily interpreted as a reflection of ellipsoidal collapse (as opposed to spherical collapse assumed in the Press-Schechter formalism). We show that the average mass accretion histories follow a simple, universal profile, and we present a recipe for computing the two scale-parameters which is applicable to a wide range of halo masses and cosmologies. Together with the universal profiles for the density and angular momentum distributions of CDM haloes, these universal mass accretion histories provide a simple but accurate framework for modeling the structure and formation of dark matter haloes. In particular, they can be used as a backbone for modeling various aspects of galaxy formation where one is not interested in the detailed effects of merging. As an example we use the universal mass accretion history to compute the rate at which haloes accrete mass, which we compare to the cosmic star formation history of the Universe.

We present new models for the formation of disk galaxies that improve upon previous models by following the detailed accretion and cooling of the baryonic mass, and by using realistic distributions of specific angular momentum. Under the assumption of detailed angular momentum conservation the disks that form have density distributions that are more centrally concentrated than an exponential. We examine the influence of star formation, bulge formation, and feedback on the outcome of the surface brightness distributions of the stars. Whereas the models are succesful in reproducing high surface brightness disks, they fail to produce exponential disks of low surface brightness. In addition, the models predict truncation radii in the gas that occur at too high surface density. This signals a new problem for the standard model of disk formation: if the baryonic component of the protogalaxies out of which disk galaxies form have the same angular momentum distribution as the dark matter, disks are too compact.

We compute the specific angular momentum distributions of a sample of low mass disk galaxies observed by Swaters. We compare these distributions to those of dark matter haloes obtained by Bullock et al. from high resolution N-body simulations of structure formation in a LCDM Universe. We find that although the disk mass fractions are significantly smaller than the Universal baryon fraction, the total specific angular momenta of the disks are in good agreement with those of dark matter haloes. This suggests that disks form out of only a small fraction of the available baryons, but yet manage to draw most of the available angular momentum. In addition we find that the angular momentum distributions of disks are clearly distinct from those of the dark matter; disks lack predominantly both low and high specific angular momentum. Understanding these findings in terms of a coherent picture for disk formation is challenging. Cooling, feedback and stripping, which are the main mechanisms to explain the small disk mass fractions found, seem unable to simultaneously explain the disk's angular momentum distribution. In fact, it seems that the baryons that make up the disks must have been born out of angular momentum distributions that are clearly distinct from those of LCDM haloes. Therefore, understanding the angular momentum content of disk galaxies remains an important challenge for our picture of galaxy formation.

The standard cold dark matter (CDM) model has recently been challenged by the claim that dwarf galaxies have dark matter halos with constant density cores. Consequently, numerous alternative dark matter candidates have recently been proposed. In this paper, we scrutinize the observational evidence for the incongruity between dwarf galaxies and the CDM model. To this end, we analyze the rotation curves of 20 late-type dwarf galaxies studied by Swaters (1999). Taking the effects of beam-smearing and adiabatic contraction into account, we fit mass models to these rotation curves with dark matter halos with different cusp slopes, ranging from constant density cores to r^{-2} cusps. Uncertainties in the stellar mass-to-light ratio and the limited spatial sampling of the halo's density distribution hamper a unique mass decomposition. Consequently, the rotation curves in our sample cannot be used to discriminate between dark halos with constant density cores and r^{-1} cusps. We show that the dwarf galaxies analyzed here are consistent with cold dark matter halos in a LCDM cosmology, and that there is thus no need to abandon the idea that dark matter is cold and collisionless. However, the data is also consistent with any alternative dark matter model that produces dark matter halos with central cusps less steep than r^{-1.5}. In fact, we argue that based on existing rotation curves alone at best weak limits can be obtained on cosmological parameters and/or the nature of the dark matter.

We present high resolution R-band images of the central regions of 67 early-type galaxies obtained with the Wide Field and Planetary Camera 2 (WFPC2) aboard the Hubble Space Telescope (HST). Our sample strikingly confirms the complex morphologies of the central regions of early-type galaxies. In particular, we detect dust in 43 percent of all galaxies, and evidence for embedded stellar disks in a remarkably large fraction of 51 percent. In 14 of those galaxies the disk-like structures are misaligned with the main galaxy, suggesting that they correspond to stellar bars in S0 galaxies. We analyze the luminosity profiles of the galaxies in our sample, and classify galaxies according to their central cusp slope. To a large extent we confirm the clear dichotomy found in previous HST surveys: bright, boxy ellipticals with shallow inner cusps (`core' galaxies) on one hand and faint, disky ellipticals with steep central cusps (`power-law' galaxies) on the other hand. The advantages and shortcomings of classification schemes utilizing the extrapolated central cusp slope are discussed, and it is shown that this cusp slope might be an inadequate representation for galaxies whose luminosity profile slope changes smoothly with radius rather than resembling a broken power-law. In fact, we find evidence for an `intermediate' class of galaxies, that cannot unambiguously be classified as either core or power-law galaxies, and which have central cusp slopes and absolute magnitudes intermediate between those of core and power-law galaxies.

We examine the dust properties of a nearby distance-limited sample of early type galaxies using the WFPC2 of the Hubble Space Telescope. Dust is detected in 29 out of 67 galaxies (43%), including 12 with small nuclear dusty disks. In a separate sample of 40 galaxies biased for the detection of dust by virtue of their detection in the IRAS 100 micron band, dust is found in ~78% of the galaxies, 15 of which contain dusty disks. In those galaxies with detectable dust, the apparent mass of the dust correlates with radio and far infrared luminosity, becoming more significant for systems with filamentary dust. A majority of IRAS and radio detections are also associated with dusty galaxies rather than dustless galaxies. This indicates that thermal emission from clumpy, filamentary dust is the main source of the far-IR radiation in early type galaxies. Dust in small disk-like morphology tends to be well aligned with the major axis of the host galaxies, while filamentary dust appears to be more randomly distributed with no preference for alignment with any major galactic structure. This suggests that, if the dusty disks and filaments have a common origin, the dust originates externally and requires time to dynamically relax and settle in the galaxy potential in the form of compact disks. More galaxies with visible dust than without dust display emission lines, indicative of ionized gas, although such nuclear activity does not show a preference for dusty disk over filamentary dust. There appears to be a weak relationship between the mass of the dusty disks and central velocity dispersion of the galaxy, suggesting a connection with a similar recently recognized relationship between the latter and the black hole mass.

We discuss a comparison of models for the formation of disk galaxies both in a Universe dominated by cold dark matter (CDM) and one in which the force law is given by modified Newtonian dynamics (MOND). Our main aim is to address the claim made by McGaugh & de Blok that CDM suffers from severe fine-tuning problems, which are circumvented under MOND. As we show, CDM indeed requires some amount of tuning of the feedback efficiencies to obtain a Tully-Fisher relation (TFR) as steep as observed. However, that same model is in excellent agreement with a wide variety of additional observations. Therefore, the modest amount of feedback needed should not be regarded a fine-tuning problem. Instead, its requirement should be considered a generic prediction for CDM, which might be tested with future observations and with detailed modeling of feedback processes. We also show that galaxy formation in a MOND universe can not simultaneously reproduce the TFR and the lack of high surface brightness dwarf galaxies. We thus conclude that CDM is a more viable theory for the formation of disk galaxies than MOND.

We present detailed semi-analytical models for the formation of disk galaxies both in a Universe dominated by dark matter (DM), and in one for which the force law is given by modified Newtonian dynamics (MOND). We tune the models to fit the observed near-infrared Tully-Fisher (TF) relation, and compare numerous predictions of the resulting models with observations. The DM and MOND models are almost indistinguishable. They both yield gas mass fractions and dynamical mass-to-light ratios which are in good agreement with observations. Both models reproduce the narrow relation between global mass-to-light ratio and central surface brightness, and reveal a characteristic acceleration, contrary to claims that these relations are not predicted by DM models. Both models require SN feedback in order to reproduce the lack of high surface brightness dwarf galaxies. However, the introduction of feedback to the MOND models steepens the TF relation and increases the scatter, making MOND only marginally consistent with observations. The most serious problem for the DM models is their prediction of steep central rotation curves. However, the DM rotation curves are only slightly steeper than those of MOND, and are only marginally inconsistent with the poor resolution data on LSB galaxies.

We re-examine the disk-halo decompositions of the rotation curves of low surface brightness (LSB) galaxies with V_max > 80 km/s, taking full account of the effects of beam smearing. We show that the rotation curves of these LSB galaxies are, contrary to previous claims, consistent with dark halos with steep central cusps in their density distribution. In fact, the spatial resolution of the data is not sufficient to put any meaningful constraints on the density profiles of their dark halos, or on cosmological parameters. This has important implications for numerous claims made in the literature regarding the halos of LSB galaxies, such as the self-similarity of their rotation curves, and their inconsistency with certain cosmological models or with cold dark matter altogether. Only in one case is the data of sufficient spatial resolution to obtain reliable constraints on the slope of the central density distribution of the dark matter halo. For this single case, we find a central density cusp with a slope 0.55 < alpha < 1.26 at the 99.73 percent confidence levels. This contrasts strongly with the results for two dwarf galaxies (V_max < 70 km/s) that we analyze, which yield alpha < 0.5 at the same level of confidence. This possibly suggests that halos with constant density cores are restricted to low-mass systems. We show that violent outflows of baryonic matter by supernovae feedback can reproduce this mass-dependence of halo cusp slopes.

New models for the formation of disk galaxies are presented. I discuss the contraints on galaxy formation that follow from fitting the model to the near-infrared Tully-Fisher (TF) relation, with an emphasis on reproducing the small amount of scatter observed. Once the parameters that describe the supernova feedback are tuned to fit the slope of the observed TF relation, the model reproduces the correct amount of TF scatter, and yields gas mass fractions, mass-to-light ratios, and characteristic accelerations that are all in excellent agreement with observations.

We present new semi-analytical models for the formation of disk galaxies with the purpose of investigating the origin of the near-infrared Tully-Fisher (TF) relation. The models assume that disks are formed by cooling of the baryons inside dark halos with realistic density profiles, and that the baryons conserve their specific angular momentum. Only gas with densities above the critical density given by Toomre's stability criterion is considered eligible for star formation, and a simple recipe for supernovae feedback is included. We emphasize the importance of extracting the proper luminosity and velocity measures from the models, something that has often been ignored in the past. The observed K-band TF relation has a slope that is steeper than simple predictions based on dynamical arguments suggest. Taking the stability related star formation threshold densities into account steepens the TF relation, decreases its scatter, and yields gas mass fractions that are in excellent agreement with observations. In order for the TF slope to be as steep as observed, further physics are required. We argue that the characteristics of the observed near-infrared TF relation do not reflect systematic variations in stellar populations, or cosmological initial conditions, but are governed by feedback. Finally we show that our models provide a natural explanation for the small amount of scatter that makes the TF relation useful as a cosmological distance indicator.

A simple, inside-out formation scenario for bulges is presented. Stability arguments are used to suggest that bulges form out of the low-angular momentum baryons that cool inside a virialized dark halo: a self-regulating mechanism transforms part of the baryonic matter into a bulge component, such that the remainder of the baryons can settle in a stable disk component. We couple this scenario to the Fall & Efstathiou theory of disk formation and search for the parameters and physical processes that determine the disk-to-bulge ratio, and therefore explain to a large extent the origin of the Hubble sequence. This simple galaxy formation scenario can account for both spirals and S0s, but fails to incorporate more bulge dominated systems.

We examine the distributions of eccentricities of orbits within mass distributions like those we see for galaxies and clusters. A comprehensive understanding of these orbital properties is essential to calculate the rates of physical processes relevant to the formation and evolution of galaxies and clusters. We derive the orbital eccentricity distributions for a number of spherical potentials. These distributions depend strongly on the velocity anisotropy, but only slightly on the shape of the potential. We also present high resolution N-body simulations of the orbital decay of satellite systems on eccentric orbits in an isothermal halo. The dynamical friction timescales are found to decrease with increasing orbital eccentricity, but the dependence is weaker than previously suggested. The orbital eccentricity stays remarkably constant throughout the decay; although the eccentricity decreases near pericenter, it increases again near apocenter, such that there is no net circularization. Finally, we discuss a number of applications. We find that the distribution of orbital eccentricities of the globular clusters in the Milky Way is consistent with a close to isotropic velocity distribution. The limited data on the galactic system of satellites, however, appears to be different and warrants further investigation as a clue to the formation and evolution of our Milky Way and its halo substructure.

We have constructed axisymmetric dynamical models of the edge-on S0 galaxy NGC 4342: simple two-integral Jeans models as well as fully general, three-integral models using a modified version of Schwarzschild's orbit superposition technique. The two-integral models suggest a black hole (BH) of 3 or $6\times 10^8 M_\odot$, depending on the data set. The three-integral models can fit all ground-based and HST data simultaneously, but only when a central BH is included. Models without BH are ruled out at better than 99.73% confidence level. We determine a BH mass of $3.0^{+1.7}_{-1.0} \times 10^8 M_\odot$. This corresponds to 2.6% of the bulge mass, making NGC 4342 one of the galaxies with the highest BH mass to bulge mass ratio currently known.

Hubble Space Telescope images of two early-type galaxies harboring both nuclear and outer stellar discs are studied in detail. By means of a photometric decomposition, the images of NGC 4342 and NGC 4570 are analyzed and the photometric properties of the nuclear discs investigated. We find a continuity of properties in the parameter space defined by the central surface brightness and the scalelength of discs in spirals, S0s and embedded discs in ellipticals, in the sense that the nuclear discs extend the observed disc properties even further towards smaller scalelengths and brighter central surface brightnesses. When including the nuclear discs, disc properties span more than four orders of magnitude in both scalelength and central surface brightness. The nuclear discs studied here are the smallest and brightest stellar discs known, and as such, they are as extreme in their photometric properties as Malin I, when compared to typical galactic discs that obey Freeman's law. We discuss a possible formation scenario in which the double-disc structure observed in these galaxies has been shaped by now dissolved bars. Based on the fact that the black holes known to exist in some of these galaxies have masses comparable to those of the nuclear discs, we explore a possible link between the black holes and the nuclear discs.

We present axisymmetric dynamical models of the edge-on S0 galaxy NGC 4342. A combination of observations from the ground and with the Hubble Space Telescope (HST) has shown that NGC 4342 rotates rapidly and has a strong central increase in velocity dispersion. We construct simple two-integral Jeans models as well as fully general, three-integral models. The latter are built using a modified version of Schwarzschild's orbit-superposition technique. The two-integral Jeans models suggest a black hole (BH) mass between 3 and 6x10^8 Msun, depending on the data set used to constrain the model, but they fail to fit the details of the observed kinematics. The three-integral models can fit all ground-based and HST data simultaneously, but only when a central BH is included. Models without BH are ruled out to a confidence level better than 99.73 per cent. We determine a BH mass of 3x10^8 Msun. This corresponds to 2.6 per cent of the total mass of the bulge, making NGC 4342 one of the galaxies with the highest BH mass to bulge mass ratio currently known. The models that best fit the data do not have a two-integral phase-space distribution function. They have rather complex dynamical structures: the velocity anisotropies are strong functions of radius reflecting the multi-component structure of this galaxy. The best fit model without BH tries to fit the high central velocity dispersion by placing stars on radial orbits. The measured rotation velocities, however, restrict the amount of radial anisotropy such that the central velocity dispersion measured with the HST can only be fit when a massive BH is included in the models.

We investigate the formation of disk-bulge-halo systems by including bulges in the Fall & Efstathiou theory of disk formation. This allows an investigation of bulge dominated disk galaxies, such as S0s and disky ellipticals. The aim of this paper is to examine whether spirals, S0s, and disky elllipticals all can be incorporated in one simple galaxy formation scenario. We investigate an inside-out formation scenario in which the inner, low angular momentum material is assumed to form the bulge. Stability arguments are used to suggest that this bulge formation is a self-regulating process in which the bulge grows until it is massive enough to allow the remaining gas to form a stable disk component. We search for the parameters and physical processes that determine the disk-to-bulge ratio, and therewith explain to a large extent the origin of the Hubble sequence. The spread in halo angular momenta coupled with a spread in the formation redshifts can explain the observed spread in disk properties and disk-to-bulge ratios from spirals to S0s. Disky ellipticals, however, are too compact to be incorporated. Upon introducing a galaxy formation efficiency that depends on redshift, as suggested by the small amount of scatter in the observed Tully-Fisher relation, this scenario naturally predicts a density-morphology relation, the amplitude of which depends on the baryon fraction of the Universe.

We present an HST study of the nuclear region of the E4 radio galaxy NGC 7052, which has a nuclear disk of dust and gas. The WFPC2 was used to obtain B, V and I broad-band images and an H_alpha+[NII] narrow-band image. The FOS was used to obtain H_alpha+[NII] spectra along the major axis, using a 0.26 arcsec diameter circular aperture. The observed rotation velocity of the ionized gas is V = 155 +/- 17 km/s at r = 0.2 arcsec from the nucleus. The Gaussian dispersion of the emission lines increases from sigma = 70 km/s at r=1 arcsec, to sigma = 400 km/s on the nucleus. To interpret the gas kinematics we construct axisymmetric models in which the gas and dust reside in a disk in the equatorial plane of the stellar body. It is assumed that the gas moves on circular orbits, with an intrinsic velocity dispersion due to turbulence. The circular velocity is calculated from the combined gravitational potential of the stars and a possible nuclear black hole (BH). Models without a BH predict a rotation curve that is shallower than observed (V_pred = 92 km/s at r = 0.2 arcsec), and are ruled out at > 99% confidence. Models with a BH of 3.3^{+2.3}_{-1.3} x 10^8 solar masses provide an acceptable fit. NGC 7052 can be added to the list of active galaxies for which HST spectra of a nuclear gas disk provide evidence for the presence of a central BH. The BH masses inferred for M87, M84, NGC 6251, NGC 4261 and NGC 7052 span a range of a factor 10, with NGC 7052 falling on the low end. By contrast, the luminosities of these galaxies are identical to within 25%. Any relation between BH mass and luminosity, as suggested by independent arguments, must therefore have a scatter of at least a factor 10.

We present circumstantial evidence that the central region of the edge-on S0 galaxy NGC 4570, which harbors a 150 pc scale nuclear disc in addition to its main outer disc, has been shaped under the influence of a small (approx. 500pc) bar. This is based on the finding of two edge-on rings, whose locations are consistent with the Inner Lindblad and Ultra Harmonic Resonances of a rapidly tumbling triaxial potential. Observed features in the photometry and rotation curve nicely correspond with the positions of the main resonances, strengthening the case for a tumbling bar potential. The relative blue colour of the ILR ring, and the complete absence of any detected ISM, indicates that the nuclear ring is made of relatively young (< 2 Gyr) stars. We discuss a possible secular-evolution scenario for this complex multi-component galaxy, which may also apply to many other S0-galaxies with observed rings and/or multiple disc-components.

We present multi-colour (U, V and I) photometry obtained with the second Wide Field and Planetary Camera (WFPC2) on the Hubble Space Telescope (HST), and spectra taken with the William Herschel Telescope (WHT) and the HST Faint Object Spectrograph (FOS), for the early-type galaxies NGC 4342 and NGC 4570. These galaxies are intermediate between ellipticals and lenticulars, and they both have a small nuclear stellar disc in addition to their main outer disc. Colour images reveal no colour differences between the nuclear discs and the bulges. Comparison of the U-V and V-I colours with stellar population models indicates that the central regions of both galaxies are of intermediate age (approx. 8 Gyr) and of high metallicity. The long-slit WHT spectra and the FOS spectra (of four times higher spatial resolution) reveal that both galaxies are rapidly rotating, and they both have velocity dispersions that increase strongly towards the centre. The FOS spectra of NGC 4342 indicate a central velocity dispersion of 420 km/s, higher than the 320 km/s measured from the WHT spectra. Also, the nuclear rotation gradient measured with the FOS is steeper than that measured with the WHT; it reaches V = 200 km/s at 0.25''. The rapid stellar motions seen in the centre of NGC 4342 suggest a large central mass concentration, possibly a massive black hole. The kinematics of the more massive NGC 4570 are less spectacular, with a central velocity dispersion of 250 km/s and a central rotation curve that reaches only 60 km/s at 0.25''.

We discuss photometric and spectroscopic data of the E/S0 galaxy NGC 4342, obtained with the Hubble Space Telescope (HST). This galaxy harbors, in addition to its outer disk, a very small, stellar disk in its nucleus. Jeans-modeling suggests the presence of a 3-6$\times 10^8 \Msun$ black hole (BH) in the nucleus. This galaxy therefore deviates significantly from the correlation between BH mass and bulge mass, found for other BH candidate galaxies.

The deprojection of the surface brightness of axisymmetric galaxies is indeterminate unless the galaxy is seen edge-on. In practice, this problem is often circumvented by making ad hoc assumptions about the density distribution. However, one can redistribute the density and still project to the same surface brightness. This is similar to adding so-called konus densities to the assumed density distribution. In this paper we investigate the maximum konus density that one can add to elliptical galaxies. In particular we focus on the uncertainties in the central densities of axisymmetric, (cusped) elliptical galaxies due to the non-uniqueness of the deprojection. In order to estimate the dynamical influence of konus densities on the central region, we calculate the mass fraction they can add to the center. We show that this mass fraction is at most a few percent. We also investigate the dynamical effect of cusped konus densities We show that konus densities can only be moderately cusped, and that an increase in cusp slope results in a decrease of the mass fraction added to the center by the konus density. We illustrate all this by the specific example of M32, and show that the uncertainty in the central mass due to deprojection is negligible compared to the inferred mass of the central black hole (BH) in this galaxy.

Recently, observations with the Hubble Space Telescope have revealed small stellar discs embedded in the nuclei of a number of ellipticals and S0s. In this paper we construct two-integral axisymmetric models for such systems. We calculate the even part of the phase-space distribution function, and specify the odd part by means of a simple parameterization. We investigate the photometric as well as the kinematic signatures of nuclear discs, including their velocity profiles (VPs), and study the influence of seeing convolution. The rotation curve of a nuclear disc gives an excellent measure of the central mass-to-light ratio whenever the VPs clearly reveal the narrow, rapidly rotating component associated with the nuclear disc. Steep cusps and seeing convolution both result in central VPs that are dominated by the bulge light, and these VPs barely show the presence of the nuclear disc, impeding measurements of the central rotation velocities of the disc stars. However, if a massive BH is present, the disc component of the VP can be seen in the wing of the bulge part, and measurements of its mean rotation provide a clear signature of the presence of the BH. This signature is insensitive to the uncertainties in the velocity anisotropy, which often lead to ambiguity in the interpretation of a central rise in velocity dispersion as due to a central BH.

We present high spatial resolution ground-based broad-band imaging, H-alpha +[NII] narrow-band imaging and long-slit spectroscopy for the E4 radio galaxy NGC 7052, which has a nuclear dust disc. We detect ionized gas with a LINER spectrum, residing also in a nuclear disc. The gas rotates rapidly and the emission line widths increase towards the nucleus. The images are well fit by an axisymmetric model with the gas and dust in a disc of 1.5" radius (340 pc) in the equatorial plane of the stellar body, viewed at an inclination of 70 degrees. We assume the gas to be on circular orbits in the equatorial plane, with in addition a local turbulent velocity dispersion. The circular velocity is calculated from the combined gravitational potential of the stars and a possible nuclear black hole. The observed gas rotation curve is well fit, either with or without a black hole. Turbulent velocities >300 km/s must be present at radii <0.5" to fit the observed nuclear line widths. Models with a black hole and no turbulence can also fit the line widths, but these models cannot fit the observed nuclear line shapes. Models with both a black hole and gas turbulence can fit the data well, but the black hole is not required by the data, and if present, its mass must be < 5 x 10^8 solar masses. This upper limit is 5 times smaller than the black hole mass inferred for M87 from HST data. Our modelling and analysis techniques will be useful also for the study of the kinematics of nuclear gas discs in other galaxies. Such data will yield important information on the presence of massive black holes.