Gravitational lensing of extended high-redshift sources by dark matter haloes

High-redshift galaxies and quasi-stellar objects (QSOs) are most likely to be strongly lensed by intervening haloes between the source and the observer. In addition, a large fraction of lensed sources is expected to be seen in the submillimetre region, as a result of the enhanced magnification bias on the steep intrinsic number counts. We extend in three directions Blain's earlier study of this effect.
First, we use a modification of the Press–Schechter mass function and detailed lens models to compute the magnification probability distribution. We compare the magnification cross-sections of populations of singular isothermal spheres and Navarro, Frenk & White (NFW) haloes and find that they are very similar, in contrast to the image-splitting statistics which were recently investigated in other studies. The distinction between the two types of density profile is therefore irrelevant for our purposes.
Secondly, we discuss quantitatively the maximum magnification, μ _max, that can be achieved for extended sources (galaxies) with realistic luminosity profiles, taking into account the possible ellipticity of the lensing potential. We find that μ _max plausibly falls into the range     for sources of     effective radius at redshifts within     .
Thirdly, we apply our model for the lensing magnification to a class of sources following the luminosity evolution typical for a unified scheme of QSO formation. As a result of the peculiar steepness of their intrinsic number counts, we find that the lensed source counts at a fiducial wave length of 850 μm can exceed the unlensed counts by several orders of magnitude at flux densities ≳100 mJy, even with a conservative choice of the maximum magnification.     

Identification of bright lenses from the astrometric observations of gravitational microlensing events

When a source star is gravitationally microlensed by a dark lens, the centroid of the source star image is displaced relative to the position of the unlensed source star, with an elliptical trajectory. Recently, routine astrometric follow-up measurements of these source star image centroid shifts by using high-precision interferometers have been proposed to measure the lens proper motion, which can resolve the lens parameter degeneracy in the photometrically determined Einstein time-scale. When an event is caused by a bright lens, on the other hand, the astrometric shift is affected by the light from the lens, but one cannot identify the existence of the bright lens from the observed trajectory because the resulting trajectory of the bright lens event is also an ellipse. As results, lensing parameters determined from the trajectory differ from those of a dark lens event, causing an incorrect identification of the lens population. In this paper, we show that although the shape and size of the astrometric centroid shift trajectory are changed because of the bright lens, the angular speed of centroid shifts around the apparent position of the unlensed source star is not affected by the lens brightness. Therefore, one can identify the existence of a bright lens and determine its brightness by comparing the lens parameters determined from the 'angular speed curve' with those determined from the trajectory of observed centroid shifts. Once the lens brightness is determined, one can correct for the lens proper motion. As the proposed method provides information about both the lens brightness (dark or bright) and the corrected values of the physical parameters of the lens, one can constrain the nature of massive compact halo objects (MACHOs) significantly better.     

Strong lensing of submillimetre galaxies: a tracer of foreground structure?

The steep source counts and negative K -corrections of bright submillimetre galaxies (SMGs) suggest that a significant fraction of those observed at high flux densities may be gravitationally lensed, and that the lensing objects may often lie at redshifts above 1, where clusters of galaxies are difficult to detect through other means. In this case, follow-up of bright SMGs may be used to identify dense structures along the line-of-sight. Here, we investigate the probability for SMGs to experience strong lensing, using the latest N -body simulations and observed source flux and redshift distributions. We find that almost all high-redshift sources with a flux density above 100 mJy will be lensed if they are not relatively local galaxies. We also give estimates of the fraction of sources experiencing strong lensing as a function of observed flux density. This has implications for planning follow-up observations for bright SMGs discovered in future surveys with the Submillimetre Common-User Bolometer Array 2 and other instruments. The largest uncertainty in these calculations is the maximum allowed lensing amplification, which is dominated by the presently unknown spatial extent of SMGs.     

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.     

Infrared observations of gravitational lensing in Abell 2219 with CIRSI

