Noise properties of gravitational lens mass reconstruction

Gravitational lensing is potentially able to observe mass-selected haloes, and to measure the projected cluster mass function. An optimal mass selection requires a quantitative understanding of the noise behaviour in mass maps. This paper is an analysis of the noise properties in mass maps reconstructed from a maximum-likelihood method.
The first part of this work is the derivation of the noise power spectrum and the mass error bars as a straightforward extension of the Kaiser & Squires algorithm for the case of a correlated noise. Very good agreement is found between these calculations and the noise properties measured in the mass reconstructions limited to non-critical clusters of galaxies. It demonstrates that Kaiser & Squires and maximum-likelihood methods have similar noise properties and that the weak lensing approximation is valid for describing these properties .
In a second stage I show that the statistics of peaks in the noise follows accurately the peak statistics of a two-dimensional Gaussian random field (using the BBKS techniques) if the smoothing aperture contains enough galaxies. This analysis provides a full procedure for deriving the significance of any convergence peak as a function of its amplitude and profile.
I demonstrate that a detailed quantitative analysis of the structures in mass maps can be carried out, and that, to a very good approximation, a mass map is the sum of the lensing signal and known two-dimensional Gaussian random noise. A straightforward application is the measurement of the projected mass function in wide-field lensing surveys, down to small mass overdensities that are individually undetectable.     

Non-parametric reconstruction of cluster mass distribution from strong lensing: modelling Abell 370

We describe a new non-parametric technique for reconstructing the mass distribution in galaxy clusters with strong lensing, i.e. from multiple images of background galaxies. The observed positions and redshifts of the images are considered as rigid constraints, and through the lens (ray-trace) equation they provide us with linear constraint equations. These constraints confine the mass distribution to some allowed region, which is then found by linear programming. Within this allowed region we study in detail the mass distribution with minimum mass-to-light variation, and also some other distributions, such as the smoothest mass distribution.
The method is applied to the extensively studied cluster Abell 370, which hosts a giant luminous arc and several other multiply imaged background galaxies. Our mass maps are constrained by the observed positions and redshifts (spectroscopic or model-inferred by previous authors) of the giant arc and multiple-image systems. The reconstructed maps obtained for Abell 370 reveal a detailed mass distribution, with substructure quite different from the light distribution. The method predicts the bimodal nature of the cluster, and that the projected mass distribution is indeed elongated along the axis defined by the two dominant cD galaxies. However, the peaks in the mass distribution appear to be offset from the centres of the cDs.
We also present an estimate for the total mass of the central region of the cluster. This is in good agreement with previous mass determinations. The total mass of the central region is M =(2.0–2.7)×10¹⁴ M⊙ h ⁻¹₅₀, depending on the solution chosen.     

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.     

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.     

Gravitational lens magnification by Abell 1689: distortion of the background galaxy luminosity function

Gravitational lensing magnifies the observed flux of galaxies behind the lens. We use this effect to constrain the total mass in the cluster Abell 1689 by comparing the lensed luminosities of background galaxies with the luminosity function of an undistorted field. Under the assumption that these galaxies are a random sample of luminosity space, this method is not limited by clustering noise. We use photometric redshift information to estimate galaxy distance and intrinsic luminosity. Knowing the redshift distribution of the background population allows us to lift the mass/background degeneracy common to lensing analysis. In this paper we use nine filters observed over 12 h with the Calar Alto 3.5-m telescope to determine the redshifts of 1000 galaxies in the field of Abell 1689. Using a complete sample of 146 background galaxies we measure the cluster mass profile. We find that the total projected mass interior to 0.25  h ⁻¹ Mpc is M _2D(<0.25  h ⁻¹ Mpc)=(0.48±0.16)×10¹⁵  h ⁻¹ M_⊙, where our error budget includes uncertainties from the photometric redshift determination, the uncertainty in the offset calibration and finite sampling. This result is in good agreement with that found by number-count and shear-based methods and provides a new and independent method to determine cluster masses.     

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.     

