Einstein's theory of gravity, General Relativity, has been precisely tested on Solar System scales, but the long-range nature of gravity is still poorly constrained. The nearby strong gravitational lens, ESO 325-G004, provides a laboratory to probe the weak-field regime of gravity and measure the spatial curvature generated per unit mass, $\gamma$. By reconstructing the observed light profile of the lensed arcs and the observed spatially resolved stellar kinematics with a single self-consistent model, we conclude that $\gamma = 0.97 \pm 0.09$ at 68% confidence. Our result is consistent with the prediction of 1 from General Relativity and provides a strong extragalactic constraint on the weak-field metric of gravity.

We present models for the dark and luminous mass structure of 12 strong lensing early-type galaxies (ETGs). We combine pixel-based modelling of multiband HST/ACS imaging with Jeans modelling of kinematics obtained from Keck/ESI spectra to disentangle the dark and luminous contributions to the mass. Assuming a gNFW profile for the dark matter halo and a spatially constant stellar-mass-to-light ratio $\Upsilon_{\star}$ for the baryonic mass, we infer distributions for $\Upsilon_{\star}$ consistent with IMFs that are heavier than the Milky Way's (with a global mean mismatch parameter relative to a Chabrier IMF $\mu_{\alpha c} = 1.80 \pm 0.14$) and halo inner density slopes which span a large range but are generally cuspier than the dark-matter-only prediction ($\mu_{\gamma'} = 2.01_{-0.22}^{+0.19}$). We investigate possible reasons for overestimating the halo slope, including the neglect of spatially varying stellar-mas-to-light ratios and/or stellar orbital anisotropy, and find that a quarter of the systems prefer radially declining stellar-mass-to-light ratio gradients, but that the overall effect on our inference on the halo slope is small. We suggest a coherent explanation of these results in the context of inside-out galaxy growth, and that the relative importance of different baryonic processes in shaping the dark halo may depend on halo environment.

We present the first constraints on stellar mass-to-light ratio gradients in an early-type galaxy (ETG) using multiple dynamical tracer populations to model the dark and luminous mass structure simultaneously. We combine the kinematics of the central starlight, two globular cluster populations and satellite galaxies in a Jeans analysis to obtain new constraints on M87's mass structure, employing a flexible mass model which allows for radial gradients in the stellar mass-to-light ratio. We find that, in the context of our model, a radially declining stellar-mass-to-light ratio is strongly favoured. Modelling the stellar mass-to-light ratio as following a power law, $\Upsilon_{\star} \sim R^{-\mu}$, we infer a power-law slope $\mu = -0.54 \pm 0.05$; equally, parameterising the stellar-mass-to-light ratio via a central mismatch parameter relative to a Salpeter IMF, $\alpha$, and scale radius $R_M$, we find $\alpha > 1.48$ at $95\%$ confidence and $R_M = 0.35 \pm 0.04$ kpc. We use stellar population modelling of high-resolution 11-band HST photometry to show that such a steep gradient cannot be achieved by variations in only the metallicity, age, dust extinction and star formation history if the stellar initial mass function (IMF) remains spatially constant. On the other hand, the stellar mass-to-light ratio gradient that we find is consistent with an IMF whose inner slope changes such that it is Salpeter-like in the central $\sim 0.5$ kpc and becomes Chabrier-like within the stellar effective radius. This adds to recent evidence that the non-universality of the IMF in ETGs may be confined to their core regions, and points towards a picture in which the stars in these central regions may have formed in fundamentally different physical conditions.

We present the discovery of 3 quasar lenses in the Sloan Digital Sky Survey (SDSS), selected using two novel photometry-based selection techniques. The J0941+0518 system, with two point sources separated by 5.46" on either side of a galaxy, has source and lens redshifts $z_s = 1.54$ and $z_l = 0.343$. The AO-assisted images of J2211+1929 show two point sources separated by 1.04", corresponding to the same quasar at $z_s = 1.07,$ besides the lens galaxy and Einstein ring. Images of J2257+2349 show two point sources separated by 1.67" on either side of an E/S0 galaxy. The extracted spectra show two images of the same quasar at redshift $z_s = 2.10$. In total, the two selection techniques identified 309 lens candidates, including 47 known lenses, and 6 previously ruled out candidates. 55 of the remaining candidates were observed using NIRC2 and ESI at Keck Observatory, EFOSC2 at the ESO-NTT (La Silla), and SAM and the Goodman spectrograph at SOAR. Of the candidates observed, 3 were confirmed as lenses, 36 were ruled out, and 16 remain inconclusive. Taking into account that we recovered known lenses, this gives us a success rate of at least 50/309 (16%). This initial campaign demonstrates the power of purely photometric selection techniques in finding lensed quasars. Developing and refining these techniques is essential for efficient identification of these rare lenses in ongoing and future photometric surveys.