We present the first detection of a gravitational depletion signal at near-infrared wavelengths, based on deep panoramic images of the cluster Abell 2219 ( z =0.22) taken with the Cambridge Infrared Survey Instrument (CIRSI) at the prime focus of the 4.2-m William Herschel Telescope. Infrared studies of gravitational depletion offer a number of advantages over similar techniques applied at optical wavelengths, and can provide reliable total masses for intermediate-redshift clusters. Using the maximum-likelihood technique developed by Schneider, King & Erben, we detect the gravitational depletion at the 3 confidence level. By modelling the mass distribution as a singular isothermal sphere and ignoring the uncertainty in the unlensed number counts, we find an Einstein radius of (66 per cent confidence limit). This corresponds to a projected velocity dispersion of _v 800 km s¹, in agreement with constraints from strongly lensed features. For a Navarro, Frenk & White mass model, the radial dependence observed indicates a best-fitting halo scalelength of 125 h ¹ kpc. We investigate the uncertainties arising from the observed fluctuations in the unlensed number counts, and show that clustering is the dominant source of error. We extend the maximum-likelihood method to include the effect of incompleteness, and discuss the prospects of further systematic studies of lensing in the near-infrared band.     

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.     

Probing subparsec structure in the Lyman α forest with gravitational microlensing

We present the results of microlens ray-tracing simulations showing the effect of absorbing material between a source quasar and a lensing galaxy in a gravitational lens system. We find that, in addition to brightness fluctuations due to microlensing, the strength of the absorption line relative to the continuum varies with time, with the properties of the variations depending on the structure of the absorbing material. We conclude that such variations will be measurable via ultraviolet spectroscopy of image A of the gravitationally lensed quasar Q2237+0305 if the Lyman α clouds between the quasar and the lensing galaxy possess structure on scales smaller than ∼0.1 pc. The time-scale for the variations is on the order of years to decades, although very short-term variability can occur. While the Lyman α lines may not be accessible at all wavelengths, this approach is applicable to any absorption system, including metal lines.     

CLASS B 1359 + 154: Modelling lensing by a group of galaxies

The recently discovered gravitationally lensed system CLASS B1359 +154 appears to have six detectable images of a single background source at a redshift of 3.235. A group of galaxies acts as the lens, at a redshift of ∼ 1. The present work identifies two distinct, physically plausible image configurations, a 7-image one and a 9-image one. Mass models are constructed corresponding to realizations of these two configurations. Both models call for, in addition to non-singular galaxy-type lenses, a larger scale mass component that resembles the extended dark matter distributions seen in relatively low-redshift galaxy groups. It is presently observationally impossible to study the extended X-ray emission from a group at such a high redshift, hence lensing studies are of some interest. A lensed system with a high image multiplicity does not necessarily admit of a unique lensing interpretation; discrimination is possible with additional observable details (e.g., the image parities, which are uncommon among even the simpler systems).     

Smooth matter and source size in microlensing simulations of gravitationally lensed quasars

Several gravitationally lensed quasars are observed with anomalous magnifications in pairs of images that straddle a critical curve. Simple theoretical arguments suggest that the magnification of these images should be approximately equivalent, whereas one image is observed to be significantly demagnified. Microlensing provides a possible explanation for this discrepancy. There are two key parameters when modelling this effect. The first, the fraction of smooth matter in the lens at the image positions, has been explored by Schechter & Wambsganss. They have shown that the anomalous flux ratio observed in the lensed quasar MG 0414+0534 is a priori a factor of 5 more likely if the assumed smooth matter content in the lens model is increased from 0 to 93 per cent. The second parameter, the size of the emission region, is explored in this paper, and shown to be more significant. We find that the broadening of the magnification probability distributions due to smooth matter content is washed out for source sizes that are predicted by standard models for quasars. We apply our model to the anomalous lensed quasar MG 0414+0534, and find a 95 per cent upper limit of  2.62 × 10¹⁶  h ^−1/2₇₀ ( M /M_⊙)^1/2 cm  on the radius of the I -band emission region. The smooth matter percentage in the lens is unconstrained.     

The Cosmic Lens All-Sky Survey – II. Gravitational lens candidate selection and follow-up

