Search

We present the MOdel for the Rise of GAlaxies aNd Active nuclei (MORGANA). Starting from the merger trees of dark matter halos and a model for the evolution of substructure within the halos, the complex physics of baryons is modeled with a set of state-of-the-art models that describe the mass, metal and energy flows between the various components and phases of a galaxy. The processes of shock-heating and cooling, star formation, feedback, galactic winds and super-winds, accretion onto BHs and AGN feedback are described by new models. In particular, the evolution of the halo gas explicitly follows the thermal and kinetic energies of the hot and cold phases, while star formation and feedback follow the results of the multi-phase model by Monaco (2004a). The increased level of sophistication allows to move from a phenomenological description of gas physics, based on simple scalings with the depth of the DM halo potential, toward a fully physically motivated one. The comparison of the predictions of MORGANA with a basic set of galactic data reveals from the one hand an overall rough agreement, and from the other hand highlights a number of well- or less-known problems: (i) producing the cutoff of the luminosity function requires to force the quenching of the late cooling flows by AGN feedback, (ii) the normalization of the Tully-Fisher relation of local spirals cannot be recovered unless the dark matter halos are assumed to have a very low concentration, (iii) the mass function of HI gas is not easily fitted at small masses, unless a similarly low concentration is assumed, (iv) there is an excess of small elliptical galaxies at z=0. These discrepancies, more than the points of agreement with data, give important clues on the missing ingredients of galaxy formation. (ABRIDGED), Comment: 35 pages, figures included, uses mn2e.cls. Revised cooling model, results are slightly changed, conclusions are unchanged. MNRAS, in press

