Search

Gravitational lensing may render individual high-mass stars detectable out to cosmological distances, and several extremely magnified stars have in recent years been detected out to redshifts |$z\approx 6$|⁠. Here, we present Muspelheim, a model for the evolving spectral energy distributions of both metal-enriched and metal-free stars at high redshifts. Using this model, we argue that lensed stars will form a highly biased sample of the intrinsic distribution of stars across the Hertzsprung–Russell diagram, and that this bias will typically tend to favour the detection of lensed stars in evolved stages characterized by low effective temperatures, even though stars only spend a minor fraction of their lifetimes in such states. We also explore the prospects of detecting individual, lensed metal-free (Population III) stars at high redshifts using the James Webb Space Telescope (JWST). We find that very massive (⁠|$\gtrsim 100\ \mathrm{M}_\odot$|⁠) Population III stars at |$z\gtrsim 6$| may potentially be detected by JWST in surveys covering large numbers of strong-lensing clusters, provided that the Population III stellar initial mass function is sufficiently top-heavy, that these stars evolve to effective temperatures |$\le 15000$| K, and that the cosmic star formation rate density of Pop III stars reaches |$\gtrsim 10^{-4}\ \mathrm{M}_\odot$| cMpc |$^{-3}$| yr |$^{-1}$| at |$z\approx 6$| –10. Various ways to distinguish metal-free lensed stars from metal-enriched ones are also discussed. [ABSTRACT FROM AUTHOR]

MACS J0600.1-2008 (MACS0600) is an X-ray-luminous, massive galaxy cluster at |$z_{\mathrm{d}}=0.43$|⁠ , studied previously by the REionization LensIng Cluster Survey and ALMA Lensing Cluster Survey projects which revealed a complex, bimodal mass distribution and an intriguing high-redshift object behind it. Here, we report on the results of a combined analysis of the extended strong lensing (SL), X-ray, Sunyaev–Zeldovich (SZ), and galaxy luminosity-density properties of this system. Using new JWST and ground-based Gemini-N and Keck data, we obtain 13 new spectroscopic redshifts of multiply-imaged galaxies and identify 12 new photometric multiple-image systems and candidates, including two multiply-imaged |$z\sim 7$| objects. Taking advantage of the larger areal coverage, our analysis reveals an additional bimodal, massive SL structure which we measure spectroscopically to lie adjacent to the cluster and whose existence was implied by previous SL-modelling analyses. While based in part on photometric systems identified in ground-based imaging requiring further verification, our extended SL model suggests that the cluster may have the second-largest critical area and effective Einstein radius observed to date, |$A_{\mathrm{crit}}\simeq 2.16\, \mathrm{arcmin}^2$| and |$\theta _{\mathrm{E}}=49.7^{\prime \prime }\pm 5.0^{\prime \prime }$| for a source at |$z_{\mathrm{s}}=2$|⁠ , enclosing a total mass of |$M(\lt \theta _{\mathrm{E}})=(4.7\pm 0.7)\times 10^{14}\, \mathrm{M}_{\odot }$|⁠. These results are also supported by the galaxy luminosity distribution, and the SZ and X-ray data. Yet another, probably related massive cluster structure, discovered in X-rays 5 arcmin (1.7 Mpc) further north, suggests that MACS0600 is part of an even larger filamentary structure. This discovery adds to several recent detections of massive structures around SL galaxy clusters and establishes MACS0600 as a prime target for future high-redshift surveys with JWST. [ABSTRACT FROM AUTHOR]

In the context of upcoming large-scale surveys like Euclid, the necessity for the automation of strong lens detection is essential. While existing machine learning pipelines heavily rely on the classification probability (P), this study intends to address the importance of integrating additional metrics, such as Information Content (IC) and the number of pixels above the segmentation threshold (⁠|$\rm {\mathit{n}_{s}}$|⁠), to alleviate the false positive rate in unbalanced data-sets. In this work, we introduce a segmentation algorithm (U-Net) as a supplementary step in the established strong gravitational lens identification pipeline (Denselens), which primarily utilizes |$\rm {\mathit{P}_{mean}}$| and |$\rm {IC_{mean}}$| parameters for the detection and ranking. The results demonstrate that the inclusion of segmentation enables significant reduction of false positives by approximately 25 per cent in the final sample extracted from DenseLens, without compromising the identification of strong lenses. The main objective of this study is to automate the strong lens detection process by integrating these three metrics. To achieve this, a decision tree-based selection process is introduced, applied to the Kilo Degree Survey (KiDS) data. This process involves rank-ordering based on classification scores (⁠|$\rm {\mathit{P}_{mean}}$|⁠), filtering based on Information Content (⁠|$\rm {IC_{mean}}$|⁠), and segmentation score (⁠|$\rm {n_{s}}$|⁠). Additionally, the study presents 14 newly discovered strong lensing candidates identified by the U-Denselens network using the KiDS DR4 data. [ABSTRACT FROM AUTHOR]

Galaxy–galaxy strong gravitational lensing (GGSL) is a powerful probe for the formation and evolution of galaxies and cosmology, while the sample size of GGSLs leads to considerable uncertainties and potential bias. The China Space Station Telescope (CSST, to be launched in late 2026) will conduct observations across 17 500 square degrees of the sky, capturing images in the |$ugriz$| bands with a spatial resolution comparable to that of the Hubble Space Telescope. We ran a set of Monte Carlo simulations to predict that the CSST's wide-field survey will observe |$\sim$| 160 000 galaxy–galaxy strong lenses over its lifespan, increasing the number of existing galaxy–galaxy strong lens samples by three orders of magnitude. This is comparable to the capabilities of the |$\it Euclid$| telescope but with the added benefit of additional colour information. Specifically, the CSST can detect strong lenses with Einstein radii about |$0.64\pm 0.42 ^{\prime \prime }$|⁠ , corresponding to the velocity dispersions of |$217.19 \pm 50.55 \, \text{km/s}$|⁠. These lenses exhibit a median magnification of |$\sim$| 5. The apparent magnitude of the unlensed sources in the g band is |$25.87 \pm 1.19$|⁠. The signal-to-noise ratio of the lensed images covers a range of |${\sim} 20$| to |${\sim} 1000$|⁠ , allowing us to determine the Einstein radius with an accuracy ranging from |${\sim} 1~{{\ \rm per\ cent}}$| to |${\sim} 0.1~{{\ \rm per\ cent}}$|⁠ , ignoring various modelling systematics. Our estimates indicate that CSST can observe rare systems like double source-plane and spiral galaxy lenses. The above selection functions of the CSST strong lensing observation help optimize the strategy of finding and modelling GGSLs. [ABSTRACT FROM AUTHOR]

Strong gravitational lensing provides a purely gravitational means to infer properties of dark matter haloes and thereby constrain the particle nature of dark matter. Strong lenses sometimes appear as four lensed images of a background quasar accompanied by spatially resolved emission from the quasar host galaxy encircling the main deflector (lensed arcs). We present methodology to simultaneously reconstruct lensed arcs and relative image magnifications (flux ratios) in the presence of full populations of subhaloes and line-of-sight haloes. To this end, we develop a new approach for multiplane ray tracing that accelerates lens mass and source light reconstruction by factors of |$\sim\!\! 100\!\!-\!\!1000$|⁠. Using simulated data, we show that simultaneous reconstruction of lensed arcs and flux ratios isolates small-scale perturbations to flux ratios by dark matter substructure from uncertainties associated with the main deflector mass profile on larger angular scales. Relative to analyses that use only image positions and flux ratios to constrain the lens model, incorporating arcs strengthens likelihood ratios penalizing warm dark matter with a suppression scale |$m_{\rm {hm}} / {\rm M}_{\odot }$| in the ranges of |$\left[10^7 \!\!-\!\! 10^{7.5}\right]$|⁠ , |$\left[10^{7.5} \!\!-\!\! 10^{8}\right]$|⁠ , |$\left[10^8 \!\!-\!\! 10^{8.5}\right]$|⁠ , and |$\left[10^{8.5} \!\!-\!\! 10^{9}\right]$| by factors of 1.3, 2.5, 5.6, and 13.1, respectively, for a cold dark matter ground truth. The 95 per cent exclusion limit improves by 0.5 dex in |$\log _{10} m_{\rm {hm}}$|⁠. The enhanced sensitivity to low-mass haloes enabled by these methods pushes the observational frontier of substructure lensing to the threshold of galaxy formation, enabling stringent tests of any theory that alters the properties of dark matter haloes. [ABSTRACT FROM AUTHOR]

