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.     

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.     

Microlensing by halo MACHOs with a spatially varying mass function

The main aim of microlensing experiments is to evaluate the mean mass of massive compact halo objects (MACHOs) and the mass fraction of the Galactic halo made by this type of dark matter. Statistical analysis shows that by considering a Dirac-Delta mass function (MF) for the MACHOs, their mean mass is about that of a white dwarf star. This result is, however, in discrepancy with other observations such as those of non-observed expected white dwarfs in the Galactic halo which give rise to metal abundance, polluting the interstellar medium by their evolution. Here we use the hypothesis of the spatially varying MF of MACHOs, proposed by Kerins and Evans to interpret microlensing events. In this model, massive lenses with a lower population contribute to the microlensing events more frequently than do dominant brown dwarfs. This effect causes the mean mass of the observed lenses to be larger than the mean mass of all the lenses. A likelihood analysis is performed to find the best parameters of the spatially varying MF that are compatible with the duration distribution of Large Magellanic Cloud microlensing candidates of the MACHO experiment.     

Microlensing in dark matter haloes

Using eight dark matter haloes extracted from fully self-consistent cosmological N -body simulations, we perform microlensing experiments. A hypothetical observer is placed at a distance of 8.5 kpc from the centre of the halo measuring optical depths, event durations and event rates towards the direction of the Large Magellanic Cloud. We simulate 1600 microlensing experiments for each halo. Assuming that the whole halo consists of massive astronomical compact halo objects (MACHOs),   f = 1.0  , and a single MACHO mass is   m _M= 1.0 M_⊙  , the simulations yield mean values of  τ= 4.7^+5.0_−2.2× 10⁻⁷  and  Γ= 1.6^+1.3_−0.6× 10⁻⁶  events star⁻¹ yr⁻¹. We find that triaxiality and substructure can have major effects on the measured values so that τ and Γ values of up to three times the mean can be found. If we fit our values of τ and Γ to the MACHO collaboration observations, we find   f = 0.23^+0.15_−0.13  and   m _M= 0.44^+0.24_−0.16  . Five out of the eight haloes under investigation produce f and m _M values mainly concentrated within these bounds.     

The observed concentration–mass relation for galaxy clusters

The properties of clusters of galaxies offer key insights into the assembly process of structure in the universe. Numerical simulations of cosmic structure formation in a hierarchical, dark matter dominated universe suggest that galaxy cluster concentrations, which are a measure of a halo's central density, decrease gradually with virial mass. However, cluster observations have yet to confirm this correlation. The slopes of the run of measured concentrations with virial mass are often either steeper or flatter than that predicted by simulations. In this work, we present the most complete sample of observed cluster concentrations and masses yet assembled, including new measurements for 10 strong-lensing clusters, thereby more than doubling the existing number of strong-lensing concentration estimates. We fit a power law to the observed concentrations as a function of virial mass, and find that the slope is consistent with the slopes found in simulations, though our normalization factor is higher. Observed lensing concentrations appear to be systematically larger than X-ray concentrations, a more pronounced effect than that found in simulations. We also find that at a fixed mass, the bulk of observed cluster concentrations are distributed lognormally, with the exception of a few anomalously high concentration clusters. We examine the physical processes likely responsible for the discrepancy between lensing and X-ray concentrations, and for the anomalously high concentrations in particular. The forthcoming Millennium simulation results will offer the most comprehensive comparison set to our findings of an observed concentration–mass power law relation.     

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.     

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.     

Microlensing of tidal debris on the Magellanic great circle

Increasing evidence suggests that the Galactic halo is lumpy on kpc scales as a result of the accretion of at least a dozen small galaxies [Large and Small Magellanic Clouds (LMC/SMC), Sgr, Fornax, etc.]. Faint stars in such lumpy structures can significantly microlense a background star with an optical depth of 10⁻⁷–10⁻⁶, which is comparable to the observed value to the LMC. The observed microlensing events towards the LMC can be explained by a tidal debris tail from the progenitor of the Magellanic Clouds and Magellanic Stream. The LMC stars can either lense stars in the debris tail a few kpc behind the LMC, or be lensed by stars in the part of the debris tail in front of the LMC. The models are consistent with an elementary particle dominated Galactic halo without massive compact halo objects (MACHOs). They also differ from Sahu's LMC-self-lensing model by predicting a higher optical depth and event rate and lower concentration of events to the LMC centre.     

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.     

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.     

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.