Abstract

Modeling the mass distributions of strong gravitational lenses is often necessary in order to use them as astrophysical and cosmological probes. With the large number of lens systems (≳105) expected from upcoming surveys, it is timely to explore efficient modeling approaches beyond traditional Markov chain Monte Carlo techniques that are time consuming. We train a convolutional neural network (CNN) on images of galaxy-scale lens systems to predict the five parameters of the singular isothermal ellipsoid (SIE) mass model (lens center x and y, complex ellipticity ex and ey, and Einstein radius θE). To train the network we simulate images based on real observations from the Hyper Suprime-Cam Survey for the lens galaxies and from the Hubble Ultra Deep Field as lensed galaxies. We tested different network architectures and the effect of different data sets, such as using only double or quad systems defined based on the source center and using different input distributions of θE. We find that the CNN performs well, and with the network trained on both doubles and quads with a uniform distribution of θE > 0.5″ we obtain the following median values with 1σ scatter: Δx = (0.00−0.30+0.30)″, Δy = (0.00−0.29+0.30)″, ΔθE = (0.07−0.12+0.29)″, Δex = −0.01−0.09+0.08, and Δey = 0.00−0.09+0.08. The bias in θE is driven by systems with small θE. Therefore, when we further predict the multiple lensed image positions and time-delays based on the network output, we apply the network to the sample limited to θE > 0.8″. In this case the offset between the predicted and input lensed image positions is (0.00−0.29+0.29)″ and (0.00−0.31+0.32)″ for the x and y coordinates, respectively. For the fractional difference between the predicted and true time-delay, we obtain 0.04−0.05+0.27. Our CNN model is able to predict the SIE parameter values in fractions of a second on a single CPU, and with the output we can predict the image positions and time-delays in an automated way, such that we are able to process efficiently the huge amount of expected galaxy-scale lens detections in the near future.

Highlights

  • Strong gravitational lensing has become a very powerful tool for probing various properties of the Universe

  • We give an overview of the different data set assumptions in Table 1, as well as the best hyperparameter values that depend on the assumed data set

  • We presented a convolutional neural network to model in a fully automated way and very quickly the mass distribution of galaxy-scale strong lens systems by assuming a singular isothermal ellipsoid (SIE) profile

Read more

Summary

Introduction

Strong gravitational lensing has become a very powerful tool for probing various properties of the Universe. To shed light on this debate, one possibility is to observe the SN Ia spectroscopically at very early stages, which is normally difficult because SN detections are often close to peak luminosity, past the early phase If this SN is lensed, we can use the position of the first appearing image, together with a mass model of the underlying lens galaxy, to predict the position and time when the images will appear. Based on newer upcoming surveys like the LSST, which will target around 20 000 deg of the southern hemisphere in six different filters (u, g, r, i, z, y), together with the Euclid imaging survey from space operated by the European Space Agency (ESA; Laureijs et al 2011), we expect billions of galaxy images containing on the order of one hundred thousand lenses (Collett 2015) To deal with this huge amount of images there are ongoing efforts to develop fast and automated algorithms to find lenses in the first place. Each quoted parameter estimate is the median of its 1D marginalized posterior probability density function, and the quoted uncertainties show the 16th and 84th percentiles (i.e., the bounds of a 68% credible interval)

Simulation of strongly lensed images
Lens galaxies from HSC
Sources from HUDF
Mock lens systems
Neural networks and their architecture
Results
Data set containing double or quads only
Comparison to other modeling codes
Summary and conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call