The Hubble constant (⁠|$H_0$|⁠), a crucial parameter in cosmology, quantifies the expansion rate of the universe so its precise measurement is important to understand the fundamental dynamics of our evolving universe. One of the major limitations of measuring |$H_0$| using time-delay cosmography is the presence of the mass-sheet degeneracy (MSD) in the lens mass modelling. We propose and quantitatively assess the use of galaxy–galaxy shear measurements to break the MSD in the strong lensing mass modelling. We use stacked galaxy–galaxy lensing profiles and corresponding covariance matrices from Huang et al. to constrain the MSD in lens mass modelling with a highly flexible mass profile. Our analyses show that if ideally all galaxy–galaxy lensing measurements from the Hyper Suprime-Cam survey can be used to constrain the MSD, we can achieve |$\sim 10~{{\ \rm per\ cent}}$| precision on the MSD constraint. We forecast that galaxy–galaxy lensing measurements from Legacy Survey of Space and Time (LSST)-like surveys can in general constrain the MSD with |$\sim 1\,\mathrm{ per\,cent}-3~{{\ \rm per\ cent}}$| precision. Furthermore, if we push weak lensing measurements to a lower angular scale of |$\sim 0.04\,\rm Mpc$|⁠ , a survey like LSST can provide |$\sim 1~{{\ \rm per\ cent}}$| precision on the MSD constraint, enabling a measurement of |$H_0$| at the 1 per cent level. We demonstrate that galaxy–galaxy weak lensing can robustly constrain the MSD independent of stellar kinematics of the deflector, with wide-field survey data alone. [ABSTRACT FROM AUTHOR]

Efficient algorithms are being developed to search for strong gravitational lens systems owing to increasing large imaging surveys. Neural networks have been successfully used to discover galaxy-scale lens systems in imaging surveys such as the Kilo Degree Survey, Hyper-Suprime Cam (HSC) Survey, and Dark Energy Survey over the last few years. Thus, it has become imperative to understand how some of these networks compare, their strengths and the role of the training data sets which are essential in supervised learning algorithms used commonly in neural networks. In this work, we present the first-of-its-kind systematic comparison and benchmarking of networks from four teams that have analysed the HSC Survey data. Each team has designed their training samples and developed neural networks independently but coordinated a priori in reserving specific data sets strictly for test purposes. The test sample consists of mock lenses, real (candidate) lenses, and real non-lenses gathered from various sources to benchmark and characterize the performance of each of the network. While each team's network performed much better on their own constructed test samples compared to those from others, all networks performed comparable on the test sample with real (candidate) lenses and non-lenses. We also investigate the impact of swapping the training samples among the teams while retaining the same network architecture. We find that this resulted in improved performance for some networks. These results have direct implications on measures to be taken for lens searches with upcoming imaging surveys such as the Rubin-Legacy Survey of Space and Time, Roman, and Euclid. [ABSTRACT FROM AUTHOR]

Gravitational lensing is the phenomenon where the presence of matter (called a lens) bends the path of light-like trajectories travelling nearby. Similar to the geometric optics limit of electromagnetic waves, gravitational lensing of gravitational waves (GWs) can occur in geometric optics condition when GW wavelength is much smaller than the Schwarzschild radius of the lens, that is, |$\lambda _{\mathrm{ GW}} \ll$| R |$^{\rm s}_{\rm lens}$|⁠. This is known as the strong lensing regime for which a multiple-image system with different magnifications and phase shifts is formed. We developed GLANCE , Gravitational Lensing Authenticator using Non-modelled Cross-correlation Exploration, a novel technique to detect strongly lensed GW signals. We demonstrate that cross-correlation between two noisy reconstruction of polarized GW signals shows a non-zero value when the signals are lensed counterparts. The relative strength between the signal cross-correlation and noise cross-correlation can quantify the significance of the event(s) being lensed. Since lensing biases the inference of source parameters, primarily the luminosity distance, a joint parameter estimation of the source and lens-induced parameters is incorporated using a Bayesian framework. We applied GLANCE to synthetic strong lensing data and showed that it can detect lensed GW signals and correctly constrain the injected source and lens parameters, even when one of the signals is below match-filtered threshold signal-to-noise ratio. This demonstrates GLANCE 's capability as a robust detection technique for strongly lensed GW signals and can distinguish between lensed and unlensed events. [ABSTRACT FROM AUTHOR]

This work studies microlensing effects in strongly lensed gravitational wave (GW) signals corresponding to global minima images in galaxy-scale lenses. We find that stellar microlenses alone are unable to introduce noticeable wave effects in the global minima GW signals at strong lensing magnification |$({\mu})\lt 50$| with match value between unlensed and lensed GW signals being above |${\sim }99.5~{{\ \rm per \, cent}}$| in |${\sim }90~{{\ \rm per \, cent}}$| of systems implying that GW signals corresponding to global minima can be treated as reference signal to determine the amount of microlensing in other strongly lensed counterparts. Since the stellar microlenses introduce negligible wave effects in global minima, they can be used to probe the intermediate-mass black hole (IMBH) lenses in the galaxy lens. We show that the presence of an IMBH lens with mass in the range |$[50,10^3]~{\rm M_\odot }$| such that the global minima lies within five Einstein radius of it, the microlensing effects at |$f\lt 10^2$| Hz are mainly determined by the IMBH lens for |${\mu} \lt 50$|⁠. Assuming that a typical strong lensing magnification of 3.8 and high enough signal-to-noise ratio (in the range |${\simeq }[10, 30]$|⁠) to detect the microlensing effect in GW signals corresponding to global minima, with non-detection of IMBH-led microlensing effects in |${\simeq }15~({\simeq }150)$| lensed GW signals, we can rule out dark matter fraction |$\gt 10~{{\ \rm per \, cent}}~(\gt 1~{{\ \rm per \, cent}})$| made of IMBH population inside galaxy lenses with mass values |$\gt 150~{\rm M_\odot }$| with |${\sim }$| 90 per cent confidence. Although we have specifically used IMBHs as an example, the same analysis applies to any subhalo (or compact objects) with lensing masses (i.e. the total mass inside Einstein radius) satisfying the above criterion. [ABSTRACT FROM AUTHOR]

Cold gas evolution ties the formation of dark matter haloes to the star formation history of the universe. A primary component of cold gas, neutral atomic hydrogen (HI), can be traced by its 21-cm emission line. However, the faintness of this emission typically limits individual detections to low redshifts (⁠|$z\lesssim 0.2$|⁠). To address this limitation, we investigate the potential of targeting gravitationally lensed systems. Building on our prior galaxy–galaxy simulations, we have developed a ray-tracing code to simulate lensed HI images for known galaxies situated behind the massive hubble frontier field galaxy clusters. Our findings reveal the existence of high HI mass, high HI magnification systems in these cluster-lensing scenarios. Through simulations of hundreds of sources, we have identified compelling targets within the redshift range |$z\approx 0.7 - 1.5$|⁠. The most promising candidate from our simulations is the Great Arc at z = 0.725 in Abell 370, which should be detectable by MeerKAT in approximately 50 h. Importantly, the derived HI mass is predicted to be relatively insensitive to systematic uncertainties in the lensing model, and should be constrained within a factor of |${\sim }2.5$| for a 95 per cent confidence interval. [ABSTRACT FROM AUTHOR]

The population of strong lensing galaxies is a subset of intermediate-redshift massive galaxies, whose population-level properties are not yet well understood. In the near future, thousands of multiply imaged systems are expected to be discovered by wide-field surveys like Rubin Observatory's Legacy Survey of Space and Time and Euclid. With the soon-to-be robust population of quadruply lensed quasars, or quads, in mind, we introduce a novel technique to elucidate the empirical distribution of the galaxy population properties. Our re-imagining of the prevailing strong lensing analysis does not fit mass models to individual lenses, but instead starts with parametric models of many galaxy populations, which include generally ignored mass distribution complexities and exclude external shear for now. We construct many mock galaxy populations with different properties and obtain populations of quads from each of them. The mock 'observed' population of quads is then compared to those from the mocks using a model-free analysis based on a three-dimensional subspace of directly observable quad image properties. The distance between two quad populations in the space of image properties is measured by a metric |$\eta$|⁠ , and the distance between their parent galaxy populations in the space of galaxy properties is measured by |$\zeta$|⁠. We find a well-defined relation between |$\eta$| and |$\zeta$|⁠. The discovered relation between the space of image properties and the space of galaxy properties allows for the observed galaxy population properties to be estimated from the properties of their quads, which will be conducted in a future paper. [ABSTRACT FROM AUTHOR]

