No Arabic abstract
We implement a detailed dust model into the L-Galaxies semi-analytical model which includes: injection of dust by type II and type Ia supernovae (SNe) and AGB stars; grain growth in molecular clouds; and destruction due to supernova-induced shocks, star formation, and reheating. Our grain growth model follows the dust content in molecular clouds and the inter-cloud medium separately, and allows growth only on pre-existing dust grains. At early times, this can make a significant difference to the dust growth rate. Above $zsim8$, type II SNe are the primary source of dust, whereas below $zsim8$, grain growth in molecular clouds dominates, with the total dust content being dominated by the latter below $zsim6$. However, the detailed history of galaxy formation is important for determining the dust content of any individual galaxy. We introduce a fit to the dust-to-metal (DTM) ratio as a function of metallicity and age, which can be used to deduce the DTM ratio of galaxies at any redshift. At $zlesssim3$, we find a fairly flat mean relation between metallicity and the DTM, and a positive correlation between metallicity and the dust-to-gas (DTG) ratio, in good agreement with the shape and normalisation of the observed relations. We also match the normalisation of the observed stellar mass -- dust mass relation over the redshift range of $0-4$, and to the dust mass function at $z=0$. Our results are important in interpreting observations on the dust content of galaxies across cosmic time, particularly so at high redshift.
We investigate the dynamical evolution of galaxies in groups with different formation epochs. Galaxy groups have been selected to be in different dynamical states, namely dynamically old and dynamically young, which reflect their early and late formation times, respectively, based on their halo mass assembly. Brightest galaxies in dynamically young groups have suffered their last major galaxy merger typically $sim 2$ Gyr more recently than their counterparts in dynamically old groups. Furthermore, we study the evolution of velocity dispersion in these two classes and compare them with the analytic models of isolated halos. The velocity dispersion of dwarf galaxies in high mass, dynamically young groups increases slowly in time, while the analogous dispersion in dynamically old high-mass groups is constant. In contrast, the velocity dispersion of giant galaxies in low mass groups decreases rapidly at late times. This increasing velocity bias is caused by dynamical friction, and starts much earlier in the dynamically old groups. The recent {sc Radio-SAGE} model of galaxy formation suggests that radio luminosities of central galaxies, considered to be tracers of AGN activity, are enhanced in halos that assembled more recently, independent of the time since the last major merger.
We introduce a new physical recipe into the De Lucia and Blaizot version of the Munich semi-analytic model built upon the Millennium dark matter simulation: the tidal stripping of stellar material from satellite galaxies during mergers. To test the significance of the new physical process we apply a Monte Carlo Markov Chain parameter estimation technique constraining the model with the $K$-band luminosity function, $B-V$ colours and the black hole-bulge mass relation. The differences in parameter correlations, and in the allowed regions in likelihood space, reveal the impact of the new physics on the basic ingredients of the model, such as the star-formation laws, feedback recipes and the black hole growth model. With satellite disruption in place, we get a model likelihood four times higher than in the original model, indicating that the new process seems to be favoured by observations. This is achieved mainly due to a reduction in black hole growth that produces a better agreement between the properties of central black holes and host galaxies. Compared to the best-fit model without disruption, the new model removes the excess of dwarf galaxies in the original recipe with a more modest supernova heating. The new model is now consistent with the three observational data sets used to constrain it, while significantly improving the agreement with observations for the distribution of metals in stars. Moreover, the model now follows the build up of intra-cluster light.
(Abridged) This paper predicts self-consistent faint galaxy counts from the UV to the submm wavelength range. The STARDUST spectral energy distributions described in Devriendt et al. (1999) are embedded within the explicit cosmological framework of a simple semi-analytic model of galaxy formation and evolution. We build a class of models which capture the luminosity budget of the universe through faint galaxy counts and redshift distributions in the whole wavelength range spanned by our spectra. In contrast with a rather stable behaviour in the optical and even in the far-IR, the submm counts are dramatically sensitive to variations in the cosmological parameters and changes in the star formation history. Faint submm counts are more easily accommodated within an open universe with a low value of $Omega_0$, or a flat universe with a non-zero cosmological constant. This study illustrates the implementation of multi-wavelength spectra into a semi-analytic model. In spite of its simplicity, it already provides fair fits of the current data of faint counts, and a physically motivated way of interpolating and extrapolating these data to other wavelengths and fainter flux levels.
We present a variation of the recently updated Munich semi-analytical galaxy formation model, L-Galaxies, with a new gas stripping method. Extending earlier work, we directly measure the local environmental properties of galaxies to formulate a more accurate treatment of ram-pressure stripping for all galaxies. We fully re-calibrate the modified L-Galaxies model using a Markov Chain Monte Carlo (MCMC) method with the stellar mass function and quenched fraction of galaxies at $0leq zleq2$ as constraints. Due to this re-calibration, global galaxy population relations, including the stellar mass function, quenched fractions versus galaxy mass and HI mass function are all largely unchanged and remain consistent with observations. By comparing to data on galaxy properties in different environments from the SDSS and HSC surveys, we demonstrate that our modified model improves the agreement with the quenched fractions and star formation rates of galaxies as a function of environment, stellar mass, and redshift. Overall, in the vicinity of haloes with total mass $10^{12}$ to $10^{15}rm M_{odot}$ at $z=0$, our new model produces higher quenched fractions and stronger environmental dependencies, better recovering observed trends with halocentric distance up to several virial radii. By analysing the actual amount of gas stripped from galaxies in our model, we show that those in the vicinity of massive haloes lose a large fraction of their hot halo gas before they become satellites. We demonstrate that this affects galaxy quenching both within and beyond the halo boundary. This is likely to influence the correlations between galaxies up to tens of megaparsecs.
It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.