Strong Gravitational Lensing and Microlensing of Supernovae

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


Brief history
In an insightful and pioneering publication, Refsdal (1964) pointed out that supernovae (SNe) would be particularly interesting sources for studies involving strong gravitational lensing.This was arguably among the biggest early milestones in the field of gravitational lensing, following the realisation by Zwicky (1937) that the scenario, first proposed by Einstein in 1936, had realistic applications in extragalactic astronomy.Refsdal pointed out that the measurement of the time delay between the arrival of SN images could be used to infer the Hubble constant (H 0 ).

History of searches for lensed SNe behind clusters
The idea to use "gravitational telescopes", i.e., known lensing galaxy clusters, to boost the faint signals from distant supernovae started to get traction about thirty years ago (Kovner and Paczynski, 1988;Sullivan et al., 2000;Gunnarsson and Goobar, 2003).Since the lensing magnification boosts the signal from the faint distant source behind the lens by a factor of µ, but leaves the dominant foreground sky noise unaffected, the signal-to-noise ratio scales as SN R ∝ µ √ t, hence leading to a gain factor µ 2 in exposure length.However, the solid angle at the source planes shrinks by a factor µ behind the lens.Thus, the net gain/loss of searching for supernovae behind massive clusters is therefore a non-trivial combination of field of view, limiting depth, and supernova luminosity functions.Although not spectroscopically confirmed, several SN candidates were eventually found in cadenced observations in the Near-IR at the Very Large Telescopes (VLT) of the European Southern Observatory (ESO) (Stanishev et al., 2009;Goobar et al., 2009), including a z = 1.7 core-collapse SN behind Abell 1689, with a lens model inferred magnification µ = 4.3 ± 0.3 (Amanullah et al., 2011).In spite of the limitations of the survey, Petrushevska et al. (2016) demonstrated the power of gravitational lensing to place meaningful limits on the rate of core-collapse supernovae at very high redshifts.
1.2 Discoveries!PS1-10afx was first reported by Chornock et al. (2013) as an unusual superluminous SN.However, shortly thereafter, Quimby et al. (2013) showed that it was a perfect match to a lensed SN Ia at redshift z = 1.388, with a large amplification, µ ∼ 30.Quimby et al. (2014) eventually also identified a foreground lens, at z = 1.117.Even if the putative lensing object was too close to the SN, it was consistent with lensing models and existing data.Since the lensed SN Ia classification only gained acceptance four years after the explosion, further investigation of the lensing nature of the SN or its type was not possible.Moreover, a subsequent single-band HST imaging (PI: Chornock) proved insufficient to verify the presence of any lensed images of the SN host galaxy.
The first detection of multiple images from a supernova came with SN Refsdal at redshift z = 1.49(Kelly et al., 2015), found in Hubble Space Telescope (HST) surveys of the massive lensing galaxy cluster MACS J1149.6+2223 at z = 0.54.Based on both light curve shapes and spectroscopy, SN Refsdal was classified as a Type II SN, resembling the iconic SN 1987A in the LMC (Kelly et al., 2016b).This is a peculiar class of faint SNe, quite rare in the local universe.At the time of discovery, 4 multiple images of the SN were visible with HST (images S1-S4, see Fig. 1 taken from Grillo et al. 2018), in a cross configuration around one of the foreground cluster galaxies.There are two other images of the background spiral galaxy hosting the SN, labelled as SX (inset (a) in Fig. 1) and SY (the northern-most image, just outside the top-right corner of inset (a)).
Using mass models of the cluster, the discovery team (Kelly et al., 2015) predicted that SY happened before S1-S4, whereas SX would appear after S1-S4.However, the time of SX reappearance was uncertain from their model, ranging from 1 to 10 years.To refine the delay prediction, Grillo et al. (2016) used the Multi Unit Spectroscopic Explorer (MUSE; Bacon et al., 2010) on ESO's VLT to obtain spectroscopic observations of the field.The spectroscopic redshift measurements helped to separate foreground lenses from background sources.This spectroscopic data set was shared with multiple modelling teams, who attempted to predict the reappearance of image SX using the new observational data.Teams predicted their delays before the reappearance (Treu et al., 2016;Jauzac et al., 2016;Grillo et al., 2016;Kawamata et al., 2016), providing a true blind test of the predictions and the modelling capabilities.Most teams predicted a short time delay of SX with respect to S1, within a year's time.Kelly et al. (2016a) detected the reappearance of SX in December 2015, and the resulting constraint on the time delay and magnification of SX agreed well with the predictions from two teams (Grillo et al., 2016;Kawamata et al., 2016).
Through the monitoring of the multiple images of SN Refsdal, Rodney et al. (2016) have measured the time delays and magnification ratios among images S1, Fig. 1 Hubble Space Telescope image of SN Refsdal, the first strongly lensed SN system with spatially resolved images.Inset (a) shows the image SX that was detected in December 2015 (Kelly et al., 2016a), and inset (b) shows the multiple images S1, S2, S3 and S4 that were first discovered in November 2014 (Kelly et al., 2015).There is another image SY located in the northern most image of the spiral host galaxy, next to the top-right corner of the inset (a).Image taken from Grillo et al. (2018).Original image credit: NASA, ESA/Hubble.S2, S3 and S4.The detection of the reappearance of image SX provided an estimate of its time delay with respect to the other images (Kelly et al., 2016a), and the timedelay measurement from subsequent monitoring is forthcoming (P.Kelly, private communications).The SX time delay is expected to be precise with uncertainties of at most a few percent given its long delay of ∼1 year, providing a great opportunity of measuring the Hubble constant with gravitational lens time delays (e.g., Grillo et al., 2018, 2020, Kelly et al. in prep.).
The first spatially resolved multiply-imaged Type Ia supernova, iPTF16geu, was detected by the intermediate Palomar Transient Factory (Goobar et al., 2017).As shown in Fig. 2, the ground-based imager at iPTF could not spatially resolve the very compact system with Einstein radius θ E = 0.3 .However, because of the "standard candle" nature of SNe Ia, it became clear that strong lensing was the most likely explanation, since the SN was more than 30 standard deviations brighter than the expectations for a SN Ia at its measured redshift, z = 0.409.The spectra used to classify the supernova showed spectral lines from both the host galaxy at the same redshift, but also the deflecting galaxy at z = 0.216.Thanks to multi-band follow-up with HST, an accurate (model independent) measurement of the magnification was made, µ = 67.8+2.6  −2.9 (Dhawan et al., 2020), after correction for non-negligible extinction by dust in both the host and lens galaxies.The time delays between the SN images for this system were very small, about a day or less (More et al., 2017;Dhawan et al., 2020).The flux ratios between the supernova images (see Fig. 2) were not consistent with expectations from lensing, hinting contributions from stellar microlensing (More et al., 2017).This was further confirmed after differential dust extinction from within the lensing galaxy was also accounted for (Mörtsell et al., 2020).
In 2021, another strongly lensed SN, SN Requiem, was discovered in archival HST imaging of the galaxy cluster MACSJ0138.0−2155by Rodney et al. (2021).A bright quiescent and evolved elliptical galaxy at redshift z = 1.95 is lensed into giant arcs by the galaxy cluster, making this system one of the brightest NIR objects (e.g., Newman et al., 2018).Near the giant arcs are three point sources that were present in the single-epoch image in 2016 (Newman et al., 2018), but absent in the 2019 HST image (HST Program -REQUIEM; HST-GO-15663; PI:Akhshik).These 3 point sources are identified as multiple images of a SN hosted in the z = 1.95 galaxy, and the SN is likely to be of Type Ia.Based on the differences in the colour of the 3 SN images and SN Ia templates, the time delays could be estimated from the single-epoch observation in 2016, although with large uncertainties.The predicted time delays suggest that the fourth image of SN Requiem will appear in the future (circa 2037).
Recently, Goobar et al. (2022) discovered another strongly lensed SNe Ia, named Supernova Zwicky (a.k.a.SN 2022qmx), in the Zwicky Transient Facility (ZTF; Bellm et al., 2019).The discovery of this system is similar to the case of iPTF16geu -the multiple lensed SN Ia images are closely separated (with an Einstein radius of θ E ∼ 0.17 ) and not spatially resolved by ZTF, thus resulting in a brightness through lensing magnification that is substantially higher than expected for a SN Ia.Pierel et al. (2022)   These first few strongly lensed SNe have opened a new window to probe cosmology and supernovae through these spectacular phenomena.Two reviews (Oguri, 2019;Liao et al., 2022) on gravitationally lensed transients have illustrated the exciting utilities of lensed transients.Compared to these two reviews, we focus and dive deeper into the specific case of strongly lensed supernovae.In Section 2, we review the potential of lensed SNe as a cosmological probe.In Section 3, we describe the use of lensed SNe as an astrophysical probe, such as constraining the SN progenitors and dust in galaxies.In Section 4, we present various searches of lensed SNe and the expected rates in current and upcoming surveys.We summarise in Section 5.

