Showing 1–38 of 38 results for author: Schuldt, S

Lens Modeling of STRIDES Strongly Lensed Quasars using Neural Posterior Estimation

Authors: Sydney Erickson, Sebastian Wagner-Carena, Phil Marshall, Martin Millon, Simon Birrer, Aaron Roodman, Thomas Schmidt, Tommaso Treu, Stefan Schuldt, Anowar Shajib, Padma Venkatraman, The LSST Dark Energy Science Collaboration

Abstract: Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory's Legacy Survey of Space and Time (LSST) and the Euclid Wide Survey, we re… ▽ More Strongly lensed quasars can be used to constrain cosmological parameters through time-delay cosmography. Models of the lens masses are a necessary component of this analysis. To enable time-delay cosmography from a sample of $\mathcal{O}(10^3)$ lenses, which will soon become available from surveys like the Rubin Observatory's Legacy Survey of Space and Time (LSST) and the Euclid Wide Survey, we require fast and standardizable modeling techniques. To address this need, we apply neural posterior estimation (NPE) for modeling galaxy-scale strongly lensed quasars from the Strong Lensing Insights into the Dark Energy Survey (STRIDES) sample. NPE brings two advantages: speed and the ability to implicitly marginalize over nuisance parameters. We extend this method by employing sequential NPE to increase precision of mass model posteriors. We then fold individual lens models into a hierarchical Bayesian inference to recover the population distribution of lens mass parameters, accounting for out-of-distribution shift. After verifying our method using simulated analogs of the STRIDES lens sample, we apply our method to 14 Hubble Space Telescope single-filter observations. We find the population mean of the power-law elliptical mass distribution slope, $γ_{\text{lens}}$, to be $\mathcal{M}_{γ_{\text{lens}}}=2.13 \pm 0.06$. Our result represents the first population-level constraint for these systems. This population-level inference from fully automated modeling is an important stepping stone towards cosmological inference with large samples of strongly lensed quasars. △ Less

Submitted 13 October, 2024; originally announced October 2024.

The next step in galaxy cluster strong lensing: modeling the surface brightness of multiply-imaged sources

Authors: Ana Acebron, Claudio Grillo, Sherry H. Suyu, Giuseppe Angora, Pietro Bergamini, Gabriel B. Caminha, Sebastian Ertl, Amata Mercurio, Mario Nonino, Piero Rosati, Han Wang, Andrea Bolamperti, Massimo Meneghetti, Stefan Schuldt, Eros Vanzella

Abstract: Overcoming both modeling and computational challenges, we present, for the first time, the extended surface-brightness distribution model of a strongly-lensed source in a complex galaxy-cluster-scale system. We exploit the high-resolution Hubble Space Telescope (HST) imaging and extensive Multi Unit Spectroscopic Explorer spectroscopy to build an extended strong-lensing model, in a full multi-plan… ▽ More Overcoming both modeling and computational challenges, we present, for the first time, the extended surface-brightness distribution model of a strongly-lensed source in a complex galaxy-cluster-scale system. We exploit the high-resolution Hubble Space Telescope (HST) imaging and extensive Multi Unit Spectroscopic Explorer spectroscopy to build an extended strong-lensing model, in a full multi-plane formalism, of SDSS J1029+2623, a lens cluster at $z = 0.588$ with three multiple images of a background quasar ($z = 2.1992$). Going beyond typical cluster strong-lensing modeling techniques, we include as observables both the positions of 26 pointlike multiple images from seven background sources, spanning a wide redshift range between 1.02 and 5.06, and the extended surface-brightness distribution of the strongly-lensed quasar host galaxy, over $\sim78000$ HST pixels. In addition, we model the light distribution of seven objects, angularly close to the strongly-lensed quasar host, over $\sim9300$ HST pixels. Our extended lens model reproduces well both the observed intensity and morphology of the quasar host galaxy in the HST F160W band (with a 0''.03 pixel scale). The reconstructed source shows a single, compact, and smooth surface-brightness distribution, for which we estimate an intrinsic magnitude of 23.3 $\pm$ 0.1 in the F160W band and a half-light radius of (2.39 $\pm$ 0.03) kpc. The increased number of observables enables the accurate determination of the total mass of line-of-sight halos lying angularly close to the extended arc. This work paves the way for a new generation of galaxy cluster strong-lens models, where additional, complementary lensing observables are directly incorporated as model constraints. △ Less

Submitted 2 October, 2024; originally announced October 2024.

Comments: 19 pages, 9 Figures, 2 Tables. Accepted for publication in ApJ

Using Convolutional Neural Networks to Search for Strongly Lensed Quasars in KiDS DR5

Authors: Zizhao He, Rui Li, Yiping Shu, Crescenzo Tortora, Xinzhong Er, Raoul Canameras, Stefan Schuldt, Nicola R. Napolitano, Bharath Chowdhary N, Qihang Chen, Nan Li, Haicheng Feng, Limeng Deng, Guoliang Li, L. V. E. Koopmans, Andrej Dvornik

Abstract: Gravitationally strongly lensed quasars (SL-QSO) offer invaluable insights into cosmological and astrophysical phenomena. With the data from ongoing and next-generation surveys, thousands of SL-QSO systems can be discovered expectedly, leading to unprecedented opportunities. However, the challenge lies in identifying SL-QSO from enormous datasets with high recall and purity in an automated and eff… ▽ More Gravitationally strongly lensed quasars (SL-QSO) offer invaluable insights into cosmological and astrophysical phenomena. With the data from ongoing and next-generation surveys, thousands of SL-QSO systems can be discovered expectedly, leading to unprecedented opportunities. However, the challenge lies in identifying SL-QSO from enormous datasets with high recall and purity in an automated and efficient manner. Hence, we developed a program based on a Convolutional Neural Network (CNN) for finding SL-QSO from large-scale surveys and applied it to the Kilo-degree Survey Data Release 5 (KiDS DR5). Our approach involves three key stages: firstly, we pre-selected ten million bright objects (with $r$-band $\tt{MAG\_AUTO} < 22$), excluding stars from the dataset; secondly, we established realistic training and test sets to train and fine-tune the CNN, resulting in the identification of 4195 machine candidates, and the false positive rate (FPR) of $\sim$1/2000 and recall of 0.8125 evaluated by using the real test set containing 16 confirmed lensed quasars; thirdly, human inspections were performed for further selections, and then 272 SL-QSO candidates were eventually found in total, including 16 high-score, 118 median-score, and 138 lower-score candidates, separately. Removing the systems already confirmed or identified in other papers, we end up with 229 SL-QSO candidates, including 7 high-score, 95 median-score, and 127 lower-score candidates, and the corresponding catalog is publicly available online. △ Less

Submitted 25 September, 2024; originally announced September 2024.

Comments: 11 Figures, 4 Tables. Submitted to ApJ. Comments Welcome!

Euclid: The Early Release Observations Lens Search Experiment

Authors: J. A. Acevedo Barroso, C. M. O'Riordan, B. Clément, C. Tortora, T. E. Collett, F. Courbin, R. Gavazzi, R. B. Metcalf, V. Busillo, I. T. Andika, R. Cabanac, H. M. Courtois, J. Crook-Mansour, L. Delchambre, G. Despali, L. R. Ecker, A. Franco, P. Holloway, N. Jackson, K. Jahnke, G. Mahler, L. Marchetti, P. Matavulj, A. Melo, M. Meneghetti , et al. (182 additional authors not shown)

Abstract: We investigate the ability of the Euclid telescope to detect galaxy-scale gravitational lenses. To do so, we perform a systematic visual inspection of the $0.7\,\rm{deg}^2$ Euclid ERO data towards the Perseus cluster using both the high-resolution VIS $I_{\scriptscriptstyle\rm E}$ band, and the lower resolution NISP bands. We inspect every extended source brighter than magnitude $23$ in… ▽ More We investigate the ability of the Euclid telescope to detect galaxy-scale gravitational lenses. To do so, we perform a systematic visual inspection of the $0.7\,\rm{deg}^2$ Euclid ERO data towards the Perseus cluster using both the high-resolution VIS $I_{\scriptscriptstyle\rm E}$ band, and the lower resolution NISP bands. We inspect every extended source brighter than magnitude $23$ in $I_{\scriptscriptstyle\rm E}$ with $41$ expert human classifiers. This amounts to $12\,086$ stamps of $10^{\prime\prime}\,\times\,10^{\prime\prime}$. We find $3$ grade A and $13$ grade B candidates. We assess the validity of these $16$ candidates by modelling them and checking that they are consistent with a single source lensed by a plausible mass distribution. Five of the candidates pass this check, five others are rejected by the modelling and six are inconclusive. Extrapolating from the five successfully modelled candidates, we infer that the full $14\,000\,{\rm deg}^2$ of the Euclid Wide Survey should contain $100\,000^{+70\,000}_{-30\,000}$ galaxy-galaxy lenses that are both discoverable through visual inspection and have valid lens models. This is consistent with theoretical forecasts of $170\,000$ discoverable galaxy-galaxy lenses in Euclid. Our five modelled lenses have Einstein radii in the range $0.\!\!^{\prime\prime}68\,<\,θ_\mathrm{E}\,<1.\!\!^{\prime\prime}24$, but their Einstein radius distribution is on the higher side when compared to theoretical forecasts. This suggests that our methodology is likely missing small Einstein radius systems. Whilst it is implausible to visually inspect the full Euclid data set, our results corroborate the promise that Euclid will ultimately deliver a sample of around $10^5$ galaxy-scale lenses. △ Less

Submitted 12 August, 2024; originally announced August 2024.

Comments: 21 pages, 20 figures, submitted to A&A

