Showing 1–50 of 50 results for author: Vernardos, G

Exploiting the diversity of modeling methods to probe systematic biases in strong lensing analyses

Authors: A. Galan, G. Vernardos, Q. Minor, D. Sluse, L. Van de Vyvere, M. Gomer

Abstract: In the past decade, the diversity of strong lens modeling methods has exploded, from being purely analytical to pixelated and non-parametric, or based on deep learning. We embrace this diversity by selecting different software packages and use them to blindly model independently simulated Hubble Space Telescope imaging data. To overcome the difficulties arising from using different codes and conve… ▽ More In the past decade, the diversity of strong lens modeling methods has exploded, from being purely analytical to pixelated and non-parametric, or based on deep learning. We embrace this diversity by selecting different software packages and use them to blindly model independently simulated Hubble Space Telescope imaging data. To overcome the difficulties arising from using different codes and conventions, we use the COde-independent Organized LEns STandard (COOLEST) to store, compare and release all models in a self-consistent and human-readable manner. From an ensemble of six modeling methods, we study the recovery of the lens potential parameters and properties of the reconstructed source. In particular, we simulate and infer parameters of an elliptical power-law mass distribution with external shear for the lens while each modeling method reconstructs the source differently. We find that overall, both lens and source properties are recovered reasonably well, but systematic biases arise in all methods. Interestingly, we do not observe that a single method is significantly more accurate than others, and the amount of bias largely depends on the specific lens or source property of interest. By combining posterior distributions from individual methods using equal weights, the maximal systematic biases on lens model parameters inferred from individual models are reduced by a factor of 5.4 on average. We investigate a selection of modeling effects that partly explain the observed biases, such as the cuspy nature of the background source and the accuracy of the point spread function. This work introduces, for the first time, a generic framework to compare and ease the combination of models obtained from different codes and methods, which will be key to retain accuracy in future strong lensing analyses. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 24 pages, 16 figures. Submitted to A&A

Measuring the Substructure Mass Power Spectrum of 23 SLACS Strong Galaxy-Galaxy Lenses with Convolutional Neural Networks

Authors: Joshua Fagin, Georgios Vernardos, Grigorios Tsagkatakis, Yannis Pantazis, Anowar J. Shajib, Matthew O'Dowd

Abstract: 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 per… ▽ More 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 favoring a high frequency power spectrum across our lens population. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 23 pages, 22 figures

How to Break the Mass Sheet Degeneracy with the Lightcurves of Microlensed Type Ia Supernovae

Authors: Luke Weisenbach, Thomas Collett, Ana Sainz de Murieta, Coleman Krawczyk, Georgios Vernardos, Wolfgang Enzi, Andrew Lundgren

Abstract: 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… ▽ More 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 Ia population. In this work, we show how the temporal microlensing variations as the supernova 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 $\approx$6 doubly imaged systems and $\approx$0.3 quadruply imaged systems per year in LSST. At the end of the ten year 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}$\% level, or $1.8\pm0.2$\% if we can only extract time delays from the third of systems with counter images brighter than $i=24$ mag. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 14 pages, 14 figures. As submitted to MNRAS. Comments welcome

Microlensing of strongly lensed quasars

Authors: G. Vernardos, D. Sluse, D. Pooley, R. W. Schmidt, M. Millon, L. Weisenbach, V. Motta, T. Anguita, P. Saha, M. O'Dowd, A. Peel, P. L. Schechter

Abstract: Strong gravitational lensing of quasars has the potential to unlock the poorly understood physics of these fascinating objects, as well as serve as a probe of the lensing mass distribution and of cosmological parameters. In particular, gravitational microlensing by compact bodies in the lensing galaxy can enable mapping of quasar structure to $\lt 10^{-6}$ arcsec scales. Some of this potential has… ▽ More Strong gravitational lensing of quasars has the potential to unlock the poorly understood physics of these fascinating objects, as well as serve as a probe of the lensing mass distribution and of cosmological parameters. In particular, gravitational microlensing by compact bodies in the lensing galaxy can enable mapping of quasar structure to $\lt 10^{-6}$ arcsec scales. Some of this potential has been realized over the past few decades, however the upcoming era of large sky surveys promises to bring this to full fruition. Here we review the theoretical framework of this field, describe the prominent current methods for parameter inference from quasar microlensing data across different observing modalities, and discuss the constraints so far derived on the geometry and physics of quasar inner structure. We also review the application of strong lensing and microlensing to constraining the granularity of the lens potential, i.e. the contribution of the baryonic and dark matter components, and the local mass distribution in the lens, i.e. the stellar mass function. Finally, we discuss the future of the field, including the new possibilities that will be opened by the next generation of large surveys and by new analysis methods now being developed. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: To be submitted to Space Science Reviews, Topical Collection "Strong Gravitational Lensing", eds. J. Wambsganss et al

Accelerating galaxy dynamical modeling using a neural network for joint lensing and kinematics analyses

Authors: Matthew R. Gomer, Sebastian Ertl, Luca Biggio, Han Wang, Aymeric Galan, Lyne Van de Vyvere, Dominique Sluse, Georgios Vernardos, Sherry H. Suyu

Abstract: Strong gravitational lensing is a powerful tool to provide constraints on galaxy mass distributions and cosmological parameters, such as the Hubble constant, $H_0$. Nevertheless, inference of such parameters from images of lensing systems is not trivial as parameter degeneracies can limit the precision in the measured lens mass and cosmological results. External information on the mass of the lens… ▽ More Strong gravitational lensing is a powerful tool to provide constraints on galaxy mass distributions and cosmological parameters, such as the Hubble constant, $H_0$. Nevertheless, inference of such parameters from images of lensing systems is not trivial as parameter degeneracies can limit the precision in the measured lens mass and cosmological results. External information on the mass of the lens, in the form of kinematic measurements, is needed to ensure a precise and unbiased inference. Traditionally, such kinematic information has been included in the inference after the image modeling, using spherical Jeans approximations to match the measured velocity dispersion integrated within an aperture. However, as spatially resolved kinematic measurements become available via IFU data, more sophisticated dynamical modeling is necessary. Such kinematic modeling is expensive, and constitutes a computational bottleneck which we aim to overcome with our Stellar Kinematics Neural Network (SKiNN). SKiNN emulates axisymmetric modeling using a neural network, quickly synthesizing from a given mass model a kinematic map which can be compared to the observations to evaluate a likelihood. With a joint lensing plus kinematic framework, this likelihood constrains the mass model at the same time as the imaging data. We show that SKiNN's emulation of a kinematic map is accurate to considerably better precision than can be measured (better than $1\%$ in almost all cases). Using SKiNN speeds up the likelihood evaluation by a factor of $\sim 200$. This speedup makes dynamical modeling economical, and enables lens modelers to make effective use of modern data quality in the JWST era. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: (13 pages, 9 figures, submitted to Astronomy & Astrophysics)

Euclid Preparation XXXIII. Characterization of convolutional neural networks for the identification of galaxy-galaxy strong lensing events

Authors: Euclid Collaboration, L. Leuzzi, M. Meneghetti, G. Angora, R. B. Metcalf, L. Moscardini, P. Rosati, P. Bergamini, F. Calura, B. Clément, R. Gavazzi, F. Gentile, M. Lochner, C. Grillo, G. Vernardos, N. Aghanim, A. Amara, L. Amendola, S. Andreon, N. Auricchio, S. Bardelli, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia , et al. (194 additional authors not shown)

