A genetic algorithm for the non-parametric inversion of strong lensing systems

We present a non-parametric technique to infer the projected mass distribution of a gravitational lens system with multiple strong-lensed images. The technique involves a dynamic grid in the lens plane on which the mass distribution of the lens is approximated by a sum of basis functions, one per grid cell. We used the projected mass densities of Plummer spheres as basis functions. A genetic algorithm then determines the mass distribution of the lens by forcing images of a single source, projected back on to the source plane, to coincide as well as possible. Averaging several tens of solutions removes the random fluctuations that are introduced by the reproduction process of genomes in the genetic algorithm and highlights those features common to all solutions. Given the positions of the images and the redshifts of the sources and the lens, we show that the mass of a gravitational lens can be retrieved with an accuracy of a few percent and that, if the sources sufficiently cover the caustics, the mass distribution of the gravitational lens can also be reliably retrieved. A major advantage of the algorithm is that it makes full use of the information contained in the radial images, unlike methods that minimize the residuals of the lens equation, and is thus able to accurately reconstruct also the inner parts of the lens.     

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⁻⁴.     

ZEN2: a narrow J-band search for z∼ 9 Lyα emitting galaxies directed towards three lensing clusters

We present the results of a continuing survey to detect Lyα emitting galaxies at redshifts   z ∼ 9  : the ' z equals nine' (ZEN) survey. We have obtained deep VLT Infrared Spectrometer and Array Camera observations in the narrow J -band filter NB119 directed towards three massive lensing clusters: Abell clusters 1689, 1835 and 114. The foreground clusters provide a magnified view of the distant Universe and permit a sensitive test for the presence of very high redshift galaxies. We search for   z ∼ 9 Lyα  emitting galaxies displaying a significant narrow-band excess relative to accompanying J -band observations that remain undetected in Hubble Space Telescope ( HST )/Advanced Camera for Surveys (ACS) optical images of each field. No sources consistent with this criterion are detected above the unlensed 90 per cent point-source flux limit of the narrow-band image,   F _NB= 3.7 × 10⁻¹⁸ erg s⁻¹ cm⁻²  . To date, the total coverage of the ZEN survey has sampled a volume at   z ∼ 9  of approximately 1700 comoving Mpc³ to a Lyα emission luminosity of  10⁴³ erg s⁻¹  . We conclude by considering the prospects for detecting   z ∼ 9 Lyα  emitting galaxies in light of both observed galaxy properties at   z < 7  and simulated populations at   z > 7  .     

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.     

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.     

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.     

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.     

A search for the third lensed image in JVAS B1030+074

Central gravitational image detection is very important for the study of the mass distribution of the inner parts (∼100 pc) of lens galaxies. However, the detection of such images is extremely rare and difficult. We present a 1.7-GHz High Sensitivity Array (HSA) observation of the double-image radio lens system B1030+074. The data are combined with archive Very Long Baseline Array and global very long baseline interferometry (VLBI) observations, and careful consideration is given to the effects of noise, clean ing and self-calibration. An upper limit is derived for the strength of the central image of 180 μJy (90 per cent confidence level), considerably greater than would have been expected on the basis of a simple analysis. This gives a lower limit of ∼10³ for the ratio of the brightest image to the central image. For cusped models of lens mass distributions, we have made use of this non-detection to constrain the relation between inner power-law slope β of the lensing galaxy mass profile, and its break radius r _b. For   r _b > 130 pc  the power-law slope is required to be close to isothermal  (β > 1.8)  . A flatter inner slope is allowed if a massive black hole is present at the centre of the lensing galaxy, but the effect of the black hole is small unless it is ∼10 times more massive than that implied by the relation between black hole mass and stellar velocity dispersion. By comparing four epochs of VLBI observations, we also detected possible superluminal motion in the jet in the brighter image A. The B jet remains unresolved, as expected from a simple lens model of the system.