Early-Time Observations of SN 2023wrk: A Luminous Type Ia Supernova with Significant Unburned Carbon in the Outer Ejecta

Authors: Jialian Liu, Xiaofeng Wang, Cristina Andrade, Pierre-Alexandre Duverne, Jujia Zhang, Liping Li, Zhenyu Wang, Felipe Navarete, Andrea Reguitti, Stefan Schuldt, Yongzhi Cai, Alexei V. Filippenko, Yi Yang, Thomas G. Brink, WeiKang Zheng, Ali Esamdin, Abdusamatjan Iskandar, Chunhai Bai, Jinzhong Liu, Xin Li, Maokai Hu, Gaici Li, Wenxiong Li, Xiaoran Ma, Shengyu Yan , et al. (22 additional authors not shown)

Abstract: We present extensive photometric and spectroscopic observations of the nearby Type Ia supernova (SN) 2023wrk at a distance of about 40 Mpc. The earliest detection of this SN can be traced back to a few hours after the explosion. Within the first few days the light curve shows a bump feature, while the B - V color is blue and remains nearly constant. The overall spectral evolution is similar to tha… ▽ More We present extensive photometric and spectroscopic observations of the nearby Type Ia supernova (SN) 2023wrk at a distance of about 40 Mpc. The earliest detection of this SN can be traced back to a few hours after the explosion. Within the first few days the light curve shows a bump feature, while the B - V color is blue and remains nearly constant. The overall spectral evolution is similar to that of an SN 1991T/SN 1999aa-like SN Ia, while the C II $\lambda6580$ absorption line appears to be unusually strong in the first spectrum taken at $t \approx -$15.4 days after the maximum light. This carbon feature disappears quickly in subsequent evolution but it reappears at around the time of peak brightness. The complex evolution of the carbon line and the possible detection of Ni III absorption around 4700 Å and 5300 Å in the earliest spectra indicate macroscopic mixing of fuel and ash. The strong carbon lines is likely related to collision of SN ejecta with unbound carbon, consistent with the predictions of pulsational delayed-detonation or carbon-rich circumstellar-matter interaction models. Among those carbon-rich SNe Ia with strong C II $\lambda6580$ absorption at very early times, the line-strength ratio of C II to Si II and the B-V color evolution are found to exhibit large diversity, which may be attributed to different properties of unbound carbon and outward-mixing $^{56}$Ni. △ Less

Submitted 22 July, 2024; originally announced July 2024.

Comments: Accepted for publication in the Astrophysical Journal (27 pages, 14 figures, 7 tables)

TDCOSMO. XVI. Measurement of the Hubble Constant from the Lensed Quasar WGD$\,$2038$-$4008

Authors: Kenneth C. Wong, Frédéric Dux, Anowar J. Shajib, Sherry H. Suyu, Martin Millon, Pritom Mozumdar, Patrick R. Wells, Adriano Agnello, Simon Birrer, Elizabeth J. Buckley-Geer, Frédéric Courbin, Christopher D. Fassnacht, Joshua Frieman, Aymeric Galan, Huan Lin, Philip J. Marshall, Jason Poh, Stefan Schuldt, Dominique Sluse, Tommaso Treu

Abstract: Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038-4008 using new time-delay measurements and previous mass models by TDCOSMO. This is… ▽ More Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant ($H_{0}$). The TDCOSMO collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038-4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain $D_{Δt} = 1.68^{+0.40}_{-0.38}$ Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat $Λ$CDM cosmology, we constrain the Hubble constant to be $H_{0} = 65^{+23}_{-14}\, \rm km\ s^{-1}\,Mpc^{-1}$. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory's Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future. △ Less

Submitted 5 July, 2024; v1 submitted 4 June, 2024; originally announced June 2024.

Comments: 8 pages, 5 figures, 3 tables; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 689, A168 (2024)

HOLISMOKES XIII: Strong-lens candidates at all mass scales and their environments from the Hyper-Suprime Cam and deep learning

Authors: Stefan Schuldt, Raoul Canameras, Irham T. Andika, Satadru Bag, Alejandra Melo, Yiping Shu, Sherry H. Suyu, Stefan Taubenberger, Claudio Grillo

Abstract: We have performed a systematic search for galaxy-scale strong lenses using Hyper Suprime-Cam imaging data, focusing on lenses in overdense environments. To identify these lens candidates, we exploit our neural network from HOLISMOKES VI, which is trained on realistic gri mock-images as positive examples, and real images as negative examples. Compared to our previous work, we lower the i-Kron radiu… ▽ More We have performed a systematic search for galaxy-scale strong lenses using Hyper Suprime-Cam imaging data, focusing on lenses in overdense environments. To identify these lens candidates, we exploit our neural network from HOLISMOKES VI, which is trained on realistic gri mock-images as positive examples, and real images as negative examples. Compared to our previous work, we lower the i-Kron radius limit to >0.5". This results in an increase by around 73 million sources to more than 135 million images. During our visual multi-stage grading of the network candidates, we now inspect simultaneously larger stamps (80"x80") to identify large, extended arcs cropped in the 10"x10" cutouts, and classify additionally their overall environment. Here we also reinspect our previous lens candidates and classify their environment. Using these 546 visually identified lens candidates, we further define various criteria by exploiting extensive and complementary photometric redshift catalogs, to select the candidates in overdensities. In total, we identified 24 grade-A and 138 grade-B candidates with either spatially-resolved multiple images or extended, distorted arcs in the new sample. Furthermore, with our different techniques, we identify in total 237/546 lens candidates in a cluster-like or overdense environment, containing only 49 group- or cluster-scale re-discoveries. These results demonstrate the feasibility of downloading and applying network classifiers to hundreds of million cutouts, necessary in the upcoming era of big data from deep, wide-field imaging surveys like Euclid and the Rubin Observatory Legacy Survey of Space and Time, while leading to a sample size that can be inspected by humans. These networks, with false-positive rates of ~0.01%, are very powerful tools to identify such rare galaxy-scale strong lensing systems, while also aiding in the discovery of new strong lensing clusters. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 12+11 pages, 5+9 figures, 3+1 Tables, submitted to A&A, comments are welcome

Systematic comparison of neural networks used in discovering strong gravitational lenses

Authors: Anupreeta More, Raoul Canameras, Anton T. Jaelani, Yiping Shu, Yuichiro Ishida, Kenneth C. Wong, Kaiki Taro Inoue, Stefan Schuldt, Alessandro Sonnenfeld

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

Submitted 21 May, 2024; originally announced May 2024.

Comments: 13 pages, 8 figures

Improved model of the Supernova Refsdal cluster MACS J1149.5+2223 thanks to VLT/MUSE

Authors: S. Schuldt, C. Grillo, G. B. Caminha, A. Mercurio, P. Rosati, T. Morishita, M. Stiavelli, S. H. Suyu, P. Bergamini, M. Brescia, F. Calura, M. Meneghetti

Abstract: We present new VLT/MUSE observations of the Hubble Frontier Field (HFF) galaxy cluster MACS J1149.5+2223, lensing the well-known supernova "Refsdal" into multiple images, which enabled the first cosmological applications with a strongly lensed supernova. Thanks to these data, targeting a northern region of the cluster and thus complementing our previous MUSE program on the cluster core, we release… ▽ More We present new VLT/MUSE observations of the Hubble Frontier Field (HFF) galaxy cluster MACS J1149.5+2223, lensing the well-known supernova "Refsdal" into multiple images, which enabled the first cosmological applications with a strongly lensed supernova. Thanks to these data, targeting a northern region of the cluster and thus complementing our previous MUSE program on the cluster core, we release a new catalog containing 162 secure spectroscopic redshifts. We confirm 22 cluster members, which were previously only photometrically selected, and detect ten additional ones, resulting in a total of 308 secure members, of which 63% are spectroscopically confirmed. We further identify 17 new spectroscopic multiple images belonging to 6 different background sources. By exploiting MUSE data, in combination with the deep HFF images, we develop an improved total mass model of MACS J1149.5+2223. This model includes 308 total mass components for the member galaxies and requires four additional mass profiles, one of which is associated with a cluster galaxy overdensity identified in the North, representing the DM mass distribution on larger scales. The values of the resulting 34 free parameters are optimized based on the observed positions of 106 multiple images from 34 different families, that cover the redshift range between 1.240 and 5.983. Our final model has a multiple image position rms value of 0.39", which is well in agreement with that of other cluster lens models. With this refined mass model, we pave the way towards even better strong-lensing analyses that will exploit the deep and high resolution observations with HST and JWST on a pixel level in the region of the supernova Refsdal host. This will increase the number of observables by around two orders of magnitudes, thus offering us the opportunity of carrying out more precise and accurate cosmographic measurements. △ Less

Submitted 4 September, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 14 pages, 7 figures, 4+1 tables, published with A&A

Detecting unresolved lensed SNe Ia in LSST using blended light curves

Authors: Satadru Bag, Simon Huber, Sherry H. Suyu, Nikki Arendse, Irham Taufik Andika, Raoul Canameras, Alex Kim, Eric Linder, Kushal Lodha, Alejandra Melo, Anupreeta More, Stefan Schuldt, Arman Shafieloo

