Analytic relations for magnifications and time delays in gravitational lenses with fold and cusp configurations

Gravitational lensing provides a unique and powerful probe of the mass distributions of distant galaxies. Four-image lens systems with fold and cusp configurations have two or three bright images near a critical point. Within the framework of singularity theory, we derive analytic relations that are satisfied for a light source that lies a small but finite distance from the astroid caustic of a four-image lens. Using a perturbative expansion of the image positions, we show that the time delay between the close pair of images in a fold lens scales with the cube of the image separation, with a constant of proportionality that depends on a particular third derivative of the lens potential. We also apply our formalism to cusp lenses, where we develop perturbative expressions for the image positions, magnifications and time delays of the images in a cusp triplet. Some of these results were derived previously for a source asymptotically close to a cusp point, but using a simplified form of the lens equation whose validity may be in doubt for sources that lie at astrophysically relevant distances from the caustic. Along with the work of Keeton, Gaudi & Petters, this paper demonstrates that perturbation theory plays an important role in theoretical lensing studies.     

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.     

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.     

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.     

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.     

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.     

Modelling the first probable two-plane lens system B2114+022: reproducing two compact radio cores A and D

We test possible lensing scenarios of the JVAS system B2114+022, in which two galaxies at different redshifts ('G1' at z ₁=0.3157 and 'G2' at z ₂=0.5883) are found within 2 arcsec of quadruple radio sources. For our investigation, we use possible lensing constraints derived from a wealth of data on the radio sources obtained with VLA, MERLIN, VLBA and EVN as well as HST imaging data on the two galaxies, which were presented recently in Augusto et al. In the present study, we focus on reproducing the widest separated, observationally similar radio components A and D as lensed images. We first treat G2 (which is the more distant one from the geometric centre) as a shear term, and then consider two-plane lensing explicitly including G2's potential at the z ₂ plane as the first case of two-plane lens modelling. Our modelling results not only support the hypothesis that the system includes gravitationally lensed images of a higher-redshift extragalactic object, but they also show that the explicit inclusion of G2's potential at the second lens plane is necessary in order to fit the data with astrophysically plausible galaxy parameters. Finally, we illustrate a natural consequence of a two-plane lens system, namely the prediction of distortion as well as shift and stretching of G2's isophotes by G1's potential, which can in principle be measured by subtracting out G1's light distribution in an image of high signal-to-noise ratio and good angular resolution, especially a multicolour one.     

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.     

Gravitational lensing of anisotropic sources

In strong gravitational lensing, the multiple images we see correspond to light rays that leave the source in slightly different directions. If the source emission is anisotropic, the images may differ from conventional lensing predictions (which assume isotropy). To identify scales on which source anisotropy may be important, we study the angle δ between the light rays emerging from the source, for different lensing configurations. If the lens has a power-law profile   M ∝ R ^γ  , the angle δ initially increases with lens redshift and then either diverges (for a steep profile  γ < 1  ), remains constant (for an isothermal profile  γ= 1  ), or vanishes (for a shallow profile  γ > 1  ) as   z _l→ z _s  . The scaling with lens mass is roughly  δ∝ M ^1/(2−γ)  . The results for an Navarro–Frenk–White (NFW) profile are qualitatively similar to those for a shallow power law, with δ peaking at about half the redshift of the source (not half the distance). In practice, beaming could modify the statistics of beamed sources lensed by massive clusters: for an opening angle  θ_jet  , there is a probability as high as   P ∼ 0.02–0.07(θ_jet/0.5°)⁻¹  that one of the lensed images may be missed (for  2 ≲ z _s≲ 6  ). Differential absorption within active galactic nuclei (AGNs) could modify the flux ratios of AGNs lensed by clusters; a sample of AGNs lensed by clusters could provide further constraints on the sizes of absorbing regions. Source anisotropy is not likely to be a significant effect in galaxy-scale strong lensing.     

Analytic relations between the observed gravitational microlensing parameters with and without the effect of blending

When a microlensing light curve is contaminated by blended light from unresolved stars near the line of sight to the lensed star, the light curve shape and corresponding parametrization for the event will differ from the values expected when the event is not affected by blending. As a result, blending makes it difficult to identify the major lens population and to estimate the amount of lensing matter. In order to estimate the effect of blending on the result of lensing experiments, it is, therefore, essential to know how the observed lensing parameters change depending on the fraction of blended light. Previously, the changed lensing parameters were obtained with a statistical method that not only required a large amount of computation time but also was prone to uncertainty. In this paper, we derive analytic relations between the lensing parameters with and without the effect of blending. By using these relations, we investigate the dependence of the observed lensing parameters on the amount of blended light, the impact parameter and the threshold amplification for event detection.     

Discs in early-type lensing galaxies: effects on magnification ratios and measurements of H0

Observations of early-type galaxies, both in the local Universe and in clusters at medium redshifts, suggest that these galaxies often contain discs or disc-like structures. Using the results of Kelson et al. for the incidence of disc-components among the galaxies in the redshift   z = 0.33  cluster CL 1358+62, we investigate the effect of disc structures on the lensing properties of early-type galaxies. Statistical properties, like magnification cross-sections and the expected number of quad (four-image) lens systems, are not affected greatly by the inclusion of discs that contain less than ∼10 per cent of the total stellar mass. However, the properties of individual lens systems are affected. We estimate that 10–30 per cent of all quad-lens systems, with early-type deflector galaxies, would be affected measurably by the presence of disc components. Intriguingly, the image magnification ratios are altered significantly. The amplitude of the predicted change is sufficient to explain the observed magnification ratios in systems like B1422+231 without requiring compact substructure. Furthermore, time-delays between images also change; fitting a bulge-only model to early-type lenses that in fact contain a disc would yield a value of the Hubble constant H ₀ that is systematically too low by ∼25 per cent.     

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.     

