Search

We consider gravitational lensing of a background source by a finite system of point-masses. The problem of determining the maximum possible number of lensed images has been completely resolved in the single-plane setting (where the point masses all reside in a single lens plane), but this problem remains open in the multiplane setting. We construct examples of K-plane point-mass gravitational lens ensembles that produce ∏ i = 1 K (5 g i − 5) images of a single background source, where gi is the number of point masses in the ith plane. This gives asymptotically (for large gi with K fixed) 5K times the minimal number of lensed images. Our construction uses Rhie's single-plane examples and a structured parameter-rescaling algorithm to produce preliminary systems of equations with the desired number of solutions. Utilizing the stability principle from the differential topology, we then show that preliminary (nonphysical) examples can be perturbed to produce physically meaningful examples while preserving the number of solutions. We provide numerical simulations illustrating the result of our construction, including positions of lensed images and the structure of critical curves and caustics. We observe an interesting "caustic of multiplicity" phenomenon that occurs in the nonphysical case and has a noticeable effect on the caustic structure in the physically meaningful perturbative case. [ABSTRACT FROM AUTHOR]

The Hubble Frontier Fields data, along with multiple data sets obtained by other telescopes, have provided some of the most extensive constraints on cluster lenses to date. Multiple lens modelling teams analyzed the fields and made public a number of deliverables. By comparing these results, we can then undertake a unique and vital test of the state of cluster lens modelling. Specifically, we see how well the different teams can reproduce similar magnifications and mass profiles. We find that the circularly averaged mass profiles of the fields are remarkably constrained (scatter |$\lt 5{{\ \rm per\ cent}}$|⁠) at distances of 1 arcmin from the cluster core, yet magnifications can vary significantly. Averaged across the six fields, we find a bias of −6 per  cent (−17 per cent) and a scatter of ∼40 per cent (∼65 per cent) at a modest magnification of 3 (10). Statistical errors reported by individual teams are often significantly smaller than the differences among all the teams, indicating the importance of continued systematics studies in cluster lensing. [ABSTRACT FROM AUTHOR]

Cluster lensing has become an important tool in the search for high-redshift galaxies through its ability to magnify sources. In order to determine the intrinsic properties of these galaxies, lensing mass models must be constructed to determine the magnification of the images. These models are traditionally 2D, focusing on the mass within the cluster and either ignoring or approximating any contribution from line-of-sight galaxies. In this paper, we present the first full set of 3D mass models of the six Hubble Frontier Fields and use them to test for systematic biases in magnifications due to using the traditional 2D approach. We find that omitting foreground or background galaxies causes image position offsets between 0.1 and 0.4 arcsec, a non-negligible fraction of the typical 0.3–0.7 arcsec residuals of current state-of-the-art models. We also find that median image magnifications can shift by up to 6 per cent, though it is dependent on the field. This can be alleviated in some cases by approximating the mass in the lensing plane, but a 5 per cent magnification bias still exists in other cases; image position offsets are also improved, but are still present at 0.10 arcsec. [ABSTRACT FROM AUTHOR]

In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhaloes through lens modelling, but recent work has suggested that the full population of subhaloes could be probed using a power spectrum analysis. In this paper, we quantify the power spectrum of small-scale structure in simulated galaxies, with the goal of understanding theoretical predictions and setting the stage for using measurements of the power spectrum to test dark matter models. We use a sample of simulated galaxies generated from the galacticus semi-analytic model to determine the power spectrum distribution first in the CDM paradigm and then in a warm dark matter scenario. We find that a measurement of the slope and amplitude of the power spectrum on galaxy strong lensing scales (k ∼ 1 kpc−1) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhaloes can be directly detected via gravitational imaging. [ABSTRACT FROM AUTHOR]

