A maximum-entropy method for reconstructing the projected mass distribution of gravitational lenses

The maximum-entropy method is applied to the problem of reconstructing the projected mass density of a galaxy cluster using its gravitational lensing effects on background galaxies. We demonstrate the method by reconstructing the mass distribution in a model cluster using simulated shear and magnification data to which Gaussian noise is added. The mass distribution is reconstructed directly and the inversion is regularized using the entropic prior for this positive additive distribution. For realistic noise levels, we find that the method faithfully reproduces the main features of the cluster mass distribution not only within the observed field but also slightly beyond it. We estimate the uncertainties in the reconstruction by calculating an analytic approximation to the covariance matrix of the reconstruction values of each pixel. This result is compared with error estimates derived from Monte Carlo simulations for different noise realizations and found to be in good agreement.     

Spectroscopic confirmation of redshifts predicted by gravitational lensing

We present deep spectroscopic measurements of 18 distant field galaxies identified as gravitationally lensed arcs in a Hubble Space Telescope image of the cluster Abell 2218. Redshifts of these objects were predicted by Kneib et al. using a lensing analysis constrained by the properties of two bright arcs of known redshift and other multiply imaged sources. The new spectroscopic identifications were obtained using long exposures with the LDSS-2 spectrograph on the William Herschel Telescope, and demonstrate the capability of that instrument to reach new limits, R ≃24; the lensing magnification implies true source magnitudes as faint as R ≃25. Statistically, our measured redshifts are in excellent agreement with those predicted from Kneib et al.'s lensing analysis, and this gives considerable support to the redshift distribution derived by the lensing inversion method for the more numerous and fainter arclets extending to R ≃25.5. We explore the remaining uncertainties arising from both the mass distribution in the central regions of Abell 2218 and the inversion method itself, and conclude that the mean redshift of the faint field population at R ≃25.5 ( B ∼26–27) is low, 〈 z 〉=0.8–1. We discuss this result in the context of redshift distributions estimated from multicolour photometry. Although such comparisons are not straightforward, we suggest that photometric techniques may achieve a reasonable level of agreement, particularly when they include near-infrared photometry with discriminatory capabilities in the 1< z <2 range.     

Resolving the Mass Discrepancy between Strong and Weak Lensing Methods

We readdress the outstanding cluster mass discrepancy between strong and weak gravitational lensing techniques utilizing updated data of both giant arcs and weak lensing measurements from the literature.We find that the systematically higher values of cluster masses revealed by strong lensing can be attributed to the oversimplification of the lensing model when estimating the cluster mass enclosed within the giant arcs.This arises because inhomogeneities and substructures in the central cores of clusters may invalidate the spherical symmetry assumption used widely in previous applications.When a more realistic modeling of the arcs is used,then the masses by strong lensing agree fairly well with those given by weak lensing when both are extrapolated to the same cluster regions.We conclude that as of now no significant discrepancy has been found among different cluster mass estimators including optical galaxies,X-ray gas and lensing.     

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.     

Gravitational lensing of type Ia supernovae by galaxy clusters

We propose a method to remove the mass-sheet degeneracy that arises when the mass of galaxy clusters is inferred from gravitational shear. The method utilizes high-redshift standard candles that undergo weak lensing. Natural candidates for such standard candles are type Ia supernovae (SNe Ia).
When corrected with the light-curve shape (LCS), the peak magnitude of SNe Ia provides a standard candle with an uncertainty in apparent magnitude of Δ m ≃0.1–0.2. Gravitational magnification of a background SN Ia by an intervening cluster would cause a mismatch between the observed SN Ia peak magnitude compared with that expected from its LCS and redshift. The average detection rate for SNe Ia with a significant mismatch of ≥2Δ m behind a cluster at z ≃0.05–0.15 is about 1–2 supernovae per cluster per year at J , I , R ≲25–26.
Since SNe are point-like sources for a limited period, they can experience significant microlensing by massive compact halo objects (MACHOs) in the intracluster medium. Microlensing events caused by MACHOs of ∼10⁻⁴ M⊙ are expected to have time-scales similar to that of the SN light curve. Both the magnification curve by a MACHO and the light curve of a SN Ia have characteristic shapes that allow us to separate them. Microlensing events caused by MACHOs of smaller mass can unambiguously be identified in the SN light curve if the latter is continuously monitored. The average number of identifiable microlensing events per nearby cluster ( z ≲0.05) per year is ∼0.02 ( f /0.01), where f is the fraction of the cluster mass in MACHOs of masses 10⁻⁷< M _macho/M⊙<10⁻⁴.     

