Search

Observations suggest that the dark matter and stars in early-type galaxies 'conspire' to produce a surprisingly simple distribution of total mass, ρ(r) ∝ ρ−γ, with γ ≈ 2. 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 an 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 γ 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 Sloan Lens ACS and GALaxy-Lyα EmitteR sYstems (GALLERY) lenses. Interestingly, the consistency breaks down when we measure the dependence of γ 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. [ABSTRACT FROM AUTHOR]

We study the highly magnified arc SGAS J122651.3+215220 caused by a star-forming galaxy at zs = 2.93 crossing the lensing caustic cast by the galaxy cluster SDSS J1226+2152 (zl = 0.43), using Hubble Space Telescope observations. We report in the arc several asymmetric surface brightness features whose angular separations are a fraction of an arcsecond from the lensing critical curve and appear to be highly but unequally magnified image pairs of underlying compact sources, with one brightest pair having clear asymmetry consistently across four filters. One explanation of unequal magnification is microlensing by intracluster stars, which induces independent flux variations in the images of individual or groups of source stars in the lensed galaxy. For a second possibility, intracluster dark matter subhaloes invisible to telescopes effectively perturb lensing magnifications near the critical curve and give rise to persistently unequal image pairs. Our modelling suggests, at least for the most prominent identified image pair, that the microlensing hypothesis is in tension with the absence of notable asymmetry variation over a six-year baseline, while subhaloes of ∼106-$10^8\, \mathrm{ M}_\odot$ anticipated from structure formation with cold dark matter typically produce stationary and sizable asymmetries. We judge that observations at additional times and more precise lens models are necessary to stringently constrain temporal variability and robustly distinguish between the two explanations. The arc under this study is a scheduled target of a Director's Discretionary Early Release Science program of the James Webb Space Telescope, which will provide deep images and a high-resolution view with integral field spectroscopy.

We report the results of the STRong lensing Insights into the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates with ESO Faint Object Spectrograph and Camera 2 on the NTT and Echellette Spectrograph and Imager onKeck, confirming 10 newlensed quasars and 10 quasar pairs. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged (DESJ0345.2545, z = 1.68), and one is quadruply imaged (quad: DESJ0053.2012, z = 3.8). Singular isothermal ellipsoid models for the doubles, based on high-resolution imaging from SAMI on Southern Astrophysical Research Telescope or Near InfraRed Camera 2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcsec. After spectroscopic follow-up, we extract multi-epoch grizY photometry of confirmed lensed quasars and contaminant quasar+star pairs from DES data using parametric multiband modelling, and compare variability in each system's components. By measuring the reduced χ2 associated with fitting all epochs to the samemagnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems. Based on our spectroscopic follow-up, this variability information improves selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with most remaining contaminants being star-forming galaxies. Using mock lensed quasar light curves we demonstrate that selection based only on variability will over-represent the quad fraction by 10 per cent over a complete DES magnitude-limited sample, explained by the magnification bias and hence lower luminosity/more variable sources in quads.

We report lensing magnifications, extinction, and time-delay estimates for the first resolved, multiply imaged Type Ia supernova iPTF16geu, at z = 0.409, using Hubble Space Telescope (HST) observations in combination with supporting ground-based data. Multiband photometry of the resolved images provides unique information about the differential dimming due to dust in the lensing galaxy. Using HST and Keck AO reference images taken after the SN faded, we obtain a total lensing magnification for iPTF16geu of μ = 67.8-2.9+2.6, accounting for extinction in the host and lensing galaxy. As expected from the symmetry of the system, we measure very short time-delays for the three fainter images with respect to the brightest one: -0.23 ± 0.99, -1.43 ± 0.74, and 1.36 ± 1.07 d. Interestingly, we find large differences between the magnifications of the four supernova images, even after accounting for uncertainties in the extinction corrections: Δm1 = - 3.88-0.06+0.07, Δm2 = -2.99-0.08+0.09, Δm3 = - 2.19-0.05+0.14, and Δm4 = - 2.40-0.12+0.14 mag, discrepant with model predictions suggesting similar image brightnesses. A possible explanation for the large differences is gravitational lensing by substructures, micro-or millilensing, in addition to the large-scale lens causing the image separations. We find that the inferred magnification is insensitive to the assumptions about the dust properties in the host and lens galaxy.

Supermassive black holes (SMBHs) are a key catalyst of galaxy formation and evolution, leading to an observed correlation between SMBH mass M BH and host galaxy velocity dispersion σe. Outside the local Universe, measurements of M 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 ∼1 kpc from the galaxy centre. Using multiband Hubble Space Telescope imaging and the lens modelling software PYAUTOLENS , we reconstruct the distribution of mass along this line of sight. Bayesian model comparison favours a point mass with M BH = 3.27 ± 2.12 × 1010 M⊙ (3σ confidence limit); an ultramassive black hole. One model gives a comparable Bayesian evidence without an SMBH; however, we argue this model is nonphysical given its base assumptions. This model still provides an upper limit of M BH ≤ 5.3 × 1010 M⊙, because an SMBH above this mass deforms the lensed image ∼1 kpc from Abell 1201's centre. This builds on previous work using central images to place upper limits on M 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 BH−σe relation. Results are available at https://github.com/Jammy2211/autolens_abell_1201. [ABSTRACT FROM AUTHOR]

We search for gravitational wave (GW) events from LIGO-Virgo's third run that may have been affected by gravitational lensing. Gravitational lensing delays the arrival of GWs, and alters their amplitude – thus biasing the inferred progenitor masses. This would provide a physically well-understood interpretation of GW detections in the 'mass gap' between neutron stars and black holes, as gravitationally lensed binary neutron star (BNS) mergers. We selected three GW detections in LIGO-Virgo's third run for which the probability of at least one of the constituent compact objects being in the mass gap was reported as high with low latency – i.e. candidate lensed BNS mergers. Our observations of powerful strong lensing clusters located adjacent to the peak of their sky localization error maps reached a sensitivity |$\rm AB\simeq 25.5$| in the z ′ band with the GMOS instruments on the Gemini telescopes, and detected no candidate lensed optical counterparts. We combine recent kilonova light-curve models with recent predictions of the lensed BNS population and the properties of the objects that we followed up to show that realistic optical counterparts were detectable in our observations. Further detailed analysis of two of the candidates suggests that they are a plausible pair of images of the same low-mass binary black hole merger, lensed by a local galaxy or small group of galaxies. This further underlines that access to accurate mass information with low latency would improve the efficiency of candidate lensed BNS selection. [ABSTRACT FROM AUTHOR]

We introduce a machine learning method for estimating the sensitivity of strong lens observations to dark matter subhaloes in the lens. Our training data include elliptical power-law lenses, Hubble Deep Field sources, external shear, and noise and PSF for the Euclid VIS instrument. We set the concentration of the subhaloes using a v max– r max relation. We then estimate the dark matter subhalo sensitivity in 16 000 simulated strong lens observations with depth and resolution resembling Euclid VIS images. We find that with a 3σ detection threshold, 2.35 per cent of pixels inside twice the Einstein radius are sensitive to subhaloes with a mass M max ≤ 1010 M⊙, 0.03 per cent are sensitive to M max ≤ 109 M⊙, and the limit of sensitivity is found to be M max = 108.8 ± 0.2 M⊙. Using our sensitivity maps and assuming CDM, we estimate that Euclid-like lenses will yield |$1.43^{+0.14}_{-0.11}[f_\mathrm{sub}^{-1}]$| detectable subhaloes per lens in the entire sample, but this increases to |$35.6^{+0.9}_{-0.9}[f_\mathrm{sub}^{-1}]$| per lens in the most sensitive lenses. Estimates are given in units of the inverse of the substructure mass fraction |$f_\mathrm{sub}^{-1}$|⁠. Assuming f sub = 0.01, one in every 70 lenses in general should yield a detection, or one in every ∼ three lenses in the most sensitive sample. From 170 000 new strong lenses detected by Euclid, we expect ∼2500 new subhalo detections. We find that the expected number of detectable subhaloes in warm dark matter models only changes relative to cold dark matter for models which have already been ruled out, i.e. those with half-mode masses M hm > 108 M⊙. [ABSTRACT FROM AUTHOR]