The free streaming length of dark matter particles determines the abundance of structure on sub-galactic scales. We present a statistical technique, amendable to any parametrization of subhalo density profile and mass function, to probe dark matter on these scales with quadrupole image lenses. We consider a warm dark matter particle with a mass function characterized by a normalization and free streaming scale m hm. We forecast bounds on dark matter warmth for 120–180 lenses, attainable with future surveys, at typical lens (source) redshifts of 0.5 (1.5) in early-type galaxies with velocity dispersions of  220–270 km s-1. We demonstrate that limits on m hm deteriorate rapidly with increasing uncertainty in image fluxes, underscoring the importance of precise measurements and accurate lens models. For our forecasts, we assume the deflectors in the lens sample do not exhibit complex morphologies, so we neglect systematic errors in their modelling. Omitting the additional signal from line-of-sight haloes, our constraints underestimate the true power of the method. Assuming cold dark matter, for a low normalization, corresponding to the destruction of all subhaloes within the host scale radius, we forecast 2σ bounds on m hm (thermal relic mass) of 107.5 (5.0), 108 (3.6), and |$10^{8.5} \ (2.7) \ {\rm M}_{\odot } \left(\rm {keV}\right)$| for flux errors of |$2{{\ \rm per\ cent}}$|⁠, |$4{{\ \rm per\ cent}}$|⁠, and |$8{{\ \rm per\ cent}}$|⁠. With a higher normalization, these constraints improve to 107.2 (6.6), 107.5 (5.3), and |$10^{7.8} \ (4.3) \ {\rm M}_{\odot } \left(\rm {keV}\right)$| with 120 systems. We are also able to measure the normalization of the mass function, which has implications for baryonic feedback models and tidal stripping. [ABSTRACT FROM AUTHOR]

Gravitational lens modelling of spatially resolved sources is a challenging inverse problem with many observational constraints and model parameters. We examine established pixel-based source reconstruction algorithms for de-lensing the source and constraining lens model parameters. Using test data for four canonical lens configurations, we explore statistical and systematic uncertainties associated with gridding, source regularization, interpolation errors, noise, and telescope pointing. Specifically, we compare two gridding schemes in the source plane: a fully adaptive grid that follows the lens mapping but is irregular, and an adaptive Cartesian grid. We also consider regularization schemes that minimize derivatives of the source (using two finite difference methods) and introduce a scheme that minimizes deviations from an analytic source profile. Careful choice of gridding and regularization can reduce ‘discreteness noise’ in the χ2 surface that is inherent in the pixel-based methodology. With a gridded source, some degree of interpolation is unavoidable, and errors due to interpolation need to be taken into account (especially for high signal-to-noise data). Different realizations of the noise and telescope pointing lead to slightly different values for lens model parameters, and the scatter between different ‘observations’ can be comparable to or larger than the model uncertainties themselves. The same effects create scatter in the lensing magnification at the level of a few per cent for a peak signal-to-noise ratio of 10, which decreases as the data quality improves. [ABSTRACT FROM PUBLISHER]

In strong gravitational lens systems, the light bending is usually dominated by one main galaxy, but may be affected by other mass along the line of sight (LOS). Shear and convergence can be used to approximate the contributions from less significant perturbers (e.g. those that are projected far from the lens or have a small mass), but higher order effects need to be included for objects that are closer or more massive. We develop a framework for multiplane lensing that can handle an arbitrary combination of tidal planes treated with shear and convergence and planes treated exactly (i.e. including higher order terms). This framework addresses all of the traditional lensing observables including image positions, fluxes, and time delays to facilitate lens modelling that includes the non-linear effects due to mass along the LOS. It balances accuracy (accounting for higher order terms when necessary) with efficiency (compressing all other LOS effects into a set of matrices that can be calculated up front and cached for lens modelling). We identify a generalized multiplane mass sheet degeneracy, in which the effective shear and convergence are sums over the lensing planes with specific, redshift-dependent weighting factors. [ABSTRACT FROM PUBLISHER]

ABSTRACT We investigate the properties of dark matter substructure in the gravitational lens HE 0435−1223 ( zl=0.455) via its effects on the positions and flux ratios of the quadruply imaged background quasar ( zs= 1.689). We start with a smooth mass model, add individual, truncated isothermal clumps near the lensed images and use the Bayesian evidence to compare the quality of different models. Compared with smooth models, models with at least one clump near image A are strongly favoured. The mass of this clump within its Einstein radius is log10( MAEin) = 7.65+0.87− 0.84 (in units of h−170 M⊙). The Bayesian evidence provides weaker support for a second clump near image B, with log10( MBEin) = 6.55+1.01− 1.51. We also examine models with a full population of substructure, and find the mass fraction in substructure at the Einstein radius to be fsub> rsim 0.000 77, assuming the total clump masses follow a mass function d N/d M∝ M−1.9 over the range M= 107-1010 M⊙. Few-clump and population models produce similar Bayesian evidence values, so neither type of model is objectively favoured. [ABSTRACT FROM AUTHOR]

