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

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

In the local Universe, the existence of very young galaxies (VYGs), having formed at least half their stellar mass in the last 1 Gyr, is debated. We predict the present-day fraction of VYGs among central galaxies as a function of galaxy stellar mass. For this, we apply to high mass resolution Monte-Carlo halo merger trees (MCHMTs) three (one) analytical models of galaxy formation, where the ratio of stellar to halo mass (mass growth rate) is a function of halo mass and redshift. Galaxy merging is delayed until orbital decay by dynamical friction. With starbursts associated with halo mergers, our models predict typically one percent of VYGs up to galaxy masses of $10^{10}$ M$_\odot$, falling rapidly at higher masses, and VYGs are usually associated with recent major mergers of their haloes. Without these starbursts, two of the models have VYG fractions reduced by 1 or 2 dex at low or intermediate stellar masses, and VYGs are rarely associated with major halo mergers. In comparison, the state-of-the-art semi-analytical model (SAM) of Henriques et al. produces only 0.01% of VYGs at intermediate masses. Finally, the Menci et al. SAM run on MCMHTs with Warm Dark Matter cosmology generates 10 times more VYGs at masses below $10^8$ M$_\odot$ than when run with Cold Dark Matter. The wide range in these VYG fractions illustrates the usefulness of VYGs to constrain both galaxy formation and cosmological models.

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.

Halo bias is the one of the key ingredients of the halo models. It was shown at a given redshift to be only dependent, to the first order, on the halo mass. In this study, four types of cosmic web environments: clusters, filaments, sheets and voids are defined within a state of the art high resolution $N$-body simulation. Within those environments, we use both halo-dark matter cross-correlation and halo-halo auto correlation functions to probe the clustering properties of halos. The nature of the halo bias differs strongly among the four different cosmic web environments we describe. With respect to the overall population, halos in clusters have significantly lower biases in the {$10^{11.0}\sim 10^{13.5}h^{-1}\rm M_\odot$} mass range. In other environments however, halos show extremely enhanced biases up to a factor 10 in voids for halos of mass {$\sim 10^{12.0}h^{-1}\rm M_\odot$}. Such a strong cosmic web environment dependence in the halo bias may play an important role in future cosmological and galaxy formation studies. Within this cosmic web framework, the age dependency of halo bias is found to be only significant in clusters and filaments for relatively small halos $\la 10^{12.5}\msunh$.

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

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

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

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

Based on the star formation histories (SFH) of galaxies in halos of different masses, we develop an empirical model to grow galaxies in dark mattet halos. This model has very few ingredients, any of which can be associated to observational data and thus be efficiently assessed. By applying this model to a very high resolution cosmological $N$-body simulation, we predict a number of galaxy properties that are a very good match to relevant observational data. Namely, for both centrals and satellites, the galaxy stellar mass function (SMF) up to redshift $z\simeq4$ and the conditional stellar mass functions (CSMF) in the local universe are in good agreement with observations. In addition, the 2-point correlation is well predicted in the different stellar mass ranges explored by our model. Furthermore, after applying stellar population synthesis models to our stellar composition as a function of redshift, we find that the luminosity functions in $^{0.1}u$, $^{0.1}g$, $^{0.1}r$, $^{0.1}i$ and $^{0.1}z$ bands agree quite well with the SDSS observational results down to an absolute magnitude at about -17.0. The SDSS conditional luminosity functions (CLF) itself is predicted well. Finally, the cold gas is derived from the star formation rate (SFR) to predict the HI gas mass within each mock galaxy. We find a remarkably good match to observed HI-to-stellar mass ratios. These features ensure that such galaxy/gas catalogs can be used to generate reliable mock redshift surveys.

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

We propose a common terminology for use in describing both temporal merger trees and spatial structure trees for dark-matter halos. We specify a unified data format in HDF5 and provide example I/O routines in C, FORTRAN and PYTHON.

