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.     

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.     

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.     

Observations of lensed relativistic jets as a tool for constraining lens galaxy parameters

The possibility of using lensed relativistic jets on very small angular scales to construct proper models of spiral lens galaxies and to independently determine the Hubble constant is considered. The system B0218+357 is used as an example to illustrate that there exists a great choice of model parameters adequately reproducing its observed large-scale properties but leading to a significant spread in the Hubble constant. The jet image position angle is suggested as an additional parameter that allows the range of models under consideration to be limited. It is shown that themodels for which the jet image position angles differ by at least 40° can be distinguished between themselves during observations on very small angular scales. The possibility of observing the geometric properties of lensed relativistic jets and measuring the superluminal velocities of knot images on time scales of several months with very long baseline space interferometers is discussed.     

Degeneracies and scaling relations in general power-law models for gravitational lenses

The time-delay in gravitational lenses can be used to derive the Hubble constant in a relatively simple way. The results of this method are less dependent on astrophysical assumptions than in many other methods. For systems with accurately measured positions and time-delays, the most important uncertainty is related to the mass model used. Simple parametric models like isothermal ellipsoidal mass distributions seem to provide consistent results with a reasonably small scatter when applied to several lens systems. We discuss a family of models with a separable radial power law and an arbitrary angular dependence for the potential   ψ = r ^β F ( θ )  . Isothermal potentials are a special case of these models with   β =1  . An additional external shear is used to take into account perturbations from other galaxies. Using a simple linear formalism for quadruple lenses, we can derive H ₀ as a function of the observables and the shear. If the latter is fixed, the result depends on the assumed power-law exponent according to   H ₀∝(2- β )/ β   . The effect of external shear is quantified by introducing a 'critical shear' γ _c as a measure for the amount of shear that changes the result significantly. The analysis shows that in the general case H ₀ and γ _c do not depend on the position of the lens galaxy. Spherical lens models with images close to the Einstein radius with fitted external shear differ by a factor of   β /2  from shearless models, leading to   H ₀∝2- β   in this case. We discuss these results and compare them with numerical models for a number of real lens systems.     

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.     

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

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 anatomy of a quadruply imaged gravitational lens system

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

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.     

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.     

Non-parametric strong lens inversion of SDSS J1004+4112

In this article, we study the well-known strong lensing system SDSS J1004+4112. Not only does it host a large-separation lensed quasar with measured time-delay information, but several other lensed galaxies have been identified as well. A previously developed strong lens inversion procedure that is designed to handle a wide variety of constraints is applied to this lensing system and compared to results reported in other works. Without the inclusion of a tentative central image of one of the galaxies as a constraint, we find that the model recovered by the other constraints indeed predicts an image at that location. An inversion which includes the central image provides tighter constraints on the shape of the central part of the mass map. The resulting model also predicts a central image of a second galaxy where indeed an object is visible in the available Advanced Camera for Surveys images. We find masses of  2.5 × 10¹³  and  6.1 × 10¹³ M_⊙  within a radius of 60 and 110 kpc, respectively, confirming the results from other authors. The resulting mass map is compatible with an elliptical generalization of a projected NFW profile, with   r _s= 58⁺²¹₋₁₃  arcsec and   c _vir= 3.91 ± 0.74  . The orientation of the elliptical NFW profile closely follows the orientation of the central cluster galaxy and the overall distribution of cluster members.     

Chandra observations of the galaxy cluster Abell 1835

We present the analysis of 30 ks of Chandra observations of the galaxy cluster Abell 1835. Overall, the X-ray image shows a relaxed morphology, although we detect substructure in the inner 30-kpc radius. Spectral analysis shows a steep drop in the X-ray gas temperature from ∼12 keV in the outer regions of the cluster to ∼4 keV in the core. The Chandra data provide tight constraints on the gravitational potential of the cluster which can be parametrized by a Navarro, Frenk & White model. The X-ray data allow us to measure the X-ray gas mass fraction as a function of radius, leading to a determination of the cosmic matter density of
   