Abstract: Forthcoming imaging surveys will potentially increase the number of known galaxy-scale strong lenses by several orders of magnitude. For this to happen, images of tens of millions of galaxies will have to be inspected to identify potential candidates. In this context, deep learning techniques are particularly suitable for the finding patterns in large data sets, and convolutional neural networks (… ▽ More Forthcoming imaging surveys will potentially increase the number of known galaxy-scale strong lenses by several orders of magnitude. For this to happen, images of tens of millions of galaxies will have to be inspected to identify potential candidates. In this context, deep learning techniques are particularly suitable for the finding patterns in large data sets, and convolutional neural networks (CNNs) in particular can efficiently process large volumes of images. We assess and compare the performance of three network architectures in the classification of strong lensing systems on the basis of their morphological characteristics. We train and test our models on different subsamples of a data set of forty thousand mock images, having characteristics similar to those expected in the wide survey planned with the ESA mission \Euclid, gradually including larger fractions of faint lenses. We also evaluate the importance of adding information about the colour difference between the lens and source galaxies by repeating the same training on single-band and multi-band images. Our models find samples of clear lenses with $\gtrsim 90\%$ precision and completeness, without significant differences in the performance of the three architectures. Nevertheless, when including lenses with fainter arcs in the training set, the three models' performance deteriorates with accuracy values of $\sim 0.87$ to $\sim 0.75$ depending on the model. Our analysis confirms the potential of the application of CNNs to the identification of galaxy-scale strong lenses. We suggest that specific training with separate classes of lenses might be needed for detecting the faint lenses since the addition of the colour information does not yield a significant improvement in the current analysis, with the accuracy ranging from $\sim 0.89$ to $\sim 0.78$ for the different models. △ Less

Submitted 26 January, 2024; v1 submitted 17 July, 2023; originally announced July 2023.

Comments: 23 pages,12 figures

Journal ref: Astronomy & Astrophysics,2024, 681, A68

Strong gravitational lensing as a probe of dark matter

Authors: S. Vegetti, S. Birrer, G. Despali, C. D. Fassnacht, D. Gilman, Y. Hezaveh, L. Perreault Levasseur, J. P. McKean, D. M. Powell, C. M. O'Riordan, G. Vernardos

Abstract: Dark matter structures within strong gravitational lens galaxies and along their line 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 a way that is independent of the baryonic content of matter structures on subgalactic scales. In this chapter, we describe how galaxy-sca… ▽ More Dark matter structures within strong gravitational lens galaxies and along their line 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 a way that is independent of the baryonic content of matter structures on subgalactic scales. In this chapter, we describe how galaxy-scale strong gravitational lensing observations are sensitive to the physical nature of dark matter. We provide a 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 near future. 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. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Strong gravitational lensing and microlensing of supernovae

Authors: Sherry H. Suyu, Ariel Goobar, Thomas Collett, Anupreeta More, Giorgos Vernardos

Abstract: Strong gravitational lensing and microlensing of supernovae (SNe) are emerging as a new probe of cosmology and astrophysics in recent years. We provide an overview of this nascent research field, starting with a summary of the first discoveries of strongly lensed SNe. We describe the use of the time delays between multiple SN images as a way to measure cosmological distances and thus constrain cos… ▽ More Strong gravitational lensing and microlensing of supernovae (SNe) are emerging as a new probe of cosmology and astrophysics in recent years. We provide an overview of this nascent research field, starting with a summary of the first discoveries of strongly lensed SNe. We describe the use of the time delays between multiple SN images as a way to measure cosmological distances and thus constrain cosmological parameters, particularly the Hubble constant, whose value is currently under heated debates. New methods for measuring the time delays in lensed SNe have been developed, and the sample of lensed SNe from the upcoming Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to provide competitive cosmological constraints. Lensed SNe are also powerful astrophysical probes. We review the usage of lensed SNe to constrain SN progenitors, acquire high-z SN spectra through lensing magnifications, infer SN sizes via microlensing, and measure properties of dust in galaxies. The current challenge in the field is the rarity and difficulty in finding lensed SNe. We describe various methods and ongoing efforts to find these spectacular explosions, forecast the properties of the expected sample of lensed SNe from upcoming surveys particularly the LSST, and summarize the observational follow-up requirements to enable the various scientific studies. We anticipate the upcoming years to be exciting with a boom in lensed SN discoveries. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 48 pages, 18 figures; to be submitted to Space Science Reviews, Topical Collection "Strong Gravitational Lensing", eds. J. Wambsganss et al

Journal ref: 2024, Space Science Reviews, Volume 220, Issue 1, article id.13 (https://link.springer.com/article/10.1007/s11214-024-01044-7)

Strong Lensing by Galaxies

Authors: A. J. Shajib, G. Vernardos, T. E. Collett, V. Motta, D. Sluse, L. L. R. Williams, P. Saha, S. Birrer, C. Spiniello, T. Treu

Abstract: Strong gravitational lensing at the galaxy scale is a valuable tool for various applications in astrophysics and cosmology. The primary uses of galaxy-scale lensing are to study elliptical galaxies' mass structure and evolution, constrain the stellar initial mass function, and measure cosmological parameters. Since the discovery of the first galaxy-scale lens in the 1980s, this field has made sign… ▽ More Strong gravitational lensing at the galaxy scale is a valuable tool for various applications in astrophysics and cosmology. The primary uses of galaxy-scale lensing are to study elliptical galaxies' mass structure and evolution, constrain the stellar initial mass function, and measure cosmological parameters. Since the discovery of the first galaxy-scale lens in the 1980s, this field has made significant advancements in data quality and modeling techniques. In this review, we describe the most common methods for modeling lensing observables, especially imaging data, as they are the most accessible and informative source of lensing observables. We then summarize the primary findings from the literature on the astrophysical and cosmological applications of galaxy-scale lenses. We also discuss the current limitations of the data and methodologies and provide an outlook on the expected improvements in both areas in the near future. △ Less

Submitted 1 May, 2024; v1 submitted 19 October, 2022; originally announced October 2022.

Comments: Submitted to Space Science Reviews, Topical Collection "Strong Gravitational Lensing", eds. J. Wambsganss et al. This version: updated after peer review

Resolving the Vicinity of Supermassive Black Holes with Gravitational Microlensing

Authors: Henry Best, Joshua Fagin, Georgios Vernardos, Matthew O'Dowd

Abstract: In the near future, wide field surveys will discover 1000's of new strongly lensed quasars, and these will be monitored with unprecedented cadence by the Legacy Survey of Space and Time (LSST). Many of these will undergo caustic-crossing microlensing events over the 10-year LSST survey, in which a sharp caustic feature from a stellar body in the lensing galaxy crosses the inner accretion disk. Cau… ▽ More In the near future, wide field surveys will discover 1000's of new strongly lensed quasars, and these will be monitored with unprecedented cadence by the Legacy Survey of Space and Time (LSST). Many of these will undergo caustic-crossing microlensing events over the 10-year LSST survey, in which a sharp caustic feature from a stellar body in the lensing galaxy crosses the inner accretion disk. Caustic-crossing events offer the unique opportunity to probe the vicinity of the central supermassive black hole for 100s of quasars with multi-platform follow-up triggered by LSST monitoring. To prepare for these observations, we have developed detailed simulations of caustic-crossing light curves. These employ a realistic analytic model of the inner accretion disk that reveals the strong surface brightness asymmetries introduced when fully accounting for both special- and general-relativistic effects. We demonstrate that an inflection in the caustic-crossing light curve due to the innermost stable circular orbit (ISCO) can be detected in reasonable follow-up observations and can be analyzed to constrain ISCO size. We also demonstrate that a convolutional neural network can be trained to predict ISCO size more reliably than traditional approaches and can also recover source orientation with high accuracy. △ Less

Submitted 11 June, 2024; v1 submitted 19 October, 2022; originally announced October 2022.

Comments: 18 pages, 17 figures

Modeling lens potentials with continuous neural fields in galaxy-scale strong lenses

Authors: Luca Biggio, Georgios Vernardos, Aymeric Galan, Austin Peel

Abstract: Strong gravitational lensing is a unique observational tool for studying the dark and luminous mass distribution both within and between galaxies. Given the presence of substructures, current strong lensing observations demand more complex mass models than smooth analytical profiles, such as power-law ellipsoids. In this work, we introduce a continuous neural field to predict the lensing potential… ▽ More Strong gravitational lensing is a unique observational tool for studying the dark and luminous mass distribution both within and between galaxies. Given the presence of substructures, current strong lensing observations demand more complex mass models than smooth analytical profiles, such as power-law ellipsoids. In this work, we introduce a continuous neural field to predict the lensing potential at any position throughout the image plane, allowing for a nearly model-independent description of the lensing mass. We apply our method on simulated Hubble Space Telescope imaging data containing different types of perturbations to a smooth mass distribution: a localized dark subhalo, a population of subhalos, and an external shear perturbation. Assuming knowledge of the source surface brightness, we use the continuous neural field to model either the perturbations alone or the full lensing potential. In both cases, the resulting model is able to fit the imaging data, and we are able to accurately recover the properties of both the smooth potential and of the perturbations. Unlike many other deep learning methods, ours explicitly retains lensing physics (i.e., the lens equation) and introduces high flexibility in the model only where required, namely, in the lens potential. Moreover, the neural network does not require pre-training on large sets of labelled data and predicts the potential from the single observed lensing image. Our model is implemented in the fully differentiable lens modeling code Herculens. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Using wavelets to capture deviations from smoothness in galaxy-scale strong lenses

Authors: Aymeric Galan, Georgios Vernardos, Austin Peel, Frédéric Courbin, Jean-Luc Starck

Abstract: Modeling the mass distribution of galaxy-scale strong gravitational lenses is a task of increasing difficulty. The high-resolution and depth of imaging data now available render simple analytical forms ineffective at capturing lens structures spanning a large range in spatial scale, mass scale, and morphology. In this work, we address the problem with a novel multiscale method based on wavelets. W… ▽ More Modeling the mass distribution of galaxy-scale strong gravitational lenses is a task of increasing difficulty. The high-resolution and depth of imaging data now available render simple analytical forms ineffective at capturing lens structures spanning a large range in spatial scale, mass scale, and morphology. In this work, we address the problem with a novel multiscale method based on wavelets. We tested our method on simulated Hubble Space Telescope (HST) imaging data of strong lenses containing the following different types of mass substructures making them deviate from smooth models: (1) a localized small dark matter subhalo, (2) a Gaussian random field (GRF) that mimics a nonlocalized population of subhalos along the line of sight, and (3) galaxy-scale multipoles that break elliptical symmetry. We show that wavelets are able to recover all of these structures accurately. This is made technically possible by using gradient-informed optimization based on automatic differentiation over thousands of parameters, which also allow us to sample the posterior distributions of all model parameters simultaneously. By construction, our method merges the two main modeling paradigms - analytical and pixelated - with machine-learning optimization techniques into a single modular framework. It is also well-suited for the fast modeling of large samples of lenses. All methods presented here are publicly available in our new Herculens package. △ Less

Submitted 21 November, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 24 pages, 12 figures, accepted for publication in A&A

Journal ref: A&A 668, A155 (2022)

The Very Knotty Lenser: exploring the role of regularization in source and potential reconstructions using Gaussian Process Regression

Authors: Georgios Vernardos, Leon V. E. Koopmans

Abstract: Reconstructing lens potentials and lensed sources can easily become an underconstrained problem, even when the degrees of freedom are low, due to degeneracies, particularly when potential perturbations superimposed on a smooth lens are included. Regularization has traditionally been used to constrain the solutions where the data failed to do so, e.g. in unlensed parts of the source. In this explor… ▽ More Reconstructing lens potentials and lensed sources can easily become an underconstrained problem, even when the degrees of freedom are low, due to degeneracies, particularly when potential perturbations superimposed on a smooth lens are included. Regularization has traditionally been used to constrain the solutions where the data failed to do so, e.g. in unlensed parts of the source. In this exploratory work, we go beyond the usual choices of regularization and adopt observationally motivated priors for the source brightness. We also perform a similar comparison when reconstructing lens potential perturbations, which are assumed to be stationary, i.e. permeate the entire field of view. We find that physically motivated priors lead to lower residuals, avoid overfitting, and are decisively preferred within a Bayesian quantitative framework in all the examples considered. For the perturbations, choosing the wrong regularization can have a detrimental effect that even high-quality data cannot correct for, while using a purely smooth lens model can absorb them to a very high degree and lead to biased solutions. Finally, our new implementation of the semi-linear inversion technique provides the first quantitative framework for measuring degeneracies between the source and the potential perturbations. △ Less

Submitted 18 February, 2022; originally announced February 2022.

Comments: 25 pages, 21 figures, 5 tables, submitted to MNRAS

Rubin-Euclid Derived Data Products: Initial Recommendations

Authors: Leanne P. Guy, Jean-Charles Cuillandre, Etienne Bachelet, Manda Banerji, Franz E. Bauer, Thomas Collett, Christopher J. Conselice, Siegfried Eggl, Annette Ferguson, Adriano Fontana, Catherine Heymans, Isobel M. Hook, Éric Aubourg, Hervé Aussel, James Bosch, Benoit Carry, Henk Hoekstra, Konrad Kuijken, Francois Lanusse, Peter Melchior, Joseph Mohr, Michele Moresco, Reiko Nakajima, Stéphane Paltani, Michael Troxel , et al. (95 additional authors not shown)

Abstract: This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum… ▽ More This report is the result of a joint discussion between the Rubin and Euclid scientific communities. The work presented in this report was focused on designing and recommending an initial set of Derived Data products (DDPs) that could realize the science goals enabled by joint processing. All interested Rubin and Euclid data rights holders were invited to contribute via an online discussion forum and a series of virtual meetings. Strong interest in enhancing science with joint DDPs emerged from across a wide range of astrophysical domains: Solar System, the Galaxy, the Local Volume, from the nearby to the primaeval Universe, and cosmology. △ Less

Submitted 13 October, 2022; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: Report of the Rubin-Euclid Derived Data Products Working Group, 78 pages, 11 figures

TDCOSMO. VII. Boxyness/discyness in lensing galaxies : Detectability and impact on $H_0$

Authors: Lyne Van de Vyvere, Matthew R. Gomer, Dominique Sluse, Dandan Xu, Simon Birrer, Aymeric Galan, Georgios Vernardos

Abstract: In the context of gravitational lensing, the density profile of lensing galaxies is often considered to be perfectly elliptical. Potential angular structures are generally ignored, except to explain flux ratios anomalies. Surprisingly, the impact of azimuthal structures on extended images of the source has not been characterized, nor its impact on the H0 inference. We address this task by creating… ▽ More In the context of gravitational lensing, the density profile of lensing galaxies is often considered to be perfectly elliptical. Potential angular structures are generally ignored, except to explain flux ratios anomalies. Surprisingly, the impact of azimuthal structures on extended images of the source has not been characterized, nor its impact on the H0 inference. We address this task by creating mock images of a point source embedded in an extended source, lensed by an elliptical galaxy on which multipolar components are added to emulate boxy/discy isodensity contours. Modeling such images with a density profile free of angular structure allow us to explore the detectability of image deformation induced by the multipoles in the residual frame. Multipole deformations are almost always detectable for our highest signal-to-noise mock data. However the detectability depends on the lens ellipticity and Einstein radius, on the S/N of the data, and on the specific lens modeling strategy. Multipoles also introduce small changes to the time delays. We therefore quantify how undetected multipoles would impact H0 inference. When no multipoles are detected in the residuals, the impact on H0 for a given lens is in general less than a few km/s/Mpc, but in the worst case scenario, combining low S/N in the ring and large intrinsic boxyness/discyness, the bias on H0 can reach 10-12 km/s/Mpc. If we now look at the inference on H0 from a population of lensing galaxies, having a distribution of multipoles representative of what is found in the light-profile of elliptical galaxies, we then find a systematic bias on H0 < 1%. The comparison of our mock systems to the state-of-the-art time delay lens sample studied by the H0LiCOW and TDCOSMO collaborations, indicates that multipoles are currently unlikely to be a source of substantial systematic bias on the inferred value of H0 from time-delay lenses. △ Less

Submitted 3 March, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: Accepted for publication in A&A; 21 pages ; content indentical to previous version, only form improvements thanks to A&A language editing ; link to the notebook associated https://github.com/TDCOSMO/TD_data_public

Journal ref: A&A 659, A127 (2022)

Strong lensing in UNIONS: Toward a pipeline from discovery to modeling

Authors: E. Savary, K. Rojas, M. Maus, B. Clément, F. Courbin, R. Gavazzi, J. H. H. Chan, C. Lemon, G. Vernardos, R. Cañameras, S. Schuldt, S. H. Suyu, J. -C. Cuillandre, S. Fabbro, S. Gwyn, M. J. Hudson, M. Kilbinger, D. Scott, C. Stone

Abstract: We present a search for galaxy-scale strong gravitational lenses in the initial 2 500 square degrees of the Canada-France Imaging Survey (CFIS). We designed a convolutional neural network (CNN) committee that we applied to a selection of 2 344 002 exquisite-seeing $r$-band images of color-selected luminous red galaxies (LRGs). Our classification uses a realistic training set where the lensing gala… ▽ More We present a search for galaxy-scale strong gravitational lenses in the initial 2 500 square degrees of the Canada-France Imaging Survey (CFIS). We designed a convolutional neural network (CNN) committee that we applied to a selection of 2 344 002 exquisite-seeing $r$-band images of color-selected luminous red galaxies (LRGs). Our classification uses a realistic training set where the lensing galaxies and the lensed sources are both taken from real data, namely the CFIS $r$-band images themselves and the Hubble Space Telescope (HST). A total of 9 460 candidates obtain a score above 0.5 with the CNN committee. After a visual inspection of the candidates, we find a total of 133 lens candidates, of which 104 are completely new. The set of false positives mainly contains ring, spiral, and merger galaxies, and to a lesser extent galaxies with nearby companions. We classify 32 of the lens candidates as secure lenses and 101 as maybe lenses. For the 32 highest quality lenses, we also fit a singular isothermal ellipsoid mass profile with external shear along with an elliptical Sersic profile for the lens and source light. This automated modeling step provides distributions of properties for both sources and lenses that have Einstein radii in the range $0.5\arcsec<θ_E<2.5\arcsec$. Finally, we introduce a new lens and/or source single-band deblending algorithm based on auto-encoder representation of our candidates. This is the first time an end-to-end lens-finding and modeling pipeline is assembled together, in view of future lens searches in a single band, as will be possible with Euclid. △ Less

Submitted 24 September, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 29 pages, 21 figures, accepted by A&A, in press

Journal ref: A&A 666, A1 (2022)

Discovery of Strongly Lensed Quasars in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS)

Authors: J. H. H. Chan, C. Lemon, F. Courbin, R. Gavazzi, B. Clément, M. Millon, E. Paic, K. Rojas, E. Savary, G. Vernardos, J. -C. Cuillandre, S. Fabbro, S. Gwyn, M. J. Hudson, M. Kilbinger, A. McConnachie

Abstract: We report the discovery of five new doubly-imaged lensed quasars from the first 2500 square degrees of the ongoing Canada-France Imaging Survey (CFIS), which is a component of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), selected from initial catalogues of either Gaia pairs or MILLIQUAS quasars. We take advantage of the deep, 0.6'' median-seeing $r$-band imaging of CFIS to confi… ▽ More We report the discovery of five new doubly-imaged lensed quasars from the first 2500 square degrees of the ongoing Canada-France Imaging Survey (CFIS), which is a component of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), selected from initial catalogues of either Gaia pairs or MILLIQUAS quasars. We take advantage of the deep, 0.6'' median-seeing $r$-band imaging of CFIS to confirm the presence of multiple point sources with similar colour of $u-r$, via convolution of the Laplacian of the point spread function. Requiring similar-colour point sources with flux ratios less than 2.5 mag in $r$-band, reduces the number of candidates from 256314 to 7815. After visual inspection we obtain 30 high-grade candidates, and prioritise spectroscopic follow-up for those showing signs of a lensing galaxy upon subtraction of the point sources. We obtain long-slit spectra for 18 candidates with ALFOSC on the 2.56-m Nordic Optical Telescope (NOT), confirming five new doubly lensed quasars with $1.21<z<3.36$ and angular separations from 0.8'' to 2.5''. One additional system is a probable lensed quasar based on the CFIS imaging and existing SDSS spectrum. We further classify six objects as nearly identical quasars -- still possible lenses but without the detection of a lensing galaxy. Given our recovery rate ($83\%$) of existing optically bright lenses within the CFIS footprint, we expect that a similar strategy, coupled with $u-r$ colour-selection from CFIS alone, will provide an efficient and complete discovery of small-separation lensed quasars of source redshifts below $z=2.7$ within the CFIS $r$-band magnitude limit of 24.1 mag. △ Less