We use cosmological simulations to study a characteristic evolution pattern of high redshift galaxies. Early, stream-fed, highly perturbed, gas-rich discs undergo phases of dissipative contraction into compact, star-forming systems (blue nuggets) at z~4-2. The peak of gas compaction marks the onset of central gas depletion and inside-out quenching into compact ellipticals (red nuggets) by z~2. These are sometimes surrounded by gas rings or grow extended dry stellar envelopes. The compaction occurs at a roughly constant specific star-formation rate (SFR), and the quenching occurs at a constant stellar surface density within the inner kpc ($\Sigma_1$). Massive galaxies quench earlier, faster, and at a higher $\Sigma_1$ than lower-mass galaxies, which compactify and attempt to quench more than once. This evolution pattern is consistent with the way galaxies populate the SFR-radius-mass space, and with gradients and scatter across the main sequence. The compaction is triggered by an intense inflow episode, involving (mostly minor) mergers, counter-rotating streams or recycled gas, and is commonly associated with violent disc instability. The contraction is dissipative, with the inflow rate >SFR, and the maximum $\Sigma_1$ anti-correlated with the initial spin parameter, as predicted by Dekel & Burkert (2014). The central quenching is triggered by the high SFR and stellar/supernova feedback (possibly also AGN feedback) due to the high central gas density, while the central inflow weakens as the disc vanishes. Suppression of fresh gas supply by a hot halo allows the long-term maintenance of quenching once above a threshold halo mass, inducing the quenching downsizing.

We address the formation of massive stellar spheroids between redshifts $z=4$ and 1 using a suite of AMR hydro-cosmological simulations. The spheroids form as bulges, and the spheroid mass growth is partly driven by violent disc instability (VDI) and partly by mergers. A kinematic decomposition to disc and spheroid yields that the mass fraction in the spheroid is between 50\% and 90\% and is roughly constant in time, consistent with a cosmological steady state of VDI discs that are continuously fed from the cosmic web. The density profile of the spheroid is typically "classical", with a Sersic index $n = 4.5\pm 1$, independent of whether it grew by mergers or VDI and independent of the feedback strength. The disc is characterized by $n=1.5\pm 0.5$, and the whole galaxy by $n=3\pm 1$. The high-redshift spheroids are compact due to the dissipative inflow of gas and the high universal density. The stellar surface density within the effective radius of each galaxy as it evolves remains roughly constant in time after its first growth. For galaxies of a fixed stellar mass, the surface density is higher at higher redshifts.

A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z=0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the Universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z=0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant.

We study the properties of giant clumps and their radial gradients in high-$z$ disc galaxies using AMR cosmological simulations. Our sample consists of 770 snapshots in the redshift range $z=4-1$ from 29 galaxies that at $z=2$ span the stellar mass range $(0.2-3)\times 10^{11}M_{\odot}$. Extended gas discs exist in 83% of the snapshots. Clumps are identified by gas density in 3D and their stellar and dark matter components are considered thereafter. While most of the overdensities are diffuse and elongated, 91% of their mass and 83% of their star-fromation rate (SFR) are in compact round clumps. Nearly all galaxies have a central, massive bulge clump, while 70% of the discs show off-center clumps, 3-4 per galaxy. The fraction of clumpy discs peaks at intermediate disc masses. Clumps are divided based on dark-matter content into $\textit{in-situ}$ and $\textit{ex-situ}$, originating from violent disc instability (VDI) and minor mergers respectively. 60% of the discs are in a VDI phase showing off-center $\textit{in-situ}$ clumps, which contribute 1-7% of the disc mass and 5-45% of its SFR. The $\textit{in-situ}$ clumps constitute 75% of the off-center clumps in terms of number and SFR but only half the mass, each clump containing on average 1% of the disc mass and 6% of its SFR. They have young stellar ages, $100-400 {\rm Myr}$, and high specific SFR (sSFR), $1-10 {\rm Gyr}^{-1}$. They exhibit gradients resulting from inward clump migration, where the inner clumps are somewhat more massive and older, with lower gas fraction and sSFR and higher metallicity. Similar observed gradients indicate that clumps survive outflows. The \exsitu clumps have stellar ages $0.5-3 {\rm Gyr}$ and sSFR $\sim 0.1-2 {\rm Gyr}^{-1}$, and they exhibit weaker gradients. Massive clumps of old stars at large radii are likely \exsitu mergers, though half of them share the disc rotation.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

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

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

