Context. Gravitationally lensed type Ia supernovae (glSNe Ia) are unique astronomical tools that can be used to study cosmological parameters, distributions of dark matter, the astrophysics of the supernovae, and the intervening lensing galaxies themselves. A small number of highly magnified glSNe Ia have been discovered by ground-based telescopes such as the Zwicky Transient Facility (ZTF), but simulations predict that a fainter, undetected population may also exist. Aims. We present a systematic search for glSNe Ia in the ZTF archive of alerts distributed from June 1 2019 to September 1 2022. Methods. Using the AMPEL platform, we developed a pipeline that distinguishes candidate glSNe Ia from other variable sources. Initial cuts were applied to the ZTF alert photometry (with constraints on the peak absolute magnitude and the distance to a catalogue-matched galaxy, as examples) before forced photometry was obtained for the remaining candidates. Additional cuts were applied to refine the candidates based on their light curve colours, lens galaxy colours, and the resulting parameters from fits to the SALT2 SN Ia template. The candidates were also cross-matched with the DESI spectroscopic catalogue. Results. Seven transients were identified that passed all the cuts and had an associated galaxy DESI redshift, which we present as glSN Ia candidates. Although superluminous supernovae (SLSNe) cannot be fully rejected as contaminants, two events, ZTF19abpjicm and ZTF22aahmovu, are significantly different from typical SLSNe and their light curves can be modelled as two-image glSN Ia systems. From this two-image modelling, we estimate time delays of 22 ± 3 and 34 ± 1 days for the two events, respectively, which suggests that we have uncovered a population of glSNe Ia with longer time delays. Conclusions. The pipeline is efficient and sensitive enough to parse full alert streams. It is currently being applied to the live ZTF alert stream to identify and follow-up future candidates while active. This pipeline could be the foundation for glSNe Ia searches in future surveys, such as the Rubin Observatory Legacy Survey of Space and Time.

For approximately a decade, we have finally entered the era of discoveries of multiply imaged gravitationally lensed supernovae. To date, all cluster-lensed supernovae, very distant, faint and spatially resolved, have been found from space. In contrast, those deflected by individual galaxies have been very compact and bright enough to be identified with wide-field ground-based surveys through the magnification of the ‘standard candles’ method, i.e. without the need to spatially resolve the individual images. We review the challenges in identifying these extremely rare events, as well as the unique opportunities they offer for two major applications: time-delay cosmography and the study of the properties of deflecting bodies acting as lenses.

This article is part of the Theo Murphy meeting issue ‘Multi-messenger gravitational lensing (Part 1)’.

This thesis focuses on two important topics in astrophysics: the detection of kilonovae (KNe) and gravitationally lensed supernovae (glSNe) in optical surveys. In the first part, the study quantifies the impact of survey depth and choice of filters on the detection probability of KNe. The results highlight the importance of accounting for asymmetries expected for KNe, and despite several search campaigns, no KNe were detected by the Zwicky Transient Facility (ZTF). Nonetheless, non-detection studies provided meaningful constraints on the luminosity function and on the rates of KNe. The findings contribute to advancing our understanding of these rare, fast, and faint transients. I also discuss the improvements in measuring the Hubble constant with follow-up data of KNe, including broadband photometry and spectrophotometric data from the upcoming IFU instrument MAAT. 

The second part of the thesis focuses on gravitationally lensed supernovae. The ZTF survey was expected to detect more than one strongly lensed supernova per year, but only one was identified in the first five years. The study presents simulations of lightcurves for lensed supernovae and new rates based on realistic survey simulations for ZTF. Optimal cuts to distinguish lensed supernovae from normal unlensed supernovae are also provided. The thesis discusses time delay and lightcurve modeling for the one event found during ZTF, SN Zwicky, and the lessons learned from it.

The techniques developed in this thesis can be applied to future surveys to increase the detection rate of KNe and glSNe. These events and their underlying physics provide valuable insights in cosmology.

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

Fulfilling the rich promise of rapid advances in time-domain astronomy is only possible through confronting our observations with physical models and extracting the parameters that best describe what we see. Here, we introduce REDBACK; a Bayesian inference software package for electromagnetic transients. REDBACK provides an object-orientated PYTHON interface to over 12 different samplers and over 100 different models for kilonovae, supernovae, gamma-ray burst afterglows, tidal disruption events, engine-driven transients among other explosive transients. The models range in complexity from simple analytical and semi-analytical models to surrogates built upon numerical simulations accelerated via machine learning. REDBACK also provides a simple interface for downloading and processing data from various catalogues such as Swift and FINK. The software can also serve as an engine to simulate transients for telescopes such as the Zwicky Transient Facility and Vera Rubin with realistic cadences, limiting magnitudes, and sky coverage or a hypothetical user-constructed survey or a generic transient for target-of-opportunity observations with different telescopes. As a demonstration of its capabilities, we show how REDBACK can be used to jointly fit the spectrum and photometry of a kilonova, enabling a more powerful, holistic probe into the properties of a transient. We also showcase general examples of how REDBACK can be used as a tool to simulate transients for realistic surveys, fit models to real, simulated, or private data, multimessenger inference with gravitational waves, and serve as an end-to-end software toolkit for parameter estimation and interpreting the nature of electromagnetic transients.

