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