We report the final results of the search for gravitationally lensed flat-spectrum radio sources found in the combination of CLASS (Cosmic Lens All-Sky Survey) and JVAS (Jodrell Bank VLA Astrometric Survey). VLA (Very Large Array) observations of 16 503 sources have been made, resulting in the largest sample of arcsec-scale lens systems available. Contained within the 16 503 sources is a complete sample of 11 685 sources which have two-point spectral indices between 1.4 and 5 GHz flatter than −0.5, and 5-GHz flux densities  ≥30 mJy  . A subset of 8958 sources form a well-defined statistical sample suitable for analysis of the lens statistics. We describe the systematic process by which 149 candidate lensed sources were picked from the statistical sample on the basis of possessing multiple compact components in the 0.2-arcsec resolution VLA maps. Candidates were followed up with 0.05-arcsec resolution MERLIN and 0.003-arcsec VLBA observations at 5 GHz and rejected as lens systems if they failed well-defined surface brightness and/or morphological tests. To illustrate the candidate elimination process, we show examples of sources representative of particular morphologies that have been ruled out by the follow-up observations. 194 additional candidates, not in the well-defined sample, were also followed up. Maps for all the candidates can be found on the World Wide Web at http://www.jb.man.ac.uk/research/gravlens/index.html . We summarize the properties of each of the 22 gravitational lens systems in JVAS/CLASS. 12 are double-image systems, nine are four-image systems and one is a six-image system. 13 constitute a statistically well-defined sample giving a point-source lensing rate of  1:690 ± 190  . The interpretation of the results in terms of the properties of the lensing galaxy population and cosmological parameters will be published elsewhere.     

Exploiting magnification bias in ultradeep submillimetre-wave surveys using ALMA

The surface density of populations of galaxies with steep/shallow source counts is increased/decreased by gravitational lensing magnification. These effects are usually called 'magnification bias' and 'depletion', respectively. However, if sources are demagnified by lensing, then the situation is reversed, and the detectable surface density of galaxies with a shallow source count, as expected at the faintest flux densities, is increased. In general, demagnified sources are difficult to detect and study: exquisite subarcsec angular resolution and surface brightness sensitivity are required, and emission from the lensing object must not dominate the image. These unusual conditions are expected to be satisfied for observations made of the dense swarm of demagnified images that could form very close to the line of sight through the centre of a rich cluster of galaxies using the forthcoming submillimetre-wave Atacama Large Millimeter Array (ALMA) interferometer. The demagnified images of most of the background galaxies lying within about 1 arcmin of a rich cluster of galaxies could be detected in a single 18-arcsec-diameter ALMA field centred on the cluster core, providing an effective increase in the ALMA field of view. This technique could allow a representative sample of faint,  10–100 μJy  submillimetre galaxies to be detected several times more rapidly than in a blank field.     

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.     

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.     

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.     

Search for dark matter using the phenomenon of strong gravitational lensing

We provide a brief overview of the methods for estimating the dark matter content in the Universe based on the phenomenon of strong gravitational lensing—the method of macroimage flux ratio anomalies and the method based on an analysis of the probability distribution for the magnification of macroimages due to microlensing events. Both methods require the specification of a macrolens model, knowledge of the spatial structure or at least the effective size of the source, and numerical simulation of microlensing events followed by a comparison of the simulation results with observational data. Using the quadruply lensed quasar Q2237+0305 as an example, we show the effect of the spatial source structure on the shape of the magnification probability distribution in microlensing events. We also point out the need to take into account the contribution from the intrinsic quasar variability to the observed light curve and to develop a physically justified algorithm to fit the observational data. For the first time, based on all the available observations of Q2237+0305, we have constructed the magnification probability histograms for all four macroimages. We analyze the possibility of using them to estimate the content of continuously distributed (dark) matter in the galaxy Q2237+0305 at a distance from its nucleus that corresponds to the macroimage locations.     

On the streaming motions of haloes and galaxies

Little is known about the statistics of gravitationally lensed quasars at large (7–30 arcsec) image separations, which probe masses on the scale of galaxy clusters. We have carried out a survey for gravitationally lensed objects, among sources in the FIRST 20-cm radio survey that have unresolved optical counterparts in the digitizations of the Palomar Observatory Sky Survey. From the statistics of ongoing surveys that search for quasars among FIRST sources, we estimate that there are about 9100 quasars in this source sample, making this one of the largest lensing surveys to date. Using broad-band imaging, we have isolated all objects with double radio components separated by 5–30 arcsec that have unresolved optical counterparts with similar BVI colours. Our criteria for similar colours conservatively allow for observational error and for colour variations due to time delays between lensed images. Spectroscopy of these candidates shows that none of the pairs are lensed quasars. This sets an upper limit (95 per cent confidence) on the lensing fraction in this survey of 3.3×10⁻⁴, assuming 9100 quasars. Although the source redshift distribution is poorly known, a rough calculation of the expected lensing frequency and the detection efficiencies and biases suggests that simple theoretical expectations are of the same order of magnitude as our observational upper limit. Our procedure is novel in that our exhaustive search for lensed objects does not require prior identification of the quasars in the sample as such. Characterization of the FIRST-selected quasar population will enable use of our result to constrain quantitatively the mass properties of clusters.     

