VLA 8.4-GHz monitoring observations of the CLASS gravitational lens B1933+503

The complex 10-component gravitational lens system B1933+503 has been monitored with the Very Large Array (VLA) during the period 1998 February–June with a view to measuring the time delay between the four compact components and hence to determine the Hubble parameter H ₀. Here we present the results of an 'A' configuration 8.4-GHz monitoring campaign that consists of 37 epochs with an average spacing of 2.8 d. The data have yielded light curves for the four flat-spectrum radio components (components 1, 3, 4 and 6). We observe only small flux density changes in the four flat-spectrum components that we do not believe are predominantly intrinsic to the source. Therefore the variations do not allow us to determine the independent time delays in this system. However, the data do allow us to accurately determine the flux density ratios between the four flat-spectrum components. These will prove important as modelling constraints and could prove crucial in future monitoring observations should these data show only a monotonic increase or decrease in the flux densities of the flat-spectrum components.     

Models of the Cosmic Horseshoe gravitational lens J1004+4112

We model the extremely massive and luminous lens galaxy in the Cosmic Horseshoe Einstein ring system J1004+4112, recently discovered in the Sloan Digital Sky Survey. We use the semilinear method of Warren & Dye, which pixelizes the source surface brightness distribution, to invert the Einstein ring for sets of parametrized lens models. Here, the method is refined by exploiting Bayesian inference to optimise adaptive pixelization of the source plane and to choose between three differently parametrized models: a singular isothermal ellipsoid, a power-law model and a Navarro, Frenk & White (NFW) profile. The most probable lens model is the power law with a volume mass density  ρ∝ r ^{−1.96±0.02}  and an axis ratio of ∼0.8. The mass within the Einstein ring (i.e. within a cylinder with projected distance of ∼30 kpc from the centre of the lens galaxy) is  (5.02 ± 0.09) × 10¹² M _⊙  , and the mass-to-light ratio is ∼30. Even though the lens lies in a group of galaxies, the preferred value of the external shear is almost zero. This makes the Cosmic Horseshoe unique amongst large separation lenses, as almost all the deflection comes from a single, very massive galaxy with little boost from the environment.     

The edge-on spiral gravitational lens B1600+434

We present new observations of the gravitational lens (GL) system B1600+434, strongly suggesting that the lens is an edge-on spiral galaxy. These observations are used to constrain the mass model of the system, in particular the oblateness and velocity dispersion of the dark matter halo around the lensing galaxy. From an analytical model we find a lower limit on the halo oblateness q _halo=( c/a )_ρ≳0.4; more detailed numerical models give a lower limit of q _halo≳0.5. We determine an average halo velocity dispersion of σ_halo=190±15 km s⁻¹ over all non-singular isothermal elliptical (NIE) halo models. Constraining the models to larger and more massive discs decreases this average by only 10 km s⁻¹. A lower limit of σ_halo≳150 km s⁻¹ is found, even for disc masses larger than the mass inside the Einstein radius. This lower limit indicates the need for a massive dark matter halo, contributing at least half of the mass inside the Einstein radius. Time-delay calculations give (54±3)/ h ₅₀ d for the NIE halo model and (70±4)/ h ₅₀ d for the modified Hubble profile (MHP) halo model. Although the time delay for both NIE and MHP halo models is well constrained on our parameter grid, it strongly depends on the halo surface density profile. We furthermore find that the presence of a flat luminous mass distribution can severely alter the statistical properties of the lens.     

The complex gravitational lens system B1933+503

We report the discovery of the most complex arcsec-scale radio gravitational lens system yet known. B1933+503 was found during the course of the CLASS survey and MERLIN and VLA radio maps reveal up to 10 components. Four of these are compact and have flat spectra; the rest are more extended and have steep spectra. The background lensed object appears to consist of a flat-spectrum core (quadruply imaged) and two compact 'lobes' symmetrically disposed relative to the core. One of the lobes is quadruply imaged while the other is doubly imaged. An HST observation of the system with the WFPC2 shows a galaxy with an axial ratio of 0.5, but none of the images of the background object is detected. A redshift of 0.755 has been measured for the lens galaxy.     

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.     

B0712+472: a new radio four-image gravitational lens

