Search

While quiescent galaxies have comparable amounts of cool gas in their outer circumgalactic medium (CGM) compared to star-forming galaxies, they have significantly less interstellar gas. However, open questions remain on the processes causing galaxies to stop forming stars and stay quiescent . Theories suggest dynamical interactions with the hot corona prevent cool gas from reaching the galaxy, therefore predicting the inner regions of quiescent galaxy CGMs are devoid of cool gas. However, there is a lack of understanding of the inner regions of CGMs due to the lack of spatial information in quasar-sightline methods. We present integral-field spectroscopy probing 10--20~kpc (2.4--4.8 R\textsubscript{e}) around a massive quiescent galaxy using a gravitationally lensed star-forming galaxy. We detect absorption from Magnesium (MgII) implying large amounts of cool atomic gas (10\textsuperscript{8.4} -- 10\textsuperscript{9.3} M\textsubscript{$\odot$} with T$\sim$10\textsuperscript{4} Kelvin), in comparable amounts to star-forming galaxies. Lens modeling of Hubble imaging also reveals a diffuse asymmetric component of significant mass consistent with the spatial extent of the MgII absorption, and offset from the galaxy light profile. This study demonstrates the power of galaxy-scale gravitational lenses to not only probe the gas around galaxies, but to also independently probe the mass of the CGM due to it's gravitational effect., Comment: accepted nature communications physics

The surface brightness distribution of massive early-type galaxies (ETGs) often deviates from a perfectly elliptical shape. To capture these deviations in their isophotes during an ellipse fitting analysis, Fourier modes of order $m = 3, 4$ are often used. In such analyses the centre of each ellipse is treated as a free parameter, which may result in offsets from the centre of light, particularly for ellipses in the outer regions. This complexity is not currently accounted for in the mass models used in either strong gravitational lensing or galactic dynamical studies. In this work, we adopt a different approach, using the $m=1$ Fourier mode to account for this complexity while keeping the centres of all perturbed ellipses fixed, showing that it fits the data equally well. We applied our method to the distribution of light emission to a sample of ETGs from the MASSIVE survey and found that the majority have low $m_1$ amplitudes, below 2 percent. Five out of the 30 galaxies we analysed have high $m_1$ amplitudes, ranging from 2 to 10 percent in the outer parts ($R \gtrsim 3$ kpc), all of which have a physically associated companion. Based on our findings, we advocate the use of the $m=1$ multipole in the mass models used in strong lensing and dynamical studies, particularly for galaxies with recent or ongoing interactions., Comment: 7 pages, 6 figures. Submitted to MNRAS

Due to the space radiation environment at L2, ESA's Euclid mission will be subject to a large amount of highly energetic particles over its lifetime. These particles can cause damage to the detectors by creating defects in the silicon lattice. These defects degrade the returned image in several ways, one example being a degradation of the Charge Transfer Efficiency, which appears as readout trails in the image data. This can be problematic for the Euclid VIS instrument, which aims to measure the shapes of galaxies to a very high degree of accuracy. Using a special clocking technique called trap pumping, the single defects in the CCDs can be detected and characterised. Being the first instrument in space with this capability, it will provide novel insights into the creation and evolution of radiation-induced defects and give input to the radiation damage correction of the scientific data. We present the status of the radiation damage of the Euclid VIS CCDs and how it has evolved over the first year in space., Comment: Submitted to the Proceedings for the SPIE Space Telescopes and Instrumentation 2024: Optical, Infrared, and Millimeter Wave, paper no. 13092-24, 10 pages, 6 figures