Search for exotic matter from gravitational microlensing observations of stars

We consider small-scale spheroidal clusters of weakly interacting massive particles in our Galaxy as non-compact gravitational microlenses and predict the appearance of caustics in the plane of a lensed source. The crossing of these caustics by a lensed star can produce a large variety of light curves, including some observed in actual microlensing events that have been interpreted as manifestations of binary gravitational lenses. We consider also observable effects during the gravitational microlensing of stars of non-zero angular size with a given brightness distribution across their disks by such an exotic objects as natural wormholes and objects whose space-time environment is described with the NUT metric. We demonstrate that, under certain conditions, the microlensing light curves, chromatic and polarizational effects due to the properties of the lens and the star disk brightness distributions can differ considerably from those observed for a Schwarzschild gravitational lens, so that their analysis can facilitate the identification of such objects.     

Cluster lenses

Clusters of galaxies are the most recently assembled, massive, bound structures in the Universe. As predicted by General Relativity, given their masses, clusters strongly deform space-time in their vicinity. Clusters act as some of the most powerful gravitational lenses in the Universe. Light rays traversing through clusters from distant sources are hence deflected, and the resulting images of these distant objects therefore appear distorted and magnified. Lensing by clusters occurs in two regimes, each with unique observational signatures. The strong lensing regime is characterized by effects readily seen by eye, namely, the production of giant arcs, multiple images, and arclets. The weak lensing regime is characterized by small deformations in the shapes of background galaxies only detectable statistically. Cluster lenses have been exploited successfully to address several important current questions in cosmology: (i) the study of the lens(es)—understanding cluster mass distributions and issues pertaining to cluster formation and evolution, as well as constraining the nature of dark matter; (ii) the study of the lensed objects—probing the properties of the background lensed galaxy population—which is statistically at higher redshifts and of lower intrinsic luminosity thus enabling the probing of galaxy formation at the earliest times right up to the Dark Ages; and (iii) the study of the geometry of the Universe—as the strength of lensing depends on the ratios of angular diameter distances between the lens, source and observer, lens deflections are sensitive to the value of cosmological parameters and offer a powerful geometric tool to probe Dark Energy. In this review, we present the basics of cluster lensing and provide a current status report of the field.     

General polytropic spheres as gravitational lenses

We propose hydrostatic polytropic spheres governed by the Lane-Emden equation (LEE) of index $n$ as a novel set of physical models for axially averaged gravitational lenses anywhere in the Universe, alternative to the familiar singular isothermal sphere (SIS) and the Navarro–Frenk–White (NFW) profile, as such general polytropic spheres are conceptually simple, versatile in representing a series of equations of state, and able to address both the inner core and cusp features. As LEE is nonlinear, there exist several distinct classes of LEE solutions to serve as physical lens models. With a few scaling parameters, the complete problem can be readily reconstructed with full physical dimensions. A given mass density profile satisfying LEE produces lensing effects that are solely determined by a dimensionless parameter $q$ which contains geometric and kinematic information about the source-lens-observer system. The lens mapping and tangential shear or distortion profile are derived, first analytically for special cases and then asymptotically at the outskirts or near the edge of the lens. Numerical procedures for calculating full lensing profiles of a general lens are developed. Our results include the analytical “singular polytropic sphere” (SPS) profile which generalizes the SIS model and may outperform the latter in modeling dark matter halos among others. We further point out that dynamic models of general polytropic spheres in self-similar evolution can serve as several broad classes of gravitational lenses and produce time-dependent lensing effects slow or fast depending on the pertinent time scales. Astrophysical sources that can be lensed include electromagnetic wave sources in the entire frequency band, gravitational wave sources in the entire frequency band, gravitons even possibly with finite masses, neutrino sources of three different types, neutron sources, and ultra high energy cosmic rays (UHECRs) of electrically charged particles which can also interact with magnetic fields. We discuss and elabrate applications to dark matter halos, hypermassive black holes and supermassive black holes in the entire Universe including the early Universe, magnetized supermassive stars, static and dynamically evolving spherical and cylindrical lenses in contexts of astrophysics and cosmology.