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.     

Weak gravitational lensing in the standard cold dark matter model, using an algorithm for three-dimensional shear

We investigate the effects of weak gravitational lensing in the standard cold dark matter cosmology, using an algorithm that evaluates the shear in three dimensions. The algorithm has the advantage of variable softening for the particles, and our method allows the appropriate angular diameter distances to be applied to every evaluation location within each three-dimensional simulation box. We investigate the importance of shear in the distance–redshift relation, and find it to be very small. We also establish clearly defined values for the smoothness parameter in the relation, finding its value to be at least 0.83 at all redshifts in our simulations. From our results, obtained by linking the simulation boxes back to source redshifts of 4, we are able to observe the formation of structure in terms of the computed shear, and also note that the major contributions to the shear come from a very broad range of redshifts. We show the probability distributions for the magnification, source ellipticity and convergence, and also describe the relationships amongst these quantities for a range of source redshifts. We find a broad range of magnifications and ellipticities; for sources at a redshift of 4, 97.5 per cent of all lines of sight show magnifications up to 1.39 and ellipticities up to 0.23. There is clear evidence that the magnification is not linear in the convergence, as might be expected for weak lensing, but contains contributions from higher order terms in both the convergence and the shear. Our results for the one-point distribution functions are generally different from those obtained by other authors using two-dimensional (planar) approaches, and we suggest reasons for the differences. Our magnification distributions for sources at redshifts of 1 and 0.5 are also very different from the results used by other authors to assess the effect on the perceived value of the deceleration parameter, and we briefly address this question.     

Rotation in gravitational lenses

Gravitational lensing deflects light. A single lens deflector can only shear images, but cannot induce rotations. Multiple lens planes can induce rotations. Such rotations can be observed in quadruply imaged sources, and can be used to distinguish between two proposed solutions of the flux anomaly problem: substructures in lensing galaxies versus large-scale structure. We predict the expected amount of rotation due to large-scale structure in strong lensing systems, and show how this effect can be measured using ∼mas very long baseline interferometry astrometry of quadruple lenses with extended source structures. The magnitude of rotation is around 1°. The biggest theoretical uncertainty is the power spectrum of dark matter on very small scales. This procedure can potentially be turned around to measure the dark matter power spectrum on very small scales. We list the predicted rms rotation angles for several quadruple lenses with known lens and source redshifts.     

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

Source–lens clustering effects on the skewness of the lensing convergence

The correlation between source galaxies and lensing potentials causes a systematic effect on measurements of cosmic shear statistics, known as the source–lens clustering (SLC) effect. The SLC effect on the skewness of lensing convergence, S ₃, is examined using a non-linear semi-analytic approach and is checked against numerical simulations. The semi-analytic calculations have been performed in a wide variety of generic models for the redshift distribution of source galaxies and power-law models for the bias parameter between the galaxy and dark matter distributions. The semi-analytic predictions are tested successfully against numerical simulations. We find the relative amplitude of the SLC effect on S ₃ to be of the order of  5–40  per cent. It depends significantly on the redshift distribution of sources and on the way in which the bias parameter evolves. We discuss possible measurement strategies to minimize the SLC effects.     

FAst STatistics for weak Lensing (FASTLens): fast method for weak lensing statistics and map making

With increasingly large data sets, weak lensing measurements are able to measure cosmological parameters with ever-greater precision. However, this increased accuracy also places greater demands on the statistical tools used to extract the available information. To date, the majority of lensing analyses use the two-point statistics of the cosmic shear field. These can be either studied directly using the two-point correlation function or in Fourier space, using the power spectrum. But analysing weak lensing data inevitably involves the masking out of regions, for example to remove bright stars from the field. Masking out the stars is common practice but the gaps in the data need proper handling. In this paper, we show how an inpainting technique allows us to properly fill in these gaps with only   N log  N   operations, leading to a new image from which we can compute straightforwardly and with a very good accuracy both the power spectrum and the bispectrum. We then propose a new method to compute the bispectrum with a polar fft algorithm, which has the main advantage of avoiding any interpolation in the Fourier domain. Finally, we propose a new method for dark matter mass map reconstruction from shear observations, which integrates this new inpainting concept. A range of examples based on 3D N -body simulations illustrates the results.     