We examine the morphological and kinematical properties of SPT-2147, a strongly lensed, massive, dusty, star-forming galaxy at $z = 3.762$. Combining data from JWST, HST, and ALMA, we study the galaxy's stellar emission, dust continuum and gas properties. The imaging reveals a central bar structure in the stars and gas embedded within an extended disc with a spiral arm-like feature. The kinematics confirm the presence of the bar and of the regularly rotating disc. Dynamical modeling yields a dynamical mass, ${M}_{\rm dyn} = (9.7 \pm 2.0) \times 10^{10}$ ${\rm M}_{\odot}$, and a maximum rotational velocity to velocity dispersion ratio, $V / \sigma = 9.8 \pm 1.2$. From multi-band imaging we infer, via SED fitting, a stellar mass, ${M}_{\star} = (6.3 \pm 0.9) \times 10^{10}$ $\rm{M}_{\odot}$, and a star formation rate, ${\rm SFR} = 781 \pm 99$ ${\rm M_{\odot} yr^{-1}}$, after correcting for magnification. Combining these measurements with the molecular gas mass, we derive a baryonic-to-total mass ratio of ${M}_{\rm bar} / {M}_{\rm dyn} = 0.9 \pm 0.2$ within 4.0 kpc. This finding suggests that the formation of bars in galaxies begins earlier in the history of the Universe than previously thought and can also occur in galaxies with elevated gas fractions., Comment: 19 pages, 11 figures. Submitted to MNRAS

In a strong gravitational lensing system, the distorted light from a source is analysed to infer the properties of the lens. However, light emitted by the lens itself can contaminate the image of the source, introducing systematic errors in the analysis. We present a simple and efficient lens light model based on the well-tested multi-Gaussian expansion (MGE) method for representing galaxy surface brightness profiles, which we combine with a semi-linear inversion scheme for pixelized source modelling. Testing it against realistic mock lensing images, we show that our scheme can fit the lensed images to the noise level, with relative differences between the true input and best-fit lens light model remaining below 5%. We apply the MGE lens light model to 38 lenses from the HST SLACS sample. We find that the new scheme provides a good fit for the majority of the sample with only 3 exceptions -- these show clear asymmetric residuals in the lens light. We examine the radial dependence of the ellipticity and position angles and confirm that it is common for a typical lens galaxy to exhibit twisting, non-elliptical isophotes and boxy / disky isophotes. Our MGE lens light model will be a valuable tool for understanding the hidden complexity of the lens mass distribution., Comment: 22 pages, 15 figures. Submitted to MNRAS. Comments Welcome

With their high angular resolutions of 30--100 mas, large fields of view, and complex optical systems, imagers on next-generation optical/near-infrared space observatories, such as the Near-Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST), present both new opportunities for science and also new challenges for empirical point spread function (PSF) characterization. In this context, we introduce ShOpt, a new PSF fitting tool developed in Julia and designed to bridge the advanced features of PIFF (PSFs in the Full Field of View) with the computational efficiency of PSFEx (PSF Extractor). Along with ShOpt, we propose a suite of non-parametric statistics suitable for evaluating PSF fit quality in space-based imaging. Our study benchmarks ShOpt against the established PSF fitters PSFEx and PIFF using real and simulated COSMOS-Web Survey imaging. We assess their respective PSF model fidelity with our proposed diagnostic statistics and investigate their computational efficiencies, focusing on their processing speed relative to the complexity and size of the PSF models. We find that ShOpt can already achieve PSF model fidelity comparable to PSFEx and PIFF while maintaining competitive processing speeds, constructing PSF models for large NIRCam mosaics within minutes., Comment: 53 pages, 27 figures, submitted to Astronomical Journal. Same as v2 with a corrected typo, which contained revisions from v1