Line-of-sight effects in strong gravitational lensing have long been treated as a nuisance. However, it was recently proposed that the line-of-sight shear could be a cosmological observable in its own right, if it is not degenerate with lens model parameters. We first demonstrate that the line-of-sight shear can be accurately measured from a simple simulated strong lensing image with per cent precision. We then extend our analysis to more complex simulated images and stress test the recovery of the line-of-sight shear when using deficient fitting models, finding that it escapes from degeneracies with lens model parameters, albeit at the expense of the precision. Lastly, we check the validity of the tidal approximation by simulating and fitting an image generated in the presence of many line-of-sight dark matter haloes, finding that an explicit violation of the tidal approximation does not necessarily prevent one from measuring the line-of-sight shear. [ABSTRACT FROM AUTHOR]

We report the spectroscopic follow-up of 175 lensed quasar candidates selected using Gaia Data Release 2 observations following Paper III of this series. Systems include 86 confirmed lensed quasars and a further 17 likely lensed quasars based on imaging and/or similar spectra. We also confirm 11 projected quasar pairs and 11 physical quasar pairs, while 25 systems are left as unclassified quasar pairs – pairs of quasars at the same redshift, which could be either distinct quasars or potential lensed quasars. Especially interesting objects include eight quadruply imaged quasars of which two have BAL sources, an apparent triple, and a doubly lensed LoBaL quasar. The source redshifts and image separations of these new lenses range between 0.65–3.59 and 0.78–6.23 arcsec, respectively. We compare the known population of lensed quasars to an updated mock catalogue at image separations between 1 and 4 arcsec, showing a very good match at z < 1.5. At z > 1.5, only 47 per cent of the predicted number are known, with 56 per cent of these missing lenses at image separations below 1.5 arcsec. The missing higher redshift, small-separation systems will have fainter lensing galaxies, and are partially explained by the unclassified quasar pairs and likely lenses presented in this work, which require deeper imaging. Of the 11 new reported projected quasar pairs, 5 have impact parameters below 10 kpc, almost tripling the number of such systems, which can probe the innermost regions of quasar host galaxies through absorption studies. We also report four new lensed galaxies discovered through our searches, with source redshifts ranging from 0.62 to 2.79. [ABSTRACT FROM AUTHOR]

Cross-referencing a watch list of galaxy groups and clusters with transient detections from real-time streams of wide-field survey data is a promising method for discovering gravitationally lensed explosive transients, including supernovae, kilonovae, gravitational waves, and gamma-ray bursts in the next 10 yr. However, currently there exists no catalogue of objects with both sufficient angular extent and depth to adequately perform such a search. In this study, we develop a cluster-finding method capable of creating an all-sky list of galaxy group- and cluster-scale objects out to z ≃ 1 based on their lens-plane properties and using only existing data from wide-field infrared surveys such as VISTA Hemisphere Survey and UKIRT Hemisphere Survey, and all-sky WISE (Wide-field Infrared Survey Explorer) data. In testing this method, we recover 91 per cent of a sample containing known and candidate lensing objects with Einstein radii of θE ≥ 5 arcsec. We also search the surrounding regions of this test sample for other groups and clusters using our method and verify the existence of any significant findings by visual inspection, deriving estimates of the false positive rate that are as low as 6 per cent. The method is also tested on simulated Rubin data from their Data Preview 0 programme, which yields complementary results of a good recovery rate of ≳80 per cent for M 200 ≥ 7 × 1013 M⊙ clusters and with no false positives produced in our test region. Importantly, our method is positioned to create a watch list in advance of Rubin's LSST, as it utilizes only existing data, therefore enabling the discovery of lensed transients early within the survey's lifetime. [ABSTRACT FROM AUTHOR]

General Relativity (GR) has been successfully tested mainly at Solar system scales; however, galaxy-scale tests have become popular in the last few decades. In this work, we investigate the ηPPN parameter, which is commonly defined by the ratio of two scalar potentials that appears in the cosmological linearly perturbed metric. Under the assumption of GR and a vanish anisotropic stress tensor, ηPPN = 1. Using ALMA, HST , and VLT/MUSE data, we combine mass measurements, using gravitational lensing and galactic dynamics, for the SDP.81 lens galaxy (z = 0.299) to constrain ηPPN. By using a flexible and self-consistent mass profile, our fiducial model takes into account the contribution of the stellar mass and a dark matter halo to reconstruct the lensed galaxy and the spatially resolved stellar kinematics. We infer, after accounting for systematic uncertainties related to the mass model, cosmology, and kinematics, |$\eta _{\text{PPN}} = 1.13^{+0.03}_{-0.03}\pm 0.20\, (\text{sys})$|⁠ , which is in accordance with GR predictions. Better spectroscopy data are needed to push the systematics down and bring the uncertainty to the percentage level since our analysis shows that the main source of the systematics is related to kinematics, which heavily depends on the signal-to-noise ratio of the spectra. [ABSTRACT FROM AUTHOR]

We present the clump populations detected in 18 lensed galaxies at redshifts 1--8.5 within the lensing cluster field SMACS0723. The recent JWST Early Release Observations of this poorly known region of the sky have revealed numerous point-like sources within and surrounding their host galaxies, undetected in the shallower Hubble Space Telescope images. We use JWST multi-and photometry and the lensing model of this galaxy cluster to estimate the intrinsic sizes and magnitudes of the stellar clumps. We derive optical restframe effective radii from <10 to hundreds pc and masses ranging from ∼105 to 109 M ⊙, overlapping with massive star clusters in the local universe. Clump ages range from 1 Myr to 1 Gyr. We compare the crossing time to the age of the clumps and determine that between 45 and 60  per cent of the detected clumps are consistent with being gravitationally bound. On average, the dearth of Gyr old clumps suggests that the dissolution time scales are shorter than 1 Gyr. We see a significant increase in the luminosity (mass) surface density of the clumps with redshift. Clumps in reionization era galaxies have stellar densities higher than star clusters in the local universe. We zoom in into single galaxies at redshift <6 and find for two galaxies, the Sparkler and the Firework, that their star clusters/clumps show distinctive colour distributions and location surrounding their host galaxy that are compatible with being accredited or formed during merger events. The ages of some of the compact clusters are between 1 and 4 Gyr, e.g. globular cluster precursors formed around 9--12 Gyr ago. Our study, conducted on a small sample of galaxies, shows the potential of JWST observations for understanding the conditions under which star clusters form in rapidly evolving galaxies. [ABSTRACT FROM AUTHOR]

