We investigate how strong gravitational lensing can test contemporary models of massive elliptical (ME) galaxy formation, by combining a traditional decomposition of their visible stellar distribution with a lensing analysis of their mass distribution. As a proof of concept, we study a sample of three ME lenses, observing that all are composed of two distinct baryonic structures, a `red' central bulge surrounded by an extended envelope of stellar material. Whilst these two components look photometrically similar, their distinct lensing effects permit a clean decomposition of their mass structure. This allows us to infer two key pieces of information about each lens galaxy: (i) the stellar mass distribution (without invoking stellar populations models) and (ii) the inner dark matter halo mass. We argue that these two measurements are crucial to testing models of ME formation, as the stellar mass profile provides a diagnostic of baryonic accretion and feedback whilst the dark matter mass places each galaxy in the context of LCDM large scale structure formation. We also detect large rotational offsets between the two stellar components and a lopsidedness in their outer mass distributions, which hold further information on the evolution of each ME. Finally, we discuss how this approach can be extended to galaxies of all Hubble types and what implication our results have for studies of strong gravitational lensing.

Using ALMA, we report high angular-resolution observations of the redshift z=3.63 galaxy, G09v1.97, one of the most luminous strongly lensed galaxies discovered by the H-ATLAS survey. We present 0"2-0"4 resolution images of the rest-frame 188 and 419$\mu$m dust continuum and the CO(6-5), H2O(211-202) and J=2 H2O+ line emission. We also report the detection of H$_2^{18}$O in this source. The dust continuum and molecular gas emission are resolved into a nearly complete ~1"5 diameter Einstein ring plus a weaker image in the center, which is caused by a special dual deflector lensing configuration. The observed line profiles of the CO, H2O and H2O+ lines are strikingly similar. In the source plane, we reconstruct the dust continuum images and the spectral cubes of the line emission at sub-kpc scales. The reconstructed dust emission in the source plane is dominated by a compact disk with an effective radius of 0.7kpc plus an overlapping extended disk with a radius twice as large. While the average magnification for the dust continuum is $\mu$~10-11, the magnification of the line emission varies 5 to 22 across different velocity components. The emission lines have similar spatial and kinematic distributions. The molecular gas and dust content reveal that G09v1.97 is a gas-rich major merger in its pre-coalescence phase. Both of the merging companions are intrinsically ULIRGs with LIR reaching $\gtrsim 4\times10^{12}L_\odot$, and the total LIR of G09v1.97 is $1.4\times10^{13}L_\odot$. The approaching southern galaxy shows no obvious kinematic structure with a semi-major half-light radius a_s=0.4kpc, while the receding galaxy resembles an a_s=1.2kpc rotating disk. The two galaxies are separated by a projected distance of 1.3kpc, bridged by weak line emission that is co-spatially located with the cold-dust-emission peak, suggesting a large amount of cold ISM in the interacting region. (abridged)

This work presents AutoLens, the first entirely automated modeling suite for the analysis of galaxy-scale strong gravitational lenses. AutoLens simultaneously models the lens galaxy's light and mass whilst reconstructing the extended source galaxy on an adaptive pixel-grid. The method's approach to source-plane discretization is amorphous, adapting its clustering and regularization to the intrinsic properties of the lensed source. The lens's light is fitted using a superposition of Sersic functions, allowing AutoLens to cleanly deblend its light from the source. Single component mass models representing the lens's total mass density profile are demonstrated, which in conjunction with light modeling can detect central images using a centrally cored profile. Decomposed mass modeling is also shown, which can fully decouple a lens's light and dark matter and determine whether the two component are geometrically aligned. The complexity of the light and mass models are automatically chosen via Bayesian model comparison. These steps form AutoLens's automated analysis pipeline, such that all results in this work are generated without any user-intervention. This is rigorously tested on a large suite of simulated images, assessing its performance on a broad range of lens profiles, source morphologies and lensing geometries. The method's performance is excellent, with accurate light, mass and source profiles inferred for data sets representative of both existing Hubble imaging and future Euclid wide-field observations.

We present integral field spectroscopy of galaxy cluster Abell 3827, using ALMA and VLT/MUSE. It reveals an unusual configuration of strong gravitational lensing in the cluster core, with at least seven lensed images of a single background spiral galaxy. Lens modelling based on HST imaging had suggested that the dark matter associated with one of the cluster's central galaxies may be offset. The new spectroscopic data enable better subtraction of foreground light, and better identification of multiple background images. The inferred distribution of dark matter is consistent with being centered on the galaxies, as expected by LCDM. Each galaxy's dark matter also appears to be symmetric. Whilst we do not find an offset between mass and light (suggestive of self-interacting dark matter) as previously reported, the numerical simulations that have been performed to calibrate Abell 3827 indicate that offsets and asymmetry are still worth looking for in collisions with particular geometries. Meanwhile, ALMA proves exceptionally useful for strong lens image identifications.

We present a new pixelized method for the inversion of gravitationally lensed extended source images which we term adaptive semi-linear inversion (SLI). At the heart of the method is an h-means clustering algorithm which is used to derive a source plane pixelization that adapts to the lens model magnification. The distinguishing feature of adaptive SLI is that every pixelization is derived from a random initialization, ensuring that data discretization is performed in a completely different and unique way for every lens model parameter set. We compare standard SLI on a fixed source pixel grid with the new method and demonstrate the shortcomings of the former when modeling singular power law ellipsoid (SPLE) lens profiles. In particular, we demonstrate the superior reliability and efficiency of adaptive SLI which, by design, fixes the number of degrees of freedom (NDOF) of the optimization and thereby removes biases present with other methods that allow the NDOF to vary. In addition, we highlight the importance of data discretization in pixel-based inversion methods, showing that adaptive SLI averages over significant systematics that are present when a fixed source pixel grid is used. In the case of the SPLE lens profile, we show how the method successfully samples its highly degenerate posterior probability distribution function with a single non-linear search. The robustness of adaptive SLI provides a firm foundation for the development of a strong lens modeling pipeline, which will become necessary in the short-term future to cope with the increasing rate of discovery of new strong lens systems.

We have modelled Atacama Large Millimeter/sub-millimeter Array (ALMA) long baseline imaging of the strong gravitational lens system H-ATLAS J090311.6+003906 (SDP.81). We have reconstructed the distribution of band 6 and 7 continuum emission in the z=3.042 source and we have determined its kinematic properties by reconstructing CO(5-4) and CO(8-7) line emission in bands 4 and 6. The continuum imaging reveals a highly non-uniform distribution of dust with clumps on scales of 200 pc. In contrast, the CO line emission shows a relatively smooth, disk-like velocity field which is well fit by a rotating disk model with an inclination angle of 40+/-5 degrees and an asymptotic rotation velocity of 320 km/s. The inferred dynamical mass within 1.5kpc is (3.5+/-0.5)x10^{10} M_sol which is comparable to the total molecular gas masses of (2.7+/-0.5)x10^{10} M_sol and (3.5+/-0.6)x10^{10} M_sol from the dust continuum emission and CO emission respectively. Our new reconstruction of the lensed HST near-infrared emission shows two objects which appear to be interacting, with the rotating disk of gas and dust revealed by ALMA distinctly offset from the near-infrared emission. The clumpy nature of the dust and a low value of the Toomre parameter of Q=0.3 suggest that the disk is in a state of collapse. We estimate a star formation rate in the disk of 470+/-80 M_sol/yr with an efficiency 65 times greater than typical low-redshift galaxies. Our findings add to the growing body of evidence that the most infra-red luminous, dust obscured galaxies in the high redshift Universe represent a population of merger induced starbursts.

We present Keck-Adaptive Optics and Hubble Space Telescope high resolution near-infrared (IR) imaging for 500 um-bright candidate lensing systems identified by the Herschel Multi-tiered Extra-galactic Survey (HerMES) and Herschel Astrophysical Terahertz Survey (H-ATLAS). Out of 87 candidates with near-IR imaging, 15 (~17%) display clear near-IR lensing morphologies. We present near-IR lens models to reconstruct and recover basic rest-frame optical morphological properties of the background galaxies from 12 new systems. Sources with the largest near-IR magnification factors also tend to be the most compact, consistent with the size bias predicted from simulations and pre- vious lensing models for sub-millimeter galaxies. For four new sources that also have high-resolution sub-mm maps, we test for differential lensing between the stellar and dust components and find that the 880 um magnification factor (u_880) is ~1.5 times higher than the near-IR magnification factor (u_NIR), on average. We also find that the stellar emission is ~2 times more extended in size than dust. The rest-frame optical properties of our sample of Herschel-selected lensed SMGs are consistent with those of unlensed SMGs, which suggests that the two populations are similar.

We have determined the mass-density radial profiles of the first five strong gravitational lens systems discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We present an enhancement of the semi-linear lens inversion method of Warren & Dye which allows simultaneous reconstruction of several different wavebands and apply this to dual-band imaging of the lenses acquired with the Hubble Space Telescope. The five systems analysed here have lens redshifts which span a range, 0.22<z<0.94. Our findings are consistent with other studies by concluding that: 1) the logarithmic slope of the total mass density profile steepens with decreasing redshift; 2) the slope is positively correlated with the average total projected mass density of the lens contained within half the effective radius and negatively correlated with the effective radius; 3) the fraction of dark matter contained within half the effective radius increases with increasing effective radius and increases with redshift.