In this second paper, we develop an analytical theory of quasi-equatorial lensing by Kerr black holes. In this setting we solve perturbatively our general lens equation with displacement given in Paper I, going beyond weak-deflection Kerr lensing to third order in our expansion parameter [variant_greek_epsilon], which is the ratio of the angular gravitational radius to the angular Einstein radius. We obtain new formulas and results for the bending angle, image positions, image magnifications, total unsigned magnification, and centroid, all to third order in [variant_greek_epsilon] and including the displacement. New results on the time delay between images are also given to second order in [variant_greek_epsilon], again including displacement. For all lensing observables we show that the displacement begins to appear only at second order in [variant_greek_epsilon]. When there is no spin, we obtain new results on the lensing observables for Schwarzschild lensing with displacement. [ABSTRACT FROM AUTHOR]

We develop a unified, analytic framework for gravitational lensing by Kerr black holes. In this first paper, we present a new, general lens equation and magnification formula governing lensing by a compact object. Our lens equation assumes that the source and observer are in the asymptotically flat region and does not require a small angle approximation. Furthermore, it takes into account the displacement that occurs when the light ray's tangent lines at the source and observer do not meet on the lens plane. We then explore our lens equation in the case when the compact object is a Kerr black hole. Specifically, we give an explicit expression for the displacement when the observer is in the equatorial plane of the Kerr black hole as well as for the case of spherical symmetry. [ABSTRACT FROM AUTHOR]

Nested sampling has emerged as a valuable tool for Bayesian analysis, in particular for determining the Bayesian evidence. The method is based on a specific type of random sampling of the likelihood function and prior volume of the parameter space. I study the statistical uncertainty in the evidence computed with nested sampling. I examine the uncertainty estimator from Skilling and introduce a new estimator based on a detailed analysis of the statistical properties of nested sampling. Both perform well in test cases and make it possible to obtain the statistical uncertainty in the evidence with no additional computational cost. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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

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 , the angle δ initially increases with lens redshift and then either diverges (for a steep profile ), remains constant (for an isothermal profile ), or vanishes (for a shallow profile ) as . The scaling with lens mass is roughly . 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 , there is a probability as high as that one of the lensed images may be missed (for ). 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We study microlensing of the faint images that form close to the centres of strong gravitational lens galaxies. These central images, which have finally begun to yield to observations, naturally appear in dense stellar fields and may be particularly sensitive to fine granularity in the mass distribution. The microlensing magnification maps for overfocused (i.e. demagnified) images differ strikingly from those for magnified images. In particular, the familiar ‘fold’ and ‘cusp’ features of maps for magnified images are only present for certain values of the fraction f of the surface mass density contained in stars. For central images, the dispersion in microlensing magnifications is generally larger than for normal (minimum and saddle) images, especially when the source is comparable to or larger than the stellar Einstein radius. The dispersion depends in a complicated way on f; this behaviour may hold the key to using microlensing as a probe of the relative densities of stars and dark matter in the cores of distant galaxies. Quantitatively, we predict that the central image C in PMN J1632−0033 has a magnification dispersion of 0.6 mag for , or 0.3 mag for . For comparison, the dispersions are 0.5–0.6 mag for image B and 0.05–0.1 mag for image A, if ; and just 0.1 mag for B and 0.008 mag for A if . (The dispersions can be extrapolated to larger sources sizes as .) Thus, central images are more susceptible than other lensed images to microlensing and hence good probes for measuring source sizes. [ABSTRACT FROM AUTHOR]

Microlensing promises to be a powerful tool for studying distant galaxies and quasars. As the data and models improve, there are systematic effects that need to be explored. Quasar continuum and broad-line regions may respond differently to microlensing due to their different sizes; to understand this effect, we study microlensing of finite sources by a mass function of stars. We find that microlensing is insensitive to the slope of the mass function but does depend on the mass range. For negative-parity images, diluting the stellar population with dark matter increases the magnification dispersion for small sources and decreases it for large sources. This implies that the quasar continuum and broad-line regions may experience very different microlensing in negative-parity lensed images. We confirm earlier conclusions that the surface brightness profile and geometry of the source have little effect on microlensing. Finally, we consider non-circular sources. We show that elliptical sources that are aligned with the direction of shear have larger magnification dispersions than sources with perpendicular alignment, an effect that becomes more prominent as the ellipticity increases. Elongated sources can lead to more rapid variability than circular sources, which raises the prospect of using microlensing to probe source shape. [ABSTRACT FROM AUTHOR]