Submitted 1 December, 2021; v1 submitted 18 October, 2021; originally announced October 2021.

Comments: 15 pages, 8 figures

Journal ref: A&A 659, A140 (2022)

SEAGLE--III: Towards resolving the mismatch in the dark-matter fraction in early-type galaxies between simulations and observations

Authors: Sampath Mukherjee, Léon V. E. Koopmans, Crescenzo Tortora, Matthieu Schaller, R. Benton Metcalf, Joop Schaye, Georgios Vernardos

Abstract: The central dark-matter fraction of galaxies is sensitive to feedback processes during galaxy formation. Strong gravitational lensing has been effective in the precise measurement of the dark-matter fraction inside massive early-type galaxies. Here, we compare the projected dark-matter fraction of early-type galaxies inferred from the SLACS strong-lens survey, with those obtained from the EAGLE, I… ▽ More The central dark-matter fraction of galaxies is sensitive to feedback processes during galaxy formation. Strong gravitational lensing has been effective in the precise measurement of the dark-matter fraction inside massive early-type galaxies. Here, we compare the projected dark-matter fraction of early-type galaxies inferred from the SLACS strong-lens survey, with those obtained from the EAGLE, Illustris, and IllustrisTNG hydro-dynamical simulations. Previous comparisons with some simulations revealed a large discrepancy, with considerably higher inferred dark-matter fractions -- by factors 2-3 -- inside half of the effective radius in observed strong-lens galaxies as compared to simulated galaxies. Here, we report good agreement between EAGLE and SLACS for the dark-matter fractions inside both half of the effective radius and the effective radius as a function of the galaxy's stellar mass, effective radius, and total mass-density slope. However, for IllustrisTNG and Illustris, the dark-matter fractions are lower than observed. This work consistently assumes a Chabrier IMF, which suggests that a different IMF (although not excluded) is not necessary to resolve this mismatch. The differences in the stellar feedback model between EAGLE and Illustris and IllustrisTNG, are likely the dominant cause of the difference in their dark-matter fraction, and density slope. △ Less

Submitted 12 November, 2021; v1 submitted 14 October, 2021; originally announced October 2021.

Comments: Matches with the published version. Fixed an error in Fig 2. Updated reference. Conclusions unchanged. 7 pages, 4 figures

Constraining quasar structure using high-frequency microlensing variations and continuum reverberation

Authors: E. Paic, G. Vernardos, D. Sluse, M. Millon, F. Courbin, J. H. Chan, V. Bonvin

Abstract: Gravitational microlensing is a powerful tool to probe the inner structure of strongly lensed quasars and to constrain parameters of the stellar mass function of lens galaxies. This is done by analysing microlensing light curves between the multiple images of strongly lensed quasars, under the influence of three main variable components: 1- the continuum flux of the source, 2- microlensing by star… ▽ More Gravitational microlensing is a powerful tool to probe the inner structure of strongly lensed quasars and to constrain parameters of the stellar mass function of lens galaxies. This is done by analysing microlensing light curves between the multiple images of strongly lensed quasars, under the influence of three main variable components: 1- the continuum flux of the source, 2- microlensing by stars in the lens galaxy and 3- reverberation of the continuum by the Broad Line Region (BLR). The latter, ignored by state-of-the-art microlensing techniques, can introduce high-frequency variations which we show carry information on the BLR size. We present a new method which includes all these components simultaneously and fits the power spectrum of the data in the Fourier space, rather than the observed light curve itself. In this new framework, we analyse COSMOGRAIL light curves of the two-image system QJ0158-4325 known to display high-frequency variations. Using exclusively the low frequency part of the power spectrum our constraint on the accretion disk radius agrees with the thin disk model estimate and previous work that fit the microlensing light curves in real space. However, if we also take into account the high-frequency variations, the data favour significantly smaller disk sizes than previous microlensing measurements. In this case, our results are in agreement with the thin disk model prediction only if we assume very low mean masses for the microlens population, i.e. <M> = 0.01 $M_\odot$. Eventually, including the differentially microlensed continuum reverberation by the BLR successfully explains the high frequencies without requiring such low mass microlenses. This allows us to measure, for the first time, the size of the BLR using single-band photometric monitoring, $R_{BLR}$ = $1.6^{+1.5}_{-0.8}\times 10^{17}$cm, in agreement with estimates using the BLR size-luminosity relation. △ Less

Submitted 11 October, 2021; originally announced October 2021.

Journal ref: A&A 659, A21 (2022)

High-quality strong lens candidates in the final Kilo Degree survey footprint

Authors: R. Li, N. R. Napolitano, C. Spiniello, C. Tortora, K. Kuijken, L. V. E. Koopmans, P. Schneider, F. Getman, L. Xie, L. Long, W. Shu, G. Vernardos, Z. Huang, G. Covone, A. Dvornik, C. Heymans, H. Hildebrandt, M. Radovich, A. H. Wright