We describe simple useful toy models for key processes of galaxy formation in its most active phase, at z > 1, and test the approximate expressions against the typical behaviour in a suite of high-resolution hydro-cosmological simulations of massive galaxies at z = 4-1. We address in particular the evolution of (a) the total mass inflow rate from the cosmic web into galactic haloes based on the EPS approximation, (b) the penetration of baryonic streams into the inner galaxy, (c) the disc size, (d) the implied steady-state gas content and star-formation rate (SFR) in the galaxy subject to mass conservation and a universal star-formation law, (e) the inflow rate within the disc to a central bulge and black hole as derived using energy conservation and self-regulated Q ~ 1 violent disc instability (VDI), and (f) the implied steady state in the disc and bulge. The toy models provide useful approximations for the behaviour of the simulated galaxies. We find that (a) the inflow rate is proportional to mass and to (1+z)^5/2, (b) the penetration to the inner halo is ~50% at z = 4-2, (c) the disc radius is ~5% of the virial radius, (d) the galaxies reach a steady state with the SFR following the accretion rate into the galaxy, (e) there is an intense gas inflow through the disc, comparable to the SFR, following the predictions of VDI, and (f) the galaxies approach a steady state with the bulge mass comparable to the disc mass, where the draining of gas by SFR, outflows and disc inflows is replenished by fresh accretion. Given the agreement with simulations, these toy models are useful for understanding the complex phenomena in simple terms and for back-of-the-envelope predictions.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

A simple, 1-equation, galaxy formation model is applied to both the halo merger tree derived from a high-resolution dissipationless cosmological simulation and to 1/4 million Monte-Carlo halo merger trees. The galaxy formation model involves a sharp entropy barrier against the accretion of gas onto low-mass halos, the shock heating of infalling gas far from the central regions of massive halos, and supernova feedback that drives the gas out of shallow halo potential wells. With the first approach, we show that the large majority of galaxies within group- and cluster-mass halos, known to be mainly dwarf ellipticals, have acquired the bulk of their stellar mass through gas accretion and not via galaxy mergers. With the second approach, we qualitatively reproduce the downsizing trend of greater ages at greater masses in stars and predict an upsizing trend of greater ages as one proceeds to masses lower than 10^10 M_Sun. We find that the fraction of galaxies with very young stellar populations (more than half the stellar mass formed within the last 1.5 Gyr) is a function of present-day stellar mass, which peaks at 0.5% at m_crit=10^7.5-9.5 M_Sun, roughly corresponding to the masses of blue compact dwarfs. We predict that the stellar mass function of galaxies should not show a maximum at m_stars > 10^{5.5}, M_Sun, with a power-law stellar mass function with slope approx -1.6 if the IGM temperature in the outskirts of halos before reionization is set by H2 cooling. We speculate on the nature of the lowest mass galaxies.

We apply a simple, one-equation, galaxy formation model on top of the halos and subhalos of a high-resolution dark matter cosmological simulation to study how dwarf galaxies acquire their mass and, for better mass resolution, on over 10^5 halo merger trees, to predict when they form their stars. With the first approach, we show that the large majority of galaxies within group- and cluster-mass halos have acquired the bulk of their stellar mass through gas accretion and not via galaxy mergers. We deduce that most dwarf ellipticals are not built up by galaxy mergers. With the second approach, we constrain the star formation histories of dwarfs by requiring that star formation must occur within halos of a minimum circular velocity set by the evolution of the temperature of the IGM, starting before the epoch of reionization. We qualitatively reproduce the downsizing trend of greater ages at greater masses and predict an upsizing trend of greater ages as one proceeds to masses lower than m_crit. We find that the fraction of galaxies with very young stellar populations (more than half the mass formed within the last 1.5 Gyr) is a function of present-day mass in stars and cold gas, which peaks at 0.5% at m_crit=10^6-8 M_Sun, corresponding to blue compact dwarfs such as I Zw 18. We predict that the baryonic mass function of galaxies should not show a maximum at masses above 10^5.5, M_Sun, and we speculate on the nature of the lowest mass galaxies.

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

Although a fair amount of work has been devoted to growing Monte-Carlo merger trees which resemble those built from an N-body simulation, comparatively little effort has been invested in quantifying the caveats one necessarily encounters when one extracts trees directly from such a simulation. To somewhat revert the tide, this paper seeks to provide its reader with a comprehensive study of the problems one faces when following this route. The first step to building merger histories of dark matter haloes and their subhaloes is to identify these structures in each of the time outputs (snapshots) produced by the simulation. Even though we discuss a particular implementation of such an algorithm (called AdaptaHOP) in this paper, we believe that our results do not depend on the exact details of the implementation but extend to most if not all (sub)structure finders. We then highlight different ways to build merger histories from AdaptaHOP haloes and subhaloes, contrasting their various advantages and drawbacks. We find that the best approach to (sub)halo merging histories is through an analysis that goes back and forth between identification and tree building rather than one which conducts a straightforward sequential treatment of these two steps. This is rooted in the complexity of the merging trees which have to depict an inherently dynamical process from the partial temporal information contained in the collection of instantaneous snapshots available from the N-body simulation.