Search

We present multi-wavelength data of SN2020acct, a double-peaked stripped-envelope supernova (SN) in NGC2981 at ~150 Mpc. The two peaks are temporally distinct, with maxima separated by 58 rest-frame days, and a factor of 20 reduction in flux between. The first is luminous (M$_{r}$ = -18.00 $\pm$ 0.02 mag), blue (g - r = 0.27 $\pm$ 0.03 mag), and displays spectroscopic signatures of interaction with hydrogen-free circumstellar material. The second peak is fainter (M$_{r}$ = -17.29 $\pm$ 0.03 mag), and spectroscopically similar to an evolved stripped-envelope SNe, with strong blended forbidden [Ca II] and [O II] features. No other known double-peak SN exhibits a light curve similar to that of SN 2020acct. We find the likelihood of two individual SNe occurring in the same star-forming region within that time to be highly improbable, while an implausibly fine-tuned configuration would be required to produce two SNe from a single binary system. We find that the peculiar properties of SN2020acct match models of pulsational pair instability (PPI), in which the initial peak is produced by collisions of shells of ejected material, shortly followed by a terminal explosion. Pulsations from a star with a 72 M$_{\odot}$ helium core provide an excellent match to the double-peaked light curve. The local galactic environment has a metallicity of 0.4 Z$_{\odot}$, a level where massive single stars are not expected retain enough mass to encounter the PPI. However, late binary mergers or a low-metallicity pocket may allow the required core mass. We measure the rate of SN 2020acct-like events to be $<3.3\times10^{-8}$ Mpc$^{-3}$ yr$^{-1}$ at z = 0.07, or <0.1% of the total core-collapse SN rate., Comment: 31 pages, 14 figures. Submitted to ApJL, comments welcome

We present UV/optical/NIR observations and modeling of supernova (SN) 2024ggi, a type II supernova (SN II) located in NGC 3621 at 7.2 Mpc. Early-time ("flash") spectroscopy of SN 2024ggi within +0.8 days of discovery shows emission lines of H I, He I, C III, and N III with a narrow core and broad, symmetric wings (i.e., IIn-like) arising from the photoionized, optically-thick, unshocked circumstellar material (CSM) that surrounded the progenitor star at shock breakout. By the next spectral epoch at +1.5 days, SN 2024ggi showed a rise in ionization as emission lines of He II, C IV, N IV/V and O V became visible. This phenomenon is temporally consistent with a blueward shift in the UV/optical colors, both likely the result of shock breakout in an extended, dense CSM. The IIn-like features in SN 2024ggi persist on a timescale of $t_{\rm IIn} = 3.8 \pm 1.6$ days at which time a reduction in CSM density allows the detection of Doppler broadened features from the fastest SN material. SN 2024ggi has peak UV/optical absolute magnitudes of $M_{\rm w2} = -18.7$ mag and $M_{\rm g} = -18.1$ mag that are consistent with the known population of CSM-interacting SNe II. Comparison of SN 2024ggi with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium (nLTE) radiative-transfer simulations suggests a progenitor mass-loss rate of $\dot{M} = 10^{-2}$M$_{\odot}$ yr$^{-1}$ ($v_w$ = 50 km/s), confined to a distance of $r < 5\times 10^{14}$ cm. Assuming a wind velocity of $v_w$ = 50 km/s, the progenitor star underwent an enhanced mass-loss episode in the last ~3 years before explosion., Comment: Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:2306.04721, arXiv:2403.02382

We present an ultraviolet-to-infrared search for the electromagnetic (EM) counterpart to GW190425, the second-ever binary neutron star (BNS) merger discovered by the LIGO-Virgo-KAGRA Collaboration (LVK). GW190425 was more distant and had a larger localization area than GW170817, therefore we use a new tool teglon to redistribute the GW190425 localization probability in the context of galaxy catalogs within the final localization volume. We derive a 90th percentile area of 6,688 deg$^{2}$, a $\sim$1.5$\times$ improvement relative to the LIGO/Virgo map, and show how teglon provides an order of magnitude boost to the search efficiency of small ($\leq$1 deg$^{2}$) field-of-view instruments. We combine our data with all publicly reported imaging data, covering 9,078.59 deg$^2$ of unique area and 48.13% of the LIGO/Virgo-assigned localization probability, to calculate the most comprehensive kilonova, short gamma-ray burst (sGRB) afterglow, and model-independent constraints on the EM emission from a hypothetical counterpart to GW190425 to date under the assumption that no counterpart was found in these data. If the counterpart were similar to AT 2017gfo, there was a 28.4% chance that it would have been detected in the combined dataset. We are relatively insensitive to an on-axis sGRB, and rule out a generic transient with a similar peak luminosity and decline rate as AT 2017gfo to 30% confidence. Finally, across our new imaging and all publicly-reported data, we find 28 candidate optical counterparts that we cannot rule out as being associated with GW190425, finding that 4 such counterparts discovered within the localization volume and within 5 days of merger exhibit luminosities consistent with a kilonova., Comment: 41 pages, 11 figures, Submitted to ApJ

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., Comment: 44 pages (68 pages with Appendix), 15 figures, accepted to ApJ

We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early-time ($\delta t < 2$ days) spectra show transient, narrow emission lines from shock ionization of confined ($r < 10^{15}$ cm) circumstellar material (CSM). The observed electron-scattering broadened line profiles (i.e., IIn-like) of HI, He I/II, C III/IV, and N III/IV/V from the CSM persist on a characteristic timescale ($t_{\rm IIn}$) that marks a transition to a lower-density CSM and the emergence of Doppler-broadened features from the fast-moving SN ejecta. Our sample, the largest to date, consists of 39 SNe with early-time IIn-like features in addition to 35 "comparison" SNe with no evidence of early-time IIn-like features, all with ultraviolet observations. The total sample consists of 50 unpublished objects with 474 previously unpublished spectra and 50 multiband light curves, collected primarily through the Young Supernova Experiment and Global Supernova Project collaborations. For all sample objects, we find a significant correlation between peak ultraviolet brightness and both $t_{\rm IIn}$ and the rise time, as well as evidence for enhanced peak luminosities in SNe II with IIn-like features. We quantify mass-loss rates and CSM density for the sample through matching of peak multiband absolute magnitudes, rise times, $t_{\rm IIn}$ and optical SN spectra with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium (nLTE) radiative-transfer simulations. For our grid of models, all with the same underlying explosion, there is a trend between the duration of the electron-scattering broadened line profiles and inferred mass-loss rate: $t_{\rm IIn} \approx 3.8[\dot{M}/(0.01 \textrm{M}_{\odot} \textrm{yr}^{-1})]$ days., Comment: 58 pages, 24 figures, submitted to ApJ. Supplementary figures available on Github (https://github.com/wynnjacobson-galan/Flash_Spectra_Sample). Data release following publication

We present the full Hubble diagram of photometrically-classified Type Ia supernovae (SNe Ia) from the Dark Energy Survey supernova program (DES-SN). DES-SN discovered more than 20,000 SN candidates and obtained spectroscopic redshifts of 7,000 host galaxies. Based on the light-curve quality, we select 1635 photometrically-identified SNe Ia with spectroscopic redshift 0.10$< z <$1.13, which is the largest sample of supernovae from any single survey and increases the number of known $z>0.5$ supernovae by a factor of five. In a companion paper, we present cosmological results of the DES-SN sample combined with 194 spectroscopically-classified SNe Ia at low redshift as an anchor for cosmological fits. Here we present extensive modeling of this combined sample and validate the entire analysis pipeline used to derive distances. We show that the statistical and systematic uncertainties on cosmological parameters are $\sigma_{\Omega_M,{\rm stat+sys}}^{\Lambda{\rm CDM}}=$0.017 in a flat $\Lambda$CDM model, and $(\sigma_{\Omega_M},\sigma_w)_{\rm stat+sys}^{w{\rm CDM}}=$(0.082, 0.152) in a flat $w$CDM model. Combining the DES SN data with the highly complementary CMB measurements by Planck Collaboration (2020) reduces uncertainties on cosmological parameters by a factor of 4. In all cases, statistical uncertainties dominate over systematics. We show that uncertainties due to photometric classification make up less than 10% of the total systematic uncertainty budget. This result sets the stage for the next generation of SN cosmology surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time., Comment: 39 pages, 19 figures; Submitted to ApJ; companion paper Dark Energy Collaboration et al. on consecutive arxiv number 2401.02929

We present cosmological constraints from the sample of Type Ia supernovae (SN Ia) discovered during the full five years of the Dark Energy Survey (DES) Supernova Program. In contrast to most previous cosmological samples, in which SN are classified based on their spectra, we classify the DES SNe using a machine learning algorithm applied to their light curves in four photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host galaxies. After accounting for the likelihood of each SN being a SN Ia, we find 1635 DES SNe in the redshift range $0.100.5$ SNe compared to the previous leading compilation of Pantheon+, and results in the tightest cosmological constraints achieved by any SN data set to date. To derive cosmological constraints we combine the DES supernova data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning $0.025

We present the Keck Infrared Transient Survey (KITS), a NASA Key Strategic Mission Support program to obtain near-infrared (NIR) spectra of astrophysical transients of all types, and its first data release, consisting of 105 NIR spectra of 50 transients. Such a data set is essential as we enter a new era of IR astronomy with the James Webb Space Telescope (JWST) and the upcoming Nancy Grace Roman Space Telescope (Roman). NIR spectral templates will be essential to search JWST images for stellar explosions of the first stars and to plan an effective Roma} SN Ia cosmology survey, both key science objectives for mission success. Between 2022 February and 2023 July, we systematically obtained 274 NIR spectra of 146 astronomical transients, representing a significant increase in the number of available NIR spectra in the literature. The first data release includes data from the 2022A semester. We systematically observed three samples: a flux-limited sample that includes all transients $<$17 mag in a red optical band (usually ZTF r or ATLAS o bands); a volume-limited sample including all transients within redshift $z < 0.01$ ($D \approx 50$ Mpc); and an SN Ia sample targeting objects at phases and light-curve parameters that had scant existing NIR data in the literature. The flux-limited sample is 39% complete (60% excluding SNe Ia), while the volume-limited sample is 54% complete and is 79% complete to $z = 0.005$. All completeness numbers will rise with the inclusion of data from other telescopes in future data releases. Transient classes observed include common Type Ia and core-collapse supernovae, tidal disruption events (TDEs), luminous red novae, and the newly categorized hydrogen-free/helium-poor interacting Type Icn supernovae. We describe our observing procedures and data reduction using Pypeit, which requires minimal human interaction to ensure reproducibility.

