MoMaF: the Mock Map Facility

We present the Mock Map Facility, a powerful tool for converting theoretical outputs of hierarchical galaxy formation models into catalogues of virtual observations. The general principle is straightforward: mock observing cones can be generated using semi-analytically post-processed snapshots of cosmological N -body simulations. These cones can then be projected to synthesize mock sky images. To this end, the paper describes in detail an efficient technique for creating such mock cones and images from the galaxies in cosmological simulations ( galics ) semi-analytic model, providing the reader with an accurate quantification of the artefacts it introduces at every step. We show that replication effects introduce a negative bias on the clustering signal – typically peaking at less than 10 per cent 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 approximately one-tenth of the simulation box size. For the purpose of analysing our method, we show that number counts and redshift distributions obtained with galics / momaf compare well with K -band observations and the two-degree field galaxy redshift survey. Given finite-volume effects, we also show that the model can reproduce the automatic plate measuring machine angular correlation function. The momaf results discussed here are made publicly available to the astronomical community through a public data base. 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.     

QSO lensing magnification associated with galaxy groups

We simulated both the matter and light (galaxy) distributions in a wedge of the Universe and calculated the gravitational lensing magnification caused by the mass along the line-of-sight of galaxies and galaxy groups identified in sky surveys. A large volume redshift cone containing cold dark matter particles mimics the expected cosmological matter distribution in a flat universe with low matter density and a cosmological constant. We generate a mock galaxy catalogue from the matter distribution and identify thousands of galaxy groups in the luminous sky projection. We calculate the expected magnification around galaxies and galaxy groups and then the induced quasi-stellar object (QSO)–lens angular correlation due to magnification bias. This correlation is observable and can be used both to estimate the average mass of the lens population and to make cosmological inferences. We also use analytical calculations and various analyses to compare the observational results with theoretical expectations for the cross-correlation between faint QSOs from the 2dF Survey and nearby galaxies and groups from the Automated Plate Measurement and Sloan Digital Sky Survey Early Data Release. The observed QSO–lens anticorrelations are stronger than the predictions for the cosmological model used. This suggests that there could be unknown systematic errors in the observations and data reduction, or that the model used is not adequate. If the observed signal is assumed to be solely due to gravitational lensing, then the lensing is stronger than expected, due to more massive galactic structures or more efficient lensing than simulated.     

Topology and geometry of the CfA2 redshift survey

We analyse the redshift space topology and geometry of the nearby Universe by computing the Minkowski functionals of the Updated Zwicky Catalogue (UZC). The UZC contains the redshifts of almost 20 000 galaxies, is 96 per cent complete to the limiting magnitude m _Zw=15.5, and includes the Center for Astrophysics (CfA) Redshift Survey (CfA2). From the UZC we can extract volume-limited samples reaching a depth of 70  h ⁻¹ Mpc before sparse sampling dominates. We quantify the shape of the large-scale galaxy distribution by deriving measures of planarity and filamentarity from the Minkowski functionals. The nearby Universe shows a large degree of planarity and a small degree of filamentarity. This quantifies the sheet-like structure of the Great Wall, which dominates the northern region (CfA2N) of the UZC. We compare these results with redshift space mock catalogues constructed from high-resolution N -body simulations of two cold dark matter (CDM) models with either a decaying massive neutrino ( τ CDM) or a non-zero cosmological constant (ΛCDM). We use semi-analytic modelling to form and evolve galaxies in these dark matter‐only simulations. We are thus able, for the first time, to compile redshift space mock catalogues which contain galaxies, along with their observable properties, rather than dark matter particles alone. In both models the large-scale galaxy distribution is less coherent than the observed distribution, especially with regard to the large degree of planarity of the real survey. However, given the small volume of the region studied, this disagreement can still be a result of cosmic variance, as shown by the agreement between the ΛCDM model and the southern region of CfA2.     

Faint Lyman-break galaxies as a crucial test for galaxy formation models

It has recently been shown that galaxy formation models within the Λ cold dark matter cosmology predict that, compared to the observed population, small galaxies (with stellar masses  <10¹¹ M_⊙  ) form too early, are too passive since   z ∼ 3  and host too old stellar populations at   z = 0  . We then expect an overproduction of small galaxies at   z ≳ 4  that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the morgana galaxy formation model and grasil spectrophotometric  +  radiative transfer code to generate mock catalogues of deep fields observed with Hubble Space Telescope Advanced Camera for Surveys. We add observational noise and the effect of Lyman α emission, and perform colour–colour selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i -dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of   z ₈₅₀∼ 27 V   -dropouts at   z ∼ 5  . We quantify the properties of these 'excess' galaxies and discuss the implications: these galaxies are hosted in dark matter haloes with circular velocities in excess of 100 km s⁻¹, and their suppression may require a deep rethinking of stellar feedback processes taking place in galaxy formation.     