Gravitational lensing in quasar samples

The first cosmic mirage was discovered approximately 20 years ago as the double optical counterpart of a radio source. This phenomenon had been predicted some 70 years earlier as a consequence of General Relativity. We present here a summary of what we have learnt since. The applications are so numerous that we had to concentrate on a few selected aspects of this new field of research. This review is focused on strong gravitational lensing, i.e. the formation of multiple images, in QSO samples. It is intended to give the reader an up-to-date status of the observations and to present an overview of its most interesting potential applications in cosmology and astrophysics, as well as numerous important results achieved so far. The first section follows an intuitive approach to the basics of gravitational lensing and is developed in view of our interest in multiply imaged quasars. The astrophysical and cosmological applications of gravitational lensing are outlined in Sect. 2 and the most important results are presented in Sect. 5. Sections 3 and 4 are devoted to the observations. Finally, conclusions are summarized in the last section. We have tried to avoid duplication with existing (and excellent) introductions to the field of gravitational lensing. For this reason, we did not concentrate on the individual properties of specific lens models, as these are already well presented in Narayan and Bartelmann (1996) and on a more intuitive ground in Refsdal and Surdej (1994). Wambsganss (1998) proposes a broad view on gravitational lensing in astronomy; the reviews by Fort and Mellier (1994) and Hattori et al. (1999) deal with lensing by galaxy clusters; microlensing in the Galaxy and the local group is reviewed by Paczyński (1996) and a general panorama on weak lensing is given by Bartelmann and Schneider (1999) and Mellier (1999). The monograph on the theory of gravitational lensing by Schneider, Ehlers and Falco (1992) also remains a reference in the field. Received 4 April 2000 / Published online 9 August 2000     

Strong gravitational lensing with SKA

The advent of new observational facilities in the last two decades has allowed the rapid discovery and high-resolution optical imaging of many strong lens systems from galaxy to cluster scales, as well as their spectroscopic follow-up. Radio telescopes have played the dominant role in the systematic detection of dozens of new arcsec-scale lens systems. For the future, we expect nothing less! The next major ground- and space-based facilities, especially the Square Kilometer Array can discover tens of thousands of new lens systems in large sky surveys. For optical imaging and spectroscopic follow-up a strong synergy with planned optical facilities is needed. Here, we discuss the field where strong gravitational lensing is expected to play the dominant role and where SKA can have a major impact: The study of the internal mass structure and evolution of galaxies and clusters to z 1. In addition, studies of more exotic phenomena are contemplated. For example, milli- and micro-lensing can provide a way to measure the mass-functions of stars and CDM substructure at cosmological distances. All-sky radio monitoring will also rapidly develop the field of time-domain lensing.     

Strong lensing optical depths in a ΛCDM universe

We investigate strong gravitational lensing in the concordance ΛCDM cosmology by carrying out ray tracing along past light cones through the Millennium Simulation, the largest simulation of cosmic structure formation ever carried out. We extend previous ray-tracing methods in order to take full advantage of the large volume and the excellent spatial and mass resolution of the simulation. As a function of source redshift we evaluate the probability that an image will be highly magnified, will be highly elongated or will be one of a set of multiple images. We show that such strong lensing events can almost always be traced to a single dominant lensing object and we study the mass and redshift distribution of these primary lenses. We fit analytic models to the simulated dark haloes in order to study how our optical depth measurements are affected by the limited resolution of the simulation and of the lensing planes that we construct from it. We conclude that such effects lead us to underestimate total strong lensing cross-sections by about 15 per cent. This is smaller than the effects expected from our neglect of the baryonic components of galaxies. Finally we investigate whether strong lensing is enhanced by material in front of or behind the primary lens. Although strong lensing lines of sight are indeed biased towards higher than average mean densities, this additional matter typically contributes only a few per cent of the total surface density.     

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.     

基于Shapelets基函数的引力透镜质量重构方法

引力透镜效应是探测星系团物质分布的有效方法之一.目前,利用引力透镜数据重构星系团质量分布的主流方法可以分为两大类,即参数法和非参数法.在实际研究工作中,受限于质量模型假设和计算分辨率等方面的影响,现有的重构算法仍有诸多亟需解决的问题.基于Shapelets基函数的引力透镜质量重构方法通过基函数来实现引力透镜质量重构,使用Shapelets基函数分解引力透镜势,以引力透镜中多重像的位置和背景星系椭率畸变为限制条件来迭代求解基函数系数从而得到透镜体的质量分布.通过拟合一个模拟的NFW (Navarro,Frenk and White)透镜系统测试了新方法的可行性,结果表明新方法可以在整体上重构出透镜体的质量分布,并拟合出接近真实的源位置,能够为星系团质量测量提供一套灵活且高效的重构算法.     

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.     

Random gravitational lensing

A simple three-dimensional approach to calculate the probability of observing gravitational lensing, caused by a chance superposition of lens and source, is presented. The formalism is then used to estimate the lensing probabilities in a number of astrophysical situations. As indicated in earlier works, we find that the chance of random gravitational lensing being detected locally is small, and the chance at extra-galactic distances is large.     

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.