The Legacy Survey of Space and Time (LSST) will start in late-summer 2026, revolutionizing transient astronomy. Here, we present the Dark Energy Camera (DECam) Shadow Survey, which is designed to maximize the science potential of LSST by shadowing LSST observations of local galaxy-cluster fields, producing a nightly cadence of these fields. The Shadow Survey will discover extremely young supernovae (SNe), SN precursors, as well as other explosive transients and exotic phenomena, helping to characterize such transients at unprecedented cadence and depth when combined with LSST. We describe our workflow, pipeline, public data releases, and candidate vetting. As an early result of Shadow, we present the fitful luminous blue variable (LBV) eruptions of AT2017des in the Virgo-Cluster galaxy NGC4532. AT2017des has short-timescale variability (of order 10 days), peaking at around M_r=-12.5mag, brighter than normal LBVs, and similar to the more extreme flaring of hot LBVs/SN impostors such as SN2000ch, AT2016blu, and the precursor activity of SN2009ip. Our spectral time-series reveals features typical of these hot LBVs and SN impostors/precursors. Combining our data with long-baseline photometry from additional observatories, we find that the peaks of the outbursts of AT2017des are getting brighter over time, with 2026 peak fluxes being up to 5 times greater than in 2023 and an average brightening of ∼0.05 mag yr^-1. The peaks of AT2017des are more luminous than those of most other LBVs, only being fainter than bright precursors such as SN2009ip, and extreme SN impostors such as AT2016blu. AT2017des may therefore be “ramping up” to a terminal explosion.

Time-domain and multimessenger astronomy (MMA/TDA) targets demand rapid-response follow-up observations. In many cases, it is the only way to make discoveries and advance our understanding of the astrophysical phenomena, for example, kilonovae accompanying gravitational waves from compact object mergers, shock breakout in supernovae, prompt emission from GRBs, etc. Presently the MMA/TDA follow-up workflow requires wrangling disparate software packages and user interfaces. We present an end-to-end software tool for the community, the Gemini Observation and Analysis of Targets System (GOATS), which unifies and simplifies the workflow, particularly for Gemini follow-up observations. GOATS achieves this by integrating services from Gemini Observatory and its parent organization, NSF NOIRLab. From a single platform, GOATS enables enhanced target selection via NOIRLab's ANTARES alert broker, triggering of Gemini (and other facilities within the Astronomical Event Observatory Network), automated data retrieval from the Gemini Observatory Archive, and interactive data reduction and analysis through Gemini's DRAGONS software and NOIRLab's Astro Data Lab science platform. GOATS was successfully deployed in an end-to-end demonstration of real-time follow-up of Rubin/LSST alerts with NOIRLab facilities. As part of this demonstration, we selected targets from the Rubin alert stream and triggered follow-up observations within minutes of the Rubin detections. We obtained spectra for several targets and classified them as supernova of various types (Ia, IIP, Ib/c) with redshifts ranging from 0.05 to 0.35. By eliminating the need to manually connect tools and automating repetitive tasks, GOATS lowers the entry barrier and allows users to focus on the scientific interpretation of the observation results.

We present a synopsis of the project to establish 32 new faint (16.5 <= V <= 19.8) DA white dwarf (DAWD) spectrophotometric standards distributed over the whole sky. Our results validate the use of fully radiative pure hydrogen model fluxes for hot DA white dwarfs to predict the observed broadband fluxes from near-ultraviolet through the near-infrared to accuracies of a few parts per thousand. After fitting the line of sight reddenings simultaneously with the model spectral energy distributions of these stars against spectroscopic and multiband photometric observations, we have shown that residuals have an rms of typically 0.4%. This indicates that the complications from interstellar dust extinction have been adequately mitigated. Our stars supplement the three brighter DAWDs that define the flux scale of CALSPEC. The consequent photometric accuracy, their all-sky coverage, and their brightness range that matches the dynamic range of large telescopes, constitute an unprecedented ensemble of standard stars for both ground as well as space based use. This paper targets readers who may wish to use these as standard stars and provides them with the essential content to understand their strengths and limitations, without traversing the technical details of analysis that are already captured in a series of papers since 2016. The narrative here describes the motivation, justification, and evolution of the analysis methods; the input data that constrain the modeling; as well as the stability of our results in the face of future improvements in models.