Statistics of weak lensing at small angular scales: analytical predictions for lower order moments

Weak lensing surveys are expected to provide direct measurements of the statistics of the projected dark matter distribution. Most analytical studies of weak lensing statistics have been limited to quasi-linear scales as they relied on perturbative calculations. On the other hand, observational surveys are likely to probe angular scales less than 10 arcmin, for which the relevant physical length-scales are in the non-linear regime of gravitational clustering. We use the hierarchical ansatz to compute the multipoint statistics of the weak lensing convergence for these small smoothing angles. We predict the multipoint cumulants and cumulant correlators up to fourth order and compare our results with high-resolution ray-tracing simulations. Averaging over a large number of simulation realizations for four different cosmological models, we find close agreement with the analytical calculations. In combination with our work on the probability distribution function, these results provide accurate analytical models for the full range of weak lensing statistics. The models allow for a detailed exploration of cosmological parameter space and of the dependence on angular scale and the redshift distribution of source galaxies. We compute the dependence of the higher moments of the convergence on the parameters Ω and Λ.     

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.     

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.     

Weak gravitational lensing by voids

We consider the prospects for detecting weak gravitational lensing by underdensities (voids) in the large-scale matter distribution. We derive the basic expressions for magnification and distortion by spherical voids. Clustering of the background sources and cosmic variance are the main factors that limit in principle the detection of lensing by voids. We conclude that only voids with radii larger than ∼100  h ⁻¹ Mpc have lensing signal-to-noise ratio larger than unity.     

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.     

MOND plus classical neutrinos are not enough for cluster lensing

Clusters of galaxies offer a robust test bed for probing the nature of dark matter that is insensitive to the assumption of the gravity theories. Both Modified Newtonian Dynamics (MOND) and General Relativity (GR) would require similar amounts of non-baryonic matter in clusters as MOND boosts the gravity only mildly on cluster scales. Gravitational lensing allows us to estimate the enclosed mass in clusters on small (∼20–50 kpc) and large (∼several 100 kpc) scales independent of the assumptions of equilibrium. Here, we show for the first time that a combination of strong and weak gravitational lensing effects can set interesting limits on the phase-space density of dark matter in the centres of clusters. The phase-space densities derived from lensing observations are inconsistent with neutrino masses ranging from 2–7 eV, and hence do not support the 2 eV-range particles required by MOND. To survive, the most plausible modification for MOND may be an additional degree of dynamical freedom in a covariant incarnation.     

Designing weak lensing surveys: a generalized eigenmode analysis

We study the estimators of various second-order weak lensing statistics such as the shear correlation functions  ξ_±  and the aperture mass dispersion  〈 M ²_ap〉  which can directly be constructed from weak lensing shear maps. We compare the efficiency with which these estimators can be used to constrain cosmological parameters. To this end we introduce the Karhunen–Loève (KL) eigenmode analysis techniques for weak lensing surveys. These tools are shown to be very effective as a diagnostics for optimizing survey strategies. The usefulness of these tools to study the effect of angular binning, the depth and width of the survey and noise contributions due to intrinsic ellipticities and number density of source galaxies on the estimation of cosmological parameters is demonstrated. Results from independent analysis of various parameters and joint estimations are compared. We also study how degeneracies among various cosmological and survey parameters affect the eigenmodes associated with these parameters.     

Angular cross-correlation of galaxies: a probe of gravitational lensing by large-scale structure

The angular cross-correlation between two galaxy samples separated in redshift is shown to be a useful measure of weak lensing by large-scale structure. Angular correlations in faint galaxies arise as a result of spatial clustering of the galaxies as well as gravitational lensing by dark matter along the line of sight. The lensing contribution to the two-point autocorrelation function is typically small compared with the gravitational clustering. However, the cross-correlation between two galaxy samples is almost unaffected by gravitational clustering provided that their redshift distributions do not overlap. The cross-correlation is then induced by magnification bias resulting from lensing by large-scale structure. We compute the expected amplitude of the cross-correlation for popular theoretical models of structure formation. For two populations with mean redshifts of ≃0.3 and 1, we find a cross-correlation signal of ≃1 per cent on arcmin scales and ≃3 per cent on scales of a few arcsec. The dependence on the cosmological parameters Ω and Λ, the dark matter power spectrum and the bias factor of the foreground galaxy population is explored.     