Abstract: Strong-gravitationally lensed supernovae (LSNe) are promising probes for providing absolute distance measurements using gravitational lens time delays. Spatially unresolved LSNe offer an opportunity to enhance the sample size for precision cosmology. We predict that there will be approximately $3$ times more unresolved than resolved LSNe Ia in the Legacy Survey of Space and Time (LSST) by the Rubi… ▽ More Strong-gravitationally lensed supernovae (LSNe) are promising probes for providing absolute distance measurements using gravitational lens time delays. Spatially unresolved LSNe offer an opportunity to enhance the sample size for precision cosmology. We predict that there will be approximately $3$ times more unresolved than resolved LSNe Ia in the Legacy Survey of Space and Time (LSST) by the Rubin Observatory. In this article, we explore the feasibility of detecting unresolved LSNe Ia from the shape of the observed blended light curves using deep learning techniques, and we find that $\sim 30\%$ can be detected with a simple 1D CNN using well-sampled $rizy$-band light curves (with a false-positive rate of $\sim 3\%$). Even when the light curve is well-observed in only a single band among $r$, $i$, and $z$, detection is still possible with false-positive rates ranging from $\sim 4-7\%$, depending on the band. Furthermore, we demonstrate that these unresolved cases can be detected at an early stage using light curves up to $\sim20$ days from the first observation, with well-controlled false-positive rates, providing ample opportunities for triggering follow-up observations. Additionally, we demonstrate the feasibility of time-delay estimations using solely LSST-like data of unresolved light curves, particularly for doubles, when excluding systems with low time delay and magnification ratio. However, the abundance of such systems among those unresolved in LSST poses a significant challenge. This approach holds potential utility for upcoming wide-field surveys, and overall results could significantly improve with enhanced cadence and depth in the future surveys. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 16 pages, 9 figures, submitted to A&A

Dissecting the Interstellar Media of A Wolf-Rayet Galaxy at $z=2.76$

Authors: Takahiro Morishita, Massimo Stiavelli, Stefan Schuldt, Claudio Grillo

Abstract: We report JWST/NIRSpec observations of a star-forming galaxy at $z=2.76$, MACSJ1149-WR1. We securely detect two temperature-sensitive auroral lines, [SIII]6312 (7.4$σ$) and [OII]7320+7331 doublets (10$σ$), and tentatively [NII]5755 ($2.3σ$) for the first time in an individual galaxy at $z>1$. We perform a detailed analysis of its interstellar media (ISM), and derive electron temperatures, various… ▽ More We report JWST/NIRSpec observations of a star-forming galaxy at $z=2.76$, MACSJ1149-WR1. We securely detect two temperature-sensitive auroral lines, [SIII]6312 (7.4$σ$) and [OII]7320+7331 doublets (10$σ$), and tentatively [NII]5755 ($2.3σ$) for the first time in an individual galaxy at $z>1$. We perform a detailed analysis of its interstellar media (ISM), and derive electron temperatures, various heavy element abundances (O/H, N/O, S/O, and Ar/O) in the hot ionized region, and the neutral fraction in the warm ionized region. MACSJ1149-WR1 shows a broad feature at the wavelength of HeII 4686, which consists of a broad ($\sim1000$km/s), blue-shifted ($\sim-110$km/s) line component. Taken together with its mildly elevated N/O abundance, we conclude that MACSJ1149-WR1 is experiencing a young starburst ($<10$Myr), likely hosting a large number of Wolf-Rayet stars. None of its spectral features support the presence of AGN, including: $(i)$ the absence of broad components and velocity shifts in Hydrogen recombination lines, $(ii)$ low [FeII]${\rm 1.257 μm}$/Pa$β$ ratio, and $(iii)$ the absence of high-ionization lines. Our analysis using HeI lines reveals a higher electron temperature and a higher attenuation value, indicating that HeI may probe a smaller spatial scale than HI, presumably the region dominated by the aforementioned Wolf-Rayet stars. The star formation rates derived from various HeI lines broadly agree with those from Hydrogen recombination lines. We thus advocate that HeI can be an excellent, independent probe of multi-phase ISM in the era of JWST. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

Lensed Type Ia Supernova "Encore" at z=2: The First Instance of Two Multiply-Imaged Supernovae in the Same Host Galaxy

Authors: J. D. R. Pierel, A. B. Newman, S. Dhawan, M. Gu, B. A. Joshi, T. Li, S. Schuldt, L. G. Strolger, S. H. Suyu, G. B. Caminha, S. H. Cohen, J. M. Diego, J. C. J. Dsilva, S. Ertl, B. L. Frye, G. Granata, C. Grillo, A. M. Koekemoer, J. Li, A. Robotham, J. Summers, T. Treu, R. A. Windhorst, A. Zitrin, S. Agarwal , et al. (38 additional authors not shown)

Abstract: A bright ($m_{\rm F150W,AB}$=24 mag), $z=1.95$ supernova (SN) candidate was discovered in JWST/NIRCam imaging acquired on 2023 November 17. The SN is quintuply-imaged as a result of strong gravitational lensing by a foreground galaxy cluster, detected in three locations, and remarkably is the second lensed SN found in the same host galaxy. The previous lensed SN was called "Requiem", and therefore… ▽ More A bright ($m_{\rm F150W,AB}$=24 mag), $z=1.95$ supernova (SN) candidate was discovered in JWST/NIRCam imaging acquired on 2023 November 17. The SN is quintuply-imaged as a result of strong gravitational lensing by a foreground galaxy cluster, detected in three locations, and remarkably is the second lensed SN found in the same host galaxy. The previous lensed SN was called "Requiem", and therefore the new SN is named "Encore". This makes the MACS J0138.0$-$2155 cluster the first known system to produce more than one multiply-imaged SN. Moreover, both SN Requiem and SN Encore are Type Ia SNe (SNe Ia), making this the most distant case of a galaxy hosting two SNe Ia. Using parametric host fitting, we determine the probability of detecting two SNe Ia in this host galaxy over a $\sim10$ year window to be $\approx3\%$. These observations have the potential to yield a Hubble Constant ($H_0$) measurement with $\sim10\%$ precision, only the third lensed SN capable of such a result, using the three visible images of the SN. Both SN Requiem and SN Encore have a fourth image that is expected to appear within a few years of $\sim2030$, providing an unprecedented baseline for time-delay cosmography. △ Less

Submitted 22 July, 2024; v1 submitted 2 April, 2024; originally announced April 2024.

Comments: Accepted, ApJL

JWST Photometric Time-Delay and Magnification Measurements for the Triply-Imaged Type Ia "Supernova H0pe" at z = 1.78

Authors: J. D. R. Pierel, B. L. Frye, M. Pascale, G. B. Caminha, W. Chen, S. Dhawan, D. Gilman, M. Grayling, S. Huber, P. Kelly, S. Thorp, N. Arendse, S. Birrer, M. Bronikowski, R. Canameras, D. Coe, S. H. Cohen, C. J. Conselice, S. P. Driver, J. C. J. Dsilva, M. Engesser, N. Foo, C. Gall, N. Garuda, C. Grillo , et al. (38 additional authors not shown)

Abstract: Supernova (SN) H0pe is a gravitationally lensed, triply-imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply-imaged SNe provide a rare opportunity to constrain the Hubble constant ($H_0$), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is locate… ▽ More Supernova (SN) H0pe is a gravitationally lensed, triply-imaged, Type Ia SN (SN Ia) discovered in James Webb Space Telescope imaging of the PLCK G165.7+67.0 cluster of galaxies. Well-observed multiply-imaged SNe provide a rare opportunity to constrain the Hubble constant ($H_0$), by measuring the relative time delay between the images and modeling the foreground mass distribution. SN H0pe is located at $z=1.783$, and is the first SN Ia with sufficient light curve sampling and long enough time delays for an $H_0$ inference. Here we present photometric time-delay measurements and SN properties of SN H0pe. Using JWST/NIRCam photometry we measure time delays of $Δt_{ab}=-116.6^{+10.8}_{-9.3}$ and $Δt_{cb}=-48.6^{+3.6}_{-4.0}$ observer-frame days relative to the last image to arrive (image 2b; all uncertainties are $1σ$), which corresponds to a $\sim5.6\%$ uncertainty contribution for $H_0$ assuming $70 \rm{km s^{-1} Mpc^{-1}}$. We also constrain the absolute magnification of each image to $μ_{a}=4.3^{+1.6}_{-1.8}$, $μ_{b}=7.6^{+3.6}_{-2.6}$, $μ_{c}=6.4^{+1.6}_{-1.5}$ by comparing the observed peak near-IR magnitude of SN H0pe to the non-lensed population of SNe Ia. △ Less

Submitted 22 July, 2024; v1 submitted 27 March, 2024; originally announced March 2024.

Comments: Accepted in ApJ

Diverse Oxygen Abundance in Early Galaxies Unveiled by Auroral Line Analysis with JWST

Authors: Takahiro Morishita, Massimo Stiavelli, Claudio Grillo, Piero Rosati, Stefan Schuldt, Michele Trenti, Pietro Bergamini, Kristan N. Boyett, Ranga-Ram Chary, Nicha Leethochawalit, Guido Roberts-Borsani, Tommaso Treu, Eros Vanzella