We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas accretion onto supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted onto the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disc instabilities in star forming discs. We compare our predictions with the luminosity function (LF) observed for Active Galactic Nuclei (AGN) and we confirm that a non-instantaneous accretion timescale (either in the form of a low-angular momentum reservoir or as an assumed light curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies or from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies ($> 10^{11} M_\odot$) without changing their overall assembly history. These winds also affect the SFR of lower mass galaxies, resulting in a too large fraction of passive galaxies at $< 10^{10} M_\odot$. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than $10^8 M_\odot$ are already in a self-regulated state as inferred from observations., 19 pages, 12 figures, 1 table, replaced with MNRAS accepted version

We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs., 9 pages; 3 figures

We present results of cosmological simulations of disc galaxies carried out with the GADGET-3 TreePM+SPH code, where star formation and stellar feedback are described using our MUlti Phase Particle Integrator model. This description is based on a simple multiphase model of the interstellar medium at unresolved scales, where mass and energy flows among the components are explicitly followed by solving a system of ordinary differential equations. Thermal energy from supernovae is injected into the local hot phase, so as to avoid that it is promptly radiated away. A kinetic feedback prescription generates the massive outflows needed to avoid the overproduction of stars. We use two sets of zoomed-in initial conditions of isolated cosmological haloes with masses (2-3) × 1012 M☉, both available at several resolution levels. In all cases we obtain spiral galaxies with small bulge-over-total stellar mass ratios (B/T ̃ 0.2), extended stellar and gas discs, flat rotation curves and realistic values of stellar masses. Gas profiles are relatively flat, molecular gas is found to dominate at the centre of galaxies, with star formation rates following the observed Schmidt-Kennicutt relation. Stars kinematically belonging to the bulge form early, while disc stars show a clear inside-out formation pattern and mostly form after redshift z = 2. However, the baryon conversion efficiencies in our simulations differ from the relation given by Moster et al. at a 3σ level, thus indicating that our stellar discs are still too massive for the dark matter halo in which they reside. Results are found to be remarkably stable against resolution. This further demonstrates the feasibility of carrying out simulations producing a realistic population of galaxies within representative cosmological volumes, at a relatively modest resolution.

Recent inspections of local available data suggest that the almost linear relation between the stellar mass of spheroids ($M_{\rm sph}$) and the mass of the super massive Black Holes (BHs) residing at their centres, shows a break below $M_{\rm sph} \sim 10^{10}\ {\rm M}_\odot$, with a steeper, about quadratic relation at smaller masses. We investigate the physical mechanisms responsible for the change in slope of this relation, by comparing data with the results of the semi-analytic model of galaxy formation MORGANA, which already predicted such a break in its original formulation. We find that the change of slope is mostly induced by effective stellar feedback in star-forming bulges. The shape of the relation is instead quite insensitive to other physical mechanisms connected to BH accretion such as disc instabilities, galaxy mergers, Active Galactic Nucleus (AGN) feedback, or even the exact modelling of accretion onto the BH, direct or through a reservoir of low angular momentum gas. Our results support a scenario where most stars form in the disc component of galaxies and are carried to bulges through mergers and disc instabilities, while accretion onto BHs is connected to star formation in the spheroidal component. Therefore, a model of stellar feedback that produces stronger outflows in star-forming bulges than in discs will naturally produce a break in the scaling relation. Our results point to a form of co-evolution especially at lower masses, below the putative break, mainly driven by stellar feedback rather than AGN feedback., MNRAS accepted, 10 pages, 6 figures, 1 table

[Abridged] In this paper, we carry out a detailed analysis of the performance of two different methods to identify the diffuse stellar light in cosmological hydrodynamical simulations of galaxy clusters. One method is based on a dynamical analysis of the stellar component. The second method is closer to techniques commonly employed in observational studies. Both the dynamical method and the method based on the surface brightness limit criterion are applied to the same set of hydrodynamical simulations for a large sample about 80 galaxy clusters. We find significant differences between the ICL and DSC fractions computed with the two corresponding methods, which amounts to about a factor of two for the AGN simulations, and a factor of four for the CSF set. We also find that the inclusion of AGN feedback boosts the DSC and ICL fractions by a factor of 1.5-2, respectively, while leaving the BCG+ICL and BCG+DSC mass fraction almost unchanged. The sum of the BCG and DSC mass stellar mass fraction is found to decrease from ~80 per cent in galaxy groups to ~60 per cent in rich clusters, thus in excess of what found from observational analysis. We identify the average surface brightness limits that yields the ICL fraction from the SBL method close to the DSC fraction from the dynamical method. These surface brightness limits turn out to be brighter in the CSF than in the AGN simulations. This is consistent with the finding that AGN feedback makes BCGs to be less massive and with shallower density profiles than in the CSF simulations. The BCG stellar component, as identified by both methods, are slightly older and more metal-rich than the stars in the diffuse component., 18 Pages, 15 figures. Matches to MNRAS published version

We use a sample of 1669 QSOs ($r, Comment: 8 pages, 9 figures, MNRAS accepted

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., La lista completa de autores que integran el documento puede consultarse en el archivo., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 at z~0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from the CLASH-VLT survey. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as part of our ESO Large Programme, sampling the cluster out to a radius ~2R200 (4 Mpc). We classify 412 galaxies as passive, with strong Hdelta absorption (red and blue galaxies, and with emission lines from weak to very strong. A number of tests for substructure detection are applied to analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hdelta galaxies trace the cluster center and the WNW-ESE elongated structure. The red strong Hdelta galaxies also mark a secondary, dense peak ~2 Mpc at ESE. The emission line galaxies cluster in several loose structures, mostly outside R200. The observational scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red strong Hdelta galaxies. The red strong Hdelta galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions., Comment: A&A accepted, 19 pages, 30 figures, minor language changes

A well calibrated method to describe the environment of galaxies at all redshifts is essential for the study of structure formation. Such a calibration should include well understood correlations with halo mass, and the possibility to identify galaxies which dominate their potential well (centrals), and their satellites. Focusing on z = 1 and 2 we propose a method of environmental calibration which can be applied to the next generation of low to medium resolution spectroscopic surveys. Using an up-to-date semi-analytic model of galaxy formation, we measure the local density of galaxies in fixed apertures on different scales. There is a clear correlation of density with halo mass for satellite galaxies, while a significant population of low mass centrals is found at high densities in the neighbourhood of massive haloes. In this case the density simply traces the mass of the most massive halo within the aperture. To identify central and satellite galaxies, we apply an observationally motivated stellar mass rank method which is both highly pure and complete, especially in the more massive haloes where such a division is most meaningful. Finally we examine a test case for the recovery of environmental trends: the passive fraction of galaxies and its dependence on stellar and halo mass for centrals and satellites. With careful calibration, observationally defined quantities do a good job of recovering known trends in the model. This result stands even with reduced redshift accuracy, provided the sample is deep enough to preserve a wide dynamic range of density., Comment: 19 pages, 12 figures, accepted for publication in MNRAS

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., 27 pages, 14 figures, 2 tables. MNRAS, in press

We compare predictions of cooled masses and cooling rates from three stripped-down Semi-Analytic Models (SAMs) of galaxy formation with the results of N-body+SPH simulations with gas particle mass of 3.9x10^6 Msun, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ~2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z=0 to within ~20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved halos is subtracted out, SAMs tend to under-predict by ~30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses SAMs tend to over-predict cooling rates, though the numerical result may be affected by the use of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z~2, cold gas accumulates very quickly in the simulated halos. This accumulation is reproduced by SAMs with varying degrees of accuracy., Comment: 20 pages, 12 figures, accepted by MNRAS

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8, Comment: 32 pages, 20 figures, 9 tables. Published in MNRAS. Replaced Fig. 8 (small scaling bug in the previous version)

We present results from high--resolution cosmological hydrodynamical simulations of a Milky--Way-sized halo, aimed at studying the effect of feedback on the nature of gas accretion. Simulations include a model of inter-stellar medium and star formation, in which SN explosions provide effective thermal feedback. We distinguish between gas accretion onto the halo, which occurs when gas particles cross the halo virial radius, and gas accretion onto the central galaxy, which takes place when gas particles cross the inner one-tenth of the virial radius. Gas particles can be accreted through three different channels, depending on the maximum temperature value, $T_{\rm max}$, reached during the particles' past evolution: a cold channel for $T_{\rm max}10^6$K, and a warm one for intermediate values of $T_{\rm max}$. We find that the warm channel is at least as important as the cold one for gas accretion onto the central galaxy. This result is at variance with previous findings that the cold mode dominates gas accretion at high redshift. We ascribe this difference to the different supernova feedback scheme implemented in our simulations. While results presented so far in the literature are based on uneffective SN thermal feedback schemes and/or the presence of a kinetic feedback, our simulations include only effective thermal feedback. We argue that observational detections of a warm accretion mode in the high--redshift circum-galactic medium would provide useful constraints on the nature of the feedback that regulates star formation in galaxies., 6 pages, 3 figures, accepted for publication in ApJL

We present evidence that 'bona fide' disks and starburst systems occupy distinct regions in the gas mass versus star formation (SF) rate plane, both for the integrated quantities and for the respective surface densities. This result is based on CO observations of galaxy populations at low and high redshifts, and on the current consensus for the CO luminosity to gas mass conversion factors. The data suggest the existence of two different star formation regimes: a long-lasting mode for disks and a more rapid mode for starbursts, the latter probably occurring during major mergers or in dense nuclear SF regions. Both modes are observable over a large range of SF rates. The detection of CO emission from distant near-IR selected galaxies reveals such bimodal behavior for the first time, as they allow us to probe gas in disk galaxies with much higher SF rates than are seen locally. The different regimes can potentially be interpreted as the effect of a top-heavy IMF in starbursts. However, we favor a different physical origin related to the fraction of molecular gas in dense clouds. The IR luminosity to gas mass ratio (i.e., the SF efficiency) appears to be inversely proportional to the dynamical (rotation) timescale. Only when accounting for the dynamical timescale, a universal SF law is obtained, suggesting a direct link between global galaxy properties and the local SF rate., 5 pages, 4 figures. ApJ Letters in press

[abridged] The properties of galaxies with the reddest observed R-K colors (Extremely Red Objects, EROs), including their apparent division into passive and obscured active objects with roughly similar number densities, are a known challenge for models of galaxy formation. We produce mock catalogues generated by interfacing the predictions of the semi-analytical MORGANA model for the evolution of galaxies in a Lambda-CDM cosmology with the spectro-photometric + radiative transfer code GRASIL and Infrared (IR) template library to show that the model correctly reproduces number counts, redshift distributions and active fractions of R-K>5 sources. We test the robustness of our results against different dust attenuations and, most importantly, against the inclusion of TP-AGB stars in Simple Stellar Populations used to generate galaxy spectra, and find that the inclusion of TP-AGBs has a relevant effect, in that it allows to increase by a large factor the number of very red active objects at all color cuts. We find that though the most passive and the most obscured active galaxies have a higher probability of being selected as EROs, many EROs have intermediate properties and the population does not show bimodality in specific star formation rate (SSFR). We predict that deep observations in the Far-IR, from 100 to 500 micron, are the most efficient way to constrain the SSFR of these objects; we give predictions for future Herschel observations. Finally, we test whether a simple evolutionary sequence for the formation of z=0 massive galaxies, going through a sub-mm-bright phase and then a ERO phase, are typical in this galaxy formation model. We find that this sequence holds for ~25 per cent of z=0 massive galaxies, while the model typically shows a more complex connection between sub-mm, ERO and massive galaxies. [abridged], 14 pages, 12 figures, MNRAS in press

By comparing galaxy formation models within the ΛCDM cosmology with available data on high redshift galaxies, we show that small galaxies , with stellar masses

[abridged] It has been widely claimed that several lines of observational evidence point towards a "downsizing" (DS) of the process of galaxy formation over cosmic time. This behavior is sometimes termed "anti-hierarchical", and contrasted with the "bottom-up" assembly of the dark matter structures in Cold Dark Matter models. In this paper we address three different kinds of observational evidence that have been described as DS: the stellar mass assembly, star formation rate and the ages of the stellar populations in local galaxies. We compare a broad compilation of available data-sets with the predictions of three different semi-analytic models of galaxy formation within the Lambda-CDM framework. In the data, we see only weak evidence at best of DS in stellar mass and in star formation rate. We find that, when observational errors on stellar mass and SFR are taken into account, the models acceptably reproduce the evolution of massive galaxies, over the entire redshift range that we consider. However, lower mass galaxies are formed too early in the models and are too passive at late times. Thus, the models do not correctly reproduce the DS trend in stellar mass or the archaeological DS, while they qualitatively reproduce the mass-dependent evolution of the SFR. We demonstrate that these discrepancies are not solely due to a poor treatment of satellite galaxies but are mainly connected to the excessively efficient formation of central galaxies in high-redshift haloes with circular velocities ~100-200 km/s. [abridged], Comment: MNRAS accepted, 16 pages, 10 figures

The high end of the stellar mass function of galaxies is observed to have little evolution since z~1. This represents a stringent constraint for merger--based models, aimed at explaining the evolution of the most massive galaxies in the concordance LambdaCDM cosmology. In this Letter we show that it is possible to remove the tension between the above observations and model predictions by allowing a fraction of stars to be scattered to the Diffuse Stellar Component (DSC) of galaxy clusters at each galaxy merger, as recently suggested by the analysis of N-body hydrodynamical simulations. To this purpose, we use the MORGANA model of galaxy formation in a minimal version, in which gas cooling and star formation are switched off after z=1. In this way, any predicted evolution of the galaxy stellar mass function is purely driven by mergers. We show that, even in this extreme case, the predicted degree of evolution of the high end of the stellar mass function is larger than that suggested by data. Assuming instead that a significant fraction, ~30 per cent, of stars are scattered in the DSC at each merger event, leads to a significant suppression of the predicted evolution, in better agreement with observational constraints, while providing a total amount of DSC in clusters which is consistent with recent observational determinations., Comment: 5 pages, figures included; ApJ Letters, in press. Revision: reference added