The first James Webb Space Telescope (JWST) data on the massive colliding cluster El Gordo allow for 23 known families of multiply lensed images to be confirmed and for eight new members of these families to be identified. Based on these families, which have been confirmed spectroscopically by MUSE, we derived an initial lens model. This model guided the identification of 37 additional families of multiply lensed galaxies, among which 28 are entirely new systems, and nine were previously known. The initial lens model determined geometric redshifts for the 37 new systems. The geometric redshifts agree reasonably well with spectroscopic or photometric redshifts when those are available. The geometric redshifts enable two additional models that include all 60 families of multiply lensed galaxies spanning a redshift range 2 < z < 6. The derived dark-matter distribution confirms the double-peak configuration of mass found by earlier work with the southern and northern clumps having similar masses. We confirm that El Gordo is the most massive known cluster at z > 0:8 and has an estimated virial mass close the maximum mass allowed by standard cosmological models. The JWST images also reveal the presence of small-mass perturbers that produce small lensing distortions. The smallest of these is consistent with being a dwarf galaxy at z = 0:87 and has an estimated mass of 3:8 109 M, making it the smallest substructure found at z > 0:5. The JWST images also show several candidate caustic-crossing events. One of them is detected at high significance at the expected position of the critical curve and is likely a red supergiant star at z = 2:1878. This would be the first red supergiant found at cosmological distances. The cluster lensing should magnify background objects at z > 6, making more of them visible than in blank fields of a similar size, but there appears to be a deficiency of such objects. [ABSTRACT FROM AUTHOR]

We have measured the redshifts and single-aperture velocity dispersions of eight lens galaxies using the data collected by the Echellette Spectrograph and Imager (ESI) and Low Resolution Imaging Spectrometer (LRIS) at W.M. Keck observatory on different observing nights spread over three years (2018-2020). These results, combined with other ancillary data, such as high-resolution images of the lens systems, and time delays, are necessary to increase the sample size of the quasar-galaxy lens systems for which the Hubble constant can be measured, using the time-delay strong lensing method, hence increasing the precision of its inference. Typically, the 2D spectra of the quasar-galaxy lens systems get spatially blended due to seeing by ground-based observations. As a result, the extracted lensing galaxy (deflector) spectra become significantly contaminated by quasar light, which affects the ability to extract meaningful information about the deflector. To account for spatial blending and extract less contaminated and higher signal-to-noise ratio (S/N) 1D spectra of the deflectors, a forward modeling method has been implemented. From the extracted spectra, we have measured redshifts using prominent absorption lines and single aperture velocity dispersions using the penalized pixel fitting code ppxf. In this paper, we report the redshifts and single aperture velocity dispersions of eight lens galaxies - J0147+4630, B0445+123, B0631+519, J0659+1629, J08182613, J0924+0219, J1433+6007, and J1817+2729. Among these systems, six do not have previously measured velocity dispersions; for the other two, our measurements are consistent with previously reported values. Additionally, we have measured the previously unknown redshifts of the deflectors in J08182613 and J1817+2729 to be 0:866 ± 0:002 and 0:408 ± 0:002, respectively. [ABSTRACT FROM AUTHOR]

When strong gravitational lenses are to be used as an astrophysical or cosmological probe, models of their mass distributions are often needed. We present a new, time-efficient automation code for the uniform modeling of strongly lensed quasars with GLEE, a lens-modeling software for multiband data. By using the observed positions of the lensed quasars and the spatially extended surface brightness distribution of the host galaxy of the lensed quasar, we obtain a model of the mass distribution of the lens galaxy. We applied this uniform modeling pipeline to a sample of nine strongly lensed quasars for which images were obtained with the Wide Field Camera 3 of the Hubble Space Telescope. The models show well-reconstructed light components and a good alignment between mass and light centroids in most cases. We find that the automated modeling code significantly reduces the input time during the modeling process for the user. The time for preparing the required input files is reduced by a factor of 3 from ~3 h to about one hour. The active input time during the modeling process for the user is reduced by a factor of 10 from ~10 h to about one hour per lens system. This automated uniform modeling pipeline can efficiently produce uniform models of extensive lens-system samples that can be used for further cosmological analysis. A blind test that compared our results with those of an independent automated modeling pipeline based on the modeling software Lenstronomy revealed important lessons. Quantities such as Einstein radius, astrometry, mass flattening, and position angle are generally robustly determined. Other quantities, such as the radial slope of the mass density profile and predicted time delays, depend crucially on the quality of the data and on the accuracy with which the point spread function is reconstructed. Better data and/or a more detailed analysis are necessary to elevate our automated models to cosmography grade. Nevertheless, our pipeline enables the quick selection of lenses for follow-up and further modeling, which significantly speeds up the construction of cosmography-grade models. This important step forward will help us to take advantage of the increase in the number of lenses that is expected in the coming decade, which is an increase of several orders of magnitude. [ABSTRACT FROM AUTHOR]

We study the total and baryonic mass distributions of the deflector SDSS J0100+1818 through a full strong lensing analysis. The system is composed of an ultra-massive early-type galaxy at z = 0.581, with a total stellar mass of (1.5 ± 0.3)×1012 M⊙ and a stellar velocity dispersion of (450 ± 40) km s−1, surrounded by ten multiple images of three background sources, two of which are spectroscopically confirmed at z = 1.880. We took advantage of high-resolution HST photometry and VLT/X-shooter spectroscopy to measure the positions of the multiple images and performed a strong lensing study with the software GLEE. We tested different total mass profiles for the lens and modeled the background sources first as point-like and then as extended objects. We successfully predict the positions of the observed multiple images and reconstruct over approximately 7200 HST pixels the complex surface brightness distributions of the sources. We measured the cumulative total mass profile of the lens and find a total mass value of (9.1 ± 0.1)×1012 M⊙, within the Einstein radius of approximately 42 kpc, and stellar-over-total mass fractions ranging from (49 ± 12)%, at the half-light radius (Re = 9.3 kpc) of the lens galaxy, to (10 ± 2)%, in the outer regions (R = 70 kpc). These results suggest that the baryonic mass component of SDSS J0100+1818 is very concentrated in its core and that the lens early-type galaxy (or group) is immersed in a massive dark matter halo, which allows it to act as a powerful gravitational lens, creating multiple images with exceptional angular separations. This is consistent with what has been found in other ultra-high-mass candidates at intermediate redshift. We also measured the physical sizes of the distant sources, resolving them down to a few hundred parsecs. Finally, we quantify and discuss a relevant source of systematic uncertainties on the reconstructed sizes of background galaxies, associated with the adopted lens total mass model. [ABSTRACT FROM AUTHOR]

Modeling of strong gravitational lenses is a necessity for further applications in astrophysics and cosmology. With the large number of detections in current and upcoming surveys, such as the Rubin Legacy Survey of Space and Time (LSST), it is pertinent to investigate automated and fast analysis techniques beyond the traditional and time-consuming Markov chain Monte Carlo sampling methods. Building upon our (simple) convolutional neural network (CNN), we present here another CNN, specifically a residual neural network (ResNet), that predicts the five mass parameters of a singular isothermal ellipsoid (SIE) profile (lens center x and y, ellipticity ex and ey, Einstein radius θE) and the external shear (γext, 1, γext, 2) from ground-based imaging data. In contrast to our previous CNN, this ResNet further predicts the 1σ uncertainty for each parameter. To train our network, we use our improved pipeline to simulate lens images using real images of galaxies from the Hyper Suprime-Cam Survey (HSC) and from the Hubble Ultra Deep Field as lens galaxies and background sources, respectively. We find very good recoveries overall for the SIE parameters, especially for the lens center in comparison to our previous CNN, while significant differences remain in predicting the external shear. From our multiple tests, it appears that most likely the low ground-based image resolution is the limiting factor in predicting the external shear. Given the run time of milli-seconds per system, our network is perfectly suited to quickly predict the next appearing image and time delays of lensed transients. Therefore, we use the network-predicted mass model to estimate these quantities and compare to those values obtained from our simulations. Unfortunately, the achieved precision allows only a first-order estimate of time delays on real lens systems and requires further refinement through follow-up modeling. Nonetheless, our ResNet is able to predict the SIE and shear parameter values in fractions of a second on a single CPU, meaning that we are able to efficiently process the huge amount of galaxy-scale lenses expected in the near future. [ABSTRACT FROM AUTHOR]