Cosmological probe
While the idea of using lensed SNe for measuring the Hubble constant dates back to Refsdal (1964), this time-delay method was first realised with lensed quasars that were discovered decades before lensed SNe, starting with the first lensed quasar system found by Walsh et al. (1979).The use of lensed quasars as a cosmological probe is reviewed by Birrer et al. (2022) in the same series (see also, e.g., Treu and Marshall, 2016;Suyu et al., 2018;Treu et al., 2022).We briefly summarise the approach of time-delay cosmography in Section 2.1, and focus on SNe as lensed background sources, which are "standardisable candles" and a relative distance indicator.In Section 2.2, we describe approaches to measure time delays of lensed SNe, particularly novel ones through spectroscopic observations of SNe, which complement conventional techniques using photometric light curves.In Section 2.3, we show the present cosmological constraints and future forecasts from lensed SNe.

Jackpot: two cosmological probes in one
We describe the various distance measurements that we can obtain from strongly lensed SNe, through the lensing time delays and the SN light curves.We briefly summarise the determination of the lensing distances that are covered in detail in the parallel reviews by, e.g., Birrer et al. (2022), and focus particularly on the new aspects that SNe bring.
The expression for the time delay between images i and j of a lensed SN is where D ∆t is the time-delay distance, c is the speed of light, and τ is the Fermat potential.The time-delay distance (e.g., Refsdal, 1964;Suyu et al., 2010) is defined by where z d is the deflector (lens) redshift, and D d , D s , and D ds are the angular diameter distance to the deflector, to the source, and between the deflector and the source, respectively.The Fermat potential difference between two SN image positions θ i and θ j (with corresponding source position β) is where the Fermat potential is and ψ(θ) is the lens potential.The dimensionless surface mass density κ(θ), a.k.a.lensing convergence, is related to the lens potential via the Poisson equation, From Equation (1), we see that by measuring the time delays ∆t ij and modelling the deflection and line-of-sight mass distributions to get ∆τ ij , we can infer the time-delay distance D ∆t .Since D ∆t is a combination of three angular diameter distances (Equation 2) and is thus proportional to H −1 0 , time-delay lenses allow us to measure directly H 0 with weak dependence on other cosmological parameters.
When the lens galaxy stellar velocity dispersion is measured, then the combination of time delays, velocity dispersion and high-resolution imaging of the lens system allows D d to be measured (e.g., Paraficz and Hjorth, 2009;Jee et al., 2015Jee et al., , 2019)), in addition to D ∆t .The joint constraint of D ∆t and D d provides even more leverage on determining cosmological parameters (e.g., Jee et al., 2016;Yıldırım et al., 2020).
It is clear from Equation (4) that in order to constrain H 0 with strong lens time delays, one needs to know both the 2D lens potential and the unlensed source position, neither of which are directly observable.The use of lens modelling is required to infer these quantities.However, strong lens models are subject to degeneracies, which are a major source of uncertainty for time-delay cosmography (Schneider and Sluse, 2014).The largest problem is the mass-sheet degeneracy: re-scaling any model of the lensing convergence κ and adding a constant sheet of surface mass density leaves the predicted images unchanged, but alters the predicted time delay between the images (Falco et al., 1985).Breaking the mass-sheet degeneracy is therefore necessary to constrain H 0 from any lens with observations of image positions and time delays.In order to break the mass-sheet degeneracy, additional information is required, either non-lensing information (e.g., the dynamical mass of the system), the presence of sources at multiple redshifts (though see Schneider, 2014) or knowledge of the intrinsic luminosity or size of the unlensed source (Kolatt and Bartelmann, 1998).It is on this final point -a known luminosity of the source -that strongly lensed SNe shine as a powerful potential cosmological probe.
In this context, Type Ia SNe play a special role.Thanks to the homogeneity of their luminosity, these thermonuclear explosions could be used as sharp distance estimators to measure the accelerated expansion of the universe (Riess et al., 1998;Perlmutter et al., 1999), leading to the discovery of dark energy.State-of-the-art SN Ia cosmological surveys have demonstrated that the intrinsic luminosity of SNe Ia only varies by about 0.1 mag (Betoule et al., 2014;Scolnic et al., 2018).Birrer et al. (2021) have shown the great benefit of these sharp distance indicators to break the mass-sheet degeneracy to improve the accuracy of the measurements of the Hubble constant from time delays.
Despite the great opportunity presented by having a lensed standard candle, the exploitation of lensed SNe poses an additional challenge: microlensing by stars (Dobler and Keeton, 2006).The intrinsic size of a SN is comparable to the Einstein radius of an individual star in the lens galaxy, a case analogous to microlensing of lensed quasars (Vernardos et al. 2023, in prep.).This means that observed magnification is sensitive to the -essentially unobservable -positions of stars in the lensing galaxy.Worse, as the supernova expands, the SN atmosphere crosses microlensing caustics.This causes a change in total magnification of the supernova and differential magnification across the atmosphere (Bagherpour et al., 2006;Goldstein et al., 2018;Huber et al., 2019;Mörtsell et al., 2020); this means that the light curves of each image can look quite different even though they are formed from the same background source, as illustrated in Fig. 4.This makes it a challenge for measuring accurate time delays between images, though early colour curves can overcome this due to an achromatic phase in the microlensing at early times (Goldstein et al., 2018;Huber et al., 2021b) as can spectroscopic observations of absorption lines in the SN atmosphere (Johansson et al., 2021;Bayer et al., 2021).The SNe never reach a large enough size for the microlensing effect of multiple stars to average out, so microlensing makes it hard to standardise lensed SNe.The scale of this problem depends in detail on the magnification of the image, the type of the image (time-delay minimum, maximum or saddle) and fraction of density in stars at the location of the image.For high macro-magnification images, the scatter can be as large as 1 magnitude in flux (Schechter and Wambsganss, 2002;Yahalomi et al., 2017;Diego et al., 2022), entirely washing out the possibility of using lensed SNe as a standard candle.For images forming further from the Einstein Radius, and particularly for large mass systems where the dark matter fraction is higher at the Einstein radius, the scatter is much smaller and can be less than the intrinsic scatter of a Type Ia SN.Foxley-Marrable et al. (2018) found that 20 percent of lensed SNe will contain an image with scatter comparable to the intrinsic scatter of a type Ia supernova luminosity, assuming the dark matter fractions of elliptical galaxies derived in Auger et al. (2010) and a Saltpeter IMF.However, these standardisable images only form far beyond the Einstein radius of the lens.This means that the counter images are very close to the centre of the lens and demagnified.Measuring time delays will therefore be challenging, although the faint images arrive second so follow-up with a larger telescope may ameliorate this problem.
The 20 percent standardisable fraction is subject to two key caveats: the stellar initial mass function is assumed to be Salpeter-like and the ellipticity of the lensing mass follows that of the light.The assumption of a particular IMF sets the stellar-to-dark matter fraction.Fig. 5 shows the importance of this assumption: at very high stellar fraction there are no strongly lensed images with minimal microlensing scatter (Weisenbach et al., 2021).With a population of lensed SNe, Foxley-Marrable et al. ( 2018) find that observing the amount of microlensing induced scatter can be used to constrain the IMF of the lens galaxy, with strong sensitivity for lenses with Einstein radius between 0.2 and 0.5 .The assumption that mass follows light produces more elliptical lens mass distributions than is seen in cold-dark-matter-only simulations.Increasing the ellipticity of the lens moves the dashed lines in Fig. 5 to the left.Strongly lensed images can form at lower magnification, where microlensing is less significant.If real lenses are more spherical, then the standardisable fraction will fall.Therefore, the standardisation of lensed supernovae is subject to somewhat uncertain astrophysics.Exactly what we observe from a large sample of lensed supernovae will inform not just cosmology but also the astrophysics of matter in the strong lensing galaxy.
Despite microlensing distortions on the light curves of SNe and reductions in the number of standardisable SN, lensed SNe still have various advantages over the more conventional lensed quasars for time-delay cosmography.The drastically varying brightness of lensed SNe enable shorter monitoring campaigns (months) to obtain the light curves and thus more efficient measurements of time delays.Simulations of microlensed SN light curves with realistic photometric uncertainties showed that time delays can still be recovered accurately and precisely (Huber et al., 2019;Pierel and Rodney, 2019), and the dominant source of uncertainty in the time delays is typically photometric uncertainties rather than microlensing distortions (Huber et al., 2021a).Another advantage of lensed SNe is that SNe fade after several months, revealing both the lens galaxy and the lensed SN host galaxy more clearly that enable more accurate lens mass modelling (e.g., Ding et al., 2021).In particular, spatially-resolved stellar kinematics of the foreground lens galaxy can be more readily acquired after the SN images fade, and the combination of lensing and kinematic data allows one to break the mass-sheet degeneracy (e.g., Barnabè et al., 2012;Yıldırım et al., 2020Yıldırım et al., , 2021;;Birrer et al., 2020;Shajib et al., 2023).Stellar kinematic maps of the lensed SN host galaxy would further constrain the lens mass distribution (Chirivì et al., 2020).We therefore anticipate lensed SNe to be an efficient cosmological probe in the upcoming era of time-domain astronomy when hundreds of such lensed SN events are expected to happen (e.g., Oguri and Marshall, 2010;Quimby et al., 2014;Goldstein et al., 2019;Wojtak et al., 2019).

