Dekel & Birnboim06 proposed that the mass-scale that separates late-type and early-type galaxies is linked to the critical halo mass for the propagation of a stable shock and showed that they could reproduce the observed bimodality scale for plausible values of the metallicity of the accreted gas and the shock radius. Here, we take their analysis one step further and present a new semianalytic model that computes the shock radius from first principles. This advancement allows us to compute the critical mass individually for each halo. Separating cold-mode and hot-mode accretion has little effect on the final galaxy masses if feedback does not preferentially couple to the hot gas. We also present an improved model for stellar feedback where 70% of the wind mass is in a cold galactic fountain with a shorter reaccretion timescale at high masses. The latter is the key mechanism that allows us to reproduce the low-mass end of the mass function of galaxies over the entire redshift range 0<z<2.5. Cooling must be mitigated to avoid overpredicting the number density of galaxies with stellar mass greater than a hundred billion Solar masses but is important to form intermediate-mass galaxies. At virial masses greater than three hundred billion Solar masses, cold accretion is more important at high z, where gas is accreted from smaller solid angles, but this is not true at lower masses because high-z filaments have lower metallicities. Our predictions are consistent with the observed metallicity evolution of the intergalactic medium at 0<z<5.

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.

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

GalICS 2.0 is a new semianalytic code to model the formation and evolution of galaxies in a cosmological context. N-body simulations based on a Planck cosmology are used to construct halo merger trees, track subhaloes, compute spins and measure concentrations. The accretion of gas onto galaxies and the morphological evolution of galaxies are modelled with prescriptions derived from hydrodynamic simulations. Star formation and stellar feedback are described with phenomenological models (as in other semianalytic codes). GalICS 2.0 computes rotation speeds from the gravitational potential of the dark matter, the disc and the central bulge. As the rotation speed depends not only on the virial velocity but also on the ratio of baryons to dark matter within a galaxy, our calculation predicts a different Tully-Fisher relation from models in which the rotation speed is proportional to the virial velocity. This is why GalICS 2.0 is able to reproduce the galaxy stellar mass function and the Tully-Fisher relation simultaneously. Our results are also in agreement with halo masses from weak lensing and satellite kinematics, gas fractions, the relation between star formation rate (SFR) and stellar mass, the evolution of the cosmic SFR density, bulge-to-disc ratios, disc sizes and the Faber-Jackson relation.

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.