Millions of supernovae will be discovered with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). As a result, spectrographs around the world will have to make difficult decisions about which supernova candidates receive spectroscopic follow-ups. This work identifies the minimum spectral resolution, R_λ = Δλ, as a function of signal-to-noise ratio (SNR) at which spectral classification of supernova subtypes becomes impossible. We include supernova types Ia, Ia-91T, Ia-91bg, Iax, Ib, Ic, broad-lined Ic, IIb, IIP, and Ibn in this work. We produce a definition of SNR based on specific lines for each SN subtype that allows us to generate homogeneous datasets at 16 different values of R_λ and 14 different SNR's and we tested the classification performance of a recently developed deep-learning classifier, ABC-SN, on each R_λ and SNR combination. We find that classification of supernova spectra into a refined taxonomy that separates, for example, between different subtypes of stripped envelope supernovae, is possible at low resolution and low SNR with no loss in model performance down to R_λ = 50 and SNR = 5. Classification performance is only minimally impacted even as low as R_λ = 25. We hope that astronomers using the LSST alert stream, as well as designers of future instruments and observatories, will benefit from knowing what spectral resolution is necessary to classify a supernova for arbitrary .

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (16.5 < V < 19.5) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zeropoints and WD parameters (surface gravity log g, effective temperature T_eff, extinction A_V, dust relation parameter R_V) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy and ground-based optical spectroscopy to sub-percent precision. Photometric residuals for the sample are the lowest yet yielding <0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory's Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

Fainter standard stars are essential for the calibration of larger telescopes. This work adds to the calibration spectra (CALSPEC) database of 19 faint white dwarfs (WDs) with all-sky coverage and V magnitudes between 16.5 and 18.7. Included for these stars is a new ultraviolet (UV) Hubble Space Telescope (HST) Space Telescope Imaging Spectrograph between 1150 and 3000 & Aring; with a resolution of similar to 500. Pure hydrogen WD models are fit to these UV spectra and to six-band HST/Wide Field Camera 3 photometry at 0.28-1.6 mu m to construct predicted model spectral energy distributions covering wavelengths from 900 & Aring; to the James Webb Space Telescope limit of 30 mu m using well-established CALSPEC procedures for producing flux standards with the goal of 1% accuracy.

Light echoes give us a unique perspective on the nature of supernovae and nonterminal stellar explosions. Spectroscopy of light echoes can reveal details on the kinematics of the ejecta, probe asymmetry, and reveal details of ejecta interaction with circumstellar matter, thus expanding our understanding of these transient events. However, the spectral features arise from a complex interplay between the source photons, the reflecting dust geometry, and the instrumental setup and observing conditions. In this work, we present an improved method for modeling these effects in light echo spectra, one that relaxes the simplifying assumption of a light-curve-weighted sum, and instead estimates the true relative contribution of each phase of a transient event to the observed spectrum. We discuss our logic, the gains we obtain over light echo analysis methods used in the past, and prospects for further improvements. Lastly, we show how the new method improves our analysis of echoes from Tycho's supernova (SN 1572) as an example.

We introduce a pipeline that performs rapid image subtraction and source selection to detect transients, with a focus on identifying gravitational-wave optical counterparts using the Dark Energy Camera (DECam). In this work, we present the pipeline steps from processing raw data to identification of astrophysical transients on individual exposures. We process DECam data and build difference images using the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Science Pipelines software, and we use flags and principal component analysis to select transients on a per-exposure basis, without associating the results from different exposures. Those candidates will be sent to brokers for further classification and alert distribution. We validate our pipeline using archival exposures that cover various types of objects, and the tested targets include a kilonova (GW170817), supernovae, stellar flares, variable stars (in a resolved galaxy or the Milky Way Bulge), and serendipitous objects. Overall, the data processing produces clean light curves that are comparable with published results, demonstrating the photometric quality of our pipeline. Real transients can be well selected by our pipeline when sufficiently bright (signal-to-noise ratio greater than or similar to 15). This pipeline is intended to serve as a tool for the broader research community. Although this pipeline is designed for DECam, our method can be easily applied to other instruments and future LSST observations.

