No Arabic abstract
This is the first paper of a series that describes the methods and basic results of the GalICS model (for Galaxies In Cosmological Simulations). GalICS is a hybrid model for hierarchical galaxy formation studies, combining the outputs of large cosmological N-body simulations with simple, semi-analytic recipes to describe the fate of the baryons within dark matter halos. The simulations produce a detailed merging tree for the dark matter halos including complete knowledge of the statistical properties arising from the gravitational forces. We intend to predict the overall statistical properties of galaxies, with special emphasis on the panchromatic spectral energy distribution emitted by galaxies in the UV/optical and IR/submm wavelength ranges. In this paper, we outline the physically motivated assumptions and key free parameters that go into the model, comparing and contrasting with other parallel efforts. We specifically illustrate the success of the model in comparison to several datasets, showing how it is able to predict the galaxy disc sizes, colours, luminosity functions from the ultraviolet to far infrared, the Tully--Fisher and Faber--Jackson relations, and the fundamental plane in the local universe. We also identify certain areas where the model fails, or where the assumptions needed to succeed are at odds with observations, and pay special attention to understanding the effects of the finite resolution of the simulations on the predictions made. Other papers in this series will take advantage of different data sets available in the literature to extend the study of the limitations and predictive power of GalICS, with particular emphasis put on high-redshift galaxies.
We construct the Numerical Galaxy Catalog ($ u$GC), based on a semi-analytic model of galaxy formation combined with high-resolution N-body simulations in a $Lambda$-dominated flat cold dark matter ($Lambda$CDM) cosmological model. The model includes several essential ingredients for galaxy formation, such as merging histories of dark halos directly taken from N-body simulations, radiative gas cooling, star formation, heating by supernova explosions (supernova feedback), mergers of galaxies, population synthesis, and extinction by internal dust and intervening HI clouds. As the first paper in a series using this model, we focus on basic photometric, structural and kinematical properties of galaxies at present and high redshifts. Two sets of model parameters are examined, strong and weak supernova feedback models, which are in good agreement with observational luminosity functions of local galaxies in a range of observational uncertainty. Both models agree well with many observations such as cold gas mass-to-stellar luminosity ratios of spiral galaxies, HI mass functions, galaxy sizes, faint galaxy number counts and photometric redshift distributions in optical pass-bands, isophotal angular sizes, and cosmic star formation rates. In particular, the strong supernova feedback model is in much better agreement with near-infrared (K-band) faint galaxy number counts and redshift distribution than the weak feedback model and our previous semi-analytic models based on the extended Press-Schechter formalism. (Abridged)
High-resolution N-body re-simulations of 15 massive (10^{14}-10^{15} Msun) dark matter haloes have been combined with the hybrid galaxy formation model GalICS (Hatton et al. 2003), to study the formation and evolution of galaxies in clusters, within the framework of the hierarchical merging scenario. New features in GalICS include a better description of galaxy positioning within dark matter haloes, a more reliable computation of the temperature of the inter-galactic medium as a function of redshift, and a description of the ram pressure stripping process. We focus on the luminosity functions, morphological fractions and colour distributions of galaxies in clusters and in cluster outskirts, at z=0. No systematic dependency on cluster richness is found either for the galaxy luminosity functions, morphological mixes, or colour distributions. Moving from higher density (cluster cores), to lower density environments (cluster outskirts), we detect a progressive flattening of the luminosity functions, an increase of the fraction of spirals and a decrease of that of ellipticals and S0s, and the progressive emergence of a bluer tail in the distributions of galaxy colours, especially for spirals. As compared to cluster spirals, early-type galaxies show a flatter luminosity function, and more homogeneous and redder colours. An overall good agreement is found between our results and the observations, particularly in terms of the cluster luminosity functions and morphological mixes. However, some discrepancies are also apparent, with too faint magnitudes of the brightest cluster members, especially in the B band, and galaxy colours tendentially too red (or not blue enough) in the model, with respect to the observations. Finally, ram pressure stripping appears to affect very little our results.
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 generate mock galaxy catalogues for a grid of different cosmologies, using rescaled N-body simulations in tandem with a semi-analytic model run using consistent parameters. Because we predict the galaxy bias, rather than fitting it as a nuisance parameter, we obtain an almost pure constraint on sigma_8 by comparing the projected two-point correlation function we obtain to that from the SDSS. A systematic error arises because different semi-analytic modelling assumptions allow us to fit the r-band luminosity function equally well. Combining our estimate of the error from this source with the statistical error, we find sigma_8=0.97 +/- 0.06. We obtain consistent results if we use galaxy samples with a different magnitude threshold, or if we select galaxies by b_J-band rather than r-band luminosity and compare to data from the 2dFGRS. Our estimate for sigma_8 is higher than that obtained for other analyses of galaxy data alone, and we attempt to find the source of this difference. We note that in any case, galaxy clustering data provide a very stringent constraint on galaxy formation models.