During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network, the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star–black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (False Alarm Rate less than 1 yr⁻¹) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg² within the skymaps of the high-significance GW events, reaching a median depth of g ≈ 20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than −16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH events), along with our O4a follow-ups. Following Kasliwal et al., we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag day⁻¹ can peak at a magnitude brighter than −17.5 mag.

In the version of this article initially published, there was in an error in the third-to-last sentence of the abstract, now reading, in part, “we calculate a rate of 0.02–0.01 +0.04 Gpc–3 yr–1”, where Gpc was spelled out as gigapascals, not gigaparsecs. Also, the scale label (2″) was missing in the lower-left corner of Fig. 1b. The errors have been corrected in the HTML and PDF versions of the article.

Detecting gravitationally lensed supernovae is among the biggest challenges in astronomy. It involves a combination of two very rare phenomena: catching the transient signal of a stellar explosion in a distant galaxy and observing it through a nearly perfectly aligned foreground galaxy that deflects light towards the observer. Here we describe how high-cadence optical observations with the Zwicky Transient Facility, with its unparalleled large field of view, led to the detection of a multiply imaged type Ia supernova, SN Zwicky, also known as SN 2022qmx. Magnified nearly 25-fold, the system was found thanks to the standard candle nature of type Ia supernovae. High-spatial-resolution imaging with the Keck telescope resolved four images of the supernova with very small angular separation, corresponding to an Einstein radius of only θ_E = 0.167″ and almost identical arrival times. The small θ_E and faintness of the lensing galaxy are very unusual, highlighting the importance of supernovae to fully characterize the properties of galaxy-scale gravitational lenses, including the impact of galaxy substructures.

Tidal disruption events (TDEs) are bursts of electromagnetic energy that are released when supermassive black holes at the centres of galaxies violently disrupt a star that passes too close¹. TDEs provide a window through which to study accretion onto supermassive black holes; in some rare cases, this accretion leads to launching of a relativistic jet^{2,3,4,5,6,7,8,9}, but the necessary conditions are not fully understood. The best-studied jetted TDE so far is Swift J1644+57, which was discovered in γ-rays, but was too obscured by dust to be seen at optical wavelengths. Here we report the optical detection of AT2022cmc, a rapidly fading source at cosmological distance (redshift z = 1.19325) the unique light curve of which transitioned into a luminous plateau within days. Observations of a bright counterpart at other wavelengths, including X-ray, submillimetre and radio, supports the interpretation of AT2022cmc as a jetted TDE containing a synchrotron ‘afterglow’, probably launched by a supermassive black hole with spin greater than approximately 0.3. Using four years of Zwicky Transient Facility¹⁰ survey data, we calculate a rate of 0.02^+0.04_−0.01 Gpc⁻³ yr⁻¹ for on-axis jetted TDEs on the basis of the luminous, fast-fading red component, thus providing a measurement complementary to the rates derived from X-ray and radio observations¹¹. Correcting for the beaming angle effects, this rate confirms that approximately 1 per cent of TDEs have relativistic jets. Optical surveys can use AT2022cmc as a prototype to unveil a population of jetted TDEs.

The merger of two compact objects of which at least one is a neutron star is signalled by transient electromagnetic emission in a kilonova (KN). This event is accompanied by gravitational waves and possibly other radiation messengers such as neutrinos or cosmic rays. The electromagnetic emission arises from the radioactive decay of heavy r-process elements synthesized in the material ejected during and after the merger. In this paper we show that the analysis of KNe light curves can provide cosmological distance measurements and constrain the properties of the ejecta. In this respect, MAAT, the new Integral Field Unit in the OSIRIS spectrograph on the 10.4 m Gran Telescopio CANARIAS (GTC), is well suited for the study of KNe by performing absolute spectro-photometry over the entire 3600 − 10 000 Å spectral range. Here, we study the most representative cases regarding the scientific interest of KNe from binary neutron stars, and we evaluate the observational prospects and performance of MAAT on the GTC to do the following: (a) study the impact of the equation of state on the KN light curve, and determine to what extent bounds on neutron star (NS) radii or compactness deriving from KN peak magnitudes can be identified and (b) measure the Hubble constant, H₀, with precision improved by up to 40%, when both gravitational wave data and photometric-light curves are used. In this context we discuss how the equation of state, the viewing angle, and the distance affect the precision and estimated value of H₀.