A new four-image gravitational lens system, B0712+472, has been discovered during the Cosmic Lens All-Sky Survey. This system consists of four flat-spectrum radio images that are also seen on a Hubble Space Telescope ( HST ) image, together with the lensing galaxy. We present MERLIN, VLA and VLBA maps and WHT spectra of the system as well as the HST images. The light distribution of the lensing galaxy is highly elongated and so too is the mass distribution deduced from modelling. We suggest a redshift of ∼1.33 for the lensed object; the lens redshift will require further investigation. The discovery of this new system further increases the ratio of four-image to two-image lens systems currently known, exacerbating problems of required ellipticity of matter distributions in lensing galaxies.     

A complete infrared Einstein ring in the gravitational lens system B1938 + 666

We report the discovery, using NICMOS on the Hubble Space Telescope , of an arcsec-diameter Einstein ring in the gravitational lens system B1938+666. The lensing galaxy is also detected, and is most likely an early-type galaxy. Modelling of the ring is presented and compared with the radio structure from MERLIN maps. We show that the Einstein ring is consistent with the gravitational lensing of an extended infrared component, centred between the two radio components.     

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.     

CLASS B2108+213: a new wide-separation gravitational lens system

We present observations of CLASS B2108+213, the widest separation gravitational lens system discovered by the Cosmic Lens All-Sky Survey. Radio imaging using the VLA at 8.46 GHz and MERLIN at 5 GHz shows two compact components separated by 4.56 arcsec with a faint third component in between which we believe is emission from a lensing galaxy. 5-GHz VLBA observations reveal milliarcsecond-scale structure in the two lensed images that is consistent with gravitational lensing. Optical emission from the two lensed images and two lensing galaxies within the Einstein radius is detected in Hubble Space Telescope imaging. Furthermore, an optical gravitational arc, associated with the strongest lensed component, has been detected. Surrounding the system is a number of faint galaxies which may help explain the wide image separation. A plausible mass distribution model for CLASS B2108+213 is also presented.     

IRAM observations of JVAS/CLASS gravitational lenses

We have searched for molecular absorption lines at millimetre wavelengths in 11 gravitational lens systems discovered in the JVAS/CLASS surveys of flat spectrum radio sources. Spectra of only one source 1030+074 were obtained in the 3-, 2- and 1.3-mm bands at the frequencies corresponding to common molecular transitions of CO and HCO⁺ as continuum emission was not found in any of the other sources. We calculated upper limits to the column density in molecular absorption for 1030+074, using an excitation temperature of 15 K, to be N _CO<6.3×10¹³ cm⁻² and N _HCO+<1.3×10¹¹ cm⁻² , equivalent to hydrogen column density of the order N _H<10¹⁸ cm⁻² , assuming standard molecular abundances. We also present the best upper limits of the continuum at the lower frequency for the other 10 gravitational lenses.     

Self-consistent gravitational lens reconstruction

We present a new method for directly determining accurate, self-consistent cluster lens mass and shear maps in the strong lensing regime from the magnification bias of background galaxies. The method relies upon pixellization of the surface mass density distribution which allows us to write down a simple, solvable set of equations. We also show how pixellization can be applied to methods of mass determination from measurements of shear and present a simplified method of application. The method is demonstrated with cluster models and applied to magnification data from the lensing cluster Abell 1689.     

NICMOS images of JVAS/CLASS gravitational lens systems

We present Hubble Space Telescope ( HST ) infrared images of four gravitational lens systems from the JVAS/CLASS gravitational lens survey and compare the new infrared HST pictures with previously published WFPC2 HST optical images and radio maps. Apart from the wealth of information that we get from the flux ratios and accurate positions and separations of the components of the lens systems, which we can use as inputs for better constraints on the lens models, we are able to discriminate between reddening and optical/radio microlensing as the possible cause of differences observed in the flux ratios of the components across the three wavelength bands. Substantial reddening has been known to be present in the lens system B1600+434 and has been further confirmed by the present infrared data. In the two systems B0712+472 and B1030+074 microlensing has been pinpointed as the main cause of the flux ratio discrepancy both in the optical/infrared and in the radio, the radio possibly caused by the substructure revealed in the lensing galaxies. In B0218+357, however, the results are still not conclusive. If we are actually seeing the two 'true' components of the lens system then the flux ratio differences are attributed to a combination of microlensing and reddening or are alternatively the result of some variability in at least one of the images. Otherwise the second 'true' component of B0218+357 may be completely absorbed by a molecular cloud and the anomalous flux density ratios and large difference in separation between the optical/infrared and radio that we see can be explained by emission either from a foreground object or from part of the lensing galaxy.