Quasar–galaxy and galaxy–galaxy cross-correlations: model predictions with realistic galaxies

Several measurements of quasi-stellar object (QSO)–galaxy correlations have reported signals much larger than predictions of magnification by large-scale structure. We find that the expected signal depends strongly on the properties of the foreground galaxy population. On arcmin scales, it can be either larger or smaller by a factor of 2 for different galaxy types in comparison with a linearly biased version of the mass distribution. Thus the resolution of some of the excess measurements may lie in examining the halo occupation properties of the galaxy population sampled by a given survey; this is also the primary information such measurements will provide.
We use the halo model of clustering and simulations to predict the magnification-induced cross-correlations and errors for forthcoming surveys. With the full Sloan Digital Sky Survey, the statistical errors will be below 1 per cent for the galaxy–galaxy correlations and significantly larger for QSO–galaxy correlations. Thus accurate constraints on parameters of the galaxy halo occupation distribution can be obtained from small-scale measurements and on the bias parameter from large scales. Since the lensing-induced cross-correlation measures the first moment of the halo occupation number of galaxies, these measurements can provide the basis for interpreting galaxy clustering measurements that measure the second- and higher-order moments.     

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.     

Infrared observations of gravitational lensing in Abell 2219 with CIRSI

We present the first detection of a gravitational depletion signal at near-infrared wavelengths, based on deep panoramic images of the cluster Abell 2219 ( z =0.22) taken with the Cambridge Infrared Survey Instrument (CIRSI) at the prime focus of the 4.2-m William Herschel Telescope. Infrared studies of gravitational depletion offer a number of advantages over similar techniques applied at optical wavelengths, and can provide reliable total masses for intermediate-redshift clusters. Using the maximum-likelihood technique developed by Schneider, King & Erben, we detect the gravitational depletion at the 3 confidence level. By modelling the mass distribution as a singular isothermal sphere and ignoring the uncertainty in the unlensed number counts, we find an Einstein radius of (66 per cent confidence limit). This corresponds to a projected velocity dispersion of _v 800 km s¹, in agreement with constraints from strongly lensed features. For a Navarro, Frenk & White mass model, the radial dependence observed indicates a best-fitting halo scalelength of 125 h ¹ kpc. We investigate the uncertainties arising from the observed fluctuations in the unlensed number counts, and show that clustering is the dominant source of error. We extend the maximum-likelihood method to include the effect of incompleteness, and discuss the prospects of further systematic studies of lensing in the near-infrared band.     

Giant cluster arcs as a constraint on the scattering cross-section of dark matter

We carry out ray tracing through five high-resolution simulations of a galaxy cluster, to study how its ability to produce giant gravitationally lensed arcs is influenced by the collision cross-section of its dark matter. In three cases typical dark matter particles in the cluster core undergo between 1 and 100 collisions per Hubble time; two more explore the long ('collisionless') and short ('fluid') mean free path limits. We study the size and shape distributions of arcs and compute the cross-section for producing 'extreme' arcs of various sizes. Even a few collisions per particle modifies the core structure enough to destroy the ability of the cluster to produce long, thin arcs. For larger collision frequencies the cluster must be scaled up to unrealistically large masses before it regains the ability to produce giant arcs. None of our models with self-interacting dark matter (except the 'fluid' limit) is able to produce radial arcs; even the case with the smallest scattering cross-section must be scaled to the upper limit of observed cluster masses before it produces radial arcs. Apparently the elastic collision cross-section of dark matter in clusters must be very small, below 0.1 cm² g⁻¹, to be compatible with the observed ability of clusters to produce both radial arcs and giant arcs.     

Non-parametric inversion of gravitational lensing systems with few images using a multi-objective genetic algorithm

