Strong-lensing optical depths in a ΛCDM universe – II. The influence of the stellar mass in galaxies

We investigate how strong gravitational lensing in the concordance ΛCDM cosmology is affected by the stellar mass in galaxies. We extend our previous studies, based on ray tracing through the Millennium Simulation, by including the stellar components predicted by galaxy formation models. We find that the inclusion of these components greatly enhances the probability for strong lensing compared to a 'dark matter only' universe. The identification of the 'lenses' associated with strong-lensing events reveals that the stellar mass of galaxies (i) significantly enhances the strong-lensing cross-sections of group and cluster haloes and (ii) gives rise to strong lensing in smaller haloes, which would not produce noticeable effects in the absence of the stars. Even if we consider only image splittings ≳10 arcsec, the luminous matter can enhance the strong-lensing optical depths by up to a factor of 2.     

The image separation distribution of strong lenses: halo versus subhalo populations

We present a halo model prediction of the image separation distribution of strong lenses. Our model takes into account the subhalo population, which has been ignored in previous studies, as well as the conventional halo population. Haloes and subhaloes are linked to central and satellite galaxies by adopting a universal scaling relation between masses of (sub)haloes and luminosities of galaxies. Our model predicts that 10–20 per cent of lenses should be caused by the subhalo population. The fraction of lensing by satellite galaxies (subhaloes) peaks at ∼1 arcsec and decreases rapidly with increasing image separations. We compute fractions of lenses which lie in groups and clusters and find them to be ∼14 and ∼4 per cent, respectively; nearly half of such lenses are expected to be produced by satellite galaxies, rather than central parts of haloes. We also study mass distributions of lensing haloes and find that, even at image separations of ∼3 arcsec, the deviation of lens mass distributions from isothermal profiles is large; at or beyond ∼3 arcsec, image separations are enhanced significantly by surrounding haloes. Our model prediction agrees reasonably well with observed image separation distributions from galaxy to cluster scales.     

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.     

Degeneracy of semi-analytic solution for axially symmetric lens equations

We derive a simple semi-analytical approximation for lens equations with an arbitrary radially symmetric mass density ρ( r ), when   r /ξ₀≪ 1  and ξ₀ is the scalelength of the density profile. At the strong lensing regime, which is mostly constrained by the inner part of the mass density profile, we assume ρ∝ r ^α.
A dark matter (DM) haloes (GNFW model) are parametrized through a shape parameter α, a concentration parameter c ₁ and the total mass M . We apply our semi-analytical model to show how the solutions of the axially symmetric lens equations are degenerated in respect to the parameters α and c ₁.
In the case of an asymmetric dual image lens system, similar effective degeneracy is produced when the geometry of the lens is relaxed. Because it is impossible to determine the exact location of the source image, a family of solutions is acquired when the mass of the lens object and location of the observed images are fixed.
Our results indicate that the amount of degeneration is only weakly affected by the asymmetry in the lensing geometry set-up, e.g. the observational effective degeneracy is very close to the true physical degeneracy of the Einstein ring solutions. Basically with high-enough values for the concentration parameter, the degeneracy spawns the whole range for the shape parameter  α=[−2.0, −1.0]  .     

Constraining the Mass Distribution of Cluster Galaxies by Weak Lensing

Analysing the weak lensing distortions of the images of faint background galaxies provides a means to constrain the average mass distribution of cluster galaxies and potentially to test the extent of their dark matter haloes as a function of the density of their environment. The observable image distortions are a consequence of the interplay between the effects of a global cluster mass distribution and the perturbations resulting from individual cluster galaxies. Starting from a reconstruction of the cluster mass distribution with conventional techniques, we apply a maximum likelihood method to infer the average properties of an ensemble of cluster galaxies. From simulations this approach is found to be reliable as long as the galaxies including their dark matter haloes only contribute a small fraction to the total mass of the system. If their haloes are extended, the galaxies contain a substantial mass fraction. In this case our method is still applicable in the outer regions of clusters, where the surface mass density is low, but yields biased estimates of the parameters describing the mass profiles of the cluster galaxies in the central part of the cluster. In that case it will be necessary to resort to more sophisticated strategies by modelling cluster galaxies and an underlying global mass distribution simultaneously. We conclude that galaxy–galaxy lensing in clusters provides a unique means to probe the presence and extent of dark haloes of cluster galaxies.     