Secure confirmation that a gravitational wave (GW) has been gravitationally lensed would bring together these two pillars of General Relativity for the first time. This breakthrough is challenging for many reasons, including: GW sky localization uncertainties dwarf the angular scale of gravitational lensing, the mass and structure of gravitational lenses is diverse, the mass function of stellar remnant compact objects is not yet well constrained, and GW detectors do not operate continuously. We introduce a new approach that is agnostic to the mass and structure of the lenses, compare the efficiency of different methods for lensed GW discovery, and explore detection of lensed kilonova counterparts as a direct method for localizing candidates. Our main conclusions are: (1) lensed neutron star mergers (NS–NS) are magnified into the 'mass gap' between NS and black holes, therefore selecting candidates from public GW alerts with high mass gap probability is efficient, (2) the rate of detectable lensed NS–NS will approach one per year in the mid-2020s, (3) the arrival time difference between lensed NS–NS images is |$1\, \rm s\lesssim \Delta \mathit{ t}\lesssim 1\, yr$|⁠ , and thus well-matched to the operations of GW detectors and optical telescopes, (4) lensed kilonova counterparts are faint at peak (e.g. r AB ≃ 24–26 in the mid-2020s), fade quickly (⁠|$d\lt 2\, \rm d$|⁠), and are detectable with target of opportunity observations with large wide-field telescopes. For example, just ≲ 0.25 per cent of Vera C. Rubin Observatory's observing time will be sufficient to follow up one well-localized candidate per year. Our predictions also provide a physically well-defined basis for exploring electromagnetically the exciting new 'mass gap' discovery space. [ABSTRACT FROM AUTHOR]

Despite the success of galaxy-scale strong gravitational lens studies with Hubble-quality imaging, a number of well-studied strong lenses remains small. As a result, robust comparisons of the lens models to theoretical predictions are difficult. This motivates our application of automated Bayesian lens modelling methods to observations from public data releases of overlapping large ground-based imaging and spectroscopic surveys: Kilo-Degree Survey (KiDS) and Galaxy and Mass Assembly (GAMA), respectively. We use the open-source lens modelling software pyautolens to perform our analysis. We demonstrate the feasibility of strong lens modelling with large-survey data at lower resolution as a complementary avenue to studies that utilize more time-consuming and expensive observations of individual lenses at higher resolution. We discuss advantages and challenges, with special consideration given to determining background source redshifts from single-aperture spectra and to disentangling foreground lens and background source light. High uncertainties in the best-fitting parameters for the models due to the limits of optical resolution in ground-based observatories and the small sample size can be improved with future study. We give broadly applicable recommendations for future efforts, and with proper application, this approach could yield measurements in the quantities needed for robust statistical inference. [ABSTRACT FROM AUTHOR]

Current observations favour that the massive ultraviolet-bright clumps with a median stellar mass of |$\sim 10^7\, {\rm M}_{\odot }$|⁠ , ubiquitously observed in z ∼ 1–3 galaxies, are star-forming regions formed in situ in galaxies. It has been proposed that they result from gas fragmentation due to gravitational instability of gas-rich, turbulent, and high-redshift discs. We bring support to this scenario by reporting the new discovery of giant molecular clouds (GMCs) in the strongly lensed, clumpy, main-sequence galaxy, A521-sys1, at z = 1.043. Its CO(4–3) emission was mapped with the Atacama Large Millimetre/submillimetre Array (ALMA) at an angular resolution of 0.19 × 0.16 arcsec2, reading down to 30 pc, thanks to gravitational lensing. We identified 14 GMCs, most being virialized, with |$10^{5.9}-10^{7.9}\, {\rm M}_{\odot }$| masses and a median |$800\, {\rm M}_{\odot }~\mathrm{pc}^{-2}$| molecular gas mass surface density, that are, respectively, 100 and 10 times higher than for nearby GMCs. They are also characterized by 10 times higher supersonic turbulence with a median Mach number of 60. They end up to fall above the Larson scaling relations, similarly to the GMCs in another clumpy z ≃ 1 galaxy, the Cosmic Snake, although differences between the two sets of high-redshift GMCs exist. Altogether they support that GMCs form with properties that adjust to the ambient interstellar medium conditions prevalent in the host galaxy whatever its redshift. The detected A521-sys1 GMCs are massive enough to be the parent gas clouds of stellar clumps, with a relatively high star formation efficiency per free-fall time of ∼11 per cent. [ABSTRACT FROM AUTHOR]

We report the first 5σ detection of H  i 21 cm emission from a star-forming galaxy at redshift z ∼ 1.3 (nearly 9 billion years ago) using upgraded Giant Metrewave Radio Telescope (uGMRT). This is the highest redshift H  i detection in emission from an individual galaxy to date. The emission is strongly boosted by the gravitational lens, an early-type elliptical galaxy, at redshift z ∼ 0.13. The measured H  i mass of the galaxy is |$M_{\rm H\, \small {\rm I}} = (0.90 \pm 0.14 \pm 0.05) \times 10^{10}\, \mathrm{M}_{\odot }$|⁠ , which is almost twice the inferred stellar mass of the galaxy, indicating an extended structure of the H  i gas inside the galaxy. By fitting 2D Gaussian to the H  i signal at the peak of the spectral line, we find the source to be marginally resolved with the position angle consistent with the emission being tangential to the critical curve of the lens mass distribution. This indicates that the solid angle of the approaching H  i line flux comes very close to the inner lens caustic and results in very high magnification. These results, for the first time, demonstrate the feasibility of observing high-redshift H  i in a lensed system with the modest amount of telescope time and open up exciting new possibilities for probing the cosmic evolution of neutral gas with existing and upcoming low-frequency radio telescopes in the near future. [ABSTRACT FROM AUTHOR]

Galaxies and galaxy groups located along the line of sight towards gravitationally lensed quasars produce high-order perturbations of the gravitational potential at the lens position. When these perturbation are too large, they can induce a systematic error on H0 of a few per cent if the lens system is used for cosmological inference and the perturbers are not explicitly accounted for in the lens model. In this work, we present a detailed characterization of the environment of the lens system WFI 2033−4723 (zsrc = 1.662, zlens = 0.6575), one of the core targets of the H0LiCOW project for which we present cosmological inferences in a companion paper. We use the Gemini and ESO-Very Large telescopes to measure the spectroscopic redshifts of the brightest galaxies towards the lens, and use the ESO-MUSE integral field spectrograph to measure the velocity-dispersion of the lens (σlos = 250+−1521 km s−1) and of several nearby galaxies. In addition, we measure photometric redshifts and stellar masses of all galaxies down to i < 23 mag, mainly based on Dark Energy Survey imaging (DR1). Our new catalogue, complemented with literature data, more than doubles the number of known galaxy spectroscopic redshifts in the direct vicinity of the lens, expanding to 116 (64) the number of spectroscopic redshifts for galaxies separated by less than 3 arcmin (2 arcmin) from the lens. Using the flexion-shift as a measure of the amplitude of the gravitational perturbation, we identify two galaxy groups and three galaxies that require specific attention in the lens models. The ESO MUSE data enable us to measure the velocity-dispersions of three of these galaxies. These results are essential for the cosmological inference analysis presented in Rusu et al.