Abstract: We present 97 new high-quality strong lensing candidates found in the final $\sim 350\,\rm deg^2$, that completed the full $\sim 1350\,\rm deg^2$ area of the Kilo-Degree Survey (KiDS). Together with our previous findings, the final list of high-quality candidates from KiDS sums up to 268 systems. The new sample is assembled using a new Convolutional Neural Network (CNN) classifier applied to $r$-b… ▽ More We present 97 new high-quality strong lensing candidates found in the final $\sim 350\,\rm deg^2$, that completed the full $\sim 1350\,\rm deg^2$ area of the Kilo-Degree Survey (KiDS). Together with our previous findings, the final list of high-quality candidates from KiDS sums up to 268 systems. The new sample is assembled using a new Convolutional Neural Network (CNN) classifier applied to $r$-band (best seeing) and $g,~r,~i$ color-composited images separately. This optimizes the complementarity of the morphology and color information on the identification of strong lensing candidates. We apply the new classifiers to a sample of luminous red galaxies (LRGs) and a sample of bright galaxies (BGs) and select candidates that received a high probability to be a lens from the CNN ($P_{\rm CNN}$). In particular, setting $P_{\rm CNN}>0.8$ for the LRGs, the $1$-band CNN predicts 1213 candidates, while the $3$-band classifier yields 1299 candidates, with only $\sim$30\% overlap. For the BGs, in order to minimize the false positives, we adopt a more conservative threshold, $P_{\rm CNN} >0.9$, for both CNN classifiers. This results in 3740 newly selected objects. The candidates from the two samples are visually inspected by 7 co-authors to finally select 97 "high-quality" lens candidates which received mean scores larger than 6 (on a scale from 0 to 10). We finally discuss the effect of the seeing on the accuracy of CNN classification and possible avenues to increase the efficiency of multi-band classifiers, in preparation of next-generation surveys from ground and space. △ Less

Submitted 30 December, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: Published by APJ

Time delay estimation in unresolved lensed quasars

Authors: L. Biggio, A. Domi, S. Tosi, G. Vernardos, D. Ricci, L. Paganin, G. Bracco

Abstract: Time-delay cosmography can be used to infer the Hubble parameter $H_0$ by measuring the relative time delays between multiple images of gravitationally-lensed quasars. A few of such systems have already been used to measure $H_0$: their time delays were determined from the multiple images light curves obtained by regular, years long, monitoring campaigns. Such campaigns can hardly be performed by… ▽ More Time-delay cosmography can be used to infer the Hubble parameter $H_0$ by measuring the relative time delays between multiple images of gravitationally-lensed quasars. A few of such systems have already been used to measure $H_0$: their time delays were determined from the multiple images light curves obtained by regular, years long, monitoring campaigns. Such campaigns can hardly be performed by any telescope: many facilities are often over-subscribed with a large amount of observational requests to fulfill. While the ideal systems for time-delay measurements are lensed quasars whose images are well resolved by the instruments, several lensed quasars have a small angular separation between the multiple images, and would appear as a single, unresolved, image to a large number of telescopes featuring poor angular resolutions or located in not privileged geographical locations. Methods allowing to infer the time delay also from unresolved light curves would boost the potential of such telescopes and greatly increase the available statistics for $H_0$ measurements. This work presents a study of unresolved lensed quasar systems to estimate the time delay using a deep learning-based approach that exploits the capabilities of one-dimensional convolutional neural networks. Experiments on state-of-the-art simulations of unresolved light curves show the potential of the proposed method and pave the way for future applications in time-delay cosmography. △ Less

Submitted 3 February, 2022; v1 submitted 3 October, 2021; originally announced October 2021.

Strong lens systems search in the Dark Energy Survey using Convolutional Neural Networks

Authors: K. Rojas, E. Savary, B. Clément, M. Maus, F. Courbin, C. Lemon, J. H. H. Chan, G. Vernardos, R. Joseph, R. Cañameras, A. Galan

Abstract: We performed a search for strong lens galaxy-scale systems in the first data release of the Dark Energy Survey (DES), from a color-selected parent sample of 18~745~029 Luminous Red Galaxies (LRGs). Our search was based on a Convolutional Neural Network (CNN) to grade our LRG selection with values between 0 (non-lens) and 1 (lens). Our training set was data-driven, i.e. using lensed sources taken f… ▽ More We performed a search for strong lens galaxy-scale systems in the first data release of the Dark Energy Survey (DES), from a color-selected parent sample of 18~745~029 Luminous Red Galaxies (LRGs). Our search was based on a Convolutional Neural Network (CNN) to grade our LRG selection with values between 0 (non-lens) and 1 (lens). Our training set was data-driven, i.e. using lensed sources taken from HST COSMOS images and where the light distribution of the lens plane was taken directly from DES images of our LRGs. A total of 76~582 cutouts obtained a score above 0.9. These were visually inspected and resulted in two catalogs. The first one contains 405 lens candidates, where 90 present clear lensing features and counterparts, while the others 315 require more evidence, such as higher resolution images or spectra to be conclusive. A total of 186 candidates were totally new identified in this search. The second catalog includes 539 ring galaxy candidates that will be useful to train CNNs against this type of false positives. For the 90 best lens candidates we carried out color-based deblending of the lens and source light without fitting any analytical profile to the data. The method turned out to be very efficient in the deblending, even for very compact objects and for objects with very complex morphology. Finally, from the 90 best lens candidates we selected 52 systems having one single deflector, to test an automated modeling pipeline which successfully modeled 79\% of the sample within an acceptable amount of computing time. △ Less

Submitted 31 August, 2021; originally announced September 2021.

Comments: 37 pages, 31 figures, submitted to A&A

Journal ref: A&A 668, A73 (2022)

Simulating time-varying strong lenses

Authors: Georgios Vernardos

Abstract: We present a self-consistent and versatile forward modelling software package that can produce time series and pixel-level simulations of time-varying strongly lensed systems. The time dimension, which needs to take into account different physical mechanisms for variability such as microlensing, has been missing from existing approaches and it is of direct relevance to time delay, and consequently… ▽ More We present a self-consistent and versatile forward modelling software package that can produce time series and pixel-level simulations of time-varying strongly lensed systems. The time dimension, which needs to take into account different physical mechanisms for variability such as microlensing, has been missing from existing approaches and it is of direct relevance to time delay, and consequently H0, measurements and caustic crossing event predictions. Such experiments are becoming more streamlined, especially with the advent of time domain surveys, and understanding their systematic and statistical uncertainties in a model-aware and physics-driven way can help improve their accuracy and precision. Here we demonstrate the software's capabilities by exploring the effect of measuring time delays from lensed quasars and supernovae in many wavelengths and under different microlensing and intrinsic variability assumptions. In this initial application, we find that the cadence of the observations and combining information from different wavelengths plays an important role in the correct recovery of the time delays. The Mock Lenses in Time (MOLET) software package is available at: \url{https://github.com/gvernard/molet} △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: 14 pages, 8 figures, submitted to MNRAS

Discovery of two Einstein crosses from massive post--blue nugget galaxies at z>1 in KiDS

Authors: N. R. Napolitano, R. Li, C. Spiniello, C. Tortora, A. Sergeyev, G. D'Ago, X. Guo, L. Xie, M. Radovich, N. Roy, L. V. E. Koopmans, K. Kuijken, M. Bilicki, T. Erben, F. Getman, C. Heymans, H. Hildebrandt, C. Moya, H. Y. Shan, G. Vernardos, A. H. Wright

Abstract: We report the discovery of two Einstein Crosses (ECs) in the footprint of the Kilo-Degree Survey (KiDS): KIDS J232940-340922 and KIDS J122456+005048. Using integral field spectroscopy from MUSE@VLT, we confirm their gravitational-lens nature. In both cases, the four spectra of the source clearly show a prominence of absorption features, hence revealing an evolved stellar population with little sta… ▽ More We report the discovery of two Einstein Crosses (ECs) in the footprint of the Kilo-Degree Survey (KiDS): KIDS J232940-340922 and KIDS J122456+005048. Using integral field spectroscopy from MUSE@VLT, we confirm their gravitational-lens nature. In both cases, the four spectra of the source clearly show a prominence of absorption features, hence revealing an evolved stellar population with little star formation. The lensing model of the two systems, assuming a singular isothermal ellipsoid (SIE) with external shear, shows that: 1) the two crosses, located at redshift $z=0.38$ and 0.24, have Einstein radius $R_{\rm E}=5.2$ kpc and 5.4 kpc, respectively; 2) their projected dark matter fractions inside the half effective radius are 0.60 and 0.56 (Chabrier IMF); 3) the sources are ultra-compact galaxies, $R_{\rm e}\sim0.9$ kpc (at redshift $z_{\rm s}=1.59$) and $R_{\rm e}\sim0.5$ kpc ($z_{\rm s}=1.10$), respectively. These results are unaffected by the underlying mass density assumption. Due to size, blue color and absorption-dominated spectra, corroborated by low specific star-formation rates derived from optical-NIR spectral energy distribution fitting, we argue that the two lensed sources in these ECs are blue nuggets migrating toward their quenching phase. △ Less

Submitted 19 November, 2020; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted for publication on APJL

Projected Cosmological Constraints from Strongly Lensed Supernovae with the Roman Space Telescope

Authors: J. D. R. Pierel, S. Rodney, G. Vernardos, M. Oguri, R. Kessler, T. Anguita

Abstract: One of the primary mission objectives for the Roman Space Telescope is to investigate the nature of dark energy with a variety of methods. Observations of Type Ia supernovae (SNIa) will be one of the principal anchors of the Roman cosmology program, through traditional luminosity distance measurements. This SNIa cosmology program can provide another valuable cosmological probe, without altering th… ▽ More One of the primary mission objectives for the Roman Space Telescope is to investigate the nature of dark energy with a variety of methods. Observations of Type Ia supernovae (SNIa) will be one of the principal anchors of the Roman cosmology program, through traditional luminosity distance measurements. This SNIa cosmology program can provide another valuable cosmological probe, without altering the mission strategy: time delay cosmography with gravitationally lensed SN. In this work, we forecast lensed SN cosmology constraints with the Roman Space Telescope, while providing useful tools for future work. Using anticipated characteristics of the Roman SNIa survey, we have constructed mock catalogs of expected resolved lensing systems, as well as strongly lensed Type Ia and core-collapse (CC) SN light curves, including microlensing effects. We predict Roman will find ~11 lensed SNIa and ~20 CCSN, dependent on the survey strategy. Next, we estimate the time delay precision obtainable with Roman (Ia: ~2 days, CC: ~3 days), and use a Fisher matrix analysis to derive projected constraints on $H_0$,$Ω_m$, and the dark energy Equation of State (EOS), $w$, for each SNIa survey strategy. A strategy optimized for high-redshift SNIa discovery is preferred when considering the constraints possible from both SNIa and lensed SN cosmology, also delivering ~1.5 times more lensed SN than other proposed survey strategies. △ Less

Submitted 24 February, 2021; v1 submitted 23 October, 2020; originally announced October 2020.

Comments: Accepted Version (ApJ, 1/4/21)

Quantifying the structure of strong gravitational lens potentials with uncertainty-aware deep neural networks