Gravitational lensing of extended high-redshift sources by dark matter haloes

High-redshift galaxies and quasi-stellar objects (QSOs) are most likely to be strongly lensed by intervening haloes between the source and the observer. In addition, a large fraction of lensed sources is expected to be seen in the submillimetre region, as a result of the enhanced magnification bias on the steep intrinsic number counts. We extend in three directions Blain's earlier study of this effect.
First, we use a modification of the Press–Schechter mass function and detailed lens models to compute the magnification probability distribution. We compare the magnification cross-sections of populations of singular isothermal spheres and Navarro, Frenk & White (NFW) haloes and find that they are very similar, in contrast to the image-splitting statistics which were recently investigated in other studies. The distinction between the two types of density profile is therefore irrelevant for our purposes.
Secondly, we discuss quantitatively the maximum magnification, μ _max, that can be achieved for extended sources (galaxies) with realistic luminosity profiles, taking into account the possible ellipticity of the lensing potential. We find that μ _max plausibly falls into the range     for sources of     effective radius at redshifts within     .
Thirdly, we apply our model for the lensing magnification to a class of sources following the luminosity evolution typical for a unified scheme of QSO formation. As a result of the peculiar steepness of their intrinsic number counts, we find that the lensed source counts at a fiducial wave length of 850 μm can exceed the unlensed counts by several orders of magnitude at flux densities ≳100 mJy, even with a conservative choice of the maximum magnification.     

Strong lensing optical depths in a ΛCDM universe

We investigate strong gravitational lensing in the concordance ΛCDM cosmology by carrying out ray tracing along past light cones through the Millennium Simulation, the largest simulation of cosmic structure formation ever carried out. We extend previous ray-tracing methods in order to take full advantage of the large volume and the excellent spatial and mass resolution of the simulation. As a function of source redshift we evaluate the probability that an image will be highly magnified, will be highly elongated or will be one of a set of multiple images. We show that such strong lensing events can almost always be traced to a single dominant lensing object and we study the mass and redshift distribution of these primary lenses. We fit analytic models to the simulated dark haloes in order to study how our optical depth measurements are affected by the limited resolution of the simulation and of the lensing planes that we construct from it. We conclude that such effects lead us to underestimate total strong lensing cross-sections by about 15 per cent. This is smaller than the effects expected from our neglect of the baryonic components of galaxies. Finally we investigate whether strong lensing is enhanced by material in front of or behind the primary lens. Although strong lensing lines of sight are indeed biased towards higher than average mean densities, this additional matter typically contributes only a few per cent of the total surface density.     

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.     

On Einstein clusters as galactic dark matter haloes

We consider global and gravitational lensing properties of the recently suggested Einstein clusters of weakly interacting massive particles (WIMPs) as galactic dark matter haloes. Being tangential pressure dominated, Einstein clusters are strongly anisotropic systems which can describe any galactic rotation curve by specifying the anisotropy. Due to this property, Einstein clusters may be considered as dark matter candidates. We analyse the stability of the Einstein clusters against both radial and non-radial pulsations, and we show that the Einstein clusters are dynamically stable. With the use of the Buchdahl type inequalities for anisotropic bodies, we derive upper limits on the velocity of the particles defining the cluster. These limits are consistent with those obtained from stability considerations. The study of light deflection shows that the gravitational lensing effect is slightly smaller for the Einstein clusters as compared to the singular isothermal density sphere model for dark matter. Therefore, lensing observations may discriminate, at least, in principle, between Einstein cluster and the other dark matter models.     

Statistical distribution of gravitational-lensing excursion angles: winding ways to us from the deep Universe

