We present the Mock Map Facility, a powerful tool to generate mock catalogues or images from semi-analytically post-processed snapshots of cosmological N-body simulations. The paper describes in detail an efficient technique to create such mocks from the GALICS semi-analytic model, providing the reader with an accurate quantification of the artifacts it introduces at every step. We show that replication effects introduce a negative bias on the clustering signal -- typically peaking at less than 10 percent around the correlation length. We also thoroughly discuss how the clustering signal is affected by finite volume effects, and show that it vanishes at scales larger than about a tenth of the simulation box size. For the purpose of analysing our method, we show that number counts and redshift distributions obtained with GALICS and MOMAF compare well to K-band observations and to the 2dFGRS. Given finite volume effects, we also show that the model can reproduce the APM angular correlation function. The MOMAF results discussed here are made publicly available to the astronomical community through a public database. Moreover, a user-friendly Web interface (http://galics.iap.fr) allows any user to recover her/his own favourite galaxy samples through simple SQL queries. The flexibility of this tool should permit a variety of uses ranging from extensive comparisons between real observations and those predicted by hierarchical models of galaxy formation, to the preparation of observing strategies for deep surveys and tests of data processing pipelines.

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.

This paper illustrates how mock observational samples of high-redshift galaxies with sophisticated selection criteria can be extracted from the predictions of GALICS, a hybrid model of hierarchical galaxy formation that couples the outputs of large cosmological simulations and semi-analytic recipes to describe dark matter collapse and the physics of baryons respectively. As an example of this method, we focus on the properties of Lyman Break Galaxies at redshift 3. With the MOMAF software package described in a companion paper, we generate a mock observational sample with selection criteria as similar as possible to those implied in the actual observations of z = 3 LBGs by Steidel et al.(1995). Our model predictions are in good agreement with the observed number density and 2D correlation function. We investigate the optical/IR luminosity budget as well as several other physical properties of LBGs and find them to be in general agreement with observed values. Looking into the future of these LBGs we predict that 75% of them end up as massive ellipticals today, even though only 35% of all our local ellipticals are predicted to have a LBG progenitor. In spite of some shortcomings, this new 'mock observation' method clearly represents a necessary first step toward a more accurate comparison between hierarchical models of galaxy formation and real observational surveys.

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.

The next generation of submillimetre/millimetre instruments will provide us with a deeper insight into the mechanisms that rule galaxy formation. As the brightest starbursts are thought to be heavily obscured at optical wavelengths, the opening of this new window will complement the present observations, and enable a detailed investigation of the hierarchical merging of galaxies at remote epochs. In this context, we will present preliminary results of simulations currently developed for deep galaxy surveys with FIRST. N-body simulations of galaxy mergers are being developed to produce realistic morphologies and star formation histories. A Schmidt law is used to account for the evolution of the stellar and gas content of each particle. These simulations will be gathered in a library, including the temporal information with a timestep of 10 Myr. We then build consistent spectra accounting for this star formation history at each resolution element of the simulations. These morphologies will be included in the framework of the GalICS hybrid model of hierarchical galaxy formation, which reproduces the main observational constraints. We intend to present synthetic maps with the characteristics of the FIRST PACS/SPIRE instruments and discuss the optimal strategy for deep surveys complemented with the ground-based ALMA project.

Given the failure of existing models for redshift-space distortions to provide a highly accurate measure of the beta-parameter, and the ability of forthcoming surveys to obtain data with very low random errors, it becomes necessary to develop better models for these distortions. Here we review the failures of the commonly-used velocity dispersion models and present an empirical method for extracting beta from the quadrupole statistic that has little systematic offset over a wide range of beta and cosmologies. This empirical model is then applied to an ensemble of mock 2dF southern strip surveys, to illustrate the technique and see how accurately we can recover the true value of beta. We compare this treatment with the error we expect to find due only to the finite volume of the survey. We find that non-linear effects reduce the range of scales over which beta can be fitted, and introduce covariances between nearby modes in addition to those introduced by the convolution with the survey window function. The result is that we are only able to constrain beta to a 1-sigma accuracy of 25% (beta=0.55 +/- 0.14 for the cosmological model considered). We explore one possible means of reducing this error, that of cluster collapse, and show that accurate application of this method can greatly reduce the effect of non-linearities, improving the determination of beta. We conclude by demonstrating that, when the contaminating effects of clusters are dealt with, this simple analysis of the full 2dF survey yields beta=0.55 +/- 0.04. For this model this represents a determination of beta to an accuracy of 8% and hence an important constraint on the cosmological density parameter Omega_0.

Recent results from a number of redshift surveys suggest that the Universe is well described by an inhomogeneous, fractal distribution on the largest scales probed. This distribution has been found to have fractal dimension, D, approximately equal to 2.1, in contrast to a homogeneous distribution in which the dimension should approach the value 3 as the scale is increased. In this paper we demonstrate that estimates of D, based on the conditional density of galaxies, are prone to bias from several sources. These biases generally result in a smaller measured fractal dimension than the true dimension of the sample. We illustrate this behaviour in application to the Stromlo-APM redshift survey, showing that this dataset in fact provides evidence for fractal dimension increasing with survey depth. On the largest scale probed, r approximately 60 h^{-1}Mpc, we find evidence for a distribution with dimension D = 2.76 +/- 0.10. A comparison between this sample and mock Stromlo-APM catalogues taken from Nbody simulations (which assume a CDM cosmology) reveals a striking similarity in the behaviour of the fractal dimension. Thus we find no evidence for inhomogeneity in excess of that expected from conventional cosmological theory. We consider biases affecting future large surveys and demonstrate, using mock SDSS catalogues, that this survey will be able to measure the fractal dimension on scales at which we expect to see full turn-over to homogeneity, in an accurate and unbiased way.

We present a comprehensive set of mock 2dF and SDSS galaxy redshift surveys constructed from a set of large, high-resolution cosmological N-body simulations. The radial selection functions and geometrical limits of the catalogues mimic those of the genuine surveys. The catalogues span a wide range of cosmologies, including both open and flat universes. In all the models the galaxy distributions are biased so as to approximately reproduce the observed galaxy correlation function on scales of 1--10 Mpc/h In some cases models with a variety of different biasing prescriptions are included. All the mock catalogues are publically available at http://star-www.dur.ac.uk/~cole/mocks/main.html . We expect these catalogues to be a valuable aid in the development of the new algorithms and statistics that will be used to analyse the 2dF and SDSS surveys when they are completed in the next few years. Mock catalogues of the PSCZ survey of IRAS galaxies are also available at the same WWW location.

We use a set of large, high-resolution cosmological N-body simulations to examine the redshift-space distortions of galaxy clustering on scales of order 10-200h^{-1} Mpc. Galaxy redshift surveys currently in progress will, on completion, allow us to measure the quadrupole distortion in the 2-point correlation function, \xi(\sigma,\pi), or its Fourier transform, the power spectrum, P(k,\mu), to a high degree of accuracy. On these scales we typically find a positive quadrupole, as expected for coherent infall onto overdense regions and outflow from underdense regions, but the distortion is substantially weaker than that predicted by pure linear theory. We assess two models that may be regarded as refinements to linear theory, the Zel'dovich approximation and a dispersion model in which the non-linear velocities generated by the formation of virialized groups and clusters are treated as random perturbations to the velocities predicted by linear theory. We find that neither provides an adequate physical description of the clustering pattern. If used to model redshift spacedistortions on scales for 10<\lambda <200 h^{-1}Mpc the estimated value of \beta (\beta=f(\Omega_0)/b where f(\Omega_0) ~ \Omega_0^{0.6} and b is the galaxy bias parameter) is liable to systematic errors of order ten per cent or more. We discuss how such systematics can be avoided by i) development of a more complete model of redshift distortions and ii) the direct use of galaxy catalogues generated from non-linear N-body simulations.