Supermassive black holes (SMBHs) are a key catalyst of galaxy formation and evolution, leading to an observed correlation between SMBH mass $M_{\rm BH}$ and host galaxy velocity dispersion $\sigma_{\rm e}$. Outside the local Universe, measurements of $M_{\rm BH}$ are usually only possible for SMBHs in an active state: limiting sample size and introducing selection biases. Gravitational lensing makes it possible to measure the mass of non-active SMBHs. We present models of the $z=0.169$ galaxy-scale strong lens Abell~1201. A cD galaxy in a galaxy cluster, it has sufficient `external shear' that a magnified image of a $z = 0.451$ background galaxy is projected just $\sim 1$ kpc from the galaxy centre. Using multi-band Hubble Space Telescope imaging and the lens modeling software $\texttt{PyAutoLens}$ we reconstruct the distribution of mass along this line of sight. Bayesian model comparison favours a point mass with $M_{\rm BH} = 3.27 \pm 2.12\times10^{10}\,$M$_{\rm \odot}$ (3$\sigma$ confidence limit); an ultramassive black hole. One model gives a comparable Bayesian evidence without a SMBH, however we argue this model is nonphysical given its base assumptions. This model still provides an upper limit of $M_{\rm BH} \leq 5.3 \times 10^{10}\,$M$_{\rm \odot}$, because a SMBH above this mass deforms the lensed image $\sim 1$ kpc from Abell 1201's centre. This builds on previous work using central images to place upper limits on $M_{\rm BH}$, but is the first to also place a lower limit and without a central image being observed. The success of this method suggests that surveys during the next decade could measure thousands more SMBH masses, and any redshift evolution of the $M_{\rm BH}$--$\sigma_{\rm e}$ relation. Results are available at https://github.com/Jammy2211/autolens_abell_1201., Comment: Accepted in MNRAS, 27 pages, 22 figures

The distribution of mass in galaxy-scale strong gravitational lenses is often modelled as an elliptical power law plus 'external shear', which notionally accounts for neighbouring galaxies and cosmic shear. We show that it does not. Except in a handful of rare systems, the best-fit values of external shear do not correlate with independent measurements of shear: from weak lensing in 45 Hubble Space Telescope images, or in 50 mock images of lenses with complex distributions of mass. Instead, the best-fit shear is aligned with the major or minor axis of 88% of lens galaxies; and the amplitude of the external shear increases if that galaxy is disky. We conclude that 'external shear' attached to a power law model is not physically meaningful, but a fudge to compensate for lack of model complexity. Since it biases other model parameters that are interpreted as physically meaningful in several science analyses (e.g. measuring galaxy evolution, dark matter physics or cosmological parameters), we recommend that future studies of galaxy-scale strong lensing should employ more flexible mass models., Comment: 13 pages, 11 figures, submitted to MNRAS and Review Response Submitted

The fast growth of supermassive black holes and their feedback to the host galaxies play an important role in regulating the evolution of galaxies, especially in the early Universe. However, due to cosmological dimming and the limited angular resolution of most observations, it is difficult to resolve the feedback from the active galactic nuclei (AGN) to their host galaxies. Gravitational lensing, for its magnification, provides a powerful tool to spatially differentiate emission originated from AGN and host galaxy at high redshifts. Here we report a discovery of a radio lobe in a strongly lensed starburst quasar, H1413+117 or Cloverleaf at redshift $z= 2.56$, based on observational data at optical, sub-millimetre, and radio wavelengths. With both parametric and non-parametric lens models and with reconstructed images on the source plane, we find a differentially lensed, kpc scaled, single-sided radio lobe, located at ${\sim}1.2\,\mathrm{kpc}$ to the north west of the host galaxy on the source plane. From the spectral energy distribution in radio bands, we find that the radio lobe has an energy turning point residing between 1.5 GHz and 8 GHz, indicating an age of 20--50 Myr. This could indicate a feedback switching of Cloverleaf quasar from the jet mode to the quasar mode.

We present the survey design, implementation, and outlook for COSMOS-Web, a 255 hour treasury program conducted by the James Webb Space Telescope in its first cycle of observations. COSMOS-Web is a contiguous 0.54 deg$^2$ NIRCam imaging survey in four filters (F115W, F150W, F277W, and F444W) that will reach 5$\sigma$ point source depths ranging $\sim$27.5-28.2 magnitudes. In parallel, we will obtain 0.19 deg$^2$ of MIRI imaging in one filter (F770W) reaching 5$\sigma$ point source depths of $\sim$25.3-26.0 magnitudes. COSMOS-Web will build on the rich heritage of multiwavelength observations and data products available in the COSMOS field. The design of COSMOS-Web is motivated by three primary science goals: (1) to discover thousands of galaxies in the Epoch of Reionization ($64$ and place constraints on the formation of the Universe's most massive galaxies ($M_\star>10^{10}$\,M$_\odot$), and (3) directly measure the evolution of the stellar mass to halo mass relation using weak gravitational lensing out to $z\sim2.5$ and measure its variance with galaxies' star formation histories and morphologies. In addition, we anticipate COSMOS-Web's legacy value to reach far beyond these scientific goals, touching many other areas of astrophysics, such as the identification of the first direct collapse black hole candidates, ultracool sub-dwarf stars in the Galactic halo, and possibly the identification of $z>10$ pair-instability supernovae. In this paper we provide an overview of the survey's key measurements, specifications, goals, and prospects for new discovery., Comment: 46 pages, 16 figures, ApJ accepted

The cold dark matter (DM) model predicts that every galaxy contains thousands of DM subhalos; almost all other DM models include a physical process that smooths away the subhalos. The subhalos are invisible, but could be detected via strong gravitational lensing, if they lie on the line of sight to a multiply-imaged background source, and perturb its apparent shape. We present a predominantly automated strong lens analysis framework, and scan for DM subhalos in Hubble Space Telescope imaging of 54 strong lenses. We identify five DM subhalo candidates, including two especially compelling candidates (one previously known in SLACS0946+1006) where a subhalo is favoured after all of our tests for systematics. We find that the detectability of subhalos depends upon the assumed parametric form for the lens galaxy's mass distribution, especially its degree of azimuthal freedom. Using separate components for dark matter and stellar mass reveals two DM subhalo candidates and removes four false-positives compared to the single power-law mass model that is common in the literature. We identify 45 lenses without substructures, the number of which is key to statistical tests able to rule out models of e.g. warm or self-interacting DM. Our full analysis results are available at https://github.com/Jammy2211/autolens_subhalo., Comment: 25 Pages, 15 Figures

Observations suggest that the dark matter and stars in early-type galaxies `conspire' to produce a surprisingly simple distribution of total mass, $\rho(r)\propto\rho^{-\gamma}$, with $\gamma\approx2$. We measure the distribution of mass in 48 early-type galaxies that gravitationally lens a resolved background source. By fitting the source light in every pixel of images from the Hubble Space Telescope, we find a mean $\langle\gamma\rangle=2.075_{-0.024}^{+0.023}$ with intrinsic scatter between galaxies of $\sigma_\gamma=0.172^{+0.022}_{-0.032}$ for the overall sample. This is consistent with, and has similar precision to traditional techniques that employ spectroscopic observations to supplement lensing with mass estimates from stellar dynamics. Comparing measurements of $\gamma$ for individual lenses using both techniques, we find a statistically insignificant correlation of $-0.150^{+0.223}_{-0.217}$ between the two, indicating a lack of statistical power or deviations from a power-law density in certain lenses. At fixed surface mass density, we measure a redshift dependence, $\partial\langle\gamma\rangle/\partial z=0.345^{+0.322}_{-0.296}$, that is consistent with traditional techniques for the same sample of SLACS and GALLERY lenses. Interestingly, the consistency breaks down when we measure the dependence of $\gamma$ on the surface mass density of a lens galaxy. We argue that this is tentative evidence for an inflection point in the total-mass density profile at a few times the galaxy effective radius -- breaking the conspiracy., Comment: 14 pages, 11 figures, 5 tables, submitted to MNRAS