We present photometric and spectroscopic data for SN 2022joj, a nearby peculiar Type Ia supernova (SN Ia) with a fast decline rate ($\rm{\Delta m_{15,B}=1.4}$ mag). SN 2022joj shows exceedingly red colors, with a value of approximately ${B-V \approx 1.1}$ mag during its initial stages, beginning from $11$ days before maximum brightness. As it evolves the flux shifts towards the blue end of the spectrum, approaching ${B-V \approx 0}$ mag around maximum light. Furthermore, at maximum light and beyond, the photometry is consistent with that of typical SNe Ia. This unusual behavior extends to its spectral characteristics, which initially displayed a red spectrum and later evolved to exhibit greater consistency with typical SNe Ia. We consider two potential explanations for this behavior: double detonation from a helium shell on a sub-Chandrasekhar-mass white dwarf and Chandrasekhar-mass models with a shallow distribution of $\rm{^{56}Ni}$. The shallow nickel models could not reproduce the red colors in the early light curves. Spectroscopically, we find strong agreement between SN 2022joj and double-detonation models with white dwarf masses around 1 $\rm{M_{\odot}}$ and thin He-shell between 0.01 and 0.02 $\rm{M_{\odot}}$. Moreover, the early red colors are explained by line-blanketing absorption from iron-peak elements created by the double detonation scenario in similar mass ranges. However, the nebular spectra composition in SN 2022joj deviates from expectations for double detonation, as we observe strong [Fe III] emission instead of [Ca II] lines as anticipated from double detonation models. More detailed modeling, e.g., including viewing angle effects, is required to test if double detonation models can explain the nebular spectra.

We present UV/optical observations and models of supernova (SN) 2023ixf, a type II SN located in Messier 101 at 6.9 Mpc. Early-time ("flash") spectroscopy of SN 2023ixf, obtained primarily at Lick Observatory, reveals emission lines of H I, He I/II, C IV, and N III/IV/V with a narrow core and broad, symmetric wings arising from the photo-ionization of dense, close-in circumstellar material (CSM) located around the progenitor star prior to shock breakout. These electron-scattering broadened line profiles persist for $\sim$8 days with respect to first light, at which time Doppler broadened features from the fastest SN ejecta form, suggesting a reduction in CSM density at $r \gtrsim 10^{15}$ cm. The early-time light curve of SN2023ixf shows peak absolute magnitudes (e.g., $M_{u} = -18.6$ mag, $M_{g} = -18.4$ mag) that are $\gtrsim 2$ mag brighter than typical type II supernovae, this photometric boost also being consistent with the shock power supplied from CSM interaction. Comparison of SN 2023ixf to a grid of light curve and multi-epoch spectral models from the non-LTE radiative transfer code CMFGEN and the radiation-hydrodynamics code HERACLES suggests dense, solar-metallicity, CSM confined to $r = (0.5-1) \times 10^{15}$ cm and a progenitor mass-loss rate of $\dot{M} = 10^{-2}$ M$_{\odot}$yr$^{-1}$. For the assumed progenitor wind velocity of $v_w = 50$ km s$^{-1}$, this corresponds to enhanced mass-loss (i.e., ``super-wind'' phase) during the last $\sim$3-6 years before explosion., Comment: 19 pages, 8 figures. Accepted for publication in ApJL

The modern study of astrophysical transients has been transformed by an exponentially growing volume of data. Within the last decade, the transient discovery rate has increased by a factor of ~20, with associated survey data, archival data, and metadata also increasing with the number of discoveries. To manage the data at this increased rate, we require new tools. Here we present YSE-PZ, a transient survey management platform that ingests multiple live streams of transient discovery alerts, identifies the host galaxies of those transients, downloads coincident archival data, and retrieves photometry and spectra from ongoing surveys. YSE-PZ also presents a user with a range of tools to make and support timely and informed transient follow-up decisions. Those subsequent observations enhance transient science and can reveal physics only accessible with rapid follow-up observations. Rather than automating out human interaction, YSE-PZ focuses on accelerating and enhancing human decision making, a role we describe as empowering the human-in-the-loop. Finally, YSE-PZ is built to be flexibly used and deployed; YSE-PZ can support multiple, simultaneous, and independent transient collaborations through group-level data permissions, allowing a user to view the data associated with the union of all groups in which they are a member. YSE-PZ can be used as a local instance installed via Docker or deployed as a service hosted in the cloud. We provide YSE-PZ as an open-source tool for the community., Comment: 23 pages, 9 figures, submitted to PASP

We present extensive optical photometry of the afterglow of GRB~221009A. Our data cover $0.9 - 59.9$\,days from the time of \textit{Swift} and \textit{Fermi} GRB detections. Photometry in $rizy$-band filters was collected primarily with Pan-STARRS and supplemented by multiple 1- to 4-meter imaging facilities. We analyzed the Swift X-ray data of the afterglow and found a single decline rate power-law $f(t) \propto t^{-1.556\pm0.002}$ best describes the light curve. In addition to the high foreground Milky Way dust extinction along this line of sight, the data favour additional extinction to consistently model the optical to X-ray flux with optically thin synchrotron emission. We fit the X-ray-derived power-law to the optical light curve and find good agreement with the measured data up to $5-6$\,days. Thereafter we find a flux excess in the $riy$ bands which peaks in the observer frame at $\sim20$\,days. This excess shares similar light curve profiles to the type Ic broad-lined supernovae SN~2016jca and SN~2017iuk once corrected for the GRB redshift of $z=0.151$ and arbitrarily scaled. This may be representative of a supernova emerging from the declining afterglow. We measure rest-frame absolute peak AB magnitudes of $M_g=-19.8\pm0.6$ and $M_r=-19.4\pm0.3$ and $M_z=-20.1\pm0.3$. If this is an SN component, then Bayesian modelling of the excess flux would imply explosion parameters of $M_{\rm ej}=7.1^{+2.4}_{-1.7}$ M$_{\odot}$, $M_{\rm Ni}=1.0^{+0.6}_{-0.4}$ M$_{\odot}$, and $v_{\rm ej}=33,900^{+5,900}_{-5,700} kms^{-1}$, for the ejecta mass, nickel mass and ejecta velocity respectively, inferring an explosion energy of $E_{\rm kin}\simeq 2.6-9.0\times10^{52}$ ergs., Comment: Submitted to ApJL on 23rd January 2023, for consideration for publication in the special issue on GRB 221009A. Accepted on 2nd March 2023. The results of this paper are under press embargo until 28th March 2023. 15 pages, 7 figures, 1 table