Authors: Georgios Vernardos, Grigorios Tsagkatakis, Yannis Pantazis

Abstract: Gravitational lensing is a powerful tool for constraining substructure in the mass distribution of galaxies, be it from the presence of dark matter sub-halos or due to physical mechanisms affecting the baryons throughout galaxy evolution. Such substructure is hard to model and is either ignored by traditional, smooth modelling, approaches, or treated as well-localized massive perturbers. In this w… ▽ More Gravitational lensing is a powerful tool for constraining substructure in the mass distribution of galaxies, be it from the presence of dark matter sub-halos or due to physical mechanisms affecting the baryons throughout galaxy evolution. Such substructure is hard to model and is either ignored by traditional, smooth modelling, approaches, or treated as well-localized massive perturbers. In this work, we propose a deep learning approach to quantify the statistical properties of such perturbations directly from images, where only the extended lensed source features within a mask are considered, without the need of any lens modelling. Our training data consist of mock lensed images assuming perturbing Gaussian Random Fields permeating the smooth overall lens potential, and, for the first time, using images of real galaxies as the lensed source. We employ a novel deep neural network that can handle arbitrary uncertainty intervals associated with the training dataset labels as input, provides probability distributions as output, and adopts a composite loss function. The method succeeds not only in accurately estimating the actual parameter values, but also reduces the predicted confidence intervals by 10 per cent in an unsupervised manner, i.e., without having access to the actual ground truth values. Our results are invariant to the inherent degeneracy between mass perturbations in the lens and complex brightness profiles for the source. Hence, we can quantitatively and robustly quantify the smoothness of the mass density of thousands of lenses, including confidence intervals, and provide a consistent ranking for follow-up science. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Comments: 13 pages, 10 figures, accepted for publication in MNRAS

A Quasar Microlensing Light Curve Generator for LSST

Authors: Favio Neira, Timo Anguita, Georgios Vernardos

Abstract: We present a tool to generate mock quasar microlensing light curves and sample them according to any observing strategy. An updated treatment of the fixed and random velocity components of observer, lens, and source is used, together with a proper alignment with the external shear defining the magnification map caustic orientation. Our tool produces quantitative results on high magnification event… ▽ More We present a tool to generate mock quasar microlensing light curves and sample them according to any observing strategy. An updated treatment of the fixed and random velocity components of observer, lens, and source is used, together with a proper alignment with the external shear defining the magnification map caustic orientation. Our tool produces quantitative results on high magnification events and caustic crossings, which we use to study three lensed quasars known to display microlensing, viz. RX J1131-1231, HE 0230-2130, and Q 2237+0305, as they would be monitored by The Rubin Observatory Legacy Survey of Space and Time (LSST). We conclude that depending on the location on the sky, the lens and source redshift, and the caustic network density, the microlensing variability may deviate significantly than the expected $\sim$20-year average time scale (Mosquera & Kochanek 2011, arXiv:1104.2356). We estimate that $\sim300$ high magnification events with $Δ$mag$>1$ mag could potentially be observed by LSST each year. The duration of the majority of high magnification events is between 10 and 100 days, requiring a very high cadence to capture and resolve them. Uniform LSST observing strategies perform the best in recovering microlensing high magnification events. Our web tool can be extended to any instrument and observing strategy, and is freely available as a service at http://gerlumph.swin.edu.au/tools/lsst_generator/, along with all the related code. △ Less

Submitted 29 October, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: 10 pages, 6 figures, 2 tables. Published in MNRAS. Updated Table 1

New high-quality strong lens candidates with deep learning in the Kilo Degree Survey

Authors: R. Li, N. R. Napolitano, C. Tortora, C. Spiniello, L. V. E. Koopmans, Z. Huang, N. Roy, G. Vernardos, S. Chatterjee, B. Giblin, F. Getman, M. Radovich, G. Covone, K. Kuijken

Abstract: We report new high-quality galaxy scale strong lens candidates found in the Kilo Degree Survey data release 4 using Machine Learning. We have developed a new Convolutional Neural Network (CNN) classifier to search for gravitational arcs, following the prescription by \cite{2019MNRAS.484.3879P} and using only $r-$band images. We have applied the CNN to two "predictive samples": a Luminous red galax… ▽ More We report new high-quality galaxy scale strong lens candidates found in the Kilo Degree Survey data release 4 using Machine Learning. We have developed a new Convolutional Neural Network (CNN) classifier to search for gravitational arcs, following the prescription by \cite{2019MNRAS.484.3879P} and using only $r-$band images. We have applied the CNN to two "predictive samples": a Luminous red galaxy (LRG) and a "bright galaxy" (BG) sample ($r<21$). We have found 286 new high probability candidates, 133 from the LRG sample and 153 from the BG sample. We have then ranked these candidates based on a value that combines the CNN likelihood to be a lens and the human score resulting from visual inspection (P-value) and we present here the highest 82 ranked candidates with P-values $\ge 0.5$. All these high-quality candidates have obvious arc or point-like features around the central red defector. Moreover, we define the best 26 objects, all with scores P-values $\ge 0.7$ as a "golden sample" of candidates. This sample is expected to contain very few false positives and thus it is suitable for follow-up observations. The new lens candidates come partially from the the more extended footprint adopted here with respect to the previous analyses, partially from a larger predictive sample (also including the BG sample). These results show that machine learning tools are very promising to find strong lenses in large surveys and more candidates that can be found by enlarging the predictive samples beyond the standard assumption of LRGs. In the future, we plan to apply our CNN to the data from next-generation surveys such as the Large Synoptic Survey Telescope, Euclid, and the Chinese Space Station Optical Survey. △ Less

Submitted 17 June, 2020; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: Accepted by APJ

The Most Powerful Lenses in the Universe: Quasar Microlensing as a Probe of the Lensing Galaxy

Authors: David Pooley, Timo Anguita, Saloni Bhatiani, George Chartas, Matthew Cornachione, Xinyu Dai, Carina Fian, Evencio Mediavilla, Christopher Morgan, Verónica Motta, Leonidas A. Moustakas, Sampath Mukherjee, Matthew J. O'Dowd, Karina Rojas, Dominique Sluse, Georgios Vernardos, Rachel Webster

Abstract: Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will… ▽ More Optical and X-ray observations of strongly gravitationally lensed quasars (especially when four separate images of the quasar are produced) determine not only the amount of matter in the lensing galaxy but also how much is in a smooth component and how much is composed of compact masses (e.g., stars, stellar remnants, primordial black holes, CDM sub-halos, and planets). Future optical surveys will discover hundreds to thousands of quadruply lensed quasars, and sensitive X-ray observations will unambiguously determine the ratio of smooth to clumpy matter at specific locations in the lensing galaxies and calibrate the stellar mass fundamental plane, providing a determination of the stellar $M/L$. A modest observing program with a sensitive, sub-arcsecond X-ray imager, combined with the planned optical observations, can make those determinations for a large number (hundreds) of the lensing galaxies, which will span a redshift range of $\sim$$0.25<z<1.5$ △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: Astro2020 Science White Paper

Astro2020 Science White Paper - Quasar Microlensing: Revolutionizing our Understanding of Quasar Structure and Dynamics

Authors: Leonidas Moustakas, Matthew O'Dowd, Timo Anguita, Rachel Webster, George Chartas, Matthew Cornachione, Xinyu Dai, Carina Fian, Damien Hutsemekers, Jorge Jimenez-Vicente, Kathleen Labrie, Geraint Lewis, Chelsea Macleod, Evencio Mediavilla, Christopher W Morgan, Veronica Motta, Anna Nierenberg, David Pooley, Karina Rojas, Dominique Sluse, Georgios Vernardos, Joachim Wambsganss, Suk Yee Yong

Abstract: Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide infor… ▽ More Microlensing by stars within distant galaxies acting as strong gravitational lenses of multiply-imaged quasars, provides a unique and direct measurement of the internal structure of the lensed quasar on nano-arcsecond scales. The measurement relies on the temporal variation of high-magnification caustic crossings which vary on timescales of days to years. Multiwavelength observations provide information from distinct emission regions in the quasar. Through monitoring of these strong gravitational lenses, a full tomographic view can emerge with Astronomical-Unit scale resolution. Work to date has demonstrated the potential of this technique in about a dozen systems. In the 2020s there will be orders of magnitude more systems to work with. Monitoring of lens systems for caustic-crossing events to enable triggering of multi-platform, multi-wavelength observations in the 2020s will fulfill the potential of quasar microlensing as a unique and comprehensive probe of active black hole structure and dynamics. △ Less

Submitted 29 April, 2019; originally announced April 2019.

Comments: White paper submitted to Astro2020 decadal survey; 7 pages, 3 figures

Quasar microlensing light curve analysis using deep machine learning

Authors: Georgios Vernardos, Grigorios Tsagkatakis

