-
Predicting High Magnification Events in Microlensed Quasars in the Era of LSST using Recurrent Neural Networks
Authors:
Joshua Fagin,
Eric Paic,
Favio Neira,
Henry Best,
Timo Anguita,
Martin Millon,
Matthew O'Dowd,
Dominique Sluse,
Georgios Vernardos
Abstract:
Upcoming wide field surveys such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) will monitor thousands of strongly lensed quasars over a 10-year period. Many of these monitored quasars will undergo high magnification events (HMEs) through microlensing as the accretion disk crosses a caustic, places of infinite magnification. Microlensing allows us to map the inner regions of the…
▽ More
Upcoming wide field surveys such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) will monitor thousands of strongly lensed quasars over a 10-year period. Many of these monitored quasars will undergo high magnification events (HMEs) through microlensing as the accretion disk crosses a caustic, places of infinite magnification. Microlensing allows us to map the inner regions of the accretion disk as it crosses a caustic, even at large cosmological distances. The observational cadences of LSST are not ideal for probing the inner regions of the accretion disk, so there is a need to predict HMEs as early as possible to trigger high-cadence multi-band or spectroscopic follow-up observations. Here we simulate a diverse and realistic sample of 10-year quasar microlensing light curves to train a recurrent neural network (RNN) to predict HMEs before they occur by classifying the location of the peaks at each time step. This is the first deep learning approach to predict HMEs. We give estimates at how well we expect to predict HME peaks during LSST and benchmark how our metrics change with different cadence strategies. With LSST-like observations, we can predict approximately 55% of HME peaks corresponding to tens to hundreds per year and a false positive rate of around 20% compared to the number of HMEs. Our network can be continuously applied throughout the LSST survey, providing crucial alerts to optimize follow-up resources.
△ Less
Submitted 13 September, 2024;
originally announced September 2024.
-
Probing compact dark matter objects with microlensing in gravitationally lensed quasars
Authors:
Petra Awad,
James H. H. Chan,
Martin Millon,
Frederic Courbin,
Eric Paic
Abstract:
The microlensing signal in the light curves of gravitationally lensed quasars can shed light on the dark matter (DM) composition in their lensing galaxies. Here, we investigate a sample of six lensed quasars from the most recent and best COSMOGRAIL observations: HE~1104$-$1805, HE~0435$-$1223, RX~J1131$-$1231, WFI~2033$-$4723, PG~1115$+$080, and J1206$+$4332, yielding a total of eight microlensing…
▽ More
The microlensing signal in the light curves of gravitationally lensed quasars can shed light on the dark matter (DM) composition in their lensing galaxies. Here, we investigate a sample of six lensed quasars from the most recent and best COSMOGRAIL observations: HE~1104$-$1805, HE~0435$-$1223, RX~J1131$-$1231, WFI~2033$-$4723, PG~1115$+$080, and J1206$+$4332, yielding a total of eight microlensing light curves, when combining independent image pairs and typically spanning ten years. We explore the microlensing signals to determine whether the standard assumptions on the stellar populations are sufficient to account for the amplitudes of the measured signals. We use the most detailed lens models to date from the H0LiCOW/TDCOSMO collaboration to generate simulated microlensing light curves. Finally, we propose a methodology based on the Kolmogorov-Smirnov test to verify whether the observed microlensing amplitudes in our data are compatible with the most standard scenario, whereby galaxies are composed of stars as compact bodies and smoothly distributed DM. Given our current sample, we show that the standard scenario cannot be rejected, in contrast with previous results by Hawkins (2020a), claiming that a population of stellar mass primordial black holes (PBHs) is necessary to explain the amplitude of the microlensing signals in lensed quasar light curves. We further estimate the number of microlensing light curves needed to distinguish between the standard scenario with stellar microlensing and a scenario that describes all the DM contained in galaxies in the form of compact objects such as PBHs, with a mean mass of $0.2M_{\odot}$. We find that about 900 microlensing curves from the Rubin Observatory will be sufficient to discriminate between the two extreme scenarios at a 95\% confidence level.
△ Less
Submitted 3 April, 2023;
originally announced April 2023.
-
The impact of human expert visual inspection on the discovery of strong gravitational lenses
Authors:
Karina Rojas,
Thomas E. Collett,
Daniel Ballard,
Mark R. Magee,
Simon Birrer,
Elizabeth Buckley-Geer.,
James H. H. Chan,
Benjamin Clément,
José M. Diego,
Fabrizio Gentile,
Jimena González,
Rémy Joseph,
Jorge Mastache,
Stefan Schuldt,
Crescenzo Tortora,
Tomás Verdugo,
Aprajita Verma,
Tansu Daylan,
Martin Millon,
Neal Jackson,
Simon Dye,
Alejandra Melo,
Guillaume Mahler,
Ricardo L. C. Ogando,
Frédéric Courbin
, et al. (31 additional authors not shown)
Abstract:
We investigate the ability of human 'expert' classifiers to identify strong gravitational lens candidates in Dark Energy Survey like imaging. We recruited a total of 55 people that completed more than 25$\%$ of the project. During the classification task, we present to the participants 1489 images. The sample contains a variety of data including lens simulations, real lenses, non-lens examples, an…
▽ More
We investigate the ability of human 'expert' classifiers to identify strong gravitational lens candidates in Dark Energy Survey like imaging. We recruited a total of 55 people that completed more than 25$\%$ of the project. During the classification task, we present to the participants 1489 images. The sample contains a variety of data including lens simulations, real lenses, non-lens examples, and unlabeled data. We find that experts are extremely good at finding bright, well-resolved Einstein rings, whilst arcs with $g$-band signal-to-noise less than $\sim$25 or Einstein radii less than $\sim$1.2 times the seeing are rarely recovered. Very few non-lenses are scored highly. There is substantial variation in the performance of individual classifiers, but they do not appear to depend on the classifier's experience, confidence or academic position. These variations can be mitigated with a team of 6 or more independent classifiers. Our results give confidence that humans are a reliable pruning step for lens candidates, providing pure and quantifiably complete samples for follow-up studies.
△ Less
Submitted 25 April, 2023; v1 submitted 9 January, 2023;
originally announced January 2023.
-
Evidence for a milliparsec-separation Supermassive Binary Black Hole with quasar microlensing
Authors:
M. Millon,
C. Dalang,
C. Lemon,
D. Sluse,
E. Paic,
J. H. H. Chan,
F. Courbin
Abstract:
We report periodic oscillations in the 15-year long optical light curve of the gravitationally lensed quasar QJ0158-4325. The signal is enhanced during a high magnification microlensing event undergone by the fainter lensed image of the quasar, between 2003 and 2010. We measure a period of $P_{o}=172.6\pm0.9$ days. We explore four scenarios to explain the origin of the periodicity: 1- the high mag…
▽ More
We report periodic oscillations in the 15-year long optical light curve of the gravitationally lensed quasar QJ0158-4325. The signal is enhanced during a high magnification microlensing event undergone by the fainter lensed image of the quasar, between 2003 and 2010. We measure a period of $P_{o}=172.6\pm0.9$ days. We explore four scenarios to explain the origin of the periodicity: 1- the high magnification microlensing event is due to a binary star in the lensing galaxy, 2- QJ0158-4325 contains a massive binary black hole system in its final dynamical stage before merging, 3- the quasar accretion disk contains a bright inhomogeneity in Keplerian motion around the black hole, and 4- the accretion disk is in precession. Among these four scenarios, only a supermassive binary black hole can account for both the short observed period and the amplitude of the signal, through the oscillation of the accretion disk towards and away from high-magnification regions of a microlensing caustic. The short measured period implies that the semi-long axis of the orbit is $\sim10^{-3}$pc, and the coalescence timescale is $t_{coal}\sim1000$ years, assuming that the decay of the orbit is solely powered by the emission of gravitational waves. The probability of observing a system so close to coalescence suggests either a much larger population of supermassive black hole binaries than predicted, or, more likely, that some other mechanism significantly increases the coalescence timescale. Three tests of the binary black hole hypothesis include: i) the recurrence of oscillations in photometric monitoring during any future microlensing events in either image, ii) spectroscopic detection of Doppler shifts (up to 0.01$c$), and iii) the detection of gravitational waves through Pulsar Timing Array experiments, such as the SKA, which will have the sensitivity to detect the $\sim$100 nano-hertz emission.
△ Less
Submitted 7 December, 2022; v1 submitted 1 July, 2022;
originally announced July 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.
-
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.
-
J1721$+$8842: a gravitationally lensed binary quasar with a proximate damped Lyman-$α$ absorber
Authors:
C. Lemon,
M. Millon,
D. Sluse,
F. Courbin,
M. Auger,
J. Chan,
E. Paic,
A. Agnello
Abstract:
High-redshift binary quasars provide key insights into mergers and quasar activity, and are useful tools for probing the spatial kinematics and chemistry of galaxies along the line-of-sight. However, only three sub-10-kpc binaries have been confirmed above $z=1$. Gravitational lensing would provide a way to easily resolve such binaries, study them in higher resolution, and provide more sightlines,…
▽ More
High-redshift binary quasars provide key insights into mergers and quasar activity, and are useful tools for probing the spatial kinematics and chemistry of galaxies along the line-of-sight. However, only three sub-10-kpc binaries have been confirmed above $z=1$. Gravitational lensing would provide a way to easily resolve such binaries, study them in higher resolution, and provide more sightlines, though the required alignment with a massive foreground galaxy is rare. Through image deconvolution of StanCam Nordic Optical Telescope (NOT) monitoring data, we reveal two further point sources in the known, $z \approx 2.38$, quadruply lensed quasar (quad), J1721+8842. An ALFOSC/NOT long-slit spectrum shows that the brighter of these two sources is a quasar with $z = 2.369 \pm 0.007$ based on the C III] line, while the C III] redshift of the quad is $z = 2.364 \pm 0.003$. Lens modelling using point source positions rules out a single source model, favouring an isothermal lens mass profile with two quasar sources separated by $\sim6.0$ kpc (0.73$^{\prime \prime}$) in projection. Given the resolving ability from lensing and current lensed quasar statistics, this discovery suggests a large population of undiscovered, unlensed sub-10-kpc binaries. We also analyse spectra of two images of the quad, showing narrow Ly$α$ emission within the trough of a proximate damped Ly$α$ absorber (PDLA). An apparent mismatch between the continuum and narrow line flux ratios provides a new potential tool for simultaneously studying microlensing and the quasar host galaxy. Signs of the PDLA are also seen in the second source, however a deeper spectrum is still required to confirm this. Thanks to the multiple lines-of-sight from lensing and two quasar sources, this system offers simultaneous sub-parsec and kpc-scale probes of a PDLA.
△ Less
Submitted 2 September, 2021;
originally announced September 2021.
-
TDCOSMO II: 6 new time delays in lensed quasars from high-cadence monitoring at the MPIA 2.2m telescope
Authors:
M. Millon,
F. Courbin,
V. Bonvin,
E. Buckley-Geer,
C. D. Fassnacht,
J. Frieman,
P. J. Marshall,
S. H. Suyu,
T. Treu,
T. Anguita,
V. Motta,
A. Agnello,
J. H. H. Chan,
D. C. -Y Chao,
M. Chijani,
D. Gilman,
K. Gilmore,
C. Lemon,
J. R. Lucey,
A. Melo,
E. Paic,
K. Rojas,
D. Sluse,
P. R. Williams,
A. Hempel
, et al. (3 additional authors not shown)
Abstract:
We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the Max Planck Institute for Astrophysics (MPIA) 2.2 m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain…
▽ More
We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the Max Planck Institute for Astrophysics (MPIA) 2.2 m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measured time delays between all pairs of multiple images with only one or two seasons of monitoring with the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate was obtained for the delay between image A and image B of DES 0407-5006, where $τ_{AB} = -128.4^{+3.5}_{-3.8}$ d (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combined our high-cadence data with measurements from decade-long light curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 d on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily.
△ Less
Submitted 22 October, 2020; v1 submitted 17 June, 2020;
originally announced June 2020.
-
COSMOGRAIL XIX: Time delays in 18 strongly lensed quasars from 15 years of optical monitoring
Authors:
M. Millon,
F. Courbin,
V. Bonvin,
E. Paic,
G. Meylan,
M. Tewes,
D. Sluse,
P. Magain,
J. H. H. Chan,
A. Galan,
R. Joseph,
C. Lemon,
O. Tihhonova,
R. I. Anderson,
M. Marmier,
B. Chazelas,
M. Lendl,
A. H. M. J. Triaud,
A. Wyttenbach
Abstract:
We present the results of 15 years of monitoring lensed quasars, which was conducted by the COSMOGRAIL programme at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed systems are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which nine reach a relative p…
▽ More
We present the results of 15 years of monitoring lensed quasars, which was conducted by the COSMOGRAIL programme at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed systems are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which nine reach a relative precision better than 15% for at least one time delay. To achieve this, we developed an automated version of the curve-shifting toolbox PyCS to ensure robust estimation of the time delay in the presence of microlensing, while accounting for the errors due to the imperfect representation of microlensing. We also re-analysed the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, by adding six and two new seasons of monitoring, respectively, and confirming the previous time-delay measurements. When the time delay measurement is possible, we corrected the light curves of the lensed images from their time delay and present the difference curves to highlight the microlensing signal contained in the data. To date, this is the largest sample of decade-long lens monitoring data, which is useful to measure $H_0$ and the size of quasar accretion discs with microlensing as well as to study quasar variability.
△ Less
Submitted 25 August, 2020; v1 submitted 13 February, 2020;
originally announced February 2020.
-
COSMOGRAIL XVII: Time delays for the quadruply imaged quasar PG 1115+080
Authors:
V. Bonvin,
J. H. H. Chan,
M. Millon,
K. Rojas,
F. Courbin,
G. C. -F. Chen,
C. D. Fassnacht,
E. Paic,
M. Tewes,
D. C. -Y. Chao,
M. Chijani,
D. Gilman,
K. Gilmore,
P. Williams,
E. Buckley-Geer,
J. Frieman,
P. J. Marshall,
S. H. Suyu,
T. Treu,
A. Hempel,
S. Kim,
R. Lachaume,
M. Rabus,
T. Anguita,
G. Meylan
, et al. (2 additional authors not shown)
Abstract:
We present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our resuls are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyse existing light curves from the literature that we complete with an additional three seasons of monitoring wi…
▽ More
We present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our resuls are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyse existing light curves from the literature that we complete with an additional three seasons of monitoring with the Mercator telescope at La Palma Observatory. When exploring the possible source of bias we consider the so-called microlensing time delay, a potential source of systematic error so far never directly accounted for in previous time-delay publications. In fifteen years of data on PG 1115+080, we find no strong evidence of microlensing time delay. Therefore not accounting for this effect, our time-delay estimates on the individual data sets are in good agreement with each other and with the literature. Combining the data sets, we obtain the most precise time-delay estimates to date on PG 1115+080, with Dt(AB) = 8.3+1.5-1.6 days (18.7% precision), Dt(AC) = 9.9+1.1-1.1 days (11.1%) and Dt(BC) = 18.8+1.6-1.6 days (8.5%). Turning these time delays into cosmological constraints is done in a companion paper that makes use of ground-based Adaptive Optics (AO) with the Keck telescope.
△ Less
Submitted 24 April, 2018;
originally announced April 2018.
-
COSMOGRAIL XVI: Time delays for the quadruply imaged quasar DES J0408-5354 with high-cadence photometric monitoring
Authors:
F. Courbin,
V. Bonvin,
E. Buckley-Geer,
C. D. Fassnacht,
J. Frieman,
H. Lin,
P. J. Marshall,
S. H. Suyu,
T. Treu,
T. Anguita,
V. Motta,
G. Meylan,
E. Paic,
M. Tewes,
A. Agnello,
D. C. -Y. Chao,
M. Chijani,
D. Gilman,
K. Rojas,
P. Williams,
A. Hempel,
S. Kim,
R. Lachaume,
M. Rabus,
T. M. C. Abbott
, et al. (52 additional authors not shown)
Abstract:
We present time-delay measurements for the new quadruply imaged quasar DES J0408-5354, the first quadruply imaged quasar found in the Dark Energy Survey (DES). Our result is made possible by implementing a new observational strategy using almost daily observations with the MPIA 2.2m telescope at La Silla observatory and deep exposures reaching a signal-to-noise ratio of about 1000 per quasar image…
▽ More
We present time-delay measurements for the new quadruply imaged quasar DES J0408-5354, the first quadruply imaged quasar found in the Dark Energy Survey (DES). Our result is made possible by implementing a new observational strategy using almost daily observations with the MPIA 2.2m telescope at La Silla observatory and deep exposures reaching a signal-to-noise ratio of about 1000 per quasar image. This data quality allows us to catch small photometric variations (a few mmag rms) of the quasar, acting on temporal scales much shorter than microlensing, hence making the time delay measurement very robust against microlensing. In only 7 months we measure very accurately one of the time delays in DES J0408-5354: Dt(AB) = -112.1 +- 2.1 days (1.8%) using only the MPIA 2.2m data. In combination with data taken with the 1.2m Euler Swiss telescope, we also measure two delays involving the D component of the system Dt(AD) = -155.5 +- 12.8 days (8.2%) and Dt(BD) = -42.4 +- 17.6 days (41%), where all the error bars include systematics. Turning these time delays into cosmological constraints will require deep HST imaging or ground-based Adaptive Optics (AO), and information on the velocity field of the lensing galaxy.
△ Less
Submitted 10 October, 2017; v1 submitted 28 June, 2017;
originally announced June 2017.