Time-delay measurements
The time delays between multiple lensed images are a key ingredient for time-delay cosmography.The fractional uncertainties in the time delays translate directly to the fractional uncertainties in the D ∆t and D d measurements (Equation 1).
Traditionally, the measurements of time delays involve monitoring the lens system to acquire the light curves of either lensed quasars (e.g., Fassnacht et al., 2002;Courbin et al., 2018;Millon et al., 2020a,b) or lensed supernovae (e.g., Rodney et al., 2016;Dhawan et al., 2020).Fig. 6 shows an example of the light curves of the four images (S1, S2, S3 and S4) of SN Refsdal obtained by Rodney et al. (2016).Curve-shifting techniques that account for microlensing distortions are then applied to these light curves to extract the time delays (e.g., Tewes et al., 2013;Pierel and Rodney, 2019;Millon et al., 2020c).By making use of the characteristic SN Ia SEDs, Dhawan et al. (2020) used SN Ia template light curves to fit to the HST images of iPTF16geu and measured the 3 independent time delays between the four SN images.New machine learning approaches such as random forests are also being developed to infer the delays of microlensed SN Ia systems from the light curves (e.g., Huber et al., 2021a).
SNe are not only drastically changing their brightness, but their spectra also evolve substantially.The spectroscopic evolution of supernovae offers a new avenue to measure the time delays that are complementary and competitive to the light-curve techniques, as demonstrated by Johansson et al. (2021) for the case of iPTF16geu.The photospheric velocity in SNe Ia evolves with time as v(t) ∝ t −0.22 (Piro and Nakar, 2014), which allows the use of spectral dating to constrain the time delay between images, especially at early times when the velocities of the spectral features is high and the changes from day-to-day are significant.For the most common type of core-collapse supernovae, SNe IIP, where the light curves rise fast to peak and then reach a plateau phase that could last for hundreds of days, spectral time delays may be the best way forward.Fig. 7 shows an example of the spectra of Type IIP SN1999em from the TARDIS simulation (Kerzendorf and Sim, 2014;Vogl et al., 2019Vogl et al., , 2020) ) at multiple epochs after explosion (Bayer et al., 2021).These epochs (in rest-frame days) are in the plateau phase of this Type IIP SN.The spectra show prominent absorption lines, especially Hβ, FeII and Hα that are indicated on the figure, and the absorption wavelength increases as a function of time.It is precisely such evolutions in the spectral features that allow us to infer the time delays between two SN images.
In more detail, a sequence of spectra of the first appearing SN image, say SN image A (such as those shown in Fig. 7) yields the absorption wavelength of each spectral feature as a function of the SN phase (days after explosion).When we obtain a single spectra of a trailing SN image, say SN image B, and measure the absorption wavelength of the same spectral feature, we can determine the phase of SN B by using the relation between the absorption wavelength and phase from SN A. The determined phase of SN B relative to SN A, together with the known observational times of the epochs, then allows us to compute the time delays between SN images A and B. Microlensing of SN introduces scatter in the relation between the absorption wavelength and phase, and this effect can be quantified using microlensing maps such as those from the GERLUMPH software (Vernardos and Fluke, 2014;Vernardos et al., 2014Vernardos et al., , 2015)).Accounting for the effect of microlensing that distorts spectra, Bayer et al. (2021) have shown that this technique can yield time delays with uncertainties of ∼2 days per spectral feature with signal-to-noise in the spectra of ∼20 per wavelength bin of 3Å width.With multiple spectral features and epochs, the inferred time delays from such spectral approach will be even more precise and accurate.