Type Ia supernovae (SNe Ia) are standardizable candles that must be modeled empirically to yield cosmological constraints. To understand the robustness of this modeling to variations in the model training procedure, we build an end-to-end pipeline to test the recently developed SALT3 model. We explore the consequences of removing pre-2000s low-$z$ or poorly calibrated $U$-band data, adjusting the amount and fidelity of SN Ia spectra, and using a model-independent framework to simulate the training data. We find the SALT3 model surfaces are improved by having additional spectra and $U$-band data, and can be shifted by $\sim 5\%$ if host galaxy contamination is not sufficiently removed from SN spectra. We find that resulting measurements of $w$ are consistent to within $2.5\%$ for all training variants explored in this work, with the largest shifts coming from variants that add color-dependent calibration offsets or host galaxy contamination to the training spectra, and those that remove pre-2000s low-$z$ data. These results demonstrate that the SALT3 model training procedure is largely robust to reasonable variations in the training data, but that additional attention must be paid to the treatment of spectroscopic data in the training process. We also find that the training procedure is sensitive to the color distributions of the input data; the resulting $w$ measurement can be biased by $\sim2\%$ if the color distribution is not sufficiently wide. Future low-$z$ data, particularly $u$-band observations and high signal-to-noise ratio SN Ia spectra, will help to significantly improve SN Ia modeling in the coming years., Comment: 17 pages, 10 figures. Published in ApJS

We present the Young Supernova Experiment Data Release 1 (YSE DR1), comprised of processed multi-color Pan-STARRS1 (PS1) griz and Zwicky Transient Facility (ZTF) gr photometry of 1975 transients with host-galaxy associations, redshifts, spectroscopic/photometric classifications, and additional data products from 2019 November 24 to 2021 December 20. YSE DR1 spans discoveries and observations from young and fast-rising supernovae (SNe) to transients that persist for over a year, with a redshift distribution reaching z~0.5. We present relative SN rates from YSE's magnitude- and volume-limited surveys, which are consistent with previously published values within estimated uncertainties for untargeted surveys. We combine YSE and ZTF data, and create multi-survey SN simulations to train the ParSNIP and SuperRAENN photometric classification algorithms; when validating our ParSNIP classifier on 472 spectroscopically classified YSE DR1 SNe, we achieve 82% accuracy across three SN classes (SNe Ia, II, Ib/Ic) and 90% accuracy across two SN classes (SNe Ia, core-collapse SNe). Our classifier performs particularly well on SNe Ia, with high (>90%) individual completeness and purity, which will help build an anchor photometric SNe Ia sample for cosmology. We then use our photometric classifier to characterize our photometric sample of 1483 SNe, labeling 1048 (~71%) SNe Ia, 339 (~23%) SNe II, and 96 (~6%) SNe Ib/Ic. YSE DR1 provides a training ground for building discovery, anomaly detection, and classification algorithms, performing cosmological analyses, understanding the nature of red and rare transients, exploring tidal disruption events and nuclear variability, and preparing for the forthcoming Vera C. Rubin Observatory Legacy Survey of Space and Time., Comment: Accepted to ApJS; 64 pages; 35 figures; 10 tables

There is growing evidence that Type Ia supernovae (SNe Ia) may originate from multiple explosion channels. Previous studies have indicated that the ejecta velocity of SNe Ia is one powerful tool to discriminate between different channels. In this work, we study ~400 confirmed SNe Ia discovered by the Pan-STARRS1 Medium Deep Survey (PS1-MDS), and obtain a sample of ~50 SNe Ia that have near-peak Si II 6355 velocity (Vsi) measurements. We investigate the relationships between Vsi and various parameters, including SN light-curve width, color, host-galaxy properties, and redshift. No significant trends are identified between Vsi and light-curve parameters. Regarding the host-galaxy properties, we see a significant trend that high-velocity (HV) SNe Ia (Vsi > 12000 km/s) tend to reside in more massive galaxies compared to normal-velocity (NV) SNe Ia (Vsi < 12000 km/s) when combining both the PS1-MDS dataset and those from previous low-z studies. While we do not see a significant trend between Vsi and redshift, HV SNe Ia appear to be more prevalent in low-z samples than in high-z samples. We discuss several possibilities that could potentially contribute to this trend. Furthermore, we investigate the potential bias on SN Ia distances and find no significant difference in Hubble residuals between HV and NV subgroups., Comment: 15 pages, 13 figures, 3 tables, accepted for publication in MNRAS

We present optical and near-infrared (NIR) observations of the Type Icn supernova (SN Icn) 2022ann, the fifth member of its newly identified class of SNe. Its early optical spectra are dominated by narrow carbon and oxygen P-Cygni features with absorption velocities of 800 km/s; slower than other SNe Icn and indicative of interaction with a dense, H/He-poor circumstellar medium (CSM) that is outflowing slower than a typical Wolf-Rayet wind velocity of $>$1000 km/s. We identify helium in NIR spectra obtained two weeks after maximum and in optical spectra at three weeks, demonstrating that the CSM is not fully devoid of helium. We never detect broad spectral features from SN ejecta, including in spectra extending to the nebular phase, a unique characteristic among SNe~Icn. Compared to other SNe Icn, SN 2022ann has a low luminosity, with a peak o-band absolute magnitude of -17.7, and evolves slowly. We model the bolometric light curve and find it is well-described by 1.7 M_Sun of SN ejecta interacting with 0.2 M_sun of CSM. We place an upper limit of 0.04 M_Sun of Ni56 synthesized in the explosion. The host galaxy is a dwarf galaxy with a stellar mass of 10^7.34 M_Sun (implied metallicity of log(Z/Z_Sun) $\approx$ 0.10) and integrated star-formation rate of log(SFR) = -2.20 M_sun/yr; both lower than 97\% of the galaxies observed to produce core-collapse supernovae, although consistent with star-forming galaxies on the galaxy Main Sequence. The low CSM velocity, nickel and ejecta masses, and likely low-metallicity environment disfavour a single Wolf-Rayet progenitor star. Instead, a binary companion star is likely required to adequately strip the progenitor before explosion and produce a low-velocity outflow. The low CSM velocity may be indicative of the outer Lagrangian points in the stellar binary progenitor, rather than from the escape velocity of a single Wolf-Rayet-like massive star., Comment: 21 pages, 13 figures, to be submitted to MNRAS

We present pre- and post-explosion observations of the Type II-P supernova (SN~II-P) 2019mhm located in NGC~6753. Based on optical spectroscopy and photometry, we show that SN\,2019mhm exhibits broad lines of hydrogen with a velocity of $-8500\pm200$~km~s$^{-1}$ and a $111\pm2$~day extended plateau in its luminosity, typical of the Type II-P subclass. We also fit its late-time bolometric light curve and infer that it initially produced a ${}^{56}$Ni mass of $1.3 \times 10^{-2}$~M$_{\odot} \pm 5.5 \times 10^{-4}$~$M_\odot$. Using imaging from the Wide Field Planetary Camera 2 on the {\it Hubble Space Telescope} obtained 19~years before explosion, we aligned to a post-explosion Wide Field Camera 3 image and demonstrate that there is no detected counterpart to the SN to a limit of $>$24.53~mag in F814W, corresponding to an absolute magnitude limit of $M_{\rm F814W} < -7.7$~mag. Comparing to massive-star evolutionary tracks, we determine that the progenitor star had a maximum zero-age main sequence mass $<$17.5~M$_{\odot}$, consistent with other SN~II-P progenitor stars. SN\,2019mhm can be added to the growing population of SNe~II-P with both direct constraints on the brightness of their progenitor stars and well-observed SN properties., Comment: 11 pages, 7 figures, updated to match the version ApJ accepted

A large fraction of Type Ia supernova (SN Ia) observations over the next decade will be in the near-infrared (NIR), at wavelengths beyond the reach of the current standard light-curve model for SN Ia cosmology, SALT3 ($\sim 2800$--8700$A$ central filter wavelength). To harness this new SN Ia sample and reduce future light-curve standardization systematic uncertainties, we train SALT3 at NIR wavelengths (SALT3-NIR) up to 2 $\mu$m with the open-source model-training software SALTShaker, which can easily accommodate future observations. Using simulated data we show that the training process constrains the NIR model to $\sim 2$--3% across the phase range ($-20$ to $50$ days). We find that Hubble residual (HR) scatter is smaller using the NIR alone or optical+NIR compared to optical alone, by up to $\sim 30$% depending on filter choice (95% confidence). There is significant correlation between NIR light-curve stretch measurements and luminosity, with stretch and color corrections often improving HR scatter by up to $\sim20%$. For SN Ia observations expected from the \textit{Roman Space Telescope}, SALT3-NIR increases the amount of usable data in the SALT framework by $\sim 20$% at redshift $z\lesssim0.4$ and by $\sim 50$% at $z\lesssim0.15$. The SALT3-NIR model is part of the open-source {\tt SNCosmo} and {\tt SNANA} SN Ia cosmology packages.

