Search

With high-quality data from programs like the Hubble Frontier Fields, cluster lensing has reached the point that models are dominated by systematic rather than statistical uncertainties. We introduce a Bayesian framework to quantify systematic effects by determining how different lens modelling choices affect the results. Our framework includes a new two-sample test for quantifying the difference between posterior probability distributions that are sampled by methods like Monte Carlo Markov chains. We use the framework to examine choices related to the selection and treatment of cluster member galaxies in two of the Frontier Field clusters: Abell 2744 and MACS J0416.1–2403. When selecting member galaxies, choices about depth and area affect the models; we find that model results are robust for an I-band magnitude limit of mlim ≥ 22.5 mag and a radial cut of rlim ≥ 90 arcsec (from the centre of the field), although the radial limit likely depends on the spatial extent of lensed images. Mass is typically assigned to galaxies using luminosity/mass scaling relations. We find that the slopes of the scaling relations can have significant effects on lens model parameters but only modest effects on lensing magnifications. Interestingly, scatter in the scaling relations affects the two fields differently. This analysis illustrates how our framework can be used to analyse lens modelling choices and guide future cluster lensing programs.

Cluster lens models are affected by a variety of choices in the lens modelling process. We have begun a programme to develop a systematic error budget for cluster lens modelling. Here, we examine the selection of image constraints as a potential systematic effect. For constraining the mass model, we find that it is more important to have images be spatially distributed around the cluster than to have them distributed in redshift. We also find that some image sets appear to be more important than others in terms of how well they constrain the models; the ‘important’ image sets typically include an image that lies close to a lensing critical curve as well as an image that is relatively isolated from other images (providing constraints in a region that would otherwise lack lensing information). These conclusions can help guide observing programmes that seek follow-up data for candidate lensed images.

We report X-ray observations of the most distant known gravitationally lensed quasar, J0439+1634 at $z=6.52$, which is also a broad absorption line (BAL) quasar, using the XMM-Newton Observatory. With a 130 ks exposure, the quasar is significantly detected as a point source at the optical position with a total of 358$^{+19}_{-19}$ net counts using the EPIC instrument. By fitting a power-law plus Galactic absorption model to the observed spectra, we obtain a spectral slope of $\Gamma=1.45^{+0.10}_{-0.09}$. The derived optical-to-X-ray spectral slope $\alpha_{\rm{ox}}$ is $-2.07^{+0.01}_{-0.01}$, suggesting that the X-ray emission of J0439+1634 is weaker by a factor of 18 than the expectation based on its 2500 Angstrom luminosity and the average $\alpha_{\rm{ox}}$ vs. luminosity relationship. This is the first time that an X-ray weak BAL quasar at $z>6$ has been observed spectroscopically. Its X-ray weakness is consistent with the properties of BAL quasars at lower redshift. By fitting a model including an intrinsic absorption component, we obtain intrinsic column densities of $N_{\rm{H}}=2.8^{+0.7}_{-0.6}\times10^{23}\,\rm{cm}^{-2}$ and $N_{\rm{H}}= 4.3^{+1.8}_{-1.5}\times10^{23}\,\rm{cm}^{-2}$, assuming a fixed $\Gamma$ of 1.9 and a free $\Gamma$, respectively. The intrinsic rest-frame 2--10 keV luminosity is derived as $(9.4-15.1)\times10^{43}\,\rm{erg\,s}^{-1}$, after correcting for lensing magnification ($\mu=51.3$). The absorbed power-law model fitting indicates that J0439+1634 is the highest redshift obscured quasar with a direct measurement of the absorbing column density. The intrinsic high column density absorption can reduce the X-ray luminosity by a factor of $3-7$, which also indicates that this quasar could be a candidate of intrinsically X-ray weak quasar., Comment: 8 pages, 3 figures, 1 table; Accepted for publication in ApJL

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 $\prod_{i=1}^K (5g_i-5)$ images of a single background source, where $g_i$ is the number of point masses in the $i^\text{th}$ plane. This gives asymptotically (for large $g_i$ with $K$ fixed) $5^K$ 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 differential topology, we then show that the 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 the positions of lensed images as well as the structure of the 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., 20 pages, 8 figures. The paper will appear in the Journal of Mathematical Physics

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 two-dimensional, 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 three-dimensional 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-0.4", a non-negligible fraction of the typical 0.3-0.7" residuals of current state-of-the-art models. We also find that median image magnifications can shift by up to 6%, 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% magnification bias still exists in other cases; image position offsets are also improved, but are still present at 0.10"., Comment: 26 pages, 9 figures, 14 tables. Accepted for publication in MNRAS