The continuous limit of the multiple lens effect and the optical scalar equation

We study the continuous limit of the multiple gravitational lensing theory based on the thin lens approximation. We define a new, light-path dependent angular diameter distance     and show that it satisfies the optical scalar equation. The distance provides relations between quantities used in gravitational lensing theory (the convergence, the shear and the twist terms) and those used in scalar optics theory (the rates of expansion, shear and rotation).     

Weak gravitational lensing in different cosmologies, using an algorithm for shear in three dimensions

We present the results of weak gravitational lensing statistics in four different cosmological N -body simulations. The data have been generated using an algorithm for the three-dimensional shear, which makes use of a variable softening facility for the N -body particle masses, and enables a physical interpretation for the large-scale structure to be made. Working in three dimensions also allows the correct use of the appropriate angular diameter distances.
Our results are presented on the basis of the filled-beam approximation in view of the variable particle softening scheme in our algorithm. The importance of the smoothness of matter in the Universe for the weak lensing results is discussed in some detail.
The low-density cosmology with a cosmological constant appears to give the broadest distributions for all the statistics computed for sources at high redshifts. In particular, the range in magnification values for this cosmology has implications for the determination of the cosmological parameters from high-redshift type Ia supernovae. The possibility of determining the density parameter from the non-Gaussianity in the probability distribution for the convergence is discussed.     

Testing dark energy with the Advanced Liquid-mirror Probe of Asteroids, Cosmology and Astrophysics

The Advanced Liquid-mirror Probe of Asteroids, Cosmology and Astrophysics (ALPACA) is a proposed 8-m liquid-mirror telescope surveying  ∼1000 deg²  of the Southern hemisphere sky. It will be a remarkably simple and inexpensive telescope that none the less will deliver a powerful sample of optical data for studying dark energy. The bulk of the cosmological data consist of nightly, high signal-to-noise ratio, multiband light curves of Type Ia supernovae (SNe Ia). At the end of the 3-yr run, ALPACA is expected to collect  ≳100 000  SNe Ia up to   z ∼ 1  . This will allow us to reduce present systematic uncertainties affecting the standard-candle relation. The survey will also provide several other data sets such as the detection of baryon acoustic oscillations in the matter power spectrum and shear weak-lensing measurements. In this preliminary analysis, we forecast constraints on dark energy parameters from SNe Ia and baryon acoustic oscillations. The combination of these two data sets will provide competitive constraints on the dark energy parameters under minimal prior assumptions. Further studies are needed to address the accuracy of weak-lensing measurements.     

A Remark on Using Gravitational Lensing Probability as a Probe of the Central Regions of CDM Halos

We calculate the gravitational lensing probabilities by cold dark matter (CDM) halos with different density profiles, and compare them with current observations from the Cosmic Lens All-Sky Survey (CLASS) and the Jodrell-Bank VLA Astrometric Survey (JVAS). We find that the lensing probability is dramatically sensitive to the clumping of the dark matter, or quantitatively, the concentration parameter. We also find that our predicted lensing probabilities in most cases show inconsistency with the observations. It is argued that high lensing probability may not be an effective tool for probing the statistical properties of inner structures of dark matter halos.     

A new quadruple gravitational lens system: CLASS B0128+437

High-resolution MERLIN observations of a newly discovered four-image gravitational lens system, B0128+437, are presented. The system was found after a careful re-analysis of the entire CLASS data set. The MERLIN observations resolve four components in a characteristic quadruple-image configuration; the maximum image separation is 542 mas and the total flux density is 48 mJy at 5 GHz. A best-fitting lens model with a singular isothermal ellipsoid results in large errors in the image positions. A significantly improved fit is obtained after the addition of a shear component, suggesting that the lensing system is more complex and may consist of multiple deflectors. The integrated radio spectrum of the background source indicates that it is a gigahertz peaked spectrum source. It may therefore be possible to resolve structure within the radio images with deep VLBI observations and thus to constrain the lensing mass distribution better.