Galaxy density profiles and shapes – I. Simulation pipeline for lensing by realistic galaxy models

Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, which we aim to investigate in a consistent way. In this first paper in a series, we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii.     

Galaxy density profiles and shapes – II. Selection biases in strong lensing surveys

Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio and image separation distribution, with a flexible treatment of magnification bias to mimic different survey strategies. We present our results for two families of density profiles: cusped and deprojected Sérsic models. While we use fixed lens and source redshifts for most of the analysis, we show that the results are applicable to other redshift combinations, and we also explore the physics of how our results change for very different redshifts. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sérsic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population.     

Discovery of four gravitational lensing systems by clusters in the SDSS DR6

We report the discovery of 4 strong gravitational lensing systems by visual inspections of the Sloan Digital Sky Survey images of galaxy clusters in Data Release 6 (SDSS DR6). Two of the four systems show Einstein rings while the others show tangen-tial giant arcs. These arcs or rings have large angular separations ( 8") from the bright central galaxies and show bluer color compared with the red cluster galaxies. In addition,we found 5 probable and 4 possible lenses by galaxy clusters.     

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.     

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.     

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.     

Weak gravitational flexion

Flexion is the significant third-order weak gravitational lensing effect responsible for the weakly skewed and arc-like appearance of lensed galaxies. Here we demonstrate how flexion measurements can be used to measure galaxy halo density profiles and large-scale structure on non-linear scales, via galaxy–galaxy lensing, dark matter mapping and cosmic flexion correlation functions. We describe the origin of gravitational flexion, and discuss its four components, two of which are first described here. We also introduce an efficient complex formalism for all orders of lensing distortion. We proceed to examine the flexion predictions for galaxy–galaxy lensing, examining isothermal sphere and Navarro–Frenk–White (NFW) profiles and both circularly symmetric and elliptical cases. We show that in combination with shear we can precisely measure galaxy masses and NFW halo concentrations. We also show how flexion measurements can be used to reconstruct mass maps in two-dimensional projection on the sky, and in three dimensions in combination with redshift data. Finally, we examine the predictions for cosmic flexion, including convergence–flexion cross-correlations, and we find that the signal is an effective probe of structure on non-linear scales.     

Angular cross-correlation of galaxies: a probe of gravitational lensing by large-scale structure

The angular cross-correlation between two galaxy samples separated in redshift is shown to be a useful measure of weak lensing by large-scale structure. Angular correlations in faint galaxies arise as a result of spatial clustering of the galaxies as well as gravitational lensing by dark matter along the line of sight. The lensing contribution to the two-point autocorrelation function is typically small compared with the gravitational clustering. However, the cross-correlation between two galaxy samples is almost unaffected by gravitational clustering provided that their redshift distributions do not overlap. The cross-correlation is then induced by magnification bias resulting from lensing by large-scale structure. We compute the expected amplitude of the cross-correlation for popular theoretical models of structure formation. For two populations with mean redshifts of ≃0.3 and 1, we find a cross-correlation signal of ≃1 per cent on arcmin scales and ≃3 per cent on scales of a few arcsec. The dependence on the cosmological parameters Ω and Λ, the dark matter power spectrum and the bias factor of the foreground galaxy population is explored.     

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.     

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.     

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.     

An atlas of predicted exotic gravitational lenses

We present an investigation of the relationships between the radio properties of a giant radio galaxy MRC B0319−454 and the surrounding galaxy distribution with the aim of examining the influence of intergalactic gas and gravity associated with the large-scale structure on the evolution in the radio morphology. Our new radio continuum observations of the radio source, with high surface brightness sensitivity, images the asymmetries in the megaparsec-scale radio structure in total intensity and polarization. We compare these with the three-dimensional galaxy distribution derived from galaxy redshift surveys. Galaxy density gradients are observed along and perpendicular to the radio axis: the large-scale structure is consistent with a model wherein the galaxies trace the ambient intergalactic gas and the evolution of the radio structures are ram-pressure limited by this associated gas. Additionally, we have modelled the off-axis evolution of the south-west radio lobe as deflection of a buoyant jet backflow by a transverse gravitational field: the model is plausible if entrainment is small. The case study presented here is a demonstration that giant radio galaxies may be useful probes of the warm-hot intergalactic medium believed to be associated with moderately over dense galaxy distributions.     

QSO lensing magnification associated with galaxy groups

We simulated both the matter and light (galaxy) distributions in a wedge of the Universe and calculated the gravitational lensing magnification caused by the mass along the line-of-sight of galaxies and galaxy groups identified in sky surveys. A large volume redshift cone containing cold dark matter particles mimics the expected cosmological matter distribution in a flat universe with low matter density and a cosmological constant. We generate a mock galaxy catalogue from the matter distribution and identify thousands of galaxy groups in the luminous sky projection. We calculate the expected magnification around galaxies and galaxy groups and then the induced quasi-stellar object (QSO)–lens angular correlation due to magnification bias. This correlation is observable and can be used both to estimate the average mass of the lens population and to make cosmological inferences. We also use analytical calculations and various analyses to compare the observational results with theoretical expectations for the cross-correlation between faint QSOs from the 2dF Survey and nearby galaxies and groups from the Automated Plate Measurement and Sloan Digital Sky Survey Early Data Release. The observed QSO–lens anticorrelations are stronger than the predictions for the cosmological model used. This suggests that there could be unknown systematic errors in the observations and data reduction, or that the model used is not adequate. If the observed signal is assumed to be solely due to gravitational lensing, then the lensing is stronger than expected, due to more massive galactic structures or more efficient lensing than simulated.     

Weighing a galaxy bar in the lens Q2237 +=0305

In the gravitational lens system Q2237+0305 the cruciform quasar image geometry is twisted by 10c by the lens effect of a bar in the lensing galaxy. This effect can be used to measure the mass of the bar. We construct a new lensing model for this system with a power-law elliptical bulge and a Ferrers bar. The observed ellipticity of the optical isophotes of the galaxy leads to a nearly isothermal elliptical profile for the bulge, with a total quasar magnification of 16⁺⁵₋₄. We measure a bar mass of (7.5 ∼ 1.5) −10⁸ h ⁻¹₇₅ M⊙ in the region inside the quasar images.     

Testing Bekenstein's relativistic Modified Newtonian Dynamics with lensing data

We propose to use multiple-imaged gravitational lenses to set limits on gravity theories without dark matter, specifically tensor–vector–scalar (TeVeS) theory, a theory which is consistent with fundamental relativistic principles and the phenomenology of Modified Newtonian Dynamics (MOND) theory. After setting the framework for lensing and cosmology, we analytically derive the deflection angle for the point lens and the Hernquist galaxy profile, and study their patterns in convergence, shear and amplification. Applying our analytical lensing models, we fit galaxy-quasar lenses in the CfA-Arizona Space Telescope Lens Survey (CASTLES) sample. We do this with three methods, fitting the observed Einstein ring sizes, the image positions, or the flux ratios. In all the cases, we consistently find that stars in galaxies in MOND/TeVeS provide adequate lensing. Bekenstein's toy μ function provides more efficient lensing than the standard MOND μ function. But for a handful of lenses, a good fit would require a lens mass orders of magnitude larger/smaller than the stellar mass derived from luminosity unless the modification function μ and modification scale a ₀ for the universal gravity were allowed to be very different from what spiral galaxy rotation curves normally imply. We discuss the limitation of present data and summarize constraints on the MOND μ function. We also show that the simplest TeVeS 'minimal-matter' cosmology, a baryonic universe with a cosmological constant, can fit the distance–redshift relation from the supernova data, but underpredicts the sound horizon size at the last scattering. We conclude that lensing is a promising approach to differentiate laws of gravity.     

Extragalactic integral field spectroscopy on the Gemini telescopes

We have been undertaking a programme on the Gemini 8‐m telescopes to demonstrate the power of integral field spectroscopy, using the optical GMOS spectrograph, and the new CIRPASS instrument in the near‐infrared. Here we present some preliminary results from 3D spectroscopy of extra‐galactic objects, mapping the emission lines in a 3CR radio galaxy and in a gravitationally lensed arc, exploring dark matter sub‐structure through observations of an Einstein Cross gravitational lens, and the star formation time‐scales of young massive clusters in the starburst galaxy NGC 1140. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Depletion of background galaxies owing to the cluster lens CL0024+1654: U- and R-band observations

We have obtained U - and R -band observations of the depletion of background galaxies resulting from the gravitational lensing of the galaxy cluster CL0024+1654 ( z =0.39). The radial depletion curves show a significant depletion in both bands within a radius of 40–70 arcsec from the cluster centre. This is the first time that depletion is detected in the U band. This gives independent evidence for a break in the slope of the U -band luminosity function at faint magnitudes. The radially averaged R -band depletion curve is broader and deeper than in the U band. The differences can be attributed to the wavelength dependence of the slope of the luminosity function and to the different redshift distribution of the objects probed in the two bands. We estimate the Einstein radius, r _E, of a singular isothermal sphere lens model using maximum-likelihood analysis. Adopting a slope of the number counts of α =0.2 and using the background density found beyond r =150 arcsec, we find r _E=17±3 and 25±3 arcsec in the U and R bands, respectively. When combined with the redshift of the single background galaxy at z =1.675 seen as four giant arcs around 30 arcsec from the cluster centre, these values indicate a median redshift in the range 〈 z _S〉≈0.7 to 1.1 for the U _AB≥24 mag and R _AB≥24 mag populations.     

An analytic model for non-spherical lenses in covariant MOdified Newtonian Dynamics

Strong gravitational lensing by galaxies in MOdified Newtonian Dynamics (MOND) has until now been restricted to spherically symmetric models. These models were able to account for the size of the Einstein ring of observed lenses, but were unable to account for double-imaged systems with collinear images, as well as four-image lenses. Non-spherical models are generally cumbersome to compute numerically in MOND, but we present here a class of analytic non-spherical models that can be applied to fit double-imaged and quadruple-imaged systems. We use them to obtain a reasonable MOND fit to 10 double-imaged systems, as well as to the quadruple-imaged system Q2237+030 which is an isolated bulge-disc lens producing an Einstein cross. However, we also find five double-imaged systems and three quadruple-imaged systems for which no reasonable MOND fit can be obtained with our models. We argue that this is mostly due to the intrinsic limitation of the analytic models, even though the presence of small amounts of additional dark mass on galaxy scales in MOND is also plausible.     

SDSS J120923.7+264047: a new massive galaxy cluster with a bright giant arc

Highly magnified lensed galaxies allow us to probe the morphological and spectroscopic properties of high-redshift stellar systems in great detail. However, such objects are rare, and there are only a handful of lensed galaxies that are bright enough for a high-resolution spectroscopic study with current instrumentation. We report the discovery of a new massive lensing cluster, SDSS J120923.7+264047, at z = 0.558. Present around the cluster core, at angular distances of up to ∼40 arcsec, are many arcs and arc candidates, presumably due to lensing of background galaxies by the cluster gravitational potential. One of the arcs, 21 arcsec long, has an r -band magnitude of 20, making it one of the brightest known lensed galaxies. We obtained a low-resolution spectrum of this galaxy, using the Keck-I telescope, and found it is at redshift of z = 1.018.