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.     

Testing Bekenstein's relativistic Modified Newtonian Dynamics with lensing data

We propose to use multiple-imaged gravitational lenses to set limits on gravity theories without dark matter, specifically tensor–vector–scalar (TeVeS) theory, a theory which is consistent with fundamental relativistic principles and the phenomenology of Modified Newtonian Dynamics (MOND) theory. After setting the framework for lensing and cosmology, we analytically derive the deflection angle for the point lens and the Hernquist galaxy profile, and study their patterns in convergence, shear and amplification. Applying our analytical lensing models, we fit galaxy-quasar lenses in the CfA-Arizona Space Telescope Lens Survey (CASTLES) sample. We do this with three methods, fitting the observed Einstein ring sizes, the image positions, or the flux ratios. In all the cases, we consistently find that stars in galaxies in MOND/TeVeS provide adequate lensing. Bekenstein's toy μ function provides more efficient lensing than the standard MOND μ function. But for a handful of lenses, a good fit would require a lens mass orders of magnitude larger/smaller than the stellar mass derived from luminosity unless the modification function μ and modification scale a ₀ for the universal gravity were allowed to be very different from what spiral galaxy rotation curves normally imply. We discuss the limitation of present data and summarize constraints on the MOND μ function. We also show that the simplest TeVeS 'minimal-matter' cosmology, a baryonic universe with a cosmological constant, can fit the distance–redshift relation from the supernova data, but underpredicts the sound horizon size at the last scattering. We conclude that lensing is a promising approach to differentiate laws of gravity.     

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.     

The impact of non-Gaussian errors on weak lensing surveys

The weak lensing power spectrum carries cosmological information via its dependence on the growth of structure and on geometric factors. Since much of the cosmological information comes from scales affected by non-linear clustering, measurements of the lensing power spectrum can be degraded by non-Gaussian covariances. Recently, there have been conflicting studies about the level of this degradation. We use the halo model to estimate it and include new contributions related to the finite size of lensing surveys, following Rimes and Hamilton's study of three-dimensional simulations. We find that non-Gaussian correlations between different multipoles can degrade the cumulative signal-to-noise ratio (S/N) for the power spectrum amplitude by up to a factor of 2 (or 5 for a worst-case model that exceeds current N -body simulation predictions). However, using an eight-parameter Fisher analysis, we find that the marginalized errors on individual parameters are degraded by less than 10 per cent (or 20 per cent for the worst-case model). The smaller degradation in parameter accuracy is primarily because: individual parameters in a high-dimensional parameter space are degraded much less than the volume of the full Fisher ellipsoid; lensing involves projections along the line of sight, which reduce the non-Gaussian effect; some of the cosmological information comes from geometric factors which are not degraded at all. We contrast our findings with those of Lee and Pen who suggested a much larger degradation in information content. Finally, our results give a useful guide for exploring survey design by giving the cosmological information returns for varying survey area, depth and the level of some systematic errors.     

Sources of contamination to weak lensing three-point statistics: constraints from N-body simulations

We investigate the impact of the observed correlation between a galaxy's shape and its surrounding density field on the measurement of third-order weak lensing shear statistics. Using numerical simulations, we estimate the systematic error contribution to a measurement of the third-order moment of the aperture mass statistic (GGG) from three-point intrinsic ellipticity correlations (III), and the three-point coupling between the weak lensing shear experienced by distant galaxies and the shape of foreground galaxies (GGI and GII). We find that third-order weak lensing statistics are typically more strongly contaminated by these physical systematics compared to second-order shear measurements, contaminating the measured three-point signal for moderately deep surveys with a median redshift   z _m∼ 0.7  by ∼15 per cent. It has been shown that accurate photometric redshifts will be crucial to correct for this effect, once a model and the redshift dependence of the effect can be accurately constrained. To this end we provide redshift-dependent fitting functions to our results and propose a new tool for the observational study of intrinsic galaxy alignments. For a shallow survey with   z _m∼ 0.4  we find III to be an order of magnitude larger than the expected cosmological GGG shear signal. Compared to the two-point intrinsic ellipticity correlation which is similar in amplitude to the two-point shear signal at these survey depths, third-order statistics therefore offer a promising new way to constrain models of intrinsic galaxy alignments. Early shallow data from the next generation of very wide weak lensing surveys will be optimal for this type of study.     

