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.

The difference in stellar metallicity between red and blue galaxies with the same mass constrains the timescale over which red galaxies ceased to form stars. Here we investigate this constraint with the GalICS 2.0 semi-analytic model of galaxy formation. The advantage of this approach is that the time of pericentric passages for satellite galaxies and the mass-loading factor for galactic winds are not free parameters of the chemical evolution model. The former is determined by the N-body simulation used to construct the merger trees, the latter by the requirement that GalICS 2.0 should reproduce the stellar mass function of galaxies. When we compare our theoretical predictions with observations, we find that GalICS 2.0 can reproduce the observed metallicity difference only if quenching is preceded by a burst of star formation, which contributes to the chemical enrichment of the stellar population. Physically, this burst can be explained as tidally-induced star formation or as an effect of ram pressure, which not only strips gas from galaxies but also compresses it, accelerating its conversion into stars.

We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component S\'ersic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers.

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 present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift.

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works.

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

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

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.

For most galaxies, the shutdown of star formation was a slow process that took four billion years. An analysis of thousands of galaxies suggests that 'strangulation' by their environment was the most likely cause. See Letter on page 192 of Peng, Maiolino and Cochrane (2015), Nature, vol. 521.

Roughly half of all stars reside in galaxies without significant ongoing star formation. However, galaxy formation models indicate that it is energetically challenging to suppress the cooling of gas and the formation of stars in galaxies that lie at the centers of their dark matter halos. In this Letter, we show that the dependence of quiescence on black hole and stellar mass is a powerful discriminant between differing models for the mechanisms that suppress star formation. Using observations of 91 star-forming and quiescent central galaxies with directly-measured black hole masses, we find that quiescent galaxies host more massive black holes than star-forming galaxies with similar stellar masses. This observational result is in qualitative agreement with models that assume that effective, more-or-less continuous AGN feedback suppresses star formation, strongly suggesting the importance of the black hole in producing quiescence in central galaxies.