We investigate statistical distributions of differences in gravitational-lensing deflections between two light rays, the so-called lensing excursion angles. A probability distribution function of the lensing excursion angles, which plays a key role in estimates of lensing effects on angular clustering of objects (such as galaxies, quasi-stellar objects and also the cosmic microwave background temperature map), is known to consist of two components: a Gaussian core and an exponential tail. We use numerical gravitational-lensing experiments in a ΛCDM cosmology for quantifying these two components. We especially focus on the physical processes responsible for generating these two components. We develop a simple empirical model for the exponential tail which allows us to explore its origin. We find that the tail is generated by the coherent lensing scatter by massive haloes with   M > 10¹⁴  h ⁻¹ M_⊙  at   z < 1  and that its exponential shape arises due to the exponential cut-off of the halo mass function at that mass range. On scales larger than 1 arcmin, the tail does not have a practical influence on the lensing effects on the angular clustering. Our model predicts that the coherent scatter may have non-negligible effects on angular clustering at subarcminute scales.     

Numerical action reconstruction of the dynamical history of dark matter haloes in N-body simulations

We test the ability of the numerical action method (NAM) to recover the individual orbit histories of mass tracers in an expanding universe, given the masses and redshift-space coordinates at the present epoch. The mass tracers are represented by dark matter (DM) haloes identified in a region of radius  26  h ⁻¹ Mpc  of a high-resolution N -body simulation of the standard Λ cold dark matter (CDM) cosmology. Since previous tests of NAM at this scale have traced the underlying distribution of DM particles rather than extended haloes, our study offers an assessment of the accuracy of NAM in a scenario which more closely approximates the complex dynamics of actual galaxy haloes. We show that NAM can recover the present-day halo distances with typical distance errors of less than 3 per cent and radial peculiar velocities with a dispersion of  ∼130 km s⁻¹  . The accuracy of individual orbit reconstructions was limited by the inability of NAM, in some instances, to correctly model the positions of haloes at early times solely on the basis of the redshifts, angular positions and masses of the haloes at the present epoch. Improvements in the quality of NAM reconstructions may be possible using the present-day three-dimensional halo velocities and distances to further constrain the dynamics. This velocity data is expected to become available for nearby galaxies in the coming generations of observations by Space Interferometry Mission ( SIM ) and Global Astrometric Interferometer for Astrophysics ( GAIA ).     

The Sunyaev–Zel'dovich effect and Faraday rotation contributions of galaxy groups to the CMB angular power spectrum

The Sunyaev–Zel'dovich (SZ) effect and the Faraday rotation from haloes are examined over a wide mass range, including gas condensation and magnetic field evolution. Contributions to the cosmic microwave background (CMB) angular power spectrum are evaluated for galaxy clusters, galaxy groups and galaxies. Smaller mass haloes are found to play a more important role than massive haloes for the B -mode polarization associated with the SZ CMB anisotropies. The B modes from the Faraday rotation dominate the secondary B modes caused by gravitational lensing at  ℓ > 3000  . Measurement of B -mode polarization in combination with the SZ power spectrum can potentially provide important constraints on intracluster magnetic field and gas evolution at early epochs.     

Gravitational lensing effects on high-redshift Type II supernova studies with the Next Generation Space Telescope

We derive the expected Type II supernova (SN) differential number counts, N ( m ), and Hubble diagram for SCDM and LCDM cosmological models, taking into account the effects of gravitational lensing (GL) produced by the intervening cosmological mass. The mass distribution of dark matter haloes (i.e. the lenses) is obtained by means of a Monte Carlo method applied to the PressSchechter mass function. The haloes are assumed to have a Navarro, Frenk & White (NFW) density profile, in agreement with recent simulations of hierarchical cosmological models. Up to z =15, the (SCDM, LCDM) models predict a total number of (857, 3656) SNII yr¹ in 100 surveyed 44 arcmin² fields of the Next Generation Space Telescope ( NGST ). NGST will be able to reach the peak of the N ( m ) curve, located at AB 30(31) for SCDM (LCDM) in J and K wavelength bands, and detect (75 per cent, 51 per cent) of the above SN events. This will allow a detailed study of the early cosmic star formation history, as traced by SNIIe. N ( m ) is only very mildly affected by the inclusion of lensing effects. In addition, GL introduces a moderate uncertainty in the determination of cosmological parameters from Hubble diagrams, when these are pushed to higher z . For example, for a 'true' LCDM with (_M=0.4, =0.6), without proper account of GL, one would instead derive We briefly compare our results with previous similar work and discuss the limitations of the model.     