We follow our first paper with an analysis of the ensemble of the extensive preexplosion ground- and space-based infrared observations of the red supergiant (RSG) progenitor candidate for the nearby core-collapse supernova SN 2023ixf in Messier 101, together with optical data prior to the explosion obtained with the Hubble Space Telescope (HST). We have confirmed the association of the progenitor candidate with the supernova (SN), as well as constrained the metallicity at the SN site, based on SN observations with instruments at Gemini-North. The internal host extinction to the SN has also been confirmed from a high-resolution Keck spectrum. We fit the observed spectral energy distribution (SED) for the star, accounting for its intrinsic variability, with dust radiative-transfer modeling, which assumes a silicate-rich dust shell ahead of the underlying stellar photosphere. The star is heavily dust obscured, likely the dustiest progenitor candidate yet encountered. We found median estimates of the star’s effective temperature and luminosity of 2770 K and 9.0 × 10 ^4 L _⊙ , with 68% credible intervals of 2340–3150 K and (7.5–10.9) × 10 ^4 L _⊙ , respectively. The candidate may have a Galactic RSG analog, IRC −10414, with a strikingly similar SED and luminosity. Via comparison with single-star evolutionary models we have constrained the initial mass of the progenitor candidate from 12 M _⊙ to as high as 14 M _⊙ . We have had available to us an extraordinary view of the SN 2023ixf progenitor candidate, which should be further followed up in future years with HST and the James Webb Space Telescope.

We present LAISS (Lightcurve Anomaly Identification and Similarity Search), an automated pipeline to detect anomalous astrophysical transients in real-time data streams. We deploy our anomaly detection model on the nightly ZTF Alert Stream via the ANTARES broker, identifying a manageable $\sim$1-5 candidates per night for expert vetting and coordinating follow-up observations. Our method leverages statistical light-curve and contextual host-galaxy features within a random forest classifier, tagging transients of rare classes (spectroscopic anomalies), of uncommon host-galaxy environments (contextual anomalies), and of peculiar or interaction-powered phenomena (behavioral anomalies). Moreover, we demonstrate the power of a low-latency ($\sim$ms) approximate similarity search method to find transient analogs with similar light-curve evolution and host-galaxy environments. We use analogs for data-driven discovery, characterization, (re-)classification, and imputation in retrospective and real-time searches. To date we have identified $\sim$50 previously known and previously missed rare transients from real-time and retrospective searches, including but not limited to: SLSNe, TDEs, SNe IIn, SNe IIb, SNe Ia-CSM, SNe Ia-91bg-like, SNe Ib, SNe Ic, SNe Ic-BL, and M31 novae. Lastly, we report the discovery of 325 total transients, all observed between 2018-2021 and absent from public catalogs ($\sim$1% of all ZTF Astronomical Transient reports to the Transient Name Server through 2021). These methods enable a systematic approach to finding the "needle in the haystack" in large-volume data streams. Because of its integration with the ANTARES broker, LAISS is built to detect exciting transients in Rubin data.

In this White Paper, we present recommendations for the scientific community and funding agencies to foster the infrastructure for a collaborative multi-messenger and time-domain astronomy (MMA/TDA) ecosystem. MMA/TDA is poised for breakthrough discoveries in the coming decade. In much the same way that expanding beyond the optical bandpass revealed entirely new and unexpected discoveries, cosmic messengers beyond light (i.e., gravitational waves, neutrinos, and cosmic rays) open entirely new windows to answer some of the most fundamental questions in (astro)physics: heavy element synthesis, equation of state of dense matter, particle acceleration, etc. This field was prioritized as a frontier scientific pursuit in the 2020 Decadal Survey on Astronomy and Astrophysics via its "New Windows on the Dynamic Universe" theme. MMA/TDA science presents technical challenges distinct from those experienced in other disciplines. Successful observations require coordination across myriad boundaries -- different cosmic messengers, ground vs. space, international borders, etc. -- all for sources that may not be well localized, and whose brightness may be changing rapidly with time. Add that all of this work is undertaken by real human beings, with distinct backgrounds, experiences, cultures, and expectations, that often conflict. To address these challenges and help MMA/TDA realize its full scientific potential in the coming decade (and beyond), the second in a series of community workshops sponsored by the U.S. National Science Foundation (NSF) and NASA titled "Windows on the Universe: Establishing the Infrastructure for a Collaborative Multi-Messenger Ecosystem" was held on October 16-18, 2023 in Tucson, AZ. Here we present the primary recommendations from this workshop focused on three key topics -- hardware, software, and people and policy. [abridged]