[abridged] We quantify the morphological evolution of z~0 massive galaxies ($M*/M_\odot\sim10^{11}$) from z~3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. The morphologies strongly evolve from z~3. At z<1, the population matches the massive end of the Hubble sequence, with 30% of spheroids, 50% of galaxies with equally dominant disk and bulge components and 20% of disks. At z~2-3 there is a majority of irregular systems (~60-70%) with still 30% of spheroids. We then analyze the SFRs, gas fractions and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around 30-40% were already bulges at z~2.5, with low average SFRs and gas-fractions (10-15%), high Sersic indices (n>3-4) and small effective radii ($R_e$~1 kpc) pointing towards an early formation through gas-rich mergers or VDI. Between z~ 2.5 and z~0, they rapidly increase their size by a factor of ~4-5, become all passive but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex-situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disks systems, essentially happening at z>1. It results in the growth of a significant bulge component (n~3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1/3 keeps a rather small bulge (n~1.5-2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease of the star formation activity and the gas fractions. The growth of the effective radii scales roughly with $H(z)^{-1}$ and it is therefore consistent with the expected growth of disks in galaxy haloes.

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

We studied the chronology of galactic bulge and disc formation by analysing the relative contributions of these components to the B-band rest-frame luminosity density at different epochs. We present the first estimate of the evolution of the fraction of rest-frame B-band light in galactic bulges and discs since redshift z~0.8. We performed a bulge-to-disc decomposition of HST/ACS images of 3266 galaxies in the zCOSMOS-bright survey with spectroscopic redshifts in the range 0.7 < z < 0.9. We find that the fraction of B-band light in bulges and discs is $(26 \pm 4)%$ and $(74 \pm 4)%$, respectively. When compared with rest-frame B-band measurements of galaxies in the local Universe in the same mass range ($10^{9} M_{\odot}\lessapprox M \lessapprox 10^{11.5} M_{\odot}$), we find that the B-band light in discs decreases by ~30% from z~0.7-0.9 to z~0, while the light from the bulge increases by ~30% over the same period of time. We interpret this evolution as the consequence of star formation and mass assembly processes, as well as morphological transformation, which gradually shift stars formed at half the age of the Universe from star-forming late-type/irregular galaxies toearlier types and ultimately into spheroids.

The relation between galaxy luminosity L and halo virial velocity v_vir required to fit the galaxy luminosity function differs from the observed Tully-Fisher relation between L and disc speed v_rot. Hence the problem of reproducing the galaxy luminosity function and the Tully-Fisher relation simultaneously has plagued semianalytic models since their inception. Here we study the relation between v_rot and v_vir by fitting observational average rotation curves of disc galaxies binned in luminosity. We show that the v_rot - v_vir relation that we obtain in this way can fully account for this seeming inconsistency. Therefore, the reconciliation of the luminosity function with the Tully-Fisher relation rests on the complex dependence of v_rot on v_vir, which arises because the ratio of stellar mass to dark matter mass is a strong function of halo mass.

We compare state-of-the-art semi-analytic models of galaxy formation as well as advanced sub-halo abundance matching models with a large sample of early-type galaxies from SDSS at z < 0.3. We focus our attention on the dependence of median sizes of central galaxies on host halo mass. The data do not show any difference in the structural properties of early-type galaxies with environment, at fixed stellar mass. All hierarchical models considered in this work instead tend to predict a moderate to strong environmental dependence, with the median size increasing by a factor of about 1.5-3 when moving from low to high mass host haloes. At face value the discrepancy with the data is highly significant, especially at the cluster scale, for haloes above log Mhalo > 14. The convolution with (correlated) observational errors reduces some of the tension. Despite the observational uncertainties, the data tend to disfavour hierarchical models characterized by a relevant contribution of disc instabilities to the formation of spheroids, strong gas dissipation in (major) mergers, short dynamical friction timescales, and very short quenching timescales in infalling satellites. We also discuss a variety of additional related issues, such as the slope and scatter in the local size-stellar mass relation, the fraction of gas in local early-type galaxies, and the general predictions on satellite galaxies.

The link between massive red galaxies in the local Universe and star-forming galaxies at high redshift is investigated with a semi-analytic model that has proven successful in many ways, e.g. explaining the galaxy colour-magnitude bimodality and the stellar mass-age relation for red-sequence galaxies. The model is used to explore the processes that drive star formation in different types of galaxies as a function of stellar mass and redshift. We find that most z=2-4 star-forming galaxies with M_*>10^10 M_Sun evolve into red-sequence galaxies. Also, most of the massive galaxies on the red-sequence today have passed through a phase of intense star formation at z>2. Specifically, ~ 90% of today's red galaxies with M_*>10^11 M_Sun were fed during this phase by cold streams including minor mergers. Gas-rich major mergers are rare and the effects of merger-driven starbursts are ephemeral. On the other hand, major mergers are important in powering the most extreme starbursts. Gas-rich mergers also explain the tail of intermediate-mass red galaxies that form relatively late, after the epoch of peak star formation. In two thirds of the currently red galaxies that had an intense star formation event at $z<1$, this event was triggered by a merger.

We use both an HI-selected and an optically-selected galaxy sample to directly measure the abundance of galaxies as a function of their "baryonic" mass (stars + atomic gas). Stellar masses are calculated based on optical data from the Sloan Digital Sky Survey (SDSS) and atomic gas masses are calculated using atomic hydrogen (HI) emission line data from the Arecibo Legacy Fast ALFA (ALFALFA) survey. By using the technique of abundance matching, we combine the measured baryonic function (BMF) of galaxies with the dark matter halo mass function in a LCDM universe, in order to determine the galactic baryon fraction as a function of host halo mass. We find that the baryon fraction of low-mass halos is much smaller than the cosmic value, even when atomic gas is taken into account. We find that the galactic baryon deficit increases monotonically with decreasing halo mass, in contrast with previous studies which suggested an approximately constant baryon fraction at the low-mass end. We argue that the observed baryon fractions of low mass halos cannot be explained by reionization heating alone, and that additional feedback mechanisms (e.g. supernova blowout) must be invoked. However, the outflow rates needed to reproduce our result are not easily accommodated in the standard picture of galaxy formation in a LCDM universe.

We introduce a toy model that describes (in a single equation) the mass in stars as a function of halo mass and redshift. Our model includes the suppression of gas accretion from gravitational shock heating and AGN jets mainly for M_halo > M_shock ~ 10^12 M_Sun and from a too hot IGM onto haloes with v_circ < 40 km/s, as well as stellar feedback that drives gas out of haloes mainly with v_circ < 120 km/s. We run our model on the merger trees of the haloes and subhaloes of a high-resolution dark matter cosmological simulation. The galaxy mass is taken as the maximum between the mass given by the model and the sum of the masses of its progenitors (reduced by tidal stripping). Designed to reproduce the present-day stellar mass function of galaxies, our model matches fairly well the evolution of the cosmic stellar density. It leads to the same z=0 relation between central galaxy stellar and halo mass as the one found by abundance matching and also as that previously measured at high mass on SDSS centrals. Our model also predicts a bimodal distribution (centrals and satellites) of stellar masses for given halo mass, in good agreement with SDSS observations. The relative importance of mergers depends much more on stellar than halo mass. Galaxies with m_stars > 10^11 M_Sun/h acquire most of their mass through mergers (mostly major and gas-poor), as expected from our model's shutdown of gas accretion at high M_halo. However, mergers are rare for m_stars < 10^11 M_Sun/h (greater than our mass resolution), a consequence of the curvature of the stellar vs. halo mass relation. So gas accretion must be the dominant growth mechanism for intermediate and low mass galaxies, e.g. dwarf ellipticals in clusters, except that gas-rich galaxy mergers account for the bulk of the growth of ellipticals with m_stars ~ 10^10.5 M_Sun/h, which we predict must be the typical mass of ULIRGs.

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.

Virtually all massive galaxies, including our own, host central black holes ranging in mass from millions to billions of solar masses. The growth of these black holes releases vast amounts of energy that powers quasars and other weaker active galactic nuclei. A tiny fraction of this energy, if absorbed by the host galaxy, could halt star formation by heating and ejecting ambient gas. A central question in galaxy evolution is the degree to which this process has caused the decline of star formation in large elliptical galaxies, which typically have little cold gas and few young stars, unlike spiral galaxies.

We analyze a sample of galaxies chosen to have F24>0.5 mJy and satisfy a certain IRAC color criterion. IRS spectra yield redshifts, spectral types, and PAH luminosities, to which we add multi-wavelength broadband photometry. Stellar population modeling and IRS spectra together demonstrate that the double criteria used to select this sample have efficiently isolated massive star-forming galaxies at z~1.9. This is the first starburst-dominated ULIRG sample at high redshift with total infrared luminosity measured directly from FIR and millimeter photometry, and as such gives us the first accurate view of broadband SEDs for starburst galaxies at extremely high luminosity and at all wavelengths. Similar broadband data are assembled for three other galaxy samples -- local starburst galaxies, local AGN/ULIRGS, and a second 24mu-luminous z~1.9 sample dominated by AGN. L(PAH)/L(IR) for the new z~1.9 starburst sample is the highest ever seen, some three times higher than in local starbursts, whereas in AGNs this ratio is depressed below the starburst trend, often severely. Several pieces of evidence imply that AGNs exist in this starburst dominated sample, except two of which even host very strong AGN, while they still have very strong PAH emission. The ACS images show most objects have very extended morphologies in the rest-frame UV band, thus extended distribution of PAH molecules. Such an extended distribution prevents further destruction PAH molecules by central AGNs. We conclude that objects in this sample are ULIRGs powered mainly by starburst; and the total infrared luminosity density contributed by this type of objects is 0.9-2.6x 10^7 Lsol/Mpc^3.

One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.

Black holes release energy via the production of photons in their accretion discs but also via the acceleration of jets. We investigate the relative importance of these two paths over cosmic time by determining the mechanical luminosity function (LF) of radio sources and by comparing it to a previous determination of the bolometric LF of active galactic nuclei (AGN) from X-ray, optical and infrared observations. The mechanical LF of radio sources is computed in two steps: the determination of the mechanical luminosity as a function of the radio luminosity and its convolution with the radio LF of radio sources. Even with the large uncertainty deriving from the former, we can conclude that the contribution of jets is unlikely to be much larger than 10% of the AGN energy budget at any cosmic epoch.

We have used GADGET2 to simulate the formation of an elliptical galaxy in a cosmological dark matter halo with mass 3x10^12M_Sun/h. Using a stellar population synthesis model has allowed us to compute magnitudes, colours and surface brightness profiles. We have included a model to follow the growth of a central black hole and we have compared the results of simulations with and without feedback from AGNs. We have studied the interplay between cold gas accretion and merging in the development of galactic morphologies, the link between colour and morphology evolution, the effect of AGN feedback on the photometry of early type galaxies, the redshift evolution in the properties of quasar hosts, and the impact of AGN winds on the chemical enrichment of the intergalactic medium (IGM). We have found that the early phases of galaxy formation are driven by the accretion of cold filamentary flows, which form a disc at the centre of the dark matter halo. When the dark matter halo is sufficiently massive to support the propagation of a stable shock, cold accretion is shut down, and the star formation rate begins to decline. Mergers transform the disc into an elliptical galaxy, but also bring gas into the galaxy. Without a mechanism that removes gas from the merger remnants, the galaxy ends up with blue colours, atypical for its elliptical morphology. AGN feedback can solve this problem even with a fairly low heating efficiency. We have also demonstrated that AGN winds are potentially important for the metal enrichment of the IGM a high redshift.(abridged)

We address the origin of the `downsizing' of elliptical galaxies, according to which the stars in more massive galaxies formed earlier and over a shorter period than those in less massive galaxies. We show that this could be the natural result of a shutdown of star formation in dark matter haloes above a critical mass of 10^12MSun. This is demonstrated using a semianalytic simulation of galaxy formation within the standard hierarchical scenario of structure formation. The assumed threshold mass is motivated by the prediction of stable shock heating above this mass and the finding that such a shutdown reproduces the observed distribution of galaxies in luminosity and colour.The shutdown at a critical halo mass introduces a characteristic stellar mass for the transition of galaxies into the `red sequence' of the galaxy colour-magnitude diagram. Central galaxies of haloes that are more massive today have reached this mass earlier and can therefore grow further along the red sequence by dry mergers, ending up more massive and containing older stars. Small galaxies formed in haloes below the critical mass can shutdown late, when they fall into haloes above the critical mass and become satellites. While our semianalytic simulation that incorporates an explicit shutdown reproduces downsizing as inferred from the stellar ages of ellipticals, we explain why it is much harder to detect downsizing using the mass functions of different galaxy types.

We use three-dimensional high-resolution adaptive-mesh-refinement simulations to investigate if mechanical feedback from active galactic nucleus jets can halt a massive cooling flow in a galaxy cluster and give rise to a self-regulated accretion cycle. We start with a 3 x 10^9MSun black hole at the centre of a spherical halo with the mass of the Virgo cluster. Initially, all the baryons are in a hot intracluster medium in hydrostatic equilibrium within the dark matter's gravitational potential. The black hole accretes the surrounding gas at the Bondi rate and a fraction of the accretion power is returned into the intracluster medium mechanically through the production of jets. The accretion, initially slow (~0.0002MSun/yr), becomes catastrophic, as the gas cools and condenses in the dark matter's potential. Therefore, it cannot prevent the cooling catastrophe at the centre of the cluster. However, after this rapid phase, where the accretion rate reaches a peak of ~0.2MSun/yr, the cavities inflated by the jets become highly turbulent. The turbulent mixing of the shock-heated gas with the rest of the intracluster medium puts a quick end to this short-lived rapid-growth phase. After dropping by almost two orders of magnitudes, the black hole accretion rate stabilises at ~0.006MSun/yr, without significant variations for several billions of years, indicating that a self-regulated steady-state has been reached. This accretion rate corresponds to a negligible increase of the black hole mass over the age of the Universe, but is sufficient to create a quasi-equilibrium state in the cluster core.

We compare the galaxy population of an SPH simulation to those predicted by the GalICS semi-analytic model and a stripped down version without supernova and AGN feedback. The SPH simulation and the no-feedback GalICS model make similar predictions for the baryonic mass functions of galaxies and for the dependence of these mass functions on environment and redshift. The two methods also make similar predictions for the galaxy content of dark matter haloes as a function of halo mass and for the gas accretion history of galaxies. Both the SPH and no-feedback GalICS models predict a bimodal galaxy population at z=0. The "red'' sequence of gas poor, old galaxies is populated mainly by satellite systems while, contrary to observations, the central galaxies of massive haloes lie on the "blue'' star-forming sequence as a result of continuing hot gas accretion at late times. Furthermore, both models overpredict the observed baryonic mass function, especially at the high mass end. In the full GalICS model, supernova-driven outflows reduce the masses of low and intermediate mass galaxies by about a factor of two. AGN feedback suppresses gas cooling in large haloes, producing a sharp cut-off in the baryonic mass function and moving the central galaxies of these massive haloes to the red sequence. Our results imply that the observational failings of the SPH simulation and the no-feedback GalICS model are a consequence of missing input physics rather than computational inaccuracies, that truncating gas accretion by satellite galaxies automatically produces a bimodal galaxy distribution with a red sequence, but that explaining the red colours of the most massive galaxies requires a mechanism like AGN feedback that suppresses the accretion onto central galaxies in large haloes.

I consider a toy model of self-regulated black hole accretion. The black hole grows through Bondi accretion and a fraction of the accretion power is distributed as thermal feedback into the surrounding gas. The gas expands or contracts until AGN heating and radiative cooling balance each other. The balance of heating and cooling is used to determine a quasi-equilibrium temperature at which the black hole accretes in self-regulated equilibrium with the surrounding intergalactic medium. This temperature grows with the black hole mass. The temperature increase is very steep around a critical black hole mass due to the shape of the cooling function. The quasi-equilibrium temperature cannot exceed the virial temperature or the AGN will drive a thermal wind. This limits the black hole mass to a maximum value determined by the depth of the potential well. In the regime in which cooling is dominated by bremsstrahlung, this model determines a relation between black hole mass and halo characteristic velocity of the form M_BH ~ v^4. The predictions of the model are consistent with the observed black hole mass -- bulge velocity dispersion relation.

We reproduce the blue and red sequences in the observed joint distribution of colour and magnitude for galaxies at low and high redshifts using hybrid N-body/semi-analytic simulations of galaxy formation. The match of model and data is achieved by mimicking the effects of cold flows versus shock heating coupled to feedback from active galactic nuclei (AGNs), as predicted by Dekel & Birnboim (2006). After a critical epoch z=3, only haloes below a critical shock-heating mass 10^12MSun enjoy gas supply by cold flows and form stars, while cooling and star formation are shut down abruptly above this mass. The shock-heated gas is kept hot because being dilute it is vulnerable to feedback from energetic sources such as AGNs in their self-regulated mode. The shutdown explains in detail the bright-end truncation of the blue sequence at ~L*, the appearance of luminous red-and-dead galaxies on the red sequence starting already at z~2, the colour bimodality, its strong dependence on environment density and its correlations with morphology and other galaxy properties. Before z~2-3, even haloes above the shock-heating mass form stars by cold streams penetrating through the hot gas. This explains the bright star-forming galaxies at z~3-4, the early appearance of massive galaxies on the red sequence, the high cosmological star-formation rate at high redshifts and the subsequent low rate at low redshifts.

The Active Galactic Nuclei In Cosmological Simulations (AGNICS) project includes AGN into the GalICS hybrid (N-body + semi-analytic) model to explore the quasar-galaxy link in a cosmological perspective. The key problems are the quasar fuelling mechanism, the origin of the BH to bulge mass relation, the causal and chronological link between BH growth and galaxy formation, the properties of quasar hosts and the role of AGN feedback in galaxy formation. This first paper has two goals: to describe the general structure and the assumptions of the method, and to apply AGNICS to studying the joint formation of BHs and spheroids in galaxy mergers. We investigate under what conditions this scenario can reproduce the local BH distribution and the evolution of the quasar population. AGNICS contains two star formation modes: a quiescent one in discs and a starburst one in protospheroids, the latter triggered by mergers and disc instabilities. We assume that BH growth is linked to the starburst mode. The simplest version of this scenario, in which the BH accretion rate and the star formation rate in the starburst component are simply related by a constant of proportionality, does not to reproduce the cosmic evolution of the quasar population. A model in which this relation is modulated by a scaling with the square root of the starburst gas density can explain the evolution of the quasar luminosity function in B-band and X-rays (accounting for the obscuration inferred from X-ray studies). The scatter and the tilt that this model introduces in the BH-to-bulge mass relation are within the observational constraints. The quasar contribution grows with bolometric luminosity and for a given bulge mass the most massive BHs are in the bulges with the oldest stars.

We present a method for generating virtual observations from smoothed-particle-hydrodynamics (SPH) simulations. This method includes stellar population synthesis models and the reprocessing of starlight by dust to produce realistic galaxy images. We apply this method and simulate the merging of two identical giant Sa galaxies. The merger remnant is an elliptical galaxy. The merger concentrates the gas content of the two galaxies into the nuclear region. The gas that flows into the nuclear region refuels the central black holes of the merging galaxies. We follow the refuelling of the black holes during the merger semi-analytically. In the simulation presented in this article, the black holes grow from 3 x 10^7 to 1.8X 10^8 Solar masses, with a peak AGN luminosity of M_B ~ -23.7. We study how the morphological and spectral properties of the system evolve during the merger and work out the predictions of this scenario for the properties of host galaxies during the active phase. The peak of AGN activity coincides with the merging of the two galactic nuclei and occurs at a stage when the remnant looks like a lenticular galaxy. The simulation predicts the formation of a circumnuclear starburst ring/dusty torus with an opening angle of 30-40 degrees and made of clouds with n_H=10^24 cm^-2. The average optical depth of the torus is quite high, but the obscuring medium is patchy, so that there still exist lines of sight where the AGN is visible in a nearly edge-on view. For the same reason, there are lines of sight where the AGN is completely obscured in the face-on view.

We investigate the hypothesis that quasars formed together with the stellar populations of early-type galaxies. This hypothesis, in conjunction with the stellar ages of early-type galaxies from population synthesis models, the relation of black hole mass to bulge velocity dispersion, and the velocity dispersion distribution from the Sloan Digital Sky Survey, completely determines the cosmic accretion hystory of supermassive black holes and the redshift evolution of the characteristic luminosity. On the other hand the precise shape of the luminosity function of quasars depends on the light curve of quasars and, in the optical, but not so much in X-rays, on the covering factor of the dust surrounding the AGN. We find a plausible set of assumptions for which the coeval formation of supermassive black holes and elliptical galaxies is in good agreement with the observed B-band and X-ray luminosity functions of quasars.

Quasars are widely believed to be powered by accretion onto supermassive black holes and there is now considerable evidence for a link between mergers, quasars and the formation of spheroids. Cattaneo, Haehnelt & Rees (1999) have demonstrated that a very simple model in which supermassive black holes form and accrete most of their mass in mergers of galaxies of comparable masses can reproduce the observed relation of black hole mass to bulge luminosity. Here we show that this simple model can account for the luminosity function of quasars and for the redshift evolution of the quasar population provided a few additional assumptions are made. We use the extended Press-Schechter formalism to simulate the formation of galaxies in hierarchical models of the formation of structures and we assume that, when two galaxies of comparable masses merge, their central black holes coalesce and a fraction of the gas in the merger remnant is accreted by the supermassive black hole over a time-scale of about 10^7 yr. We find that the decrease in the merging rate with cosmic time and the depletion in the amount of cold gas available due to the formation of stars are not sufficient to explain the strong decline in the space density of bright quasars between z=2 and z=0, since larger and larger structures form, which can potentially host brighter and brighter quasars. To explain the redshift evolution of the space density of bright quasars between z=2 and z=0 we need to assume that there is a dependence on redshift either in the fraction of available gas accreted or in the time-scale for accretion.

The growth of supermassive black holes by merging and accretion in hierarchical models of galaxy formation is studied by means of Monte Carlo simulations. A tight linear relation between masses of black holes and masses of bulges arises if if the mass accreted by supermassive black holes scales linearly with the mass forming stars and if the redshift evolution of mass accretion tracks closely that of star formation. Differences in redshift evolution between black hole accretion and star formation introduce considerable scatter in this relation. A non-linear relation between black hole accretion and star formation results in a non-linear relation between masses of remnant black holes and masses of bulges. The relation of black hole mass to bulge luminosity obseved in nearby galaxies and its scatter are reproduced reasonably well by models in which black hole accretion and star formation are linearly related but do not track each other in redshift. This suggests that a common mechanism determines the efficiency for black hole accretion and the efficiency for star formation, especially for bright bulges.