The unknown cause of the correlation between Type Ia supernova (SN Ia) Hubble residuals and their host-galaxy masses (the "mass step") may bias cosmological parameter measurements. To better understand the mass step, we develop a SALT3 light-curve model for SN cosmology that uses the host-galaxy masses of 296 low-redshift SNe Ia to derive a spectral-energy distribution--host-galaxy mass relationship. The resulting model has larger Ca II H&K, Ca II near-infrared triplet, and Si II equivalent widths for SNe in low-mass host galaxies at 2.2-2.7$\sigma$ significance; this indicates higher explosion energies per unit mass in low-mass-hosted SNe. The model has phase-dependent changes in SN Ia colors as a function of host mass, indicating intrinsic differences in mean broadband light curves. Although the model provides a better fit to the SN data overall, it does not substantially reduce data--model residuals for a typical light curve in our sample nor does it significantly reduce Hubble residual dispersion. This is because we find that previous SALT models parameterized most host-galaxy dependencies with their first principal component, although they failed to model some significant spectral variations. Our new model is luminosity and cosmology independent, and applying it to data reduces the mass step by $0.021\pm0.002$ mag (uncertainty accounts for correlated data sets); these results indicate that $\sim$35% of the mass step can be attributed to luminosity-independent effects. This SALT model version could be trained using alternative host-galaxy properties and at different redshifts, and therefore will be a tool for understanding redshift-dependent correlations between SNe Ia and their host properties as well as their impact on cosmological parameter measurements., Comment: 19 pages, 14 figures, updated to match the version accepted by ApJ

Massive black holes (BHs) at the centres of massive galaxies are ubiquitous. The population of BHs within dwarf galaxies, on the other hand, is evasive. Dwarf galaxies are thought to harbour BHs with proportionally small masses, including intermediate mass BHs, with masses $10^{2} < M_{BH} < 10^{6} M_{\odot}$. Identification of these systems has historically relied upon the detection of light emitted from accreting gaseous discs close to the BHs. Without this light, they are difficult to detect. Tidal disruption events (TDEs), the luminous flares produced when a star strays close to a BH and is shredded, are a direct way to probe massive BHs. The rise times of these flares theoretically correlate with the BH mass. Here we present AT2020neh, a fast rising TDE candidate, hosted by a dwarf galaxy. AT2020neh can be described by the tidal disruption of a main sequence star by a 10$^{4.7} - 10^{5.9} M_{\odot}$ BH. We find the observable rate of fast rising nuclear transients like AT2020neh to be rare, at $\lesssim 2 \times 10^{-8}$ events Mpc$^{-3}$ yr$^{-1}$. Finding non-accreting BHs in dwarf galaxies is important to determine how prevalent BHs are within these galaxies, and constrain models of BH formation. AT2020neh-like events may provide a galaxy-independent method of measuring IMBH masses., Comment: Accepted for publication in Nature Astronomy

We present photometric and spectroscopic data on three extragalactic luminous red novae (LRNe): AT2018bwo, AT2021afy, and AT2021blu. AT2018bwo was discovered in NGC45 (at 6.8 Mpc) a few weeks after the outburst onset. During the monitoring period, the transient reached a peak luminosity of 10^40 erg/s. AT2021afy, hosted by UGC10043 (49.2 Mpc), showed a double-peaked light curve, with the two peaks reaching a similar luminosity of 2.1(+-0.6)x10^41 erg/s. For AT2021blu in UGC5829, (8.6 Mpc), the pre-outburst phase was well-monitored by several photometric surveys, and the object showed a slow luminosity rise before the outburst. The light curve of AT2021blu was sampled with an unprecedented cadence until the object disappeared behind the Sun, and it was then recovered at late phases. The light curve of AT2021blu shows a double peak, with a prominent early maximum reaching a luminosity of 6.5x10^40 erg/s, which is half of that of AT2021afy. The spectra of AT2021afy and AT2021blu display the expected evolution for LRNe: a blue continuum dominated by prominent Balmer lines in emission during the first peak, and a redder continuum consistent with that of a K-type star with narrow absorption metal lines during the second, broad maximum. The spectra of AT2018bwo are markedly different, with a very red continuum dominated by broad molecular features in absorption. As these spectra closely resemble those of LRNe after the second peak, AT2018bwo was probably discovered at the very late evolutionary stages. This would explain its fast evolution and the spectral properties compatible with that of an M-type star. From the analysis of deep frames of the LRN sites years before the outburst, and considerations of the light curves, the quiescent progenitor systems of the three LRNe were likely massive, with primaries ranging from 13Mo for AT2018bwo, to 13-18Mo for AT2021blu, and over 40Mo for AT2021afy., Comment: 33 pages, 19 figures, 7 tables (plus 3 available at the CDS). Accepted for publication in A&A

Type Ia supernovae (SNe Ia) are more precise standardizable candles when measured in the near-infrared (NIR) than in the optical. With this motivation, from 2012-2017 we embarked on the RAISIN program with the Hubble Space Telescope (HST) to obtain rest-frame NIR light curves for a cosmologically distant sample of 37 SN Ia ($0.2 \lesssim z \lesssim 0.6$) discovered by Pan-STARRS and the Dark Energy Survey. By comparing higher-$z$ HST data with 42 SN Ia at $z<0.1$ observed in the NIR by the Carnegie Supernova Project, we construct a Hubble diagram from NIR observations (with only time of maximum light and some selection cuts from optical data) to pursue a unique avenue to constrain the dark energy equation of state parameter, $w$. We analyze the dependence of the full set of Hubble residuals on the SN Ia host galaxy mass and find Hubble residual steps of size $\sim$0.06-0.1~mag with 1.5- to 2.5-$\sigma$ significance depending on the method and step location. Combining our NIR sample with CMB constraints, we find $1+w=-0.17\pm0.12$ (stat$+$syst). The largest systematic errors are the redshift-dependent SN selection biases and the properties of the NIR mass step. We also use these data to measure $H_0=75.9\pm 2.2$ km s$^{-1}$ Mpc$^{-1}$ from stars with geometric distance calibration in the hosts of 8 SNe Ia observed in the NIR versus $H_0=71.2\pm3.8$ km s$^{-1}$ Mpc$^{-1}$ using an inverse distance ladder approach tied to Planck. Using optical data we find $1+w=-0.10\pm0.09$ and with optical and NIR data combined, we find $1+w=-0.06\pm0.07$; these shifts of up to 0.11 in $w$ could point to inconsistency in optical versus NIR SN models. There will be many opportunities to improve this NIR measurement and better understand systematic uncertainties through larger low-$z$ samples, new light-curve models, calibration improvements, and by building high-$z$ samples from the Roman Space Telescope., Comment: Updated cosmology constraints and text to match published version, and with data release at https://github.com/djones1040/RAISIN_DataRelease

The final explosive fate of massive stars, and the nature of the compact remnants they leave behind (black holes and neutron stars), are major open questions in astrophysics. Many massive stars are stripped of their outer hydrogen envelopes as they evolve. Such Wolf-Rayet (W-R) stars emit strong and rapidly expanding (v_wind>1000 km/s) winds indicating a high escape velocity from the stellar surface. A fraction of this population is also helium depleted, with spectra dominated by highly-ionized emission lines of carbon and oxygen (Types WC/WO). Evidence indicates that the most commonly-observed supernova (SN) explosions that lack hydrogen and helium (Types Ib/Ic) cannot result from massive WC/WO stars, leading some to suggest that most such stars collapse directly into black holes without a visible supernova explosions. Here, we present observations of supernova SN 2019hgp, discovered about a day after explosion. The short rise time and rapid decline place it among an emerging population of rapidly-evolving transients (RETs). Spectroscopy reveals a rich set of emission lines indicating that the explosion occurred within a nebula composed of carbon, oxygen, and neon. Narrow absorption features show that this material is expanding at relatively high velocities (>1500 km/s) requiring a compact progenitor. Our observations are consistent with an explosion of a massive WC/WO star, and suggest that massive W-R stars may be the progenitors of some rapidly evolving transients., Comment: Unedited author version, Nature in press

