We present cosmological zoom-in hydro-dynamical 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 sub-resolution 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.

No abstract

We investigate the impact of galactic outflow modelling 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.

We study the impact of relatively weak AGN feedback on the interstellar medium of intermediate and massive elliptical galaxies. We find that the AGN activity, while globally heating the ISM, naturally stimulates some degree of hot gas cooling on scales of several kpc. This process generates the persistent presence of a cold ISM phase, with mass ranging between 10 4 and 5 × 10 7 M , where the latter value is appropriate for group centered, massive galaxies. Widespread cooling occurs where the ratio of cooling to free-fall time before the activation of the AGN feedback satisfies t cool /t ff 70, that is we find a less restrictive threshold than commonly quoted in the literature. This process helps explaining the body of observations of cold gas (both ionized and neutral/molecular) in Ellipticals and, perhaps, the residual star formation detected in many early-type galaxies. The amount and distribution of the off-center cold gas vary irregularly with time.The cold ISM velocity field is irregular, initially sharing the (outflowing) turbulent hot gas motion. Typical velocity dispersions of the cold gas lie in the range 100 − 200 km s −1 . Freshly generated cold gas often forms a cold outflow and can appear kinematically misaligned with respect to the stars.We also follow the dust evolution in the hot and cold gas. We find that the internally generated cold ISM has a very low dust content, with representative values of the dust-to-gas ratio of 10 −4 − 10 −5 . Therefore, this cold gas can escape detection in the traditional dust-absorption maps.

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 AGN feedback. Our sub-resolution model describes a multiphase 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 halo, DM 2 · 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 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.

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 halo, DM 2 · 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 et al. (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.

We use 47 gravitational wave sources from the Third LIGO–Virgo–Kamioka Gravitational Wave Detector Gravitational Wave Transient Catalog (GWTC–3) to estimate the Hubble parameter H(z), including its current value, the Hubble constant H 0. Each gravitational wave (GW) signal provides the luminosity distance to the source, and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and H(z). The source mass distribution displays a peak around 34 M ⊙, followed by a drop-off. Assuming this mass scale does not evolve with the redshift results in a H(z) measurement, yielding H 0 = 68 − 8 + 12 km s − 1 Mpc − 1 (68% credible interval) when combined with the H 0 measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the H 0 estimate from GWTC–1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+, statistically marginalizing over the redshifts of each event’s potential hosts. Assuming a fixed BBH population, we estimate a value of H 0 = 68 − 6 + 8 km s − 1 Mpc − 1 with the galaxy catalog method, an improvement of 42% with respect to our GWTC–1 result and 20% with respect to recent H 0 studies using GWTC–2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about H 0) is the well-localized event GW190814.

The properties of quasar-host galaxies might be determined by the growth and feedback of their supermassive (SMBH, 108 − 10 M⊙) black holes. We investigate such connection with a suite of cosmological simulations of massive (halo mass ≈1012 M⊙) galaxies at z ≃ 6 which include a detailed sub-grid multiphase gas and accretion model. BH seeds of initial mass 105 M⊙ grow mostly by gas accretion, and become SMBH by z = 6 setting on the observed MBH − M⋆ relation without the need for a boost factor. Although quasar feedback crucially controls the SMBH growth, its impact on the properties of the host galaxy at z = 6 is negligible. In our model, quasar activity can both quench (via gas heating) or enhance (by ISM over-pressurization) star formation. However, we find that the star formation history is insensitive to such modulation as it is largely dominated, at least at z > 6, by cold gas accretion from the environment that cannot be hindered by the quasar energy deposition. Although quasar-driven outflows can achieve velocities $> 1000~\rm km~s^{-1}$, only ≈4 per cent of the outflowing gas mass can actually escape from the host galaxy. These findings are only loosely constrained by available data, but can guide observational campaigns searching for signatures of quasar feedback in early galaxies.