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.     

The anatomy of a quadruply imaged gravitational lens system

The key to using a strong gravitational lens system to measure the Hubble constant is to obtain an accurate model of the lens potential. In this paper, we investigate the properties of gravitational lens B1608+656, a quadruply imaged lens system with an extended source intensity distribution. Our analysis is valid for generic quadruply lensed systems. Limit curves and isophotal separatrices are defined for such systems, and we show that the isophotal separatrices must intersect at the critical curves and the satellite isophotes must be tangent to the limit curves. The most recent model of B1608+656 by Koopmans et al. satisfies these criteria for some, but not all, of the isophotal separatrices within the observational uncertainty. We study a non-parametric method of potential reconstruction proposed by Blandford, Surpi & Kundic and demonstrate that although the method works in principle and elucidates image formation, the initial potential only converges to the true model when it is within ∼1 per cent of the true model.     

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.     

The lens and source of the optical Einstein ring gravitational lens ER 0047–2808

We perform a detailed analysis of the optical gravitational lens ER 0047–2808 imaged with the Wide Field Planetary Camera 2 on the Hubble Space Telescope . Using software specifically designed for the analysis of resolved gravitational lens systems, we focus on how the image alone can constrain the mass distribution in the lens galaxy. We find that the data are of sufficient quality to strongly constrain the lens model with no a priori assumptions about the source. Using a variety of mass models, we find statistically acceptable results for elliptical isothermal-like models with an Einstein radius of 1.17 arcsec. An elliptical power-law model  (Σ∝ R ^−β)  for the surface mass density favours a slope slightly steeper than isothermal with  β= 1.08 ± 0.03  . Other models including a constant mass-to-light ratio (M/L), pure Navarro, Frenk & White halo and (surprisingly) an isothermal sphere with external shear are ruled out by the data. We find the galaxy light profile can only be fit with a Sérsic plus point-source model. The resulting total  M/L _B   contained within the images is  4.7  h ₆₅± 0.3  . In addition, we find the luminous matter is aligned with the total mass distribution within a few degrees. This is the first time a resolved optical gravitational lens image has been quantitatively reproduced using a non-parametric source.
The source, reconstructed by the software, is revealed to have two bright regions, with an unresolved component inside the caustic and a resolved component straddling a fold caustic. The angular size of the entire source is ∼0.1 arcsec and its (unlensed) Lyα flux is  3 × 10⁻¹⁷ erg s⁻¹ cm⁻²  .     

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.     

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.     

MERLIN/VLA imaging of the gravitational lens system B0218+357

Gravitational lenses offer the possibility of accurately determining the Hubble parameter ( H ₀) over cosmological distances, and B0218+357 is one of the most promising systems for an application of this technique. In particular this system has an accurately measured time delay (10.5±0.4 d) and preliminary mass modelling has given a value for H ₀ of 69₋₁₉⁺¹³ km s⁻¹ Mpc⁻¹. The error on this estimate is now dominated by the uncertainty in the mass modelling. As this system contains an Einstein ring it should be possible to constrain the model better by imaging the ring at high resolution. To achieve this we have combined data from MERLIN and the Very Large Array (VLA) at a frequency of 5 GHz. In particular MERLIN has been used in multifrequency mode in order to improve substantially the aperture coverage of the combined data set. The resulting map is the best that has been made of the ring and contains many new and interesting features. Efforts are currently underway to exploit the new data for lensing constraints using the LensClean algorithm.     

Modelling the 10-image lensed system B1933+503

A gravitational lens model is presented for the newly discovered 10-image system B1933+503. The underlying object, revealed by modelling, is a triple radio source on the scale of a couple of hundred mas that is well-aligned along the line of sight with a foreground and somewhat flattened lensing galaxy, the orientation and location of which match those of an observed galaxy, known to be at a redshift of 0.755. Uncertainties in the modelling are obtained by a Monte Carlo exercise. Observational tests of the lens model are proposed, and the time delays between various pairs of images are determined, as the core of the source is known to be significantly variable. Future observations of the lens hold the key to using B1933+503 to constrain Hubble's constant. Despite the absence of a source redshift, the utility of the system as a probe of the structure of the lens galaxy is unparalleled as it provides a surfeit of easily identifiable constraints for modelling the system.     

The new gravitational lens system B1030 + 074

We report the discovery of a new double-image gravitational lens system, B1030 + 074, which was found during the Jodrell Bank–VLA Astrometric Survey (JVAS). We have collected extensive radio data on the system using the VLA, MERLIN, the EVN and the VLBA, and optical observations using WFPC2 on the HST . The lensed images are separated by 1.56 arcsec and their flux density ratio at centimetric wavelengths is approximately 14:1, although the ratio is slightly frequency-dependent and the images appear to be time-variable. The HST pictures show both the lensed images and the lensing galaxy close to the weaker image. The lensing galaxy has substructure which could be a spiral arm or an interacting galaxy. We have modelled B1030 + 074 using a singular isothermal ellipsoid which yields a time delay of 156/ h ₅₀ d. This lens is likely to be suitable for the measurement of the Hubble constant.     

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 new gravitational lens from the MUSCLES survey: ULAS J082016.1+081216

We present observations of a new double-image gravitational lens system, ULAS J082016.1+081216, of image separation 2.3 arcsec and high (∼6) flux ratio. The system is selected from the Sloan Digital Sky Survey (SDSS) spectroscopic quasar list using new high-quality images from the UKIRT (United Kingdom Infrared Telescope) Deep Sky Survey (UKIDSS). The lensed quasar has a source redshift of 2.024, and we identify the lens galaxy as a faint red object of redshift  0.803 ± 0.001  . Three other objects from the UKIDSS survey, selected in the same way, were found not to be lens systems. Together with the earlier lens found using this method, the SDSS–UKIDSS lenses have the potential to significantly increase the number of quasar lenses found in SDSS, to extend the survey to higher flux ratios and lower separations, and to give greater completeness which is important for statistical purposes.     

High-resolution observations and mass modelling of the CLASS gravitational lens B1152+199

We present a series of high-resolution radio and optical observations of the CLASS gravitational lens system B1152+199 obtained with the Multi-Element Radio-Linked Interferometer Network, Very Long Baseline Array and Hubble Space Telescope . Based on the milliarcsecond-scale substructure of the lensed radio components and precise optical astrometry for the lensing galaxy, we construct models for the system and place constraints on the galaxy mass profile. For a single galaxy model with surface mass density  Σ(r)∝r^{− β}   , we find that  0.95 β 1.21  at 2 σ confidence. Including a second deflector to represent a possible satellite galaxy of the primary lens leads to slightly steeper mass profiles.     

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.