We report the first measurements with sub-kiloparsec spatial resolution of strongly inverted gas-phase metallicity gradients in two dwarf galaxies at z ∼ 2. The galaxies have stellar masses ∼109 , specific star formation rate ∼20 Gyr-1, and global metallicity (1/4 solar), assuming the strong-line calibrations of [O iii]/Hβ and [O ii]/Hβ from Maiolino et al. Their radial metallicity gradients are measured to be highly inverted, i.e., 0.122 ± 0.008 and 0.111 ± 0.017 dex kpc-1, which is hitherto unseen at such small masses in similar redshift ranges. From the Hubble Space Telescope observations of the source nebular emission and stellar continuum, we present two-dimensional spatial maps of star formation rate surface density, stellar population age, and gas fraction, which show that our galaxies are currently undergoing rapid mass assembly via disk inside-out growth. More importantly, using a simple chemical evolution model, we find that the gas fractions for different metallicity regions cannot be explained by pure gas accretion. Our spatially resolved analysis based on a more advanced gas regulator model results in a spatial map of net gaseous outflows, triggered by active central starbursts, that potentially play a significant role in shaping the spatial distribution of metallicity by effectively transporting stellar nucleosynthesis yields outwards. The relation between wind mass loading factors and stellar surface densities measured in different regions of our galaxies shows that a single type of wind mechanism, driven by either energy or momentum conservation, cannot explain the entire galaxy. These sources present a unique constraint on the effects of gas flows on the early phase of disk growth from the perspective of spatially resolved chemical evolution within individual systems.

We test the ability of strong lensing data to constrain the size of a central core in the dark matter haloes of galaxy clusters, using Abell 611 as a prototype. Using simulated data, we show that modelling a cluster halo with ellipticity in the gravitational potential can bias the inferred mass and concentration, which may bias the inferred central density when weak lensing or X-ray data are added. We also highlight the possibility for spurious constraints on the core size if the radial density profile is different from the assumed model. These systematics can be ameliorated if central images are present in the data. Applying our methodology to Abell 611 and imposing a reasonable prior on the stellar mass-to-light ratio restricts the core size to be less than about 4 kpc, with a minimum reduced χ2 of 0.28, and with RMS error of 0″.12 in image positions. Such small cores imply a constraint on the dark matter self-interaction cross section of the order of 0.1 cm2, g-1 at relative velocities of about 1500 km s-1. In addition, our results imply a bottom-heavy (super-Salpeter) initial mass function for the central galaxy in the Abell 611.

Supernovae that are strongly gravitationally lensed (gLSNe) by elliptical galaxies are powerful probes of astrophysics and cosmology that will be discovered systematically by wide-field, high-cadence imaging surveys such as the Zwicky Transient Facility (ZTF) and the Large Synoptic Survey Telescope (LSST). Here we use pixel-level simulations that include observing strategy, target selection, supernova properties, and dust to forecast the rates and properties of gLSNe that ZTF and LSST will find. Applying the resolution-insensitive discovery strategy of Goldstein et al., we forecast that ZTF (LSST) can discover 0.02 (0.79) 91bg-like, 0.17 (5.92) 91T-like, 1.22 (47.84) Type Ia, 2.76 (88.51) Type IIP, 0.31 (12.78) Type IIL, and 0.36 (15.43) Type Ib/c gLSNe per year, with uncertainties dominated by uncertainties in the supernova rate. We also forecast that the surveys can discover at least 3.75 (209.32) Type IIn gLSNe per year, for a total of at least 8.60 (380.60) gLSNe per year under fiducial observing strategies. ZTF gLSNe have a median z s = 0.9, z l = 0.35, , Δt max = 10 days, min(θ) = 0.″25, and N img = 4. LSST gLSNe are less compact and less magnified, with a median z s = 1.0, z l = 0.4, , Δt max = 25 days, min(θ) = 0.″6, and N img = 2. We develop a model of the supernova-host galaxy connection and find that the vast majority of gLSN host galaxies will be multiply imaged, enabling detailed constraints on lens models with sufficiently deep high-resolution imaging taken after the supernova has faded. We release the results of our simulations as catalogs at http://portal.nersc.gov/project/astro250/glsne/.

We search Dark Energy Survey (DES) Year 3 imaging data for galaxy-galaxy strong gravitational lenses using convolutional neural networks. We generate 250 000 simulated lenses at redshifts > 0.8 from which we create a data set for training the neural networks with realistic seeing, sky and shot noise. Using the simulations as a guide, we build a catalogue of 1.1 million DES sources with 1.8 < g − i < 5, 0.6 < g − r < 3, r mag > 19, g mag > 20, and i mag > 18.2. We train two ensembles of neural networks on training sets consisting of simulated lenses, simulated non-lenses, and real sources. We use the neural networks to score images of each of the sources in our catalogue with a value from 0 to 1, and select those with scores greater than a chosen threshold for visual inspection, resulting in a candidate set of 7301 galaxies. During visual inspection, we rate 84 as 'probably' or 'definitely' lenses. Four of these are previously known lenses or lens candidates. We inspect a further 9428 candidates with a different score threshold, and identify four new candidates. We present 84 new strong lens candidates, selected after a few hours of visual inspection by astronomers. This catalogue contains a comparable number of high-redshift lenses to that predicted by simulations. Based on simulations, we estimate our sample to contain most discoverable lenses in this imaging and at this redshift range.

We present spatially resolved analysis of a lensed galaxy, SDSS1958+5950 at z = 2.225, from the Cambridge Sloan Survey ofWide Arcs in the Sky. We use our new high-resolution imaging data to construct a robust lens model for the galaxy group at z = 0.214.We employ the updated lens model to combine the Integral Field Spectrographic observations on two highly distorted images of the lensed target. We adopt a forward-modelling approach to deconvolve the effects of point spread function from the combined source-plane reconstruction. The approach is adapted to the lens model magnification and enables a resolution of ~170 pc in the galaxysource plane. We propose an ongoing merger as the origin of the lensed system on the basis of its source-plane morphology, kinematics and rest-frame emission-line ratios. Using our novel technique of adaptive coadded source plane reconstruction, we are able to detect different components in the velocity gradient that were not seen in previous studies of this object, plausibly belonging to different components in the merging system.

We present spatially resolved analysis of a lensed galaxy, SDSS1958+5950 at z = 2.225, from the Cambridge Sloan Survey ofWide Arcs in the Sky. We use our new high-resolution imaging data to construct a robust lens model for the galaxy group at z = 0.214.We employ the updated lens model to combine the Integral Field Spectrographic observations on two highly distorted images of the lensed target. We adopt a forward-modelling approach to deconvolve the effects of point spread function from the combined source-plane reconstruction. The approach is adapted to the lens model magnification and enables a resolution of ~170 pc in the galaxysource plane. We propose an ongoing merger as the origin of the lensed system on the basis of its source-plane morphology, kinematics and rest-frame emission-line ratios. Using our novel technique of adaptive coadded source plane reconstruction, we are able to detect different components in the velocity gradient that were not seen in previous studies of this object, plausibly belonging to different components in the merging system.