We report early-time ultraviolet (UV) and optical spectroscopy of the young, nearby Type II supernova (SN) 2022wsp obtained by the Hubble Space Telescope (HST)/STIS at about 10 and 20 days after the explosion. The SN 2022wsp UV spectra are compared to those of other well-observed Type II/IIP SNe, including the recently studied Type IIP SN 2021yja. Both SNe exhibit rapid cooling and similar evolution during early phases, indicating a common behavior among SNe II. Radiative-transfer modeling of the spectra of SN 2022wsp with the TARDIS code indicates a steep radial density profile in the outer layer of the ejecta, a supersolar metallicity, and a relatively high total extinction of E(B-V) = 0.35 mag. The early-time evolution of the photospheric velocity and temperature derived from the modeling agree with the behavior observed from other previously studied cases. The strong suppression of hydrogen Balmer lines in the spectra suggests interaction with a pre-existing circumstellar environment could be occurring at early times. In the SN 2022wsp spectra, the absorption component of the Mg II P Cygni profile displays a double-trough feature on day +10 that disappears by day +20. The shape is well reproduced by the model without fine-tuning the parameters, suggesting that the secondary blueward dip is a metal transition that originates in the SN ejecta.

This paper presents a new optical imaging survey of four deep drilling fields (DDFs), two Galactic and two extragalactic, with the Dark Energy Camera (DECam) on the 4 meter Blanco telescope at the Cerro Tololo Inter-American Observatory (CTIO). During the first year of observations in 2021, $>$4000 images covering 21 square degrees (7 DECam pointings), with $\sim$40 epochs (nights) per field and 5 to 6 images per night per filter in $g$, $r$, $i$, and/or $z$, have become publicly available (the proprietary period for this program is waived). We describe the real-time difference-image pipeline and how alerts are distributed to brokers via the same distribution system as the Zwicky Transient Facility (ZTF). In this paper, we focus on the two extragalactic deep fields (COSMOS and ELAIS-S1), characterizing the detected sources and demonstrating that the survey design is effective for probing the discovery space of faint and fast variable and transient sources. We describe and make publicly available 4413 calibrated light curves based on difference-image detection photometry of transients and variables in the extragalactic fields. We also present preliminary scientific analysis regarding Solar System small bodies, stellar flares and variables, Galactic anomaly detection, fast-rising transients and variables, supernovae, and active galactic nuclei.

Observational evidence points to a red supergiant (RSG) progenitor for SN 2023ixf. The progenitor candidate has been detected in archival images at wavelengths (>0.6 micron) where RSGs typically emit profusely. This object is distinctly variable in the infrared (IR). We characterize the variability using pre-explosion mid-IR (3.6 and 4.5 micron) Spitzer and ground-based near-IR (JHKs) archival data jointly covering 19 yr. The IR light curves exhibit significant variability with RMS amplitudes in the range of 0.2-0.4 mag, increasing with decreasing wavelength. From a robust period analysis of the more densely sampled Spitzer data, we measure a period of 1091+/-71 days. We demonstrate using Gaussian Process modeling that this periodicity is also present in the near-IR light curves, thus indicating a common physical origin, which is likely pulsational instability. We use a period-luminosity relation for RSGs to derive a value of M_K=-11.58+/-0.31 mag. Assuming a late M spectral type, this corresponds to log(L/L_sun)=5.27+/-0.12 at T_eff=3200 K and to log(L/L_sun)=5.37+/-0.12 at T_eff=3500 K. This gives an independent estimate of the progenitor's luminosity, unaffected by uncertainties in extinction and distance. Assuming the progenitor candidate underwent enhanced dust-driven mass-loss during the time of these archival observations, and using an empirical period-luminosity-based mass-loss prescription, we obtain a mass-loss rate of around (2-4)x10^-4 M_sun/yr. Comparing the above luminosity with stellar evolution models, we infer an initial mass for the progenitor candidate of 20+/-4 M_sun, making this one of the most massive progenitors for a Type II SN detected to-date.