Lensing magnification effects on the cosmic shear statistics

Gravitational lensing causes a correlation between a population of foreground large-scale structures and the observed number density of the background distant galaxies as a consequence of the flux magnification and the lensing area distortion. This correlation has not been taken into account in calculations of the theoretical predictions of the cosmic shear statistics but may cause a systematic error in a cosmic shear measurement. We examine its impact on the cosmic shear statistics using the semi-analytic approach. We find that the lensing magnification has no practical influence on the cosmic shear variance. Exploring the possible shapes of the redshift distribution of source galaxies, we find that the relative amplitude of the effect on the convergence skewness is 3 per cent at most.     

QSO lensing magnification associated with galaxy groups

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

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

Detecting mass substructure in galaxy clusters: an aperture mass statistic for gravitational flexion

Gravitational flexion has been introduced as a technique by which one can map out and study substructure in clusters of galaxies. Previous analyses involving flexion have measured the individual galaxy–galaxy flexion signal, or used either parametric techniques or a Kaiser, Squires and Broadhurst (KSB)-type inversion to reconstruct the mass distribution in Abell 1689. In this paper, we present an aperture mass statistic for flexion, and apply it to the lensed images of background galaxies obtained by ray-tracing simulations through a simple analytic mass distribution and through a galaxy cluster from the Millennium Simulation. We show that this method is effective at detecting and accurately tracing structure within clusters of galaxies on subarcminute scales with high signal to noise even using a moderate background source number density and image resolution. In addition, the method provides much more information about both the overall shape and the small-scale structure of a cluster of galaxies than can be achieved through a weak lensing mass reconstruction using gravitational shear data. Lastly, we discuss how the zero-points of the aperture mass might be used to infer the masses of structures identified using this method.     

Investigating cosmological weak lensing with the ray-bundle method

Using the ray-bundle method for calculating gravitational lens magnifications, we outline a method by which the magnification probability may be determined specifically in the weak lensing limit for cosmological models obtained from N -body simulations.
16 different models are investigated, which are variations on three broad classes of cold dark matter model: the standard model with  (Ω₀, λ ₀)=(1.0,0.0)  , the open model with  (Ω₀, λ ₀)=(0.3,0.0)  and the lambda model, which is a flat model with a cosmological constant  (Ω₀, λ ₀)=(0.3,0.7)  .
The effects of varying the Hubble parameter, H ₀, the power spectrum shape parameter, Γ, and the cluster mass normalization, σ ₈, are studied. It is shown that there is no signature of these parameters in the weak lensing magnification distributions. The magnification probability distributions are also shown to be independent of the numerical parameters such as the lens mass and simulation box size in the N -body simulations.     

X-ray and strong lensing mass estimate of MS2137.3–2353

We present new mass estimates of the cluster of galaxies MS2137.3–2353, inferred from X-ray and strong lensing analyses. This cluster exhibits an outstanding strong lensing configuration and indicates a well-relaxed dynamical state, being most suitable for a mass reconstruction which combines both techniques. Despite this, several previous studies have claimed a significant discrepancy between the X-ray and the strong lensing mass estimates. The primary aim of this paper is to address and explain this mismatch. For this purpose, we have analysed Chandra observations to recover the profiles of the intracluster medium properties and, assuming a functional form for the matter density, the total mass distribution. The notable strong-lensing features of MS2137.3 allow us to reconstruct its projected mass in the central regions with good accuracy, by taking advantage of the lensing inversion code lenstool . We compare the results obtained for both methods. Our mass estimates for MS2137.3 are in agreement within errors, leading to a mean, extrapolated value of   M ₂₀₀≃ 4.4 ± 0.3 × 10¹⁴ M_⊙  , under the assumption of the Navarro–Frenk–White (NFW) mass profile. However, the strong lensing mass estimate is affected by the details of the brightest cluster galaxy mass modelling, since the radial arc is a very sensitive probe of the total mass derivative in the central region. In particular, we do not find evidence for a high concentration for the NFW density profile, as reported in some earlier works.     

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.     

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.     

