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 imprint of cosmic reionization on galaxy clustering

We consider the effect of reionization on the clustering properties of galaxy samples at intermediate redshifts ( z ∼ 0.3–5.5). Current models for the reionization of intergalactic hydrogen predict that overdense regions will be reionized early, thus delaying the build-up of stellar mass in the progenitors of massive lower redshift galaxies. As a result, the stellar populations observed in intermediate-redshift galaxies are somewhat younger and hence brighter in overdense regions of the Universe. Galaxy surveys would therefore be sensitive to galaxies with a somewhat lower dark matter mass in overdense regions. The corresponding increase in the observed number density of galaxies can be parametrized as a galaxy bias due to reionization. We model this process using merger trees combined with a stellar synthesis code. Our model demonstrates that reionization has a significant effect on the clustering properties of galaxy samples that are selected based on their star formation properties. The bias correction in Lyman-break galaxies (including those in proposed baryonic oscillation surveys at z < 1) is at the level of 10–20 per cent for a halo mass of  10¹² M_⊙  , leading to corrections factors of 1.5–2 in the halo mass inferred from measurements of clustering length. The reionization of helium could also lead to a sharp increase in the amplitude of the galaxy correlation function at z ∼ 3. We find that the reionization bias is approximately independent of scale and halo mass. However, since the traditional galaxy bias is mass dependent, the reionization bias becomes relatively more important for lower mass systems. The correction to the bias due to reionization is very small in surveys of luminous red galaxies at z < 1.     

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 .     

The dependence of galaxy formation on cosmological parameters: can we distinguish between the WMAP1 and WMAP3 parameter sets?

We combine N -body simulations of structure growth with physical modelling of galaxy evolution to investigate whether the shift in cosmological parameters between the first- and third-year results from the Wilkinson Microwave Anisotropy Probe ( WMAP ) affects predictions for the galaxy population. Structure formation is significantly delayed in the WMAP3 cosmology, because the initial matter fluctuation amplitude is lower on the relevant scales. The decrease in dark matter clustering strength is, however, almost entirely offset by an increase in halo bias, so predictions for galaxy clustering are barely altered. In both cosmologies, several combinations of physical parameters can reproduce observed, low-redshift galaxy properties; the star formation, supernova feedback and active galactic nucleus feedback efficiencies can be played off against each other to give similar results. Models which fit observed luminosity functions predict projected two-point correlation functions which scatter by about 10–20 per cent on large scale and by larger factors on small scale, depending both on cosmology and on details of galaxy formation. Measurements of the pairwise velocity distribution prefer the WMAP1 cosmology, but careful treatment of the systematics is needed. Given present modelling uncertainties, it is not easy to distinguish between the WMAP1 and WMAP3 cosmologies on the basis of low-redshift galaxy properties. Model predictions diverge more dramatically at high redshift. Better observational data at   z > 2  will better constrain galaxy formation and perhaps also cosmological parameters.     

SDSS galaxy clustering: luminosity and colour dependence and stochasticity

Differences in clustering properties between galaxy subpopulations complicate the cosmological interpretation of the galaxy power spectrum, but can also provide insights about the physics underlying galaxy formation. To study the nature of this relative clustering, we perform a counts-in-cells analysis of galaxies in the Sloan Digital Sky Survey in which we measure the relative bias between pairs of galaxy subsamples of different luminosities and colours. We use a generalized  χ²  test to determine if the relative bias between each pair of subsamples is consistent with the simplest deterministic linear bias model, and we also use a maximum likelihood technique to further understand the nature of the relative bias between each pair. We find that the simple, deterministic model is a good fit for the luminosity-dependent bias on scales above  ∼2  h ⁻¹ Mpc  , which is good news for using magnitude-limited surveys for cosmology. However, the colour-dependent bias shows evidence for stochasticity and/or non-linearity which increases in strength towards smaller scales, in agreement with previous studies of stochastic bias. Also, confirming hints seen in earlier work, the luminosity-dependent bias for red galaxies is significantly different from that of blue galaxies: both luminous and dim red galaxies have higher bias than moderately bright red galaxies, whereas the biasing of blue galaxies is not strongly luminosity dependent. These results can be used to constrain galaxy formation models and also to quantify how the colour and luminosity selection of a galaxy survey can impact measurements of the cosmological matter power spectrum.     