Hot DA white dwarfs have fully radiative pure hydrogen atmospheres that are the least complicated to model. Pulsationally stable, they are fully characterized by their effective temperature Teff, and surface gravity log g, which can be deduced from their optical spectra and used in model atmospheres to predict their spectral energy distribution (SED). Based on this, three bright DAWDs have defined the spectrophotometric flux scale of the CALSPEC system of HST. In this paper we add 32 new fainter (16.5 < V < 19.5) DAWDs spread over the whole sky and within the dynamic range of large telescopes. Using ground based spectra and panchromatic photometry with HST/WFC3, a new hierarchical analysis process demonstrates consistency between model and observed fluxes above the terrestrial atmosphere to < 0.004 mag rms from 2700 Å to 7750 Å and to 0.008 mag rms at 1.6μm for the total set of 35 DAWDs. These DAWDs are thus established as spectrophotometric standards with unprecedented accuracy from the near ultraviolet to the near-infrared, suitable for both ground and space based observatories. They are embedded in existing surveys like SDSS, PanSTARRS and GAIA, and will be naturally included in the LSST survey by Rubin Observatory. With additional data and analysis to extend the validity of their SEDs further into the IR, these spectrophotometric standard stars could be used for JWST, as well as for the Roman and Euclid observatories.

We verified for photometric stability a set of DA white dwarfs with Hubble Space Telescope magnitudes from the near-ultraviolet to the near-infrared and ground-based spectroscopy by using time-spaced observations from the Las Cumbres Observatory network of telescopes. The initial list of 38 stars was whittled to 32 final ones which comprise a high quality set of spectrophotometric standards. These stars are homogeneously distributed around the sky and are all fainter than r 16.5 mag. Their distribution is such that at least two of them would be available to be observed from any observatory on the ground at any time at airmass less than two. Light curves and different variability indices from the Las Cumbres Observatory data were used to determine the stability of the candidate standards. When available, Pan-STARRS1, Zwicky Transient Facility and TESS data were also used to confirm the star classification. Our analysis showed that four DA white dwarfs may exhibit evidence of photometric variability, while a fifth is cooler than our established lower temperature limit, and a sixth star might be a binary. In some instances, due to the presence of faint nearby red sources, care should be used when observing a few of the spectrophotometric standards with ground-based telescopes. Light curves and finding charts for all the stars are provided.

Ongoing large-scale optical time-domain surveys, such as the Zwicky Transient Facility (ZTF), are producing alerts at unprecedented rates. Analysis of transient sources has so far followed two distinct paths: archival analysis of data on transient sources at a time when they are no longer observable and real-time analysis at the time when the sources are first detected. The latter is the realm of alert brokers such as the Arizona-NOIRLab Temporal Analysis and Response to Events System (ANTARES). In this paper, we synthesize the two analysis paths and present a first systematic study of archival alert-broker data, focusing on extragalactic transients with multipeaked light curves identified in the ANTARES archive of ZTF alerts. Our analysis yields a sample of 37 such sources, including core-collapse supernovae (with two analogs of iPTF14hls), thermonuclear supernovae interacting with their surrounding circumstellar medium, tidal disruption events, luminous blue variables, and as yet unclassified objects. A large fraction of the identified sources is currently active, warranting allocation of follow-up resources in the immediate future to further constrain their nature and the physical processes at work.

We present three epochs of early-time ultraviolet (UV) and optical HST/STIS spectroscopy of the young, nearby Type IIP supernova (SN) 2021yja. We complement the HST data with two earlier epochs of Swift UVOT spectroscopy. The HST and Swift UVOT spectra are consistent with those of other well-studied Type IIP SNe. The UV spectra exhibit rapid cooling at early times, while less dramatic changes are seen in the optical. We also present Lick/KAIT optical photometry up to the late-time tail phase, showing a very long plateau and shallow decline compared with other SNe IIP. Our modeling of the UV spectrum with the TARDIS radiative transfer code produces a good fit for a high-velocity explosion, a low total extinction E(B - V) = 0.07 mag, and a subsolar metallicity. We do not find a significant contribution to the UV flux from an additional heating source, such as interaction with the circumstellar medium, consistent with the observed flat plateau. Furthermore, the velocity width of the Mg ii lambda 2798 line is comparable to that of the hydrogen Balmer lines, suggesting that the UV emission is confined to a region close to the photosphere.