Galaxies acting as gravitational lenses are surrounded by, at most, a handful of images. This apparent paucity of information forces one to make the best possible use of what information is available to invert the lens system. In this paper, we explore the use of a genetic algorithm to invert in a non-parametric way strong lensing systems containing only a small number of images. Perhaps the most important conclusion of this paper is that it is possible to infer the mass distribution of such gravitational lens systems using a non-parametric technique. We show that including information about the null space (i.e. the region where no images are found) is prerequisite to avoid the prediction of a large number of spurious images, and to reliably reconstruct the lens mass density. While the total mass of the lens is usually constrained within a few per cent, the fidelity of the reconstruction of the lens mass distribution depends on the number and position of the images. The technique employed to include null space information can be extended in a straightforward way to add additional constraints, such as weak-lensing data or time-delay information.     

Effects of cluster galaxies on arc statistics

We present the results of a set of numerical simulations evaluating the effect of cluster galaxies on arc statistics.
We perform a first set of gravitational lensing simulations using three independent projections for each of nine different galaxy clusters obtained from N -body simulations. The simulated clusters consist of dark matter only. We add a population of galaxies to each cluster, mimicking the observed luminosity function and the spatial galaxy distribution, and repeat the lensing simulations including the effects of cluster galaxies, which themselves act as individual lenses. Each galaxy is represented by a spherical Navarro, Frenk & White density profile.
We consider the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of lengths, widths, curvature radii and length-to-width ratios of long arcs. We find some changes to the distribution of short-arc properties in the presence of cluster galaxies. The differences appear in the distribution of curvature radii for arc lengths smaller than 12 arcsec, while the distributions of lengths, widths and length-to-width ratios are significantly changed only for arcs shorter than 4 arcsec.     

The JVAS/CLASS search for 6-arcsec to 15-arcsec image separation lensing

The Jodrell Bank–VLA Astrometric Survey (JVAS) and the Cosmic Lens All Sky Survey (CLASS) have been systematically searched for multiple gravitational imaging of sources with image separations between 6 arcsec and 15 arcsec, associated with galaxy group and cluster lensing masses. The radio and optical follow-up observations of all candidates are presented. From a total of ∼15 000 sources only one weak candidate remains and this is not contained in the statistically complete sample of flat-spectrum JVAS/CLASS sources of 11 670 sources. A simple Press–Schechter analysis is performed. For singular isothermal sphere lenses the lack of multiple image systems is inconsistent with the currently favoured cosmologies with     at the 4.2 σ level. Cored isothermal lenses reduce the expected number of lens systems and we suggest that the most probable interpretation of our results is that the surface mass density of groups and clusters of galaxies is not high enough to cause multiple imaging and the presence of the mass concentrations associated with individual galaxies is required to produce image separations such as those in B0957+561.     

Chandra observations of the galaxy cluster Abell 1835

We present the analysis of 30 ks of Chandra observations of the galaxy cluster Abell 1835. Overall, the X-ray image shows a relaxed morphology, although we detect substructure in the inner 30-kpc radius. Spectral analysis shows a steep drop in the X-ray gas temperature from ∼12 keV in the outer regions of the cluster to ∼4 keV in the core. The Chandra data provide tight constraints on the gravitational potential of the cluster which can be parametrized by a Navarro, Frenk & White model. The X-ray data allow us to measure the X-ray gas mass fraction as a function of radius, leading to a determination of the cosmic matter density of
   
. The projected mass within a radius of ∼150 kpc implied by the presence of gravitationally lensed arcs in the cluster is in good agreement with the mass models preferred by the Chandra data. We find a radiative cooling time of the X-ray gas in the centre of Abell 1835 of about
   
. Cooling-flow model fits to the Chandra spectrum and a deprojection analysis of the Chandra image both indicate the presence of a young cooling flow (∼     with an integrated mass deposition rate of     within a radius of 30 kpc. We discuss the implications of our results in the light of recent Reflection Grating Spectrograph (RGS) observations of Abell 1835 with XMM-Newton .