No Arabic abstract
Massive early-type galaxies have higher metallicities and higher ratios of $alpha$ elements to iron than their less massive counterparts. Reproducing these correlations has long been a problem for hierarchical galaxy formation theory, both in semi-analytic models and cosmological hydrodynamic simulations. We show that a simulation in which gas cooling in massive dark haloes is quenched by radio-mode active galactic nuclei (AGNs) feedback naturally reproduces the observed trend between $alpha$/Fe and the velocity dispersion of galaxies, $sigma$. The quenching occurs earlier for more massive galaxies. Consequently, these galaxies complete their star formation before $alpha$/Fe is diluted by the contribution from type Ia supernovae. For galaxies more massive than $sim 10^{11}~M_odot$ whose $alpha$/Fe correlates positively with stellar mass, we find an inversely correlated mass-metallicity relation. This is a common problem in simulations in which star formation in massive galaxies is quenched either by quasar- or radio-mode AGN feedback. The early suppression of gas cooling in progenitors of massive galaxies prevents them from recapturing enriched gas ejected as winds. Simultaneously reproducing the [$alpha$/Fe]-$sigma$ relation and the mass-metallicity relation is, thus, difficult in the current framework of galaxy formation.
We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.
Over the last decades, cosmological simulations of galaxy formation have been instrumental for advancing our understanding of structure and galaxy formation in the Universe. These simulations follow the non-linear evolution of galaxies modeling a variety of physical processes over an enormous range of scales. A better understanding of the physics relevant for shaping galaxies, improved numerical methods, and increased computing power have led to simulations that can reproduce a large number of observed galaxy properties. Modern simulations model dark matter, dark energy, and ordinary matter in an expanding space-time starting from well-defined initial conditions. The modeling of ordinary matter is most challenging due to the large array of physical processes affecting this matter component. Cosmological simulations have also proven useful to study alternative cosmological models and their impact on the galaxy population. This review presents a concise overview of the methodology of cosmological simulations of galaxy formation and their different applications.
Observations have established that the diffuse intergalactic medium (IGM) at z ~ 3 is enriched to ~0.1-1% solar metallicity and that the hot gas in large clusters of galaxies (ICM) is enriched to 1/3-1/2 solar metallicity at z=0. Metals in the IGM may have been removed from galaxies (in which they presumably form) during dynamical encounters between galaxies, by ram-pressure stripping, by supernova-driven winds, or as radiation-pressure driven dust efflux. This study develops a method of investigating the chemical enrichment of the IGM and of galaxies, using already completed cosmological simulations. To these simulations, we add dust and (gaseous) metals, distributing the dust and metals in the gas according to three simple parameterized prescriptions, one for each enrichment mechanism. These prescriptions are formulated to capture the basic ejection physics, and calibrated when possible with empirical data. Our results indicate that dynamical removal of metals from >~ 3*10^8 solar mass galaxies cannot account for the observed metallicity of low-column density Ly-alpha absorbers, and that dynamical removal from >~ 3*10^10 solar mass galaxies cannot account for the ICM metallicities. Dynamical removal also fails to produce a strong enough mass-metallicity relation in galaxies. In contrast, either wind or radiation-pressure ejection of metals from relatively large galaxies can plausibly account for all three sets of observations (though it is unclear whether metals can be distributed uniformly enough in the low-density regions without overly disturbing the IGM, and whether clusters can be enriched quite as much as observed). We investigate in detail how our results change with variations in our assumed parameters, and how results for the different ejection processes compare. (Abridged)
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-big 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.
The fate of metals ejected by young OB associations into the Interstellar Medium (ISM) is investigated numerically. In particular, we study the enrichment of the cold gas phase, which is the material that forms molecular clouds. Following previous work, the expansion and collision of two supershells in a diffuse ISM is simulated, in this case also introducing an advected quantity which represents the metals expelled by the young stars. We adopt the simplest possible approach, not differentiating between metals coming from stellar winds and those coming from supernovae. Even though the hot, diffuse phase of the ISM receives a significant amount of metals from the stars, the cold phase is efficiently shielded, with very little metal enrichment. Significant enrichment of the cold ISM will therefore be delayed by at least the cooling time of this hot phase. No variations in cloud metallicity with distance from the OB association or with direction are found, which means that the shell collision does little to enhance the metallicity of the cold clumps. We conclude that the stellar generation that forms out of molecular structures, triggered by shell collisions cannot be significantly enriched.