We report ALMA observations of the dust continuum and {\cii} emission of the host galaxy of J0439+1634, a gravitationally lensed quasar at $z=6.5$. Gravitational lensing boosts the source-plane resolution to $\sim0\farcs15$ $(\sim0.8\text{ kpc})$. The lensing model derived from the ALMA data is consistent with the fiducial model in \citet{fan19} based on {\it HST} imaging. The host galaxy of J0439+1634 can be well-fitted by a S\'ersic profile consistent with an exponential disk, both in the far-infrared (FIR) continuum and the {\cii} emission. The overall magnification is $4.53\pm0.05$ for the continuum and $3.44\pm0.05$ for the {\cii} line. The host galaxy of J0439+1634 is a compact ultra-luminous infrared galaxy, with a total star formation rate (SFR) of $1.56\times10^{3}M_\odot/\text{year}$ after correcting for lensing and an effective radius of $0.74$ kpc. The resolved regions in J0439+1634 follow the ``{\cii} deficit," where the {\cii}-to-FIR ratio decreases with FIR surface brightness. The reconstructed velocity field of J0439+1634 appears to be rotation-like. The maximum line-of-sight rotation velocity of 130 km/s at a radius of 2 kpc. However, our data cannot be fit by an axisymmetric thin rotating disk, and the inclination of the rotation axis, $i$, remains unconstrained. We estimate the dynamical mass of the host galaxy to be $7.9\sin^{-2}(i)\times10^{9}M_\odot$. J0439+1634 is likely to have a high gas-mass fraction and an oversized SMBH compared to local relations. The SFR of J0439+1634 reaches the maximum possible values, and the SFR surface density is close to the highest value seen in any star-forming galaxy currently known in the universe., Comment: 14 pages, 7 figures. Accepted by ApJ

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 modeling 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 modeling. 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, 25 pages, 21 figures. Accepted for publication in MNRAS

Magnetic fields play a pivotal role in the physics of interstellar medium in galaxies, but there are few observational constraints on how they evolve across cosmic time. Spatially resolved synchrotron polarization maps at radio wavelengths reveal well-ordered large-scale magnetic fields in nearby galaxies that are believed to grow from a seed field via a dynamo effect. To directly test and characterize this theory requires magnetic field strength and geometry measurements in cosmologically distant galaxies, which are challenging to obtain due to the limited sensitivity and angular resolution of current radio telescopes. Here, we report the cleanest measurements yet of magnetic fields in a galaxy beyond the local volume, free of the systematics traditional techniques would encounter. By exploiting the scenario where the polarized radio emission from a background source is gravitationally lensed by a foreground galaxy at z = 0.439 using broadband radio polarization data, we detected coherent $\mu$G magnetic fields in the lensing disk galaxy as seen 4.6 Gyrs ago, with similar strength and geometry to local volume galaxies. This is the highest redshift galaxy whose observed coherent magnetic field property is compatible with a mean-field dynamo origin., Comment: 29 pages, 5 figures (including Supplementary Information). Published in Nature Astronomy on August 28, 2017

We report Northern Extended Millimeter Array (NOEMA) CO($J = 3 - 2$) observations of the dusty star-forming galaxy ACT-S\,J020941+001557 at $z = 2.5528$, which was detected as an unresolved source in the Atacama Cosmology Telescope (ACT) equatorial survey. Our spatially resolved spectral line data support the derivation of a gravitational lens model from 37 independent velocity channel maps using a pixel-based algorithm, from which we infer a velocity-dependent magnification factor $\mu \approx 7-22$ with a luminosity-weighted mean $\left\approx 13$. The resulting source-plane reconstruction is consistent with a rotating disk, although other scenarios cannot be ruled out by our data. After correction for lensing, we derive a line luminosity $L^{\prime}_{\rm CO(3-2)}= (5.53\pm 0.69) \times 10^{10}\,{\rm \,K\,km\,s^{-1}\,pc^{2}}$, a cold gas mass $M_{{\rm gas}}= (3.86 \pm 0.33) \times 10^{10}\,M_{\odot}$, a dynamical mass $M_{\rm dyn}\,{\rm sin}^2\,i = 3.9^{+1.8}_{-1.5} \times 10^{10}\,M_{\odot}$, and a gas mass fraction $f_{\rm gas}\,{\rm csc}^2\,i = 1.0^{+0.8}_{-0.4}$. The line brightness temperature ratio of $r_{3,1}\approx 1.6$ relative to a Green Bank Telescope CO($J=1-0$) detection may be elevated by a combination of external heating of molecular clouds, differential lensing, and/or pointing errors., Comment: 8 pages, 5 figures, accepted to ApJ