Flux ratio anomalies in quasar lenses can be attributed to dark matter substructure surrounding the lensing galaxy and, thus, used to constrain the substructure mass fraction. Previous applications of this approach infer a substructure abundance that potentially in tension with the predictions of a $\Lambda$CDM cosmology. However, the assumption that all flux ratio anomalies are due to substructure is a strong one, and alternative explanations have not been fully investigated. Here, we use new high-resolution near-IR Keck~II adaptive optics imaging for the lens system CLASS B0712+472 to perform pixel-based lens modelling for this system and, in combination with new VLBA radio observations, show that the inclusion of the disc in the lens model can explain the flux ratio anomalies without the need for dark matter substructures. The projected disc mass comprises 16% of the total lensing mass within the Einstein radius and the total disc mass is $1.79 \times 10^{10} M_{sun}$. The case of B0712+472 adds to the evidence that not all flux ratio anomalies are due to dark subhaloes, and highlights the importance of taking the effects of baryonic structures more fully into account in order to obtain an accurate measure of the substructure mass fraction.

We present an exploration of the mass structure of a sample of 12 strongly lensed massive, compact early-type galaxies at redshifts $z\sim0.6$ to provide further possible evidence for their inside-out growth. We obtain new ESI/Keck spectroscopy and infer the kinematics of both lens and source galaxies, and combine these with existing photometry to construct (a) the fundamental plane (FP) of the source galaxies and (b) physical models for their dark and luminous mass structure. We find their FP to be tilted towards the virial plane relative to the local FP, and attribute this to their unusual compactness, which causes their kinematics to be totally dominated by the stellar mass as opposed to their dark matter; that their FP is nevertheless still inconsistent with the virial plane implies that both the stellar and dark structure of early-type galaxies is non-homologous. We also find the intrinsic scatter of their FP to be comparable to the local value, indicating that variations in the stellar mass structure outweight variations in the dark halo in the central regions of early-type galaxies. Finally, we show that inference on the dark halo structure -- and, in turn, the underlying physics -- is sensitive to assumptions about the stellar initial mass function (IMF), but that physically-motivated assumptions about the IMF imply haloes with sub-NFW inner density slopes, and may present further evidence for the inside-out growth of compact early-type galaxies via minor mergers and accretion.

We present a new sample of strong gravitational lens systems where both the foreground lenses and background sources are early-type galaxies. Using imaging from HST/ACS and Keck/NIRC2, we model the surface brightness distributions and show that the sources form a distinct population of massive, compact galaxies at redshifts $0.4 \lesssim z \lesssim 0.7$, lying systematically below the size-mass relation of the global elliptical galaxy population at those redshifts. These may therefore represent relics of high-redshift red nuggets or their partly-evolved descendants. We exploit the magnifying effect of lensing to investigate the structural properties, stellar masses and stellar populations of these objects with a view to understanding their evolution. We model these objects parametrically and find that they generally require two S\'ersic components to properly describe their light profiles, with one more spheroidal component alongside a more envelope-like component, which is slightly more extended though still compact. This is consistent with the hypothesis of the inside-out growth of these objects via minor mergers. We also find that the sources can be characterised by red-to-blue colour gradients as a function of radius which are stronger at low redshift -- indicative of ongoing accretion -- but that their environments generally appear consistent with that of the general elliptical galaxy population, contrary to recent suggestions that these objects are predominantly associated with clusters.