Abstract: We introduce a deep machine learning approach to studying quasar microlensing light curves for the first time by analyzing hundreds of thousands of simulated light curves with respect to the accretion disc size and temperature profile. Our results indicate that it is possible to successfully classify very large numbers of diverse light curve data and measure the accretion disc structure. The detai… ▽ More We introduce a deep machine learning approach to studying quasar microlensing light curves for the first time by analyzing hundreds of thousands of simulated light curves with respect to the accretion disc size and temperature profile. Our results indicate that it is possible to successfully classify very large numbers of diverse light curve data and measure the accretion disc structure. The detailed shape of the accretion disc brightness profile is found to play a negligible role, in agreement with Mortonson et al. (2005). The speed and efficiency of our deep machine learning approach is ideal for quantifying physical properties in a `big-data' problem setup. This proposed approach looks promising for analyzing decade-long light curves for thousands of microlensed quasars, expected to be provided by the Large Synoptic Survey Telescope. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Spectroscopic confirmation and modelling of two lensed quadruple quasars in the Dark Energy Survey public footprint

Authors: C. Spiniello, A. V. Sergeyev, L. Marchetti, C. Tortora, N. R. Napolitano, V. Shalyapin, A. Agnello, F. I. Getman, M. Vaccari, S. Serjeant, L. V. E. Koopmans, A. J. Baker, T. H. Jarrett, G. Covone, G. Vernardos

Abstract: Quadruply lensed quasars are extremely rare objects, but incredibly powerful cosmological tools. Only few dozen are known in the whole sky. Here we present the spectroscopic confirmation of two new quadruplets WG0214-2105 and WG2100-4452 discovered by Agnello & Spiniello (2018) within the Dark Energy Survey (DES) public footprints. We have conducted spectroscopic follow-up of these systems with th… ▽ More Quadruply lensed quasars are extremely rare objects, but incredibly powerful cosmological tools. Only few dozen are known in the whole sky. Here we present the spectroscopic confirmation of two new quadruplets WG0214-2105 and WG2100-4452 discovered by Agnello & Spiniello (2018) within the Dark Energy Survey (DES) public footprints. We have conducted spectroscopic follow-up of these systems with the Southern African Large Telescope as part of a program that aims at confirming the largest possible number of optically selected strong gravitational lensing systems in the Equatorial and Southern Hemisphere. For both systems, we present the spectra for the sources and deflectors that allowed us to estimate the source redshifts and unambiguously confirm their lensing nature. For the brighter deflector (WG2100-4452), we measure the stellar velocity dispersion from the spectrum. We also obtain photometry for both lenses, directly from DES multi-band images, isolating the lens galaxies from the quasar images. One of the quadruplets, WG0214-2105, was also observed by Pan-STARRS, allowing us to estimate the apparent brightness of each quasar image at two different epochs, and thus to find evidence for flux variability. This result could suggest a microlensing event for the faintest components, although intrinsic variability cannot be excluded with only two epochs. Finally, we present simple lens models for both quadruplets, obtaining Einstein radii, SIE velocity dispersions, ellipticities, and position angles of the lens systems, as well as time delay predictions assuming a concordance cosmological model. △ Less

Submitted 16 January, 2019; v1 submitted 14 January, 2019; originally announced January 2019.

Comments: Submitted to MNRAS, 10 pages, 6 figures, 4 tables. Some typos have been corrected with respect to the previous version

Microlensing flux ratio predictions for Euclid

Authors: Georgios Vernardos

Abstract: Quasar microlensing flux ratios are used to unveil properties of the lenses in large collections of lensed quasars, like the ones expected to be produced by the Euclid survey. This is achieved by using the direct survey products, without any (expensive) follow-up observations or monitoring. First, the theoretical flux ratio distribution of samples of hundreds of mock quasar lenses is calculated fo… ▽ More Quasar microlensing flux ratios are used to unveil properties of the lenses in large collections of lensed quasars, like the ones expected to be produced by the Euclid survey. This is achieved by using the direct survey products, without any (expensive) follow-up observations or monitoring. First, the theoretical flux ratio distribution of samples of hundreds of mock quasar lenses is calculated for different Initial Mass Functions (IMFs) and Sersic radial profiles for the lens compact matter distribution. Then, mock observations are created and compared to the models to recover the underlying one. The most important factor for determining the flux ratio properties of such samples is the value of the smooth matter fraction at the location of the multiple images. Doubly lensed CASTLES-like quasars are the most promising systems to constrain the IMF and the mass components for a sample of lenses. △ Less

Submitted 18 December, 2018; originally announced January 2019.

Comments: 14 pages, 12 figures, 3 tables, accepted for publication in MNRAS

SEAGLE--II: Constraints on feedback models in galaxy formation from massive early type strong lens galaxies

Authors: Sampath Mukherjee, Leon V. E. Koopmans, R. Benton Metcalf, Crescenzo Tortora, Matthieu Schaller, Joop Schaye, Giorgos Vernardos, Fabio Bellagamba

Abstract: We use nine different galaxy formation scenarios in ten cosmological simulation boxes from the EAGLE suite of ΛCDM hydrodynamical simulations to assess the impact of feedback mechanisms in galaxy formation and compare these to observed strong gravitational lenses. To compare observations with simulations, we create strong lenses with $M_\star$ > $10^{11}$ $M_\odot$ with the appropriate resolution… ▽ More We use nine different galaxy formation scenarios in ten cosmological simulation boxes from the EAGLE suite of ΛCDM hydrodynamical simulations to assess the impact of feedback mechanisms in galaxy formation and compare these to observed strong gravitational lenses. To compare observations with simulations, we create strong lenses with $M_\star$ > $10^{11}$ $M_\odot$ with the appropriate resolution and noise level, and model them with an elliptical power-law mass model to constrain their total mass density slope. We also obtain the mass-size relation of the simulated lens-galaxy sample. We find significant variation in the total mass density slope at the Einstein radius and in the projected stellar mass-size relation, mainly due to different implementations of stellar and AGN feedback. We find that for lens selected galaxies, models with either too weak or too strong stellar and/or AGN feedback fail to explain the distribution of observed mass-density slopes, with the counter-intuitive trend that increasing the feedback steepens the mass density slope around the Einstein radius ($\approx$ 3-10 kpc). Models in which stellar feedback becomes inefficient at high gas densities, or weaker AGN feedback with a higher duty cycle, produce strong lenses with total mass density slopes close to isothermal (i.e. -d log(ρ)/d log(r) $\approx$ 2.0) and slope distributions statistically agreeing with observed strong lens galaxies in SLACS and BELLS. Agreement is only slightly worse with the more heterogeneous SL2S lens galaxy sample. Observations of strong-lens selected galaxies thus appear to favor models with relatively weak feedback in massive galaxies. △ Less

Submitted 4 March, 2021; v1 submitted 4 January, 2019; originally announced January 2019.

Comments: Accepted in MNRAS

LinKS: Discovering galaxy-scale strong lenses in the Kilo-Degree Survey using Convolutional Neural Networks

Authors: C. E. Petrillo, C. Tortora, G. Vernardos, L. V. E. Koopmans, G. Verdoes Kleijn, M. Bilicki, N. R. Napolitano, S. Chatterjee, G. Covone, A. Dvornik, T. Erben, F. Getman, B. Giblin, C. Heymans, J. T. A. de Jong, K. Kuijken, P. Schneider, H. Shan, C. Spiniello, A. H. Wright

Abstract: We present a new sample of galaxy-scale strong gravitational-lens candidates, selected from 904 square degrees of Data Release 4 of the Kilo-Degree Survey (KiDS), i.e., the "Lenses in the Kilo-Degree Survey" (LinKS) sample. We apply two Convolutional Neural Networks (ConvNets) to $\sim88\,000$ colour-magnitude selected luminous red galaxies yielding a list of 3500 strong-lens candidates. This list… ▽ More We present a new sample of galaxy-scale strong gravitational-lens candidates, selected from 904 square degrees of Data Release 4 of the Kilo-Degree Survey (KiDS), i.e., the "Lenses in the Kilo-Degree Survey" (LinKS) sample. We apply two Convolutional Neural Networks (ConvNets) to $\sim88\,000$ colour-magnitude selected luminous red galaxies yielding a list of 3500 strong-lens candidates. This list is further down-selected via human inspection. The resulting LinKS sample is composed of 1983 rank-ordered targets classified as "potential lens candidates" by at least one inspector. Of these, a high-grade subsample of 89 targets is identified with potential strong lenses by all inspectors. Additionally, we present a collection of another 200 strong lens candidates discovered serendipitously from various previous ConvNet runs. A straightforward application of our procedure to future Euclid or LSST data can select a sample of $\sim3000$ lens candidates with less than 10 per cent expected false positives and requiring minimal human intervention. △ Less

Submitted 16 January, 2019; v1 submitted 7 December, 2018; originally announced December 2018.

Comments: 19 pages, 11 figures, accepted for publication in MNRAS

A joint microlensing analysis of lensing mass and accretion disc models

Authors: Georgios Vernardos

Abstract: Microlensing of multiply imaged quasars is a unique probe of quasar structure, down to the size of the accretion disc and the central black hole. Flux ratios between close pairs of images of lensed quasars can be used to constrain the accretion disc size and temperature profile. The starting point of any microlensing model is the macromodel of the lens, which provides the convergence and shear val… ▽ More Microlensing of multiply imaged quasars is a unique probe of quasar structure, down to the size of the accretion disc and the central black hole. Flux ratios between close pairs of images of lensed quasars can be used to constrain the accretion disc size and temperature profile. The starting point of any microlensing model is the macromodel of the lens, which provides the convergence and shear values at the location of the multiple images. Here I present a new approach of microlensing modelling independently of the macromodel of the lens. The technique is applied to the close pair of images A1 and A2 of MG 0414+0534, for a set of flux ratios with large variation with respect to wavelength. The inferred accretion disc size and temperature profile measurements, as well as the smooth matter fraction at the location of the images, are quite robust under a wide range of macromodel variations. A case of using purely microlensing data (flux ratios) to constrain the macromodel is also presented. This is a first application of the technique on a fiducial system and set of flux ratios; the method is readily applicable to collections of such objects and can be extended to light curve and/or imaging data. △ Less

Submitted 28 August, 2018; originally announced August 2018.

Comments: 10 pages, 4 figures, accepted for publication in MNRAS

Testing Convolutional Neural Networks for finding strong gravitational lenses in KiDS

Authors: C. E. Petrillo, C. Tortora, S. Chatterjee, G. Vernardos, L. V. E. Koopmans, G. Verdoes Kleijn, N. R. Napolitano, G. Covone, L. S. Kelvin, A. M. Hopkins

Abstract: Convolutional Neural Networks (ConvNets) are one of the most promising methods for identifying strong gravitational lens candidates in survey data. We present two ConvNet lens-finders which we have trained with a dataset composed of real galaxies from the Kilo Degree Survey (KiDS) and simulated lensed sources. One ConvNet is trained with single \textit{r}-band galaxy images, hence basing the class… ▽ More Convolutional Neural Networks (ConvNets) are one of the most promising methods for identifying strong gravitational lens candidates in survey data. We present two ConvNet lens-finders which we have trained with a dataset composed of real galaxies from the Kilo Degree Survey (KiDS) and simulated lensed sources. One ConvNet is trained with single \textit{r}-band galaxy images, hence basing the classification mostly on the morphology. While the other ConvNet is trained on \textit{g-r-i} composite images, relying mostly on colours and morphology. We have tested the ConvNet lens-finders on a sample of 21789 Luminous Red Galaxies (LRGs) selected from KiDS and we have analyzed and compared the results with our previous ConvNet lens-finder on the same sample. The new lens-finders achieve a higher accuracy and completeness in identifying gravitational lens candidates, especially the single-band ConvNet. Our analysis indicates that this is mainly due to improved simulations of the lensed sources. In particular, the single-band ConvNet can select a sample of lens candidates with $\sim40\%$ purity, retrieving 3 out of 4 of the confirmed gravitational lenses in the LRG sample. With this particular setup and limited human intervention, it will be possible to retrieve, in future surveys such as Euclid, a sample of lenses exceeding in size the total number of currently known gravitational lenses. △ Less

Submitted 1 October, 2018; v1 submitted 12 July, 2018; originally announced July 2018.

Comments: 16 pages, 10 figures. Accepted for publication in MNRAS

HST imaging of four gravitationally lensed quasars

Authors: Nicholas F. Bate, Georgios Vernardos, Matthew J. O'Dowd, Daniel M. Neri-Larios, Rachel L. Webster, David J. E. Floyd, Robert L. Barone-Nugent, Kathleen Labrie, Anthea L. King, Suk Yee Yong

Abstract: We present new HST WFPC3 imaging of four gravitationally lensed quasars: MG 0414+0534; RXJ 0911+0551; B 1422+231; WFI J2026-4536. In three of these systems we detect wavelength-dependent microlensing, which we use to place constraints on the sizes and temperature profiles of the accretion discs in each quasar. Accretion disc radius is assumed to vary with wavelength according to the power-law rela… ▽ More We present new HST WFPC3 imaging of four gravitationally lensed quasars: MG 0414+0534; RXJ 0911+0551; B 1422+231; WFI J2026-4536. In three of these systems we detect wavelength-dependent microlensing, which we use to place constraints on the sizes and temperature profiles of the accretion discs in each quasar. Accretion disc radius is assumed to vary with wavelength according to the power-law relationship $r\propto λ^p$, equivalent to a radial temperature profile of $T\propto r^{-1/p}$. The goal of this work is to search for deviations from standard thin disc theory, which predicts that radius goes as wavelength to the power $p=4/3$. We find a wide range of power-law indices, from $p=1.4^{+0.5}_{-0.4}$ in B 1422+231 to $p=2.3^{+0.5}_{-0.4}$ in WFI J2026-4536. The measured value of $p$ appears to correlate with the strength of the wavelength-dependent microlensing. We explore this issue with mock simulations using a fixed accretion disc with $p=1.5$, and find that cases where wavelength-dependent microlensing is small tend to under-estimate the value of $p$. This casts doubt on previous ensemble single-epoch measurements which have favoured low values using samples of lensed quasars that display only moderate chromatic effects. Using only our systems with strong chromatic microlensing we prefer $p>4/3$, corresponding to shallower temperature profiles than expected from standard thin disc theory. △ Less

Submitted 10 July, 2018; originally announced July 2018.

Comments: 20 pages, 11 figures, accepted for publication in MNRAS

KiDS-SQuaD: The KiDS Strongly lensed Quasar Detection project

Authors: C. Spiniello, A. Agnello, N. R. Napolitano, A. V. Sergeyev, F. I. Getman, C. Tortora, M. Spavone, M. Bilicki, H. Buddelmeijer, L. V. E. Koopmans, K. Kuijken, G. Vernardos, E. Bannikova, M. Capaccioli

Abstract: New methods have been recently developed to search for strong gravitational lenses, in particular lensed quasars, in wide-field imaging surveys. Here, we compare the performance of three different, morphology- and photometry- based methods to find lens candidates over the Kilo-Degree Survey (KiDS) DR3 footprint (440 deg$^2$). The three methods are: i) a multiplet detection in KiDS-DR3 and/or Gaia-… ▽ More New methods have been recently developed to search for strong gravitational lenses, in particular lensed quasars, in wide-field imaging surveys. Here, we compare the performance of three different, morphology- and photometry- based methods to find lens candidates over the Kilo-Degree Survey (KiDS) DR3 footprint (440 deg$^2$). The three methods are: i) a multiplet detection in KiDS-DR3 and/or Gaia-DR1, ii) direct modeling of KiDS cutouts and iii) positional offsets between different surveys (KiDS-vs-Gaia, Gaia-vs-2MASS), with purpose-built astrometric recalibrations. The first benchmark for the methods has been set by the recovery of known lenses. We are able to recover seven out of ten known lenses and pairs of quasars observed in the KiDS DR3 footprint, or eight out of ten with improved selection criteria and looser colour pre-selection. This success rate reflects the combination of all methods together, which, taken individually, performed significantly worse (four lenses each). One movelty of our analysis is that the comparison of the performances of the different methods has revealed the pros and cons of the approaches and, most of all, the complementarities. We finally provide a list of high-grade candidates found by one or more methods, awaiting spectroscopic follow-up for confirmation. Of these, KiDS 1042+0023 is to our knowledge the first confirmed lensed quasar from KiDS, exhibiting two quasar spectra at the same source redshift at either sides of a red galaxy, with uniform flux-ratio $f\approx1.25$ over the wavelength range $0.45μ\mathrm{m}<λ<0.75μ\mathrm{m}.$ △ Less

Submitted 11 July, 2018; v1 submitted 31 May, 2018; originally announced May 2018.

Comments: 12 pages, 4 figures, 4 tables, accepted for publication in MNRAS

The Impact of Microlensing on the Standardisation of Strongly Lensed Type Ia Supernovae

Authors: Max Foxley-Marrable, Thomas E. Collett, Georgios Vernardos, Daniel A. Goldstein, David Bacon

Abstract: We investigate the effect of microlensing on the standardisation of strongly lensed Type Ia supernovae (GLSNe Ia). We present predictions for the amount of scatter induced by microlensing across a range of plausible strong lens macromodels. We find that lensed images in regions of low convergence, shear and stellar density are standardisable, where the microlensing scatter is < 0.15 magnitudes, co… ▽ More We investigate the effect of microlensing on the standardisation of strongly lensed Type Ia supernovae (GLSNe Ia). We present predictions for the amount of scatter induced by microlensing across a range of plausible strong lens macromodels. We find that lensed images in regions of low convergence, shear and stellar density are standardisable, where the microlensing scatter is < 0.15 magnitudes, comparable to the intrinsic dispersion of for a typical SN Ia. These standardisable configurations correspond to asymmetric lenses with an image located far outside the Einstein radius of the lens. Symmetric and small Einstein radius lenses (< 0.5 arcsec) are not standardisable. We apply our model to the recently discovered GLSN Ia iPTF16geu and find that the large discrepancy between the observed flux and the macromodel predictions from More et al. (2017) cannot be explained by microlensing alone. Using the mock GLSNe Ia catalogue of Goldstein et al. (2017), we predict that ~ 22% of GLSNe Ia discovered by LSST will be standardisable, with a median Einstein radius of 0.9 arcseconds and a median time-delay of 41 days. By breaking the mass-sheet degeneracy the full LSST GLSNe Ia sample will be able to detect systematics in H0 at the 0.5% level. △ Less

Submitted 25 May, 2018; v1 submitted 21 February, 2018; originally announced February 2018.

Comments: 11 pages, 8 Figures. Accepted by MNRAS May 17 2018

SEAGLE - I: A pipeline for simulating and modeling strong lenses from cosmological hydrodynamic simulations

Authors: Sampath Mukherjee, Léon. V. E. Koopmans, Robert Benton Metcalf, Nicolas Tessore, Crescenzo Tortora, Matthieu Schaller, Joop Schaye, Robert A. Crain, Giorgos Vernardos, Fabio Bellagamba, Tom Theuns

Abstract: In this paper we introduce the SEAGLE (i.e. Simulating EAGLE LEnses) program, that approaches the study of galaxy formation through strong gravitational lensing, using a suite of high-resolution hydrodynamic simulations, Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. We introduce the simulation and analysis pipeline and present the first set of results from our analysis… ▽ More In this paper we introduce the SEAGLE (i.e. Simulating EAGLE LEnses) program, that approaches the study of galaxy formation through strong gravitational lensing, using a suite of high-resolution hydrodynamic simulations, Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. We introduce the simulation and analysis pipeline and present the first set of results from our analysis of early-type galaxies. We identify and extract an ensemble of simulated lens galaxies and use the GLAMER ray-tracing lensing code to create mock lenses similar to those observed in the SLACS and SL2S surveys, using a range of source parameters and galaxy orientations, including observational effects such as the Point-Spread-Function (PSF), pixelization and noise levels, representative of single-orbit observations with the Hubble Space Telescope (HST) using the ACS-F814W filter. We subsequently model these mock lenses using the code LENSED, treating them in the same way as observed lenses. We also estimate the mass model parameters directly from the projected surface mass density of the simulated galaxy, using an identical mass model family. We perform a three-way comparison of all the measured quantities with real lenses. We find the average total density slope of EAGLE lenses, $t=2.26\; (0.25\; \rm{rms})$ to be higher than SL2S, $t=2.16$ or SLACS, $t=2.08$. We find a very strong correlation between the external shear ($γ$) and the complex ellipticity ($ε$), with $γ\sim ε/4$. This correlation indicates a degeneracy in the lens mass modeling. We also see a dispersion between lens modeling and direct fitting results, indicating systematical biases. △ Less

Submitted 29 June, 2018; v1 submitted 19 February, 2018; originally announced February 2018.

Comments: Accepted in MNRAS, Highlights: Figure- 6, 12, 13

The Strong Gravitational Lens Finding Challenge

Authors: R. Benton Metcalf, M. Meneghetti, Camille Avestruz, Fabio Bellagamba, Clécio R. Bom, Emmanuel Bertin, Rémi Cabanac, F. Courbin, Andrew Davies, Etienne Decencière, Rémi Flamary, Raphael Gavazzi, Mario Geiger, Philippa Hartley, Marc Huertas-Company, Neal Jackson, Eric Jullo, Jean-Paul Kneib, Léon V. E. Koopmans, François Lanusse, Chun-Liang Li, Quanbin Ma, Martin Makler, Nan Li, Matthew Lightman , et al. (11 additional authors not shown)

Abstract: Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives a… ▽ More Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged) △ Less

Submitted 20 March, 2019; v1 submitted 10 February, 2018; originally announced February 2018.

Comments: 22 pages, 16 figures, 4 tables, accepted version for A&A

Journal ref: A&A 625, A119 (2019)

Finding Strong Gravitational Lenses in the Kilo Degree Survey with Convolutional Neural Networks

Authors: C. E. Petrillo, C. Tortora, S. Chatterjee, G. Vernardos, L. V. E. Koopmans, G. Verdoes Kleijn, N. R. Napolitano, G. Covone, P. Schneider, A. Grado, J. McFarland

Abstract: The volume of data that will be produced by new-generation surveys requires automatic classification methods to select and analyze sources. Indeed, this is the case for the search for strong gravitational lenses, where the population of the detectable lensed sources is only a very small fraction of the full source population. We apply for the first time a morphological classification method based… ▽ More The volume of data that will be produced by new-generation surveys requires automatic classification methods to select and analyze sources. Indeed, this is the case for the search for strong gravitational lenses, where the population of the detectable lensed sources is only a very small fraction of the full source population. We apply for the first time a morphological classification method based on a Convolutional Neural Network (CNN) for recognizing strong gravitational lenses in $255$ square degrees of the Kilo Degree Survey (KiDS), one of the current-generation optical wide surveys. The CNN is currently optimized to recognize lenses with Einstein radii $\gtrsim 1.4$ arcsec, about twice the $r$-band seeing in KiDS. In a sample of $21789$ colour-magnitude selected Luminous Red Galaxies (LRG), of which three are known lenses, the CNN retrieves 761 strong-lens candidates and correctly classifies two out of three of the known lenses. The misclassified lens has an Einstein radius below the range on which the algorithm is trained. We down-select the most reliable 56 candidates by a joint visual inspection. This final sample is presented and discussed. A conservative estimate based on our results shows that with our proposed method it should be possible to find $\sim100$ massive LRG-galaxy lenses at $z\lsim 0.4$ in KiDS when completed. In the most optimistic scenario this number can grow considerably (to maximally $\sim$2400 lenses), when widening the colour-magnitude selection and training the CNN to recognize smaller image-separation lens systems. △ Less

Submitted 23 August, 2017; v1 submitted 24 February, 2017; originally announced February 2017.

Comments: 24 pages, 17 figures. Published in MNRAS

The Herschel-ATLAS: a sample of 500μm-selected lensed galaxies over 600 square degrees

Authors: M. Negrello, S. Amber, A. Amvrosiadis, Z. -Y. Cai, A. Lapi, J. Gonzalez-Nuevo, G. De Zotti, C. Furlanetto, S. Maddox, M. Allen, T. Bakx, R. S. Bussmann, A. Cooray, G. Covone, L. Danese, H. Dannerbauer, H. Fu, J. Greenslade, M. Gurwell, R. Hopwood, L. V. E. Koopmans, N. Napolitano, H. Nayyeri, A. Omont, C. E. Petrillo , et al. (21 additional authors not shown)

Abstract: We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500μm extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly support… ▽ More We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500μm extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500μm and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: 25 pages, 17 figures, accepted for publication in MNRAS

Data Compression in the Petascale Astronomy Era: a GERLUMPH case study

Authors: Dany Vohl, Christopher J. Fluke, Georgios Vernardos

Abstract: As the volume of data grows, astronomers are increasingly faced with choices on what data to keep -- and what to throw away. Recent work evaluating the JPEG2000 (ISO/IEC 15444) standards as a future data format standard in astronomy has shown promising results on observational data. However, there is still a need to evaluate its potential on other type of astronomical data, such as from numerical… ▽ More As the volume of data grows, astronomers are increasingly faced with choices on what data to keep -- and what to throw away. Recent work evaluating the JPEG2000 (ISO/IEC 15444) standards as a future data format standard in astronomy has shown promising results on observational data. However, there is still a need to evaluate its potential on other type of astronomical data, such as from numerical simulations. GERLUMPH (the GPU-Enabled High Resolution cosmological MicroLensing parameter survey) represents an example of a data intensive project in theoretical astrophysics. In the next phase of processing, the ~27 terabyte GERLUMPH dataset is set to grow by a factor of 100 -- well beyond the current storage capabilities of the supercomputing facility on which it resides. In order to minimise bandwidth usage, file transfer time, and storage space, this work evaluates several data compression techniques. Specifically, we investigate off-the-shelf and custom lossless compression algorithms as well as the lossy JPEG2000 compression format. Results of lossless compression algorithms on GERLUMPH data products show small compression ratios (1.35:1 to 4.69:1 of input file size) varying with the nature of the input data. Our results suggest that JPEG2000 could be suitable for other numerical datasets stored as gridded data or volumetric data. When approaching lossy data compression, one should keep in mind the intended purposes of the data to be compressed, and evaluate the effect of the loss on future analysis. In our case study, lossy compression and a high compression ratio do not significantly compromise the intended use of the data for constraining quasar source profiles from cosmological microlensing. △ Less

Submitted 4 May, 2016; v1 submitted 21 May, 2015; originally announced May 2015.

Comments: 15 pages, 9 figures, 5 tables. Published in the Special Issue of Astronomy & Computing on The future of astronomical data formats

Journal ref: Astronomy & Computing, 12, 200 (2015)

GERLUMPH Data Release 2: 2.5 billion simulated microlensing light curves

Authors: Georgios Vernardos, Christopher J. Fluke, Nicholas F. Bate, Darren Croton, Dany Vohl

Abstract: In the upcoming synoptic all--sky survey era of astronomy, thousands of new multiply imaged quasars are expected to be discovered and monitored regularly. Light curves from the images of gravitationally lensed quasars are further affected by superimposed variability due to microlensing. In order to disentangle the microlensing from the intrinsic variability of the light curves, the time delays bet… ▽ More In the upcoming synoptic all--sky survey era of astronomy, thousands of new multiply imaged quasars are expected to be discovered and monitored regularly. Light curves from the images of gravitationally lensed quasars are further affected by superimposed variability due to microlensing. In order to disentangle the microlensing from the intrinsic variability of the light curves, the time delays between the multiple images have to be accurately measured. The resulting microlensing light curves can then be analyzed to reveal information about the background source, such as the size of the quasar accretion disc. In this paper we present the most extensive and coherent collection of simulated microlensing light curves; we have generated $>2.5$ billion light curves using the GERLUMPH high resolution microlensing magnification maps. Our simulations can be used to: train algorithms to measure lensed quasar time delays, plan future monitoring campaigns, and study light curve properties throughout parameter space. Our data are openly available to the community and are complemented by online eResearch tools, located at http://gerlumph.swin.edu.au . △ Less

Submitted 2 March, 2015; originally announced March 2015.

Comments: 9 pages, 4 figures, accepted for publication in ApJSS, data available from http://gerlumph.swin.edu.au

The effect of macromodel uncertainties on microlensing modelling of lensed quasars

Authors: Georgios Vernardos, Christopher J. Fluke

Abstract: Cosmological gravitational microlensing has been proven to be a powerful tool to constrain the structure of multiply imaged quasars, especially the accretion disc and central supermassive black-hole system. However, the derived constraints on models may be affected by large systematic errors introduced in the various stages of modelling, namely, the macromodels, the microlensing magnification maps… ▽ More Cosmological gravitational microlensing has been proven to be a powerful tool to constrain the structure of multiply imaged quasars, especially the accretion disc and central supermassive black-hole system. However, the derived constraints on models may be affected by large systematic errors introduced in the various stages of modelling, namely, the macromodels, the microlensing magnification maps, and the convolution with realistic disc profiles. In particular, it has been known that different macromodels of the galaxy lens that fit the observations equally well, can lead to different values of convergence, $κ$, and shear, $γ$, required to generate magnification maps. So far, 25 microlensed quasars have been studied using microlensing techniques, where each system has been modelled and analyzed individually, or in small samples. This is about to change due to the upcoming synoptic all-sky surveys, which are expected to discover thousands of quasars suitable for microlensing studies. In the present study we investigate the connection between macromodels of the galaxy lens and microlensing magnification maps throughout the parameter space in preparation for future studies of large statistical samples of systems displaying microlensing. In particular, we use 55,900 maps produced by the GERLUMPH parameter survey (available online at http://gerlumph.swin.edu.au ) and identify regions of parameter space where macromodel uncertainties ($Δκ$,$Δγ$) lead to statistically different magnification maps. Strategies for mitigating the effect of $Δκ$,$Δγ$ uncertainties are discussed in order to understand and control this potential source of systematic errors in accretion disc constraints derived from microlensing. △ Less

Submitted 4 September, 2014; originally announced September 2014.

Comments: 14 pages, 10 figures, 1 table, accepted for publication in MNRAS

Adventures in the microlensing cloud: large datasets, eResearch tools, and GPUs

Authors: Georgios Vernardos, Christopher J Fluke

Abstract: As astronomy enters the petascale data era, astronomers are faced with new challenges relating to storage, access and management of data. A shift from the traditional approach of combining data and analysis at the desktop to the use of remote services, pushing the computation to the data, is now underway. In the field of cosmological gravitational microlensing, future synoptic all--sky surveys are… ▽ More As astronomy enters the petascale data era, astronomers are faced with new challenges relating to storage, access and management of data. A shift from the traditional approach of combining data and analysis at the desktop to the use of remote services, pushing the computation to the data, is now underway. In the field of cosmological gravitational microlensing, future synoptic all--sky surveys are expected to bring the number of multiply imaged quasars from the few tens that are currently known to a few thousands. This inflow of observational data, together with computationally demanding theoretical modelling via the production of microlensing magnification maps, requires a new approach. We present our technical solutions to supporting the GPU-Enabled, High Resolution cosmological MicroLensing parameter survey (GERLUMPH). This extensive dataset for cosmological microlensing modelling comprises over 70,000 individual magnification maps and ${\sim}10^6$ related results. We describe our approaches to hosting, organizing, and serving ${\sim}$30 Terabytes of data and metadata products. We present a set of online analysis tools developed with PHP, JavaScript and WebGL to support access and analysis of GELRUMPH data in a Web browser. We discuss our use of graphics processing units (GPUs) to accelerate data production, and we release the core of the GPU-D direct inverse ray--shooting code (Thompson et al., 2010; Astrophysics Source Code Library, record ascl:1403.001) used to generate the magnification maps. All of the GERLUMPH data and tools are available online from http://gerlumph.swin.edu.au . This project made use of gSTAR, the GPU Supercomputer for Theoretical Astrophysical Research. △ Less

Submitted 2 June, 2014; originally announced June 2014.

Comments: 24 pages, 14 figures, 6 tables, Astronomy & Computing in press

GERLUMPH Data Release 1: High-resolution cosmological microlensing magnification maps and eResearch tools

Authors: Georgios Vernardos, Christopher J Fluke, Nicholas F Bate, Darren Croton

Abstract: As synoptic all-sky surveys begin to discover new multiply lensed quasars, the flow of data will enable statistical cosmological microlensing studies of sufficient size to constrain quasar accretion disc and supermassive black hole properties. In preparation for this new era, we are undertaking the GPU-Enabled, High Resolution cosmological MicroLensing parameter survey (GERLUMPH). We present here… ▽ More As synoptic all-sky surveys begin to discover new multiply lensed quasars, the flow of data will enable statistical cosmological microlensing studies of sufficient size to constrain quasar accretion disc and supermassive black hole properties. In preparation for this new era, we are undertaking the GPU-Enabled, High Resolution cosmological MicroLensing parameter survey (GERLUMPH). We present here the GERLUMPH Data Release 1, which consists of 12342 high resolution cosmological microlensing magnification maps and provides the first uniform coverage of the convergence, shear and smooth matter fraction parameter space. We use these maps to perform a comprehensive numerical investigation of the mass-sheet degeneracy, finding excellent agreement with its predictions. We study the effect of smooth matter on microlensing induced magnification fluctuations. In particular, in the minima and saddle-point regions, fluctuations are enhanced only along the critical line, while in the maxima region they are always enhanced for high smooth matter fractions (~0.9). We describe our approach to data management, including the use of an SQL database with a Web interface for data access and online analysis, obviating the need for individuals to download large volumes of data. In combination with existing observational databases and online applications, the GERLUMPH archive represents a fundamental component of a new microlensing eResearch cloud. Our maps and tools are publicly available at http://gerlumph.swin.edu.au/. △ Less

Submitted 29 January, 2014; originally announced January 2014.

Comments: 9 pages, 9 figures, 2 tables, accepted for publication in ApJ

Journal ref: Vernardos et al., 2014, ApJS, 211, 16

A new parameter space study of cosmological microlensing

Authors: Georgios Vernardos, Christopher J. Fluke

Abstract: Cosmological gravitational microlensing is a useful technique for understanding the structure of the inner parts of a quasar, especially the accretion disk and the central supermassive black hole. So far, most of the cosmological microlensing studies have focused on single objects from ~90 currently known lensed quasars. However, present and planned all-sky surveys are expected to discover thousan… ▽ More Cosmological gravitational microlensing is a useful technique for understanding the structure of the inner parts of a quasar, especially the accretion disk and the central supermassive black hole. So far, most of the cosmological microlensing studies have focused on single objects from ~90 currently known lensed quasars. However, present and planned all-sky surveys are expected to discover thousands of new lensed systems. Using a graphics processing unit (GPU) accelerated ray-shooting code, we have generated 2550 magnification maps uniformly across the convergence (κ) and shear (γ) parameter space of interest to microlensing. We examine the effect of random realizations of the microlens positions on map properties such as the magnification probability distribution (MPD). It is shown that for most of the parameter space a single map is representative of an average behaviour. All of the simulations have been carried out on the GPU-Supercomputer for Theoretical Astrophysics Research (gSTAR). △ Less

Submitted 16 June, 2013; originally announced June 2013.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Journal ref: Vernardos & Fluke, 2013, MNRAS, 434, 832