Abstract: We present deep JWST NIRSpec observations in the sightline of MACS J1149.5+2223, a massive cluster of galaxies at $z=0.54$. We report the spectroscopic redshift of 28 sources at $3<z<9.1$, including 9 sources with the detection of the [OIII]4363 auroral line. Combining these with 16 [OIII]4363-detected sources from publicly available JWST data, our sample consists of 25 galaxies with robust gas-ph… ▽ More We present deep JWST NIRSpec observations in the sightline of MACS J1149.5+2223, a massive cluster of galaxies at $z=0.54$. We report the spectroscopic redshift of 28 sources at $3<z<9.1$, including 9 sources with the detection of the [OIII]4363 auroral line. Combining these with 16 [OIII]4363-detected sources from publicly available JWST data, our sample consists of 25 galaxies with robust gas-phase metallicity measurements via the direct method. We observe a positive correlation between stellar mass and metallicity, with a $\sim0.5$\,dex offset down below the local relation. Interestingly, we find a larger than expected scatter of $\sim0.3$\,dex around the relation, which cannot be explained by redshift evolution among our sample or other third parameter. The scatter increases at higher redshift, and we attribute this to the enrichment process having higher stochasticity due to shallower potential wells, more intense feedback processes, and a higher galaxy merger rate. Despite reaching to a considerably low-mass regime ($\log M_*/M_\odot \sim7.3$), our samples have metallicity of $\log$(O/H)$+12>7$, i.e. comparable to the most metal poor galaxies in the local Universe. The search of primordial galaxies may be accomplished by extending toward a lower mass and/or by investigating inhomogeneities at smaller spatial scales. Lastly, we investigate potential systematics caused by the limitation of JWST's MSA observations. Caution is warranted when the target exceeds the slit size, as this situation could allow an overestimation of "global" metallicity, especially under the presence of strong negative metallicity gradient. △ Less

Submitted 1 August, 2024; v1 submitted 21 February, 2024; originally announced February 2024.

Comments: Accepted for publication on May 29, 2024 to ApJ

Augmenting the power of time-delay cosmography in lens galaxy clusters by probing their member galaxies. II. Cosmic chronometers

Authors: P. Bergamini, S. Schuldt, A. Acebron, C. Grillo, U. Mestric, G. Granata, G. B. Caminha, M. Meneghetti, A. Mercurio, P. Rosati, S. H. Suyu, E. Vanzella

Abstract: We present a novel approach to measuring the expansion rate and the geometry of the Universe, which combine time-delay cosmography in lens galaxy clusters with pure samples of 'cosmic chronometers' (CCs) by probing the member galaxies. The former makes use of the measured time delays between the multiple images of time-varying sources strongly lensed by galaxy clusters, while the latter exploits t… ▽ More We present a novel approach to measuring the expansion rate and the geometry of the Universe, which combine time-delay cosmography in lens galaxy clusters with pure samples of 'cosmic chronometers' (CCs) by probing the member galaxies. The former makes use of the measured time delays between the multiple images of time-varying sources strongly lensed by galaxy clusters, while the latter exploits the most massive and passive cluster member galaxies to measure the differential time evolution of the Universe. We applied two different statistical techniques, adopting realistic errors on the measured quantities, to assess the accuracy and the gain in precision on the values of the cosmological parameters. We demonstrate that the proposed combined method allows for a robust and accurate measurement of the value of the Hubble constant. In addition, this provides valuable information on the other cosmological parameters thanks to the complementarity between the two different probes in breaking parameter degeneracies. Finally, we showcase the immediate observational feasibility of the proposed joint method by taking advantage of the existing high-quality spectro-photometric data for several lens galaxy clusters. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 6 pages, 2 figures. Accepted for publication in A&A Letters

Augmenting the power of time-delay cosmography in lens galaxy clusters by probing their member galaxies I. Type Ia supernovae

Authors: A. Acebron, S. Schuldt, C. Grillo, P. Bergamini, G. Granata, U. Mestric, G. B. Caminha, M. Meneghetti, A. Mercurio, P. Rosati, S. H. Suyu, E. Vanzella

Abstract: We present a simple and promising new method to measure the expansion rate and the geometry of the universe that combines observations related to the time delays between the multiple images of time-varying sources, strongly lensed by galaxy clusters, and Type Ia supernovae, exploding in galaxies belonging to the same lens clusters. By means of two different statistical techniques that adopt realis… ▽ More We present a simple and promising new method to measure the expansion rate and the geometry of the universe that combines observations related to the time delays between the multiple images of time-varying sources, strongly lensed by galaxy clusters, and Type Ia supernovae, exploding in galaxies belonging to the same lens clusters. By means of two different statistical techniques that adopt realistic errors on the relevant quantities, we quantify the accuracy of the inferred cosmological parameter values. We show that the estimate of the Hubble constant is robust and competitive, and depends only mildly on the chosen cosmological model. Remarkably, the two probes separately produce confidence regions on the cosmological parameter planes that are oriented in complementary ways, thus providing in combination valuable information on the values of the other cosmological parameters. We conclude by illustrating the immediate observational feasibility of the proposed joint method in a well-studied lens galaxy cluster, with a relatively small investment of telescope time for monitoring from a 2 to 3m class ground-based telescope. △ Less

Submitted 5 February, 2024; v1 submitted 6 November, 2023; originally announced November 2023.

Comments: 6 pages, 2 figures, 1 table. Published in A&A Letters

The puzzling properties of the MACS1149-JD1 galaxy at z=9.11

Authors: Massimo Stiavelli, Takahiro Morishita, Marco Chiaberge, Claudio Grillo, Nicha Leethochawalit, Piero Rosati, Stefan Schuldt, Michele Trenti, Tommaso Treu

Abstract: We analyze new JWST NIRCam and NIRSpec data on the redshift 9.11 galaxy MACS1149-JD1. Our NIRCam imaging data reveal that JD1 comprises three spatially distinct components. Our spectroscopic data indicate that JD1 appears dust-free but is already enriched, $12 + \log {\rm (O/H) } = 7.90^{+0.04}_{-0.05}$. We also find that the Carbon and Neon abundances in JD1 are below the solar abundance ratio. P… ▽ More We analyze new JWST NIRCam and NIRSpec data on the redshift 9.11 galaxy MACS1149-JD1. Our NIRCam imaging data reveal that JD1 comprises three spatially distinct components. Our spectroscopic data indicate that JD1 appears dust-free but is already enriched, $12 + \log {\rm (O/H) } = 7.90^{+0.04}_{-0.05}$. We also find that the Carbon and Neon abundances in JD1 are below the solar abundance ratio. Particularly the Carbon under-abundance is suggestive of recent star formation where Type~II supernovae have already enriched the ISM in Oxygen but intermediate mass stars have not yet enriched the ISM in Carbon. A recent burst of star formation is also revealed by the star formation history derived from NIRCam photometry. Our data do not reveal the presence of a significant amount of old populations, resulting in a factor of $\sim7\times$ smaller stellar mass than previous estimates. Thus, our data support the view that JD1 is a young galaxy. △ Less

Submitted 13 October, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: Accepted for publication

The missing quasar image in the gravitationally lensed quasar HE0230$-$2130: Implications for the cored lens mass distribution and dark satellites

Authors: S. Ertl, S. Schuldt, S. H. Suyu, P. L. Schechter, A. Halkola, J. Wagner

Abstract: Strongly lensed systems with peculiar configurations allow us to probe the local properties of the deflecting lens mass while simultaneously testing general profile assumptions. The quasar HE0230$-$2130 is lensed by two galaxies at similar redshifts ($Δz \sim 0.003$) into four observed images. Using modeled quasar positions from fitting the brightness of the quasar images in ground-based imaging d… ▽ More Strongly lensed systems with peculiar configurations allow us to probe the local properties of the deflecting lens mass while simultaneously testing general profile assumptions. The quasar HE0230$-$2130 is lensed by two galaxies at similar redshifts ($Δz \sim 0.003$) into four observed images. Using modeled quasar positions from fitting the brightness of the quasar images in ground-based imaging data from the Magellan telescope, we find that lens-mass models where each of these two galaxies is parametrized with a singular power-law (PL) profile predict five quasar images. To interpret the quad configuration of the system, we tested 12 different profile assumptions with the aim of obtaining lens-mass models that correctly predict only four observed images. We tested the effects of adopting: cored profiles for the lensing galaxies; external shear; and additional profiles to represent a dark matter clump. We find that half of our model classes can produce the correct image multiplicity. By comparing the Bayesian evidence of different model parametrizations, we favor two model classes: (i) one that incorporates two singular PL profiles for the lensing galaxies and a cored isothermal sphere in the region of the previously predicted fifth image (rNIS profile), and (ii) one with a bigger lensing galaxy parametrized by a singular PL profile and the smaller galaxy by a cored PL profile with external shear. We estimated the mass of the rNIS clump for each candidate model of our final Markov chain Monte Carlo sample, and find that only 2\% are in the range of $10^6 M_{\odot} \leq M_{\rm rNIS}\leq 10^9 M_{\odot}$, which is the predicted mass range of dark matter subhalos in cold dark matter simulations, or the mass of dark-matter-dominated and low-surface-brightness galaxies. We therefore favor the models with a cored mass distribution for the lens galaxy close to the predicted fifth image. △ Less

Submitted 29 February, 2024; v1 submitted 9 August, 2023; originally announced August 2023.

Comments: 22 pages, 13 figures, 7 tables; accepted by A&A

Streamlined Lensed Quasar Identification in Multiband Images via Ensemble Networks

Authors: Irham Taufik Andika, Sherry H. Suyu, Raoul Cañameras, Alejandra Melo, Stefan Schuldt, Yiping Shu, Anna-Christina Eilers, Anton Timur Jaelani, Minghao Yue

