Search

We present cosmological zoom-in hydrodynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain-size distribution, with the MUPPI star formation and feedback subresolution model. Our dust evolution model follows carbon and silicate dust separately. To distinguish differences induced by the chaotic behaviour of simulations from those genuinely due to different simulation set-up, we run each model six times, after introducing tiny perturbations in the initial conditions. With this method, we discuss the role of various dust-related physical processes and the effect of a few possible approximations adopted in the literature. Metal depletion and dust cooling affect the evolution of the system, causing substantial variations in its stellar, gas, and dust content. We discuss possible effects on the Spectral Energy Distribution of the significant variations of the size distribution and chemical composition of grains, as predicted by our simulations during the evolution of the galaxy. We compare dust surface density, dust-to-gas ratio, and small-to-large grain mass ratio as a function of galaxy radius and gas metallicity predicted by our fiducial run with recent observational estimates for three disc galaxies of different masses. The general agreement is good, in particular taking into account that we have not adjusted our model for this purpose. [ABSTRACT FROM AUTHOR]

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 on to 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 on to 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 (AGNs) and we confirm that a non-instantaneous accretion time-scale (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 and from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies (>1011 M⊙) 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 <1010 M⊙. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than 108 M⊙ are already in a self-regulated state as inferred from observations. [ABSTRACT FROM AUTHOR]

We present simulations of galaxy formation, based on the gadget -3 code, in which a sub-resolution model for star formation and stellar feedback is interfaced with a new model for active galactic nucleus (AGN) feedback. Our sub-resolution model describes a multiphase interstellar medium (ISM), accounting for hot and cold gas within the same resolution element: we exploit this feature to investigate the impact of coupling AGN feedback energy to the different phases of the ISM over cosmic time. Our fiducial model considers that AGN feedback energy coupling is driven by the covering factors of the hot and cold phases. We perform a suite of cosmological hydrodynamical simulations of disc galaxies (⁠|$M_{\rm halo,DM} \simeq 2 \times 10^{12}$|  M⊙, at z  = 0), to investigate: (i) the effect of different ways of coupling AGN feedback energy to the multiphase ISM; (ii) the impact of different prescriptions for gas accretion (i.e. only cold gas, both cold and hot gas, with the additional possibility of limiting gas accretion from cold gas with high angular momentum); (iii) how different models of gas accretion and coupling of AGN feedback energy affect the coevolution of supermassive black holes (BHs) and their host galaxy. We find that at least a share of the AGN feedback energy has to couple with the diffuse gas, in order to avoid an excessive growth of the BH mass. When the BH only accretes cold gas, it experiences a growth that is faster than in the case in which both cold and hot gas are accreted. If the accretion of cold gas with high angular momentum is reduced, the BH mass growth is delayed, the BH mass at z  = 0 is reduced by up to an order of magnitude, and the BH is prevented from accreting below z ≲ 2, when the galaxy disc forms. [ABSTRACT FROM AUTHOR]

We perform a suite of cosmological hydrodynamical simulations of disc galaxies, with zoomed-in initial conditions leading to the formation of a halo of mass |$M_{\rm halo, \, DM} \simeq 2 \times 10^{12}$|  M⊙ at redshift z  = 0. These simulations aim at investigating the chemical evolution and the distribution of metals in a disc galaxy, and at quantifying the effect of (i) the assumed IMF, (ii) the adopted stellar yields, and (iii) the impact of binary systems originating SNe Ia on the process of chemical enrichment. We consider either a Kroupa, Tout & Gilmore (1993) or a more top-heavy Kroupa (2001) IMF, two sets of stellar yields and different values for the fraction of binary systems suitable to give rise to SNe Ia. We investigate stellar ages, SN rates, stellar and gas metallicity gradients, and stellar α-enhancement in simulations, and compare predictions with observations. We find that a Kroupa et al. (1993) IMF has to be preferred when modelling late-type galaxies in the local Universe. On the other hand, the comparison of stellar metallicity profiles and α-enhancement trends with observations of Milky Way stars shows a better agreement when a Kroupa (2001) IMF is assumed. Comparing the predicted SN rates and stellar α-enhancement with observations supports a value for the fraction of binary systems producing SNe Ia of 0.03, at least for late-type galaxies and for the considered IMFs. Adopted stellar yields are crucial in regulating cooling and star formation, and in determining patterns of chemical enrichment for stars, especially for those located in the galaxy bulge. [ABSTRACT FROM AUTHOR]

We investigate the impact of galactic outflowmodelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions (ICs) of a Milky Way-sized halo. We verify how sensitive the general properties of the simulated galaxy are to the way in which stellar feedback triggered outflows are implemented, keeping ICs, simulation code and star formation (SF) model all fixed. We present simulations that are based on a version of the GADGET3 code where our sub-resolution model is coupled with an advanced implementation of smoothed particle hydrodynamics that ensures a more accurate fluid sampling and an improved description of gas mixing and hydrodynamical instabilities. We quantify the strong interplay between the adopted hydrodynamic scheme and the sub-resolution model describing SF and feedback. We consider four different galactic outflow models, including the one introduced by Dalla Vecchia & Schaye (2012) and a scheme that is inspired by the Springel & Hernquist (2003) model. We find that the sub-resolution prescriptions adopted to generate galactic outflows are the main shaping factor of the stellar disc component at low redshift. The key requirement that a feedback model must have to be successful in producing a disc-dominated galaxy is the ability to regulate the high-redshift SF (responsible for the formation of the bulge component), the cosmological infall of gas from the large-scale environment, and gas fall-back within the galactic radius at low redshift, in order to avoid a too high SF rate at z = 0. [ABSTRACT FROM AUTHOR]