Substructure and the cusp and fold relations

Gravitational lensing of a background source by a foreground galaxy lens occasionally produces four images of the source. The cusp and the fold relations impose conditions on the ratios of magnifications of these four-image lenses. In this theoretical investigation, we explore the sensitivity of these relations to the presence of substructure in the lens. Starting with a smooth lens potential, we add varying amounts of substructure, while keeping the source position fixed, and find that the fold relation is a more robust indicator of substructure than the cusp relation for the images. This robustness is independent of the detailed spatial distribution of the substructure, as well as of the ellipticity of the lensing potential and the presence of external shear.     

Perturbation theory of N point-mass gravitational lens systems without symmetry: small mass-ratio approximation

This paper makes the first systematic attempt to determine using perturbation theory the positions of images by gravitational lensing due to arbitrary number of coplanar masses without any symmetry on a plane, as a function of lens and source parameters. We present a method of Taylor-series expansion to solve the lens equation under a small mass-ratio approximation. First, we investigate perturbative structures of a single-complex-variable polynomial, which has been commonly used. Perturbative roots are found. Some roots represent positions of lensed images, while the others are unphysical because they do not satisfy the lens equation. This is consistent with a fact that the degree of the polynomial, namely the number of zeros, exceeds the maximum number of lensed images if   N = 3  (or more). The theorem never tells which roots are physical (or unphysical). In this paper, unphysical ones are identified. Secondly, to avoid unphysical roots, we re-examine the lens equation. The advantage of our method is that it allows a systematic iterative analysis. We determine image positions for binary lens systems up to the third order in mass ratios and for arbitrary N point masses up to the second order. This clarifies the dependence on parameters. Thirdly, the number of the images that admit a small mass-ratio limit is less than the maximum number. It is suggested that positions of extra images could not be expressed as Maclaurin series in mass ratios. Magnifications are finally discussed.     

The Milky Way Galaxy as a strong gravitational lens

We study the gravitational lensing effects of spiral galaxies by taking a model of the Milky Way and computing its lensing properties. The model is composed of a spherical Hernquist bulge, a Miyamoto–Nagai disc and an isothermal halo. As a strong lens, a spiral galaxy like the Milky Way can give rise to four different imaging geometries. They are (i) three images on one side of the galaxy centre ('disc triplets'), (ii) three images with one close to the centre ('core triplets'), (iii) five images and (iv) seven images. Neglecting magnification bias, we show that the core triplets, disc triplets and fivefold imaging are roughly equally likely. Even though our models contain edge-on discs, their image multiplicities are not dominated by disc triplets. The halo is included for completeness, but it has a small effect on the caustic structure, the time delays and brightnesses of the images.
The Milky Way model has a maximum disc (i.e. the halo is not dynamically important in the inner parts). Strong lensing by nearly edge-on disc galaxies breaks the degeneracy between the relative contributions of the disc and halo to the overall rotation curve. If a spiral galaxy has a submaximum disc, then the astroid caustic shrinks dramatically in size, whilst the radial caustic shrinks more modestly. This causes changes in the relative likelihood of the image geometries, specifically (i) core triplets are now ∼9/2 times more likely than disc triplets, (ii) the cross-section for threefold imaging is reduced by a factor of ∼2/3, whilst (iii) the cross-section for fivefold imaging is reduced by ∼1/2. Although multiple imaging is less likely (the cross-sections are smaller), the average total magnification is greater. The time delays are smaller, as the total projected lensing mass is reduced.     

Systematic uncertainties in gravitational lensing models: a semi-analytical study of PG1115+080

While the Hubble constant can be derived from observable time delays between images of lensed quasars, the result is often highly sensitive to assumptions and systematic uncertainties in the lensing model. Unlike most previous authors, we put minimal restrictions on the radial profile of the lens and allow for non-elliptical lens potentials. We explore these effects using a broad class of models with a lens potential     which has an unrestricted radial profile but self-similar iso-potential contours defined by     For these potentials, the lens equations can be solved semi-analytically. The axis ratio and position angle of the lens can be determined from the image positions of quadruple gravitational lensed systems directly, independent of the radial profile. We give simple equations for estimating the power-law slope of the lens density directly from the image positions and for estimating the time delay ratios. Our method greatly simplifies the numerics for fitting observations and is fast in exploring the model parameter space. As an illustration, we apply the model to PG1115+080. An entire one-parameter sequence of models fits the observations exactly. We show that the measured image positions and time delays do not uniquely determine the Hubble constant.     

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.     

Non-parametric strong lens inversion of Cl 0024+1654: illustrating the monopole degeneracy