Abstract: Quasars experiencing strong lensing offer unique viewpoints on subjects related to the cosmic expansion rate, the dark matter profile within the foreground deflectors, and the quasar host galaxies. Unfortunately, identifying them in astronomical images is challenging since they are overwhelmed by the abundance of non-lenses. To address this, we have developed a novel approach by ensembling cutting… ▽ More Quasars experiencing strong lensing offer unique viewpoints on subjects related to the cosmic expansion rate, the dark matter profile within the foreground deflectors, and the quasar host galaxies. Unfortunately, identifying them in astronomical images is challenging since they are overwhelmed by the abundance of non-lenses. To address this, we have developed a novel approach by ensembling cutting-edge convolutional networks (CNNs) -- for instance, ResNet, Inception, NASNet, MobileNet, EfficientNet, and RegNet -- along with vision transformers (ViTs) trained on realistic galaxy-quasar lens simulations based on the Hyper Suprime-Cam (HSC) multiband images. While the individual model exhibits remarkable performance when evaluated against the test dataset, achieving an area under the receiver operating characteristic curve of $>$97.3% and a median false positive rate of 3.6%, it struggles to generalize in real data, indicated by numerous spurious sources picked by each classifier. A significant improvement is achieved by averaging these CNNs and ViTs, resulting in the impurities being downsized by factors up to 50. Subsequently, combining the HSC images with the UKIRT, VISTA, and unWISE data, we retrieve approximately 60 million sources as parent samples and reduce this to 892,609 after employing a photometry preselection to discover $z>1.5$ lensed quasars with Einstein radii of $θ_\mathrm{E}<5$ arcsec. Afterward, the ensemble classifier indicates 3080 sources with a high probability of being lenses, for which we visually inspect, yielding 210 prevailing candidates awaiting spectroscopic confirmation. These outcomes suggest that automated deep learning pipelines hold great potential in effectively detecting strong lenses in vast datasets with minimal manual visual inspection involved. △ Less

Submitted 18 August, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: Accepted for publication in the Astronomy & Astrophysics journal. 28 pages, 11 figures, and 3 tables. We welcome comments from the reader

Journal ref: A&A 678, A103 (2023)

HOLISMOKES -- XI. Evaluation of supervised neural networks for strong-lens searches in ground-based imaging surveys

Authors: R. Canameras, S. Schuldt, Y. Shu, S. H. Suyu, S. Taubenberger, I. T. Andika, S. Bag, K. T. Inoue, A. T. Jaelani, L. Leal-Taixe, T. Meinhardt, A. Melo, A. More

Abstract: While supervised neural networks have become state of the art for identifying the rare strong gravitational lenses from large imaging data sets, their selection remains significantly affected by the large number and diversity of nonlens contaminants. This work evaluates and compares systematically the performance of neural networks in order to move towards a rapid selection of galaxy-scale strong… ▽ More While supervised neural networks have become state of the art for identifying the rare strong gravitational lenses from large imaging data sets, their selection remains significantly affected by the large number and diversity of nonlens contaminants. This work evaluates and compares systematically the performance of neural networks in order to move towards a rapid selection of galaxy-scale strong lenses with minimal human input in the era of deep, wide-scale surveys. We used multiband images from PDR2 of the HSC Wide survey to build test sets mimicking an actual classification experiment, with 189 strong lenses previously found over the HSC footprint and 70,910 nonlens galaxies in COSMOS. Multiple networks were trained on different sets of realistic strong-lens simulations and nonlens galaxies, with various architectures and data pre-processing. The overall performances strongly depend on the construction of the ground-truth training data and they typically, but not systematically, improve using our baseline residual network architecture. Improvements are found when applying random shifts to the image centroids and square root stretches to the pixel values, adding z band, or using random viewpoints of the original images, but not when adding difference images to subtract emission from the central galaxy. The most significant gain is obtained with committees of networks trained on different data sets, and showing a moderate overlap between populations of false positives. Nearly-perfect invariance to image quality can be achieved by training networks either with large number of bands, or jointly with the PSF and science frames. Overall, we show the possibility to reach a TPR0 as high as 60% for the test sets under consideration, which opens promising perspectives for pure selection of strong lenses without human input using the Rubin Observatory and other forthcoming ground-based surveys. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 21 pages, 10 figures, submitted to A&A, comments are welcome

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.

Comments: 16 pages, 20 Figures

LensWatch: I. Resolved HST Observations and Constraints on the Strongly-Lensed Type Ia Supernova 2022qmx ("SN Zwicky")

Authors: J. D. R. Pierel, N. Arendse, S. Ertl, X. Huang, L. A. Moustakas, S. Schuldt, A. J. Shajib, Y. Shu, S. Birrer, M. Bronikowski, J. Hjorth, S. H. Suyu, S. Agarwal, A. Agnello, A. S. Bolton, S. Chakrabarti, C. Cold, F. Courbin, J. M. Della Costa, S. Dhawan, M. Engesser, O. D. Fox, C. Gall, S. Gomez, A. Goobar , et al. (17 additional authors not shown)

Abstract: Supernovae (SNe) that have been multiply-imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply-image… ▽ More Supernovae (SNe) that have been multiply-imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply-imaged Type Ia SN 2022qmx (aka, "SN Zwicky") at $z=0.3544$. SN Zwicky was discovered by the Zwicky Transient Facility (ZTF) in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multi-cycle "LensWatch" program. We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, F814W) but unresolved with WFC3/IR~F160W, and present an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between lens model predicted time delays ($\lesssim1$ day), and delays estimated with the single epoch of HST colors ($\lesssim3.5$ days), including the uncertainty from chromatic microlensing ($\sim1$-$1.5$ days). Our lens models converge to an Einstein radius of $θ_E=(0.168^{+0.009}_{-0.005})\prime\prime$, the smallest yet seen in a lensed SN system. The "standard candle" nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter than predicted by $\sim1.7^{+0.8}_{-0.6}$mag and $\sim0.9^{+0.8}_{-0.6}$mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models. △ Less

Submitted 22 July, 2024; v1 submitted 7 November, 2022; originally announced November 2022.

TDCOSMO X. Automated Modeling of 9 Strongly Lensed Quasars and Comparison Between Lens Modeling Software

Authors: S. Ertl, S. Schuldt, S. H. Suyu, T. Schmidt, T. Treu, S. Birrer, A. J. Shajib, D. Sluse

Abstract: To use strong gravitational lenses as an astrophysical or cosmological probe, models of their mass distributions are often needed. We present a new, time-efficient automation code for uniform modeling of strongly lensed quasars with GLEE, a lens modeling software, for high-resolution multi-band data. By using the observed positions of the lensed quasars and the spatially extended surface brightnes… ▽ More To use strong gravitational lenses as an astrophysical or cosmological probe, models of their mass distributions are often needed. We present a new, time-efficient automation code for uniform modeling of strongly lensed quasars with GLEE, a lens modeling software, for high-resolution multi-band data. By using the observed positions of the lensed quasars and the spatially extended surface brightness distribution of the lensed quasar host galaxy, we obtain a model of the mass distribution of the lens galaxy. We apply this uniform modeling pipeline to a sample of nine strongly lensed quasars with HST WFC 3 images. The models show in most cases well reconstructed light components and a good alignment between mass and light centroids. We find that the automated modeling code significantly reduces the user input time during the modeling process. The preparation time of required input files is reduced significantly. This automated modeling pipeline can efficiently produce uniform models of extensive lens system samples which can be used for further cosmological analysis. A blind test through a comparison with the results of an independent automated modeling pipeline based on the modeling software Lenstronomy reveals important lessons. Quantities such as Einstein radius, astrometry, mass flattening and position angle are generally robustly determined. Other quantities depend crucially on the quality of the data and the accuracy of the PSF reconstruction. Better data and/or more detailed analysis will be necessary to elevate our automated models to cosmography grade. Nevertheless, our pipeline enables the quick selection of lenses for follow-up monitoring and further modeling, significantly speeding up the construction of cosmography-grade models. This is an important step forward to take advantage of the orders of magnitude increase in the number of lenses expected in the coming decade. △ Less

Submitted 27 March, 2023; v1 submitted 7 September, 2022; originally announced September 2022.

Comments: 37 pages, 14 figures, published in A&A

Journal ref: A&A 672, A2 (2023)

HOLISMOKES -- X. Comparison between neural network and semi-automated traditional modeling of strong lenses

Authors: S. Schuldt, S. H. Suyu, R. Canameras, Y. Shu, S. Taubenberger, S. Ertl, A. Halkola

Abstract: Modeling of strongly gravitationally lensed galaxies is often required in order to use them as astrophysical or cosmological probes. With current and upcoming wide-field imaging surveys, the number of detected lenses is increasing significantly such that automated and fast modeling procedures for ground-based data are urgently needed. This is especially pertinent to short-lived lensed transients i… ▽ More Modeling of strongly gravitationally lensed galaxies is often required in order to use them as astrophysical or cosmological probes. With current and upcoming wide-field imaging surveys, the number of detected lenses is increasing significantly such that automated and fast modeling procedures for ground-based data are urgently needed. This is especially pertinent to short-lived lensed transients in order to plan follow-up observations. Therefore, we present in a companion paper (submitted) a neural network predicting the parameter values with corresponding uncertainties of a Singular Isothermal Ellipsoid (SIE) mass profile with external shear. In this work, we present a newly-developed pipeline glee_auto.py to model consistently any galaxy-scale lensing system. In contrast to previous automated modeling pipelines that require high-resolution images, glee_auto.py is optimized for ground-based images such as those from the Hyper-Suprime-Cam (HSC) or the upcoming Rubin Observatory Legacy Survey of Space and Time. We further present glee_tools.py, a flexible automation code for individual modeling that has no direct decisions and assumptions implemented. Both pipelines, in addition to our modeling network, minimize the user input time drastically and thus are important for future modeling efforts. We apply the network to 31 real galaxy-scale lenses of HSC and compare the results to the traditional models. In the direct comparison, we find a very good match for the Einstein radius especially for systems with $θ_E \gtrsim 2$". The lens mass center and ellipticity show reasonable agreement. The main discrepancies are on the external shear as expected from our tests on mock systems. In general, our study demonstrates that neural networks are a viable and ultra fast approach for measuring the lens-galaxy masses from ground-based data in the upcoming era with $\sim10^5$ lenses expected. △ Less