Light curves for off-centre ignition models of Type Ia supernovae

We compare observations of the high-redshift galaxy population to the predictions of the galaxy formation model of Croton et al. and De Lucia & Blaizot. This model, implemented on the Millennium Simulation of the concordance Lambda cold dark matter cosmogony, introduces 'radio mode' feedback from the central galaxies of groups and clusters in order to obtain quantitative agreement with the luminosity, colour, morphology and clustering properties of the present-day galaxy population. Here we construct deep light cone surveys in order to compare model predictions to the observed counts and redshift distributions of distant galaxies, as well as to their inferred luminosity and mass functions out to redshift 5. With the exception of the mass functions, all these properties are sensitive to modelling of dust obscuration. A simple but plausible treatment agrees moderately well with most of the data. The predicted abundance of relatively massive  (∼ M _*)  galaxies appears systematically high at high redshift, suggesting that such galaxies assemble earlier in this model than in the real Universe. An independent galaxy formation model implemented on the same simulation matches the observed mass functions slightly better, so the discrepancy probably reflects incomplete or inaccurate galaxy formation physics rather than problems with the underlying cosmogony.     

Testing model predictions of the cold dark matter cosmology for the sizes, colours, morphologies and luminosities of galaxies with the SDSS

The huge size and uniformity of the Sloan Digital Sky Survey (SDSS) make possible an exacting test of current models of galaxy formation. We compare the predictions of the galform semi-analytical galaxy formation model for the luminosities, morphologies, colours and scalelengths of local galaxies. galform models the luminosity and size of the disc and bulge components of a galaxy, and so we can compute quantities which can be compared directly with SDSS observations, such as the Petrosian magnitude and the Sérsic index. We test the predictions of two published models set in the cold dark matter cosmology: the Baugh et al. model, which assumes a top-heavy initial mass function (IMF) in starbursts and superwind feedback, and the Bower et al. model, which uses active galactic nucleus feedback and a standard IMF. The Bower et al. model better reproduces the overall shape of the luminosity function, the morphology–luminosity relation and the colour bimodality observed in the SDSS data, but gives a poor match to the size–luminosity relation. The Baugh et al. model successfully predicts the size–luminosity relation for late-type galaxies. Both models fail to reproduce the sizes of bright early-type galaxies. These problems highlight the need to understand better both the role of feedback processes in determining galaxy sizes, in particular the treatment of the angular momentum of gas reheated by supernovae, and the sizes of the stellar spheroids formed by galaxy mergers and disc instabilities.     

Galaxy–dark matter correlations applied to galaxy–galaxy lensing: predictions from the semi-analytic galaxy formation models

We use semi-analytic models of galaxy formation combined with high-resolution N -body simulations to make predictions for galaxy–dark matter correlations and apply them to galaxy–galaxy lensing. We analyse cross-power spectra between the dark matter and different galaxy samples selected by luminosity, colour or star formation rate. We compare the predictions with the recent detection by the Sloan Digital Sky Survey (SDSS). We show that the correlation amplitude and the mean tangential shear depend strongly on the luminosity of the sample on scales below 1  h ⁻¹ Mpc, reflecting the correlation between the galaxy luminosity and the halo mass. The cross-correlation cannot, however, be used to infer the halo profile directly because different halo masses dominate on different scales and because not all galaxies are at the centres of the corresponding haloes. We compute the redshift evolution of the cross-correlation amplitude and compare it with those of galaxies and dark matter. We also compute the galaxy–dark matter correlation coefficient and show that it is close to unity on scales above 1  h ⁻¹ Mpc for all considered galaxy types. This would allow one to extract the bias and the dark matter power spectrum on large scales from the galaxy and galaxy–dark matter correlations.     

The three-dimensional power spectrum of dark and luminous matter from the VIRMOS-DESCART cosmic shear survey