. The projected mass within a radius of ∼150 kpc implied by the presence of gravitationally lensed arcs in the cluster is in good agreement with the mass models preferred by the Chandra data. We find a radiative cooling time of the X-ray gas in the centre of Abell 1835 of about
   
. Cooling-flow model fits to the Chandra spectrum and a deprojection analysis of the Chandra image both indicate the presence of a young cooling flow (∼     with an integrated mass deposition rate of     within a radius of 30 kpc. We discuss the implications of our results in the light of recent Reflection Grating Spectrograph (RGS) observations of Abell 1835 with XMM-Newton .     

A strategy for finding gravitationally lensed distant supernovae

Distant Type Ia and II supernovae (SNe) can serve as valuable probes of the history of the cosmic expansion and star formation, and provide important information on their progenitor models. At present, however, there are few observational constraints on the abundance of SNe at high redshifts. A major science driver for the Next Generation Space Telescope is the study of such very distant SNe. In this paper we discuss strategies for finding and counting distant SNe by using repeat imaging of supercritical intermediate redshift clusters whose mass distributions are well constrained via modelling of strongly lensed features. For a variety of different models for the star formation history and supernova progenitors, we estimate the likelihood of detecting lensed SNe as a function of their redshift. In the case of a survey conducted with Hubble Space Telescope ( HST ), we predict a high probability of seeing a supernova in a single return visit with either Wide Field Planetary Camera 2 or Advanced Camera for Surveys, and a much higher probability of detecting examples with     in the lensed case. Most events would represent magnified SNe II at     and a fraction will be more distant examples. We discuss various ways to classify such events using ground-based infrared photometry. We demonstrate an application of the method using the HST archival data and discuss the case of a possible event found in the rich cluster AC 114     .     

A new gravitational lens from the MUSCLES survey: ULAS J082016.1+081216

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

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.     

Submillimetre detection of the z=2.83 Lyman-break galaxy, Westphal–MM8, and implications for SCUBA2

We present confusion-limited submillimetre (submm) observations with the Submillimetre Common-User Bolometer Array (SCUBA) camera on the James Clerk Maxwell Telescope of the   z = 2.83  Lyman-break galaxy (LBG), Westphal–MM8, reaching an 850 μm sensitivity even greater than that achieved in the SCUBA map of the Hubble Deep Field region. The detection of MM8  ( S _850 μm= 1.98 ± 0.48 mJy)  , along with the literature submm detections of lensed LBGs, suggests that the LBG population may contribute significantly to the source counts of submm-selected galaxies in the 1–2 mJy regime. Additionally, submm-luminous LBGs are a viable progenitor population for the recently discovered evolved galaxies at   z ∼ 2–2.5  . These observations represent an important baseline for SCUBA2 observations which will regularly map large regions of the sky to this depth.     

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.     

The Chang–Refsdal lens revisited

This paper provides a complete theoretical treatment of the point-mass lens perturbed by constant external shear, often called the Chang–Refsdal lens. We show that simple invariants exist for the products of the (complex) positions of the four images, as well as moment sums of their signed magnifications. The image topographies and equations of the caustics and critical curves are also studied. We derive the fully analytic expressions for pre-caustics, which are the loci of non-critical points that map to the caustics under the lens mapping. They constitute boundaries of the region in the image domain that maps on to the interior of the caustics. The areas under the critical curves, caustics and pre-caustics are all evaluated, which enables us to calculate the mean magnification of the source within the caustics. Additionally, the exact analytic expression for the magnification distribution for the source in the triangular caustics is derived, as well as a useful approximate expression. Finally, we find that the Chang–Refsdal lens with additional convergence greater than unity (the 'overfocusing case') can exhibit third-order critical behaviour, if the 'reduced shear' is exactly equal to     , and that the number of images for N -point masses with non-zero constant shear cannot be greater than  5 N − 1  .