Cosmography with lensed SNe: present cosmological constraints and future forecasts
Of the five known lensed SN systems with spatially resolved SN images, SN Refsdal is the most promising system in delivering a precise and accurate H 0 measurement (with <10% uncertainty) in the near future.As mentioned in Section 1, the measurement of the time delay of image SX relative to the first image S1 is forthcoming, and is expected to have uncertainties of only a few percent (P.Kelly, private communications).For a preview of the constraining power on H 0 from SN Refsdal, Fig. 8 shows a forecast of the inferred H 0 for a range of hypothetical time-delay measurements of image SX relative to image S1, assuming the flat ΛCDM cosmological model (Grillo et al., 2020).These are based on cluster mass models of Grillo et al. (2020) that incorporate all known sources of uncertainties.We anticipate that a 3% uncertainty on SX time delay would yield approximately 6% and 40% uncertainties (1σ, including both statistical and systematic) for H 0 and Ω m , respectively.
Various studies indicate that current and future surveys would discover dozens to hundreds of lensed SNe (Goobar et al., 2002;Oguri and Marshall, 2010;Goldstein et al., 2019;Wojtak et al., 2019;Pierel et al., 2021).In particular, Oguri and Marshall (2010) anticipated thousands of lensed quasars and ∼100 lensed SNe to be detected in the Rubin Observatory Legacy Survey of Space and Time (LSST; LSST Science Collaboration et al., 2009).Using a mock sample of ∼1500 well-observed lenses (consisting of 1476 lensed quasars and 66 lensed SNe) and assuming Planck priors in a flat Universe, Oguri and Marshall (2010) expected the following 1σ uncertainties on cosmological parameters: σ(w 0 ) = 0.15, σ(w a ) = 0.41 and σ(h) = 0.017 where (w 0 , w a ) are the time-independent and time-dependent components of the dark energy equation-of-state parameter, and h is H 0 in units of 100 km s −1 Mpc −1 .Huber et al. (2019) analyzed in detail the number of lensed SNe Ia that could yield time-delay measurements with precisions better than 5% and accuracies better than 1%, with realistic microlensed SN Ia light curves and LSST observing strategies.While LSST is efficient at detecting lensed SNe given its wide survey area and depth, the observing cadence per filter is not rapid enough to map out light curves for precise delay measurements.Follow-up observations at the cadence  2020) simulated a sample of 20 lensed SNe Ia from LSST that are expected to have good time delays.Assuming that these systems would have high-resolution imaging and spatially resolved kinematic measurements to break the mass-sheet degeneracy and yield 6.6% uncertainty on D ∆t and 5% uncertainty on D d , this modest sample of 20 lensed SNe Ia could yield constraints on H 0 and Ω m of 1.3% and 19%, respectively, in flat ΛCDM, as illustrated in Fig. 9.In an open ΛCDM cosmology (i.e., allowing for a spatially curved Universe), a similar constraint on H 0 compared to flat ΛCDM is achievable, while in the flat wCDM cosmology (where the dark energy equation-of-state parameter w is allowed to vary), the H 0 constraint degrades to 3% (Suyu et al., 2020).
The forecast made by Suyu et al. (2020) did not make use of the standardisable nature of SNe Ia, due to the effects of microlensing and millilensing as mentioned in Section 2.1.Nonetheless, Foxley-Marrable et al. (2018) showed that lensed SNe Ia with asymmetric configurations, i.e., with a SN image located far outside the Einstein radius of the foreground lens, have low microlensing scatter of 0.15 mag for the outer SN image, which is comparable to the intrinsic dispersion of a typical SN Ia.These systems are thus standardisable in terms of overcoming microlensing perturbations, and Foxley-Marrable et al. ( 2018) estimated that about ∼22% of the ∼ 930 LSST systems predicted by Goldstein et al. (2018) would be standardisable.However, Goldstein et al. (2018)  ple of ∼140 lensed SNe Ia that are standardisable and with potential time-delay measurement.
Following the estimates of Foxley-Marrable et al. ( 2018) and by making use of 144 standardisable lensed SNe Ia, Birrer et al. (2021) combined non-lensed SNe Ia with lensed SNe Ia via a Bayesian hierarchical framework to infer the constraints on H 0 .The aim is to use the standardisable nature of SNe Ia to break the mass-sheet degeneracy, without using spatially resolved kinematics of the lens galaxy.Assuming optimistically that the time delays between SN image pairs can be measured with a precision of 2 days for all these systems and that the typical Einstein radius is ∼1 , Birrer et al. (2021) found that H 0 can be measured with an uncertainty of 1.5%.
The various cosmological forecasts depend sensitively on assumptions of followup observations in order to classify the SN type, measure the time delays, acquire high-resolution imaging and spectroscopy.As described in Section 2.1, these are necessary for measuring absolute distances to lensed SNe and using them as a cosmological probe.In Section 4.4, we discuss in more detail the follow-up requirements for lensed SNe.

Astrophysical probe
Lensed SNe are powerful astrophysical tools to study supernovae and galaxies, in addition to probing cosmology through their time delays.
In Section 3.1, we explain how lensed SNe provide an excellent opportunity to constrain SN progenitors.In Section 3.2, we describe how lensing magnifications enable the acquisition of high-redshift SN spectra that are crucial for SN cosmology.In Section 3.3, we outline the use of microlensing to probe SN structure.In Section 3.4, we show lensed SNe as a probe of dust in the foreground lens galaxy.

SN physics and progenitors
Whilst it is accepted that a Type Ia SN is a thermonuclear detonation of a white dwarf (WD), the cause of the detonations remains uncertain even after decades of study (Maoz et al., 2014;Livio and Mazzali, 2018).
The classic progenitor model (often referred to as the single degenerate channel) for a Type Ia SN is a close binary of a white dwarf and a post main sequence star (Whelan and Iben, 1973).As the star expands, it overflows its Roche lobe and its outer layers are accreted onto the WD until the WD reaches the Chandrasekhar mass.At this point, electron degeneracy pressure can no longer support the WD, and gravitational collapse causes thermonuclear detonation of the WD.This model is attractive, since it naturally explains the standardisability of SNe Ia, with the progenitor always having the same mass and near identical composition.However the lack of evidence for surviving companions suggest that this channel cannot account for all SNe Ia (González Hernández et al., 2012;Lundqvist et al., 2015).
An alternative progenitor model is the double degenerate channel.Here two white dwarfs merge, exceeding the Chandrasekhar mass (Iben and Tutukov, 1984;Webbink, 1984).This mechanism leaves behind no companion star, but since each pair of white dwarfs will sum to a different mass, there is no obvious reason why this channel would produce a good standard candle.In addition to these two classic channels of single-and double-degenerate systems, there are other explosion mechanisms that have been explored, including sub-Chandrasekhar explosions (e.g., Sim et al., 2010), delayed detonations (e.g., Röpke et al., 2012), and double detonations of sub-Chandrasekhar WDs (e.g., Fink et al., 2007).
Early observations of SN light curves are critical in constraining the properties of SN progenitor systems (e.g.Kasen 2010;Piro et al. 2010;Rabinak and Waxman 2011;Nugent et al. 2011;Bloom et al. 2012;Goobar et al. 2015;Piro and Morozova 2016;Noebauer et al. 2017;Kochanek 2019;Fausnaugh et al. 2019;Yao et al. 2019;Miller et al. 2020;Bulla et al. 2020).If non-degenerate matter is close to the WD, then it should be shock heated by the explosion, producing excess high energy flux in the first few hours to days of the light curve.In addition to SNe Ia, early observations of core-collapse SNe are also helpful in constraining the properties, such as the sizes, of the progenitor stars (e.g., González-Gaitán et al., 2015;Chen et al., 2022).Even with the development of wide-field optical surveys, observing these earliest moments of SNe is heavily reliant on chance.Strong lensing gives us the opportunity to predict the precise reappearance time of a SN, such that very early data can be gathered (Suwa, 2018;Suyu et al., 2020).
There are multiple challenges to this approach: early detection of the first SN image, possible image demagnification, distortions of SN light curves and spectra by microlensing, and time delay predictions of insufficient precision.We discuss each of these in turn.
Detection of a lensed SN based on only the appearance of the first image is more difficult, due to possible confusions between a lensed and an unlensed SN without the multiple SN images present.This confusion can be overcome when the lensed SN host galaxy is visible as lensed arc features in addition to the foreground lens galaxy.
Whilst lensed SNe will be discovered from highly magnified images, images occurring with significant time delay after the bright "discovery image" are typically demagnified in the case of two-image systems (Foxley-Marrable et al., 2020).The situation is better for four-image systems where some of the trailing SN images can be brighter than the discovery SN image.We expect ∼ 2/3 of the lensed SNe to be two-image systems, and ∼ 1/3 to be four-image systems (Oguri and Marshall, 2010).
As seen in Fig. 4 and shown in Huber et al. (2019), microlensing distorts the SN light curves and spectra, potentially garbling the information on progenitors.Encouragingly, Suyu et al. (2020) have investigated the microlensing effects on four different explosion models of SNe Ia and demonstrated that spectral distortions due to microlensing are 1% at the 1σ level within 10 rest-frame days after explosions.Therefore, microlensing is expected to have negligible impact on earlyphase spectra (within ∼10 rest-frame days) for deciphering SN progenitors.
Finally, time delays cannot be precisely estimated from lens modelling alone.A 5-10 percent precision is typical from the best lens models (e.g., Shajib et al., 2019), and the uncertain value of the Hubble constant given the discrepant H 0 measurements (e.g., Verde et al., 2019;Di Valentino et al., 2021) adds several more percent of uncertainty in the predicted time delays.For systems with delays 20 days that are sufficiently long to catch the trailing SN images from their beginnings (given the time it takes to detect the first SN image and to arrange follow-up observations), a 10% uncertainty would translate to 2 days.Given this uncertainty, the predicted delay can easily miss the observations within ∼ 2 days after explosion, the crucially early moments for constraining SN progenitors.
To overcome these challenges, deep and high-cadence (ideally daily) monitorings with imagers after the detection of the first appearing SN image would greatly help.Such a monitoring would acquire the early-phase light curves of trailing SN images, and spectroscopic observations of the early phases can be triggered as soon as a trailing SN image appears in the monitoring.Early-phase spectroscopic observations especially in the rest-frame UV are unprecedented and important for constraining SN progenitors (Suyu et al., 2020).The time delays of lensed SNe provide a unique and exciting avenue to acquire such spectroscopic observations for studying SN progenitors.