We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $\mu$m, [CI] 369 $\mu$m, [OI] 146 $\mu$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $\mu^{-1}$ $L_{\odot}$, where $\mu$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $\mu^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} \mu^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17., Comment: 13 pages, 6 figures. Accepted by ApJ

Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy., 8 pages, 4 figures, submitted to ApJL

Gravitational lensing has emerged as a powerful probe of the matter distribution on subgalactic scales, which itself may contain important clues about the fundamental origins and properties of dark matter. Broadly speaking, two different approaches have been taken in the literature to map the small-scale structure of the Universe using strong lensing, with one focused on measuring the position and mass of a small number of discrete massive subhalos appearing close in projection to lensed images, and the other focused on detecting the collective effect of all the small-scale structure between the lensed source and the observer. In this paper, we follow the latter approach and perform a detailed study of the sensitivity of galaxy-scale gravitational lenses to the ensemble properties of small-scale structure. As in some previous studies, we adopt the language of the substructure power spectrum to characterize the statistical properties of the small-scale density field. We present a comprehensive theory that treats lenses with extended sources as well as those with time-dependent compact sources (such as quasars) in a unified framework for the first time. Our approach uses mode functions to provide both computational advantages and insights about couplings between the lens and source. The goal of this paper is to develop the theory and gain the intuition necessary to understand how the sensitivity to the substructure power spectrum depends on the source and lens properties, with the eventual aim of identifying the most promising targets for such studies., 24 pages + appendices, 12 figures. v2: matches published version

The most distant light source that can be observed is the cosmic microwave background. Its cosmological significance is outlined in Sect. 9.1. As in the case of cosmic shear, the lensed power spectrum is the key quantity of interest. However, shear is ill-defined when the background source is the entire sky. Nevertheless, lensing “rearranges” the distribution of photons in a way that alters the observed pattern of temperature fluctuations that characterize cosmic background radiation. The spherical nature of the sky results in a power spectrum that depends on a discrete index, rather than a continuous wavenumber (Sect. 9.2). On the angular scales affected by lensing, the sky can be treated as “flat” in the vicinity of a given point. We use this approximation to derive expressions for the lensed power spectrum (Sect. 9.3.1) and the lens potential (Sect. 9.3.2). Finally, we mention a few applications in Sect. 9.4.

Broadly speaking, there are two types of strong lensing: macrolensing and microlensing. In the former case, the multiple images of the source can be spatially resolved, while in the latter case, they cannot. In this chapter, we investigate microlensing of stars within the Local Group of galaxies. In both Chaps. 2 and 4, we solved the lens equation to find the image positions, with the understanding that they could be observed individually. In this chapter, we assume only that the centroid (i.e., center of light) and total magnification of the images can be measured. We begin by considering microlensing of a point source. The simplest such case is for a lens consisting of a single point mass (Sect. 5.1). We then turn to the case of lensing by several point masses (Sect. 5.2). After formulating the problem in terms of complex variables (Sect. 5.2.1), we specialize to the binary lensing (Sect. 5.2.2) of a source by two stars close together or by a single star orbited by a planet. Extended sources are considered in Sect. 5.3. Putting microlensing theory into practice requires the statistical approach presented in Sect. 5.4. Applications of microlensing are discussed in Sect. 5.5.

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 subhalos through lens modeling, but recent work has suggested that the full population of subhalos 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 \texttt{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\sim 1$ kpc$^{-1}$) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhalos can be directly detected via gravitational imaging., Comment: 8 pages, 10 figures

Before turning to lensing by arbitrary mass distributions, we would do well to remember that the Newtonian deflection angle is off from the actual value by a factor of two. Narrowly speaking, the result of this chapter will be that factor. More broadly, though, understanding how gravitational lensing is an inevitable consequence of general relativity puts the following chapters on solid theoretical ground. We begin in Sect. 3.1 with a review of special relativity, including the mathematics of four-dimensional spacetime. The geodesic equation, whose solution describes the trajectory of a particle moving through curved spacetime, is presented in Sect. 3.2. The special case of a spherically symmetric mass distribution, which is the basis of weak-field gravitational lensing, is treated in Sects. 3.3 and 3.4. We conclude in Sect. 3.5 with a discussion of a homogeneous, isotropic universe, which is the setting for the cosmological applications of lensing given in later chapters.