Submitted 29 March, 2023; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: 15+45 pages, 7+31 figures, 1+6 tables, published with A&A

Journal ref: A&A 673, A33 (2023)

HOLISMOKES -- IX. Neural network inference of strong-lens parameters and uncertainties from ground-based images

Authors: S. Schuldt, R. Cañameras, Y. Shu, S. H. Suyu, S. Taubenberger, T. Meinhardt, L. Leal-Taixé

Abstract: Modeling of strong gravitational lenses is a necessity for further applications in astrophysics and cosmology. Especially with the large number of detections in current and upcoming surveys such as the Rubin Legacy Survey of Space and Time (LSST), it is timely to investigate in automated and fast analysis techniques beyond the traditional and time consuming Markov chain Monte Carlo sampling method… ▽ More Modeling of strong gravitational lenses is a necessity for further applications in astrophysics and cosmology. Especially with the large number of detections in current and upcoming surveys such as the Rubin Legacy Survey of Space and Time (LSST), it is timely to investigate in automated and fast analysis techniques beyond the traditional and time consuming Markov chain Monte Carlo sampling methods. Building upon our convolutional neural network (CNN) presented in Schuldt et al. (2021b), we present here another CNN, specifically a residual neural network (ResNet), that predicts the five mass parameters of a Singular Isothermal Ellipsoid (SIE) profile (lens center $x$ and $y$, ellipticity $e_x$ and $e_y$, Einstein radius $θ_E$) and the external shear ($γ_{ext,1}$, $γ_{ext,2}$) from ground-based imaging data. In contrast to our CNN, this ResNet further predicts a 1$σ$ uncertainty for each parameter. To train our network, we use our improved pipeline from Schuldt et al. (2021b) to simulate lens images using real images of galaxies from the Hyper Suprime-Cam Survey (HSC) and from the Hubble Ultra Deep Field as lens galaxies and background sources, respectively. We find overall very good recoveries for the SIE parameters, while differences remain in predicting the external shear. From our tests, most likely the low image resolution is the limiting factor for predicting the external shear. Given the run time of milli-seconds per system, our network is perfectly suited to predict the next appearing image and time delays of lensed transients in time. Therefore, we also present the performance of the network on these quantities in comparison to our simulations. Our ResNet is able to predict the SIE and shear parameter values in fractions of a second on a single CPU such that we are able to process efficiently the huge amount of expected galaxy-scale lenses in the near future. △ Less

Submitted 29 March, 2023; v1 submitted 22 June, 2022; originally announced June 2022.

Comments: 17 pages, including 11 figures, published with A&A

Journal ref: A&A 671, A147 (2023)

STRIDES: Automated uniform models for 30 quadruply imaged quasars

Authors: T. Schmidt, T. Treu, S. Birrer, A. J. Shajib, C. Lemon, M. Millon, D. Sluse, A. Agnello, T. Anguita, M. W. Auger-Williams, R. G. McMahon, V. Motta, P. Schechter, C. Spiniello, I. Kayo, F. Courbin, S. Ertl, C. D. Fassnacht, J. A. Frieman, A. More, S. Schuldt, S. H. Suyu, M. Aguena, F. Andrade-Oliveira, J. Annis , et al. (46 additional authors not shown)

Abstract: Gravitational time delays provide a powerful one step measurement of $H_0$, independent of all other probes. One key ingredient in time delay cosmography are high accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (10$^{5-6}$ CPU hours and $\sim$ 0.5-1 year, respectively). Major improvements in modeling speed are therefore necessary to… ▽ More Gravitational time delays provide a powerful one step measurement of $H_0$, independent of all other probes. One key ingredient in time delay cosmography are high accuracy lens models. Those are currently expensive to obtain, both, in terms of computing and investigator time (10$^{5-6}$ CPU hours and $\sim$ 0.5-1 year, respectively). Major improvements in modeling speed are therefore necessary to exploit the large number of lenses that are forecast to be discovered over the current decade. In order to bypass this roadblock, building on the work by Shajib et al. (2019), we develop an automated modeling pipeline and apply it to a sample of 30 quadruply imaged quasars and one lensed compact galaxy, observed by the Hubble Space Telescope in multiple bands. Our automated pipeline can derive models for 30/31 lenses with few hours of human time and <100 CPU hours of computing time for a typical system. For each lens, we provide measurements of key parameters and predictions of magnification as well as time delays for the multiple images. We characterize the cosmography-readiness of our models using the stability of differences in Fermat potential (proportional to time delay) w.r.t. modeling choices. We find that for 10/30 lenses our models are cosmography or nearly cosmography grade (<3% and 3-5% variations). For 6/30 lenses the models are close to cosmography grade (5-10%). These results are based on informative priors and will need to be confirmed by further analysis. However, they are also likely to improve by extending the pipeline modeling sequence and options. In conclusion, we show that uniform cosmography grade modeling of large strong lens samples is within reach. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 40 pages, 24 figures, 11 tables. Submitted to MNRAS

Observations of the Very Young Type Ia Supernova 2019np with Early-excess Emission

Authors: Hanna Sai, Xiaofeng Wang, Nancy Elias-Rosa, Yi Yang, Jujia Zhang, Weili Lin, Jun Mo, Anthony L. Piro, Xiangyun Zeng, Reguitti Andrea, Peter Brown, Christopher R. Burns, Yongzhi Cai, Achille Fiore, Eric Y. Hsiao, Jordi Isern, K. Itagaki, Wenxiong Li, Zhitong Li, Priscila J. Pessi, M. M. Phillips, Stefan Schuldt, Melissa Shahbandeh, Maximilian D. Stritzinger, Lina Tomasella , et al. (8 additional authors not shown)

Abstract: Early-time radiative signals from type Ia supernovae (SNe Ia) can provide important constraints on the explosion mechanism and the progenitor system. We present observations and analysis of SN 2019np, a nearby SN Ia discovered within 1-2 days after the explosion. Follow-up observations were conducted in optical, ultraviolet, and near-infrared bands, covering the phases from $\sim-$16.7 days to… ▽ More Early-time radiative signals from type Ia supernovae (SNe Ia) can provide important constraints on the explosion mechanism and the progenitor system. We present observations and analysis of SN 2019np, a nearby SN Ia discovered within 1-2 days after the explosion. Follow-up observations were conducted in optical, ultraviolet, and near-infrared bands, covering the phases from $\sim-$16.7 days to $\sim$+367.8 days relative to its $B-$band peak luminosity. The photometric and spectral evolutions of SN 2019np resembles the average behavior of normal SNe Ia. The absolute B-band peak magnitude and the post-peak decline rate are $M_{\rm max}(B)=-19.52 \pm 0.47$mag and $Δm_{\rm15}(B) =1.04 \pm 0.04$mag, respectively. No Hydrogen line has been detected in the near-infrared and nebular-phase spectra of SN 2019np. Assuming that the $^{56}$Ni powering the light curve is centrally located, we find that the bolometric light curve of SN 2019np shows a flux excess up to 5.0% in the early phase compared to the radiative diffusion model. Such an extra radiation perhaps suggests the presence of an additional energy source beyond the radioactive decay of central nickel. Comparing the observed color evolution with that predicted by different models such as interactions of SN ejecta with circumstellar matter (CSM)/companion star, a double-detonation explosion from a sub-Chandrasekhar mass white dwarf (WD), and surface $^{56}$Ni mixing, the latter one is favored. △ Less

Submitted 1 June, 2022; v1 submitted 31 May, 2022; originally announced May 2022.

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

Constraining the multi-scale dark-matter distribution in CASSOWARY 31 with strong gravitational lensing and stellar dynamics

Authors: H. Wang, R. Cañameras, G. B. Caminha, S. H. Suyu, A. Yıldırım, G. Chirivì, L. Christensen, C. Grillo, S. Schuldt

Abstract: We study the inner structure of the group-scale lens CASSOWARY 31 (CSWA 31) by adopting both strong lensing and dynamical modeling. CSWA 31 is a peculiar lens system. The brightest group galaxy (BGG) is an ultra-massive elliptical galaxy at z = 0.683 with a weighted mean velocity dispersion of $σ= 432 \pm 31$ km s$^{-1}$. It is surrounded by group members and several lensed arcs probing up to ~150… ▽ More We study the inner structure of the group-scale lens CASSOWARY 31 (CSWA 31) by adopting both strong lensing and dynamical modeling. CSWA 31 is a peculiar lens system. The brightest group galaxy (BGG) is an ultra-massive elliptical galaxy at z = 0.683 with a weighted mean velocity dispersion of $σ= 432 \pm 31$ km s$^{-1}$. It is surrounded by group members and several lensed arcs probing up to ~150 kpc in projection. Our results significantly improve previous analyses of CSWA 31 thanks to the new HST imaging and MUSE integral-field spectroscopy. From the secure identification of five sets of multiple images and measurements of the spatially-resolved stellar kinematics of the BGG, we conduct a detailed analysis of the multi-scale mass distribution using various modeling approaches, both in the single and multiple lens-plane scenarios. Our best-fit mass models reproduce the positions of multiple images and provide robust reconstructions for two background galaxies at z = 1.4869 and z = 2.763. The relative contributions from the BGG and group-scale halo are remarkably consistent in our three reference models, demonstrating the self-consistency between strong lensing analyses based on image position and extended image modeling. We find that the ultra-massive BGG dominates the projected total mass profiles within 20 kpc, while the group-scale halo dominates at larger radii. The total projected mass enclosed within $R_{eff}$ = 27.2 kpc is $1.10_{-0.04}^{+0.02} \times 10^{13}$ M$_\odot$. We find that CSWA 31 is a peculiar fossil group, strongly dark-matter dominated towards the central region, and with a projected total mass profile similar to higher-mass cluster-scale halos. The total mass-density slope within the effective radius is shallower than isothermal, consistent with previous analyses of early-type galaxies in overdense environments. △ Less