The cluster lens Cl 0024+1654 is undoubtedly one of the most beautiful examples of strong gravitational lensing, providing five large images of a single source with well-resolved substructure. Using the information contained in the positions and the shapes of the images, combined with the null space information, a non-parametric technique is used to infer the strong lensing mass map of the central region of this cluster. This yields a strong lensing mass of  1.60 × 10¹⁴ M_⊙  within a 0.5  arcmin radius around the cluster centre. This mass distribution is then used as a case study of the monopole degeneracy, which may be one of the most important degeneracies in gravitational lensing studies and which is extremely hard to break. We illustrate the monopole degeneracy by adding circularly symmetric density distributions with zero total mass to the original mass map of Cl 0024+1654. These redistribute mass in certain areas of the mass map without affecting the observed images in any way. We show that the monopole degeneracy and the mass-sheet degeneracy together lie at the heart of the discrepancies between different gravitational lens reconstructions that can be found in the literature for a given object, and that many images/sources, with an overall high image density in the lens plane, are required to construct an accurate, high-resolution mass map based on strong lensing data.     

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.     

Asymptotic formulas for the magnification of a gravitational lens system near a fold caustic

An approximate formula for the magnification of a point source near a fold caustic obtained in the first linear caustic approximation is widely used in the theory of gravitational lens systems. Here, this formula is refined to include the post-linear terms that have been found both for a point source and for an extended Gaussian source in the absence of continuous matter on the line of sight. The formulas are reduced to a form containing three additional parameters; the derivation of nontrivial corrections requires including the expansion terms in the lens equation up to the fourth order. The modified formula for an extended source is used to analyze strong microlensing events in the gravitational lens system Q2237+0305 (the Einstein Cross). For such an event on the light curve of image C (1999, OGLE data), the corrections found are statistically significant.     

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

From planes to spheres: about gravitational lens magnifications

We discuss the classic theorem according to which a gravitational lens always produces at least one image with a magnification greater than unity. This theorem seems to contradict the conservation of total flux from a lensed source. The standard solution to this paradox is based on the exact definition of the reference 'unlensed' situation, in which the lens mass can be either removed or smoothly redistributed.
We calculate magnifications and amplifications (in photon number and energy flux density) for general lensing scenarios not limited to regions close to the optical axis. In this way the formalism is naturally extended from tangential planes for the source and lensed images to complete spheres. We derive the lensing potential theory on the sphere and find that the Poisson equation is modified by an additional source term that is related to the mean density and to the Newtonian potential at the positions of observer and source. This new term generally reduces the magnification, to below unity far from the optical axis, and ensures conservation of the total photon number received on a sphere around the source.
This discussion does not affect the validity of the focusing theorem , in which the unlensed situation is defined to have an unchanged affine distance between source and observer. The focusing theorem does not contradict flux conservation, because the mean total magnification (or amplification) directly corresponds to different areas of the source (or observer) sphere in the lensed and unlensed situation. We argue that a constant affine distance does not define an astronomically meaningful reference.
By exchanging source and observer, we confirm that magnification and amplification differ according to Etherington's reciprocity law, so that surface brightness is no longer strictly conserved. At this level we also have to distinguish between different surface brightness definitions that are based on photon number, photon flux and energy flux.     

Lens models with density cusps

Lens models appropriate for representing cusped galaxies and clusters are developed. The analogue of the odd-number theorem for cusped density distributions is given. Density cusps are classified into strong, isothermal or weak, according to their lensing properties. Strong cusps cause multiple imaging for any source position, whereas isothermal and weak cusps give rise to only one image for distant sources. Isothermal cusps always possess a pseudo-caustic. When the source crosses the pseudo-caustic, the number of images changes by unity.   Two families of cusped galaxy and cluster models are examined in detail. The double power-law family has an inner cusp, followed by a transition region and an outer envelope. One member of this family — the isothermal double power-law model — possesses an exceedingly scarce property, namely the lens equation is exactly solvable for any source position. This means that the magnifications, the time delay and the lensing cross-sections are all readily available. The model has a three-dimensional density that is cusped like r ⁻² at small radii and falls off like r ⁻⁴ asymptotically. Thus, it provides a very useful representation of the lensing properties of a galaxy or cluster of finite total mass with a flat rotation curve. The second set of models studied is the single power-law family. These are single density cusps of infinite extent. The properties of the critical curves and caustics and the behaviour of the lenses in the presence of external shear are all discussed in some detail.     

The influence of central black holes on gravitational lenses