High-z SN spectra through gravitational telescopes
The precision of SNe Ia as distance indicators, and thus their use to accurately constrain the nature of dark energy, is ultimately limited by our understanding of progenitor systems, e.g., the potential evolution of the SN properties, especially their absolute magnitude, over cosmic time.Thus, confirming the standard-candle nature through spectroscopic comparisons between SNe near and far is essential.With current instruments, detailed comparisons beyond z > 1 are not possible, due to the low SNR.This will remain challenging in the James Webb Space Telescope (JWST) era.However, multiple spectra of SNe at high redshift are required to check if the progenitor population evolves with redshift.Hence, the magnification of the signal from gravitational lensing could become essential to constrain one of the biggest systematic uncertainties in SN Ia cosmology (Petrushevska et al., 2017;Johansson et al., 2021).Fig. 10 shows the noisy HST spectrum, based on six hours of observations of SN "Primo" at z = 1.55 observed by Rodney et al. (2012) compared with the expected signal if observed for shorter time through gravitational telescopes of different strengths, ∆m = 1, 1.5, 2.0 and 2.5 mag, where ∆m is the difference in magnitude coming from lensing magnification.These magnifications are expected to be typical for lensed SNe (Goldstein et al., 2019).

SN structure with microlensing
When a background source crosses a caustic network produced by microlenses, this allows one to estimate the physical dimensions of the light emitting surface, most robustly its half-light radius (Mortonson et al., 2005;Vernardos and Tsagkatakis, 2019).For the case of SNe, their expanding ejecta lead to an increasing size that is directly linked to the expansion velocity.Because the (unobservable) distance from a caustic/high magnification region on the source plane is inversely proportional to changes in brightness due to microlensing, deformed SN light curves, like the ones shown in Fig. 4, can be used to constrain the size evolution, and subsequently the expansion velocity.In combination with the measurements of the photospheric velocity of the SN ejecta from spectroscopic observations, the SN size evolution can provide information on SN explosion models.We illustrate the size measurement via microlensing by the following toy model.
Let us assume that the brightness of a lensed SN image at any given time t and wavelength λ (ignoring time delays and macromagnification for simplicity) is: where S is the two dimensional intrinsic SN brightness profile, µ is the microlensing magnification on the source plane, and the integrals are performed over the extent of the source brightness profile.Assuming, for illustrative purposes, that the intrinsic surface brightness profile (S) is a circle of radius R(t) with a known, constant surface brightness i(t, λ), we get: in polar coordinates (r, θ), where the value of the integral can be written as a radius and time-dependent factor, F (R), and whose probability distribution can be calculated numerically from magnification maps.In fact, we can drop the assumption on a uniform profile and perform the convolution in Equation ( 6) for different profile shapes.Therefore, if we know S(x, y; t, λ) from a standardised SN type and calculate F , then we can measure R(t).If the SN brightness profile shape varies over time, ratios of I(t, λ) can be considered for separations in time ∆t small enough so that the profile does not change by much.Finally, we note that Dobler and Keeton (2006) have used a similar approach, the time-weighted light curve derivative, albeit to measure the stellar mass fraction of the lensing galaxy.
In an analogous way, the relative size, or more specifically the half-light radius, of the SN ejecta at different wavelengths can be estimated from Equation ( 6).Actually, the shape of the SN profile and/or the expansion velocity may vary in different wavelengths (see Figs. 11 and 12;Huber et al., 2019;Bayer et al., 2021), leading to different half-light radii (i.e., R = R(t, λ) instead of R(t)).In the context of our toy model, the brightness ratio between two wavelengths at a given time t can be used to constrain the relative half-light radii: Such measurements of radii across wavelengths and time not only provide a check on whether the expansion is homologous (when comparing the radii to the measured velocities from spectra), but also information on the explosion models that alter the dependence of radius on wavelength.It remains to be seen whether microlensing constraints on explosion models are competitive to the existing constraints based on SN spectral evolution modelling and analysis.

Dust properties in distant galaxies
The different sightlines of multipe images offer unique opportunities to explore the properties of the interstellar medium of the deflecting galaxy, as was first shown for strongly lensed QSOs by Falco et al. (1999) and in subsequent studies (e.g., Elíasdóttir et al., 2006;Østman et al., 2008;Hjorth et al., 2013).This is important as accurate measurements of e.g., the total-to-selective extinction parameter R V ≡ Fig. 12 Specific intensity profiles from TARDIS (Kerzendorf and Sim, 2014;Vogl et al., 2019) modelling of SN 1999em, a Type IIP SN (Vogl et al., 2020) projected on the plane of the sky (Bayer et al., 2021).Left: radial intensity profiles in different filters 11 rest-frame days after explosion, in units of the Einstein radius of the microlenses (R Ein = 2.9 × 10 16 cm for this specific case).Right panel: same as the left panel but for 27 rest-frame days after explosion.A V /E(B − V ), a key number for many areas of astronomy, are very hard to establish outside the Milky Way and the Magellanic Clouds.The extinction in the V-band, A V , is defined as where m V is the apparent magnitude with dust extinction and m V,0 is the intrinsic magnitude without dust extinction.The extinction for any other wavelength band is defined in a similar way.The color excess (or reddening) between the B and V band is In the case of no dust extinction, then A V = 0, A B = 0 and E(B − V ) = 0.For a fixed R V value, higher A V values (or higher E(B − V ) values) correspond to more dust along the sight line.The observed colors in multiple bands of the images 5000 7500 10000 12500 15000 17500 20000 Observed (Å) Image 1 Image 2 Image 3 Image 4 Host Galaxy Fig. 13 The observed color excess for the resolved images from HST of iPTF16geu as a function of wavelength.The absorption from the host galaxy dust grains is plotted in black.
For Image 1 we can see that the host galaxy is the dominant source of extinction, and for images 2, 3, 4 there is a progressively larger contribution from the dust in the lens galaxy, with correspondingly higher values of the color excess.are used to infer the differential reddening as the various images travel through different regions of the lensing galaxy (i.e., infer the differences in A V between the multiple SN images).Corrections for differential extinctions in systems of strongly lensed SNe Ia is crucial to also make use of their standard candle nature to infer the lensing magnification.
Exploiting the well-known color evolution of SNe Ia, Dhawan et al. (2020) were able to infer both the extinction in the lensing galaxy for iPTF16geu, along with the common color excess from attenuation in the host galaxy using multi-band images from HST, as shown in Fig. 13.For iPTF16geu, with an impact parameter of about 1 kpc, two of the images (3 and 4) suffered significant extinction and the total magnification inferred from spatially unresolved images was significantly underestimated.Intriguingly, the best fit values of R V for the lensing galaxy at z = 0.216 were much lower than what is observed for studies of stars in the Milky-Way (Schlafly et al., 2017), but similar to what is found in studies of well-measured, highly-reddened SNe Ia (Amanullah et al., 2015).