We present the first optimal power spectrum estimation and three-dimensional deprojections for the dark and luminous matter and their cross-correlations. The results are obtained using a new optimal fast estimator, deprojected using minimum variance and Singular Value Decomposition (SVD) techniques. We show the resulting 3D power spectra for dark matter and galaxies, and their covariance for the VIRMOS-DESCART weak lensing shear and galaxy data. The survey is most sensitive to non-linear scales   k _NL∼ 1 h Mpc⁻¹  . On these scales, our 3D power spectrum of dark matter is in good agreement with the RCS 3D power spectrum found by Tegmark & Zaldarriaga. Our galaxy power is similar to that found by the 2MASS survey, and larger than that of SDSS, APM and RCS, consistent with the expected difference in galaxy population.
We find an average bias   b = 1.24 ± 0.18  for the I -selected galaxies, and a cross-correlation coefficient   r = 0.75 ± 0.23  . Together with the power spectra, these results optimally encode the entire two point information about dark matter and galaxies, including galaxy–galaxy lensing. We address some of the implications regarding galaxy haloes and mass-to-light ratios. The best-fitting 'halo' parameter   h ≡ r / b = 0.57 ± 0.16  , suggesting that dynamical masses estimated using galaxies systematically underestimate total mass.
Ongoing surveys, such as the Canada–France–Hawaii Telescope Legacy Survey, will significantly improve on the dynamic range, and future photometric redshift catalogues will allow tomography along the same principles.     

Measuring the distribution of galaxies between haloes

We develop a method to measure the probability, P ( N;   M ), of finding N galaxies in a dark matter halo of mass M from the theoretically determined clustering properties of dark matter haloes and the observationally measured clustering properties of galaxies. Knowledge of this function and the distribution of the dark matter completely specifies all clustering properties of galaxies on scales larger than the size of dark matter haloes. Furthermore, P ( N;   M ) provides strong constraints on models of galaxy formation, since it depends upon the merger history of dark matter haloes and the galaxy–galaxy merger rate within haloes. We show that measurements from a combination of the Two Micron All Sky Survey and Sloan Digital Sky Survey or Two-degree Field Galaxy Redshift Survey data sets will allow P ( N;   M ) averaged over haloes occupied by bright galaxies to be accurately measured for N =0–2 .     

lemomaf: Lensed Mock Map Facility

We present the Lensed Mock Map Facility ( lemomaf ), a tool designed to perform mock weak-lensing measurements on numerically simulated chunks of the Universe. Coupling N -body simulations to a semi-analytical model of galaxy formation, lemomaf can create realistic lensed images and mock catalogues of galaxies, at wavelengths ranging from the ultraviolet to the submillimetre. To demonstrate the power of such a tool, we compute predictions of the source–lens clustering (SLC) effect on the convergence statistics, and quantify the impact of weak lensing on galaxy counts in two different filters. We find that the SLC effect skews the probability density function of the convergence towards low values, with an intensity which strongly depends on the redshift distribution of galaxies. On the other hand, the degree of enhancement or depletion in galaxy counts due to weak lensing is independent of the SLC effect. We discuss the impact on the two-point shear statistics to be measured by future missions like SNAP and LSST . The SLC effect would bias the estimation of σ₈ from two-point statistics up to 5 per cent for a narrow redshift distribution of mean   z ∼ 0.5  , and up to 2 per cent in small angular scales for a redshift distribution of mean   z ∼ 1.5  . We conclude that accurate photometric redshifts for individual galaxies are necessary in order to quantify and isolate the SLC effect.     

Modelling galaxy–galaxy weak lensing with Sloan Digital Sky Survey groups

We use galaxy groups selected from the Sloan Digital Sky Survey (SDSS) together with mass models for individual groups to study the galaxy–galaxy lensing signals expected from galaxies of different luminosities and morphological types. We compare our model predictions with the observational results obtained from the SDSS by Mandelbaum et al. for the same samples of galaxies. The observational results are well reproduced in a Λ cold dark matter (ΛCDM) model based on the Wilkinson Microwave Anisotropy Probe ( WMAP ) 3-yr data, but a ΛCDM model with higher σ₈, such as the one based on the WMAP 1-yr data, significantly overpredicts the galaxy–galaxy lensing signal. We model, separately, the contributions to the galaxy–galaxy lensing signals from different galaxies: central versus satellite, early type versus late type and galaxies in haloes of different masses. We also examine how the predicted galaxy–galaxy lensing signal depends on the shape, density profile and the location of the central galaxy with respect to its host halo.     

The properties of Lyα emitting galaxies in hierarchical galaxy formation models

We present detailed predictions for the properties of Lyα-emitting galaxies in the framework of the Λ cold dark matter cosmology, calculated using the semi-analytical galaxy formation model galform . We explore a model that assumes a top-heavy initial mass function in starbursts and that has previously been shown to explain the sub-millimetre number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Lyα photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Lyα emitters (LAEs) over the redshift range  3 < z < 6.6  . We also examine the distribution of Lyα equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of LAEs, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of LAEs at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of LAEs at   z > 7  , a redshift range that is starting to be be probed by near-infrared surveys and using new instruments such as the Dark Ages Z Lyman Explorer (DAzLE).     

Void statistics and void galaxies in the 2dF Galaxy Redshift Survey