Strong gravitational lensing offers a compelling test of the cold dark matter paradigm, as it allows for subhaloes with masses of $\sim10^{9}$ M$_\odot$ and below to be detected. We test commonly-used techniques for detecting subhaloes superposed in images of strongly lensed galaxies. For the lens we take a simulated galaxy in a $\sim10^{13}$ M$_\odot$ halo grown in a high-resolution cosmological hydrodynamical simulation, which we view from two different directions. Though the resolution is high, we note the simulated galaxy still has an artificial core which adds additional complexity to the baryon dominated region. To remove particle noise, we represent the projected galaxy mass distribution by a series of Gaussian profiles which precisely capture the features of the projected galaxy. We first model the lens mass as a (broken) power-law density profile and then search for small haloes. Of the two projections, one has a regular elliptical shape, while the other has distinct deviations from an elliptical shape. For the former, the broken power-law model gives no false positives and correctly recovers the mass of the superposed small halo, but for the latter we find false positives and the inferred halo mass is overestimated by $\sim4-5$ times. We then use a more complex model in which the lens mass is decomposed into stellar and dark matter components. In this case, we show that we can capture the simulated galaxy's complex projected structures and correctly infer the input small halo., Comment: Accepted by MNRAS. Comments welcome!

The distribution of dark and luminous matter can be mapped around galaxies that gravitationally lens background objects into arcs or Einstein rings. New surveys will soon observe hundreds of thousands of galaxy lenses, and current, labour-intensive analysis methods will not scale up to this challenge. We instead develop a fully automatic, Bayesian method which we use to fit a sample of 59 lenses imaged by the Hubble Space Telescope in uniform conditions. We set out to \textit{leave no lens behind} and focus on ways in which automated fits fail in a small handful of lenses, describing adjustments to the pipeline that allows us to infer accurate lens models. Our pipeline ultimately fits {\em all} 59 lenses in our sample, with a high success rate key because catastrophic outliers would bias large samples with small statistical errors. Machine Learning techniques might further improve the two most difficult steps: subtracting foreground lens light and initialising a first, approximate lens model. After that, increasing model complexity is straightforward. We find a mean $\sim1\%$ measurement precision on the measurement of the Einstein radius across the lens sample which {\em does not degrade with redshift} up to at least $z=0.7$ -- in stark contrast to other techniques used to study galaxy evolution, like stellar dynamics. Our \texttt{PyAutoLens} software is open source, and is also installed in the Science Data Centres of the ESA Euclid mission., Comment: 23 pages, 20 figures, 5 tables. Submitted to MNRAS