On the interaction of jet plasma with dark matter subhaloes in active galaxies

The interaction of fragmented plasma of active galactic nuclei jets with galactic haloes via gravitational scattering and lensing by dark matter subhaloes is studied using analytical calculations and numerical Monte-Carlo method. The lensing of jet radiation by halo masses is found to be negligible and unobservable. Moving through a galactic halo jet plasma fragments are sequentially deflected on hyperbolic orbits by gravitational field of subhaloes and deviates at some angles when leaving halo, causing widening of the jet. Based on this model jet opening angles are calculated numerically for various values of jet and halo characteristics. Though these angles are very small, gravitational scattering by halo masses results in specific radial profile of jet radiation intensity, that does not depend on halo mass distribution and jet properties. The intensity of jet radiation, obeying the derived profile, decreases by reasonable observable factors giving possibility to probe the presence of dark matter subhaloes.     

Imaging the cosmic matter distribution using gravitational lensing of pre-galactic H i

21-cm emission from neutral hydrogen during and before the epoch of cosmic reionization is gravitationally lensed by material at all lower redshifts. Low-frequency radio observations of this emission can be used to reconstruct the projected mass distribution of foreground material, both light and dark. We compare the potential imaging capabilities of such 21-cm lensing with those of future galaxy lensing surveys. We use the Millennium Simulation to simulate large-area maps of the lensing convergence with the noise, resolution and redshift-weighting achievable with a variety of idealized observation programmes. We find that the signal-to-noise ratio of 21-cm lens maps can far exceed that of any map made using galaxy lensing. If the irreducible noise limit can be reached with a sufficiently large radio telescope, the projected convergence map provides a high-fidelity image of the true matter distribution, allowing the dark matter haloes of individual galaxies to be viewed directly, and giving a wealth of statistical and morphological information about the relative distributions of mass and light. For instrumental designs like that planned for the Square Kilometre Array, high-fidelity mass imaging may be possible near the resolution limit of the core array of the telescope.     

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.     

The Hubble constant from galaxy lenses: impacts of triaxiality and model degeneracies

The Hubble constant can be constrained using the time delays between multiple images of gravitationally lensed sources. In some notable cases, typical lensing analyses assuming isothermal galaxy density profiles produce low values for the Hubble constant, inconsistent with the result of the HST Key Project  (72 ± 8 km s⁻¹ Mpc⁻¹)  . Possible systematics in the values of the Hubble constant derived from galaxy lensing systems can result from a number of factors, for example, neglect of environmental effects, assumption of isothermality, or contamination by line-of-sight structures. One additional potentially important factor is the triaxial structure of the lensing galaxy halo; most lens models account for halo shape simply by perturbing the projected spherical lensing potential, an approximation that is often necessary but that is inadequate at the levels of triaxiality predicted in the cold dark matter paradigm. To quantify the potential error introduced by this assumption in estimates of the Hubble parameter, we strongly lens a distant galaxy through a sample of triaxial softened isothermal haloes and use an Markov Chain Monte Carlo method to constrain the lensing halo profile and the Hubble parameter from the resulting multiple image systems. We explore the major degeneracies between the Hubble parameter and several parameters of the lensing model, finding that without a way to accurately break these degeneracies accurate estimates of the Hubble parameter are not possible. Crucially, we find that triaxiality does not significantly bias estimates of the Hubble constant, and offer an analytic explanation for this behaviour in the case of isothermal profiles. Neglected triaxial halo shape cannot contribute to the low Hubble constant values derived in a number of galaxy lens systems.     