Photometric wide-area observations in the next decade will be capable of detecting a large number of galaxy-scale strong gravitational lenses, increasing the gravitational lens sample size by orders of magnitude. To aid in forecasting and analysis of these surveys, we construct a flexible model based on observed distributions for the lens and source properties and test it on the results of past lens searches, including SL2S, SuGOHI, and searches on the COSMOS HST and DES fields. We use this model to estimate the expected yields of some current and planned surveys, including Euclid Wide, Vera Rubin LSST, and Roman High Latitude Wide Area. The model proposed includes a set of free parameters to constrain on the identifiability of a lens in an image, allowing construction of prior probability distributions for different lens detection methods. The code used in this work is made publicly available. [ABSTRACT FROM AUTHOR]

Strong gravitational lensing can be used as a tool for constraining the substructure in the mass distribution of galaxies. In this study we investigate the power spectrum of dark matter perturbations in a population of 23 Hubble Space Telescope images of strong galaxy–galaxy lenses selected from The Sloan Lens ACS (SLACS) survey. We model the dark matter substructure as a Gaussian random field perturbation on a smooth lens mass potential, characterized by power-law statistics. We expand upon the previously developed machine learning framework to predict the power-law statistics by using a convolutional neural network (CNN) that accounts for both epistemic and aleatoric uncertainties. For the training sets, we use the smooth lens mass potentials and reconstructed source galaxies that have been previously modelled through traditional fits of analytical and shapelet profiles as a starting point. We train three CNNs with different training set: the first using standard data augmentation on the best-fitting reconstructed sources, the second using different reconstructed sources spaced throughout the posterior distribution, and the third using a combination of the two data sets. We apply the trained CNNs to the SLACS data and find agreement in their predictions. Our results suggest a significant substructure perturbation favouring a high frequency power spectrum across our lens population. [ABSTRACT FROM AUTHOR]

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-fitting lens light model remaining below 5 per cent. 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. [ABSTRACT FROM AUTHOR]