Submitted 20 October, 2022; v1 submitted 25 March, 2022; originally announced March 2022.

Comments: 22 pages, 12 figures, 5 tables, submitted to Astronomy & Astrophysics. We welcome the comments from readers

TDCOSMO. IX. Systematic comparison between lens modelling software programs: time delay prediction for WGD 2038-4008

Authors: A. J. Shajib, K. C. Wong, S. Birrer, S. H. Suyu, T. Treu, E. Buckley-Geer, H. Lin, C. E. Rusu, J. Poh, A. Palmese, A. Agnello, M. W. Auger-Williams, A. Galan, S. Schuldt, D. Sluse, F. Courbin, J. Frieman, M. Millon

Abstract: The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedur… ▽ More The importance of alternative methods for measuring the Hubble constant, such as time-delay cosmography, is highlighted by the recent Hubble tension. It is paramount to thoroughly investigate and rule out systematic biases in all measurement methods before we can accept new physics as the source of this tension. In this study, we perform a check for systematic biases in the lens modelling procedure of time-delay cosmography by comparing independent and blind time-delay predictions of the system WGD 2038-4008 from two teams using two different software programs: GLEE and LENSTRONOMY. The predicted time delays from the two teams incorporate the stellar kinematics of the deflector and the external convergence from line-of-sight structures. The un-blinded time-delay predictions from the two teams agree within $1.2σ$, implying that once the time delay is measured the inferred Hubble constant will also be mutually consistent. However, there is a $\sim$4$σ$ discrepancy between the power-law model slope and external shear, which is a significant discrepancy at the level of lens models before the stellar kinematics and the external convergence are incorporated. We identify the difference in the reconstructed point spread function (PSF) to be the source of this discrepancy. When the same reconstructed PSF was used by both teams, we achieved excellent agreement, within $\sim$0.6$σ$, indicating that potential systematics stemming from source reconstruction algorithms and investigator choices are well under control. We recommend that future studies supersample the PSF as needed and marginalize over multiple algorithms or realizations for the PSF reconstruction to mitigate the systematics associated with the PSF. A future study will measure the time delays of the system WGD 2038-4008 and infer the Hubble constant based on our mass models. △ Less

Submitted 21 September, 2022; v1 submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted version by A&A. 34 pages, 26 figures, 3 tables

Journal ref: A&A 667, A123 (2022)

HOLISMOKES. VIII. High-redshift, strong-lens search in the Hyper Suprime-Cam Subaru Strategic Program

Authors: Yiping Shu, Raoul Cañameras, Stefan Schuldt, Sherry H. Suyu, Stefan Taubenberger, Kaiki Taro Inoue, Anton T. Jaelani

Abstract: We carry out a search for strong-lens systems containing high-redshift lens galaxies with the goal of extending strong-lensing-assisted galaxy evolutionary studies to earlier cosmic time. Two strong-lens classifiers are constructed from a deep residual network and trained with datasets of different lens-redshift and brightness distributions. We classify a sample of 5,356,628 pre-selected objects f… ▽ More We carry out a search for strong-lens systems containing high-redshift lens galaxies with the goal of extending strong-lensing-assisted galaxy evolutionary studies to earlier cosmic time. Two strong-lens classifiers are constructed from a deep residual network and trained with datasets of different lens-redshift and brightness distributions. We classify a sample of 5,356,628 pre-selected objects from the Wide-layer fields in the second public data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) by applying the two classifiers to their HSC $gri$-filter cutouts. Cutting off at thresholds that correspond to a false positive rate of $10^{-3}$ on our test set, the two classifiers identify 5,468 and 6,119 strong-lens candidates. Visually inspecting the cutouts of those candidates results in 735 grade-A or B strong-lens candidates in total, of which 277 candidates are discovered for the first time. This is the single largest set of galaxy-scale strong-lens candidates discovered with HSC data to date, and nearly half of it (331/735) contains lens galaxies with photometric redshifts above 0.6. Our discoveries will serve as a valuable target list for ongoing and scheduled spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, the Subaru Prime Focus Spectrograph project, and the Maunakea Spectroscopic Explorer. △ Less

Submitted 25 May, 2022; v1 submitted 26 January, 2022; originally announced January 2022.

Comments: 15 pages, 8 figures + appendices, published on the A&A

Report number: 2022, A&A, 662, A4

Journal ref: A&A 662, A4 (2022)

Close, bright and boxy: the superluminous SN 2018hti

Authors: A. Fiore, S. Benetti, M. Nicholl, A. Reguitti, E. Cappellaro, S. Campana, S. Bose, E. Paraskeva, E. Berger, T. M. Bravo, J. Burke, Y. -Z. Cai, T. -W. Chen, P. Chen, R. Ciolfi, S. Dong, S. Gomez, M. Gromadzki, C. P. Gutiérrez, D. Hiramatsu, G. Hosseinzadeh, D. A. Howell, A. Jerkstrand, E. Kankare, A. Kozyreva , et al. (15 additional authors not shown)

Abstract: SN 2018hti was a very nearby (z=0.0614) superluminous supernova with an exceedingly bright absolute magnitude of -21.7 mag in r-band at maximum. The densely sampled pre-maximum light curves of SN 2018hti show a slow luminosity evolution and constrain the rise time to ~50 rest-frame days. We fitted synthetic light curves to the photometry to infer the physical parameters of the explosion of SN 2018… ▽ More SN 2018hti was a very nearby (z=0.0614) superluminous supernova with an exceedingly bright absolute magnitude of -21.7 mag in r-band at maximum. The densely sampled pre-maximum light curves of SN 2018hti show a slow luminosity evolution and constrain the rise time to ~50 rest-frame days. We fitted synthetic light curves to the photometry to infer the physical parameters of the explosion of SN 2018hti for both the magnetar and the CSM-interaction scenarios. We conclude that one of two mechanisms could be powering the luminosity of SN 2018hti; interaction with ~10 Msun of circumstellar material or a magnetar with a magnetic field of B_p~1.3e13 G and initial period of P_spin~1.8 ms. From the nebular spectrum modelling we infer that SN 2018hti likely results from the explosion of a ~40 Msun progenitor star. △ Less

Submitted 15 March, 2022; v1 submitted 13 November, 2021; originally announced November 2021.

Comments: 31 pages, 19 figures, replaced after acceptance by MNRAS (minor revisions compared to the previous version)

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)

HOLISMOKES. VI. New galaxy-scale strong lens candidates from the HSC-SSP imaging survey

Authors: R. Canameras, S. Schuldt, Y. Shu, S. H. Suyu, S. Taubenberger, T. Meinhardt, L. Leal-Taixé, D. C. -Y. Chao, K. T. Inoue, A. T. Jaelani, A. More

Abstract: We have carried out a systematic search for galaxy-scale strong lenses in multiband imaging from the Hyper Suprime-Cam (HSC) survey. Our automated pipeline, based on realistic strong-lens simulations, deep neural network classification, and visual inspection, is aimed at efficiently selecting systems with wide image separations (Einstein radii ~1.0-3.0"), intermediate redshift lenses (z ~ 0.4-0.7)… ▽ More We have carried out a systematic search for galaxy-scale strong lenses in multiband imaging from the Hyper Suprime-Cam (HSC) survey. Our automated pipeline, based on realistic strong-lens simulations, deep neural network classification, and visual inspection, is aimed at efficiently selecting systems with wide image separations (Einstein radii ~1.0-3.0"), intermediate redshift lenses (z ~ 0.4-0.7), and bright arcs for galaxy evolution and cosmology. We classified gri images of all 62.5 million galaxies in HSC Wide with i-band Kron radius >0.8" to avoid strict pre-selections and to prepare for the upcoming era of deep, wide-scale imaging surveys with Euclid and Rubin Observatory. We obtained 206 newly-discovered candidates classified as definite or probable lenses with either spatially-resolved multiple images or extended, distorted arcs. In addition, we found 88 high-quality candidates that were assigned lower confidence in previous HSC searches, and we recovered 173 known systems in the literature. These results demonstrate that, aided by limited human input, deep learning pipelines with false positive rates as low as ~0.01% can be very powerful tools for identifying the rare strong lenses from large catalogs, and can also largely extend the samples found by traditional algorithms. We provide a ranked list of candidates for future spectroscopic confirmation. △ Less

Submitted 7 September, 2021; v1 submitted 16 July, 2021; originally announced July 2021.

Comments: 5 pages and 4 figures (plus appendix), version published in A&A

Journal ref: A&A 653, L6 (2021)

Photometric Redshift Estimation with a Convolutional Neural Network: NetZ

Authors: S. Schuldt, S. H. Suyu, R. Cañameras, S. Taubenberger, T. Meinhardt, L. Leal-Taixé, B. C. Hsieh