The dark matter distribution in disc galaxies

We use high-quality optical rotation curves of nine low-luminosity disc galaxies to obtain the velocity profiles of the surrounding dark matter haloes. We find that they increase linearly with radius at least out to the edge of the stellar disc, implying that, over the entire stellar region, the density of the dark halo is about constant.
The properties of the mass structure of these haloes are similar to those found for a number of dwarf and low surface brightness galaxies, but provide a more substantial evidence of the discrepancy between the halo mass distribution predicted in the cold dark matter scenario and those actually detected around galaxies. We find that the density law proposed by Burkert reproduces the halo rotation curves, with halo central densities ( ρ ₀∼1–4×10⁻²⁴ g cm⁻³) and core radii ( r ₀∼5–15 kpc) scaling as ρ ₀∝ r ₀^−2/3.     

Understanding the halo-mass and galaxy-mass cross-correlation functions

We use the Millennium Simulation (MS) to measure the cross-correlation between halo centres and mass (or equivalently the average density profiles of dark haloes) in a Lambda cold dark matter (ΛCDM) cosmology. We present results for radii in the range  10  h ⁻¹ kpc < r < 30  h ⁻¹ Mpc  and for halo masses in the range  4 × 10¹⁰ < M ₂₀₀ < 4 × 10¹⁴  h ⁻¹ M_⊙  . Both at   z = 0  and at   z = 0.76  these cross-correlations are surprisingly well fitted if the inner region is approximated by a density profile of NFW or Einasto form, the outer region by a biased version of the linear mass autocorrelation function, and the maximum of the two is adopted where they are comparable. We use a simulation of galaxy formation within the MS to explore how these results are reflected in cross-correlations between galaxies and mass. These are directly observable through galaxy–galaxy lensing. Here also we find that simple models can represent the simulation results remarkably well, typically to ≲10 per cent. Such models can be used to extend our results to other redshifts, to cosmologies with other parameters, and to other assumptions about how galaxies populate dark haloes. Our galaxy formation simulation already reproduces current galaxy–galaxy lensing data quite well. The characteristic features predicted in the galaxy–galaxy lensing signal should provide a strong test of the ΛCDM cosmology as well as a route to understanding how galaxies form within it.     

Probing galactic dark matter in dense environments: on the strong lensing efficiency of galaxies in rich clusters

The recent detection by Limousin et al. of five new strong lensing events dominated by galaxy cluster members in Abell 1689, and outside the critical regime of the cluster itself, offers a way to obtain constraints on the cluster mass distribution in a region inaccessible to standard lensing analysis. In addition, modelling such systems will provide another window on the dark matter haloes of galaxies in very dense environments. Here, it is shown that the boost in image separation due to the external shear and convergence from a smooth cluster component means that more numerous, less massive galaxies have the potential to create multiple images with detectable separations, relative to isolated field galaxies. This comes in addition to a potential increase in their lensing (source plane) cross-section. To gain insight into the factors involved and as a precursor to a numerical study using N -body simulations, a simple analytic model of a cluster at   z = 0.3  lensing background galaxies at   z = 2  is considered here. The fiducial model has cluster members with isothermal density profiles and luminosities L , distributed in a Schechter function (faint-end slope  ν=−1.25  ), related to their velocity dispersions σ via the Faber–Jackson scaling L ∝σ⁴. Just outside the critical regime of the cluster, the scale of galaxy-dominated image separations is significantly increased. Folding in the fact that less massive galaxies present a lower lensing cross-section, and that the cross-section can itself be enhanced in an external field leads to a factor of a few times more detected events relative to field galaxies. These values will be higher closer to the critical curve. Given that the events in Abell 1689 were detected over a very small region of the cluster where ACS data were available, this motivates the search for such events in other clusters.