The first deep-field observations of the JWST have immediately yielded a surprisingly large number of very high redshift candidates, pushing the frontier of observability well beyond z ≳ 10. We here present a detailed SED-fitting analysis of the 10 gravitationally lensed z ∼ 9–16 galaxy candidates detected behind the galaxy cluster SMACS J0723.3−7327 in a previous paper using the BEAGLE tool. Our analysis makes use of dynamical considerations to place limits on the ages of these galaxies and of all three published SL models of the cluster to account for lensing systematics. We find the majority of these galaxies to have relatively low stellar masses |$M_{\star }\sim 10^7-10^8\, \mathrm{M}_{\odot }$| and young ages t age ∼ 10–100 Myr but with a few higher mass exceptions (⁠|$M_{\star }\sim 10^9\rm{-}10^{10}\, \mathrm{M}_{\odot }$|⁠) due to Balmer-break detections at z ∼ 9–10. Because of their very blue UV-slopes, down to β ∼ −3, all of the galaxies in our sample have extremely low dust attenuations AV ≲ 0.02. Placing the measured parameters into relation, we find a very shallow M ⋆ − M UV-slope and high sSFRs above the main sequence of star formation with no significant redshift-evolution in either relation. This is in agreement with the bright UV luminosities measured for these objects and indicates that we are naturally selecting UV-bright galaxies that are undergoing intense star formation at the time they are observed. Finally, we discuss the robustness of our high-redshift galaxy sample regarding low-redshift interlopers and conclude that low-redshift solutions can safely be ruled out for roughly half of the sample, including the highest redshift galaxies at z ∼ 12–16. These objects represent compelling targets for spectroscopic follow-up observations with JWST and ALMA. [ABSTRACT FROM AUTHOR]

The sample of time-delay gravitational lenses appropriate for studying the geometry of the Universe continues to grow as dedicated campaigns, such as the Dark Energy Survey, the (Very Large Telescope) VLT Survey Telescope (VST) ATLAS survey, and the Large Synoptic Survey Telescope, complete their census of high-redshift sources. This catalogue now includes hundreds of strong lensing systems, at least 31 of which have reasonably accurate time-delay measurements. In this paper, we use them to compare the predictions of two competing Friedmann–Lemaître–Robertson–Walker models: flat Lambda cold dark matter (ΛCDM), characterized by two adjustable parameters (H 0 and Ωm), and the R h = ct universe (with H 0 as the single free variable). Over the past decade, the latter has accounted for the data better than the standard model, most recently the emergence of well-formed galaxies discovered by JWST at cosmic dawn. Here, we show that the current sample of time-delay lenses favours R h = ct with a likelihood of ∼ |$84{{\ \rm per\ cent}}$| versus ∼ |$16{{\ \rm per\ cent}}$| for the standard model. This level of accuracy will greatly improve as the ongoing surveys uncover many thousands of additional lens systems over the next several years. [ABSTRACT FROM AUTHOR]

When travelling from their source to the observer, gravitational waves can get deflected by massive objects along their travel path. For a massive lens and a good source-lens alignment, the wave undergoes strong lensing, leading to several images with the same frequency evolution. These images are separated in time, magnified, and can undergo an overall phase shift. Searches for strongly lensed gravitational waves look for events with similar masses, spins, and sky location and linked through so-called lensing parameters. However, the agreement between these quantities can also happen by chance. To reduce the overlap between background and foreground, one can include lensing models. When doing realistic searches, one does not know which model is the correct one to be used. Using an incorrect model could lead to the non-detection of genuinely lensed events. In this work, we investigate how one can reduce the false alarm probability when searching for strongly lensed events. We focus on the impact of the addition of a model for the lens density profile and investigate the effect of potential errors in the modelling. We show that the risks of false alarm are high without the addition of a lens model. We also show that slight variations in the profile of the lens model are tolerable, but a model with an incorrect assumption about the underlying lens population causes significant errors in the identification process. We also suggest some strategies to improve confidence in the detection of strongly lensed gravitational waves. [ABSTRACT FROM AUTHOR]

In general, the cosmological redshift of an object changes with time, a phenomenon known as redshift drift. Although long known theoretically, recently interest has been renewed because of the possibility of measuring it on a reasonable time-scale. Strong gravitational lensing offers a possibility to measure it on a much shorter time-scale, by making use of a time delay of n years rather than making observations separated by n years, but, perhaps at least in part because of the expectation that the signal would be swamped by a larger change in redshift due to transverse motion of the lens, that has not attracted much interest. I present a method to extract the small signal, making use of the fact that the light-travel time through different parts of an Einstein ring is the same (and hence the difference in redshift due to redshift drift vanishes), thus enabling the measurement of redshift drift on a much shorter time-scale, and show how that can help in the measurement of the distribution of (dark) matter. [ABSTRACT FROM AUTHOR]

One of the main goals of the JWST is to study the first galaxies in the Universe. We present a systematic photometric analysis of very distant galaxies in the first JWST deep field towards the massive lensing cluster SMACS0723. As a result, we report the discovery of two galaxy candidates at z ∼ 16, only 250 million years after the big bang. We also identify two candidates at z ∼ 12 and six candidates at z ∼ 9−11. Our search extended out to z ≲ 21 by combining colour information across seven near-infrared camera and near-infrared imager and slitless spectrograph filters. By modelling the Spectral Energy Distributions (SEDs) with EAZY and BEAGLE , we test the robustness of the photometric redshift estimates. While their intrinsic (unlensed) luminosity is typical of the characteristic luminosity L* at z > 10, our high-redshift galaxies typically show small sizes and their morphologies are consistent with disks in some cases. The highest-redshift candidates have extremely blue ultraviolet-continuum slopes −3 < β < −2.4, young ages ∼10−100 Myr, and stellar masses around log (M ⋆/M⊙) = 8.8 inferred from their spectral energy distribution modelling, which indicate a rapid build-up of their stellar mass. Our search clearly demonstrates the capabilities of JWST to uncover robust photometric candidates up to very high redshifts and peer into the formation epoch of the first galaxies. [ABSTRACT FROM AUTHOR]

Future ground-based gravitational wave (GW) detectors, i.e. Einstein telescope (ET) and Cosmic Explorer (CE), are expected to detect a significant number of lensed binary neutron star (BNS) mergers, which may provide a unique tool to probe cosmology. In this paper, we investigate the detectability of the optical/infrared electromagnetic (EM) counterparts (kilonovae/afterglows) from these lensed BNS mergers by future GW detectors and EM telescopes using simple kilonova, afterglow, and lens models. ET and CE are expected to detect |${\sim}5.32^{+26.1}_{-5.10}$| and |$67.3^{+332}_{-64.7}$| lensed BNS mergers per year. We find that the EM counterparts associated with all these mergers will be detectable by an all sky-survey in the H band with the limiting magnitude m lim ≳ 27, while the detectable fraction is |${\lesssim}0.4{{\ \rm per\ cent}}$| in the g / z band if with m lim ≲ 24. Generally, it is more efficient to search the lensed EM counterparts by adopting the infrared bands than the optical/UV bands with the same m lim. Future telescopes like Vera C. Rubin Observatory, China Space Station Telescope, and Euclid can hardly detect the EM counterparts of even one lensed BNS merger. Roman Space Telescope (RST) and James Webb Space Telescope (JWST) have the capability to detect about a few or more such events per year. Moreover, the time delays and separations between the lensed image pairs are typically in the ranges from minutes to months and from 0.1 to 1 arcsec, suggesting that both the GW and EM images of most lensed BNS mergers can be well resolved by not only CE/ET in the time domain but also RST / JWST spatially. [ABSTRACT FROM AUTHOR]

Strongly lensed systems are powerful probes of the distribution of dark matter on small scales. In this paper, we show that line-of-sight haloes between the source and the observers give rise to a distinct anisotropic signature in the two-point function of the effective lensing deflection field. We show in particular that the non-linear coupling between line-of-sight haloes and the main lens plane imprints a characteristic quadrupole moment on this two-point function whose amplitude reflects the abundance of such haloes within the strongly lensed field. We discuss how, by taking ratios of different multipole moments, such observables could be made robust under the mass-sheet transform. We also demonstrate that future extremely large telescopes have the ability to detect the quadrupole moment due to this unique anisotropic signature under ideal conditions. Our approach opens the door to statistically distinguish the effect of line-of-sight haloes from that of the main-lens substructure on lensed images, hence allowing one to probe dark matter physics in a new way. [ABSTRACT FROM AUTHOR]