Dark matter structures within strong gravitational lens galaxies and along their lines of sight leave a gravitational imprint on the multiple images of lensed sources. Strong gravitational lensing provides, therefore, a key test of different dark matter models. In this article, we describe how galaxy-scale strong gravitational lensing observations are sensitive to the physical nature of dark matter. We provide an historical perspective of the field, and review its current status. We discuss the challenges and advances in terms of data, treatment of systematic errors and theoretical predictions, that will enable one to deliver a stringent and robust test of different dark matter models in the next decade. With the advent of the next generation of sky surveys, the number of known strong gravitational lens systems is expected to increase by several orders of magnitude. Coupled with high-resolution follow-up observations, these data will provide a key opportunity to constrain the properties of dark matter with strong gravitational lensing. [ABSTRACT FROM AUTHOR]

Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z ≈ 2 that we serendipitously discovered during the data reduction of the COSMOS-Web survey and that could be the most distant lens discovered to date. Methods. We extracted the visible and near-infrared photometry of the source and the lens from more than 25 bands. We combined these observations with far-infrared detections to study the dusty nature of the source and we derived the photometric redshifts and physical properties of both the lens and the source with three different spectral energy distribution (SED) fitting codes. Using JWST/NIRCam images, we also produced two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. We find the lens to be a very massive elliptical galaxy at z = 2.02 ± 0.02 with a total mass within the Einstein radius of Mtot(<θEin = (3.66 ± 0.36) × 1011M⊙ and a total stellar mass of M⋆ = 1.37−0.11+0.14 × 1011 M⊙. We also estimate it to be compact and quiescent with a specific star formation rate below 10−13 yr. Compared to stellar-to-halo mass relations from the literature, we find that the total mass of the lens within the Einstein radius is consistent with the presence of a dark matter (DM) halo of total mass Mh = 1.09−0.57+1.46 × 1013 M⊙. In addition, the background source is a M⋆ = (1.26 ± 0.17) × 1010M⊙ star-forming galaxy (SFR ≈ (78 ± 15) M⊙ yr) at z = 5.48 ± 0.06. The morphology reconstructed in the source plane shows two clear components with different colors. Dust attenuation values from SED fitting and nearby detections in the far infrared also suggest that the background source could be at least partially dust-obscured. Conclusions. We find the lens at z ≈ 2. Its total, stellar, and DM halo masses are consistent within the Einstein ring, so we do not need any unexpected changes in our description of the lens such as changing its initial mass function or including a non-negligible gas contribution. The most likely solution for the lensed source is at z ≈ 5.5. Its reconstructed morphology is complex and highly wavelength dependent, possibly because it is a merger or a main sequence galaxy with a heterogeneous dust distribution. [ABSTRACT FROM AUTHOR]

Microlensing of stars in strongly lensed galaxies can lead to temporary extreme magnification factors (μ> 1000), enabling their detection at high redshifts. Following the discovery of Icarus, several stars at cosmological distances (z > 1) have been observed using the Hubble Space Telescope (HST) and the James Webb Space Telescope (JWST). This emerging field of gravitational lensing holds promise to study individual high redshift stars. It also offers the opportunity to study the substructure in the lens plane with implications for dark matter models, as more lensed stars are detected and analysed statistically. Due to the computational demands of simulating microlensing at large magnification factors, it is important to develop fast and accurate analytical approximations for the probability of magnification in such extreme scenarios. In this study, we consider different macro-model magnification and microlensing surface mass density scenarios and study how the probability of extreme magnification factors depends on these factors. To achieve this, we created state-of-the-art large simulations of the microlensing effect in these scenarios. Through the analysis of these simulations, we derived analytical scaling relationships that can bypass the need for expensive numerical simulations. Our results are useful to interpret current observations of stars at cosmic distances which are extremely magnified and under the influence of microlenses. [ABSTRACT FROM AUTHOR]

When radio photons propagate through a non-uniform electron density volume, the plasma lensing effect can induce an extreme magnification to the observed flux at certain frequencies. Because the plasma lens acts as a diverging lens, it can extremely suppress the observed flux when aligned with the source. These two properties can theoretically cause a highly magnified fast radio burst (FRB) to become faint or even disappear for a period of time. In this paper, we interpret that the significant increase in burst counts followed by a sudden quenching in FRB 20201124A in September 2021 can be attributed to plasma lensing. Based on the one-dimensional Gaussian lens model, we search for double main-peak structures in the spectra just before its extinction on 2021 September 29. After the de-dispersion and de-scintillation procedures, we find eight bursts with double main-peaks at stable positions. There are three parameters in our modelling, the height N 0, width a of the lens and its distance D LS to the source. We reformulate them as a combined parameter |${P}_0 \propto \left(\frac{a}{\mathrm{au}}\right)\sqrt{\frac{\mathrm{kpc}}{D_{\mathrm{LS}}} \frac{\mathrm{pc}\mathrm{cm}^{-3}}{N_0} }$|⁠. The frequency spectra can give an accurate estimation of P 0 corresponding to |$\left(\frac{a}{\mathrm{au}}\right)\sqrt{\frac{\mathrm{kpc}}{D_{\mathrm{LS}}} \frac{\mathrm{pc}\mathrm{cm}^{-3}}{N_0} } \approx 28.118$|⁠ , while the time of arrival only give a relatively loose constraint on a 2/ D LS. Comparing with the observation dynamic spectra, we suggest that for a plasma lens in host galaxy, e.g. D LS ≈ 1 kpc, the width of lens can not be larger than 40 au. At last, we estimate the relative transverse motion velocity between the lens and source, |$v\approx 98\left(\frac{a}{\mathrm{au}}\right)\mathrm{km\,s^{-1}}$|⁠. [ABSTRACT FROM AUTHOR]

The standardizable nature of gravitationally lensed Type Ia supernovae (glSNe Ia) makes them an attractive target for time-delay cosmography, since a source with known luminosity breaks the mass sheet degeneracy. It is known that microlensing by stars in the lensing galaxy can add significant stochastic uncertainty to the unlensed luminosity, which is often much larger than the intrinsic scatter of the Type Ia population. In this work, we show how the temporal microlensing variations as the supernova (SN) disc expands can be used to improve the standardization of glSNe Ia. We find that SNe are standardizable if they do not cross caustics as they expand. We estimate that this will be the case for ≈6 doubly imaged systems and ≈0.3 quadruply imaged systems per year from the Vera Rubin Observatory (LSST). At the end of the 10 yr LSST survey, these systems should enable us to test for systematics in H 0 due to the mass sheet degeneracy at the |$1.00^{+0.07}_{-0.06}$| per cent level, or 1.8 ± 0.2 per cent if we can only extract time delays from the third of systems with counter-images brighter than i  = 24 mag. [ABSTRACT FROM AUTHOR]

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 along our line of sight. A small amount of external shear could come from these sources, but we show that the vast majority does not. Except in a handful of rare systems, the best-fitting values do not correlate with independent measurements of line-of-sight 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 external shear is aligned with the major or minor axis of 88 per cent of lens galaxies; and the amplitude of the external shear increases if that galaxy is discy. 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. [ABSTRACT FROM AUTHOR]

Strongly lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (SNIa) in LSST by quantifying the expected annual number of detections, the impact of stellar microlensing, follow-up feasibility, and how to best separate lensed and unlensed SNIa. We simulate SNIa lensed by galaxies, using the current LSST baseline v3.0 cadence, and find an expected number of 44 lensed SNIa detections per year. Microlensing effects by stars in the lensing galaxy are predicted to lower the lensed SNIa detections by ∼8 per cent. The lensed events can be separated from the unlensed ones by jointly considering their colours and peak magnitudes. We define a 'gold sample' of ∼10 lensed SNIa per year with time delay >10 d, >5 detections before light curve peak, and sufficiently bright (mi < 22.5 mag) for follow-up observations. In 3 yr of LSST operations, such a sample is expected to yield a 1.5 per cent measurement of the Hubble constant. [ABSTRACT FROM AUTHOR]

The flux ratios of strongly lensed quasars have previously been used to infer the properties of dark matter. In these analyses, it is crucial to separate the effect of the main lensing galaxy and the low-mass dark matter halo population. In this work, we investigate flux-ratio perturbations resulting from general third- and fourth-order multipole perturbations to the main lensing galaxy's mass profile. We simulate four lens systems, each with a different lensing configuration, without multipoles. The simulated flux ratios are perturbed by 10–40 per cent by a population of low-mass haloes consistent with cold dark matter and, in one case, also a satellite galaxy. This level of perturbation is comparable to the magnitude of flux-ratio anomalies in real data that has been previously analysed. We then attempt to fit the simulated systems using multipoles instead of low-mass haloes. We find that multipoles with amplitudes of 0.01 or less can produce flux-ratio perturbations in excess of 40 per cent. In all cases, third- or fourth-order multipoles can individually reduce the magnitude of, if not eliminate, flux-ratio anomalies. When both multipole orders are jointly included, all simulated flux ratios can be fit to within the observational uncertainty. Our results indicate that low-mass haloes and multipoles are highly degenerate when modelling quadruply imaged quasars based just on image positions and flux ratios. In the presence of this degeneracy, flux-ratio anomalies in lensed quasars alone cannot be used to place strong constraints on the properties of dark matter without additional information that can inform our priors. [ABSTRACT FROM AUTHOR]

One of the frontiers for advancing what is known about dark matter lies in using strong gravitational lenses to characterize the population of the smallest dark matter haloes. There is a large volume of information in strong gravitational lens images - the question we seek to answer is to what extent we can refine this information. To this end, we forecast the detectability of a mixed warm and cold dark matter scenario using the anomalous flux ratio method from strong gravitational lensed images. The halo mass function of the mixed dark matter scenario is suppressed relative to cold dark matter but still predicts numerous low-mass dark matter haloes relative to warm dark matter. Since the strong lensing signal receives a contribution from a range of dark matter halo masses and since the signal is sensitive to the specific configuration of dark matter haloes, not just the halo mass function, degeneracies between different forms of suppression in the halo mass function, relative to cold dark matter, can arise. We find that, with a set of lenses with different configurations of the main deflector and hence different sensitivities to different mass ranges of the halo mass function, the different forms of suppression of the halo mass function between the warm dark matter model and the mixed dark matter model can be distinguished with 40 lenses with Bayesian odds of 30:1.

Over the past year, JWST has uncovered galaxies at record-breaking distances up to z ∼ 13. The JWST UNCOVER (ultra-deep NIRSpec and NIRcam observations before the epoch of reionization) program has obtained ultra-deep multiwavelength NIRCam imaging of the massive galaxy cluster A2744 over ∼45 arcmin2 down to ∼29.5 AB mag. Here, we present a robust ultraviolet (UV) luminosity function derived through lensing clusters at 9 < z < 12. Using comprehensive end-to-end simulations, we account for all lensing effects and systematic uncertainties in deriving both the amplification factors and the effective survey volume. Our results confirm the intriguing excess of UV-bright galaxies (M UV <−20 AB mag) previously reported at z > 9 in recent JWST studies. In particular, a double power-law (DPL) describes better the bright end of the luminosity function compared to the classical Schechter form. The number density of these bright galaxies is 10–100 times larger than theoretical predictions and previous findings based on Hubble Space Telescope (HST) observations. Additionally, we measure a star formation rate density of ρSFR = 10−2.64 M⊙ yr−1 Mpc−3 at these redshifts, which is 4–10 times higher than galaxy formation models that assume a constant star formation efficiency. Future wide-area surveys and accurate modelling of lensing-assisted observations will reliably constrain both the bright and the dim end of the UV luminosity function at z > 9, which will provide key benchmarks for galaxy formation models. [ABSTRACT FROM AUTHOR]

Strong gravitational lensing is a powerful probe of the distribution of matter on sub-kpc scales. It can be used to test the existence of completely dark sub-haloes surrounding galaxies, as predicted by the standard cold dark matter model, or to test alternative dark matter models. The constraining power of the method depends strongly on photometric and astrometric precision and accuracy. We simulate and quantify the capabilities of upcoming adaptive optics systems and advanced instruments on ground-based telescopes, focusing as an illustration on the Keck Telescope (OSIRIS + KAPA, LIGER + KAPA) and the Thirty Meter Telescope (TMT; IRIS + NFIRAOS). We show that these new systems will achieve dramatic improvements over current ones in both photometric and astrometric precision. Narrow line flux ratio errors below 2 per cent, and submilliarcsecond astrometric precision will be attainable for typical quadruply imaged quasars. With TMT, the exposure times required will be of order a few minutes per system, enabling the follow-up of 100–1000 systems expected to be discovered by the Rubin, Euclid, and Roman Telescopes. [ABSTRACT FROM AUTHOR]

In this study, we investigate the impact of microlensing on gravitational wave (GW) signals in the LIGO−Virgo sensitivity band. Microlensing caused by an isolated point lens, with (redshifted) mass ranging from M Lz ∈ (1,  105) M⊙ and impact parameter y ∈ (0.01,  5), can result in a maximum mismatch of |$\sim 30~{{\ \rm per\ cent}}$| with their unlensed counterparts. When y < 1, it strongly anticorrelates with the luminosity distance enhancing the detection horizon and signal-to-noise ratio (SNR). Biases in inferred source parameters are assessed, with in-plane spin components being the most affected intrinsic parameters. The luminosity distance is often underestimated, while sky-localization and trigger times are mostly well-recovered. Study of a population of microlensed signals due to an isolated point lens primarily reveals: (i) using unlensed templates during the search causes fractional loss (20 per cent to 30 per cent) of potentially identifiable microlensed signals; (ii) the observed distribution of y challenges the notion of its high improbability at low values (y ≲ 1), especially for y ≲ 0.1; (iii) Bayes factor analysis of the population indicates that certain region in M Lz − y parameter space have a higher probability of being detected and accurately identified as microlensed. Notably, the microlens parameters for the most compelling candidate identified in previous microlensing searches, GW200208_130117, fall within a 1σ range of the aforementioned higher probability region. Identifying microlensing signatures from M Lz < 100 M⊙ remains challenging due to small microlensing effects at typical SNR values. Additionally, we also examined how microlensing from a population of microlenses influences the detection of strong lensing signatures in pairs of GW events, particularly in the posterior-overlap analysis. [ABSTRACT FROM AUTHOR]

The number of gravitational arcs systems detected is increasing quickly and should even increase at a faster rate in the near future. This wealth of new gravitational arcs requires the development of a purely automated method to reconstruct the lens and source. A general reconstruction method based on the singular perturbative approach is proposed in this paper. This method generates a lens and source reconstruction directly from the gravitational arc image. The method is fully automated and works in two steps. The first step is to generate a guess solution based on the circular solution in the singular perturbative approach. The second step is to break the sign degeneracy and to refine the solution by using a general source model. The refinement of the solution is conducted step-by-step to avoid the source-lens degeneracy issue. One important asset of this automated method is that the lens solution is written in universal terms which allows the computation of statistics. Considering the large number of lenses which should be available in the near future this ability to compute unbiased statistics is an important asset. [ABSTRACT FROM AUTHOR]

Cluster-scale strong lensing is a powerful tool for exploring the properties of dark matter and constraining cosmological models. However, due to the complex parameter space, pixelized strong lens modelling in galaxy clusters is computationally expensive, leading to the point-source approximation of strongly lensed extended images, potentially introducing systematic biases. Herein, as the first paper of the ClUsteR strong Lens modelIng for the Next-Generation observations (CURLING) program, we use lensing ray-tracing simulations to quantify the biases and uncertainties arising from the point-like image approximation for JWST -like observations. Our results indicate that the approximation works well for reconstructing the total cluster mass distribution, but can bias the magnification measurements near critical curves and the constraints on the cosmological parameters, the total matter density of the universe Ωm, and dark energy equation of state parameter w. To mitigate the biases, we propose incorporating the extended surface brightness distribution of lensed sources into the modelling. This approach reduces the bias in magnification from 46.2 per cent to 0.09 per cent for μ ∼ 1000. Furthermore, the median values of cosmological parameters align more closely with the fiducial model. In addition to the improved accuracy, we also demonstrate that the constraining power can be substantially enhanced. In conclusion, it is necessary to model cluster-scale strong lenses with pixelized multiple images, especially for estimating the intrinsic luminosity of highly magnified sources and accurate cosmography in the era of high-precision observations. [ABSTRACT FROM AUTHOR]

We report on a |$\rm {[CII]}_{158\mu \rm {m}}$| search using the Atacama Large Millimeter/submillimeter Array (ALMA) on three lensed, confirmed Ly α emitting galaxies at z ∼ 7. Our targets are ultraviolet (UV) faint systems with stellar masses on the order of M * ∼ 109 M⊙. We detect a single [C  ii ] line emission (4σ) from the brightest (L ∼ 2.4 × 1010 L ⊙) galaxy in our sample, MACS0454-1251. We determine a systemic redshift (z [C  ii ] = 6.3151 ± 0.0005) for MACS0454-1251 and measure a Ly α velocity offset of |$\Delta v \approx 300 \pm 70 \rm {km\, s}^{-1}$|⁠. The remaining two galaxies we detect no [C  ii ] but provide 3σ upper limits on their [C  ii ] line luminosities which we use to investigate the |$L_{\textrm {[CII]}} - \rm {SFR}$| relation. Overall our single [C  ii ] detection shows agreement with the relation for local dwarf galaxies. Our [C  ii ] deficient galaxies could potentially be exhibiting low metallicities (Z < Z⊙). Another possible explanation for weaker [C  ii ] emission could be strong feedback from star formation disrupting molecular clouds. We do not detect continuum emission in any of the sources, placing upper limits on their dust masses. Assuming a single dust temperature of |$T_{d}=35 \rm {K}$| dust masses (M dust) range from <4.8 × 107 M⊙ to 2.3 × 108 M⊙. Collectively, our results suggest faint reionization era sources could be metal poor and/or could have strong feedback suppressing [C  ii ] emission. [ABSTRACT FROM AUTHOR]

In recent years, there has been a remarkable development of simulation-based inference (SBI) algorithms, and they have now been applied across a wide range of astrophysical and cosmological analyses. There are a number of key advantages to these methods, centred around the ability to perform scalable statistical inference without an explicit likelihood. In this work, we propose two technical building blocks to a specific sequential SBI algorithm, truncated marginal neural ratio estimation (TMNRE). In particular, first, we develop autoregressive ratio estimation with the aim to robustly estimate correlated high-dimensional posteriors. Secondly, we propose a slice-based nested sampling algorithm to efficiently draw both posterior samples and constrained prior samples from ratio estimators, the latter being instrumental for sequential inference. To validate our implementation, we carry out inference tasks on three concrete examples: a toy model of a multidimensional Gaussian, the analysis of a stellar stream mock observation, and finally, a proof-of-concept application to substructure searches in strong gravitational lensing. In addition, we publicly release the code for both the autoregressive ratio estimator and the slice sampler. [ABSTRACT FROM AUTHOR]

We present gausSN , a Bayesian semiparametric Gaussian Process (GP) model for time-delay estimation with resolved systems of gravitationally lensed supernovae (glSNe). gausSN models the underlying light curve non-parametrically using a GP. Without assuming a template light curve for each SN type, gausSN fits for the time delays of all images using data in any number of wavelength filters simultaneously. We also introduce a novel time-varying magnification model to capture the effects of microlensing alongside time-delay estimation. In this analysis, we model the time-varying relative magnification as a sigmoid function, as well as a constant for comparison to existing time-delay estimation approaches. We demonstrate that gausSN provides robust time-delay estimates for simulations of glSNe from the Nancy Grace Roman Space Telescope and the Vera C. Rubin Observatory's Legacy Survey of Space and Time (Rubin-LSST). We find that up to 43.6 per cent of time-delay estimates from Roman and 52.9 per cent from Rubin-LSST have fractional errors of less than 5  per cent. We then apply gausSN to SN Refsdal and find the time delay for the fifth image is consistent with the original analysis, regardless of microlensing treatment. Therefore, gausSN maintains the level of precision and accuracy achieved by existing time-delay extraction methods with fewer assumptions about the underlying shape of the light curve than template-based approaches, while incorporating microlensing into the statistical error budget. gausSN is scalable for time-delay cosmography analyses given current projections of glSNe discovery rates from Rubin-LSST and Roman. [ABSTRACT FROM AUTHOR]

The flux ratios of gravitationally lensed quasars provide a powerful probe of the nature of dark matter. Importantly, these ratios are sensitive to small-scale structure, irrespective of the presence of baryons. This sensitivity may allow us to study the halo mass function even below the scales where galaxies form observable stars. For accurate measurements, it is essential that the quasar's light is emitted from a physical region of the quasar with an angular scale of milliarcseconds or larger; this minimizes microlensing effects by stars within the deflector. The warm dust region of quasars fits this criterion, as it has parsec-size physical scales and dominates the spectral energy distribution of quasars at wavelengths greater than 10 μm. The JWST Mid-Infrared Instrument is adept at detecting redshifted light in this wavelength range, offering both the spatial resolution and sensitivity required for accurate gravitational lensing flux ratio measurements. Here, we introduce our survey designed to measure the warm dust flux ratios of 31 lensed quasars. We discuss the flux-ratio measurement technique and present results for the first target, DES J0405-3308. We find that we can measure the quasar warm dust flux ratios with 3 per cent precision. Our simulations suggest that this precision makes it feasible to detect the presence of 107 M⊙ dark matter haloes at cosmological distances. Such haloes are expected to be completely dark in cold dark matter models. [ABSTRACT FROM AUTHOR]

To unlock the full spectrum of astrophysical and cosmological applications of gravitational-wave detections, it is essential to localize the associated black hole mergers to high precision inside their host galaxies. One possible method to achieve this is to compare the properties of multiple detections of gravitationally lensed binary black hole merger events with the properties of strong gravitational lens systems located in the joint sky localization of the gravitational-wave detections. In this work, we simulate the population of binary black hole mergers lensed by galaxy-scale lenses and detectable by LIGO-Virgo-Kagra in the coming decade and the population of galaxy-scale strong lenses that will be detected by Euclid. We use these simulations to investigate the prospects for localizing strongly lensed binary black hole mergers inside the lensed galaxies of 'Euclid-like' galaxy-scale strong lenses. We find that for 20– |$50\, \rm \%$| of strongly lensed gravitational-wave events the lens system is detectable with Euclid, if the event falls in its survey footprint. Of these, we expect to correctly identify the strongly lensed host galaxy as likely (with posterior probability) host galaxy – based on Bayesian evidence ranking of candidate hosts – for 34.6– |$21.9\,\mathrm{ per\,cent}$| of quadruply lensed gravitational-wave events when given an a priori 1–5 |$\deg ^{2}$| gravitational-wave-only sky localization. For triply and doubly lensed gravitational-wave events, this becomes 29.8– |$14.9\,\mathrm{ per\,cent}$| and 16.4– |$6.6\,\mathrm{ per\,cent}$| respectively. If successfully identified, however, the localization can be better than a fraction of the host-galaxy size, i.e. of order milli-arcseconds. A first detection in the coming decade, however, probably requires dedicated deep and high-resolution follow-ups and continued upgrades in the current and planned gravitational-wave detectors. [ABSTRACT FROM AUTHOR]

To better understand which sources contribute to optical passband fluxes of Q0957+561, we present and analyse light curves of the doubly imaged gravitationally lensed quasar from its discovery to 2023. After an early microlensing event, the difference light curves (describing delay-corrected flux ratios between the two quasar images) only show prominent microlensing gradients over the last 17 yr. In addition to these long time-scale gradients in the gr bands, we detect short time-scale, extrinsically induced differential variations that are highly correlated with the short time-scale intrinsic variability of the quasar in those bands. Both the accretion disc and the broad emission-line region (BELR) contribute to optical passband fluxes, and we also show that realistic contributions of the BELR account for the observed correlations in the gr bands. We would like to highlight that the BELR sources of optical passband fluxes of Q0957+561 should be taken into account when measuring accretion-disc source sizes from microlensing simulations. [ABSTRACT FROM AUTHOR]

Sagittarius A* (Sgr A*), the supermassive black hole at the centre of the Milky Way, undergoes large-amplitude near-infrared (NIR) flares that can coincide with the continuous rotation of the NIR emission region. One promising explanation for this observed NIR behaviour is a magnetic flux eruption, which occurs in three-dimensional General Relativistic Magneto-Hydrodynamic (3D GRMHD) simulations of magnetically arrested accretion flows. After running two-temperature 3D GRMHD simulations, where the electron temperature is evolved self-consistently along with the gas temperature, it is possible to calculate ray-traced images of the synchotron emission from thermal electrons in the accretion flow. Changes in the gas-dominated (σ = b 2/2ρ < 1) regions of the accretion flow during a magnetic flux eruption reproduce the NIR flaring and NIR emission region rotation of Sgr A* with durations consistent with observation. In this paper, we demonstrate that these models also predict that large (1.5x – 2x) size increases of the sub-millimeter (sub-mm) and millimeter (mm) emission region follow most NIR flares by 20–50 min. These size increases occur across a wide parameter space of black hole spin (a  = 0.3, 0.5, −0.5, and 0.9375) and initial tilt angle between the accretion flow and black hole spin axes θ0 (θ0 = 0°, 16°, and 30°). We also calculate the sub-mm polarization angle rotation and the shift of the sub-mm spectral index from zero to –0.8 during a prominent NIR flare in our high spin (a  = 0.9375) simulation. We show that, during a magnetic flux eruption, a large (∼10 r g), magnetically dominated (σ > 1), low-density, and high-temperature 'bubble' forms in the accretion flow. The drop in density inside the bubble and additional electron heating in accretion flow between 15 r g and 25 r g leads to a sub-mm size increase in corresponding images. [ABSTRACT FROM AUTHOR]

The mass distribution in massive elliptical galaxies encodes their evolutionary history, thus providing an avenue to constrain the baryonic astrophysics in their evolution. The power-law assumption for the radial mass profile in ellipticals has been sufficient to describe several observables to the noise level, including strong lensing and stellar dynamics. In this paper, we quantitatively constrained any deviation, or the lack thereof, from the power-law mass profile in massive ellipticals through joint lensing–dynamics analysis of a large statistical sample with 77 galaxy–galaxy lens systems. We performed an improved and uniform lens modelling of these systems from archival Hubble Space Telescope imaging using the automated lens modelling pipeline dolphin. We combined the lens model posteriors with the stellar dynamics to constrain the deviation from the power law after accounting for the line-of-sight lensing effects, a first for analyses on galaxy–galaxy lenses. We find that the Sloan Lens ACS Survey lens galaxies with a mean redshift of 0.2 are consistent with the power-law profile within 1.1σ (2.8σ) and the Strong Lensing Legacy Survey lens galaxies with a mean redshift of 0.6 are consistent within 0.8σ (2.1σ), for a spatially constant (Osipkov–Merritt) stellar anisotropy profile. We adopted the spatially constant anisotropy profile as our baseline choice based on previous dynamical observables of local ellipticals. However, spatially resolved stellar kinematics of lens galaxies are necessary to differentiate between the two anisotropy models. Future studies will use our lens models to constrain the mass distribution individually in the dark matter and baryonic components. [ABSTRACT FROM AUTHOR]

The arrival of the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), Euclid-Wide and Roman wide-area sensitive surveys will herald a new era in strong lens science in which the number of strong lenses known is expected to rise from |$\mathcal {O}(10^3)$| to |$\mathcal {O}(10^5)$|⁠. However, current lens-finding methods still require time-consuming follow-up visual inspection by strong lens experts to remove false positives which is only set to increase with these surveys. In this work, we demonstrate a range of methods to produce calibrated probabilities to help determine the veracity of any given lens candidate. To do this we use the classifications from citizen science and multiple neural networks for galaxies selected from the Hyper Suprime-Cam survey. Our methodology is not restricted to particular classifier types and could be applied to any strong lens classifier which produces quantitative scores. Using these calibrated probabilities, we generate an ensemble classifier, combining citizen science, and neural network lens finders. We find such an ensemble can provide improved classification over the individual classifiers. We find a false-positive rate of 10−3 can be achieved with a completeness of 46 per cent, compared to 34 per cent for the best individual classifier. Given the large number of galaxy–galaxy strong lenses anticipated in LSST, such improvement would still produce significant numbers of false positives, in which case using calibrated probabilities will be essential for population analysis of large populations of lenses and to help prioritize candidates for follow-up. [ABSTRACT FROM AUTHOR]