We rederive the number density of intervening line-of-sight haloes relative to lens subhaloes in galaxy-galaxy strong lensing observations, where these perturbers can generate detectable image fluctuations. Previous studies have calculated the detection limit of a line-of-sight small-mass dark halo by comparing the lensing deflection angles it would cause, to those caused by a subhalo within the lens. However, this overly simplifies the difference in observational consequences between a subhalo and a line-of-sight halo. Furthermore, it does not take into account degeneracies between an extra subhalo and the uncertain properties of the main lens. More in keeping with analyses of real-world observations, we regard a line-of-sight halo as detectable only if adding it to a smooth model generates a statistically significant improvement in the reconstructed image. We find that the number density of detectable line-of-sight perturbers has been overestimated by as much as a factor of two in the previous literature. For typical lensing geometries and configurations, very deep imaging is sensitive to twice as many line-of-sight perturbers as subhaloes, but moderate depth imaging is sensitive to only slightly more line-of-sight perturbers than subhaloes., Comment: Accepted by MNRAS. Comments welcome

A defining prediction of the cold dark matter (CDM) cosmological model is the existence of a very large population of low-mass haloes. This population is absent in models in which the dark matter particle is warm (WDM). These alternatives can, in principle, be distinguished observationally because halos along the line-of-sight can perturb galaxy-galaxy strong gravitational lenses. Furthermore, the WDM particle mass could be deduced because the cut-off in their halo mass function depends on the mass of the particle. We systematically explore the detectability of low-mass haloes in WDM models by simulating and fitting mock lensed images. Contrary to previous studies, we find that halos are harder to detect when they are either behind or in front of the lens. Furthermore, we find that the perturbing effect of haloes increases with their concentration: detectable haloes are systematically high-concentration haloes, and accounting for the scatter in the mass-concentration relation boosts the expected number of detections by as much as an order of magnitude. Haloes have lower concentration for lower particle masses and this further suppresses the number of detectable haloes beyond the reduction arising from the lower halo abundances alone. Taking these effects into account can make lensing constraints on the value of the mass function cut-off at least an order of magnitude more stringent than previously appreciated., Comment: MNRAS submitted. Comments welcome