In some parts of the Universe, hundreds of galaxies are bound by their mutual gravity into enormous structures known as clusters. The lensing properties of galaxy clusters differ from those of individual galaxies for two reasons: a cluster is a hundred to a thousand times more massive than a galaxy, and the inner density profile of a cluster is shallower than isothermal. The large mass creates a non-negligible probability that two or more sources will be lensed by a single cluster, while the shallow density profile allows the creation of radial arcs. Section 7.1 extends our previous discussion of strong lensing to clusters. Moreover, since they are so massive clusters can produced distortions well beyond the strong lensing region. These distortions are too small to be detected individually, but they can be extracted from statistical analysis of thousands of galaxies. Section 7.2 introduces the theory of weak lensing and describes how measurements of shape distortions can be used to constrain lens mass distributions. Section 7.3 then discusses applications of both strong and weak lensing by clusters.

We have seen the beauty and power of gravitational lensing in both the local and distant universe, but exploiting this rich phenomenon requires a quantitative treatment. Many of the most important concepts in lensing can be gleaned from the axisymmetric case, where the mass distribution of the lens has either spherical or cylindrical symmetry in three dimensions, corresponding to circular symmetry in two dimensions. Strictly speaking, we must work within the confines of general relativity, but a Newtonian derivation of the deflection angle is nevertheless instructive (Sect. 2.1) and requires only a slight mathematical tweak to bring the result in line with relativity. Because the radial extent of a lens is much smaller than either the distance from the observer to the lens or from the lens to the source, we can project along the line of sight to arrive at a two-dimensional problem (Sect. 2.2). Using this “thin lens” approximation, we derive the equation that a deflected light ray must obey and discuss how to determine the positions and magnifications of lensed images. We then use this framework to analyze a few simple yet representative lens models (Sect. 2.3). The final three sections derive conditions for multiple imaging by axisymmetric lenses and discuss how to extend our analysis to non-axisymmetric lenses.

The remainder of this book is concerned with applications of gravitational lensing of cosmological relevance. We proceed by considering ever larger mass and distance scales, starting with galaxy-scale strong lensing in the present chapter. The background sources in this context are point-like quasars or extended galaxies. As we will see in Sect. 6.1, the axisymmetric models of Chap. 2 can be readily adapted to the more general setting of Chap. 4. Only such non-axisymmetric mass distributions can be expected to describe observed lens systems. We noted in Chap. 2 that actual galaxies have finite density everywhere and offered the nonsingular isothermal sphere as an example. We extend that discussion to more general lenses in Sect. 6.2. Even all of these refinements are insufficient for understanding the full variety of observed galaxy lenses. For example, we must account for extended sources (Sect. 6.3) and for lenses that cannot be described by smooth density profiles (Sect. 6.4). Using the properties of observed images to infer the mass distribution of the lens is the subject of Sect. 6.5. Applications of the models and methods described in this chapter are outlined in Sect. 6.6.

In the remaining chapters, we apply lensing to cases where no single object, such as a galaxy or cluster of galaxies, is responsible for the deflection of light from a background source. Instead, all matter along the line of sight to the source acts collectively as a lens. Comparing large samples of observed galaxy shapes to what we would expect for particular cosmological models allows us to constrain the properties of the universe itself. The first step is to describe the statistical properties of density fluctuations in three dimensions, from which we can infer the statistics of the two-dimensional convergence field relevant for lensing. Linear perturbation theory, in which density fluctuations are assumed to be small, provides a framework for lensing by large-scale structure (Sect. 8.1.1), even though a nonlinear treatment is needed for a quantitative comparison between theory and observation. Such considerations cannot be ignored by the researcher, but they do not substantially affect the basic approach to the subject. (See Sect. 8.1.2 for a flavor of what is involved.) For a given model, whether linear or nonlinear, the correlation function and power spectrum encapsulate much of the statistics of structure formation (Sects. 8.2 and 8.3). The lensed versions of these quantities can be inferred from galaxy shapes, which are distorted by the “cosmic shear” from matter in the foreground (Sect. 8.4). Understanding this phenomenon sets the stage for the applications presented in the final section (Sect. 8.5).