Precision of diffuse 21-cm lensing

We study the limits of accuracy for weak lensing maps of dark matter using diffuse 21-cm radiation from the pre-reionization epoch using simulations. We improve on previous 'optimal' quadratic lensing estimators by using shear and convergence instead of deflection angles. This is a generalization of the deflection estimator, and is more optimal for non-Gaussian sources. The cross-power spectrum of shear and convergence is an unbiased estimator of lensing power spectrum which does not require knowledge of the source four-point function. We find that non-Gaussianity provides a limit to the accuracy of weak lensing reconstruction, even if instrumental noise is reduced to zero. The best reconstruction result is equivalent to Gaussian sources with effective independent cell of side length  2.0  h ⁻¹ Mpc  . Using a source full map from z = 10 to 20, this limiting sensitivity allows mapping of dark matter at a signal-to-noise ratio greater than 1 out to l ≲ 6000, which is better than any other proposed technique for large-area weak lensing mapping.     

An upper limit to the density inhomogeneity in the universe

An upper limit to the amplitude of the overall density fluctuation has been found by means of the gravitational lensing effect of the density inhomogeneity on the luminosities of quasars with larger redshifts. The observed differences of luminosities of quasars located at different directions are partially given by the lensing effect, therefore, a useful upper limit to the inhomogeneity can be derived if the luminosity distribution of quasars is uniform enough. We obtain that, in the case of the density parameter of the Universe =1, the overall matter should be less clustered than the luminous matter by a factor of at least 3. This result may not favour the biased dark matter scenario for the formation of large-scale structure in the Universe.     

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.     

Galaxy density profiles and shapes – I. Simulation pipeline for lensing by realistic galaxy models

Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, which we aim to investigate in a consistent way. In this first paper in a series, we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii.     

ARCS, the Arcminute Radio Cluster-lens Search – I. Selection criteria and initial results

We present the results of an unbiased radio search for gravitational lensing events with image separations between 15 and 60 arcsec, which would be associated with clusters of galaxies with masses >10^13–14 M_⊙. A parent population of 1023 extended radio sources stronger than 35 mJy with stellar optical identifications was selected using the FIRST radio catalogue at 1.4 GHz and the APM optical catalogue. The FIRST catalogue was then searched for companions to the parent sources stronger than 7 mJy and with separation in the range 15 to 60 arcsec. Higher-resolution observations of the resulting 38 lens candidates were made with the VLA at 1.4 and 5 GHz, and with MERLIN at 5 GHz in order to test the lens hypothesis in each case. None of our targets was found to be a gravitational lens system. These results provide the best current constraint on the lensing rate for this angular scale, but improved calculations of lensing rates from realistic simulations of the clustering of matter on the relevant scales are required before cosmologically significant constraints can be derived from this null result. We now have an efficient, tested observational strategy with which it will be possible to make an order-of-magnitude larger unbiased search in the near future.     

Gravitational Lensing Effects in a Modified General Relativity Model

We investigate the effects of gravitational lensing in a modified generalrelativity flat cosmological model. We show that the predicted lensingprobability for low matter densities is larger than in the Einstein-deSitter model but considerably less than in constant- models. Wealso find that the optical depth at low density rises with the sourcered-shift more slowly than in constant- models and that the lensdistribution peaks at lower values of the lens red-shift.     

Full-sky maps for gravitational lensing of the cosmic microwave background

We use the large cosmological Millennium Simulation (MS) to construct the first all-sky maps of the lensing potential and the angle, aiming at gravitational lensing of the cosmic microwave background (CMB), with the goal of properly including small-scale non-linearities and non-Gaussianity. Exploiting the Born approximation, we implement a map-making procedure based on direct ray tracing through the gravitational potential of the MS. We stack the simulation box in redshift shells up to z ∼ 11, producing continuous all-sky maps with arcmin angular resolution. A randomization scheme avoids the repetition of structures along the line of sight, and structures larger than the MS box size are added to supply the missing contribution of large-scale (LS) structures to the lensing signal. The angular power spectra of the projected lensing potential and the deflection-angle modulus agree quite well with semi-analytic estimates on scales down to a few arcmin, while we find a slight excess of power on small scales, which we interpret as being due to non-linear clustering in the MS. Our map-making procedure, combined with the LS adding technique, is ideally suited for studying lensing of CMB anisotropies, for analysing cross-correlations with foreground structures, or other secondary CMB anisotropies such as the Rees–Sciama effect.     

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.     

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.     