By constructing scaling relations for galaxies in the massive cluster MACSJ0717.5 at $z=0.545$ and comparing with those of Coma, we model the luminosity evolution of the stellar populations and the structural evolution of the galaxies. We calculate magnitudes, surface brightnesses and effective radii using HST/ACS images and velocity dispersions using Gemini/GMOS spectra, and present a catalogue of our measurements for 17 galaxies. We also generate photometric catalogues for $\sim 3000$ galaxies from the HST imaging. With these, we construct the colour-magnitude relation, the fundamental plane, the mass-to-light versus mass relation, the mass-size relation and the mass-velocity dispersion relation for both clusters. We present a new, coherent way of modelling these scaling relations simultaneously using a simple physical model in order to infer the evolution in luminosity, size and velocity dispersion as a function of redshift, and show that the data can be fully accounted for with this model. We find that (a) the evolution in size and velocity dispersion undergone by these galaxies between $z \sim 0.5$ and $z \sim 0$ is mild, with $R_e(z) \sim (1+z)^{-0.40\pm0.32}$ and $\sigma(z) \sim (1+z)^{0.09 \pm 0.27}$, and (b) the stellar populations are old, $\sim 10$ Gyr, with a $\sim 3$ Gyr dispersion in age, and are consistent with evolving purely passively since $z \sim 0.5$ with $\Delta \log M/L_B = -0.55_{-0.07}^{+0.15} z$. The implication is that these galaxies formed their stars early and subsequently grew dissipationlessly so as to have their mass already in place by $z \sim 0.5$, and suggests a dominant role for dry mergers, which may have accelerated the growth in these high-density cluster environments.

We present a general method to identify infalling substructure in discrete datasets with position and line-of-sight velocity data. We exploit the fact that galaxies falling onto a brightest cluster galaxy (BCG) in a virialised cluster, or dwarf satellites falling onto a central galaxy like the Milky Way, follow nearly radial orbits. If the orbits are exactly radial, we show how to find the probability distribution for a satellite's energy, given a tracer density for the satellite population, by solving an Abel integral equation. This is an extension of Eddington (1916)'s classical formula for the isotropic distribution function. When applied to a system of galaxies, clustering in energy space can then be quantified using the Kullback-Leibler divergence, and groups of objects can be identified which, though separated in the sky, may be falling in on the same orbit. This method is tested using mock data and applied to the satellite galaxy population around M87, the BCG in Virgo, and a number of associations are found which may represent infalling galaxy groups.

Following a number of conflicting studies of M87's mass profile, we undertake a dynamical analysis of multiple tracer populations to constrain its mass over a large radius range. We combine stellar kinematics in the central regions with the dynamics of 612 globular clusters out to 200 kpc and satellite galaxies extending to scales comparable with the virial radius. Using a spherical Jeans analysis, we are able to disentangle the mass contributions from the dark and baryonic components and set constraints on the structure of each. Assuming isotropy, we explore four different models for the dark matter halo and find that a centrally-cored dark matter distribution is preferred. We infer a stellar mass-to-light ratio $\Upsilon_{\star,v} = 6.9 \pm 0.1$ -- consistent with a Salpeter-like IMF -- and a core radius $r_c = 67 \pm 20$ kpc. We then introduce anisotropy and find that, while the halo remains clearly cored, the radial stellar anisotropy has a strong impact on both $\Upsilon_{\star,v}$ and the core's radius; here we find $\Upsilon_{\star,v} = 3.50_{-0.36}^{+0.32}$ -- consistent with a Chabrier-like IMF -- and $r_c = 19.00_{-8.34}^{+8.38}$ kpc. Thus the presence of a core at the centre of the dark halo is robust against anisotropy assumptions, while the stellar mass and core size are not. We are able to reconcile previously discrepant studies by showing that modelling the globular cluster data alone leads to the very different inference of a super-NFW cusp, thus highlighting the value of multiple-population modelling, and we point to the possible role of M87's AGN and the cluster environment in forming the central dark matter core.

We present a new photometric catalogue of the rich globular cluster (GC) system around M87, the brightest cluster galaxy in Virgo. Using archival Next Generation Virgo cluster Survey (NGVS) images in the ugriz bands, observed with CFHT/MegaPrime, we perform a careful subtraction of the galaxy's halo light in order to detect objects at small galactocentric radii as well as in the wider field, and find 17620 GC candidates over a radius range from 1.3 kpc to 445 kpc with g < 24 magnitudes. By inferring their colour, radial and magnitude distributions in a Bayesian way, we find that they are well described as a mixture of two GC populations and two distinct contaminant populations, but confirm earlier findings of radius-dependent colour gradients in both GC populations. This is consistent with a picture in which the more enriched GCs reside deeper in the galaxy's potential well, indicating a role for dissipative collapse in the formation of both the red and the blue GCs.