No Arabic abstract
We discuss a model for treating chemical enrichment by SNII and SNIa explosions in simulations of cosmological structure formation. Our model includes metal-dependent radiative cooling and star formation in dense collapsed gas clumps. Metals are returned into the diffuse interstellar medium by star particles using a local SPH smoothing kernel. A variety of chemical abundance patterns in enriched gas arise in our treatment owing to the different yields and lifetimes of SNII and SNIa progenitor stars. In the case of SNII chemical production, we adopt metal-dependent yields. Because of the sensitive dependence of cooling rates on metallicity, enrichment of galactic haloes with metals can in principle significantly alter subsequent gas infall and the build up of the stellar components. Indeed, in simulations of isolated galaxies we find that a consistent treatment of metal-dependent cooling produces 25% more stars outside the central region than simulations with a primordial cooling function. In the highly-enriched central regions, the evolution of baryons is however not affected by metal cooling, because here the gas is always dense enough to cool. A similar situation is found in cosmological simulations because we include no strong feedback processes which could spread metals over large distances and mix them into unenriched diffuse gas. We demonstrate this explicitly with test simulations which adopt super-solar cooling functions leading to large changes both in the stellar mass and in the metal distributions. We also find that the impact of metallicity on the star formation histories of galaxies may depend on their particular evolutionary history. Our results hence emphasise the importance of feedback processes for interpreting the cosmic metal enrichment.
We have developed a new scheme to treat a multiphase interstellar medium in smoothed particle hydrodynamics simulations of galaxy formation. This scheme can represent a co-spatial mixture of cold and hot ISM components, and is formulated without scale-dependent parameters. It is thus particularly suited to studies of cosmological structure formation where galaxies with a wide range of masses form simultaneously. We also present new algorithms for energy and heavy element injection by supernovae, and show that together these schemes can reproduce several important observed effects in galaxy evolution. Both in collapsing systems and in quiescent galaxies our codes can reproduce the Kennicutt relation between the surface densities of gas and of star formation. Strongly metal-enhanced winds are generated in both cases with ratios of mass-loss to star formation which are similar to those observed. This leads to a self-regulated cycle for star formation activity. The overall impact of feedback depends on galaxy mass. Star formation is suppressed at most by a factor of a few in massive galaxies, but in low-mass systems the effects can be much larger, giving star formation an episodic, bursty character. The larger the energy fraction assumed available in feedback, the more massive the outflows and the lower the final stellar masses. Winds from forming disks are collimated perpendicular to the disk plane, reach velocities up to 1000 km/s, and efficiently transport metals out of the galaxies. The asymptotically unbound baryon fraction drops from >95 per cent to ~30 per cent from the least to the most massive of our idealised galaxies, but the fraction of all metals ejected with this component exceeds 60 per cent regardless of mass. Such winds could plausibly enrich the intergalactic medium to observed levels.
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)
A stochastic model of the chemical enrichment of metal-poor systems by core-collapse (Type II) supernovae is presented, allowing for large-scale mixing of the enriched material by turbulent motions and cloud collisions in the interstellar medium. Infall of pristine material is taken into account by following the evolution of the gas density in the medium. Analytical expressions were derived for the number of stars enriched by a given number of supernovae, as well as for the amount of mass with which the ejected material from a supernova is mixed before being locked up in a subsequently formed star. It is shown that for reasonable values of the gas density (~0.1 cm-3) and of the supernova rate (~0.25 kpc-3 Myr-1) of the Galactic halo, the resulting metallicity distributions of the extreme Population II stars show a distinct cut-off at [Fe/H] ~= -4. In fact, by assuming no low-mass Population III stars were able to form out of the primordial interstellar medium, the derived fraction of stars below [Fe/H] = -4 is in agreement with observations. Moreover, the probability is high that even the most metal-poor stars observed to date have been enriched by several contributing supernovae. This partly explains the relatively small star-to-star scatter in many chemical-abundance ratios for stars down to [Fe/H] = -4, as recently found in several observational studies. Contribution from the thermonuclear (Type Ia) supernovae is found to be negligible over almost the entire extremely metal-poor regime. (***abridged***)
The distribution of metals in the intracluster medium (ICM) of galaxy clusters provides valuable information on their formation and evolution, on the connection with the cosmic star formation and on the effects of different gas processes. By analyzing a sample of simulated galaxy clusters, we study the chemical enrichment of the ICM, its evolution, and its relation with the physical processes included in the simulation and with the thermal properties of the core. These simulations, consisting of re-simulations of 29 Lagrangian regions performed with an upgraded version of the SPH GADGET-3 code, have been run including two different sets of baryonic physics: one accounts for radiative cooling, star formation, metal enrichment and supernova (SN) feedback, and the other one further includes the effects of feedback from active galactic nuclei (AGN). In agreement with observations, we find an anti-correlation between entropy and metallicity in cluster cores, and similar radial distributions of heavy-element abundances and abundance ratios out to large cluster-centric distances (~R180). In the outskirts, namely outside of ~0.2R180, we find a remarkably homogeneous metallicity distribution, with almost flat profiles of the elements produced by either SNIa or SNII. We investigated the origin of this phenomenon and discovered that it is due to the widespread displacement of metal-rich gas by early (z>2-3) AGN powerful bursts, acting on small high-redshift haloes. Our results also indicate that the intrinsic metallicity of the hot gas for this sample is on average consistent with no evolution between z=2 and z=0, across the entire radial range.
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.