Professor Stephen Smartt CBE FRS MRIA

Abstract:

We present the Virtual Research Assistant (VRA) of the ATLAS sky survey, which performs preliminary eyeballing on our clean transient data stream. The VRA uses histogram-based gradient-boosted decision tree classifiers trained on real data to score incoming alerts on two axes: “Real” and “Galactic.” The alerts are then ranked using a geometric distance such that the most “real” and “extragalactic” receive high scores; the scores are updated when new lightcurve data is obtained on subsequent visits. To assess the quality of the training we use the recall at rank K, which is more informative to our science goal than general metrics (e.g., accuracy, F1-scores). We also establish benchmarks for our metric based on the pre-VRA eyeballing strategy, to ensure our models provide notable improvements before being added to the ATLAS pipeline. Then, policies are defined on the ranked list to select the most promising alerts for humans to eyeball and to automatically remove bogus alerts. In production the VRA method has resulted in a reduction in eyeballing workload by 85% with a loss of follow-up opportunity <0.08%. It also allows us to automatically trigger follow-up observations with the Lesedi telescope, paving the way toward automated methods that will be required in the era of LSST. Finally, this is a demonstration that feature-based methods remain extremely relevant in our field, being trainable on only a few thousand samples and highly interpretable; they also offer a direct way to inject expertise into models through feature engineering.

Abstract:

Context. Type IIb supernovae (SNe IIb) often exhibit an early light curve excess (EE) preceding the main peak powered by ⁵⁶Ni decay. The physical origin of this early emission remains an open question. Among the proposed scenarios, shock cooling (SC) emission, resulting from the interaction of the shock wave with extended envelopes, is considered the most plausible mechanism. However, the occurrence rate of such events has yet to be reliably constrained. Aims. This study aims to quantify the frequency of EE in SNe IIb and investigate its physical origin by analysing optical light curves from the Asteroid Terrestrial-impact Last Alert System (ATLAS) survey, ultimately providing qualitative constraints on their progenitor systems. Methods. We identified 74 potential SNe IIb from 153 spectroscopically classified events reported in the Transient Name Server (TNS), observed by ATLAS with peak fluxes exceeding 150 μJy (18.46 mag) and explosion epoch uncertainties below six days. Using a spectral reclassification method, we selected a sample of 66 SNe IIb and a cleaned sample of 59 SNe IIb for analysis. We then applied light curve model fitting and outlier analysis to identify objects exhibiting EE emission and studied their photometric properties. Results. We identify 20 SNe IIb with EE, corresponding to a frequency of approximately 30.5% to 50%, the higher value being obtained under the most stringent observational data-quality cuts. The duration and colour evolution of the early excess support its interpretation as shock cooling in extended envelopes. We also find that EE SNe IIb exhibit faster post-peak declines than non-EE events, while both groups show similar peak absolute magnitudes and rise-time distributions. Conclusions. Our findings suggest that EE and non-EE SNe IIb likely share similar initial progenitor masses but differ in their ejecta properties, potentially due to varying degrees of binary interaction. This study constrains EE SNe frequency and photometric properties, paving the way for future theoretical work, such as hydrodynamical modelling of EE SNe light curves, which could corroborate these results and contribute to constraining the evolutionary pathways of SNe IIb progenitor systems.

Abstract:

We discuss the results of the spectroscopic and photometric monitoring of the type IIn supernova (SN) 2023ldh. Survey archive data show that the SN progenitor experienced erratic variability in the years before exploding. Beginning May 2023, the source showed a general slow luminosity rise that lasted for over four months, with some superposed luminosity fluctuations. In analogy to SN 2009ip , we call this brightening ‘Event A’. During Event A, SN 2023ldh reached a maximum absolute magnitude of M r = −15.52 ± 0.24 mag. The light curves then decreased by about 1 mag in all filters for about two weeks reaching a relative minimum, which was followed by a steep brightening (Event B) to an absolute peak magnitude of M r = −18.53 ± 0.23 mag, replicating the evolution of SN 2009ip and similar to that of type IIn SNe. The three spectra of SN 2023ldh obtained during Event A show multi-component P Cygni profiles of H I and Fe II lines. During the rise to the Event B peak, the spectrum shows a blue continuum dominated by Balmer lines in emission with Lorentzian profiles, with a full width at half maximum velocity of about 650 km s −1 . Later, in the post-peak phase, the spectrum reddens, and broader wings appear in the H α line profile. Metal lines with P Cygni profiles and velocities of about 2000 km s −1 are clearly visible. Beginning around three months past maximum and until very late phases, the Ca II lines become among the most prominent features, while H α is dominated by an intermediate-width component with a boxy profile. Although SN 2023ldh mimics the evolution of other SN 2009ip -like transients, it is slightly more luminous and has a slower photometric evolution. The surprisingly homogeneous observational properties of SN 2009ip -like events may indicate similar explosion scenarios and similar progenitor parameters.

Abstract:

Abstract With the 4-meter Multi-Object Spectroscopic Telescope (4MOST) expected to provide an influx of transient spectra when it begins observations in early 2026 we consider the potential for real-time classification of these spectra. We investigate three extant spectroscopic transient classifiers: the Deep Automated Supernova and Host classifier (DASH), Next Generation SuperFit (NGSF) and SuperNova IDentification (SNID), with a focus on comparing the completeness and purity of the transient samples they produce. We manually simulate fibre losses critical for accurately determining host-contamination and use the 4MOST Exposure Time Calculator to produce realistic, 4MOST-like, host-galaxy contaminated spectra. We investigate the three classifiers individually and in all possible combinations. We find that a combination of DASH and NGSF can produce a SN Ia sample with a purity of 99.9 per cent while successfully classifying 70 per cent of SNe Ia. However, it struggles to classify non-SN Ia transients. We investigate photometric cuts to transient magnitude and the transient’s fraction of total fibre flux, finding that both can be used to improve non-SN Ia transient classification completeness by 8–44 per cent with SNe Ibc benefitting the most and superluminous (SL) SNe the least. Finally, we present an example classification plan for live classification and the predicted purities and completeness across five transient classes: Ia, Ibc, II, SL and non-SN transients. We find that it is possible to classify 75 per cent of input spectra with >70 per cent purity in all classes except non-SN transients. Precise values can be varied using different classifiers and photometric cuts to suit the needs of a given study.

Abstract:

We present ultraviolet to infrared observations of the extraordinary Type IIn supernova 2023zkd (SN 2023zkd). Photometrically, it exhibits persistent and luminous precursor emission spanning ∼4 yr preceding discovery (Mr ≈ −15 mag, 1500 days in the observer frame), followed by a secondary stage of gradual brightening in its final year. Post-discovery, it exhibits two photometric peaks of comparable brightness (Mr ≲ −18.7 mag and Mr ≈ −18.4 mag, respectively) separated by 240 days. Spectroscopically, SN 2023zkd exhibits highly asymmetric and multicomponent Balmer and He I profiles that we attribute to ejecta interaction with fast-moving (1000–2000 km s−1) He-rich polar material and slow-moving (∼400 km s−1) equatorially distributed H-rich material. He II features also appear during the second light curve peak and evolve rapidly. Shock-driven models fit to the multiband photometry suggest that the event is powered by interaction with ∼5–6 M⊙ of CSM, with 2–3 M⊙ associated with each light curve peak, expelled during mass-loss episodes ∼3–4 yr and ∼1–2 yr prior to explosion. The observed precursor emission, combined with the extreme mass-loss rates required to power each light curve peak, favors either super-Eddington accretion onto a black hole or multiple long-lived eruptions from a massive star to luminosities that have not been previously observed. We consider multiple progenitor scenarios for SN 2023zkd, and find that the brightening optical precursor and inferred explosion properties are most consistent with a massive (MZAMS ≥ 30 M⊙) and partially stripped He star undergoing an instability-induced merger with a black hole companion.