Intrinsic correlation of galaxy shapes: implications for weak lensing measurements

Weak gravitational lensing is now established as a powerful method to measure mass fluctuations in the universe. It relies on the measurement of small coherent distortions of the images of background galaxies. Even low-level correlations in the intrinsic shapes of galaxies could however produce a significant spurious lensing signal. These correlations are also interesting in their own right, since their detection would constrain models of galaxy formation. Using     haloes found in N -body simulations, we compute the correlation functions of the intrinsic ellipticity of spiral galaxies assuming that the disc is perpendicular to the angular momentum of the dark matter halo. We also consider a simple model for elliptical galaxies, in which the shape of the dark matter halo is assumed to be the same as that of the light. For deep lensing surveys with median redshifts ∼1, we find that intrinsic correlations of ∼10⁻⁴ on angular scales     are generally below the expected lensing signal, and contribute only a small fraction of the excess signals reported on these scales. On larger scales we find limits to the intrinsic correlation function at a level ∼10⁻⁵, which gives a (model-dependent) range of separations for which the intrinsic signal is about an order of magnitude below the ellipticity correlation function expected from weak lensing. Intrinsic correlations are thus negligible on these scales for dedicated weak lensing surveys. For wider but shallower surveys such as SuperCOSMOS, APM and SDSS, we cannot exclude the possibility that intrinsic correlations could dominate the lensing signal. We discuss how such surveys could be used to calibrate the importance of this effect, as well as study spin–spin correlations of spiral galaxies.     

Scale-dependent galaxy bias in the Sloan Digital Sky Survey as a function of luminosity and colour

It has been known for a long time that the clustering of galaxies changes as a function of galaxy type. This galaxy bias acts as a hindrance to the extraction of cosmological information from the galaxy power spectrum or correlation function. Theoretical arguments show that a change in the amplitude of the clustering between galaxies and mass on large scales is unavoidable, but cosmological information can be easily extracted from the shape of the power spectrum or correlation function if this bias is independent of scale. Scale-dependent bias is generally small on large scales,   k < 0.1  h  Mpc⁻¹  , but on smaller scales can affect the recovery of  Ω_m h   from the measured shape of the clustering signal, and have a small effect on the Baryon Acoustic Oscillations. In this paper, we investigate the transition from scale-independent to scale-dependent galaxy bias as a function of galaxy population. We use the Sloan Digital Sky Survey Data Release 5 sample to fit various models, which attempt to parametrize the turn-off from scale-independent behaviour. For blue galaxies, we find that the strength of the turn-off is strongly dependent on galaxy luminosity, with stronger scale-dependent bias on larger scales for more luminous galaxies. For red galaxies, the scale dependence is a weaker function of luminosity. Such trends need to be modelled in order to optimally extract the information available in future surveys, and can help with the design of such surveys.     

Intrinsic and extrinsic galaxy alignment

We show with analytic models that the assumption of uncorrelated intrinsic ellipticities of target sources that is usually made in searches for weak gravitational lensing arising from large-scale mass inhomogeneities ('field lensing') is unwarranted. If the orientation of the galaxy image is determined either by the angular momentum or by the shape of the halo in which it forms, then the image should be aligned preferentially with the component of the tidal gravitational field perpendicular to the line of sight. Long-range correlations in the tidal field will thus lead to long-range ellipticity–ellipticity correlations that mimic the shear correlations arising from weak gravitational lensing. We calculate the ellipticity–ellipticity correlation expected if halo shapes determine the observed galaxy shape, and we discuss uncertainties (which are still considerable) in the predicted amplitude of this correlation. The ellipticity–ellipticity correlation induced by angular momenta should be smaller. We consider several methods for discriminating between the weak-lensing (extrinsic) and intrinsic correlations, including the use of redshift information. An ellipticity–tidal-field correlation also implies the existence of an alignment of images of galaxies near clusters. Although the intrinsic alignment may complicate the interpretation of field-lensing results, it is inherently interesting as it may shed light on galaxy formation as well as on structure formation.     

The mass of the Andromeda galaxy