Abstract: The redshifts of galaxies are a key attribute that is needed for nearly all extragalactic studies. Since spectroscopic redshifts require additional telescope and human resources, millions of galaxies are known without spectroscopic redshifts. Therefore, it is crucial to have methods for estimating the redshift of a galaxy based on its photometric properties, the so-called photo-$z$. We developed N… ▽ More The redshifts of galaxies are a key attribute that is needed for nearly all extragalactic studies. Since spectroscopic redshifts require additional telescope and human resources, millions of galaxies are known without spectroscopic redshifts. Therefore, it is crucial to have methods for estimating the redshift of a galaxy based on its photometric properties, the so-called photo-$z$. We developed NetZ, a new method using a Convolutional Neural Network (CNN) to predict the photo-$z$ based on galaxy images, in contrast to previous methods which often used only the integrated photometry of galaxies without their images. We use data from the Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) in five different filters as training data. The network over the whole redshift range between 0 and 4 performs well overall and especially in the high-$z$ range better than other methods on the same data. We obtain an accuracy $|z_\text{pred}-z_\text{ref}|$ of $σ= 0.12$ (68% confidence interval) with a CNN working for all galaxy types averaged over all galaxies in the redshift range of 0 to $\sim$4. By limiting to smaller redshift ranges or to Luminous Red Galaxies (LRGs), we find a further notable improvement. We publish more than 34 million new photo-$z$ values predicted with NetZ here. This shows that the new method is very simple and fast to apply, and, importantly, covers a wide redshift range limited only by the available training data. It is broadly applicable and beneficial to imaging surveys, particularly upcoming surveys like the Rubin Observatory Legacy Survey of Space and Time which will provide images of billions of galaxies with similar image quality as HSC. △ Less

Submitted 15 July, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: 10 pages, 11 figures, 1 catalog, published in A&A

Journal ref: A&A 651, A55 (2021)

HOLISMOKES -- IV. Efficient mass modeling of strong lenses through deep learning

Authors: S. Schuldt, S. H. Suyu, T. Meinhardt, L. Leal-Taixé, R. Cañameras, S. Taubenberger, A. Halkola

Abstract: Modelling the mass distributions of strong gravitational lenses is often necessary to use them as astrophysical and cosmological probes. With the high number of lens systems ($>10^5$) expected from upcoming surveys, it is timely to explore efficient modeling approaches beyond traditional MCMC techniques that are time consuming. We train a CNN on images of galaxy-scale lenses to predict the paramet… ▽ More Modelling the mass distributions of strong gravitational lenses is often necessary to use them as astrophysical and cosmological probes. With the high number of lens systems ($>10^5$) expected from upcoming surveys, it is timely to explore efficient modeling approaches beyond traditional MCMC techniques that are time consuming. We train a CNN on images of galaxy-scale lenses to predict the parameters of the SIE mass model ($x,y,e_x,e_y$, and $θ_E$). To train the network, we simulate images based on real observations from the HSC Survey for the lens galaxies and from the HUDF as lensed galaxies. We tested different network architectures, the effect of different data sets, and using different input distributions of $θ_E$. We find that the CNN performs well and obtain with the network trained with a uniform distribution of $θ_E$ $>0.5"$ the following median values with $1σ$ scatter: $Δx=(0.00^{+0.30}_{-0.30})"$, $Δy=(0.00^{+0.30}_{-0.29})" $, $Δθ_E=(0.07^{+0.29}_{-0.12})"$, $Δe_x = -0.01^{+0.08}_{-0.09}$ and $Δe_y = 0.00^{+0.08}_{-0.09}$. 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 $x$ and $y$, respectively. For the fractional difference between the predicted and true time delay, we obtain $0.04_{-0.05}^{+0.27}$. Our CNN is able to predict the SIE parameters 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 lens detections in the near future. △ Less

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

Comments: 17 pages, 14 Figures

HOLISMOKES -- II. Identifying galaxy-scale strong gravitational lenses in Pan-STARRS using convolutional neural networks

Authors: R. Canameras, S. Schuldt, S. H. Suyu, S. Taubenberger, T. Meinhardt, L. Leal-Taixe, C. Lemon, K. Rojas, E. Savary

Abstract: We present a systematic search for wide-separation (Einstein radius >1.5"), galaxy-scale strong lenses in the 30 000 sq.deg of the Pan-STARRS 3pi survey on the Northern sky. With long time delays of a few days to weeks, such systems are particularly well suited for catching strongly lensed supernovae with spatially-resolved multiple images and open new perspectives on early-phase supernova spectro… ▽ More We present a systematic search for wide-separation (Einstein radius >1.5"), galaxy-scale strong lenses in the 30 000 sq.deg of the Pan-STARRS 3pi survey on the Northern sky. With long time delays of a few days to weeks, such systems are particularly well suited for catching strongly lensed supernovae with spatially-resolved multiple images and open new perspectives on early-phase supernova spectroscopy and cosmography. We produce a set of realistic simulations by painting lensed COSMOS sources on Pan-STARRS image cutouts of lens luminous red galaxies with known redshift and velocity dispersion from SDSS. First of all, we compute the photometry of mock lenses in gri bands and apply a simple catalog-level neural network to identify a sample of 1050207 galaxies with similar colors and magnitudes as the mocks. Secondly, we train a convolutional neural network (CNN) on Pan-STARRS gri image cutouts to classify this sample and obtain sets of 105760 and 12382 lens candidates with scores pCNN>0.5 and >0.9, respectively. Extensive tests show that CNN performances rely heavily on the design of lens simulations and choice of negative examples for training, but little on the network architecture. Finally, we visually inspect all galaxies with pCNN>0.9 to assemble a final set of 330 high-quality newly-discovered lens candidates while recovering 23 published systems. For a subset, SDSS spectroscopy on the lens central regions proves our method correctly identifies lens LRGs at z~0.1-0.7. Five spectra also show robust signatures of high-redshift background sources and Pan-STARRS imaging confirms one of them as a quadruply-imaged red source at z_s = 1.185 strongly lensed by a foreground LRG at z_d = 0.3155. In the future, we expect that the efficient and automated two-step classification method presented in this paper will be applicable to the deeper gri stacks from the LSST with minor adjustments. △ Less

Submitted 7 April, 2021; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: 18 pages and 11 figures (plus appendix), version published in A&A

Journal ref: A&A 644, A163 (2020)

HOLISMOKES -- I. Highly Optimised Lensing Investigations of Supernovae, Microlensing Objects, and Kinematics of Ellipticals and Spirals

Authors: S. H. Suyu, S. Huber, R. Cañameras, M. Kromer, S. Schuldt, S. Taubenberger, A. Yıldırım, V. Bonvin, J. H. H. Chan, F. Courbin, U. Nöbauer, S. A. Sim, D. Sluse

Abstract: We present the HOLISMOKES programme on strong gravitational lensing of supernovae as a probe of supernova (SN) physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra using four different SN explosion models, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible ($<$1% distortion… ▽ More We present the HOLISMOKES programme on strong gravitational lensing of supernovae as a probe of supernova (SN) physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra using four different SN explosion models, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible ($<$1% distortion within the 1$σ$ spread, and $\lesssim$10% distortion within the 2$σ$ spread). This shows great prospects of using lensed SNe Ia to obtain intrinsic early-phase SN spectra for deciphering SN Ia progenitors. As a demonstration of the usefulness of lensed SNe Ia for cosmology, we simulate a sample of mock lensed SN Ia systems that are expected to have accurate and precise time-delay measurements in the era of the Rubin Observatory Legacy Survey of Space and Time (LSST). Adopting realistic yet conservative uncertainties on their time-delay distances and lens angular diameter distances (of 6.6% and 5%, respectively), we find that a sample of 20 lensed SNe Ia would allow a constraint on the Hubble constant ($H_0$) with 1.3% uncertainty in the flat $Λ$CDM cosmology. We find a similar constraint on $H_0$ in an open $Λ$CDM cosmology, while the constraint degrades to $3\%$ in a flat $w$CDM cosmology. We anticipate lensed SNe to be an independent and powerful probe of SN physics and cosmology in the upcoming LSST era. △ Less

Submitted 11 November, 2020; v1 submitted 19 February, 2020; originally announced February 2020.

Comments: 15 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 644, A162 (2020)

The inner dark matter distribution of the Cosmic Horseshoe (J1148+1930) with gravitational lensing and dynamics

Authors: S. Schuldt, G. Chirivì, S. H. Suyu, A. Yıldırım, A. Sonnenfeld, A. Halkola, G. F. Lewis

Abstract: We present a detailed analysis of the inner mass structure of the Cosmic Horseshoe (J1148+1930) strong gravitational lens system observed with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3). In addition to the spectacular Einstein ring, this systems shows a radial arc. We obtained the redshift of the radial arc counter image $z_\text{s,r} = 1.961 \pm 0.001$ from Gemini observations. T… ▽ More We present a detailed analysis of the inner mass structure of the Cosmic Horseshoe (J1148+1930) strong gravitational lens system observed with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3). In addition to the spectacular Einstein ring, this systems shows a radial arc. We obtained the redshift of the radial arc counter image $z_\text{s,r} = 1.961 \pm 0.001$ from Gemini observations. To disentangle the dark and luminous matter, we consider three different profiles for the dark matter distribution: a power-law profile, the NFW, and a generalized version of the NFW profile. For the luminous matter distribution, we base it on the observed light distribution that is fitted with three components: a point mass for the central light component resembling an active galactic nucleus, and the remaining two extended light components scaled by a constant M/L. To constrain the model further, we include published velocity dispersion measurements of the lens galaxy and perform a self-consistent lensing and axisymmetric Jeans dynamical modeling. Our model fits well to the observations including the radial arc, independent of the dark matter profile. Depending on the dark matter profile, we get a dark matter fraction between 60 % and 70 %. With our composite mass model we find that the radial arc helps to constrain the inner dark matter distribution of the Cosmic Hoseshoe independently of the dark matter profile. △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: 19 pages, 14 figures, 8 tables, submitted to A&A

Journal ref: A&A 631, A40 (2019)