It is well known from the Butcher-Oemler effect that galaxies in dense environment are mostly red with little star formation and the fraction of blue galaxies in galaxy groups/clusters also declines rapidly with redshifts. A recent work by Hashimoto et al. reported a local 'blue cluster' with high fraction of blue galaxies ($\sim 0.57$), higher than the model predictions. They ascribed this blue cluster to the feeding of gas along a filamentary structure around the cluster. In this work we use group catalog from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) and the state-of-art of semi-analytic model (SAM) to investigate the formation of blue clusters in local Universe. In total, we find four blue clusters with halo mass $\sim 10^{14}M_{\odot}$ at $0.02 < z < 0.082$, while only the one found by Hashimoto et al. is in a filamentary structure. The SAM predicts that blue clusters have later formation time and most blue satellite galaxies are recently accreted. We conclude that the formation of blue clusters is mainly governed by newly accreted blue satellites, rather than the effect of large-scale environment.

Quenched massive spiral galaxies have attracted great attention recently, as more data is available to constrain their environment and cold gas content. However, the quenching mechanism is still uncertain, as it depends on the mass range and baryon budget of the galaxy. In this letter, we report the identification of a rare population of very massive, quenched spiral galaxies with stellar mass $\gtrsim10^{11}{\rm~M_\odot}$ and halo mass $\gtrsim10^{13}{\rm~M_\odot}$ from the Sloan Digital Sky Survey at redshift $z\sim0.1$. Our CO observations using the IRAM-30m telescope show that these galaxies contain only a small amount of molecular gas. Similar galaxies are also seen in the state-of-the-art semi-analytical models and hydro-dynamical simulations. It is found from these theoretical models that these quenched spiral galaxies harbor massive black holes, suggesting that feedback from the central black holes has quenched these spiral galaxies. This quenching mechanism seems to challenge the popular scenario of the co-evolution between massive black holes and massive bulges.

We present a systematic analysis of the reaction of dark matter distribution to galaxy formation across more than eight orders of magnitude in stellar mass. We extend the previous work presented in the NIHAO-IV paper (Tollet et al.) by adding 46 new high resolution simulations of massive galaxies performed with the inclusion of Black Hole feedback. We show that outflows generated by the AGN are able to partially counteract the dark matter contraction due to the large central stellar component in massive haloes. The net effect is to relax the central dark matter distribution that moves to a less cuspy density profiles at halo masses larger than $\approx 3 \times 10^{12}$ M$_{\odot}$. The scatter around the mean value of the density profile slope ($\alpha$) is fairly constant ($\Delta\alpha \approx 0.3$), with the exception of galaxies with halo masses around $10^{12}$ M$_{\odot}$, at the transition from stellar to AGN feedback dominated systems, where the scatter increases by almost a factor three. We provide useful fitting formulas for the slope of the dark matter density profiles at few percent of the virial radius for the whole stellar mass range: $10^5-10^{12}$ M$_{\odot}$ ($2 \times 10^9-5 \times 10^{13}$ M$_{\odot}$ in halo mass).

A large number of new members ($\sim$150) of the Cetus Stream (CS) were identified from their clustering features in dynamical space using 6D kinematic data by combining LAMOST DR5 and Gaia DR2 surveys. They map a diffuse structure that extends over at least 100 degrees in the northern and southern Galactic hemispheres, at heliocentric distances between 20 to 50 kpc. Taking advantage of this expanded dataset, we model the stream with a suite of tailored N-body simulations. Our findings exclude the possibility that the NGC 5824 globular cluster is the core of the progenitor of the stream, as postulated by previous studies. Our best models, which successfully reproduce the features of the CS indicate that the progenitor is likely a dwarf galaxy of $\sim$ 2$\times$10$^9$M$_{\odot}$, with a diffuse disc morphology. The merger occured $\sim$ 5 Gyr ago and since then it has experienced approximately eight apo-center passages. Our results suggest that NGC 5824 was either a globular cluster situated off-centre in the dwarf progenitor or, alternatively, it was the nuclear star cluster of another dwarf galaxy that has very similar orbit as the progenitor of the CS. In both scenarios, the progenitor systems would leave streams around NGC 5824, but with distinct distance distributions. To discriminate between these scenarios, the detection and accurate distance measurements of the predicted stream around the GC are crucial, which will be possible in the upcoming LSST era. Our simulations also predict that part of the Southern Cetus stream is very likely the newly discovered Palca stream, and possibly related to another, more diffuse Southern substructure, the Eridanus-Pheonix overdensity.

We take advantage of an analytic model of galaxy formation coupled to the merger tree of an N-body simulation to study the roles of environment and stellar mass in the quenching of galaxies. The model has been originally set in order to provide the observed evolution of the stellar mass function as well as reasonable predictions of the star formation rate-stellar mass relation, from high redshift to the present time. We analyse the stellar mass and environmental quenching efficiencies and their dependence on stellar mass, halo mass (taken as a proxy for the environment) and redshift. Our analysis shows that the two quenching efficiencies are redshift, stellar and halo mass dependent, and that the halo mass is also a good proxy for the environment. The environmental quenching increases with decreasing redshift and is inefficient below $\log M_* \sim 9.5$, reaches the maximum value at $\log M_* \sim 10.5$, and decreases again, becoming poorly efficient at very high stellar mass ($\log M_* \gtrsim 11.5$). Central and satellites galaxies are mass quenched differently: for the former, the quenching efficiency depends very weakly on redshift, but strongly on stellar mass; for the latter, it strongly depends on both stellar mass and redshift in the range $10\lesssim \log M_* \lesssim 11$. According to the most recent observational results, we find that the two quenching efficiencies are not separable: intermediate mass galaxies are environmental quenched faster, as well as intermediate/massive galaxies in more massive haloes. At stellar masses lower than $\log M_* \lesssim 9.5$ both quenching mechanisms become inefficient, independently of the redshift.

Theoretical studying of the very inner structure of faint satellite galaxy requires very high-resolution hydro-dynamical simulations with realistic models for star formation, which are beginning to emerge recently. In this work we present an analytical description to model the inner kinematic of satellites in the Milky Way (MW). We use a Monte-Carlo method to produce merger trees for MW mass halo and analytical models to produce stellar mass in the satellite galaxies. We consider two important processes which can significantly modify the inner mass distribution in satellite galaxy. The first is baryonic feedback which can induce a flat inner profile depending on the star formation efficiency in the galaxy. The second is the tidal stripping to reduce and re-distribute the mass inside satellite. We apply this model to MW satellite galaxies in both CDM and thermal relic WDM models. It is found that tidal heating must be effective to produce a relatively flat distribution of the satellite circular velocities, to agree with the data. The constraint on WDM mass depends on the host halo mass. For a MW halo with dark matter mass lower than $2\times 10^{12}M_{\odot}$, a 2 keV WDM model can be safely excluded as the predicted satellite circular velocities are systematically lower than the data. For WDM with mass of 3.5 keV, it requires the MW halo mass to be larger than $1.5\times 10^{12}M_{\odot}$, otherwise the 3.5 Kev model can also be excluded. Our current model can not exclude the WDM model with mass larger than 10 Kev.

We present a Bayesian Stacking technique to directly measure the HI mass function (HIMF) and its evolution with redshift using galaxies formally below the nominal detection threshold. We generate galaxy samples over several sky areas given an assumed HIMF described by a Schechter function and simulate the HI emission lines with different levels of background noise to test the technique. We use Multinest to constrain the parameters of the HIMF in a broad redshift bin, demonstrating that the HIMF can be accurately reconstructed, using the simulated spectral cube far below the HI mass limit determined by the $5\sigma$ flux-density limit, i.e. down to $M_{\rm HI} = 10^{7.5}$ M$_{\odot}$ over the redshift range $0 < z < 0.55$ for this particular simulation, with a noise level similar to that expected for the MIGHTEE survey. We also find that the constraints on the parameters of the Schechter function, $\phi_{\star}$, $M_\star$ and $\alpha$ can be reliably fit, becoming tighter as the background noise decreases as expected, although the constraints on the redshift evolution are not significantly affected. All the parameters become better constrained as the survey area increases. In summary, we provide an optimal method for estimating the HI mass at cosmological distances that allows us to constrain the HI mass function below the detection threshold in forthcoming HI surveys. This study is a first step towards the measurement of the HIMF at high ($z>0.1$) redshifts.

We study the roles of stellar mass and environment in quenching the star formation activity of a large set of simulated galaxies by taking advantage of an analytic model coupled to the merger tree extracted from an N-body simulation. The analytic model has been set to match the evolution of the global stellar mass function since redshift $z\sim 2.3$ and give reasonable predictions of the star formation history of galaxies at the same time. We find that stellar mass and environment play different roles: the star formation rate/specific star formation rate-$M_*$ relations are independent of the environment (defined as the halo mass) at any redshift probed, $0<z<1.5$, for both star forming and quiescent galaxies, while the star formation rate-$M_{halo}$ relation strongly depends on stellar mass in the same redshift range, for both star forming and quiescent galaxies. Moreover, the star formation rate and the specific star formation rate are strongly dependent on stellar mass even when the distance from the cluster core is used as a proxy for the environment, rather than the halo mass. We then conclude that stellar mass is the main driver of galaxy quenching at any redshift probed in this study, not just at $z>1$ as generally claimed, while the environment has a minimal role. All the physical processes linked to the environment must act on very short timescales, such that they do not influence the star formation of active galaxies, but increase the probability of a given galaxy to become quiescent.

Observational work have shown that the two brightest satellite galaxies of the Milky Way (MW), the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC), are rare amongst MW analogues. It is then interesting to know whether the presence of massive satellite has any effect on the whole satellite population in MW analogues. In this article, we investigate this problem using a semi-analytical model combined with the Millennium-II Simulation. MW-analogous galaxies are defined to have similar stellar mass or dark matter halo mass to the MW. We find that, in the first case, the halo mass is larger and there are, on average, twice as many satellites in Milky Way analogs if there is a massive satellite galaxy in the system. This is mainly from the halo formation bias. The difference is smaller if MW analogues are selected using halo mass. We also find that the satellites distribution is slightly asymmetric, being more concentrated on the line connecting the central galaxy and the massive satellite and that, on average, LMC have brought in 14.7 satellite galaxies with $M_{r}<0$ at its accretion, among which 4.5 satellites are still within a distance of 50kpc from the LMC. Considering other satellites, we predict that thereare 7.8 satellites with 50kpc of the LMC. By comparing our model with the early data of Satellites Around Galactic Analogs (SAGA), a survey to observe satellite galaxies around 100 Milky Way analogues, we find that SAGA has more bright satellites and less faint satellites than our model predictions. A future comparison with the final SAGA data is needed.

We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $\alpha$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $\alpha$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $\alpha$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation.