The degeneracy between disc and halo contributions in spiral galaxy rotation curves makes it difficult to obtain a full understanding of the distribution of baryons and dark matter in disc galaxies like our own Milky Way. Using mock data, we study how constraints on dark matter profiles obtained from kinematics, strong lensing, or a combination of the two are affected by assumptions about the halo model. We compare four different models: spherical isothermal and Navarro-Frenk-White halos, along with spherical and elliptical Burkert halos. For both kinematics and lensing we find examples where different models fit the data well but give enclosed masses that are inconsistent with the true (i.e., input) values. This is especially notable when the input and fit models differ in having cored or cuspy profiles (such as fitting an NFW model when the underlying dark matter distribution follows a different profile). We find that mass biases are more pronounced with lensing than with kinematics, and using both methods can help reduce the bias and provide stronger constraints on the dark matter distributions., Comment: 13 pages, 8 figures, submitted for MNRAS publication Nov 15th 2018

Now that we have had some practice in carrying out lensing calculations and have derived the all-important factor of two, we turn to non-axisymmetric lens models in the context of the thin lens approximation. We show in Sect. 4.1 that the vector counterpart of the deflection angle is given as the gradient of a scalar potential and uses Fermat’s principle to obtain the lens equation. We then solve this equation to find the image positions for a simple but important special case. What follows thereafter is a general discussion of image magnification (Sect. 4.2) and time delay (Sect. 4.3), which can be observed and compared against specified lens models. Some general properties of the lens mapping are described in Sects. 4.4 and 4.5. We derive the conservation of surface brightness in Sect. 4.6, which allows us to consider lensing of spatially extended sources. Mathematical degeneracies in the lens equation, which constrain our ability to compare theory and observation, are discussed in Sect. 4.7. The case of a light ray deflected by multiple lenses (Sect. 4.8) rounds out our presentation of the theory of strong lensing in the weak-field limit.

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 parameterization 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_{\rm{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/sec. We demonstrate that limits on $m_{\rm{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 modeling. Omitting the additional signal from line of sight halos, our constraints underestimate the true power of the method. Assuming cold dark matter, for a low normalization, corresponding the destruction of all subhalos within the host scale radius, we forecast $2\sigma$ bounds on $m_{\rm{hm}}$ (thermal relic mass) of $10^{7.5} \ (5.0)$, $10^{8} \ (3.6)$, and $10^{8.5} \ (2.7) \ M_{\odot} \left(\rm{keV}\right)$ for flux errors of $2\%$, $4\%$, and $8\%$. With a higher normalization, these constraints improve to $10^{7.2} \ (6.6)$, $10^{7.5} \ (5.3) $, and $10^{7.8} \ (4.3) \ 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., Comment: accepted to MNRAS; added a new simulation to illustrate how constraints on free streaming length scale with mass function normalization; 17 pages, 12 figures, 1 table, 5 appendices

Strong gravitational lensing provides a powerful test of Cold Dark Matter (CDM) as it enables the detection and mass measurement of low mass haloes even if they do not contain baryons. Compact lensed sources such as Active Galactic Nuclei (AGN) are particularly sensitive to perturbing subhalos, but their use as a test of CDM has been limited by the small number of systems which have significant radio emission which is extended enough avoid significant lensing by stars in the plane of the lens galaxy, and red enough to be minimally affected by differential dust extinction. Narrow-line emission is a promising alternative as it is also extended and, unlike radio, detectable in virtually all optically selected AGN lenses. We present first results from a WFC3 grism narrow-line survey of lensed quasars, for the quadruply lensed AGN HE0435-1223. Using a forward modelling pipeline which enables us to robustly account for spatial blending, we measure the [OIII] 5007 \AA~ flux ratios of the four images. We find that the [OIII] fluxes and positions are well fit by a simple smooth mass model for the main lens. Our data rule out a $M_{600}>10^{8} (10^{7.2}) M_\odot$ NFW perturber projected within $\sim$1\farcs0 (0\farcs1) arcseconds of each of the lensed images, where $M_{600}$ is the perturber mass within its central 600 pc. The non-detection is broadly consistent with the expectations of $\Lambda$CDM for a single system. The sensitivity achieved demonstrates that powerful limits on the nature of dark matter can be obtained with the analysis of $\sim20$ narrow-line lenses., Comment: Accepted for publication in MNRAS, 15 pages, 8 figures

We present a combined strong and weak lensing analysis of the J085007.6+360428 (J0850) field, which was selected by its high projected concentration of luminous red galaxies and contains the massive cluster Zwicky 1953. Using Subaru/Suprime-Cam $BVR_{c}I_{c}i^{\prime}z^{\prime}$ imaging and MMT/Hectospec spectroscopy, we first perform a weak lensing shear analysis to constrain the mass distribution in this field, including the cluster at $z = 0.3774$ and a smaller foreground halo at $z = 0.2713$. We then add a strong lensing constraint from a multiply-imaged galaxy in the imaging data with a photometric redshift of $z \approx 5.03$. Unlike previous cluster-scale lens analyses, our technique accounts for the full three-dimensional mass structure in the beam, including galaxies along the line of sight. In contrast with past cluster analyses that use only lensed image positions as constraints, we use the full surface brightness distribution of the images. This method predicts that the source galaxy crosses a lensing caustic such that one image is a highly-magnified "fold arc", which could be used to probe the source galaxy's structure at ultra-high spatial resolution ($< 30$ pc). We calculate the mass of the primary cluster to be $\mathrm{M_{vir}} = 2.93_{-0.65}^{+0.71} \times 10^{15}~\mathrm{M_{\odot}}$ with a concentration of $\mathrm{c_{vir}} = 3.46_{-0.59}^{+0.70}$, consistent with the mass-concentration relation of massive clusters at a similar redshift. The large mass of this cluster makes J0850 an excellent field for leveraging lensing magnification to search for high-redshift galaxies, competitive with and complementary to that of well-studied clusters such as the HST Frontier Fields., Accepted for publication in The Astrophysical Journal; 14 pages, 13 figures, 3 tables

Strong gravitational lensing provides an independent measurement of the Hubble parameter ($H_0$). One remaining systematic is a bias from the additional mass due to a galaxy group at the lens redshift or along the sightline. We quantify this bias for more than 20 strong lenses that have well-sampled sightline mass distributions, focusing on the convergence $\kappa$ and shear $\gamma$. In 23% of these fields, a lens group contributes a $\ge$1% convergence bias; in 57%, there is a similarly significant line-of-sight group. For the nine time delay lens systems, $H_0$ is overestimated by 11$^{+3}_{-2}$% on average when groups are ignored. In 67% of fields with total $\kappa \ge$ 0.01, line-of-sight groups contribute $\gtrsim 2\times$ more convergence than do lens groups, indicating that the lens group is not the only important mass. Lens environment affects the ratio of four (quad) to two (double) image systems; all seven quads have lens groups while only three of 10 doubles do, and the highest convergences due to lens groups are in quads. We calibrate the $\gamma$-$\kappa$ relation: $\log(\kappa_{\rm{tot}}) = (1.94 \pm 0.34) \log(\gamma_{\rm{tot}}) + (1.31 \pm 0.49)$ with a rms scatter of 0.34 dex. Shear, which, unlike convergence, can be measured directly from lensed images, can be a poor predictor of $\kappa$; for 19% of our fields, $\kappa$ is $\gtrsim 2\gamma$. Thus, accurate cosmology using strong gravitational lenses requires precise measurement and correction for all significant structures in each lens field., Comment: 34 pages, 11 figures, accepted for publication in ApJ

With a large, unique spectroscopic survey in the fields of 28 galaxy-scale strong gravitational lenses, we identify groups of galaxies in the 26 adequately-sampled fields. Using a group finding algorithm, we find 210 groups with at least five member galaxies; the median number of members is eight. Our sample spans redshifts of 0.04 $\le z_{grp} \le$ 0.76 with a median of 0.31, including 174 groups with $0.1 < z_{grp} < 0.6$. Groups have radial velocity dispersions of 60 $\le \sigma_{grp} \le$ 1200 km s$^{-1}$ with a median of 350 km s$^{-1}$. We also discover a supergroup in field B0712+472 at $z =$ 0.29 consisting of three main groups. We recover groups similar to $\sim$ 85% of those previously reported in these fields within our redshift range of sensitivity and find 187 new groups with at least five members. The properties of our group catalog, specifically 1) the distribution of $\sigma_{grp}$, 2) the fraction of all sample galaxies that are group members, and 3) the fraction of groups with significant substructure, are consistent with those for other catalogs. The distribution of group virial masses agrees well with theoretical expectations. Of the lens galaxies, 12 of 26 (46%) (B1422+231, B1600+434, B2114+022, FBQS J0951+2635, HE0435-1223, HST J14113+5211, MG0751+2716, MGJ1654+1346, PG 1115+080, Q ER 0047-2808, RXJ1131-1231, and WFI J2033-4723) are members of groups with at least five galaxies, and one more (B0712+472) belongs to an additional, visually identified group candidate. There are groups not associated with the lens that still are likely to affect the lens model; in six of 25 (24%) fields (excluding the supergroup), there is at least one massive ($\sigma_{grp} \ge$ 500 km s$^{-1}$) group or group candidate projected within 2$^{\prime}$ of the lens., Comment: 87 pages, 8 figures, a version of this was published in ApJ

We propose a single galaxy gravitational lensing model with a mass density that has a spiral structure. Namely, we extend the arcsine gravitational lens (a truncated singular isothermal elliptical model), adding an additional parameter that controls the amount of spiraling in the structure of the mass density. An important feature of our model is that, even though the mass density is sophisticated, we succeed in integrating the deflection term in closed form using a Gauss hypergeometric function. When the spiraling parameter is set to zero, this reduces to the arcsine lens.

We present ~1" resolution (~2 kpc in the source plane) observations of the CO(1-0), CO(3-2), Halpha, and [N II] lines in the strongly-lensed z=2.26 star-forming galaxy SDSS J0901+1814. We use these observations to constrain the lensing potential of a foreground group of galaxies, and our source-plane reconstructions indicate that SDSS J0901+1814 is a nearly face-on (i~30 degrees) massive disk with r_{1/2}>~4 kpc for its molecular gas. Using our new magnification factors (mu_tot~30), we find that SDSS J0901+1814 has a star formation rate (SFR) of 268^{+63}_{-61} M_sun/yr, M_gas=(1.6^{+0.3}_{-0.2})x10^11x(alpha_CO/4.6) M_sun, and M_star=(9.5^{+3.8}_{-2.8})x10^10 M_sun, which places it on the star-forming galaxy "main sequence." We use our matched high-angular resolution gas and SFR tracers (CO and Halpha, respectively) to perform a spatially resolved (pixel-by-pixel) analysis of SDSS J0901+1814 in terms of the Schmidt-Kennicutt relation. After correcting for the large fraction of obscured star formation (SFR_Halpha/SFR_TIR=0.054^{+0.015}_{-0.014}), we find SDSS J0901+1814 is offset from "normal" star-forming galaxies to higher star formation efficiencies independent of assumptions for the CO-to-H_2 conversion factor. Our mean best-fit index for the Schmidt-Kennicutt relation for SDSS J0901+1814, evaluated with different CO lines and smoothing levels, is n=1.54+/-0.13; however, the index may be affected by gravitational lensing, and we find n=1.24+/-0.02 when analyzing the source-plane reconstructions. While the Schmidt-Kennicutt index largely appears unaffected by which of the two CO transitions we use to trace the molecular gas, the source-plane reconstructions and dynamical modeling suggest that the CO(1-0) emission is more spatially extended than the CO(3-2) emission., ApJ accepted; 40 pages, 24 figures, 8 tables

In gravitational lensing, magnification cross sections characterize the probability that a light source will have magnification greater than some fixed value, which is useful in a variety of applications. The (area) cross section is known to scale as &mu, &minus, 2 for fold caustics and &mu, 2.5 for cusp caustics. We aim to extend the results to higher-order caustic singularities, focusing on the elliptic umbilic, which can be manifested in lensing systems with two or three galaxies. The elliptic umbilic has a caustic surface, and we show that the volume cross section scales as &mu, 2.5 in the two-image region and &mu, 2 in the four-image region, where &mu, is the total unsigned magnification. In both cases our results are supported both numerically and analytically.

This chapter discusses computational and statistical methods for fitting models to strong lens data. It centres on parametric models of point-like lenses but includes extensions to composite models, free-form models, and extended sources. It describes how to use statistical tools including Monte Carlo Markov chains and nested sampling to explore the range of models that are consistent with data. Introduction Strong lensing is a versatile tool for astrophysics that can be used to study the physical properties and environments of lensing galaxies, to dissect the structure of source quasars and galaxies, to constrain cosmological parameters, and much more. Other chapters in this volume review the theory of strong lensing, the status of observations, and the variety of astrophysical applications that result. The goal of this chapter is to outline methods for fitting models to strong lens data. Since modelling is required for most applications of strong lensing, understanding the strengths and weaknesses of the analysis is key for drawing robust conclusions. When discussing methodology, we need to distinguish between point-like and extended images. Point-like images (in a lensed quasar, for example) provide constraints on the potential and its derivatives at discrete positions, which can be described with a modest number of constraint equations or a straightforward χ 2 goodness of fit statistic. Established statistical methods can then be used to find the best fit and explore the range of allowed models. In this case the barrier to entry is low in the sense that fitting basic models does not require tremendous expertise, yet the potential for growth is high in the sense that advanced analyses can combine lensing with other astrophysical probes to draw conclusions that have broad reach. Extended images, by contrast, provide many more pixels of data but require many more free parameters (associated with the unknown shape of the source). Specialized methods must be used to simultaneously fit a mass model for the lens and a light model for the source. For pedagogical purposes, I focus on analysis methods that are applicable to point-like sources but include an overview of methods for modelling extended images (Section 7.5.4).

In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Galaxy-scale strong gravitational lenses have been successfully used to study mass substructures lying close to lensed images of bright background sources. However, in typical galaxy-scale lenses, the strong lensing region only covers a small projected area of the lens's dark matter halo, implying that the vast majority of subhalos cannot be directly detected in lensing observations. In this paper, we point out that this large population of dark satellites can collectively affect gravitational lensing observables, hence possibly allowing their statistical detection. Focusing on the region of the galactic halo outside the strong lensing area, we compute from first principles the statistical properties of perturbations to the gravitational time delay and position of lensed images in the presence of a mass substructure population. We find that in the standard cosmological scenario, the statistics of these lensing observables are well approximated by Gaussian distributions. The formalism developed as part of this calculation is very general and can be applied to any halo geometry and choice of subhalo mass function. Our results significantly reduce the computational cost of including a large substructure population in lens models and enable the use of Bayesian inference techniques to detect and characterize the distributed satellite population of distant lens galaxies., 21 pages + appendices, 7 figures. v2: Some derivations streamlined, extended appendices. Matches version published in PRD

The arrival times, positions, and fluxes of multiple images in strong lens systems can be used to infer the presence of dark subhalos in the deflector, and thus test predictions of cold dark matter models. However, gravitational lensing does not distinguish between perturbations to a smooth gravitational potential arising from baryonic and non-baryonic mass. In this work, we quantify the extent to which the stellar mass distribution of a deflector can reproduce flux ratio and astrometric anomalies typically associated with the presence of a dark matter subhalo. Using Hubble Space Telescope images of nearby galaxies, we simulate strong lens systems with real distributions of stellar mass as they would be observed at redshift $z_d=0.5$. We add a dark matter halo and external shear to account for the smooth dark matter field, omitting dark substructure, and use a Monte Carlo procedure to characterize the distributions of image positions, time delays, and flux ratios for a compact background source of diameter 5 pc. By convolving high-resolution images of real galaxies with a Gaussian PSF, we simulate the most detailed smooth potential one could construct given high quality data, and find scatter in flux ratios of $\approx 10\%$, which we interpret as a typical deviation from a smooth potential caused by large and small scale structure in the lensing galaxy. We demonstrate that the flux ratio anomalies arising from galaxy-scale baryonic structure can be minimized by selecting the most massive and round deflectors, and by simultaneously modeling flux ratio and astrometric data., Comment: 26 pages, 22 figures, submitted to MNRAS

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.

I review methods for modeling gravitational lens systems comprising multiple images of a background source surrounding a foreground galaxy. In a Bayesian framework, the likelihood is driven by the nature of the data, which in turn depends on whether the source is point-like or extended. The prior encodes astrophysical expectations about lens galaxy mass distributions, either through a careful choice of model families, or through an explicit Bayesian prior applied to under-constrained free-form models. We can think about different lens modeling methods in terms of their choices of likelihoods and priors.

[Abridged] 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. 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 Sersic 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., 26 pages, 15 figures; paper I at arXiv:0808.2493; accepted for publication in MNRAS (minor revisions); PDF file with full resolution figures at http://www.sns.ias.edu/~rmandelb/paper2.pdf

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 delta between the light rays emerging from the source, for different lensing configurations. If the lens has a power law profile M propto R^gamma, the angle delta initially increases with lens redshift and then either diverges (for a steep profile gamma1) as zl approaches zs. The scaling with lens mass is roughly delta propto M^(1/(2-gamma)). The results for an NFW profile are qualitatively similar to those for a shallow power law, with delta 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 theta_jet, there is a probability as high as P ~ 0.02-0.07 (theta_jet/{0.5} deg)^-1 that one of the lensed images may be missed (for 2 < zs < 6). Differential absorption within Active Galactic Nuclei 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., minor changes, accepted to MNRAS