Methods
There are multiple ways to find lensed supernovae by making use of their image morphology, magnification, image multiplicity and/or time evolution.We briefly describe several approaches in this subsection.

Search through known lensed galaxies
A straightforward approach to find lensed SNe is to monitor known lensed galaxies and wait for a SN to explode in one of the lensed galaxies (Shu et al., 2018).Using the sample of 128 galaxy-scale strong-lens systems from the Sloan Lens ACS Survey (SLACS; Bolton et al., 2006Bolton et al., , 2008)), the SLACS for the Masses Survey (S4TM; Shu et al., 2017), and the Baryon Oscillation Spectroscopic Survey Emission-Line Lens Survey (BELLS; Brownstein et al., 2012;Bolton et al., 2012), Shu et al. (2018) estimated that the rates of strongly lensed SNe Ia and core-collapse SNe are 1.23 ± 0.12 and 10.4 ± 1.1 events per year, respectively.
One can either monitor known lens systems (especially ones with the highest star-formation and SN rates) through dedicated observing programs (Craig et al., 2021), or through wide-field imaging surveys such as the ongoing ZTF (Bellm et al., 2019).In the latter case, lens systems and lens candidates can be cross matched to transient alerts from ZTF through various brokers such as AMPEL (Nordin et al., 2019), ANTARES (Saha et al., 2014;Narayan et al., 2018;Lee et al., 2020;Matheson et al., 2021) and Lasair (Smith et al., 2019).Multi-epoch images of lens candidates can also be used to find lensed SNe based on temporal and spatial information of transients occurring near the lens candidates (Sheu et al., 2023).There are now thousands of confirmed and candidate lens systems from a wide range of lens searches, notably through machine learning approaches in recent years (e.g., Bolton et al., 2006;Limousin et al., 2009;Gavazzi et al., 2012;Brownstein et al., 2012;Vieira et al., 2013;Cañameras et al., 2015;Jacobs et al., 2017Jacobs et al., , 2019;;Marshall et al., 2016;More et al., 2016;Sonnenfeld et al., 2018Sonnenfeld et al., , 2020;;Petrillo et al., 2019;Cañameras et al., 2020Cañameras et al., , 2021;;Huang et al., 2020Huang et al., , 2021;;Savary et al., 2022;Rojas et al., 2022;Shu et al., 2022;Tran et al., 2022).However, given the relatively shallow limiting depth of r = 20.6 mag (5σ) of the ZTF survey (Bellm et al., 2019), individual lensed SN images associated with the lensed galaxies are likely below the detection threshold (Oguri and Marshall, 2010;Wojtak et al., 2019).Nonetheless, the combination of the flux from all lensed SN images can be above the detection threshold, which motivates the next search approach through magnifications.

Search through lensing magnification
While targeted searches for lensed SNe in known strongly lensed systems is straightforward, it is also constraining given the limited number of such systems.Since supernovae, especially SNe Ia, have a known narrow range of intrinsic brightness, it is feasible to identify particularly bright SNe as potential lensing candidates.Quimby et al. (2014) proposed a technique to find unresolved lensed images of SNe Ia based on their high magnification and colors.Lensed SNe are expected  to come from high redshifts where they will be in larger numbers owing to a larger volume and the lensing probability being higher for distant sources.The lensing magnification will make them appear brighter but their colors will largely be unaffected.Hence, in a color-magnitude diagram, the lensed SNe Ia will appear redder, at given magnitude, since they will be originating from higher redshifts compared to the unlensed population (see the blue circles above the black solid line in Fig. 14).
Another technique is based on having some knowledge of the redshifts of the galaxies involved, and it is interesting since it does not require multiple images to be spatially resolved.Thus it can in principle be used to detected arbitrarily compact systems, complementing other techniques, and allowing to probe the entire angular separation distribution of strongly lensed systems.Fig. 15 from Goobar et al. (2017) shows a practical realisation of detecting unresolved strongly lensed supernova: iPTF16geu was found to be a 30σ outlier when compared with other SNe Ia at the same redshift.Thus, even with the very modest ∼ 2 spatial resolution at the 48-inch telescope at the Palomar Observatory, a strongly lensed system with θ E = 0.3 could be identified.The difficulty with this approach is that spectroscopic observations are needed for classification of the transients, along with obtaining reliable redshifts of the galaxies, potentially a daunting task for large surveys with tens or even hundreds of new astrophysical transients found every night.Goldstein and Nugent (2017) proposed a methodology to lower the rate of false positives: only considering transients spatially associated with elliptical galaxies, which make up about ∼ 80% of the galaxy lenses and typically only host SNe Ia, including the sub-luminous ones.Furthermore, thanks to prominent 4000Å breaks, robust photometric redshifts can be computed for elliptical galaxies.Hence, if the transient appears to be too bright for being a SN Ia at the (photometric) redshift of the nearest galaxy, there is a good chance that it is a deflecting galaxy, as opposed to the host.Upon a secure identification, targeted high spatial-resolution imaging can be used to try to resolve the system.

Search through multiplicity of images
For imaging surveys with faint limiting depth and good angular resolution such as the upcoming LSST, the multiple SN images of a lensed SN can be resolved and detected individually.Oguri and Marshall (2010) and Wojtak et al. (2019) investigated the number of lensed SNe detectable from such an approach.Specifically, Wojtak et al. (2019) imposed the following conditions for detecting a lensed SN via its image multiplicity: (1) the maximum image separation between the multiple SN images is between 0.5 and 4 , where the lower limit is set by the expected seeing of LSST and the upper limit puts focus on galaxy-scale lens systems, (2) the flux ratio between the two images of a two-image (double) system is larger than 0.1, in order to have good contrast and clear identification of the images, and (3) at least three of the four images in a quad system are detected, and both images of a double system are detected.
For surveys with limiting magnitude that are fainter than ∼23 mag in g, r or i-band, the detection via image multiplicity is expected to detect more lensed SN systems than the magnification approach.We illustrate this in Fig. 16 for lensed SNe Ia, and refer to Wojtak et al. (2019) for other types of lensed SNe that show similar trends.A combination of the two complementary approaches (magnification and multiplicity), as labelled by "hybrid" in Fig. 16, delivers more lens systems than the individual approaches on their own.In terms of cosmological applications, the multiplicity approach tends to detect lensed SN systems with longer time delays and larger image separations, which are more suitable for timedelay cosmography.

Search through spatio-temporal images
Most of the lensed SN searches in ZTF are based on the approaches of cross matches to lens systems/candidates (Section 4.1.1)or magnification (Section 4.1.2)since the angular resolution of ZTF (with 1 pixel sizes) does not resolve the multiple SN images of typical lens systems.However, in the advent of LSST, there will be high angular-resolution and deep images with at least hundreds of epochs per filter at each sky location.Building upon the image multiplicity method described in Section 4. 1.3, Kodi Ramanah et al. (2021) proposed a machine-learning approach to classify an object based on images from cadenced surveys like the Young Supernova Experiment (YSE; Jones et al., 2021) and LSST.Machine learning techniques such as convolutional neural networks (CNN) perform well and are efficient at processing large amounts of imaging data for finding gravitational lenses (e.g., Jacobs et al., 2017Jacobs et al., , 2019;;Lanusse et al., 2018;Metcalf et al., 2019;Petrillo et al., 2019;Cañameras et al., 2020Cañameras et al., , 2021)).The spatio-temporal method of Kodi Multicolor light curve of iPTF16geu showing that the unresolved supernova was 4.3 magnitudes (30 standard deviations) brighter than expected for its redshift.The magnitudes are measured with respect to time of maximum light in the R-band at P48 and in the g-, r-, and i-bands with the P60 telescope.The solid lines show the best-fitted SN Ia model to the data, while the dashed lines indicate the expected light curves at z = 0.409 (without lensing).The bands represent the standard deviation of the brightness distribution for SNe Ia.To fit the observed light curves, a brightness boost from gravitational lensing of 4.3 magnitudes is required.Figure from Goobar et al. (2017).2021) builds upon the success of the CNN and incorporate also time-domain information.