In this work we present relensing , a package written in python whose goal is to model galaxy clusters from gravitational lensing. With relensing we extend the amount of software available, which provides the scientific community with a wide range of models that help us to compare and therefore validate the physical results that rely on them. We implement a free-form approach which computes the gravitational deflection potential on an adaptive irregular grid, from which one can characterize the cluster and its properties as a gravitational lens. Here, we use two alternative penalty functions to constrain strong lensing. We apply relensing to two toy models, in order to explore under which conditions one can get a better performance in the reconstruction. We find that by applying a smoothing to the deflection potential, we are able to increase the capability of this approach to recover the shape and size of the mass profile of galaxy clusters, as well as its magnification map. This translates into a better estimation of the critical and caustic curves. The power that the smoothing provides is also tested on the simulated clusters Ares and Hera, for which we get an rms on the lens plane of |$\sim 0.17\, {\rm arcsec}$| and |$\sim 0.16\, {\rm arcsec}$|⁠ , respectively. Our results represent an improvement with respect to reconstructions that were carried out with methods of the same nature as relensing. In its current state, relensing is available upon request. [ABSTRACT FROM AUTHOR]

Precision analysis of galaxy–galaxy strong gravitational lensing images provides a unique way of characterizing small-scale dark matter haloes, and could allow us to uncover the fundamental properties of dark matter's constituents. Recently, gravitational imaging techniques made it possible to detect a few heavy subhaloes. However, gravitational lenses contain numerous subhaloes and line-of-sight haloes, whose subtle imprint is extremely difficult to detect individually. Existing methods for marginalizing over this large population of subthreshold perturbers to infer population-level parameters are typically computationally expensive, or require compressing observations into hand-crafted summary statistics, such as a power spectrum of residuals. Here, we present the first analysis pipeline to combine parametric lensing models and a recently developed neural simulation-based inference technique called truncated marginal neural ratio estimation (TMNRE) to constrain the warm dark matter halo mass function cut-off scale directly from multiple lensing images. Through a proof-of-concept application to simulated data, we show that our approach enables empirically testable inference of the dark matter cut-off mass through marginalization over a large population of realistic perturbers that would be undetectable on their own, and over lens and source parameter uncertainties. To obtain our results, we combine the signal contained in a set of images with Hubble Space Telescope resolution. Our results suggest that TMNRE can be a powerful approach to put tight constraints on the mass of warm dark matter in the multi-keV regime, which will be relevant both for existing lensing data and in the large sample of lenses that will be delivered by near-future telescopes. [ABSTRACT FROM AUTHOR]

This study explores the gravitational lensing effects of supermassive black holes (SMBHs) in galaxy clusters. While the presence of central SMBHs in galaxies is firmly established, recent work from high-resolution simulations predict the existence of an additional population of wandering SMBHs. Though the masses of these SMBHs are a minor perturbation on the larger scale and individual galaxy scale dark matter components in the cluster, they can impact statistical lensing properties and individual lensed image configurations. Probing for these potentially observable signatures, we find that SMBHs imprint detectable signatures in rare, higher order strong lensing image configurations although they do not manifest any statistically significant detectable evidence in either the magnification distribution or the integrated shear profile. Investigating specific lensed image geometries, we report that a massive, near point-like, potential of an SMBH causes the following detectable effects: (i) image splitting leading to the generation of extra images; (ii) positional and magnification asymmetries in multiply imaged systems; and (iii) the apparent disappearance of a lensed counter image. Of these, image splitting inside the cluster tangential critical curve, is the most prevalent notable observational signature. We demonstrate these possibilities in two cases of observed giant arcs in SGAS J003341.5+024217 and RX J1347.5−1145, wherein specific image configurations seen can be reproduced with SMBHs. Future observations with high-resolution instrumentation (e.g. MAVIS-Very Large Telescope, MICADO-Extremely Large Telescope, and the upgraded ngVLA, along with data from the Euclid and Nancy Grace Roman Space Telescopes and the Rubin LSST Observatory are likely to allow us to probe these unique yet rare SMBHs lensing signatures. [ABSTRACT FROM AUTHOR]

Strong gravitational lensing offers a compelling test of the cold dark matter paradigm, as it allows for subhaloes with masses of ∼109 M⊙ 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 ∼1013 M⊙ 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; however, for the latter we find false positives and the inferred halo mass is overestimated by ∼4–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. [ABSTRACT FROM AUTHOR]

Gravitational time delays provide a powerful one-step measurement of H 0, independent of all other probes. One key ingredient in time-delay cosmography are high-accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (105–106 CPU hours and ∼0.5–1 yr, respectively). Major improvements in modelling speed are therefore necessary to exploit the large number of lenses that are forecast to be discovered over the current decade. In order to bypass this roadblock, we develop an automated modelling pipeline and apply it to a sample of 31 lens systems, observed by the Hubble Space Telescope in multiple bands. Our automated pipeline can derive models for 30/31 lenses with few hours of human time and <100 CPU hours of computing time for a typical system. For each lens, we provide measurements of key parameters and predictions of magnification as well as time delays for the multiple images. We characterize the cosmography-readiness of our models using the stability of differences in the Fermat potential (proportional to time delay) with respect to modelling choices. We find that for 10/30 lenses, our models are cosmography or nearly cosmography grade (<3 per cent and 3–5 per cent variations). For 6/30 lenses, the models are close to cosmography grade (5–10 per cent). These results utilize informative priors and will need to be confirmed by further analysis. However, they are also likely to improve by extending the pipeline modelling sequence and options. In conclusion, we show that uniform cosmography grade modelling of large strong lens samples is within reach. [ABSTRACT FROM AUTHOR]

Fast cosmological transients such as fast radio bursts (FRBs) and gamma-ray bursts (GRBs) represent a class of sources more compact than any other cosmological object. As such, they are sensitive to significant magnification via gravitational lensing from a class of lenses which are not well constrained by observations today. Low-mass primordial black holes are one such candidate which may constitute a significant fraction of the Universe's dark matter. Current observations only constrain their density in the nearby Universe, giving fast transients from cosmological distances the potential to form complementary constraints. Motivated by this, we calculate the effect that gravitational lensing from a cosmological distribution of compact objects would have on the observed rates of FRBs and GRBs. For static lensing geometries, we rule out the prospect that all FRBs are gravitationally lensed for a range of lens masses and show that lens masses greater than 10−5 M ⊙ can be constrained with 8000 unlocalized high-fluence FRBs at 1.4GHz, as might be detected by the next generation of FRB-finding telescopes. [ABSTRACT FROM AUTHOR]

Over the last two decades, strong gravitational lensing has emerged as a potential method for studying the nature and distribution of dark matter on sub-galactic scales. In addition to the main lens substructure, line-of-sight dark matter haloes contribute greatly to the subtle perturbations of lensed images. Line-of-sight haloes, unlike dark matter subhaloes, imprint distinct anisotropic and quadrupole signatures in the maps that depict the divergence and curl of the effective deflection field, respectively, giving rise to quadrupole moments of the image-plane averaged two-point correlation function of these maps. In terms of central density evolution and dark matter halo distribution, the shapes and amplitudes of the two-point function multipoles alter dramatically in the presence of warm dark matter and self-interacting dark matter. This method, in conjunction with upcoming large-scale surveys, provides the prospect of improving the constraints on dark matter at a critical time in strong gravitational lensing research. [ABSTRACT FROM AUTHOR]

One of the largest sources of systematics in time-delay cosmography arises from Mass Sheet Transformation (MST). The degeneracy associated with this transformation is often broken by an assumed profile shape, such as a power-law. A hierarchical strategy has been developed which constrains the global profile shape on a population level, constrained collectively by the kinematics measurements of the lenses. This framework allows one to include non-time-delay lenses to provide constraints to the global profile, improving the H 0 constraints. This work tests the hierarchical framework using analytical profiles, and additionally tests the capacity to combine two populations which come from the same profiles but probe different radii due to a change in source redshift. We find that the hierarchical framework is able to compensate for this effect, and the addition of non-time-delay lenses improves the H 0 constraint, even though these lenses have different Einstein radii than their time-delay counterparts. [ABSTRACT FROM AUTHOR]

Galaxy-galaxy strong lensing (SL) systems provide a unique opportunity to test modified gravity theories. Deviations from General Relativity are encoded in the post-Newtonian parameter (γ). As a preparation for the upcoming data from the Vera Rubin Observatory Legacy Survey of Space and Time (LSST), our research group collected imaging data of SL systems from ground-based telescopes and conducted spectroscopic observations of 21 systems on the Southern Astrophysical Research (SOAR) Telescope to measure the lens velocity dispersions, σ v . We briefly describe the semi-automated SL modelling of the systems in this sample and combine the results with σ v from SOAR to derive an estimate for γ. Our preliminary results yield a value of $$\gamma = 1.17_{ - 0.33}^{ + 0.29}$$ , which is consistent with General Relativity. Although the error bars are limited by the sample size, this result represents the first constraint on modified gravity obtained purely from ground-based data, with a sample completely independent from previous studies, and which allows for a self consistent end-to-end analysis. [ABSTRACT FROM AUTHOR]

Current searches for galaxy-scale strong lenses focus on massive Luminous Red Galaxies but tend to overlook late-type lenses, in part because of their smaller Einstein radii. We take advantage of the superb seeing of the UNIONS survey in the r -band to perform an imaging search for edge-on late-type lenses. We use Convolutional Neural Networks trained with simulated observations composed of images of real galaxies from UNIONS and real sources from HST. Using 3600 square degrees of the survey we test ∼7 million galaxies and find 56 systems with obvious signs of lensing. In addition, we empirically estimate the true prevalence of lenses in UNIONS by visually inspecting 120,000 randomly chosen images in the survey. We find that the number of edge-on lenses we discover with CNNs is compatible with these estimates. [ABSTRACT FROM AUTHOR]

Rare systems with multiple sources strongly lensed by a single galaxy provide the most robust way to explore the geometry of the Universe and to study the inner mass structure of lens galaxies. We present here a study of the SDSS J0100+1818 deflector, analyzing its total and baryonic mass distributions. The system comprises an ultra-massive early-type galaxy, surrounded by ten multiple images of three background sources. Exploiting high-resolution HST photometry and VLT/X-shooter spectroscopy we conduct a strong lensing analysis with the software GLEE to reconstruct the complex surface brightness distributions of the background sources over approximately 7200 HST pixels. These results are presented in our recent paper, Bolamperti et al. (2023). Finally, we present some preliminary results from new VLT/MUSE observations, that will allow us to build a new strong lensing+dynamics joint model and measure the values of the total matter density and of the cosmological constant parameters, Ω m and ΩΛ. [ABSTRACT FROM AUTHOR]

We aim to uncover the structure of the continuum and broad emission line (BEL) emitting regions in the gravitationally lensed quasar SDSS J1004+4112 through unique microlensing signatures. Analyzing 20 spectroscopic observations from 2003 to 2018, we study the striking deformations of various BEL profiles and determine the sizes of their respective emitting regions. We use the emission line cores as a baseline for no microlensing and then apply Bayesian methods to derive the sizes of the Ly α , Si IV, C IV, C III], and Mg II emitting regions, as well as of the underlying continuum-emitting sources. We find that the sizes of the emitting regions for the BELs are a few light-days across, notably smaller than in typical lensed quasars. The asymmetric distortions observed in the BELs suggest that the broad-line region lacks spherical symmetry and is likely confined to a plane. The inferred continuum emitting region sizes are larger than predictions based on standard thin-disk theory by a factor of ∼ 4. We find that the size-wavelength relation is in agreement with that of a geometrically thin and optically thick accretion disk. [ABSTRACT FROM AUTHOR]