This note presents an initial survey design for the Nancy Grace Roman High-latitude Time Domain Survey. This is not meant to be a final or exhaustive list of all the survey strategy choices, but instead presents a viable path towards achieving the desired precision and accuracy of dark energy measurements using Type Ia supernovae (SNe Ia). We describe a survey strategy that use six filters (RZYJH and F) and the prism on the Roman Wide Field Instrument. This survey has two tiers, one "wide" which targets SNe Ia at redshifts up to 1 and one "deep" targeting redshifts up to 1.7; for each, four filters are used (with Y and J used in both tiers). We propose one field each in the north and south continuous viewing zones, and expect to obtain high-quality distances of $\sim$12,000 SNe Ia with $\sim$5,000 at z > 1. We propose a wide-tier area of $\sim$19 deg$^2$ and a deep tier of $\sim$5 deg$^2$. Exposure times range from 100 s to 900 s for imaging and 900 s to 3600 s for the prism. These exposure times would reach $\sim$25.5 mag and $\sim$26.5 mag for the wide and deep tiers respectively, with deep co-add stacks reaching $\sim$28 mag and $\sim$29 mag. The total survey spans two years, with a total allocation time of six months, and a cadence of $\sim$5 days., Comment: A report to NASA from the Roman Supernova Science Investigation Teams

SN 2017jgh is a type IIb supernova discovered by Pan-STARRS during the C16/C17 campaigns of the Kepler/K2 mission. Here we present the Kepler/K2 and ground based observations of SN 2017jgh, which captured the shock cooling of the progenitor shock breakout with an unprecedented cadence. This event presents a unique opportunity to investigate the progenitors of stripped envelope supernovae. By fitting analytical models to the SN 2017jgh lightcurve, we find that the progenitor of SN 2017jgh was likely a yellow supergiant with an envelope radius of $\sim50-290~R_{\odot}$, and an envelope mass of $\sim0-1.7~M_{\odot}$. SN 2017jgh likely had a shock velocity of $\sim7500-10300$ km s$^{-1}$. Additionally, we use the lightcurve of SN 2017jgh to investigate how early observations of the rise contribute to constraints on progenitor models. Fitting just the ground based observations, we find an envelope radius of $\sim50-330~R_{\odot}$, an envelope mass of $\sim0.3-1.7~M_{\odot}$ and a shock velocity of $\sim9,000-15,000$ km s$^{-1}$. Without the rise, the explosion time can not be well constrained which leads to a systematic offset in the velocity parameter and larger uncertainties in the mass and radius. Therefore, it is likely that progenitor property estimates through these models may have larger systematic uncertainties than previously calculated., Comment: 15 pages, 9 figures. Accepted to MNRAS

We present multiwavelength observations of the Type II SN 2020pni. Classified at $\sim 1.3$ days after explosion, the object showed narrow (FWHM $<250\,\textrm{km}\,\textrm{s}^{-1}$) recombination lines of ionized helium, nitrogen, and carbon, as typically seen in flash-spectroscopy events. Using the non-LTE radiative transfer code CMFGEN to model our first high resolution spectrum, we infer a progenitor mass-loss rate of $\dot{M}=(3.5-5.3)\times10^{-3}$ M$_{\odot}$ yr$^{-1}$ (assuming a wind velocity of $v_w=200\,\textrm{km}\,\textrm{s}^{-1}$), estimated at a radius of $R_{\rm in}=2.5\times10^{14}\,\rm{cm}$. In addition, we find that the progenitor of SN 2020pni was enriched in helium and nitrogen (relative abundances in mass fractions of 0.30$-$0.40, and $8.2\times10^{-3}$, respectively). Radio upper limits are also consistent with a dense CSM, and a mass-loss rate of $\dot M>5 \times 10^{-4}\,\rm{M_{\odot}\,yr^{-1}}$. During the first 4 days after first light, we also observe an increase in velocity of the hydrogen lines (from $\sim 250\,\textrm{km}\,\textrm{s}^{-1}$ to $\sim 1000\,\textrm{km}\,\textrm{s}^{-1}$), which suggests a complex CSM. The presence of dense and confined CSM, as well as its inhomogeneous structure, suggest a phase of enhanced mass loss of the progenitor of SN 2020pni during the last year before explosion. Finally, we compare SN 2020pni to a sample of other shock-photoionization events. We find no evidence of correlations among the physical parameters of the explosions and the characteristics of the CSM surrounding the progenitors of these events. This favors the idea that the mass-loss experienced by massive stars during their final years could be governed by stochastic phenomena, and that, at the same time, the physical mechanisms responsible for this mass-loss must be common to a variety of different progenitors., Comment: 19 pages, 14 figures