Flux ratio ‘anomalies’ in quadruply imaged gravitational lenses can be explained with galactic substructure of the sort predicted by ΛCDM, but the strength and uniqueness of that hypothesis need to be further assessed. A good way to do that is to use the physical scale associated with the size of the source quasar, and its dependence on wavelength. We develop a toy model to study finite source effects in substructure lensing. Treating substructure as a singular isothermal sphere allows us to compute the images of a finite source analytically, and then to explore how the image configurations and magnifications depend on source position and size. Although simplified, our model yields instructive general principles: image positions and magnifications are basically independent of source size until the source is large enough to intersect a substructure caustic; even sources that are much larger than the substructure Einstein radius can be perturbed at a detectable level; and, most importantly, there is a tremendous amount to be learned from comparing image positions and magnifications at wavelengths that correspond to different source sizes. In a separate analysis, we carefully study four observed radio lenses to determine which of the images are anomalous. In B0712+472, the evidence for a radio flux ratio anomaly is marginal, but if the anomaly is real then image C is probably the culprit. In B1422+231, the anomaly is in image A. Interestingly, B2045+265 and B1555+375 both appear to have two anomalous images. Coincidentally, in each system one of the anomalies is in image C, and the other is in either image A or image B (both possibilities lead to acceptable models). It remains to be seen whether ΛCDM predicts enough substructure to explain multiple anomalies in multiple lenses. When we finally join our modelling results and substructure theory, we obtain lower bounds on the masses of the substructures responsible for the observed anomalies. The mass bounds are broadly consistent with expectations for ΛCDM. Perhaps more importantly, we outline various systematic effects in the mass bounds; poor knowledge of whether the substructure lies within the main lens galaxy or elsewhere along the line of sight appears to be the dominant systematic. [ABSTRACT FROM AUTHOR]

We study the impact of lens galaxy environments on the image separation distribution of lensed quasars. We account for both environmental convergence and shear, using a joint distribution derived from galaxy formation models calibrated by galaxy–galaxy lensing data and number counts of massive elliptical galaxies. We find that the external field enhances lensing probabilities, particularly, at large image separations; the increase is ∼30 per cent at and ∼200 per cent at , when we adopt a power-law source luminosity function . The enhancement is mainly driven by convergence, which boosts both the image separation and magnification bias (for a fixed lens galaxy mass). These effects have been neglected in previous studies of lens statistics. Turning the problem around, we derive the posterior convergence and shear distributions and point out that they are strong functions of image separation; lens systems with larger image separations are more likely to lie in dense environments. [ABSTRACT FROM AUTHOR]

It has been known for over a decade that many four-image gravitational lenses exhibit anomalous radio flux ratios. These anomalies can be explained by adding a clumpy cold dark matter (CDM) component to the background galactic potential of the lens. As an alternative, Evans & Witt (2003) recently suggested that smooth multipole perturbations provide a reasonable alternative to CDM substructure in some but not all cases. We generalize their method in two ways so as to determine whether multipole models can explain highly anomalous systems. We carry the multipole expansion to higher order, and also include external tidal shear as a free parameter. Fitting for the shear proves crucial to finding a physical (positive-definite density) model. For B1422+231, working to order (and including shear) yields a model that is physical but implausible. Going to higher order reduces global departures from ellipticity, but at the cost of introducing small-scale wiggles in proximity to the bright images. These localized undulations are more pronounced in B2045+265, where multipoles are required to smooth out large-scale deviations from elliptical symmetry. Such modes surely cannot be taken at face value; they must indicate that the models are trying to reproduce some other sort of structure. Our formalism naturally finds models that fit the data exactly, but we use B0712+472 to show that measurement uncertainties have little effect on our results. Finally, we consider the system B1933+503, where two sources are lensed by the same foreground galaxy. The additional constraints provided by the images of the second source render the multipole model unphysical. We conclude that external shear must be taken into account to obtain plausible models, and that a purely smooth angular structure for the lens galaxy does not provide a viable alternative to the prevailing CDM clump hypothesis. [ABSTRACT FROM AUTHOR]