It has been shown, both in simulations and observationally, that the tidal field of a large galaxy can torque its satellites such that the major axis of satellite galaxies points towards their hosts. This so-called `shape alignment' has been observed in isolated Milky Way-like galaxies but not in `Local Group'-like pairs. In this study, we investigate the shape alignment of satellite galaxies in galaxy pairs similar to the Local Group identified in the Sloan Digital Sky Survey Data Release 13 (SDSS DR13). By stacking tens of thousands of satellite galaxies around primary galaxy pairs, we find two statistically strong alignment signals. (1) The major axes of satellite galaxies located in the (projected) area between two primaries (the {\it facing} region) tend to be perpendicular to the line connecting the satellite to its host (tangential alignment), while (2) the major axes of satellite galaxies located in regions away from the other host (the {\it away} region) tend to be aligned with the line connecting the satellite to its host (radial alignment). These alignments are confirmed at $\sim5\sigma$ levels. The alignment signal increases with increasing primary brightness, decreasing pair separation, and decreasing satellite distance. The alignment signal is also found to be stronger in filamentary environments. These findings will shed light on understanding the mechanisms of how satellite galaxies are affected by the tidal field in galaxy pairs and will be useful for investigating galaxy intrinsic alignment in the analyses of weak gravitational lensing.

Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigated redshifts. We further validate the conclusion obtained in Paper I that the gas web is an unbiased tracer of total matter -- even better at high redshifts. By separating the gas mainly by temperature, we find that about 40 per cent of gas is in the so-called warm-hot intergalactic medium (WHIM). This fraction of gas mass in the WHIM decreases with redshift, especially from z = 1 (29 per cent) to z = 2.1 (10 per cent). By separating the whole WHIM gas mass into the four large-scale environments (i.e. voids, sheets, filaments, and knots), we find that about half of the WHIM gas is located in filaments. Although the total gas mass in WHIM decreases with redshift, the WHIM mass fractions in the different LSE seem unchanged.

We have used the NIHAO simulations to explore how supernovae (SNe) affect star formation in galaxies. We find that SN feedback operates on all scales from the interstellar medium (ISM) to several virial radii. SNe regulate star formation by preventing condensation of HI into H$_2$ and by moving cold neutral gas to the hot HII phase. The first effect explains why the cold neutral gas in dwarf galaxies forms stars inefficiently. The second maintains the hot ISM of massive galaxies (HII vents out at lower masses). At $v_{\rm vir}\simeq 67{\rm\,km\,s}^{-1}$, the outflow rate follows the relation: $\dot{M}_{\rm out}=23\,(v_{\rm vir}/67{\rm\,km\,s}^{-1})^{-4.6}\,{\rm SFR}$. $20\%$ to $70\%$ of the gas expelled from galaxies escapes from the halo (ejective feedback) but outflows are dominated by cold swept-up gas, most of which falls back onto the galaxy on a $\sim 1\,$Gyr timescale. This `fountain feedback' reduces the masses of galaxies by a factor of two to four, since gas spends half to three quarter of its time in the fountain. Less than $10\%$ of the ejected gas mixes with the hot circumgalactic medium and this gas is usually not reaccreted. On scales as large as $6r_{\rm vir}$, galactic winds divert the incoming gas from cosmic filaments and prevent if from accreting onto galaxies (pre-emptive feedback). This process is the main reason for the low baryon content of ultradwarves.

We study colors and metallicities of the Brightest Cluster Galaxies (BCGs) and Intra-Cluster Light (ICL) in galaxy groups and clusters, as predicted by a semi-analytic model of galaxy formation, coupled with a set of high-resolution N-body simulations. The model assumes stellar stripping and violent relaxation processes during galaxy mergers to be the main channels for the formation of the ICL. We find that BCGs are more metal-rich and redder than the ICL, at all redshifts since the ICL starts to form ($z\sim 1$). In good agreement with several observed data, our model predicts negative radial metallicity and color gradients in the BCG+ICL system. By comparing the typical colors of the ICL with those of satellite galaxies, we find that mass and metals in the ICL come from galaxies of different mass, depending on the redshift. Stripping of low mass galaxies, $9<\log M_* <10$, is the most important contributor in the early stage of the ICL formation, but the bulk of the mass/metals contents are given by intermediate/massive galaxies, $10<\log M_* <11$, at lower redshift. Our analysis supports the idea that stellar stripping is more important than galaxy mergers in building-up the ICL, and highlights the importance of colors/metallicity measurements for understanding the formation and evolution of the ICL.

We choose the bright compact group HCG 62, which was found to exhibit both excess X-ray emission and high Fe abundance to the southwest of its core, as an example to study the impact of mergers on chemical enrichment in the intragroup medium. We first reanalyze the high-quality Chandra and XMM-Newton archive data to search for the evidence for additional SN II yields, which is expected as a direct result of the possible merger-induced starburst. We reveal that, similar to the Fe abundance, the Mg abundance also shows a high value in both the innermost region and the southwest substructure, forming a high-abundance plateau, meanwhile all the SN Ia and SN II yields show rather flat distributions in $>0.1r_{200}$ in favor of an early enrichment. Then we carry out a series of idealized numerical simulations to model the collision of two initially isolated galaxy groups by using the TreePM-SPH GADGET-3 code. We find that the observed X-ray emission and metal distributions, as well as the relative positions of the two bright central galaxies with reference to the X-ray peak, can be well reproduced in a major merger with a mass ratio of 3 when the merger-induced starburst is assumed. The `best-match' snapshot is pinpointed after the third pericentric passage when the southwest substructure is formed due to gas sloshing. By following the evolution of the simulated merging system, we conclude that the effects of such a major merger on chemical enrichment are mostly restricted within the core region when the final relaxed state is reached.