A spectral-energy distribution (SED) model for Type Ia supernovae (SNe Ia) is a critical tool for measuring precise and accurate distances across a large redshift range and constraining cosmological parameters. We present an improved model framework, SALT3, which has several advantages over current models including the leading SALT2 model (SALT2.4). While SALT3 has a similar philosophy, it differs from SALT2 by having improved estimation of uncertainties, better separation of color and light-curve stretch, and a publicly available training code. We present the application of our training method on a cross-calibrated compilation of 1083 SNe with 1207 spectra. Our compilation is $2.5\times$ larger than the SALT2 training sample and has greatly reduced calibration uncertainties. The resulting trained SALT3.K21 model has an extended wavelength range $2000$-$11000$ angstroms (1800 angstroms redder) and reduced uncertainties compared to SALT2, enabling accurate use of low-$z$ $I$ and $iz$ photometric bands. Including these previously discarded bands, SALT3.K21 reduces the Hubble scatter of the low-z Foundation and CfA3 samples by 15% and 10%, respectively. To check for potential systematic uncertainties we compare distances of low ($0.01

We review the needs of the supernova community for improvements in survey coordination and data sharing that would significantly boost the constraints on dark energy using samples of Type Ia supernovae from the Vera C. Rubin Observatories, the \textit{Nancy Grace Roman Space Telescope}, and the \textit{Euclid} Mission. We discuss improvements to both statistical and systematic precision that the combination of observations from these experiments will enable. For example, coordination will result in improved photometric calibration, redshift measurements, as well as supernova distances. We also discuss what teams and plans should be put in place now to start preparing for these combined data sets. Specifically, we request coordinated efforts in field selection and survey operations, photometric calibration, spectroscopic follow-up, pixel-level processing, and computing. These efforts will benefit not only experiments with Type Ia supernovae, but all time-domain studies, and cosmology with multi-messenger astrophysics., Comment: Response to the recent DOE/NASA Request for Information. Endorsed by the Roman Supernova Science Investigation Teams and the LSST DESC Supernova Working Group

With the conclusion of the third observing run for Advanced LIGO/Virgo (O3), we present a detailed analysis of both triggered and serendipitous observations of 17 gravitational wave (GW) events (7 triggered and 10 purely serendipitous) from the Searches After Gravitational-waves Using ARizona Observatories (SAGUARO) program. We searched a total of 4935 deg$^2$ down to a median 5$\sigma$ transient detection depth of 21.1 AB mag using the Mt Lemmon 1.5 m telescope, the discovery engine for SAGUARO. In addition to triggered events within 24~hours, our transient search encompassed a time interval following GW events of $<120$~hrs, providing observations on $\sim$ 1/2 of the events accessible to the Mt Lemmon 1.5 m telescope. We covered 2.1--86\% of the LVC total probability ($P_{\rm total}$) for individual events, with a median $P_{\rm total} \approx 8\%$ within $<120$~hours. Following improvements to our pipeline and the addition of serendipitous observations, we find a total of 7 new optical candidates across 5 GW events which we are unable to rule out after searching for additional information and comparing to kilonova models. Using both publicly available and our own late-time data, we investigated a total of 252 optical candidates for these 17 events, finding only 65\% were followed up in some capacity by the community. Of the total 252 candidates, we are able to rule out an additional 12 previously reported counterpart candidates. In light of these results, we discuss lessons learned from the SAGUARO GW counterpart search. We discuss how community coordination of observations and candidate follow-up, as well as the role of archival data, are crucial to improving the efficiency of follow-up efforts and preventing unnecessary duplication of effort with limited EM resources., Comment: 44 pages, 21 figures, 2 tables, submitted to ApJ

In the next decade, transient searches from the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope will increase the sample of known Type Ia Supernovae (SN Ia) from $\sim10^3$ to $10^5$. With this reduction of statistical uncertainties on cosmological measurements, new methods are needed to reduce systematic uncertainties. Characterizing the underlying spectroscopic evolution of SN Ia remains a major systematic uncertainty in current cosmological analyses, motivating a new simulation tool for the next era of SN Ia cosmology: Build Your Own Spectral Energy Distribution (BYOSED). BYOSED is used within the SNANA framework to simulate light curves by applying spectral variations to model SEDs, enabling flexible testing of possible systematic shifts in SN Ia distance measurements. We test the framework by comparing a nominal Roman SN Ia survey simulation using a baseline SED model to simulations using SEDs perturbed with BYOSED, and investigate the impact of neglecting specific SED features in the analysis. These features include semi-empirical models of two possible, predicted relationships: between SN ejecta velocity and light curve observables, and a redshift-dependent relationship between SN Hubble residuals and host galaxy mass. We analyze each BYOSED simulation using the SALT2 & BBC framework, and estimate changes in the measured value of the dark energy equation-of-state parameter, $w$. We find a difference of $\Delta w=-0.023$ for SN velocity and $\Delta w=0.021$ for redshift-evolving host mass when compared to simulations without these features. By using BYOSED for SN Ia cosmology simulations, future analyses (e.g., Rubin and Roman SN Ia samples) will have greater flexibility to constrain or reduce such SN Ia modeling uncertainties.

The detection of AT2017gfo proved that binary neutron star mergers are progenitors of kilonovae. Using a combination of numerical-relativity and radiative-transfer simulations, the community has developed sophisticated models for these transients for a wide portion of the expected parameter space. Using these simulations and surrogate models made from them, it has been possible to perform Bayesian inference of the observed signals to infer properties of the ejected matter. It has been pointed out that combining inclination constraints derived from the kilonova with gravitational-wave measurements increases the accuracy with which binary parameters can be measured and allows a more accurate inference of the Hubble Constant. In order to not introduce biases, constraints on the inclination angle for AT2017gfo should be insensitive to the employed models. In this work, we compare different assumptions about the ejecta and radiative reprocesses used by the community and we investigate their impact on the parameter inference. While most inferred parameters agree, we find disagreement between posteriors for the inclination angle for different geometries that have been used in the literature. According to our study, the inclusion of reprocessing of the photons between different ejecta types improves the modeling fits to AT2017gfo and in some cases affects the inferred constraints. Our study motivates the inclusion of large $\sim$ 1 mag uncertainties in the kilonova models employed for Bayesian analysis to capture yet unknown systematics, especially when inferring inclination angles, although smaller uncertainties seem appropriate to capture model systematics for other parameters. We also use this method to impose soft constraints on the ejecta geometry of the kilonova AT2017gfo.

Time domain science has undergone a revolution over the past decade, with tens of thousands of new supernovae (SNe) discovered each year. However, several observational domains, including SNe within days or hours of explosion and faint, red transients, are just beginning to be explored. Here, we present the Young Supernova Experiment (YSE), a novel optical time-domain survey on the Pan-STARRS telescopes. Our survey is designed to obtain well-sampled $griz$ light curves for thousands of transient events up to $z \approx 0.2$. This large sample of transients with 4-band light curves will lay the foundation for the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope, providing a critical training set in similar filters and a well-calibrated low-redshift anchor of cosmologically useful SNe Ia to benefit dark energy science. As the name suggests, YSE complements and extends other ongoing time-domain surveys by discovering fast-rising SNe within a few hours to days of explosion. YSE is the only current four-band time-domain survey and is able to discover transients as faint $\sim$21.5 mag in $gri$ and $\sim$20.5 mag in $z$, depths that allow us to probe the earliest epochs of stellar explosions. YSE is currently observing approximately 750 square degrees of sky every three days and we plan to increase the area to 1500 square degrees in the near future. When operating at full capacity, survey simulations show that YSE will find $\sim$5000 new SNe per year and at least two SNe within three days of explosion per month. To date, YSE has discovered or observed 8.3% of the transient candidates reported to the International Astronomical Union in 2020. We present an overview of YSE, including science goals, survey characteristics and a summary of our transient discoveries to date., Comment: ApJ, in press; more information at https://yse.ucsc.edu/

We consider the effects of weak gravitational lensing on observations of 196 spectroscopically confirmed Type Ia Supernovae (SNe Ia) from years 1 to 3 of the Dark Energy Survey (DES). We simultaneously measure both the angular correlation function and the non-Gaussian skewness caused by weak lensing. This approach has the advantage of being insensitive to the intrinsic dispersion of SNe Ia magnitudes. We model the amplitude of both effects as a function of $\sigma_8$, and find $\sigma_8 = 1.2^{+0.9}_{-0.8}$. We also apply our method to a subsample of 488 SNe from the Joint Light-curve Analysis (JLA) (chosen to match the redshift range we use for this work), and find $\sigma_8 = 0.8^{+1.1}_{-0.7}$. The comparable uncertainty in $\sigma_8$ between DES-SN and the larger number of SNe from JLA highlights the benefits of homogeneity of the DES-SN sample, and improvements in the calibration and data analysis., Comment: 10 pages, 7 figures, MNRAS

We present the results from a search for the electromagnetic counterpart of the LIGO/Virgo event S190510g using the Dark Energy Camera (DECam). S190510g is a binary neutron star (BNS) merger candidate of moderate significance detected at a distance of 227$\pm$92 Mpc and localized within an area of 31 (1166) square degrees at 50\% (90\%) confidence. While this event was later classified as likely non-astrophysical in nature within 30 hours of the event, our short latency search and discovery pipeline identified 11 counterpart candidates, all of which appear consistent with supernovae following offline analysis and spectroscopy by other instruments. Later reprocessing of the images enabled the recovery of 6 more candidates. Additionally, we implement our candidate selection procedure on simulated kilonovae and supernovae under DECam observing conditions (e.g., seeing, exposure time) with the intent of quantifying our search efficiency and making informed decisions on observing strategy for future similar events. This is the first BNS counterpart search to employ a comprehensive simulation-based efficiency study. We find that using the current follow-up strategy, there would need to be 19 events similar to S190510g for us to have a 99\% chance of detecting an optical counterpart, assuming a GW170817-like kilonova. We further conclude that optimization of observing plans, which should include preference for deeper images over multiple color information, could result in up to a factor of 1.5 reduction in the total number of followups needed for discovery., Comment: Paper submitted to ApJL

We present a measurement of the Hubble constant $H_0$ using the gravitational wave (GW) event GW190814, which resulted from the coalescence of a 23 $M_\odot$ black hole with a 2.6 $M_\odot$ compact object, as a standard siren. No compelling electromagnetic counterpart has been identified for this event, thus our analysis accounts for thousands of potential host galaxies within a statistical framework. The redshift information is obtained from the photometric redshift (photo-$z$) catalog from the Dark Energy Survey. The luminosity distance is provided by the LIGO/Virgo gravitational wave sky map. Since this GW event has the second-smallest localization volume after GW170817, GW190814 is likely to provide the best constraint on cosmology from a single standard siren without identifying an electromagnetic counterpart. Our analysis uses photo-$z$ probability distribution functions and corrects for photo-$z$ biases. We also reanalyze the binary-black hole GW170814 within this updated framework. We explore how our findings impact the $H_0$ constraints from GW170817, the only GW merger associated with a unique host galaxy. From a combination of GW190814, GW170814 and GW170817, our analysis yields $H_0 = 72.0^{+ 12}_{- 8.2 }~{\rm km~s^{-1}~Mpc^{-1}}$ (68\% Highest Density Interval, HDI) for a prior in $H_0$ uniform between $[20,140]~{\rm km~s^{-1}~Mpc^{-1}}$. The addition of GW190814 and GW170814 to GW170817 improves the 68\% HDI from GW170817 alone by $\sim 18\%$, showing how well-localized mergers without counterparts can provide a significant contribution to standard siren measurements, provided that a complete galaxy catalog is available at the location of the event., Comment: 12 pages, 3 figures, replacement reflects published version

On 14 August 2019, the LIGO and Virgo Collaborations detected gravitational waves from a black hole and a 2.6 solar mass compact object, possibly the first neutron star -- black hole (NSBH) merger. In search of an optical counterpart, the Dark Energy Survey (DES) obtained deep imaging of the entire 90 percent confidence level localization area with Blanco/DECam 0, 1, 2, 3, 6, and 16 nights after the merger. Objects with varying brightness were detected by the DES Pipeline and we systematically reduced the candidate counterparts through catalog matching, light curve properties, host-galaxy photometric redshifts, SOAR spectroscopic follow-up observations, and machine-learning-based photometric classification. All candidates were rejected as counterparts to the merger. To quantify the sensitivity of our search, we applied our selection criteria to full light curve simulations of supernovae and kilonovae as they would appear in the DECam observations. Since the source class of the merger was uncertain, we utilized an agnostic, three-component kilonova model based on tidally-disrupted NS ejecta properties to quantify our detection efficiency of a counterpart if the merger included a NS. We find that if a kilonova occurred during this merger, configurations where the ejected matter is greater than 0.07 solar masses, has lanthanide abundance less than $10^{-8.56}$, and has a velocity between $0.18c$ and $0.21c$ are disfavored at the $2\sigma$ level. Furthermore, we estimate that our background reduction methods are capable of associating gravitational wave signals with a detected electromagnetic counterpart at the $4\sigma$ level in $95\%$ of future follow-up observations., Comment: Published in ApJ

We present observations of ASASSN-19dj, a nearby tidal disruption event (TDE) discovered in the post-starburst galaxy KUG 0810+227 by the All-Sky Automated Survey for Supernovae (ASAS-SN) at a distance of d $\simeq98$ Mpc. We observed ASASSN-19dj from $-$21 to 392 d relative to peak ultraviolet (UV)/optical emission using high-cadence, multiwavelength spectroscopy and photometry. From the ASAS-SN $g$-band data, we determine that the TDE began to brighten on 2019 February 6.8 and for the first 16 d the rise was consistent with a flux $\propto t^2$ power-law. ASASSN-19dj peaked in the UV/optical on 2019 March 6.5 (MJD = 58548.5) at a bolometric luminosity of $L = (6.2 \pm 0.2) \times 10^{44} \text{ erg s}^{-1}$. Initially remaining roughly constant in X-rays and slowly fading in the UV/optical, the X-ray flux increased by over an order of magnitude $\sim$225 d after peak, resulting from the expansion of the X-ray emitting region. The late-time X-ray emission is well fitted by a blackbody with an effective radius of $\sim1 \times 10^{12} \text{ cm}$ and a temperature of $\sim6 \times 10^{5} \text{ K}$. The X-ray hardness ratio becomes softer after brightening and then returns to a harder state as the X-rays fade. Analysis of Catalina Real-Time Transient Survey images reveals a nuclear outburst roughly 14.5 yr earlier with a smooth decline and a luminosity of $L_V\geq1.4 \times 10^{43}$ erg s$^{-1}$, although the nature of the flare is unknown. ASASSN-19dj occurred in the most extreme post-starburst galaxy yet to host a TDE, with Lick H$\delta_{A}$ = $7.67\pm0.17$ \AA., Comment: Updated to match accepted MNRAS version

We present a description of the Australian Dark Energy Survey (OzDES) and summarise the results from its six years of operations. Using the 2dF fibre positioner and AAOmega spectrograph on the 3.9-metre Anglo-Australian Telescope, OzDES has monitored 771 AGN, classified hundreds of supernovae, and obtained redshifts for thousands of galaxies that hosted a transient within the 10 deep fields of the Dark Energy Survey. We also present the second OzDES data release, containing the redshifts of almost 30,000 sources, some as faint as $r_{\mathrm AB}=24$ mag, and 375,000 individual spectra. These data, in combination with the time-series photometry from the Dark Energy Survey, will be used to measure the expansion history of the Universe out to $z\sim1.2$ and the masses of hundreds of black holes out to $z\sim4$. OzDES is a template for future surveys that combine simultaneous monitoring of targets with wide-field imaging cameras and wide-field multi-object spectrographs., Comment: 20 pages, 16 figures. Accepted for publication by MNRAS. Data release available at https://datacentral.org.au

Despite vast improvements in the measurement of the cosmological parameters, the nature of dark energy and an accurate value of the Hubble constant (H$_0$) in the Hubble-Lema\^itre law remain unknown. To break the current impasse, it is necessary to develop as many independent techniques as possible, such as the use of Type II supernovae (SNe II). The goal of this paper is to demonstrate the utility of SNe II for deriving accurate extragalactic distances, which will be an asset for the next generation of telescopes where more-distant SNe II will be discovered. More specifically, we present a sample from the Dark Energy Survey Supernova Program (DES-SN) consisting of 15 SNe II with photometric and spectroscopic information spanning a redshift range up to 0.35. Combining our DES SNe with publicly available samples, and using the standard candle method (SCM), we construct the largest available Hubble diagram with SNe II in the Hubble flow (70 SNe II) and find an observed dispersion of 0.27 mag. We demonstrate that adding a colour term to the SN II standardisation does not reduce the scatter in the Hubble diagram. Although SNe II are viable as distance indicators, this work points out important issues for improving their utility as independent extragalactic beacons: find new correlations, define a more standard subclass of SNe II, construct new SN II templates, and dedicate more observing time to high-redshift SNe II. Finally, for the first time, we perform simulations to estimate the redshift-dependent distance-modulus bias due to selection effects., Comment: 39 pages, 22 figures, 10 tables, Accepted for publication in MNRAS

Rapidly evolving transients (RETs), also termed fast blue optical transients (FBOTs), are a distinct class of astrophysical event. They are characterised by lightcurves that decline much faster than common classes supernovae (SNe), span vast ranges in peak luminosity and can be seen to redshifts greater than 1. Their evolution on fast timescales has hindered high quality follow-up observations, and thus their origin and explosion/emission mechanism remains unexplained. In this paper we define the largest sample of RETs to date, comprising 106 objects from the Dark Energy Survey, and perform the most comprehensive analysis of RET host galaxies. Using deep-stacked photometry and emission-lines from OzDES spectroscopy, we derive stellar masses and star-formation rates (SFRs) for 49 host galaxies, and metallicities for 37. We find that RETs explode exclusively in star-forming galaxies and are thus likely associated with massive stars. Comparing RET hosts to samples of host galaxies of other explosive transients as well as field galaxies, we find that RETs prefer galaxies with high specific SFRs, indicating a link to young stellar populations, similar to stripped-envelope SNe. RET hosts appear to show a lack of chemical enrichment, their metallicities akin to long duration gamma-ray bursts and superluminous SN host galaxies. There are no clear relationships between properties of the host galaxies and the peak magnitudes or decline rates of the transients themselves., Comment: 16 pages, 15 figures, submitted to MNRAS

After correcting for their light-curve shape and color, Type Ia supernovae (SNe Ia) are precise cosmological distance indicators. However, there remains a non-zero intrinsic scatter in the differences between measured distance and that inferred from a cosmological model (i.e., Hubble residuals or HRs), indicating that SN Ia distances can potentially be further improved. We use the open-source relational database kaepora to generate composite spectra with desired average properties of phase, light-curve shape, and HR. At many phases, the composite spectra from two subsamples with positive and negative average HRs are significantly different. In particular, in all spectra from 9 days before to 15 days after peak brightness, we find that SNe with negative HRs have, on average, higher ejecta velocities (as seen in nearly every optical spectral feature) than SNe with positive HRs. At +4 days relative to B-band maximum, using a sample of 62 SNe Ia, we measure a 0.091 +/- 0.035 mag HR step between SNe with Si II 6355 line velocities higher/lower than -11,000 km/s (the median velocity). After light-curve shape and color correction, SNe with higher velocities tend to have underestimated distance moduli relative to a cosmological model. The intrinsic scatter in our sample reduces from 0.094 mag to 0.082 mag after making this correction. Using the Si II 6355 velocity evolution of 115 SNe Ia, we estimate that a velocity difference > 500 km/s exists at each epoch between the positive-HR and negative-HR samples with 99.4% confidence. Finally at epochs later than +37 days, we observe that negative-HR composite spectra tend to have weaker spectral features in comparison to positive-HR composite spectra., Comment: 14 pages, 15 figures. Accepted for publication in MNRAS. Repository at https://github.com/msiebert1/kaepora , database and composite spectra at https://msiebert1.github.io/kaepora/ , and our documentation at https://kaepora.readthedocs.io/en/latest/index.html

While many studies have shown a correlation between properties of the light curves of Type Ia SN (SNe Ia) and properties of their host galaxies, it remains unclear what is driving these correlations. We introduce a new direct method to study these correlations by analyzing `parent' galaxies that host multiple SNe Ia 'siblings'. Here, we search the Dark Energy Survey SN sample, one of the largest samples of discovered SNe, and find 8 galaxies that hosted two likely Type Ia SNe. Comparing the light-curve properties of these SNe and recovered distances from the light curves, we find no better agreement between properties of SNe in the same galaxy as any random pair of galaxies, with the exception of the SN light-curve stretch. We show at $2.8\sigma$ significance that at least 1/2 of the intrinsic scatter of SNe Ia distance modulus residuals is not from common host properties. We also discuss the robustness with which we could make this evaluation with LSST, which will find $100\times$ more pairs of galaxies, and pave a new line of study on the consistency of Type Ia supernovae in the same parent galaxies. Finally, we argue that it is unlikely some of these SNe are actually single, lensed SN with multiple images., Comment: Submitted to ApJL. Comments welcome. Fermilab id included

We present DES16C3cje, a low-luminosity, long-lived type II supernova (SN II) at redshift 0.0618, detected by the Dark Energy Survey (DES). DES16C3cje is a unique SN. The spectra are characterized by extremely narrow photospheric lines corresponding to very low expansion velocities of $\lesssim1500$ km s$^{-1}$, and the light curve shows an initial peak that fades after 50 days before slowly rebrightening over a further 100 days to reach an absolute brightness of M$_r\sim -15.5$ mag. The decline rate of the late-time light curve is then slower than that expected from the powering by radioactive decay of $^{56}$Co but is comparable to that expected from accretion power. Comparing the bolometric light curve with hydrodynamical models, we find that DES16C3cje can be explained by either i) a low explosion energy (0.11 foe) and relatively large $^{56}$Ni production of 0.075 M$_{\odot}$ from a $\sim15$ M$_{\odot}$ red supergiant progenitor typical of other SNe II, or ii) a relatively compact $\sim40$ M$_{\odot}$ star, explosion energy of 1 foe, and 0.08 M$_{\odot}$ of $^{56}$Ni. Both scenarios require additional energy input to explain the late-time light curve, which is consistent with fallback accretion at a rate of $\sim0.5\times{10^{-8}}$ M$_{\odot}$ s$^{-1}$., Comment: Accepted in MNRAS. 17 pages, 6 figures

We present improved photometric measurements for the host galaxies of 206 spectroscopically confirmed type Ia supernovae discovered by the Dark Energy Survey Supernova Program (DES-SN) and used in the first DES-SN cosmological analysis. Fitting spectral energy distributions to the $griz$ photometric measurements of the DES-SN host galaxies, we derive stellar masses and star-formation rates. For the DES-SN sample, when considering a 5D ($z$, $x_1$, $c$, $\alpha$, $\beta$) bias correction, we find evidence of a Hubble residual `mass step', where SNe Ia in high mass galaxies ($>10^{10} \textrm{M}_{\odot}$) are intrinsically more luminous (after correction) than their low mass counterparts by $\gamma=0.040\pm0.019$mag. This value is larger by $0.031$mag than the value found in the first DES-SN cosmological analysis. This difference is due to a combination of updated photometric measurements and improved star formation histories and is not from host-galaxy misidentification. When using a 1D (redshift-only) bias correction the inferred mass step is larger, with $\gamma=0.066\pm0.020$mag. The 1D-5D $\gamma$ difference for DES-SN is $0.026\pm0.009$mag. We show that this difference is due to a strong correlation between host galaxy stellar mass and the $x_1$ component of the 5D distance-bias correction. To better understand this effect, we include an intrinsic correlation between light-curve width and stellar mass in simulated SN Ia samples. We show that a 5D fit recovers $\gamma$ with $-9$mmag bias compared to a $+2$mmag bias for a 1D fit. This difference can explain part of the discrepancy seen in the data. Improvements in modeling correlations between galaxy properties and SN is necessary to determine the implications for $\gamma$ and ensure unbiased precision estimates of the dark energy equation-of-state as we enter the era of LSST., Comment: 27 pages, 13 figures; Submitted to MNRAS

Gravitational wave (GW) events detectable by LIGO and Virgo have several possible progenitors, including black hole mergers, neutron star mergers, black hole--neutron star mergers, supernovae, and cosmic string cusps. A subset of GW events are expected to produce electromagnetic (EM) emission that, once detected, will provide complementary information about their astrophysical context. To that end, the LIGO--Virgo Collaboration (LVC) sends GW candidate alerts to the astronomical community so that searches for their EM counterparts can be pursued. The DESGW group, consisting of members of the Dark Energy Survey (DES), the LVC, and other members of the astronomical community, uses the Dark Energy Camera (DECam) to perform a search and discovery program for optical signatures of LVC GW events. DESGW aims to use a sample of GW events as standard sirens for cosmology. Due to the short decay timescale of the expected EM counterparts and the need to quickly eliminate survey areas with no counterpart candidates, it is critical to complete the initial analysis of each night's images as quickly as possible. We discuss our search area determination, imaging pipeline, and candidate selection processes. We review results from the DESGW program during the first two LIGO--Virgo observing campaigns and introduce other science applications that our pipeline enables., Comment: 15 pages, 7 figures, accepted in Astronomy and Computing, matches accepted version

The five-year Dark Energy Survey supernova programme (DES-SN) is one of the largest and deepest transient surveys to date in terms of volume and number of supernovae. Identifying and characterising the host galaxies of transients plays a key role in their classification, the study of their formation mechanisms, and the cosmological analyses. To derive accurate host galaxy properties, we create depth-optimised coadds using single-epoch DES-SN images that are selected based on sky and atmospheric conditions. For each of the five DES-SN seasons, a separate coadd is made from the other 4 seasons such that each SN has a corresponding deep coadd with no contaminating SN emission. The coadds reach limiting magnitudes of order $\sim 27$ in $g$-band, and have a much smaller magnitude uncertainty than the previous DES-SN host templates, particularly for faint objects. We present the resulting multi-band photometry of host galaxies for samples of spectroscopically confirmed type Ia (SNe Ia), core-collapse (CCSNe), and superluminous (SLSNe) as well as rapidly evolving transients (RETs) discovered by DES-SN. We derive host galaxy stellar masses and probabilistically compare stellar-mass distributions to samples from other surveys. We find that the DES spectroscopically confirmed sample of SNe Ia selects preferentially fewer high mass hosts at high redshift compared to other surveys, while at low redshift the distributions are consistent. DES CCSNe and SLSNe hosts are similar to other samples, while RET hosts are unlike the hosts of any other transients, although these differences have not been disentangled from selection effects., Comment: 18 pages, 17 figures. Submitted to MNRAS

We present an analysis of DES17X1boj and DES16E2bjy, two peculiar transients discovered by the Dark Energy Survey (DES). They exhibit nearly identical double-peaked light curves which reach very different maximum luminosities (M$_\mathrm{r}$ = -15.4 and M$_\mathrm{r}$ = -17.9, respectively). The light curve evolution of these events is highly atypical and has not been reported before. The transients are found in different host environments: DES17X1boj was found near the nucleus of a spiral galaxy, while DES16E2bjy is located in the outskirts of a passive red galaxy. Early photometric data is well fitted with a blackbody and the resulting moderate photospheric expansion velocities (1800 km/s for DES17X1boj and 4800 km/s for DES16E2bjy) suggest an explosive or eruptive origin. Additionally, a feature identified as high-velocity CaII absorption (v $\approx$ 9400km/s) in the near-peak spectrum of DES17X1boj may imply that it is a supernova. While similar light curve evolution suggests a similar physical origin for these two transients, we are not able to identify or characterise the progenitors., Comment: Accepted to MNRAS

We present the discovery of ASASSN-18jd (AT 2018bcb), a luminous optical/UV/X-ray transient located in the nucleus of the galaxy 2MASX J22434289--1659083 at $z=0.1192$. Over the year after discovery, Swift UVOT photometry shows the UV SED of the transient to be well modeled by a slowly shrinking blackbody with temperature $T \sim 2.5 \times 10^{4} \rm ~K$, a maximum observed luminosity of $L_\text{max} = 4.5^{+0.6}_{-0.3} \times 10^{44} \rm ~erg ~s^{-1}$, and a radiated energy of $E = 9.6^{+1.1}_{-0.6} \times 10^{51} \rm ~erg$. X-ray data from Swift XRT and XMM-Newton show a transient, variable X-ray flux with blackbody and power-law components that fade by nearly an order of magnitude over the following year. Optical spectra show strong, roughly constant broad Balmer emission as well as transient features attributable to He II, N III-V, O III, and coronal Fe. While ASASSN-18jd shares similarities with Tidal Disruption Events (TDEs), it is also similar to the newly-discovered nuclear transients seen in quiescent galaxies and faint Active Galactic Nuclei (AGNs)., Comment: 24 pages, 16 figures, 2 machine-readable tables (included as ancillary txt files)

Ultraviolet (UV) observations of Type Ia supernovae (SNe Ia) are crucial for constraining the properties of their progenitor systems. Theoretical studies predicted that the UV spectra, which probe the outermost layers of a SN, should be sensitive to the metal content of the progenitor. Using the largest SN Ia UV (lambda<2900 A) spectroscopic sample obtained from Neil Gehrels Swift Observatory, we investigate the dependence of UV spectra on metallicity. For the first time, our results reveal a correlation (~2 sigma) between SN Ia UV flux and host-galaxy metallicities, with SNe in more metal-rich galaxies (which are likely to have higher progenitor metallicities) having lower UV flux level. We find that this metallicity effect is only significant at short wavelengths (lambda<2700 A), which agrees well with the theoretical predictions. We produce UV spectral templates for SNe Ia at peak brightness. With our sample, we could disentangle the effect of light-curve shape and metallicity on the UV spectra. We also examine the correlation between the UV spectra and SN luminosities as parameterised by Hubble residuals. However, we do not see a significant trend with Hubble residuals. This is probably due to the large uncertainties in SN distances, as the majority of our sample members are extremely nearby (redshift z<0.01). Future work with SNe discovered in the Hubble flow will be necessary to constrain a potential metallicity bias on SN Ia cosmology., Comment: Accepted for publication in MNRAS

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches., Comment: 16 pages, 7 figures, 1 table. Accepted to ApJL