Spectro-polarimetric signatures of accretion disks in X-ray binaries and active galactic nuclei contain information on the masses and spins of their central black holes, as well as the geometry of matter in proximity to the compact objects. This information can be extracted by means of X-ray polarimetry. In this work, we present a fast analytical ray-tracing technique for polarized light (ARTPOL) that helps us to obtain the spinning black hole parameters from the observed properties. This technique can replace the otherwise time-consuming numerical ray-tracing calculations for any optically thick or geometrically thin accretion flow. For the purposes of illustration, we considered a standard optically thick, geometrically thin accretion disk in the equatorial plane of the Kerr black hole. We show that ARTPOL proves accurate for dimensionless spin parameter a ≤ 0.94 with a speed that is over four orders of magnitude faster than direct ray-tracing calculations. This approach opens up broader prospects for direct fittings of the spectro-polarimetric data from the Imaging X-ray Polarimetry Explorer. [ABSTRACT FROM AUTHOR]

We present the results of a ∼60-h multiband observational campaign with the Atacama Large Millimeter Array targeting a spectroscopically confirmed and lensed sub-L⋆ galaxy at z = 6.07, first identified during the ALMA Lensing Cluster Survey (ALCS). We sampled the dust continuum emission from rest frame 90–370 μm at six different frequencies and set constraining upper limits on the molecular gas line emission and content by targeting the CO (7 − 6) and [C I](3P2−3P1) transitions in two lensed images with μ ≳ 20. Complementing these submillimeter observations with deep optical and near-IR photometry and spectroscopy with JWST, we find this galaxy to form stars at a rate of SFR ∼ 7 M⊙ yr−1, ∼50 − 70% of which is obscured by dust. This is consistent with what one would predict for a M⋆ ∼ 7.5 × 108 M⊙ object by extrapolating the relation between the fraction of the obscured star formation rate and stellar mass at z <  2.5 and with observations of IR-detected objects at 5 <  z <  7. The light-weighted dust temperature of Tdust ∼ 50 K is similar to that of more massive galaxies at similar redshifts, although with large uncertainties and with possible negative gradients. We measure a dust mass of Mdust ∼ 1.5 × 106 M⊙ and, by combining [C I], [C II], and a dynamical estimate, a gas mass of Mgas ∼ 2 × 109 M⊙. Their ratio (δDGR) is in good agreement with predictions from models and empirical relations in the literature. The dust-to-stellar mass fraction of fdust ∼ 0.002 and the young stellar age (100 − 200 Myr) are consistent with efficient dust production via supernovae, as predicted by existing models and simulations of dust evolution. Also, the expected number density of galaxies with Mdust ∼ 106 M⊙ at z = 6 from a subset of these models is in agreement with the observational estimate that we set from the parent ALCS survey. The combination of gravitational lensing and deep multiwavelength observations allowed us to probe luminosity and mass regimes up to two orders of magnitude lower than what has been explored so far for field galaxies at similar redshifts. Our results serve as a benchmark for future observational endeavors of the high-redshift and faint sub-L⋆ galaxy population that might have driven the reionization of the Universe. [ABSTRACT FROM AUTHOR]