The recent rapid growth of the black hole (BH) catalogue from gravitational waves (GWs) has allowed us to study the substructure of black hole mass function (BHMF) beyond the simplest power–law distribution. However, the BH masses inferred from binary BH merger events, may be systematically 'brightened' or 'dimmed' by the gravitational lensing effect. In this work, we investigate the impact of gravitational lensing on the BHMF inference considering the detection of the third-generation GW detector – the Einstein Telescope (ET). We focus on high redshift, z  = 10 in order to obtain the upper limits of this effect. We use Monte Carlo (MC) method to simulate the data adopting three original BHMFs under Un-Lensed and Lensed scenarios, then recover the parameters of BHMFs from the mock data, and compare the difference of results, respectively. We found that all the parameters are well recovered within one standard deviation(std. 1σ), and all three BHMF models are reconstructed within 68 per cent credible interval, suggesting that lensing would not change the main structure drastically, even at very high redshifts and with high precision of ET. And the modest influence beyond 50 M ⊙, depends on the modeling of the high mass tail or substructure of BHMF. We conclude that the impact of lensing on BHMF inference with ET can be safely ignored in the foreseeable future. Careful handling of lensing effects is required only when focusing on an accurate estimation of the high mass end of BHMF at high redshifts. [ABSTRACT FROM AUTHOR]

Strong gravitational lensing has emerged as a promising approach for probing dark matter (DM) models on sub-galactic scales. Recent work has proposed the subhalo effective density slope as a more reliable observable than the commonly used subhalo mass function. The subhalo effective density slope is a measurement independent of assumptions about the underlying density profile and can be inferred for individual subhaloes through traditional sampling methods. To go beyond individual subhalo measurements, we leverage recent advances in machine learning and introduce a neural likelihood-ratio estimator to infer an effective density slope for populations of subhaloes. We demonstrate that our method is capable of harnessing the statistical power of multiple subhaloes (within and across multiple images) to distinguish between characteristics of different subhalo populations. The computational efficiency warranted by the neural likelihood-ratio estimator over traditional sampling enables statistical studies of DM perturbers and is particularly useful as we expect an influx of strong lensing systems from upcoming surveys. [ABSTRACT FROM AUTHOR]

Primordial gravitational waves with extremely low frequency are expected to origin from inflation in the early Universe. The detection of such kind of gravitational waves is of great significance to verify the inflationary theory and determine the energy scale of inflation. The B mode of polarization of cosmic microwave background shows to be a promising probe of extremely low frequency primordial gravitational wave. In order to find a complementary observational feature induced by primordial gravitational wave with extremely low frequency, we propose an alternative way of detection by investigating the effect of primordial gravitational wave with extremely low frequency on a gravitational lens system with a non-aligned source–deflector–observer configuration. The results show that, with a series of chosen parameters, a gravitational lens system with perturbation from extremely low frequency primordial gravitational wave could induce time delay that could deviate from the time delay deduced from the theoretical model as much as 100 per cent, meaning that a gravitational lens system with a non-aligned configuration could serve as a potential long-baseline detector of extremely low frequency primordial gravitational wave. [ABSTRACT FROM AUTHOR]