Ramanah et al. (
Rather than providing a neural network with the multi-band stacked (static) images of objects for classification, Kodi Ramanah et al. ( 2021) developed a network to take in a temporal series of images.The network architecture is a CNN that encodes long short-term memory (LSTM, a type of recurrent network; Sher-

Expected number of lensed SN events
Strongly lensed transients are very rare.Hence, deep, wide-field time-domain optical and Near-IR surveys have the best chances to find samples of lensed SNe.Among these, the (optical) LSST survey at the Vera C. Rubin Observatory to see first light in 2024, and the (NIR) Roman satellite, planned for launch some years later offer the best opportunities.In both cases, between several tens to hundreds of strongly lensed supernovae may be expected throughout the multi-year surveys (Goobar et al., 2002;Oguri and Marshall, 2010;Quimby et al., 2014;Goldstein et al., 2019;Wojtak et al., 2019;Pierel et al., 2021).The cumulative number of SNe, Type Ia and core-collapse, expected to be found every year by LSST as a  function of the detection threshold is shown in Fig. 17.A potential concern is the low cadence of the observations, which may require follow-up observations to measure time delays accurately (Huber et al., 2019).

Expected lens properties
Fig. 18 shows simulations by Goldstein et al. (2019) of the distribution of system properties for strongly lensed SNe Ia in the LSST survey.The typical source redshift is z s ∼ 0.9 with time delays between 2 and 3 weeks.With a median image separation close to 1 and the expected good seeing conditions at the Rubin Observatory coupled with the plate-scale of the camera, most events would be spatially resolved.

Requirements on follow-up observations
Full scientific exploitation of strongly lensed SNe depends on 1) early transient detection, ideally well before light curve peak; 2) rapid spectroscopic classification of the transient; 3) well-sampled, spatially resolved light curves; and 4) good wavelength coverage, necessary for accurate extinction corrections.From Fig. 17 one can see that the vast majority of the lensed SNe in forthcoming surveys will have (summed) peak magnitudes fainter than 21 mag.Hence, spectroscopic classification requires 8-m class telescopes, in most cases, potentially a major challenge.As argued in Sec.3.4, an accurate measurement of the image magnifications requires a wide lever arm in wavelength coverage of the spatially resolved images, for at least a handful of epochs along the light curve.The new space telescope, JWST, would be ideally suited for follow-up.Its Near-IR sensitivity is an excellent match to the source redshifts shown in Fig. 18.

Summary
The first discoveries of strongly lensed SNe in recent years are opening a new window of exploration for cosmological and astrophysical studies.In this review, we have an overview of the analysis and results from these first lensed SN systems, and the future prospects for this exciting new field.The main takeaway points are as follows.
-There are currently five strongly lensed SN systems with spatially resolved SN images.Three systems are lensed by galaxy clusters and two lensed by individual galaxies.-Time delays between the multiple SN images of a lensed SN allow direct measurements of (1) the time-delay distance, (2) the angular-diameter distance to the lens in the case where stellar kinematic measurements of the lens galaxy are available, and (3) the luminosity distance in the case of a Type Ia SN that is a standardisable candle, as long as the SN magnifications due to microlensing and millilensing can be accurately accounted for.These distance measurements provide competitive constraints on cosmological parameters, particularly H 0 .-Several new methods to infer the time delays of lensed SNe have been developed in recent years, employing different kinds of data such as light curves, color curves or spectral evolution of SNe.Time-delay measurements with uncertainties of ∼1 day are achievable with real and mock data.The fractional uncertainty in the delays contribute directly to the fractional uncertainty on the time-delay distance and the lens angular-diameter distance.-Cosmographic constraints, especially the Hubble constant inference, from SN Refsdal are forthcoming.We expect that a modest sample of 20 lensed SNe Ia from the upcoming LSST to yield H 0 with nearly 1% precision in flat ΛCDM cosmology.
-The time delays of lensed SNe provide a unique opportunity to acquire very early-phase observations of SNe, especially in the rest-frame UV, and probe SN progenitors.-Lensing magnifications allow the acquisition of spectra of SN above z > 1, that are crucial to reduce systematic uncertainties in SN Ia cosmology.-Microlensing of SNe provides an avenue to constrain the sizes of the SN at different wavelengths.-The dust properties in the (foreground lens) galaxies can be measured through the multiple SN sight lines in lensed SN systems.-Various methods have been developed to search for lensed SNe, through e.g., the monitoring of known lens systems, magnification, image multiplicity, and spatio-temporal evolution in the imaging.Depending on the selection criteria, we expect at least dozens of lensed SNe to explode in the upcoming LSST, with their properties summarized in Section 4.3.
We are entering a new era of lensed SNe with the upcoming wide-field cadenced imaging surveys.Rapid follow-up observations including spectroscopic typing and light curve monitoring, will be crucial and necessary to make the best use of such events for cosmological and astrophysical studies.

Fig. 2
Fig.2Wide-field image from iPTF showing the portion of the CCD camera where iPTF16geu was found in modest seeing (2 ).The insets show how the strong lensing nature, a quadruple lens with an Einstein radius of only 0.3 , could be verified using HST imaging in the optical, and laser-guide-star adaptive-optics imaging from Keck in the Near-IR.Image credit: J. Johansson.
obtained HST imaging of this system, showing clearly the four multiple images of the SN Ia in a symmetric configuration.Based on the single-epoch photometry of the SN Ia images, Pierel et al. (2022) measured short time delays of < 1 day between the multiple images, which are consistent with the predictions from their multiple lens mass models obtained using different lensing software.Both Goobar et al. (2022) and Pierel et al. (2022) find anomalous flux ratios of the SN images compared to the smooth mass model predictions, indicating the presence of dust, millilensing and/or microlensing.The fifth system of spatially-resolved lensed SN is reported very recently by Chen et al. (2022).A SN was lensed by one of the galaxies in the galaxy cluster Abell 370 in 2010 into three multiple images, and Chen et al. (2022) discovered

Fig. 3
Fig. 3 Panel a): Color image of SN Requiem showing arc-like images of the distant host galaxy (H1-H4), the three SN images (SN1-SN3) and an ellipse pointing out the expected location of SN4.Panels b-i): zoomed-in regions around the SN locations from the data taken in 2016 (b-e) and 2019 (f-i).Figure taken from Rodney et al. (2021).

Fig. 4
Fig. 4 Examples of Type Ia SNe light curves with significant microlensing.Left panel: a map of microlensing magnifications (µ, indicated by the color bar) on the source plane as a result from stars in the foreground lens galaxy.The map is for a lensed SN image with convergence κ = 0.6, shear γ = 0.6, and smooth dark matter fraction s = 0.6.In this example, the Einstein radius for the mean microlens is R Ein = 7.2 × 10 −3 pc.The two circles in solid cyan and dashed magenta indicate the size of a SN Ia at 21 rest-frame days after explosion.Middle and right panels: microlensed light curves in the g-band (middle panel) and z-band (right panel) corresponding to the SN positions shown in the left panel (solid cyan and dashed magenta).The intrinsic light curve without microlensing is shown in dotted black.When the expanding photosphere of the SN crosses a microlensing caustic with high µ, its light curve changes substantially relative to the no-microlensing case.Figure taken from Huber et al. (2019).