Weak lensing magnification reconstruction with the modified internal linear combination method

Measuring weak lensing cosmic magnification signal is very challenging due to the overwhelming intrinsic clustering in the observed galaxy distribution.In this paper,we modify the Internal Linear Combination (ILC) method to reconstruct the lensing signal with an extra constraint to suppress the intrinsic clustering.To quantify the performance,we construct a realistic galaxy catalogue for the LSST-like photometric survey,covering 20 000 deg~2 with mean source redshift at z_s~1.We find that the reconstruction performance depends on the width of the photo-z bin we choose.Due to the correlation between the lensing signal and the source galaxy distribution,the derived signal has smaller systematic bias but larger statistical uncertainty for a narrower photo-z bin.We conclude that the lensing signal reconstruction with the Modified ILC method is unbiased with a statistical uncertainty5%for bin width?z~P=0.2.     

The dark matter environment of the Abell 901/902 supercluster: a weak lensing analysis of the HST STAGES survey

We present a high-resolution dark matter reconstruction of the   z = 0.165  Abell 901/902 supercluster from a weak lensing analysis of the Hubble Space Telescope STAGES survey. We detect the four main structures of the supercluster at high significance, resolving substructure within and between the clusters. We find that the distribution of dark matter is well traced by the cluster galaxies, with the brightest cluster galaxies marking out the strongest peaks in the dark matter distribution. We also find a significant extension of the dark matter distribution of Abell 901a in the direction of an infalling X-ray group Abell 901α. We present mass, mass-to-light and mass-to-stellar mass ratio measurements of the structures and substructures that we detect. We find no evidence for variation of the mass-to-light and mass-to-stellar mass ratio between the different clusters. We compare our space-based lensing analysis with an earlier ground-based lensing analysis of the supercluster to demonstrate the importance of space-based imaging for future weak lensing dark matter 'observations'.     

Quasar–galaxy associations revisited

Gravitational lensing predicts an enhancement of the density of bright, distant quasi-stellar objects (QSOs) around foreground galaxies. We measure this QSO–galaxy correlation w _qg for two complete samples of radio-loud quasars, the southern 1 Jy and half-Jansky samples. The existence of a positive correlation between     quasars and     galaxies is confirmed at a     per cent significance level (>99.9 per cent if previous measurements on the northern hemisphere are included). A comparison with the results obtained for incomplete quasar catalogues (e.g. the Veron-Cetty and Veron compilation) suggests the existence of an 'identification bias', which spuriously increases the estimated amplitude of w _qg for incomplete samples. This effect may explain many of the surprisingly strong quasar–galaxy associations found in the literature. Nevertheless, the value of w _qg that we measure in our complete catalogues is still considerably higher than the predictions from weak lensing. We consider two effects which could help to explain this discrepancy: galactic dust extinction and strong lensing.     

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.     

Non-parametric strong lens inversion of Cl 0024+1654: illustrating the monopole degeneracy

The cluster lens Cl 0024+1654 is undoubtedly one of the most beautiful examples of strong gravitational lensing, providing five large images of a single source with well-resolved substructure. Using the information contained in the positions and the shapes of the images, combined with the null space information, a non-parametric technique is used to infer the strong lensing mass map of the central region of this cluster. This yields a strong lensing mass of  1.60 × 10¹⁴ M_⊙  within a 0.5  arcmin radius around the cluster centre. This mass distribution is then used as a case study of the monopole degeneracy, which may be one of the most important degeneracies in gravitational lensing studies and which is extremely hard to break. We illustrate the monopole degeneracy by adding circularly symmetric density distributions with zero total mass to the original mass map of Cl 0024+1654. These redistribute mass in certain areas of the mass map without affecting the observed images in any way. We show that the monopole degeneracy and the mass-sheet degeneracy together lie at the heart of the discrepancies between different gravitational lens reconstructions that can be found in the literature for a given object, and that many images/sources, with an overall high image density in the lens plane, are required to construct an accurate, high-resolution mass map based on strong lensing data.     

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.