This paper argues that the Milky Way galaxy is probably the largest member of the Local Group. The evidence comes from estimates of the total mass of the Andromeda galaxy (M31) derived from the three-dimensional positions and radial velocities of its satellite galaxies, as well as the projected positions and radial velocities of its distant globular clusters and planetary nebulae. The available data set comprises 10 satellite galaxies, 17 distant globular clusters and nine halo planetary nebulae with radial velocities. We find that the halo of Andromeda has a mass of together with a scalelength of 90 kpc and a predominantly isotropic velocity distribution. For comparison, our earlier estimate for the Milky Way halo is Although the error bars are admittedly large, this suggests that the total mass of M31 is probably less than that of the Milky Way . We verify the robustness of our results to changes in the modelling assumptions and to errors caused by the small size and incompleteness of the data set.
Our surprising claim can be checked in several ways in the near future. The numbers of satellite galaxies, planetary nebulae and globular clusters with radial velocities can be increased by ground-based spectroscopy, while the proper motions of the companion galaxies and the unresolved cores of the globular clusters can be measured using the astrometric satellites Space Interferometry Mission ( SIM ) and Global Astrometric Interferometer for Astrophysics ( GAIA ). Using 100 globular clusters at projected radii 20 R 50 kpc with both radial velocities and proper motions, it will be possible to estimate the mass within 50 kpc to an accuracy of 20 per cent. Measuring the proper motions of the companion galaxies with SIM and GAIA will reduce the uncertainty in the total mass caused by the small size of the data set to 22 per cent.     

Theoretical estimates of intrinsic galaxy alignment

It has recently been argued that the observed ellipticities of galaxies may be determined at least in part by the primordial tidal gravitational field in which the galaxy formed. Long-range correlations in the tidal field could thus lead to an ellipticity–ellipticity correlation for widely separated galaxies. We present a new model relating ellipticity to angular momentum, which can be calculated in linear theory. We use this model to calculate the angular power spectrum of intrinsic galaxy shape correlations. We show that, for low-redshift galaxy surveys, our model predicts that intrinsic correlations will dominate correlations induced by weak lensing, in good agreement with previous theoretical work and observations. We find that our model produces ' E -mode' correlations enhanced by a factor of 3.5 over B -modes on small scales, making it harder to disentangle intrinsic correlations from those induced by weak gravitational lensing.     

Tracing the nature of dark energy with galaxy distribution

Dynamical dark energy (DE) is a viable alternative to the cosmological constant. Constructing tests to discriminate between Λ and dynamical DE models is difficult, however, because the differences are not large. In this paper we explore tests based on the galaxy mass function, the void probability function (VPF), and the number of galaxy clusters. At high z , the number density of clusters shows large differences between DE models, but geometrical factors reduce the differences substantially. We find that detecting a model dependence in the cluster redshift distribution is a significant challenge. We show that the galaxy redshift distribution is potentially a more sensitive characteristic. We do this by populating dark matter haloes in N -body simulations with galaxies using well-tested halo occupation distributions. We also estimate the VPF and find that samples with the same angular surface density of galaxies, in different models, exhibition almost model-independent VPF which therefore cannot be used as a test for DE. Once again, geometry and cosmic evolution compensate each other. By comparing VPFs for samples with fixed galaxy mass limits, we find measurable differences.     

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).     

Modelling galaxy clustering by colour

We extend the mass-halo formalism for analytically generating power spectra to allow for the different clustering behaviour observed in galaxy subpopulations. Although applicable to other separations, we concentrate our methods on a simple separation by rest-frame colour into 'red' and 'blue' subpopulations through modifications to the 〈 N 〉( M ) relations and halo distribution functions for each of the subpopulations. This sort of separation is within the capabilities of the current generations of simulations as well as galaxy surveys, suggesting a potentially powerful observational constraint for current and future simulations. In anticipation of this, we demonstrate the sensitivity of the resulting power spectra to the choice of model parameters.     

Halo masses for optically selected and for radio-loud AGN from clustering and galaxy–galaxy lensing

We compute two-point correlation functions and measure the shear signal due to galaxy–galaxy lensing for 80 000 optically identified and 5700 radio-loud active galactic nuclei (AGN) from Data Release 4 of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed   M _*  . Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive haloes than optical AGN: in our samples, their mean halo masses are  1.6 × 10¹³  and  8 × 10¹¹  h ⁻¹ M_⊙  , respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter haloes of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in haloes where gas cools quasi-statically.