Recent observations indicate that many if not all galaxies host massive central black holes. In this paper we explore the influence of black holes on the lensing properties. We model the lens as an isothermal ellipsoid with a finite core radius plus a central black hole. We show that the presence of the black hole substantially changes the critical curves and caustics. If the black hole mass is above a critical value, then it will completely suppress the central images for all source positions. Realistic central black holes are likely to have masses below this critical value. Even in such subcritical cases, the black hole can suppress the central image when the source is inside a zone of influence, which depends on the core radius and black hole mass. In the subcritical cases, an additional image may be created by the black hole in some regions, which for some radio lenses may be detectable with high-resolution and large dynamic range VLBI maps. The presence of central black holes should also be taken into account when one constrains the core radius from the lack of central images in gravitational lenses.     

The Hubble constant from galaxy lenses: impacts of triaxiality and model degeneracies

The Hubble constant can be constrained using the time delays between multiple images of gravitationally lensed sources. In some notable cases, typical lensing analyses assuming isothermal galaxy density profiles produce low values for the Hubble constant, inconsistent with the result of the HST Key Project  (72 ± 8 km s⁻¹ Mpc⁻¹)  . Possible systematics in the values of the Hubble constant derived from galaxy lensing systems can result from a number of factors, for example, neglect of environmental effects, assumption of isothermality, or contamination by line-of-sight structures. One additional potentially important factor is the triaxial structure of the lensing galaxy halo; most lens models account for halo shape simply by perturbing the projected spherical lensing potential, an approximation that is often necessary but that is inadequate at the levels of triaxiality predicted in the cold dark matter paradigm. To quantify the potential error introduced by this assumption in estimates of the Hubble parameter, we strongly lens a distant galaxy through a sample of triaxial softened isothermal haloes and use an Markov Chain Monte Carlo method to constrain the lensing halo profile and the Hubble parameter from the resulting multiple image systems. We explore the major degeneracies between the Hubble parameter and several parameters of the lensing model, finding that without a way to accurately break these degeneracies accurate estimates of the Hubble parameter are not possible. Crucially, we find that triaxiality does not significantly bias estimates of the Hubble constant, and offer an analytic explanation for this behaviour in the case of isothermal profiles. Neglected triaxial halo shape cannot contribute to the low Hubble constant values derived in a number of galaxy lens systems.     

Substructures in lens galaxies: PG1115+080 and B1555+375, two fold configurations

We study the anomalous flux ratio which is observed in some four-image lens systems, where the source lies close to a fold caustic. In this case two of the images are close to the critical curve and their flux ratio should be equal to unity, instead in several cases the observed value differs significantly. The most plausible solution is to invoke the presence of substructures, as for instance predicted by the Cold Dark Matter scenario, located near the two images. In particular, we analyze the two fold lens systems PG1115+080 and B1555+375, for which there are not yet satisfactory models which explain the observed anomalous flux ratios. We add to a smooth lens model, which reproduces well the positions of the images but not the anomalous fluxes, one or two substructures described as singular isothermal spheres. For PG1115+080 we consider a smooth model with the influence of the group of galaxies described by a SIS and a substructure with mass ∼10⁵ M _⊙ as well as a smooth model with an external shear and one substructure with mass ∼10⁸ M _⊙. For B1555+375 either a strong external shear or two substructures with mass ∼10⁷ M _⊙ reproduce the data quite well.     

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.     

Rotation in gravitational lenses

Gravitational lensing deflects light. A single lens deflector can only shear images, but cannot induce rotations. Multiple lens planes can induce rotations. Such rotations can be observed in quadruply imaged sources, and can be used to distinguish between two proposed solutions of the flux anomaly problem: substructures in lensing galaxies versus large-scale structure. We predict the expected amount of rotation due to large-scale structure in strong lensing systems, and show how this effect can be measured using ∼mas very long baseline interferometry astrometry of quadruple lenses with extended source structures. The magnitude of rotation is around 1°. The biggest theoretical uncertainty is the power spectrum of dark matter on very small scales. This procedure can potentially be turned around to measure the dark matter power spectrum on very small scales. We list the predicted rms rotation angles for several quadruple lenses with known lens and source redshifts.     

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.     

Submillimetre observations of luminous z>4 radio-quiet quasars and the contribution of AGN to the submm source population

We present polarization observations of the gravitational lens system B1422+231 made at 8.4 GHz using the VLBA and the 100-m telescope at Effelsberg. All four images of the quasar show structure on the milliarcsec scale. The three bright images show tangential stretching as expected from lens models. Some basic properties of gravitational lensing are exhibited by this system. The surface brightness of images A and B are the same and the parity reversal expected in image B is revealed, for the first time, by polarization observations. There is a large differential Faraday rotation between images A and B.     

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.     

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.