Strongly lensed systems with peculiar configurations allow us to probe the local properties of the deflecting lens mass while simultaneously testing general profile assumptions. The quasar HE0230−2130 is lensed by two galaxies at similar redshifts (Δz ∼ 0.003) into four observed images. Using modeled quasar positions from fitting the brightness of the quasar images in ground-based imaging data from the Magellan telescope, we find that lens-mass models where each of these two galaxies is parametrized with a singular power-law (PL) profile predict five quasar images. One of the predicted images is unobserved despite it being distinctively offset from the lensing galaxies and likely bright enough to be observable. This missing image gives rise to new opportunities to study the mass distribution of these galaxies. To interpret the quad configuration of the system, we tested 12 different profile assumptions with the aim of obtaining lens-mass models that correctly predict only four observed images. We tested the effects of adopting: cored profiles for the lensing galaxies; external shear; and additional profiles to represent a dark matter clump. We find that half of our model classes can produce the correct image multiplicity. By comparing the Bayesian evidence of different model parametrizations, we favor two model classes: (i) one that incorporates two singular PL profiles for the lensing galaxies and a cored isothermal sphere in the region of the previously predicted fifth image (rNIS profile), and (ii) one with a bigger lensing galaxy parametrized by a singular PL profile and the smaller galaxy by a cored PL profile with external shear. We estimated the mass of the rNIS clump for each candidate model of our final Markov chain Monte Carlo sample, and find that only 2% are in the range of 106 M⊙ ≤ MrNIS ≤ 109 M⊙, which is the predicted mass range of dark matter subhalos in cold dark matter simulations, or the mass of dark-matter-dominated and low-surface-brightness galaxies. We therefore favor the models with a cored mass distribution for the lens galaxy close to the predicted fifth image. Our study further demonstrates that lensed quasar images are sensitive to the dark matter structure in the gravitational lens. We are able to describe this exotic lensing configuration with relatively simple models, which demonstrates the power of strong lensing for studying galaxies and lens substructure. [ABSTRACT FROM AUTHOR]