Strong gravitational lensing, which can make a background source galaxy appears multiple times due to its light rays being deflected by the mass of one or more foreground lens galaxies, provides astronomers with a powerful tool to study dark matter, cosmology and the most distant Universe. PyAutoLens is an open-source Python 3.6+ package for strong gravitational lensing, with core features including fully automated strong lens modeling of galaxies and galaxy clusters, support for direct imaging and interferometer datasets and comprehensive tools for simulating samples of strong lenses. The API allows users to perform ray-tracing by using analytic light and mass profiles to build strong lens systems. Accompanying PyAutoLens is the autolens workspace (see https://github.com/Jammy2211/autolens_workspace), which includes example scripts, lens datasets and the HowToLens lectures in Jupyter notebook format which introduce non experts to strong lensing using PyAutoLens. Readers can try PyAutoLens right now by going to the introduction Jupyter notebook on Binder (see https://mybinder.org/v2/gh/Jammy2211/autolens_workspace/master) or checkout the readthedocs (see https://pyautolens.readthedocs.io/en/latest/) for a complete overview of PyAutoLens's features., Comment: 2 pages, 1 figure

A major trend in academia and data science is the rapid adoption of Bayesian statistics for data analysis and modeling, leading to the development of probabilistic programming languages (PPL). A PPL provides a framework that allows users to easily specify a probabilistic model and perform inference automatically. PyAutoFit is a Python-based PPL which interfaces with all aspects of the modeling (e.g., the model, data, fitting procedure, visualization, results) and therefore provides complete management of every aspect of modeling. This includes composing high-dimensionality models from individual model components, customizing the fitting procedure and performing data augmentation before a model-fit. Advanced features include database tools for analysing large suites of modeling results and exploiting domain-specific knowledge of a problem via non-linear search chaining. Accompanying PyAutoFit is the autofit workspace (see https://github.com/Jammy2211/autofit_workspace), which includes example scripts and the HowToFit lecture series which introduces non-experts to model-fitting and provides a guide on how to begin a project using PyAutoFit. Readers can try PyAutoFit right now by going to the introduction Jupyter notebook on Binder (see https://mybinder.org/v2/gh/Jammy2211/autofit_workspace/HEAD) or checkout our readthedocs(see https://pyautofit.readthedocs.io/en/latest/) for a complete overview of PyAutoFit's features., Comment: Published in the Journal of Open Source Software

A fundamental prediction of the cold dark matter (CDM) model of structure formation is the existence of a vast population of dark matter haloes extending to subsolar masses. By contrast, other dark matter models, such as a warm thermal relic (WDM), predict a cutoff in the mass function at a mass which, for popular models, lies approximately between $10^7$ and $10^{10}~{\rm M}_\odot$. We use mock observations to demonstrate the viability of a forward modelling approach to extract information about low-mass dark haloes lying along the line-of-sight to galaxy-galaxy strong lenses. This can be used to constrain the mass of a thermal relic dark matter particle, $m_\mathrm{DM}$. With 50 strong lenses at Hubble Space Telescope resolution and a maximum pixel signal-to-noise ratio of $\sim50$, the expected median 2$\sigma$ constraint for a CDM-like model (with a halo mass cutoff at $10^{7}~{\rm M}_\odot$) is $m_\mathrm{DM} > 4.10 \, \mathrm{keV}$ (50% chance of constraining $m_{\rm DM}$ to be better than 4.10 keV). If, however, the dark matter is a warm particle of $m_\mathrm{DM}=2.2 \, \mathrm{keV}$, our 'Approximate Bayesian Computation' method would result in a median estimate of $m_\mathrm{DM}$ between 1.43 and 3.21 keV. Our method can be extended to the large samples of strong lenses that will be observed by future telescopes, and could potentially rule out the standard CDM model of cosmogony. To aid future survey design, we quantify how these constraints will depend on data quality (spatial resolution and integration time) as well as on the lensing geometry (source and lens redshifts)., Comment: Accepted by MNRAS. Comments welcome

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., Comment: 19 pages, 4 figures. source code: https://github.com/Jammy2211/PyAutoLens

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., Comment: 49 pages, 21 figures, MNRAS resubmitted after review

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., Comment: 10 pages, 5 figures. Comments welcome

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., Comment: 19 pages, 16 figures, 4 tables, Revised Submission to MNRAS