Galaxy density profiles and shapes – II. Selection biases in strong lensing surveys

Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio and image separation distribution, with a flexible treatment of magnification bias to mimic different survey strategies. We present our results for two families of density profiles: cusped and deprojected Sérsic models. While we use fixed lens and source redshifts for most of the analysis, we show that the results are applicable to other redshift combinations, and we also explore the physics of how our results change for very different redshifts. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sérsic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.     

Kinematic effect in gravitational lensing by clusters of galaxies

Gravitational lensing provides an efficient tool for the investigation of matter structures, independent of the dynamical or the hydrostatic equilibrium properties of the deflecting system. However, it depends on the kinematic status. In fact, either a translational motion or a coherent rotation of the mass distribution can affect the lensing properties. Here, light deflection by galaxy clusters in motion is considered. Even if gravitational lensing mass measurements of galaxy clusters are regarded as very reliable estimates, the kinematic effect should be considered. A typical peculiar motion with respect to the Hubble flow brings about a systematic error ≲0.3 per cent, independent of the mass of the cluster. On the other hand, the effect of the spin increases with the total mass. For cluster masses  ∼10¹⁵ M_⊙  , the effect of the gravitomagnetic term is ≲0.04 per cent on strong lensing estimates and ≲0.5 per cent in the weak-lensing analyses. The total kinematic effect on the mass estimate is then ≲1 per cent, which is negligible in current statistical studies. In the weak-lensing regime, the rotation imprints a typical angular modulation in the tangential shear distortion. This would allow, in principle, a detection of the gravitomagnetic field and a direct measurement of the angular velocity of the cluster but the required background source densities are well beyond current technological capabilities.     

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.     

Infrared constraints on the dark mass concentration observed in the cluster Abell 1942

We present a deep H -band image of the region in the vicinity of the cluster Abell 1942 containing the puzzling dark matter concentration detected in an optical weak lensing study by Erben et al. We demonstrate that our limiting magnitude, H =22 , would be sufficient to detect clusters of appropriate mass out to redshifts comparable with the mean redshift of the background sources. Despite this, our infrared image reveals no obvious overdensity of sources at the location of the lensing mass peak, nor an excess of sources in the I − H versus H colour–magnitude diagram. We use this to constrain further the luminosity and mass-to-light ratio of the putative dark clump as a function of its redshift. We find that for spatially flat cosmologies, background lensing clusters with reasonable mass-to-light ratios lying in the redshift range 0< z <1 are strongly excluded, leaving open the possibility that the mass concentration is a new type of truly dark object.     

A comparison of different cluster mass estimates: consistency or discrepancy?

Rich and massive clusters of galaxies at intermediate redshift are capable of magnifying and distorting the images of background galaxies. A comparison of different mass estimators among these clusters can provide useful information about the distribution and composition of cluster matter and its dynamical evolution. Using the hitherto largest sample of lensing clusters drawn from the literature, we compare the gravitating masses of clusters derived from the strong/weak gravitational lensing phenomena, from the X-ray measurements based on the assumption of hydrostatic equilibrium, and from the conventional isothermal sphere model for the dark matter profile characterized by the velocity dispersion and core radius of galaxy distributions in clusters. While there is excellent agreement between the weak lensing, X-ray and isothermal sphere model-determined cluster masses, these methods are likely to underestimate the gravitating masses enclosed within the central cores of clusters by a factor of 2–4 as compared with the strong lensing results. Such a mass discrepancy has probably arisen from the inappropriate applications of the weak lensing technique and the hydrostatic equilibrium hypothesis to the central regions of clusters, as well as from assuming an unreasonably large core radius for both luminous and dark matter profiles. Nevertheless, it is pointed out that these cluster mass estimators may be safely applied on scales greater than the core sizes. Namely, the overall clusters of galaxies at intermediate redshift can still be regarded as the dynamically relaxed systems, in which the velocity dispersion of galaxies and the temperature of X-ray emitting gas are good indicators of the underlying gravitational potentials of clusters.     

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.