Context. Strongly lensed quasars are fundamental sources for cosmology. The Gaia space mission covers the entire sky with the unprecedented resolution of 0.18″ in the optical, making it an ideal instrument to search for gravitational lenses down to the limiting magnitude of 21. Nevertheless, the previous Gaia Data Releases are known to be incomplete for small angular separations such as those expected for most lenses. Aims. We present the Data Processing and Analysis Consortium GravLens pipeline, which was built to analyse all Gaia detections around quasars and to cluster them into sources, thus producing a catalogue of secondary sources around each quasar. We analysed the resulting catalogue to produce scores that indicate source configurations that are compatible with strongly lensed quasars. Methods. GravLens uses the DBSCAN unsupervised clustering algorithm to detect sources around quasars. The resulting catalogue of multiplets is then analysed with several methods to identify potential gravitational lenses. We developed and applied an outlier scoring method, a comparison between the average BP and RP spectra of the components, and we also used an extremely randomised tree algorithm. These methods produce scores to identify the most probable configurations and to establish a list of lens candidates. Results. We analysed the environment of 3 760 032 quasars. A total of 4 760 920 sources, including the quasars, were found within 6″ of the quasar positions. This list is given in the Gaia archive. In 87% of cases, the quasar remains a single source, and in 501 385 cases neighbouring sources were detected. We propose a list of 381 lensed candidates, of which we identified 49 as the most promising ones. Beyond these candidates, the associate tables in this Focused Product Release allow the entire community to explore the unique Gaia data for strong lensing studies further. [ABSTRACT FROM AUTHOR]

We present 6 GHz Very Large Array radio images of 70 gravitational lens systems at 300 mas resolution, in which the source is an optically selected quasar, and nearly all of which have two lensed images. We find that about in half of the systems (40/70, with 33/70 secure), one or more lensed images are detected down to our detection limit of 20  μ Jy beam−1, similar to previous investigations and reinforcing the conclusion that typical optically selected quasars have intrinsic GHz radio flux densities of a few μ Jy (∼1023 W Hz−1 at redshifts of 1–2). In addition, for 10 cases it is likely that the lensing galaxies are detected in the radio. Available detections of, and limits on the far-infrared luminosities from the literature, suggest that nearly all of the sample lie on the radio-FIR correlation typical of star-forming galaxies, and that their radio luminosities are at least compatible with the radio emission being produced by star formation processes. One object, WISE2329−1258, has an extra radio component that is not present in optical images, and is difficult to explain using simple lens models. In-band spectral indices, where these can be determined, are generally moderately steep and consistent with synchrotron processes either from star formation/supernovae or AGNs. Comparison of the A/B image flux ratios at radio and optical wavelengths suggests a 10 per cent level contribution from finite source effects or optical extinction to the optical flux ratios, together with sporadic larger discrepancies that are likely to be due to optical microlensing. [ABSTRACT FROM AUTHOR]

Of the hundreds of z ≳ 6 quasars discovered to date, only one is known to be gravitationally lensed, despite the high lensing optical depth expected at z ≳ 6. High-redshift quasars are typically identified in large-scale surveys by applying strict photometric selection criteria, in particular by imposing non-detections in bands blueward of the Lyman-α line. Such procedures by design prohibit the discovery of lensed quasars, as the lensing foreground galaxy would contaminate the photometry of the quasar. We present a novel quasar selection methodology, applying contrastive learning (an unsupervised machine learning technique) to Dark Energy Survey imaging data. We describe the use of this technique to train a neural network which isolates an 'island' of 11 sources, of which seven are known z ∼ 6 quasars. Of the remaining four, three are newly discovered quasars (J0109−5424, z  = 6.07; J0122−4609, z  = 5.99; J0603−3923, z  = 5.94), as confirmed by follow-up and archival spectroscopy, implying a 91 per cent efficiency for our novel selection method; the final object on the island is a brown dwarf. In one case (J0109−5424), emission below the Lyman limit unambiguously indicates the presence of a foreground source, though high-resolution optical/near-infrared imaging is still needed to confirm the quasar's lensed (multiply imaged) nature. Detection in the g band has led this quasar to escape selection by traditional colour cuts. Our findings demonstrate that machine learning techniques can thus play a key role in unveiling populations of quasars missed by traditional methods. [ABSTRACT FROM AUTHOR]

We study the populations of stellar clumps in three high-redshift galaxies, at z = 4.92, 4.88, and 4.03, gravitationally lensed by the foreground galaxy clusters MS1358, RCS0224, and MACS0940, respectively. The lensed galaxies consist of multiple counter-images with large magnifications, mostly above |$\mu > 5$| and in some cases reaching |$\mu > 20$|⁠. We use rest-frame UV observations from the HST to extract and analyse their clump populations, counting 10, 3, and 11 unique sources, respectively. Most of the clumps have derived effective radii in the range |$R_{\rm eff}=10\!-\!100$|  pc, with the smallest one down to 6 pc, i.e. consistent with the sizes of individual stellar clusters. Their UV magnitudes correspond to |$\rm SFR_{UV}$| mostly in the range |$0.1\!-\!1\ {\rm M_\odot \, yr}^{-1}$|⁠ ; the most extreme ones, reaching |${\rm SFR_{UV}}=5\ {\rm M_\odot \, yr}^{-1}$| are among the UV-brightest compact (⁠|$R_{\rm eff} < 100$|  pc) star-forming regions observed at any redshift. Clump masses span a broad range from 106 to |$10^9\,{\rm M}_\odot$|⁠ ; stellar mass surface densities are comparable and in many cases larger than the ones of local stellar clusters, while being typically 10 times larger in size. By compiling published properties of clump populations at similar spatial resolution between redshifts 0 and 5, we find a tentative evolution of |$\Sigma_{\rm SFR}$| and |$\Sigma _{M_\star }$| with redshift, especially when very compact clumps (⁠|$R_{\rm eff}\leqslant 20$|  pc) are considered. We suggest that these trends with redshift reflect the changes in the host galaxy environments where clumps form. Comparisons with the local universe clumps/star clusters shows that, although rare, conditions for elevated clump |$\Sigma_{\rm SFR}$| and |$\Sigma _{M_\star }$| can be found. [ABSTRACT FROM AUTHOR]