Fig. 5
Fig. 5 Microlensing induced scatter (indicated by the color bar) for a point source lensed by a Singular Isothermal Sphere macromodel as a function of stellar fraction and macromagnification.Scatter is lower for the lensed SN image that is far from the Einstein radius (with lower convergence/shear), but there is a limit since a convergence below 0.25 produces only single images (i.e., no strong lensing) and convergence below ∼0.35 produces a faint counterimage (close to lens center) that makes time-delay measurements more challenging.Image credit: Luke Weisenbach, modified based on Weisenbach et al. (2021).

Fig. 6
Fig. 6 Light curves of SN Refsdal from Hubble Space Telescope imaging in multiple wavelength filters.Each row consists of the light curves obtained in the specific filter that is indicated in the leftmost panel.Each column shows one of the four SN images, S1-S4 (left to right), that are indicated on the top panels.Each panel shows the observed AB magnitude as a function of the observer-frame days.Figure taken from Rodney et al. (2016).

Fig. 7
Fig.7Spectral evolution of core-collapse SN1999em based on the TARDIS simulation fromVogl et al. (2019Vogl et al. ( , 2020)).Three prominent spectral features, Hβ, FeII and Hα, are labelled.The SN phases are indicated by their rest-frame days after explosion.As the SN phase increases, the spectral features become stronger and the absorption wavelengths increase.Such a sequence of spectra of the first-appearing SN image provides the wavelength-phase relation, and a measurements of the absorption wavelength of a trailing SN image therefore provides information of its SN phase and thus its time delay relative to the first-appearing SN image.With each spectrum of signal-to-noise of 20,Bayer et al. (2021) showed that the time delays can be measured with uncertainties of ∼2 days per spectral feature, even after accounting for the effects of microlensing.Figure taken from Bayer et al. (2021).

Fig. 8
Fig. 8 Forecast of H 0 from SN Refsdal as a function of the time delay between SN images SX and S1, using the reference mass model of Grillo et al. (2020).The median values of H 0 (diamonds) with the 1σ uncertainties in flat ΛCDM are shown for 9 different hypothetical SX-S1 delays, each with an assumed uncertainty of 10 days.For comparison, the blue, magenta and red bands show, respectively, the 1σ credible intervals from SH0ES (Riess et al. 2019), H0LiCOW (Wong et al. 2020) and Planck (Planck Collaboration et al. 2020).Figure taken from Grillo et al. (2020).

Fig. 10
Fig. 10 Expected signal for HST observations through grisms G102 and G141 of a SN Ia at z = 1.5 for different magnifications provided by lensing.The shaded regions in [orange, green, blue, red] indicate the 1σ uncertainty per pixel in one-hour-long exposures for ∆m = [1, 1.5, 2.0, 2.5] mag of magnification, respectively.The observed spectrum of "SN Primo" (Rodney et al., 2012, six hours exposure time) is also shown for comparison.

Fig. 11
Fig. 11 Specific intensity profile of the W7 explosion model of SNe Ia.Left: radial intensity profiles in different filters 14.9 rest-frame days after explosion, in units of the Einstein radius of the microlenses (R Ein = 2.2 × 10 16 cm for this specific case).The vertical solid cyan lines indicate the radius that encloses 99.9% of the total projected specific intensity.The vertical black dashed lines are random locations of caustics -effects of microlensing are strong when the specific intensity of the SN crosses a caustic.Right panel: same as the left panel but for a SN Ia 39.8 rest-frame days after explosion.Figure extracted and modified from Huber et al. (2019).
Fig.12Specific intensity profiles from TARDIS(Kerzendorf and Sim, 2014;Vogl et al., 2019) modelling of SN 1999em, a Type IIP SN(Vogl et al., 2020) projected on the plane of the sky(Bayer et al., 2021).Left: radial intensity profiles in different filters 11 rest-frame days after explosion, in units of the Einstein radius of the microlenses (R Ein = 2.9 × 10 16 cm for this specific case).Right panel: same as the left panel but for 27 rest-frame days after explosion.Figure taken from Bayer et al. (2021).
Fig.13The observed color excess for the resolved images from HST of iPTF16geu as a function of wavelength.The absorption from the host galaxy dust grains is plotted in black.For Image 1 we can see that the host galaxy is the dominant source of extinction, and for images 2, 3, 4 there is a progressively larger contribution from the dust in the lens galaxy, with correspondingly higher values of the color excess.Figure taken from Dhawan et al. (2020).

Figure 4 :
Figure 4: Color-magnitude diagram showing how lensed SNIa can be distinguished from un-lensed events.The blue shaded area shows the expected color-magnitude distribution for un-lensed SNIa on a log scale, and the red shaded area corresponds to core-collapse supernovae.The r − i colors for low redshift supernovae are relatively blue.However, at higher redshifts (fainter observed magnitudes), the color becomes red as the peak of the rest-frame spectral energy distribution passes through the observer-frame bands.The red limit for unlensed supernovae at a given i-band magnitude is denoted by the thick black line.Blue circles and red triangles show the distribution of lensed SNIa and core-collapse supernovae, respectively, predicted by Monte Carlo simulations (16).Filled symbols indicate objects that could be resolved from ground based observations, such as those planned by the Large Synoptic Survey Telescope (LSST).Open symbols depict objects that require high angular resolution follow-up observations to resolve spatially.The open star marks the values corresponding to the peak iband brightness of PS1-10afx, and the dash-dotted curve shows that the color evolution within one magnitude of this peak is minimal.The vertical dashed line marks the single epoch limit predicted for LSST.The arrow shows the reddening vector, assuming A V = 1.0 mag.

Fig. 17
Fig.17Cumulative number of lensed Type Ia and core-collapse supernovae detections by the LSST survey per year as a function of their peak i-band magnitude, as computed inGoldstein et al. (2019).Figure credit: Ana Sagués Carracedo.
Fig.17Cumulative number of lensed Type Ia and core-collapse supernovae detections by the LSST survey per year as a function of their peak i-band magnitude, as computed inGoldstein et al. (2019).Figure credit: Ana Sagués Carracedo.
Suyu et al. (2020)y ∼650 systems would be detected early enough by LSST to deliver reliable time delays.This implies a sam-Forecast of constraints on H 0 and Ωm in flat ΛCDM from a sample of 20 lensed SNe Ia from LSST with precise and accurate time-delay measurements.Assuming distance measurements of D ∆t and D d with 6.6% and 5% uncertainties, respectively, for each lensed SN Ia, this modest sample is expected to yield H 0 and Ωm with precisions of 1.3% and 19%, respectively.Figure taken fromSuyu et al. (2020).
The number of lensed SNe Ia expected per year over the whole sky as a function of the imaging survey limiting magnitude depth.The different colors correspond to detection in different filters (g, r, i, z or y).The dotted (dashed) lines show the expected number of lensed SNe Ia detected through the magnification (image multiplicity) approach.The hybrid approach, a combination of magnification and image multiplicity, is indicated by the solid curve.For shallow image surveys with limiting depth brighter than ∼22 such as the ZTF survey, the magnification approach dominates in providing most of the expected lensed SNe.For deep image surveys with limiting depth fainter than ∼23 such as the LSST survey, the image multiplicity approach start to dominate.Figure taken from Wojtak et al. (2019).stinsky,2020).The concept is that the temporal series of images will show the multiple SN images appearing at different epochs.Even if the multiple SN images are not well resolved, the centroids of the distribution of light would change as the multiple SN images brighten and dim at different times.Such change in the features of the object helps the neural network to distinguish lensed SN (with multiple SN images) from non-lensed SN (with single SN image).Using simulated images ofYSE, Kodi Ramanah et al. (2021)demonstrated that the spatio-temporal network improves the classification accuracy by ∼20% compared to networks that use static (e.g., single-epoch) images.The new development spatio-temporal network is very promising for application to the LSST.