We study the origin of bulge stars and their angular momentum (AM) evolution in 10 spiral galaxies with baryonic masses above $10^{10}$M$_\odot$ in the NIHAO galaxy formation simulations. The simulated galaxies are in good agreement with observations of the relation between specific AM and mass of the baryonic component and the stellar bulge-to-total ratio ($B/T$). We divide the star particles at $z=0$ into disc and bulge components using a hybrid photometric/kinematic decomposition method that identifies all central mass above an exponential disc profile as the `bulge'. By tracking the bulge star particles back in time, we find that on average 95\% of the bulge stars formed {\it in situ}, 3\% formed {\it ex situ} in satellites of the same halo, and only 2\% formed {\it ex situ} in external galaxies. The evolution of the AM distribution of the bulge stars paints an interesting picture: the higher the final $B/T$ ratio, the more the specific AM remains preserved during the bulge formation. In all cases, bulge stars migrate significantly towards the central region, reducing their average galactocentric radius by roughly a factor 2, independently of the final $B/T$ value. However, in the higher $B/T$ ($\gtrsim0.2$) objects, the velocity of the bulge stars increases and the AM of the bulge is almost conserved, whereas at lower $B/T$ values, the velocity of the bulge stars decreases and the AM of bulge reduces. The correlation between the evolution of the AM and $B/T$ suggests that bulge and disc formation are closely linked and cannot be treated as independent processes.

We analyze 88 independent high-resolution cosmological zoom-in simulations of disk galaxies in the NIHAO simulations suite to explore the connection between the atomic gas fraction and angular momentum of baryons throughout cosmic time. The study is motivated by the analytic model of \citet{obreschkow16}, which predicts a relation between the atomic gas fraction $f_{\rm atm}$ and the global atomic stability parameter $q \equiv j\sigma / (GM)$, where $M$ and $j$ are the mass and specific angular momentum of the galaxy (stars+cold gas) and $\sigma$ is the velocity dispersion of the atomic gas. We show that the simulated galaxies follow this relation from their formation ($z\simeq4$) to present within $\sim 0.5$ dex. To explain this behavior, we explore the evolution of the local Toomre stability and find that $90\%$--$100\%$ of the atomic gas in all simulated galaxies is stable at any time. In other words, throughout the entire epoch of peak star formation until today, the timescale for accretion is longer than the timescale to reach equilibrium, thus resulting in a quasi-static equilibrium of atomic gas at any time. Hence, the evolution of $f_{\rm atm}$ depends on the complex hierarchical growth history primarily via the evolution of $q$. An exception are galaxies subject to strong environmental effects.

The correlation between the spins of dark matter halos and the large-scale structure (LSS) has been studied in great detail over a large redshift range, while investigations of galaxies are still incomplete. Motivated by this point, we use the state-of-the-art hydrodynamic simulation, Illustris-1, to investigate mainly the spin--LSS correlation of galaxies at redshift of $z=0$. We mainly find that the spins of low-mass, blue, oblate galaxies are preferentially aligned with the slowest collapsing direction ($e_3$) of the large-scale tidal field, while massive, red, prolate galaxy spins tend to be perpendicular to $e_3$. The transition from a parallel to a perpendicular trend occurs at $\sim10^{9.4} M_{\odot}/h$ in the stellar mass, $\sim0.62$ in the g-r color, and $\sim0.4$ in triaxiality. The transition stellar mass decreases with increasing redshifts. The alignment was found to be primarily correlated with the galaxy stellar mass. Our results are consistent with previous studies both in N-body simulations and observations. Our study also fills the vacancy in the study of the galaxy spin--LSS correlation at $z=0$ using hydrodynamical simulations and also provides important insight to understand the formation and evolution of galaxy angular momentum.

We propose a new strategy of finding strongly-lensed supernovae (SNe) by monitoring known galaxy-scale strong-lens systems. Strongly lensed SNe are potentially powerful tools for the study of cosmology, galaxy evolution, and stellar populations, but they are extremely rare. By targeting known strongly lensed starforming galaxies, our strategy significantly boosts the detection efficiency for lensed SNe compared to a blind search. As a reference sample, we compile the 128 galaxy-galaxy strong-lens systems from the Sloan Lens ACS Survey (SLACS), the SLACS for the Masses Survey, and the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey. Within this sample, we estimate the rates of strongly-lensed Type Ia SN (SNIa) and core-collapse SN (CCSN) to be $1.23 \pm 0.12$ and $10.4 \pm 1.1$ events per year, respectively. The lensed SN images are expected to be widely separated with a median separation of 2 arcsec. Assuming a conservative fiducial lensing magnification factor of 5 for the most highly magnified SN image, we forecast that a monitoring program with a single-visit depth of 24.7 mag (5$\sigma$ point source, $r$ band) and a cadence of 5 days can detect 0.49 strongly-lensed SNIa event and 2.1 strongly-lensed CCSN events per year within this sample. Our proposed targeted-search strategy is particularly useful for prompt and efficient identifications and follow-up observations of strongly-lensed SN candidates. It also allows telescopes with small field of views and limited time to efficiently discover strongly-lensed SNe with a pencil-beam scanning strategy.

We study the growth pathways of Brightest Central Galaxies (BCGs) and Intra-Cluster Light (ICL) by means of a semi-analytic model. We assume that the ICL forms by stellar stripping of satellite galaxies and violent processes during mergers, and implement two independent models: (1) one considers both mergers and stellar stripping (named {\small STANDARD} model), and one considers only mergers (named {\small MERGERS} model). We find that BCGs and ICL form, grow and overall evolve at different times and with different timescales, but they show a clear co-evolution after redshift $z \sim 0.7-0.8$. Around 90\% of the ICL from stellar stripping is built-up in the innermost 150 Kpc from the halo centre and the dominant contribution comes from disk-like galaxies (B/T$<$0.4) through a large number of small/intermediate stripping events ($M_{strip}/M_{sat}<0.3$). The fractions of stellar mass in BCGs and in ICL over the total stellar mass within the virial radius of the halo evolve differently with time. At high redshift, the BCG accounts for the bulk of the mass, but its contribution gradually decreases with time and stays constant after $z\sim 0.4-0.5$. The ICL, instead, grows very fast and its contribution keeps increasing down to the present time. The {\small STANDARD} and the {\small MERGERS} models make very similar predictions in most of the cases, but predict different amounts of ICL associated to other galaxies within the virial radius of the group/cluster other than the BCG, at $z=0$. We then suggest that this quantity is a valid observable that can shed light on the relative importance of mergers and stellar stripping for the formation of the ICL.

We develop a new method StarGO (Stars' Galactic Origin) to identify the galactic origins of halo stars using their kinematics. Our method is based on self-organizing map (SOM), which is one of the most popular unsupervised learning algorithms. StarGO combines SOM with a novel adaptive group identification algorithm with essentially no free parameters. In order to evaluate our model, we build a synthetic stellar halo from mergers of nine satellites in the Milky Way. We construct the mock catalogue by extracting a heliocentric volume of 10 kpc from our simulations and assigning expected observational uncertainties corresponding to bright stars from Gaia DR2 and LAMOST DR5. We compare the results from StarGO against that from a Friends-of-Friends (FoF) based method in the space of orbital energy and angular momentum. We show that StarGO is able to systematically identify more satellites and achieve higher number fraction of identified stars for most of the satellites within the extracted volume. When applied to data from Gaia DR2, StarGO will enable us to reveal the origins of the inner stellar halo in unprecedented detail.

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

The convergence peaks, constructed from galaxy shape measurement in weak lensing, is a powerful probe of cosmology as the peaks can be connected with the underlined dark matter haloes. However the capability of convergence peak statistic is affected by the noise in galaxy shape measurement, signal to noise ratio as well as the contribution from the projected mass distribution from the large-scale structures along the line of sight (LOS). In this paper we use the ray-tracing simulation on a curved sky to investigate the correspondence between the convergence peak and the dark matter haloes at the LOS. We find that, in case of no noise and for source galaxies at $z_{\rm s}=1$, more than $65\%$ peaks with $\text{SNR} \geq 3$ (signal to noise ratio) are related to more than one massive haloes with mass larger than $10^{13} {\rm M}_{\odot}$. Those massive haloes contribute $87.2\%$ to high peaks ($\text{SNR} \geq 5$) with the remaining contributions are from the large-scale structures. On the other hand, the peaks distribution is skewed by the noise in galaxy shape measurement, especially for lower SNR peaks. In the noisy field where the shape noise is modelled as a Gaussian distribution, about $60\%$ high peaks ($\text{SNR} \geq 5$) are true peaks and the fraction decreases to $20\%$ for lower peaks ($ 3 \leq \text{SNR} < 5$). Furthermore, we find that high peaks ($\text{SNR} \geq 5$) are dominated by very massive haloes larger than $10^{14} {\rm M}_{\odot}$.

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.

The alignment between satellites and central galaxies has been studied in detail both in observational and theoretical works. The widely accepted fact is that the satellites preferentially reside along the major axis of their central galaxy. However, the origin and large-scale environment dependence of this alignment are still unknown. In an attempt to figure out those, we use data constructed from SDSS DR7 to investigate the large-scale environmental dependence of this alignment with emphasis on examining the alignments' dependence on the colour of the central galaxy. We find a very strong large-scale environmental dependence of the satellite-central alignment in groups with blue centrals. Satellites of blue centrals in knots are preferentially located perpendicular to the major axis of the centrals, and the alignment angle decreases with environment namely when going from knots to voids. The alignment angle strongly depend on the ${}^{0.1}(g-r)$ colour of centrals. We suggest that the satellite-central alignment is the result of a competition between satellite accretion within large scale-structure and galaxy evolution inside host haloes. For groups containing red central galaxies, the satellite-central alignment is mainly determined by the evolution effect, while for blue central dominated groups, the effect of large-scale structure plays a more important role, especially in knots. Our results provide an explanation for how the satellite-central alignment forms within different large-scale environments. The perpendicular case in groups and knots with blue centrals may also provide insight into understanding similar polar arrangements such the formation of the Milky Way and Centaurus A's satellite system.

The intra-cluster light (ICL) in observations is usually identified through the surface brightness limit method. In this paper, for the first time we produce the mock images of galaxy groups and clusters using a cosmological hydro- dynamical simulation, to investigate the ICL fraction and focus on its dependence on observational parameters, e.g., the surface brightness limit (SBL), the effects of cosmological redshift dimming, point spread function and CCD pixel size. Detailed analyses suggest that the width of point spread function has a significant effect on the measured ICL fraction, while the relatively small pixel size shows almost no influence. It is found that the measured ICL fraction depends strongly on the SBL. At a fixed SBL and redshift, the measured ICL fraction decreases with increasing halo mass, while with a much faint SBL, it does not depend on halo mass at low redshifts. In our work, the measured ICL fraction shows clear dependence on the cosmological redshift dimming effect. It is found that there are more mass locked in ICL component than light, suggesting that the use of a constant mass-to-light ratio at high surface brightness levels will lead to an underestimate of ICL mass. Furthermore, it is found that the radial profile of ICL shows a characteristic radius which is almost independent of halo mass. The current measurement of ICL from observations has a large dispersion due to different methods, and we emphasize the importance of using the same definition when observational results are compared with the theoretical predictions.

Using a series of cosmological simulations that includes one dark-matter-only (DM-only) run, one gas cooling-star formation-supernovae feedback (CSF) run and one that additionally includes feedback from active galactic nuclei (AGNs), we classify the large-scale structures with both a velocity-shear-tensor code (Vweb) and a tidal-tensor code (Pweb). We find that the baryonic processes have almost no impact on large-scale structures -- at least not when classified using aforementioned techniques. More importantly, our results confirm that the gas component alone can be used to infer the filamentary structure of the Universe practically un-biased, which could be applied to cosmology constrains. In addition, the gas filaments are classified with its velocity (Vweb) and density (Pweb) fields, which can theoretically connect to the radio observations, such as HI surveys. This will help us to bias-freely link the radio observations with DM distributions at large scale.

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

We study the connection between the observed star formation rate-stellar mass (SFR-$M_*$) relation and the evolution of the stellar mass function (SMF) by means of a subhalo abundance matching technique coupled to merger trees extracted from a N-body simulation. Our approach consist of forcing the model to match the observed SMF at redshift $z \sim 2.3$, and let it evolve down to $z \sim 0.3$ according to a $\tau$ model, an exponentially declining functional form which describes the star formation rate decay of both satellite and central galaxies. In this study, we use three different sets of SMFs: ZFOURGE data from Tomczak et al.; UltraVISTA data from Ilbert et al. and COSMOS data from Davidzon et al. We also build a mock survey combining UltraVISTA with ZFOURGE. Our modelling of quenching timescales is consistent with the evolution of the SMF down to $z \sim 0.3$, with different accuracy depending on the particular survey used for calibration. We tested our model against the observed SMFs at low redshift and it predicts residuals (observation versus model) within $1\sigma$ observed scatter along most of the stellar mass range investigated, and with mean residuals below 0.1 dex in the range $\sim [10^{8.7}-10^{11.7}] M_{\odot}$. We then compare the SFR-$M_*$ relation predicted by the model with the observed one at different redshifts. The predicted SFR-$M_*$ relation underpredicts the median SFR at fixed stellar mass relative to observations at all redshifts. Nevertheless, the shapes are consistent with the observed relations up to intermediate-mass galaxies, followed by a rapid decline for massive galaxies.

Both simulations and observations have confirmed that the spin of haloes/galaxies is correlated with the large scale structure (LSS) with a mass dependence such that the spin of low-mass haloes/galaxies tend to be parallel with the LSS, while that of massive haloes/galaxies tend to be perpendicular with the LSS. It is still unclear how this mass dependence is built up over time. We use N-body simulations to trace the evolution of the halo spin-LSS correlation and find that at early times the spin of all halo progenitors is parallel with the LSS. As time goes on, mass collapsing around massive halo is more isotropic, especially the recent mass accretion along the slowest collapsing direction is significant and it brings the halo spin to be perpendicular with the LSS. Adopting the $fractional$ $anisotropy$ (FA) parameter to describe the degree of anisotropy of the large-scale environment, we find that the spin-LSS correlation is a strong function of the environment such that a higher FA (more anisotropic environment) leads to an aligned signal, and a lower anisotropy leads to a misaligned signal. In general, our results show that the spin-LSS correlation is a combined consequence of mass flow and halo growth within the cosmic web. Our predicted environmental dependence between spin and large-scale structure can be further tested using galaxy surveys.

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

NGC\,2403, NGC\,300 and M33 are three nearby pure-disc galaxies with similar stellar mass in different environments, they are benchmarks for understanding late-type spiral galaxies in different environments. The chemical evolution and growth of their discs are investigated by using the simple chemical evolution model, in which their discs are assumed to originate and grow through the accretion of the primordial gas, and the gas outflow process is also taken into account. Through comparative study of the best-fitting model predicted star formation histories for them, we hope to derive a picture of the local environment on the evolution and star formation histories of galaxies and whether or not the isolated galaxies follow similar evolution history. Our results show that these three galaxies accumulated more than 50 percent of their stellar mass at $z\,<\,1$. It can be also found that the metallicity gradients in isolated spiral galaxies NGC\,2403 and NGC\,300 are similar and obviously steeper than that in M33, when the metallicity gradients are expressed in ${\rm dex}\,R_{\rm 25}^{-1}$. The similar metallicity gradients in NGC\,2403 and NGC\,300 indicate that they may experience similar chemical evolutionary histories. The principal epoch of star formation on the discs of NGC\,2403 and NGC\,300 is earlier than that on the disc of M33, and the mean age of stellar populations along the discs of both NGC\,2403 and NGC\,300 is older than that of M33. Our results indicates that the evolution and star formation history of a galaxy indeed depends on its local environment, at least for galaxies with stellar mass of $10^{9.2}\,\rm M_{\odot}\sim10^{9.7}\,\rm M_{\odot}$.

Both simulations and observations have found that the spin of halo/galaxy is correlated with the large scale environment, and particularly the spin of halo flips in filament. A consistent picture of halo spin evolution in different environments is still lacked. Using N-body simulation we find that halo spin with its environment evolves continuously from sheet to cluster, and the flip of halo spin happens both in filament and nodes. For the flip in filament can be explained by halo formation time and the migrating time when its environment changes from sheet to filament. For low-mass haloes, they form first in sheets and migrate into filaments later, so their mass and spin growth inside filament are lower, and the original spin is still parallel to filament. For massive haloes, they migrate into filaments first, and most of their mass and spin growth are obtained in filaments, so the resulted spin is perpendicular to filament. Our results well explain the overall evolution of cosmic web in the cold dark matter model and can be tested using high-redshift data. The scenario can also be tested against alternative models of dark matter, such as warm/hot dark matter, where the structure formation will proceed in a different way.

We use a sample of 83 high-resolution cosmological zoom-in simulations and a semi-analytic model (SAM) to study the stochasticity of galaxy formation in haloes ranging from dwarf to Milky Way masses. Our simulated galaxies reproduce the observed inefficiency of galaxy formation as expressed through the stellar, gas and baryonic Tully-Fisher relations. For HI velocities in the range (70 < V < 220 km/s), the scatter is just 0.08 to 0.14 dex, consistent with the observed intrinsic scatter at these scales. At low velocities (20 < V < 70 km/s), the simulated scatter is 0.2-0.25 dex, which could be tested with future observations. The scatter in the stellar mass versus dark halo velocity relation is constant for 30 < V < 180 km/s, and smaller (~0.17 dex) when using the maximum circular velocity of the dark matter only simulation, Vmax, compared to the virial velocity ($V_{200}$ or $V_{200}^{DMO}$). The scatter in stellar mass is correlated with halo concentration, and is minimized when using a circular velocity at a fixed fraction of the virial radius $\simeq 0.4 R_{200}$ or with $V_{\alpha}=V_{200}^{DMO} (V_{max}^{DMO} / V_{200}^{DMO})^\alpha$ with \alpha ~ 0.7, consistent with constraints from halo clustering. Using the SAM we show the correlation between halo formation time and concentration is essential in order to reproduce this result. This uniformity in galaxy formation efficiency we see in our hydrodynamical simulations and a semi-analytic model proves the simplicity and self-regulating nature of galaxy formation in a \Lambda Cold Dark Matter universe.

We study the connection between the observed star formation rate-stellar mass (SFR-$M_*$) relation and the evolution of the stellar mass function (SMF) by means of a Subhalo Abundance Matching technique coupled to merger trees extracted from a N-body simulation. Our approach, which considers both galaxy mergers and stellar stripping, is to force the model to match the observed SMF at redshift $z>2$, and let it evolve down to the present time according to the observed (SFR-$M_*$) relation. In this study, we use two different sets of SMFs and two SFR-$M_*$ relations: a simple power law and a relation with a mass-dependent slope. Our analysis shows that the evolution of the SMF is more consistent with a SFR-$M_*$ relation with a mass-dependent slope, in agreement with predictions from other models of galaxy evolution and recent observations. In order to fully and realistically describe the evolution of the SMF, both mergers and stellar stripping must be considered, and we find that both have almost equal effects on the evolution of SMF at the massive end. Taking into account the systematic uncertainties in the observed data, the high-mass end of the SMF obtained by considering stellar stripping results in good agreement with recent observational data from the Sloan Digital Sky Survey (SDSS). At $\log M_* < 11.2$, our prediction at z=0.1 is close to \citet{li-white09} data, but the high-mass end ($\log M_* > 11.2$) is in better agreement with \citet{dsouza15} data which account for more massive galaxies.

We use the NIHAO galaxy formation simulations to make predictions for the baryonic budget in present day galaxies ranging from dwarf to Milky Way masses. The sample is made of 88 independent high resolution cosmological zoom-in simulations. NIHAO galaxies reproduce key properties of observed galaxies, such as the stellar mass vs halo mass and cold gas vs stellar mass relations. Thus they make plausible predictions for the baryon budget. We present the mass fractions of stars, cold gas ($T<10^4$K), cool gas ($10^4 < T < 10^5$K), warm-hot gas ($10^5 < T < 5\times10^6$K), and hot gas (T$> 5\times10^6$K) inside the virial radius, $R_{200}$. Compared to the predicted baryon mass, using the dark halo mass and the universal baryon fraction, $f_{\rm b}\equiv \Omega_{\rm b}/\Omega_{\rm m}=0.15$, we find that all of our haloes are missing baryons. The missing mass has been relocated past 2 virial radii, and cool gas dominates the corona at low mass (M$_{200} < 3 \times 10^{11} M_\odot$) while the warm-hot gas dominates at high mass (M$_{200} > 3 \times 10^{11} M_\odot$). Haloes of mass $\sim 10^{10} M_\odot$ are missing $\sim 90\%$ of their baryons. More massive haloes ($\sim 10^{12} M_\odot$) retain a higher fraction of their baryons, with $\sim 30\%$ missing, consistent with recent observational estimates. Moreover, these more massive haloes reproduce the observed fraction of cold, warm-hot and hot gas. The fraction of cool gas we predict ($0.11\pm0.06$) is significantly lower than the observation from COS-HALOs (0.3-0.47), but agrees with the alternative analysis of Stern et al. 2016.

Semi-analytic models of galaxy formation are powerful tools to study the evolution of galaxy population in a cosmological context. However, most models over-predict the number of low-mass galaxies at high redshifts and the color of model galaxies are not right in the sense that low-mass satellite galaxies are too red and centrals are too blue. The recent version of the L-Galaxies model by Henriques et al.(H15) is a step forward to solve these problems by reproducing the evolution of stellar mass function and the overall fraction of red galaxies. In this paper we compare the two model predictions of L-Galaxies (the other is Guo et al. , G13) to the SDSS data in detail. We find that in the H15 model the red fraction of central galaxies now agrees with the data due to their implementation of strong AGN feedback, but the stellar mass of centrals in massive haloes is now slightly lower than the data. For satellite galaxies, the red fraction of low-mass galaxies ($\log M_{*}/M_{\odot} < 10$) also agrees with the data, but the color of massive satellites ($10 < \log M_{*}/M_{\odot} < 11$) is slightly bluer. The correct color of centrals and bluer color of massive satellites indicate that the quenching in massive satellites are not strong enough. We also find that there are too much red spirals and less bulge-dominated galaxies in both H15 and G13 models. Our results suggest that additional mechanisms, such as more minor merger or disk instability, are needed to slightly increase the stellar mass of central galaxy in massive galaxies, mainly in the bulge component, and the bulge dominated galaxy will be quenched or then be quenched by any other mechanisms.

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

We used 87 high resolution hydrodynamical cosmological simulations from the NIHAO suite to investigate the relation between the maximum circular velocity (Vmax) of a dark matter halo in a collisionless simulation and the velocity width of the HI gas in the same halo in the hydrodynamical simulation. These two quantities are normally used to compare theoretical and observational velocity functions and have led to a possible discrepancy between observations and predictions based on the Cold Dark Matter (CDM) model. We show that below 100 km/s, there is clear bias between HI based velocities and Vmax, that leads to an underestimation of the actual circular velocity of the halo. When this bias is taken into account the CDM model has no trouble in reproducing the observed velocity function and no lack of low velocity galaxies is actually present. Our simulations also reproduce the linewidth - stellar mass (Tully-Fisher) relation and HI sizes, indicating that the HI gas in our simulations is as extended as observed. The physical reason for the lower than expected linewidths is that, in contrast to high mass galaxies, low mass galaxies no longer have extended thin HI rotating disks, as is commonly assumed.

We use the NIHAO (Numerical Investigation of Hundred Astrophysical Objects) cosmological simulations to study the effects of galaxy formation on key properties of dark matter (DM) haloes. NIHAO consists of $\simeq 90$ high-resolution SPH simulations that include (metal-line) cooling, star formation, and feedback from massive stars and SuperNovae, and cover a wide stellar and halo mass range: $10^6 < M_* / M_{\odot} < 10^{11}$ ( $10^{9.5} < M_{\rm halo} / M_{\odot} < 10^{12.5}$). When compared to DM-only simulations, the NIHAO haloes have similar shapes at the virial radius, R_{\rm vir}, but are substantially rounder inside $\simeq 0.1R_{\rm vir}$. In NIHAO simulations $c/a$ increases with halo mass and integrated star formation efficiency, reaching $\sim 0.8$ at the Milky Way mass (compared to 0.5 in DM-only), providing a plausible solution to the long-standing conflict between observations and DM-only simulations. The radial profile of the phase-space $Q$ parameter ($\rho/\sigma^3$) is best fit with a single power law in DM-only simulations, but shows a flattening within $\simeq 0.1R_{\rm vir}$ for NIHAO for total masses $M>10^{11} M_{\odot}$. Finally, the global velocity distribution of DM is similar in both DM-only and NIHAO simulations, but in the solar neighborhood, NIHAO galaxies deviate substantially from Maxwellian. The distribution is more symmetric, roughly Gaussian, with a peak that shifts to higher velocities for Milky Way mass haloes. We provide the distribution parameters which can be used for predictions for direct DM detection experiments. Our results underline the ability of the galaxy formation processes to modify the properties of dark matter haloes.

We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(\Omega_b/\Omega_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes.

The Milky Way presents the too-big-to-fail (TBTF) problem that there are two observed satellite galaxies with maximum circular velocity larger than 55km/s, and others have velocity less than 25km/s, but the cold dark matter model predicts there should be more than 10 subhaloes with velocity larger than 25km/s. Those massive subhaloes with $25km/s<V_{max}<55km/s$ should not have failed to form stars. The TBTF problem severely challenges the CDM model. Most efforts are seeking the effects of baryonic feedback, decreasing the mass of the Milky Way, changing the properties of dark matter, so as to assign the observed low-velocity satellites into the massive subhaloes found in simulations. However, the TBTF problem can be avoided if the MW have not accreted subhaloes with velocity between $25km/s<V_{max}<55km/s$ although the probability of such a gap is lower as $\sim 1\%$ and can not be tested against observations. In this work we study the gap in stellar mass of satellite galaxies using the SDSS group catalogue and a semi-analytical model. We find that there are 1-2\% of galaxy groups with a large gap in the stellar mass of their satellites. These 'big gap' groups have accreted less massive subhaloes in their formation history and naturally display a gap between their satellite galaxies. If extrapolating our results to the Milky Way is appropriate, we conclude that it is very likely that our Milky Way has not accreted enough massive subhaloes to host those low-velocity satellites, and the TBTF problem is naturally avoided.

Detailed studies of galaxy formation require clear definitions of the structural components of galaxies. Precisely defined components also enable better comparisons between observations and simulations. We use a subsample of eighteen cosmological zoom-in simulations from the NIHAO project to derive a robust method for defining stellar kinematic discs in galaxies. Our method uses Gaussian Mixture Models in a 3D space of dynamical variables. The NIHAO galaxies have the right stellar mass for their halo mass, and their angular momenta and S\'ersic indices match observations. While the photometric disc-to-total ratios are close to 1 for all the simulated galaxies, the kinematic ratios are around ~0.5. Thus, exponential structure does not imply a cold kinematic disc. Above log(M*)~9.5, the decomposition leads to thin discs and spheroids that have clearly different properties, in terms of angular momentum, rotational support, ellipticity, [Fe/H] and [O/Fe]. At log(M*)<9.5, the decomposition selects discs and spheroids with less distinct properties. At these low masses, both the discs and spheroids have exponential profiles with high minor-to-major axes ratios, i.e. thickened discs.

The quenching of star formation in satellite galaxies is observed over a wide range of dark matter halo masses and galaxy environments. In the recent Guo et al (2011) and Fu et al (2013) semi-analytic + N-body models, the gaseous environment of the satellite galaxy is governed by the properties of the dark matter subhalo in which it resides. This quantity depends of the resolution of the N-body simulation, leading to a divergent fraction of quenched satellites in high- and low-resolution simulations. Here, we incorporate an analytic model to trace the subhaloes below the resolution limit. We demonstrate that we then obtain better converged results between the Millennium I and II simulations, especially for the satellites in the massive haloes ($\rm log M_{halo}=[14,15]$). We also include a new physical model for the ram-pressure stripping of cold gas in satellite galaxies. However, we find very clear discrepancies with observed trends in quenched satellite galaxy fractions as a function of stellar mass at fixed halo mass. At fixed halo mass, the quenched fraction of satellites does not depend on stellar mass in the models, but increases strongly with mass in the data. In addition to the over-prediction of low-mass passive satellites, the models also predict too few quenched central galaxies with low stellar masses, so the problems in reproducing quenched fractions are not purely of environmental origin. Further improvements to the treatment of the gas-physical processes regulating the star formation histories of galaxies are clearly necessary to resolve these problems.

We investigate the galaxy population in simulated proto-cluster regions using a semi-analytic model of galaxy formation, coupled to merger trees extracted from N-body simulations. We select the most massive clusters at redshift $z=0$ from our set of simulations, and follow their main progenitors back in time. The analysis shows that proto-cluster regions are dominated by central galaxies and their number decreases with time as many become satellites, clustering around the central object. In agreement with observations, we find an increasing velocity dispersion with cosmic time, the increase being faster for satellites. The analysis shows that proto-clusters are very extended regions, $\gtrsim 20 \, Mpc$ at $z \gtrsim 1$. The fraction of galaxies in proto-cluster regions that are not progenitor of cluster galaxies varies with redshift, stellar mass and area considered. It is about 20-30 per cent for galaxies with stellar mass $\sim 10^9\,{\rm M}_{\sun}$, while negligible for the most massive galaxies considered. Nevertheless, these objects have properties similar to those of progenitors. We investigate the building-up of the passive-sequence in clusters, and find that their progenitors are on average always active at any redshift of interest of proto-clusters. The main mechanism which quenches their star formation is the removal of the hot gas reservoir at the time of accretion. The later galaxies are accreted (become satellite), and the more the cold gas available, the longer the time spent as active. Central galaxies are active over all redshift range considered, although a non-negligible fraction of them become passive at redshift $z<1$, due to strong feedback from Active Galactic Nuclei.

We compare the half-light circular velocities, V_{1/2}, of dwarf galaxies in the Local Group to the predicted circular velocity curves of galaxies in the NIHAO suite of LCDM simulations. We use a subset of 34 simulations in which the central galaxy has a stellar luminosity in the range 0.5 x 10^5 < L_V < 2 x 10^8 L_{sun}. The NIHAO galaxy simulations reproduce the relation between stellar mass and halo mass from abundance matching, as well as the observed half-light size vs luminosity relation. The corresponding dissipationless simulations over-predict the V_{1/2}, recovering the problem known as too big to fail (TBTF). By contrast, the NIHAO simulations have expanded dark matter haloes, and provide an excellent match to the distribution of V_{1/2} for galaxies with L_V > 2 x 10^6 L_{sun}. For lower luminosities our simulations predict very little halo response, and tend to over predict the observed circular velocities. In the context of LCDM, this could signal the increased stochasticity of star formation in haloes below M_{halo} \sim 10^{10} M_{sun}, or the role of environmental effects. Thus, haloes that are "too big to fail", do not fail LCDM, but haloes that are "too small to pass" (the galaxy formation threshold) provide a future test of LCDM.

Since NGC300 is a bulge-less, isolated low-mass galaxy and has not experienced radial migration during its evolution history, it can be treated as an ideal laboratory to test simple galactic chemical evolution models. By assuming its disk forms gradually from continuous accretion of primordial gas and including the gas-outflow process, we construct a simple chemical evolution model for NGC300 to build a bridge between its SFH and its observed data, especially the present-day radial profiles and global observed properties (e.g., cold gas mass, star-formation rate and metallicity). By means of comparing the model predictions with the corresponding observations, we adopt the classical $\chi^{2}$ methodology to find out the best combination of free parameters $a$, $b$ and $b_{\rm out}$. Our results show that, by assuming an inside-out formation scenario and an appropriate outflow rate, our model reproduces well most of the present-day observational values, not only the radial profiles but also the global observational data for the NGC300 disk. Our results suggest that NGC300 may experience a rapid growth of its disk. Through comparing the best-fitting model predicted SFH of NGC300 with that of M33, we find that the mean stellar age of NGC300 is older than that of M33 and there is a lack of primordial gas infall onto the disk of NGC300 recently. Our results also imply that the local environment may paly a key role in the secular evolution of NGC300.

The distribution of galaxies displays anisotropy on different scales and it is often referred as galaxy alignment. To understand the origin of galaxy alignments on small scales, one must investigate how galaxies were accreted in the early universe and quantify their primordial anisotropic at the time of accretion. In this paper we use N-body simulations to investigate the accretion of dark matter subhaloes, focusing on their alignment with the host halo shape and the orientation of mass distribution on large scale, defined using the hessian matrix of the density field. The large/small (e1/e3) eigenvalues of the hessian matrix define the fast/slow collapse direction of dark matter on large scale. We find that: 1) the halo major axis is well aligned with the e3 (slow collapse) direction, and it is stronger for massive haloes; 2) subhaloes are predominately accreted along the major axis of the host halo, and the alignment increases with the host halo mass. Most importantly, this alignment is universal; 3) accretion of subhaloes with respect to the e3 direction is not universal. In massive haloes, subhaloes are accreted along the e3 (even stronger than the alignment with the halo major axis), but in low-mass haloes subhaloes are accreted perpendicular to the e3. The transit mass is lower at high redshift. The last result well explains the puzzled correlation (both in recent observations and simulations) that massive galaxies/haloes have their spin perpendicular to the filament, and the spin of low-mass galaxies/haloes is slightly aligned with the filament, under the assumption that the orbital angular momentum of subhaloes is converted to halo spin.

We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm. The haloes in our study range from dwarf to Milky Way masses, and represent an unbiased sampling of merger histories, concentrations and spin parameters. The particle masses and force softenings are chosen to resolve the mass profile to below 1% of the virial radius at all masses, ensuring that galaxy half-light radii are well resolved. Using the same treatment of star formation and stellar feedback for every object, the simulated galaxies reproduce the observed inefficiency of galaxy formation across cosmic time as expressed through the stellar mass vs halo mass relation, and the star formation rate vs stellar mass relation. We thus conclude that stellar feedback is the chief piece of physics required to limit the efficiency of star formation in galaxies less massive than the Milky Way.

Taking advantage of a high-resolution simulation coupled with a state-of-art semi-analytic model of galaxy formation, we probe the mass segregation of galaxies in groups and clusters, focusing on which physical mechanisms are driving it. We find evidence of mass segregation in groups and clusters up to the virial radius, both looking at the galaxy stellar mass and subhalo mass. The physical mechanism responsible for that is consistent with dynamical friction, a drag-force that brings more massive galaxies faster towards the innermost regions of the halo. At odds with observational results, we do not find the inclusion of low-mass galaxies in the samples, down to stellar mass $M_* = 10^9 \, M_{\odot}$, to change the overall trend shown by intermediate and massive galaxies. Moreover, stellar stripping as well as the growth of galaxies after their accretion, do not contribute either in shaping mass segregation or mixing the radial mass distribution. Beyond the virial radius we find an "anti-mass segregation" in groups that progressively weakens in clusters. The continuous accretion of new objects and recent merger events play a different role depending on the halo mass onto which accreting material is falling.

By means of our own cosmological-hydrodynamical simulation and semi-analytical model we studied galaxy population properties in clusters and groups, spanning over 10 different bands from UV to NIR, and their evolution since redshift z=2. We compare our results in terms of galaxy red/blue fractions and luminous-to-faint ratio (LFR) on the Red Sequence (RS) with recent observational data reaching beyond z=1.5. Different selection criteria were tested in order to retrieve galaxies belonging to the RS: either by their quiescence degree measured from their specific SFR ("Dead Sequence"), or by their position in a colour-colour plane which is also a function of sSFR. In both cases, the colour cut and the limiting magnitude threshold were let evolving with redshift, in order to follow the natural shift of the characteristic luminosity in the LF. We find that the Butcher-Oemler effect is wavelength-dependent, with the fraction of blue galaxies increasing steeper in optical colours than in NIR. Besides, only when applying a lower limit in terms of fixed absolute magnitude, a steep BO effect can be reproduced, while the blue fraction results less evolving when selecting samples by stellar mass or an evolving magnitude limit. We then find that also the RS-LFR behaviour, highly debated in the literature, is strongly dependent on the galaxy selection function: in particular its very mild evolution recovered when measured in terms of stellar mass, is in agreement with values reported for some of the highest redshift confirmed (proto)clusters. As to differences through environments, we find that normal groups and (to a lesser extent) cluster outskirts present the highest values of both star forming fraction and LFR at low z, while fossil groups and cluster cores the lowest: this separation among groups begins after z~0.5, while earlier all group star forming properties are undistinguishable.

NGC 5194 (M51a) is a grand-design spiral galaxy and undergoing interactions with its companion. Here we focus on investigating main properties of its star-formation history (SFH) by constructing a simple evolution model, which assumes that the disc builds up gradually by cold gas infall and the gas infall rate can be parameterizedly described by a Gaussian form. By comparing model predictions with the observed data, we discuss the probable range for free parameter in the model and then know more about the main properties of the evolution and SFH of M51a. We find that the model predictions are very sensitive to the free parameter and the model adopting a constant infall-peak time $t_{\rm p}\,=\,7.0{\rm Gyr}$ can reproduce most of the observed constraints of M51a. Although our model does not assume the gas infall time-scale of the inner disc is shorter than that of the outer disc, our model predictions still show that the disc of M51a forms inside-out. We find that the mean stellar age of M51a is younger than that of the Milky Way, but older than that of the gas-rich disc galaxy UGC 8802. In this paper, we also introduce a 'toy' model to allow an additional cold gas infall occurred recently to imitate the influence of the interaction between M51a and its companion. Our results show that the current molecular gas surface density, the SFR and the UV-band surface brightness are important quantities to trace the effects of recent interaction on galactic SF process.

Globular Clusters (GCs) are generally treated as natural validators of simple stellar population (SSP) models. However, there are still some differences between real GCs and SSPs. In this work we use a direct $N$-body simulation code {\hf Nbody6} to study the influences of dynamical interactions, metallicity and primordial binaries on Milky Way globular clusters' integrated properties. Our models start with $N=100,000$ stars, covering a metallicity range $Z=0.0001-0.02$, a subset of our models contain primordial binaries, resulting in a binary fraction as currently observed at a model age of GCs.Stellar evolution and external tidal field representative for an average Milky Way GC are taken into consideration. The integrated colours and Lick indices are calculated using BaSeL and Bluered stellar spectral libraries separately. By including dynamical interactions, our model clusters show integrated features (colours up to 0.01mag bluer, H$\beta$ up to 0.1\AA $ $ greater and [MgFe]$'$ 0.05\AA $ $ smaller) making the clusters appear slightly younger than the model clusters without dynamical interactions. This effect is caused mainly by the preferential loss of low-mass stars which have a stronger contribution to redder passbands as well as different spectral features compared to higher-mass stars. In addition, this effect is larger at lower metallicities. On the contrary, the incorporation of primordial binaries reduces this effect.

Radiation fields emitted by O, B-type stars or young stellar populations (SPs) are generally considered as significant central ionizing sources (CISs) of classic H II regions. In our previous studies, we show that the inclusion of binary interactions in stellar population synthesis models can significantly increase the ultraviolet spectrum hardness and the number of ionizing photons of intermediate-age (IA, 7<~log(t/yr)<~8) SPs. In this work, we present photoionization models of H II regions ionized by radiation fields emitted by IA SPs, and show that radiation fields of IA SPs, including binary systems, are in theory possible candidates of significant CISs of classic H II regions. When radiation fields of IA SPs comprising binary systems are used as the CISs of classic H II regions, the theoretical strengths of a number of lines (such as [O III]4959', [S II]6716', etc.), weaker than observations, are raised; the border /selection-criterion lines, between star-forming galaxies and AGNs in the diagnostic diagrams (for example, [N II]6583/Ha~[O III]5007/Hb), move into the region occupied originally by AGNs; and He II 1640 line, observed in Lyman break and high-redshift gravitationally lensed galaxies, also can be produced.

Galaxies are not randomly distributed across the universe but showing different kinds of alignment on different scales. On small scales satellite galaxies have a tendency to distribute along the major axis of the central galaxy, with dependence on galaxy properties that both red satellites and centrals have stronger alignment than their blue counterparts. On large scales, it is found that the major axes of Luminous Red Galaxies (LRGs) have correlation up to 30Mpc/h. Using hydro-dynamical simulation with star formation, we investigate the origin of galaxy alignment on different scales. It is found that most red satellite galaxies stay in the inner region of dark matter halo inside which the shape of central galaxy is well aligned with the dark matter distribution. Red centrals have stronger alignment than blue ones as they live in massive haloes and the central galaxy-halo alignment increases with halo mass. On large scales, the alignment of LRGs is also from the galaxy-halo shape correlation, but with some extent of mis-alignment. The massive haloes have stronger alignment than haloes in filament which connect massive haloes. This is contrary to the naive expectation that cosmic filament is the cause of halo alignment.

Observations have shown that the spatial distribution of satellite galaxies is not random, but rather is aligned with the major axes of central galaxies (CGs). The strength of the alignment is dependent on the properties of both the satellites and centrals. Theoretical studies using dissipationless N-body simulations are limited by their inability to directly predict the shape of CGs. Using hydrodynamical simulations including gas cooling, star formation, and feedback, we carry out a study of galaxy alignment and its dependence on the galaxy properties predicted directly from the simulations.We found that the observed alignment signal is well produced, as is the color dependence: red satellites and red centrals both show stronger alignments than their blue counterparts. The reason for the stronger alignment of red satellites is that most of them stay in the inner region of the dark matter halo where the shape of the CG better traces the dark matter distribution. The dependence of alignment on the color of CGs arises from the halo mass dependence, since the alignment between the shape of the central stellar component and the inner halo increases with halo mass. We also find that the alignment of satellites is most strongly dependent on their metallicity, suggesting that the metallicity of satellites, rather than color, is a better tracer of galaxy alignment on small scales. This could be tested in future observational studies.

Tidally enhanced stellar wind may affect horizontal branch (HB) morphology in globular clusters (GCs) by enhancing the mass loss of primary star during binary evolution. Lei et al. (2013a, 2013b) studied the effect of this kind of wind on HB morphology in details, and their results indicated that binary is a possible secondparameter (2P) candidate in GCs. Binary fraction is an very important fact in the tidally-enhanced-stellar-wind model. In this paper, we studied the effect of binary fraction on HB morphology by removing the effects of metallicity and age. Five different binary fractions (i.e., 10%, 15%, 20%, 30% and 50%) are adopted in our model calculations. The synthetic HB morphologies with different binary fractions are obtained at different metallicities and ages. We found that, due to the great influence of metallicity and age, the effect of binary fraction on HB morphology may be masked by these two parameters. However, when the effects of metallicity and age are removed, the tendency that HB morphologies become bluer with increasing of binary fractions is clearly presented. Furthermore, we compared our results with the observation by Milone et al. (2012). Our results are consistent well with the observation at metalrich and metal-poor GCs. For the GCs with intermediate metallicity, when the effect of age on HB morphology is removed, a weak tendency that HB morphologies become bluer with increasing of binary fractions is presented in all regions of GCs, which is consistent with our results obtained in this metallicity range.

Observations have shown that the spatial distribution of satellite galaxies is not random, but aligned with the major axes of central galaxies. This alignment is dependent on galaxy properties, such that red satellites are more strongly aligned than blue satellites. Theoretical work done to interpret this phenomena has found that it is due to the non-spherical nature of dark matter halos. However, most studies over-predict the alignment signal under the assumption that the central galaxy shape follows the shape of the host halo. It is also not clear whether the color dependence of alignment is due to an assembly bias or an evolution effect. In this paper we study these problems using a cosmological N-body simulation. Subhalos are used to trace the positions of satellite galaxies. It is found that the shape of dark matter halos are mis-aligned at different radii. If the central galaxy shares the same shape as the inner host halo, then the alignment effect is weaker and agrees with observational data. However, it predicts almost no dependence of alignment on the color of satellite galaxies, though the late accreted subhalos show stronger alignment with the outer layer of the host halo than their early accreted counterparts. We find that this is due to the limitation of pure N-body simulations that satellites galaxies without associated subhalos ('orphan galaxies') are not resolved. These orphan (mostly red) satellites often reside in the inner region of host halos and should follow the shape of the host halo in the inner region.

We use high-resolution Eulerian simulations to study the stability of cold gas flows in a galaxy size dark matter halo (10^12 Msun) at redshift z=2. Our simulations show that a cold stream penetrating a hot gaseous halo is stable against thermal convection and Kelvin-Helmholtz instability. We then investigate the effect of a satellite orbiting the main halo in the plane of the stream. The satellite is able to perturb the stream and to inhibit cold gas accretion towards the center of the halo for 0.5 Gyr. However, if the supply of cold gas at large distances is kept constant, the cold stream is able to re-establish itself after 0.3 Gyr. We conclude that cold streams are very stable against a large variety of internal and external perturbations.

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

Galaxy clustering sets strong constraints on the physics governing galaxy formation and evolution. However, most current models fail to reproduce the clustering of low-mass galaxies on small scales ($r<1Mpc/h$). In this paper we study the galaxy clusterings predicted from a few semi-analytical models. We firstly compare two Munich versions, Guo et al. (2011, Guo11) and De Lucia \& Blazoit (2007, DLB07). The Guo11 model well reproduces the galaxy stellar mass function, but over-predicts the clustering of low-mass galaxies on small scales. The DLB07 model provides a better fit to the clustering on small scales, but over-predicts the stellar mass function. These seem to be puzzling. We find that there is slightly more fraction of satellite galaxies residing in massive haloes in the Guo11 model, which is the dominant contribution to the clustering discrepancy between the two models. However, both models still over-predict the clustering at $0.1Mpc/h<r<10Mpc/h$ for low mass galaxies. This is because both models over-predict the number of satellites by $30\%$ in massive halos than the data. Actually, the better agreement of DLB07 model with the data on small scales comes as a coincidence as it predicts too many low-mass central galaxies which are less clustered and thus bring down the total clustering. Finally, we show the predictions from the semi-analytical of Kang et al. (2012). We find that this model can simultaneously fit the stellar mass function and galaxy clustering if the supernova feedback in satellite galaxies is stronger. We conclude that semi-analytical models are now able to solve the small-scales clustering problem, without invoking of any other new physics or changing the dark matter properties, such as the recent favored warm dark matter.

The calibrations of star formation rate (SFR) are prone to be affected by many factors, such as metallicity, initial mass function (IMF), evolutionary population synthesis (EPS) models and so on. In this paper we will discuss the effects of binary interactions, metallicity, EPS models and IMF on several widely used SFR calibrations based on the EPS models of Yunnan with and without binary interactions, BC03, SB99, PEGASE and POPSTAR. The inclusion of binary interactions makes these SFR conversion coefficients smaller (less than 0.2dex), and these differences increase with metallicity. The differences in the calibration coefficient between SFR and the luminosity of $\rm H\alpha$ recombination line (C$_{\rm H\alpha}$) and that between SFR and the ultraviolet (UV) fluxes at 1500 and 2800\,$\rm \AA$ (C$_{i, {\rm UV}}$), caused by IMF, are independent of metallicity (0.03-0.33\,dex) except $\Delta$C$_{\rm H\alpha, IMF}$ when using the POPSTAR and $\Delta$C$_{i, {\rm UV, IMF}}$ when using the PEGASE models. Moreover, we find that $L_{\rm 2800}$ is not suitable to the linear calibration of SFR at low metallicities. At last, we compare the effects of these several factors on the SFR calibrations considered in this paper. The effects of metallicity/IMF and EPS models on the C$_{\rm H\alpha}$ and C$_{\rm FIR}$ (the conversion coefficient between SFR and the far-infrared flux) are the largest among these factors, respectively. For the calibration between SFR and C$_{i, {\rm UV}}$, the effects of these several factors are comparable.

In this paper we study the tidal stripping process for satellite galaxies orbiting around a massive host galaxy, and focus on its dependence on the morphology of both satellite and host galaxy. For this purpose, we use three different morphologies for the satellites: pure disc, pure bulge and a mixture bulge+disc. Two morphologies are used for the host galaxies: bulge+disc and pure bulge. We find that while the spheroidal stellar component experiences a constant power-law like mass removal, the disc is exposed to an exponential mass loss when the tidal radius of the satellite is of the same order of the disc scale length. This dramatic mass loss is able to completely remove the stellar component on time scale of 100 Myears. As a consequence two satellites with the same stellar and dark matter masses, on the same orbit could either retain considerable fraction of their stellar mass after 10 Gyrs or being completely destroyed, depending on their initial stellar morphology. We find that there are two characteristic time scales describing the beginning and the end of the disc removal, whose values are related to the size of the disc. This result can be easily incorporated in semi-analytical models. We also find that the host morphology and the orbital parameters also have an effect on the determining the mass removal, but they are of secondary importance with respect to satellite morphology. We conclude that satellite morphology has a very strong effect on the efficiency of stellar stripping and should be taken into account in modeling galaxy formation and evolution.

In this paper we combine high resolution N-body simulations with a semi analytical model of galaxy formation to study the effects of a possible Warm Dark Matter (WDM) component on the observable properties of galaxies. We compare three WDM models with a dark matter mass of 0.5, 0.75 and 2.0 keV, with the standard Cold Dark Matter case. For a fixed set of parameters describing the baryonic physics the WDM models predict less galaxies at low (stellar) masses, as expected due to the suppression of power on small scales, while no substantial difference is found at the high mass end. However these differences in the stellar mass function, vanish when different set of parameters are used to describe the (largely unknown) galaxy formation processes. We show that is possible to break this degeneracy between DM properties and the parameterization of baryonic physics by combining observations on the stellar mass function with the Tully-Fisher relation (the relation between stellar mass and the rotation velocity at large galactic radii as probed by resolved HI rotation curves). WDM models with a too warm candidate (m<0.75 keV) cannot simultaneously reproduce the stellar mass function and the Tully-Fisher relation. We conclude that accurate measurements of the galaxy stellar mass function and the link between galaxies and dark matter haloes down to the very low-mass end can give very tight constraints on the nature of DM candidates.

An arbitrary surface mass density of gravitational lens can be decomposed into multipole components. We simulate the ray-tracing for the multipolar mass distribution of generalized SIS (Singular Isothermal Sphere) model, based on the deflection angles which are analytically calculated. The magnification patterns in the source plane are then derived from inverse shooting technique. As have been found, the caustics of odd mode lenses are composed of two overlapping layers for some lens models. When a point source traverses such kind of overlapping caustics, the image numbers change by \pm 4, rather than \pm 2. There are two kinds of images for the caustics. One is the critical curve and the other is the transition locus. It is found that the image number of the fold is exactly the average value of image numbers on two sides of the fold, while the image number of the cusp is equal to the smaller one. We also focus on the magnification patterns of the quadrupole (m = 2) lenses under the perturbations of m = 3, 4 and 5 mode components, and found that one, two, and three butterfly or swallowtail singularities can be produced respectively. With the increasing intensity of the high-order perturbations, the singularities grow up to bring sixfold image regions. If these perturbations are large enough to let two or three of the butterflies or swallowtails contact, eightfold or tenfold image regions can be produced as well. The possible astronomical applications are discussed.

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

We build the Yunnan-III evolutionary population synthesis (EPS) models by using the MESA stellar evolution code, BaSeL stellar spectra library and the initial mass functions (IMFs) of Kroupa and Salpeter, and present colours and integrated spectral energy distributions (ISEDs) of solar-metallicity stellar populations (SPs) in the range of 1Myr-15 Gyr. The main characteristic of the Yunnan-III EPS models is the usage of a set of self-consistent solar-metallicity stellar evolutionary tracks (the masses of stars are from 0.1 to 100Msun). This set of tracks is obtained by using the state-of-the-art MESA code. MESA code can evolve stellar models through thermally pulsing asymptotic giant branch (TP-AGB) phase for low- and intermediate-mass stars. By comparisons, we confirm that the inclusion of TP-AGB stars make the V-K, V-J and V-R colours of SPs redder and the infrared flux larger at ages log(t/yr)>7.6 (the differences reach the maximum at log(t/yr)~8.6, ~0.5-0.2mag for colours, ~2 times for K-band flux). The stellar evolutionary tracks, isochrones, colours and ISEDs can be obtained on request from the first author or from our website (http://www1.ynao.ac.cn/~zhangfh/). Using the isochrones, you can build your EPS models. Now the format of stellar evolutionary tracks is the same as that in the STARBURST99 code, you can put them into the STARBURST99 code and get the SP's results. Moreover, the colours involving other passbands or on other systems (for example, HST $F439W-F555W$ colour on AB system) can also be obtained on request.

We construct a parameterized model to explore the main properties of the star formation history of M33. We assume that the disk originates and grows by the primordial gas infall and adopt the simple form of gas accretion rate with one free parameter, the infall time-scale. We also include the contribution of gas outflow process. A major update of the model is that we adopt a molecular hydrogen correlated star formation law and calculate the evolution of the atomic and molecular gas separately. Comparisons between the model predictions and the observational data show that the model predictions are very sensitive to the adopted infall time-scale, while the gas outflow process mainly influences the metallicity profile. The model adopting a moderate outflow rate and an inside-out formation scenario can be in good agreement with most of observed constraints of M33 disk. We also compare the model predictions based on the molecular hydrogen correlated star formation law and that based on the Kennicutt star formation law. Our results imply that the molecular hydrogen correlated star formation law should be preferred to describe the evolution of the M33 disk, especially the radial distributions of both the cold gas and the stellar population.

Using N-body simulations and galaxy formation models, we study the galaxy stellar mass correlation and the two-point auto-correlation. The simulations are run with cosmological parameters from the WMAP first, third and seven year results, which mainly differ in the perturbation amplitude of \sigma_{8}. The stellar mass of galaxies are determined using either a semi-analytical galaxy formation model or a simple empirical abundance matching method. Compared to the SDSS DR7 data at z=0 and the DEEP2 results at z=1, we find that the predicted galaxy clusterings from the semi-analytical model are higher than the data at small scales, regardless of the adopted cosmology. Conversely, the abundance matching method predicts good agreement with the data at both z=0 and z=1 for high \sigma_8 cosmologies (WMAP1 & WMAP7), but the predictions from a low \sigma_8 cosmology (WMAP3) are significantly lower than the data at z=0. We find that the excess clustering at small-scales in the semi-analytical model mainly arises from satellites in massive haloes, indicating that either the star formation is too efficient in low-mass haloes or tidal stripping is too inefficient at high redshift. Our results show that galaxy clustering is strongly affected by the models for galaxy formation, thus can be used to constrain the baryonic physics. The weak dependence of galaxy clustering on cosmological parameters makes it difficult to constrain the WMAP1 and WMAP7 cosmologies.

The number of methods used to study the properties of galaxies is increased, and testing these methods is very important. Galactic globular clusters (GCs) provide an excellent medium for such test, because they can be considered as simple stellar populations. We present ages and metallicities for 40 Galactic GCs as determined from three publicly available techniques, including colour, Lick-index and spectrum-fitting methods, based on Bruzual & Charlot evolutionary population synthesis (EPS) models. By comparing with the ages obtained from colour-magnitude diagrams (CMDs) and metallicities obtained from spectra of stars, we are able to estimate the ability of 'these methods on determination of GCs parameters, which is absolutely necessary. As a result, we find that: (i) for the metallicity, our derived metallicities agree with those derived from the spectra of stars, Lick-index method is suitable to study metallicity for the stellar population systems in the range of -1.5=<[Fe/H]=<-0.7 and spectrum- fitting method is suitable to study metallicity for the stellar population systems in the range of -2.3=<[Fe/H]=<-1.5; (ii) for the age, these three methods have difficulties in age determination, our derived ages are smaller (about 2.0 Gyr, on average) than the results of CMDs for all these three methods. We use Vazdekis and Maraston models to analyze whether our results are dependent on EPS models, and find that the tendency of these two models is the same as that of Bruzual & Charlot models. Our results are independent of the EPS models. In addition, our test is based on the old GCs and our conclusions may hold for old stellar population systems. The whole abstract can be found in my PDF version.

Using the evolutionary population synthesis (EPS) models with and without binary interactions (BIs), we present L_Ha, L_[OII], L_{i,UV} and L_FIR for Burst, S0, Sa-Sd and Irr galaxies, and present the calibrations of star formation rate (SFR) in terms of these diagnostics. By comparison, we find that BIs lower the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors by ~0.2 dex. Moreover, BIs do not significantly vary the calibrations of SFR in terms of L_{i,UV}. In addition, BIs have little effect on the flux at 2800A. At last, the calibration of SFR from L_FIR is almost unaffected by BIs. This is caused by the fact that BIs almost do not affect the bolometric magnitudes of SPs. We also discuss the effects of initial mass function (IMF), gas-recycle assumption (GRA) and EPS models on these SFR calibrations. Comparing the results by using Salpeter (S55) IMF with those by using Miller & Scalo (MS79) IMF, we find that the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors by using S55 IMF are greater by 0.4 and 0.2 dex than those by using MS79 IMF for the models with and without BIs, respectively. The SFR.vs.L_{i,UV} and SFR.vs.L_FIR conversion factors by using S55 IMF are larger by an amount of 0.2 dex than the corresponding ones by using MS79 IMF. The inclusion of GRA only lowers these SFR calibrations at faint SFR. Moreover, comparing the results when using different EPS models, we find that the differences in the SFR.vs.L_Ha and SFR.vs.L_[OII] conversion factors reach ~ 0.7 and 0.9 dex, the difference in the SFR.vs.L_FIR conversion factor reaches 0.4 and 0.8 dex, and the differences in the SFR.vs.L_{i,UV} conversion factors reach 0.3 and 0.2 dex when using S55 and NON-S55 IMF (including Cha03, K01, K93' and MS79 IMFs, partly caused by the difference in the IMF), respectively. At last, we give the conversion coefficients between SFR and these diagnostics for all models.

We provide a comprehensive description and offer an explanation for the sizes of the faintest known galaxies in the universe, the dwarf spheroidal (dSph) satellites of the Milky Way and Andromeda. After compiling a consistent data set of half-light radii (r_{1/2}) and luminosities, we describe the size-luminosity relation of dSphs by a log-normal distribution in r_{1/2} with a mean size that varies as a function of luminosity. Accounting for modest number statistics, measurement uncertainties and surface brightness limitations, we find that the size-luminosity relations of the Milky Way and Andromeda dSph populations are statistically indistinguishable, and also very similar: their mean sizes at a given stellar luminosity differ by no more than 30%. In addition, we find that the mean size, slope and scatter of this log-normal size description of Local Group dSphs matches onto the relation of more massive low-concentration galaxies. This suggests that the stellar sizes of dSphs are ultimately related to their overall initial baryonic angular momentum. To test this hypothesis we perform a series of high resolution N-body simulations that we couple with a semi-analytic model of galaxy formation. These predict the same mean size and slope as observed in dSph satellites. At the same time, these models predict that the size-luminosity distributions for satellite galaxies around similar host-halos must be similar providing a natural explanation as to why the size distributions of Milky Way and Andromeda satellites are similar. Although strong rotation is currently not observed in dSphs, this may well be consistent with our angular-momentum-based explanation for their sizes if the disks of these galaxies have become sufficiently stirred through tidal interaction.

We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Way's halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative "close pair distribution" (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $\rm r_{gc} < 20$ kpc.

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.

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

When a satellite galaxy falls into a massive dark matter halo, it suffers the dynamical friction force which drag it into the halo center and finally it merger with the central galaxy. The time interval between entry and merger is called as the dynamical friction timescale (T_df). Many studies have been dedicated to derive T_df using analytical models or N-body simulations. These studies have obtained qualitative agreements on how T_df depends on the orbit parameters, and mass ratio between satellite and host halo. However, there are still disagreements on the accurate form of T_df . In this paper, we present a semi-analytical model to predict T_df and we focus on interpreting the discrepancies among different studies. We find that the treatment of mass loss from satellite by tidal stripping dominates the behavior of T_df . We also identify other model parameters which affect the predicted T_df.

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

The evolution of the galaxy stellar mass function is especially useful to test the current model of galaxy formation. Observational data have revealed a few inconsistencies with predictions from the $\Lambda {\rm CDM}$ model. For example, most massive galaxies have already been observed at very high redshifts, and they have experienced only mild evolution since then. In conflict with this, semi-analytical models of galaxy formation predict an insufficient number of massive galaxies at high redshift and a rapid evolution between redshift 1 and 0 . In addition, there is a strong correlation between star formation rate and stellar mass for star-forming galaxies, which can be roughly reproduced with the model, but with a normalization that is too low at high redshift. Furthermore, the stellar mass density obtained from the integral of the cosmic star formation history is higher than the measured one by a factor of 2. In this paper, we study these issues using a semi-analytical model that includes: 1) cold gas accretion in massive halos at high redshift; 2) tidal stripping of stellar mass from satellite galaxies; and 3) an evolving stellar initial mass function (bottom-light) with a higher gas recycle fraction. Our results show that the combined effects from 1) and 2) can predict sufficiently massive galaxies at high redshifts and reproduce their mild evolution at low redshift, While the combined effects of 1) and 3) can reproduce the correlation between star formation rate and stellar mass for star-forming galaxies across wide range of redshifts. A bottom-light/top-heavy stellar IMF could partly resolve the conflict between the stellar mass density and cosmic star formation history.

We study the luminosity function and the radial distribution of satellite galaxies within Milky Way sized haloes as predicted in Cold Dark Matter based models of galaxy formation, making use of numerical N-body techniques as well as three different semi-analytic model (SAMs) galaxy formation codes. We extract merger trees from very high-resolution dissipationless simulations of four Galaxy-sized DM haloes, and use these as common input for the semi-analytic models. We present a detailed comparison of our predictions with the observational data recently obtained on the Milky Way satellite luminosity function (LF). We find that semi-analytic models with rather standard astrophysical ingredients are able to reproduce the observed luminosity function over six orders of magnitude in luminosity, down to magnitudes as faint as M_V=-2. We also perform a comparison with the actual observed number of satellites as a function of luminosity, by applying the selection criteria of the SDSS survey to our simulations instead of correcting the observations for incompleteness. Using this approach we again find good agreement for both the luminosity and radial distributions of MW satellites. We investigate which physical processes in our models are responsible for shaping the predicted satellite LF, and find that tidal destruction, suppression of gas infall by a photo-ionizing background, and supernova feedback all make important contributions. We conclude that the number and luminosity of Milky Way satellites can be naturally accounted for within the (Lambda) Cold Dark Matter paradigm, and this should no longer be considered a problem.

We use galaxy catalogues constructed by combining high-resolution N-body simulations with semi-analytic models of galaxy formation to study the properties of Long Gamma-Ray Burst (LGRB) host galaxies. We assume that LGRBs originate from the death of massive young stars and analyse how results are affected by different metallicity constraints on the progenitor stars. As expected, the host sample with no metallicity restriction on the progenitor stars provides a perfect tracer of the cosmic star formation history. When LGRBs are required to be generated by low-metallicity stars, they trace a decreasing fraction of the cosmic star formation rate at lower redshift, as a consequence of the global increase in metallicity. We study the properties of host galaxies up to high redshift (~9), finding that they typically have low-metallicity (Z<0.5 Z_sun) and that they are small (M<10^9 M_sun), bluer and younger than the average galaxy population, in agreement with observational data. They are also less clustered than typical L_* galaxies in the Universe, and their descendents are massive, red and reside in groups of galaxies with halo mass between 10^{13} M_sun to 10^{14} M_sun.

It has been pointed out that the Galactic satellites all have a common mass around 1e7 Msun within 300 pc (M03), while they span almost four order of magnitudes in luminosity (Mateo et al. 1993, Strigari et al. 2008). It is argued that this may reflect a specific scale for galaxy formation or a scale for dark matter clustering. Here we use numerical simulations coupled with a semi-analytic model for galaxy formation, to predict the central mass and luminosity of galactic satellites in the LCDM model. We show that this common mass scale can be explained within the Cold Dark Matter scenario when the physics of galaxy formation is taken into account. The narrow range of M03 comes from the narrow distribution of circular velocities at time of accretion (peaking around 20 km/s) for satellites able to form stars and the not tight correlation between halo concentration and circular velocity. The wide range of satellite luminosities is due to a combination of the mass at time of accretion and the broad distribution of accretion redshifts for a given mass. This causes the satellites baryonic content to be suppressed by photo-ionization to very different extents. Our results favor the argument that the common mass M03 reflects a specific scale (circular velocity ~ 20 km/s) for star formation.

We revisit the Milky Way satellite problem using a semi-analytical model of galaxy formation and compare the predicted luminosity function to recent result from the SDSS. With cosmic photoionization, the luminosity function can be brought into broad agreement with the data between $-15< M_{V} <-2$. This improvement over previous semi-analytical model results (e.g., Benson et al.2002) is from our adoption of improved models for galaxy merger history and galaxy merging time-scales. The very faint satellites ($M_{v} > -5$) formed in halos with virial temperature over $10^{4}K$ (mass around $10^{9} M_{\odot}$ before accretion), but their baryon content are strongly suppressed by photoionization. We model the mass evolution of the subhalos, and compare the predicted mass-to-light ratio with the data. We find that the measured total mass inside the luminous radii of satellites are about 5% of their present total dark matter mass.

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

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.

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.

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

The first year Wilkinson Microwave Anisotropy Probe (WMAP) set quantitative constraints on the amplitude of any primordial non-Gaussianity. We run a series of dark matter-only N-body simulations with the WMAP constraints to investigate the effect of the presence of primordial non-Gaussianity on large scale structures. The model parameters can be constrained using the observations of protoclusters associated with Ly-$\alpha$ emitters at high redshift ($2 \leq z \leq 4$), assuming the galaxy velocity bias can be modelled properly. High redshift structure formation potentially provides a more powerful test of possible primordial non-Gaussianity than does the CMB, albeit on smaller scales. Another constraint is given by the local galaxy density probability distribution function (PDF), as mapped by the 2 degree Field Galaxy Redshift Survey (2dFGRS). The PDF of 2dFGRS \lstar galaxies is substantially higher than the standard model predictions and requires either a non-negligible bias between galaxy and dark matter on $\sim 12$~\hmpc scales or a stronger non-Gaussianity than allowed by the WMAP year one data. The latter interpretation is preferred since second-order bias corrections are negative. With a lower normalisation of the power spectrum fluctuations, sigma_8=0.74, as favoured by the WMAP 3 year data, the discrepancy between the Gaussian model and the data is even larger.

Semi-analytical models of galaxy formation can be used to predict the evolution of the number density of early-type galaxies as a function of the circular velocity at the virial radius, v_{c,vir}. Gravitational lensing probability and separation distribution on the other hand are sensitive to the velocity dispersion (or circular velocity) at about the effective radius. We adopt the singular isothermal ellipsoid (SIE) lens model to estimate the velocity dispersion at the effective radius. We use radio lenses from the Cosmic Lens All-Sky Survey and the PMN-NVSS Extragalactic Lens Survey to study how the velocity dispersions, \sigma_SIE, are related to v_{c,vir}. When we include both the lensing probability and separation distribution as our lensing constraints, we find \sigma_SIE /(200\kms) = [(1.17_{-0.26}^{+0.40}) v_{c,vir}/ (200\kms)]^{0.22^{+0.05}_{-0.04}} for 200\kms \la \sigma_SIE \la 260\kms; at \sigma_SIE = 200\kms, the ratio \sqrt{2} \sigma_SIE / v_{c,vir} is about 1.65^{+0.57}_{-0.37} (68% CL) but decreases to 0.65_{-0.12}^{+0.15} (68% CL) for \sigma_SIE = 260\kms. These results are consistent with those of Seljak (2002) obtained from galaxy-galaxy weak lensing for galaxies of around L_*. However, our results clearly suggest that the ratio must vary significantly as \sigma_SIE is varied and are marginally discrepant with the Seljak results at \sigma_SIE = 260\kms. The scaling \sigma_SIE \propto v_{c,vir}^{0.22\pm 0.05} is broadly consistent with those from galaxy occupation statistics studies and the most recent galaxy-galaxy weak lensing study. (Abridged)

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

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

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 develop an algorithm to reconstruct the triaxial shapes of dark matter halos from the anisotropic spatial distributions of their substructures for the concordance background cosmology. First, we construct an analytic model for the anisotropic spatial distribution of dark halo substructures under the assumption that the tidal field with non-zero trace in the triaxial mass distribution of the host halos generates the substructure bulk motions toward the major principal axes of the hosts. Our analytic model implies that the degree of anisotropy depends sensitively on the triaxiality of the host halos as well as the correlation between the substructure locations and the tidal shear field. Second, we set the axis ratios of the triaxial host halos as free parameters in the analytic model, fit the model to the numerical results from high-resolution N-body simulation of the concordance cosmology, and reconstruct the two axis-ratios of the host halos from the best-fit values of the free parameters. The comparison of the reconstructed axis ratios with the numerical results reveals a good agreement. Finally, we conclude that our analytic model may provide a physical understanding of the anisotropic spatial distribution of dark halo substructures and a new way to reconstruct in principle the triaxial shapes of dark matter halos from the observables.

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

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

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

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

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

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