While gravitational lens inversion holds great promise to reveal the structure of the light-deflecting mass distribution, both light and dark, the existence of various kinds of degeneracies implies that care must be taken when interpreting the resulting lens models. This article illustrates how thinking in terms of the projected potential helps to gain insight into these matters. Additionally it is shown explicitly how, when starting from a discretized version of the projected potential of one particular lens model, the technique of quadratic programming can be used to create a multitude of equivalent lens models that preserve all or a subset of lens properties. This method is applied to a number of scenarios, showing the lack of grasp on the mass outside the strong lensing region, revisiting mass redistribution in between images, and applying this to a recent model of the SDSS J1004+4112 cluster, as well as illustrating the generalized mass sheet degeneracy and source-position transformation. In the case of J1004, we show that this mass redistribution did not succeed at completely eliminating a dark mass clump recovered by grale near one of the quasar images. [ABSTRACT FROM AUTHOR]

Modern astronomical surveys are producing data sets of unprecedented size and richness, increasing the potential for high-impact scientific discovery. This possibility, coupled with the challenge of exploring a large number of sources, has led to the development of novel machine-learning-based anomaly detection approaches, such as astronomaly. For the first time, we test the scalability of astronomaly by applying it to almost 4 million images of galaxies from the Dark Energy Camera Legacy Survey. We use a trained deep learning algorithm to learn useful representations of the images and pass these to the anomaly detection algorithm isolation forest, coupled with astronomaly 's active learning method, to discover interesting sources. We find that data selection criteria have a significant impact on the trade-off between finding rare sources such as strong lenses and introducing artefacts into the data set. We demonstrate that active learning is required to identify the most interesting sources and reduce artefacts, while anomaly detection methods alone are insufficient. Using astronomaly , we find 1635 anomalies among the top 2000 sources in the data set after applying active learning, including eight strong gravitational lens candidates, 1609 galaxy merger candidates, and 18 previously unidentified sources exhibiting highly unusual morphology. Our results show that by leveraging the human–machine interface, astronomaly is able to rapidly identify sources of scientific interest even in large data sets. [ABSTRACT FROM AUTHOR]

The ΛCDM paradigm successfully explains the large-scale structure of the Universe, but is less well constrained on subgalactic scales. Gravitational lens modelling has been used to measure the imprints of dark substructures on lensed arcs, testing the small-scale predictions of ΛCDM. However, the methods required for these tests are subject to degeneracies among the lens mass model and the source light profile. We present a case study of the unique compound gravitational lens SDSSJ0946+1006, wherein a dark, massive substructure has been detected, whose reported high concentration would be unlikely in a ΛCDM universe. For the first time, we model the first two background sources in both I - and U -band HST imaging, as well as VLT-MUSE emission line data for the most distant source. We recover a lensing perturber at a 5.9σ confidence level with mass |$\log _{10}(M_\mathrm{sub}/{\rm M}_{\odot })=9.2^{+0.4}_{-0.1}$| and concentration |$\log _{10}c=2.4^{+0.5}_{-0.3}$|⁠. The concentration is more consistent with CDM subhaloes than previously reported, and the mass is compatible with that of a dwarf satellite galaxy whose flux is undetectable in the data at the location of the perturber. A wandering black hole with mass |$\log _{10}(M_\mathrm{BH}/{\rm M}_{\odot })=8.9^{+0.2}_{-0.1}$| is a viable alternative model. We systematically investigate alternative assumptions about the complexity of the mass distribution and source reconstruction; in all cases the subhalo is detected at around the ≥5σ level. However, the detection significance can be altered substantially (up to 11.3σ) by alternative choices for the source regularization scheme. [ABSTRACT FROM AUTHOR]

We present radio observations of 24 confirmed and candidate strongly lensed quasars identified by the Gaia Gravitational Lenses working group. We detect radio emission from eight systems in 5.5 and 9 GHz observations with the Australia Telescope Compact Array (ATCA), and 12 systems in 6 GHz observations with the Karl G. Jansky Very Large Array (VLA). The resolution of our ATCA observations is insufficient to resolve the radio emission into multiple lensed images, but we do detect multiple images from 11 VLA targets. We have analysed these systems using our observations in conjunction with existing optical measurements, including measuring offsets between the radio and optical positions for each image and building updated lens models. These observations significantly expand the existing sample of lensed radio quasars, suggest that most lensed systems are detectable at radio wavelengths with targeted observations, and demonstrate the feasibility of population studies with high-resolution radio imaging. [ABSTRACT FROM AUTHOR]

We have carried out Hubble Space Telescope (HST) snapshot observations at 1.1 μm of 281 candidate strongly lensed galaxies identified in the wide-area extragalactic surveys conducted with the Herschel Space Observatory. Our candidates comprise systems with flux densities at |$500\, \mu$| m, S 500 ≥ 80 mJy. We model and subtract the surface brightness distribution for 130 systems, where we identify a candidate for the foreground lens candidate. After combining visual inspection, archival high-resolution observations, and lens subtraction, we divide the systems into different classes according to their lensing likelihood. We confirm 65 systems to be lensed. Of these, 30 are new discoveries. We successfully perform lens modelling and source reconstruction on 23 systems, where the foreground lenses are isolated galaxies and the background sources are detected in the HST images. All the systems are successfully modelled as a singular isothermal ellipsoid. The Einstein radii of the lenses and the magnifications of the background sources are consistent with previous studies. However, the background source circularized radii (between 0.34 and 1.30 kpc) are ∼3 times smaller than the ones measured in the sub-millimetre/millimetre for a similarly selected and partially overlapping sample. We compare our lenses with those in the Sloan Lens Advanced Camera for Surveys (ACS) Survey confirming that our lens-independent selection is more effective at picking up fainter and diffuse galaxies and group lenses. This sample represents the first step towards characterizing the near-infrared properties and stellar masses of the gravitationally lensed dusty star-forming galaxies. [ABSTRACT FROM AUTHOR]

In this study, we investigate a recent finding based on strong lensing observations, which suggests that the sub-haloes observed in clusters exhibit greater compactness compared to those predicted by ΛCDM simulations. To address this discrepancy, we compare the cumulative sub-halo mass function and the M sub– V circ relation between observed clusters and 324 simulated clusters from |$\rm \small {The\,Three\,\,Hundred}$| project, focusing on the hydrodynamic resimulations using |$\rm \small {Gadget-X}$| and |$\rm \small {Gizmo-Simba}$| baryonic models. The cumulative sub-halo mass function of |$\rm \small {Gizmo-Simba}$| simulated clusters aligns with observations, while |$\rm \small {Gadget-X}$| simulations exhibit discrepancies in the lower sub-halo mass range, possibly due to its strong supernova feedback. Both |$\rm \small {Gadget-X}$| and |$\rm \small {Gizmo-Simba}$| simulations demonstrate a redshift evolution of the sub-halo mass function and the V circ function, with slightly fewer sub-haloes observed at lower redshifts. Neither the |$\rm \small {Gadget-X}$| nor |$\rm \small {Gizmo-Simba}$| (albeit a little closer) simulated clusters' predictions for the M sub– V circ relation align with the observational result. Further investigations on the correlation between sub-halo/halo properties and the discrepancy in the M sub– V circ relation reveal that the sub-halo's half mass radius and galaxy stellar age, the baryon fraction, and sub-halo distance from the cluster's centre, as well as the halo relaxation state, play important roles on reproducing this relation. Nonetheless, challenges persist in accurately reproducing the observed M sub– V circ relationship within our current hydrodynamic cluster simulation that adheres to the standard ΛCDM cosmology. These challenges may stem from shortcomings in our baryon modelling, numerical intricacies within the simulation, or even potential limitations of the ΛCDM framework. [ABSTRACT FROM AUTHOR]