In the 2dF Galaxy Redshift Survey, we study the properties of voids and of fainter galaxies within voids that are defined by brighter galaxies. Our results are compared with simulated galaxy catalogues from the Millennium simulation coupled with a semi-analytical galaxy formation recipe. We derive the void size distribution and discuss its dependence on the faint magnitude limit of the galaxies defining the voids. While voids among faint galaxies are typically smaller than those among bright galaxies, the ratio of the void sizes to the mean galaxy separation reaches larger values. This is well reproduced in the mock galaxy samples studied. We provide analytic fitting functions for the void size distribution. Furthermore, we study the galaxy population inside voids defined by galaxies with   B _J− 5 log  h < −20  and diameter larger than  10  h ⁻¹ Mpc  . We find a clear bimodality of galaxies inside voids and in the average field but with different characteristics. The abundance of blue cloud galaxies inside voids is enhanced. There is an indication of a slight blueshift of the blue cloud. Furthermore, galaxies in void centres have slightly higher specific star formation rates as measured by the η parameter. We determine the radial distribution of the ratio of early- and late-type galaxies through the voids. We find and discuss some differences between observations and the Millennium catalogues.     

A photometric study of the field around the candidate recoiling/binary black hole SDSS J092712.65+294344.0

We present a photometric far-ultraviolet (FUV) to K _s-band study of the field around quasar SDSS J092712.65+294344.0. The SDSS spectrum of this object shows various emission lines with two distinct redshifts, at   z = 0.699  and 0.712. Because of this peculiar spectroscopic feature, this source has been proposed as a candidate recoiling or binary black hole. A third alternative model involves two galaxies moving in the centre of a rich galaxy cluster. Here, we present a study addressing the possible presence of such a rich cluster of galaxies in the SDSS J092712.65+294344.0 field. We observed the  3.6 × 2.6  arcmin² field in the K _s band and matched the near-infrared data with the FUV and near-ultraviolet images in the Galaxy Evolution Explorer archive and the ugriz observations in the SDSS. From various colour–colour diagrams, we were able to classify the nature of 32 sources, only 6–11 of which have colours consistent with galaxies at   z ≈ 0.7  . We compare these numbers with the surface density of galaxies, stars and quasars and the expectations for typical galaxy clusters both at low and high redshift. Our study shows that the galaxy cluster scenario is in clear disagreement with the new observations.     

Are ring galaxies the ancestors of giant low surface brightness galaxies?

We simulate the collisional formation of a ring galaxy and we integrate its evolution up to 1.5 Gyr after the interaction. About 100–200 Myr after the collision, the simulated galaxy is very similar to observed ring galaxies (e.g. Cartwheel). After this stage, the ring keeps expanding and fades. Approximately 0.5–1 Gyr after the interaction, the disc becomes very large (∼100 kpc) and flat. Such extended discs have been observed only in giant low surface brightness galaxies (GLSBs). We compare various properties of our simulated galaxies (surface brightness profile, morphology, H  i spectrum and rotation curve) with the observations of four well-known GLSBs (UGC 6614, Malin 1, Malin 2 and NGC 7589). The simulations match quite well the observations, suggesting that ring galaxies could be the progenitors of GLSBs. This result is crucial for the cold dark matter (CDM) model, as it was very difficult, so far, to explain the formation of GLSBs within the CDM scenario.     

Satellite kinematics – II. The halo mass–luminosity relation of central galaxies in SDSS

The kinematics of satellite galaxies reflect the masses of the extended dark matter haloes in which they orbit, and thus shed light on the mass–luminosity relation (MLR) of their corresponding central galaxies. In this paper, we select a large sample of centrals and satellites from the Sloan Digital Sky Survey and measure the kinematics (velocity dispersions) of the satellite galaxies as a function of the r -band luminosity of the central galaxies. Using the analytical framework presented in More, van den Bosch & Cacciato, we use these data to infer both the mean and the scatter of the MLR of central galaxies, carefully taking account of selection effects and biases introduced by the stacking procedure. As expected, brighter centrals on average reside in more massive haloes. In addition, we find that the scatter in halo masses for centrals of a given luminosity,  σ_{log  M}   , also increases with increasing luminosity. As we demonstrate, this is consistent with  σ_{log  L}   , which reflects the scatter in the conditional probability function   P ( L _c| M )  , being independent of halo mass. Our analysis of the satellite kinematics yields  σ_{log  L} = 0.16  ±  0.04  , in excellent agreement with constraints from clustering and group catalogues, and with predictions from a semi-analytical model of galaxy formation. We thus conclude that the amount of stochasticity in galaxy formation, which is characterized by  σ_{log  L}   , is well constrained, independent of halo mass and in a good agreement with current models of galaxy formation.