Showing 1–50 of 117 results for author: Drout, M R

The first catalogue of candidate white dwarf-main sequence binaries in open star clusters: A new window into common envelope evolution

Authors: Steffani M. Grondin, Maria R. Drout, Jason Nordhaus, Philip S. Muirhead, Joshua S. Speagle, Ryan Chornock

Abstract: Close binary systems are the progenitors to both Type Ia supernovae and the compact object mergers that can be detected via gravitational waves. To achieve a binary with a small radial separation, it is believed that the system likely undergoes common envelope (CE) evolution. Despite its importance, CE evolution may be one of the largest uncertainties in binary evolution due to a combination of co… ▽ More Close binary systems are the progenitors to both Type Ia supernovae and the compact object mergers that can be detected via gravitational waves. To achieve a binary with a small radial separation, it is believed that the system likely undergoes common envelope (CE) evolution. Despite its importance, CE evolution may be one of the largest uncertainties in binary evolution due to a combination of computational challenges and a lack of observed benchmarks where both the post-CE and pre-CE conditions are known. Identifying post-CE systems in star clusters can partially circumvent this second issue by providing an independent age constraint on the system. For the first time, we conduct a systematic search for white dwarf (WD) and main-sequence (MS) binary systems in 299 Milky Way open star clusters. Coupling Gaia DR3 photometry and kinematics with multi-band photometry from Pan-STARRS1 and 2MASS, we apply a machine learning based approach and find 52 high-probability candidates in 38 open clusters. For a subset of our systems, we present follow-up spectroscopy from the Gemini and Lick Observatories and archival light curves from TESS, Kepler/K2 and the Zwicky Transient Facility. Examples of M-dwarfs with hot companions are spectroscopically observed, along with regular system variability. While the kinematics of our candidates are consistent with their host clusters, some systems have spatial positions offset relative to their hosts, potentially indicative of natal kicks. Ultimately, this catalogue is a first step to obtaining a set of observational benchmarks to better link post-CE systems to their pre-CE progenitors. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 35 pages, 18 figures. Submitted to ApJ

Binary Yellow Supergiants in the Magellanic Clouds I: Photometric Candidate Identification

Authors: Anna J. G. O'Grady, Maria R. Drout, Kathryn F. Neugent, Bethany Ludwig, Ylva Gotberg, Bryan M. Gaensler

Abstract: Recent works have constrained the binary fraction of evolved populations of massive stars in local galaxies such as red supergiants and Wolf-Rayet stars, but the binary fraction of yellow supergiants (YSGs) in the Hertzsprung Gap remains unconstrained. Binary evolution theory predicts that the Hertzsprung Gap is home to multiple populations of binary systems with varied evolutionary histories. In… ▽ More Recent works have constrained the binary fraction of evolved populations of massive stars in local galaxies such as red supergiants and Wolf-Rayet stars, but the binary fraction of yellow supergiants (YSGs) in the Hertzsprung Gap remains unconstrained. Binary evolution theory predicts that the Hertzsprung Gap is home to multiple populations of binary systems with varied evolutionary histories. In this paper, we develop a method to distinguish single YSGs from YSG plus O- or B-type main sequence binaries using optical and ultraviolet photometry, and then apply this method to identify candidate YSG binaries in the Magellanic Clouds. After constructing a set of combined stellar atmosphere models, we find that optical photometry is, given typical measurement and reddening uncertainties, sufficient to discern single YSGs from YSG+OB binaries if the OB-star is at least $\sim5$M$_{\odot}$ for T$_{\mathrm{eff,YSG}}\sim$ 4000 K, but requires a $\sim$20M$_{\odot}$ OB star for YSGs up to T$_{\mathrm{eff,YSG}}\sim$ 9000 K. For these hotter YSG temperatures, ultraviolet photometry allows binaries with OB companions as small as $\sim$7M$_{\odot}$ to be identified. We use color-color spaces developed from these models to search for evidence of excess blue or ultraviolet light in a set of $\sim$1000 YSG candidates in the Magellanic Clouds. We identify hundreds of candidate YSG binary systems and report a preliminary fraction of YSGs that show a blue/UV color excess of 20-60\%. Spectroscopic follow-up is now required to confirm the true nature of this population. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures, 3 tables. Submitted to The Astrophysical Journal on June 24th 2024

An Untargeted Search for Radio-Emitting Tidal Disruption Events in the VAST Pilot Survey

Authors: Hannah Dykaar, Maria R. Drout, B. M. Gaensler, David L. Kaplan, Tara Murphy, Assaf Horesh, Akash Anumarlapudi, Dougal Dobie, Laura N. Driessen, Emil Lenc, Adam Stewart

Abstract: We present a systematic search for tidal disruption events (TDEs) using radio data from the Variables and Slow Transients (VAST) Pilot Survey conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). Historically, TDEs have been identified using observations at X-ray, optical, and ultraviolet wavelengths. After discovery, a few dozen TDEs have been shown to have radio counterparts… ▽ More We present a systematic search for tidal disruption events (TDEs) using radio data from the Variables and Slow Transients (VAST) Pilot Survey conducted using the Australian Square Kilometre Array Pathfinder (ASKAP). Historically, TDEs have been identified using observations at X-ray, optical, and ultraviolet wavelengths. After discovery, a few dozen TDEs have been shown to have radio counterparts through follow-up observations. With systematic time-domain radio surveys becoming available, we can now identify new TDEs in the radio regime. A population of radio-discovered TDEs has the potential to provide several key insights including an independent constraint on their volumetric rate. We conducted a search to select variable radio sources with a single prominent radio flare and a position consistent within 2$σ$ of the nucleus of a known galaxy. While TDEs were the primary target of our search, sources identified in this search may also be consistent with active galactic nuclei exhibiting unusual flux density changes at the timescales probed, uncharacteristically bright supernovae, or a population of gamma-ray bursts. We identify a sample of 12 radio-bright candidate TDEs. The timescales and luminosities range from ~6 to 230 days and ~10$^{38}$ to 10$^{41}$ erg s$^{-1}$, consistent with models of radio emission from TDEs that launch relativistic jets. After calculating the detection efficiency of our search using a Monte Carlo simulation of TDEs, and assuming all 12 sources are jetted TDEs, we derive a volumetric rate for jetted TDEs of 0.80$^{+0.31}_{-0.23}$ Gpc$^{-3}$ yr$^{-1}$, consistent with previous empirically estimated rates. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 37 pages, 37 figures, accepted to ApJ

SN 2024ggi in NGC 3621: Rising Ionization in a Nearby, CSM-Interacting Type II Supernova

Authors: W. V. Jacobson-Galán, K. W. Davis, C. D. Kilpatrick, L. Dessart, R. Margutti, R. Chornock, R. J. Foley, P. Arunachalam, K. Auchettl, C. R. Bom, R. Cartier, D. A. Coulter, G. Dimitriadis, D. Dickinson, M. R. Drout, A. T. Gagliano, C. Gall, B. Garretson, L. Izzo, D. O. Jones, N. LeBaron, H. -Y. Miao, D. Milisavljevic, Y. -C. Pan, A. Rest , et al. (6 additional authors not shown)

Abstract: 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 circumstel… ▽ More 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. △ Less

Submitted 25 June, 2024; v1 submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:2306.04721, arXiv:2403.02382

The Gravity Collective: A Comprehensive Analysis of the Electromagnetic Search for the Binary Neutron Star Merger GW190425

Authors: D. A. Coulter, C. D. Kilpatrick, D. O. Jones, R. J. Foley, A. V. Filippenko, W. Zheng, J. J. Swift, G. S. Rahman, H. E. Stacey, A. L. Piro, C. Rojas-Bravo, J. Anais Vilchez, N. Muñoz-Elgueta, I. Arcavi, G. Dimitriadis, M. R. Siebert, J. S. Bloom, M. J. Bustamante-Rosell, K. E. Clever, K. W. Davis, J. Kutcka, P. Macias, P. McGill, P. J. Quiñonez, E. Ramirez-Ruiz , et al. (12 additional authors not shown)

Abstract: 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 catalo… ▽ More 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. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 41 pages, 11 figures, Submitted to ApJ

Final Moments II: Observational Properties and Physical Modeling of CSM-Interacting Type II Supernovae

Authors: W. V. Jacobson-Galán, L. Dessart, K. W. Davis, C. D. Kilpatrick, R. Margutti, R. J. Foley, R. Chornock, G. Terreran, D. Hiramatsu, M. Newsome, E. Padilla Gonzalez, C. Pellegrino, D. A. Howell, A. V. Filippenko, J. P. Anderson, C. R. Angus, K. Auchettl, K. A. Bostroem, T. G. Brink, R. Cartier, D. A. Coulter, T. de Boer, M. R. Drout, N. Earl, K. Ertini , et al. (30 additional authors not shown)

Abstract: We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early-time ($δ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… ▽ More We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early-time ($δ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. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 58 pages, 24 figures, submitted to ApJ. Supplementary figures available on Github (https://github.com/wynnjacobson-galan/Flash_Spectra_Sample). Data release following publication

SN2023ixf in Messier 101: the twilight years of the progenitor as seen by Pan-STARRS

Authors: Conor L. Ransome, V. Ashley Villar, Anna Tartaglia, Sebastian Javier Gonzalez, Wynn V. Jacobson-Galán, Charles D. Kilpatrick, Raffaella Margutti, Ryan J. Foley, Matthew Grayling, Yuan Qi Ni, Ricardo Yarza, Christine Ye, Katie Auchettl, Thomas de Boer, Kenneth C. Chambers, David A. Coulter, Maria R. Drout, Diego Farias, Christa Gall, Hua Gao, Mark E. Huber, Adaeze L. Ibik, David O. Jones, Nandita Khetan, Chien-Cheng Lin , et al. (6 additional authors not shown)

Abstract: The nearby type II supernova, SN2023ixf in M101 exhibits signatures of early-time interaction with circumstellar material in the first week post-explosion. This material may be the consequence of prior mass loss suffered by the progenitor which possibly manifested in the form of a detectable pre-supernova outburst. We present an analysis of the long-baseline pre-explosion photometric data in $g$,… ▽ More The nearby type II supernova, SN2023ixf in M101 exhibits signatures of early-time interaction with circumstellar material in the first week post-explosion. This material may be the consequence of prior mass loss suffered by the progenitor which possibly manifested in the form of a detectable pre-supernova outburst. We present an analysis of the long-baseline pre-explosion photometric data in $g$, $w$, $r$, $i$, $z$ and $y$ filters from Pan-STARRS as part of the Young Supernova Experiment, spanning $\sim$5,000 days. We find no significant detections in the Pan-STARRS pre-explosion light curve. We train a multilayer perceptron neural network to classify pre-supernova outbursts. We find no evidence of eruptive pre-supernova activity to a limiting absolute magnitude of $-7$. The limiting magnitudes from the full set of $gwrizy$ (average absolute magnitude $\approx$-8) data are consistent with previous pre-explosion studies. We use deep photometry from the literature to constrain the progenitor of SN2023ixf, finding that these data are consistent with a dusty red supergiant (RSG) progenitor with luminosity $\log\left(L/L_\odot\right)$$\approx$5.12 and temperature $\approx$3950K, corresponding to a mass of 14-20 M$_\odot$ △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 19 pages, 8 figures, 1 table

Spectroscopic r-Process Abundance Retrieval for Kilonovae II: Lanthanides in the Inferred Abundance Patterns of Multi-Component Ejecta from the GW170817 Kilonova

Authors: Nicholas Vieira, John J. Ruan, Daryl Haggard, Nicole M. Ford, Maria R. Drout, Rodrigo Fernández

Abstract: In kilonovae, freshly-synthesized $r$-process elements imprint features on optical spectra, as observed in AT2017gfo, the counterpart to the GW170817 binary neutron star merger. However, measuring the $r$-process compositions of the merger ejecta is computationally challenging. Vieira et al. (2023) introduced Spectroscopic $r$-Process Abundance Retrieval for Kilonovae (SPARK), a software tool to i… ▽ More In kilonovae, freshly-synthesized $r$-process elements imprint features on optical spectra, as observed in AT2017gfo, the counterpart to the GW170817 binary neutron star merger. However, measuring the $r$-process compositions of the merger ejecta is computationally challenging. Vieira et al. (2023) introduced Spectroscopic $r$-Process Abundance Retrieval for Kilonovae (SPARK), a software tool to infer elemental abundance patterns of the ejecta, and associate spectral features with particular species. Previously, we applied SPARK to the 1.4 day spectrum of AT2017gfo and inferred its abundance pattern for the first time, characterized by electron fraction $Y_e=0.31$, a substantial abundance of strontium, and a dearth of lanthanides and heavier elements. This ejecta is consistent with wind from a remnant hypermassive neutron star and/or accretion disk. We now extend our inference to spectra at 2.4 and 3.4 days, and test the need for multicomponent ejecta, where we stratify the ejecta in composition. The ejecta at 1.4 and 2.4 days is described by the same single blue component. At 3.4 days, a new redder component with lower $Y_e=0.16$ and a significant abundance of lanthanides emerges. This new redder component is consistent with dynamical ejecta and/or neutron-rich ejecta from a magnetized accretion disk. As expected from photometric modelling, this component emerges as the ejecta expands, the photosphere recedes, and the earlier bluer component dims. At 3.4 days, we find an ensemble of lanthanides, with the presence of cerium most concrete. This presence of lanthanides has important implications for the contribution of kilonovae to the $r$-process abundances observed in the Universe. △ Less

Submitted 22 December, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: Accepted to ApJ; 19 pages, 8 figures in body; 10 pages, 7 figures in appendix

Type II-P Supernova Progenitor Star Initial Masses and SN 2020jfo: Direct Detection, Light Curve Properties, Nebular Spectroscopy, and Local Environment

Authors: Charles D. Kilpatrick, Luca Izzo, Rory O. Bentley, Kenneth C. Chambers, David A. Coulter, Maria R. Drout, Thomas de Boer, Ryan J. Foley, Christa Gall, Melissa R. Halford, David O. Jones, Danial Langeroodi, Chien-Cheng Lin, Eugene A. Magnier, Peter McGill, Anna J. G. O'Grady, Yen-Chen Pan, Enrico Ramirez-Ruiz, Armin Rest, Jonathan J. Swift, Samaporn Tinyanont, V. Ashley Villar, Richard J. Wainscoat, Amanda Rose Wasserman, S. Karthik Yadavalli , et al. (1 additional authors not shown)

Abstract: We present optical, ultraviolet, and infrared data of the type II supernova (SN II) 2020jfo at 14.5 Mpc. This wealth of multiwavelength data allows to compare different metrics commonly used to estimate progenitor masses of SN II for the same object. Using its early light curve, we infer SN 2020jfo had a progenitor radius of $\approx$700 $R_{\odot}$, consistent with red supergiants of initial mass… ▽ More We present optical, ultraviolet, and infrared data of the type II supernova (SN II) 2020jfo at 14.5 Mpc. This wealth of multiwavelength data allows to compare different metrics commonly used to estimate progenitor masses of SN II for the same object. Using its early light curve, we infer SN 2020jfo had a progenitor radius of $\approx$700 $R_{\odot}$, consistent with red supergiants of initial mass $M_{\rm ZAMS}=$11-13 $M_{\odot}$. The decline in its late-time light curve is best fit by a ${}^{56}$Ni mass of 0.018$\pm$0.007 $M_{\odot}$ consistent with that ejected from SN II-P with $\approx$13 $M_{\odot}$ initial mass stars. Early spectra and photometry do not exhibit signs of interaction with circumstellar matter, implying that SN 2020jfo experienced weak mass loss within the final years prior to explosion. Our spectra at $>$250 days are best fit by models from 12 $M_{\odot}$ initial mass stars. We analyzed integral field unit spectroscopy of the stellar population near SN 2020jfo, finding its massive star population had a zero age main sequence mass of 9.7$\substack{+2.5\\-1.3} M_{\odot}$. We identify a single counterpart in pre-explosion imaging and find it has an initial mass of at most $7.2\substack{+1.2\\-0.6} M_{\odot}$. We conclude that the inconsistency between this mass and indirect mass indicators from SN 2020jfo itself is most likely caused by extinction with $A_{V}=2$-3 mag due to matter around the progenitor star, which lowered its observed optical luminosity. As SN 2020jfo did not exhibit extinction at this level or evidence for interaction with circumstellar matter between 1.6-450 days from explosion, we conclude that this material was likely confined within $\approx$3000 $R_{\odot}$ from the progenitor star. △ Less

Submitted 2 July, 2023; originally announced July 2023.

Comments: 31 pages, 12 figures, accepted to MNRAS

Stellar properties of observed stars stripped in binaries in the Magellanic Clouds

Authors: Y. Gotberg, M. R. Drout, A. P. Ji, J. H. Groh, B. A. Ludwig, P. A. Crowther, N. Smith, A. de Koter, S. E. de Mink

Abstract: Massive stars (~8-25Msun) stripped of their hydrogen-rich envelopes via binary interaction are thought to be the main progenitors for merging neutron stars and stripped-envelope supernovae. We recently presented the discovery of the first set of such stripped stars in a companion paper. Here, we fit the spectra of ten stars with new atmosphere models in order to constrain their stellar properties… ▽ More Massive stars (~8-25Msun) stripped of their hydrogen-rich envelopes via binary interaction are thought to be the main progenitors for merging neutron stars and stripped-envelope supernovae. We recently presented the discovery of the first set of such stripped stars in a companion paper. Here, we fit the spectra of ten stars with new atmosphere models in order to constrain their stellar properties precisely. We find that the stellar properties align well with the theoretical expectations from binary evolution models for helium-core burning envelope-stripped stars. The fits confirm that the stars have high effective temperatures (Teff~50-100kK), high surface gravities (log g ~5), and hydrogen-poor/helium-rich surfaces (X(H, surf)~0-0.4) while showing for the first time a range of bolometric luminosities (10^3-10^5 Lsun), small radii (~0.5-1Rsun), and low Eddington factors (Gamma_e~0.006-0.4). Using these properties, we derive intermediate current masses (~1-8Msun), which suggest that their progenitors were massive stars (~5-25Msun) and that a subset will reach core-collapse, leaving behind neutron stars or black holes. Using the model fits, we also estimate the emission rates of ionizing photons for these stars, which agree well with previous model expectations. Further, by computing models for a range of mass-loss rates, we find that the stellar winds are weaker than predicted by any existing scheme (Mdot(wind)<~ 1e-9 Msun/yr). The properties of this first sample of intermediate mass helium stars suggest they both contain progenitors of type Ib and IIb supernovae, and provide important benchmarks for binary evolution and population synthesis models. △ Less

Submitted 30 June, 2023; originally announced July 2023.

Comments: Accepted for publication in ApJ. Spectral models available via zenodo as described in text. Please also see a companion manuscript detailing the discovery of these systems in today's arXiv posting

Discovery of the missing intermediate-mass helium stars stripped in binaries

Authors: M. R. Drout, Y. Götberg, B. A. Ludwig, J. H. Groh, S. E. de Mink, A. J. G. O'Grady, N. Smith

Abstract: The theory of binary evolution predicts that many massive stars should lose their hydrogen-rich envelopes via interaction with a companion -- revealing hot helium stars with masses of $\sim$2--8M$_{\odot}$. However, only one candidate system had been identified, leaving a large discrepancy between theory and observation. Here, we present a new sample of stars -- identified via excess ultraviolet e… ▽ More The theory of binary evolution predicts that many massive stars should lose their hydrogen-rich envelopes via interaction with a companion -- revealing hot helium stars with masses of $\sim$2--8M$_{\odot}$. However, only one candidate system had been identified, leaving a large discrepancy between theory and observation. Here, we present a new sample of stars -- identified via excess ultraviolet emission -- whose luminosities, colors, and spectral morphologies are consistent with predictions for the missing population. We detect radial velocity variations indicative of binary motion and measure high temperatures ($T_{\rm eff}\sim60-100$kK), high surface gravities ($\log(g)\sim5$) and depleted surface hydrogen mass fractions ($X_{\rm{H,surf}}\lesssim0.3$), which match expectations for stars with initial masses between 8--25 M$_{\odot}$ that have been stripped via binary interaction. These systems fill the helium star mass gap between subdwarfs and Wolf-Rayet stars, and are thought to be of large astrophysical significance as ionizing sources, progenitors of stripped-envelope supernovae and merging double neutron stars. △ Less

Submitted 30 June, 2023; originally announced July 2023.

Comments: Original version of manuscript submitted in August 2022, per journal policies. Manuscript passed through referee process and was recommended for publication in Science in March 2023. Conclusions have not changed. Please contact authors for further information. Data and updated manuscript will be made available upon publication. Please also see a companion paper in today's archive posting

Luminous Radio Emission from the Superluminous Supernova 2017ens at 3.3 years after explosion

Authors: Raffaella Margutti, J. S. Bright, D. J. Matthews, D. L. Coppejans, K. D. Alexander, E. Berger, M. Bietenholz, R. Chornock, L. DeMarchi, M. R. Drout, T. Eftekhari, W. V. Jacobson-Galan, T. Laskar, D. Milisavljevic, K. Murase, M. Nicholl, C. M. B. Omand, M. Stroh, G. Terreran, A. Z. VanderLey

Abstract: We present the results from a multi-year radio campaign of the superluminous supernova (SLSN) 2017ens, which yielded the earliest radio detection of a SLSN to date at the age of $\sim$3.3 years after explosion. SN2017ens was not detected at radio frequencies in the first $\sim$300\,d of evolution but reached $L_ν\approx 10^{28}\,\rm{erg\,s^{-1}\,cm^{-2}}$ at $ν\sim 6$ GHz, $\sim1250$ days post-exp… ▽ More We present the results from a multi-year radio campaign of the superluminous supernova (SLSN) 2017ens, which yielded the earliest radio detection of a SLSN to date at the age of $\sim$3.3 years after explosion. SN2017ens was not detected at radio frequencies in the first $\sim$300\,d of evolution but reached $L_ν\approx 10^{28}\,\rm{erg\,s^{-1}\,cm^{-2}}$ at $ν\sim 6$ GHz, $\sim1250$ days post-explosion. Interpreting the radio observations in the context of synchrotron radiation from the supernova shock interaction with the circumstellar medium (CSM), we infer an effective mass-loss rate of $\approx 10^{-4}\,\rm{M_{\odot}yr^{-1}}$ at $r\sim 10^{17}$ cm from the explosion's site, for a wind speed of $v_w=50-60\,\rm{km\,s^{-1}}$ measured from optical spectra. These findings are consistent with the spectroscopic metamorphosis of SN2017ens from hydrogen-poor to hydrogen-rich $\sim190$ d after explosion reported by Chen et al., 2018. SN2017ens is thus an addition to the sample of hydrogen-poor massive progenitors that explode shortly after having lost their hydrogen envelope. The inferred circumstellar densities, implying a CSM mass up to $\sim0.5\,\rm{M_{\odot}}$, and low velocity of the ejection point at binary interactions (in the form of common envelope evolution and subsequent envelope ejection) playing a role in shaping the evolution of the stellar progenitors of SLSNe in the $\lesssim 500$ yr preceding core collapse. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: 9 pages, 4 figures, 1 table

SN2023ixf in Messier 101: A Variable Red Supergiant as the Progenitor Candidate to a Type II Supernova

Authors: Charles D. Kilpatrick, Ryan J. Foley, Wynn V. Jacobson-Galán, Anthony L. Piro, Stephen J. Smartt, Maria R. Drout, Alexander Gagliano, Christa Gall, Jens Hjorth, David O. Jones, Kaisey S. Mandel, Raffaella Margutti, Conor L. Ransome, V. Ashley Villar, David A. Coulter, Hua Gao, David Jacob Matthews, Yossef Zenati

Abstract: We present pre-explosion optical and infrared (IR) imaging at the site of the type II supernova (SN II) 2023ixf in Messier 101 at 6.9 Mpc. We astrometrically registered a ground-based image of SN 2023ixf to archival Hubble Space Telescope (HST), Spitzer Space Telescope (Spitzer), and ground-based near-IR images. A single point source is detected at a position consistent with the SN at wavelengths… ▽ More We present pre-explosion optical and infrared (IR) imaging at the site of the type II supernova (SN II) 2023ixf in Messier 101 at 6.9 Mpc. We astrometrically registered a ground-based image of SN 2023ixf to archival Hubble Space Telescope (HST), Spitzer Space Telescope (Spitzer), and ground-based near-IR images. A single point source is detected at a position consistent with the SN at wavelengths ranging from HST $R$-band to Spitzer 4.5 $μ$m. Fitting to blackbody and red supergiant (RSG) spectral-energy distributions (SEDs), we find that the source is anomalously cool with a significant mid-IR excess. We interpret this SED as reprocessed emission in a 8600 $R_{\odot}$ circumstellar shell of dusty material with a mass $\sim$5$\times10^{-5} M_{\odot}$ surrounding a $\log(L/L_{\odot})=4.74\pm0.07$ and $T_{\rm eff}=3920\substack{+200\\-160}$ K RSG. This luminosity is consistent with RSG models of initial mass 11 $M_{\odot}$, depending on assumptions of rotation and overshooting. In addition, the counterpart was significantly variable in pre-explosion Spitzer 3.6 $μ$m and 4.5 $μ$m imaging, exhibiting $\sim$70% variability in both bands correlated across 9 yr and 29 epochs of imaging. The variations appear to have a timescale of 2.8 yr, which is consistent with $κ$-mechanism pulsations observed in RSGs, albeit with a much larger amplitude than RSGs such as $α$ Orionis (Betelgeuse). △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 14 pages, 5 figures, submitted to ApJL, comments welcome

Proposed host galaxies of repeating fast radio burst sources detected by CHIME/FRB

Authors: Adaeze L. Ibik, Maria R. Drout, B. M. Gaensler, Paul Scholz, Daniele Michilli, Mohit Bhardwaj, Victoria M. Kaspi, Ziggy Pleunis, Tomas Cassanelli, Amanda M. Cook, Fengqiu A. Dong, Calvin Leung, Kiyoshi W. Masui, Jane F. Kaczmarek, Katherine J. Lu, Aaron B. Pearlman, Masoud Rafiei-Ravandi, Ketan R. Sand, Kaitlyn Shin, Kendrick M. Smith, Ingrid H. Stairs

Abstract: We present a search for host galaxy associations for the third set of repeating fast radio burst (FRB) sources discovered by the CHIME/FRB Collaboration. Using the $\sim$ 1 arcmin CHIME/FRB baseband localizations and probabilistic methods, we identify potential host galaxies of two FRBs, 20200223B and 20190110C at redshifts of 0.06024(2) and 0.12244(6), respectively. We also discuss the properties… ▽ More We present a search for host galaxy associations for the third set of repeating fast radio burst (FRB) sources discovered by the CHIME/FRB Collaboration. Using the $\sim$ 1 arcmin CHIME/FRB baseband localizations and probabilistic methods, we identify potential host galaxies of two FRBs, 20200223B and 20190110C at redshifts of 0.06024(2) and 0.12244(6), respectively. We also discuss the properties of a third marginal candidate host galaxy association for FRB 20191106C with a host redshift of 0.10775(1). The three putative host galaxies are all relatively massive, fall on the standard mass-metallicity relationship for nearby galaxies, and show evidence of ongoing star formation. They also all show signatures of being in a transitional regime, falling in the ``green valley'' which is between the bulk of star-forming and quiescent galaxies. The plausible host galaxies identified by our analysis are consistent with the overall population of repeating and non-repeating FRB hosts while increasing the fraction of massive and bright galaxies. Coupled with these previous host associations, we identify a possible excess of FRB repeaters whose host galaxies have $M_{\mathrm{u}}-M_{\mathrm{r}}$ colors redder than the bulk of star-forming galaxies. Additional precise localizations are required to confirm this trend. △ Less

Submitted 2 October, 2023; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: 23 pages, 9 figures, submitted to AAS journals

Origin of high-velocity ejecta and early red excess emission in the infant Type Ia supernova 2021aefx

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Christopher D. Matzner, Kelvin C. C. Leong, Sang Chul Kim, Hong Soo Park, Youngdae Lee

Abstract: \object{SN 2021aefx} is a normal Type Ia Supernova (SN) with red excess emission over the first $\sim$ 2 days. We present detailed analysis of this SN using our high-cadence KMTNet multi-band photometry, spectroscopy, and publicly available data. We provide the first measurements of its epochs of explosion (MJD 59529.32 $\pm$ 0.16) as well as ``first light'' (MJD 59529.85 $\pm$ 0.55) associated wi… ▽ More \object{SN 2021aefx} is a normal Type Ia Supernova (SN) with red excess emission over the first $\sim$ 2 days. We present detailed analysis of this SN using our high-cadence KMTNet multi-band photometry, spectroscopy, and publicly available data. We provide the first measurements of its epochs of explosion (MJD 59529.32 $\pm$ 0.16) as well as ``first light'' (MJD 59529.85 $\pm$ 0.55) associated with the main ejecta ${\rm{^{56}Ni}}$ distribution. This places our first detection of SN 2021aefx at $\sim -$0.5 hours since ``first light'', indicating the presence of additional power sources. Our peak-spectrum confirms its Type Ia sub-classification as intermediate between Core-Normal and Broad-Line, and we estimate the ejecta mass to be $\sim$ 1.34 $M_{\odot}$. The pre-peak spectral evolution identifies fast-expanding material reaching $>$ 40,000 km s$^{-1}$ (the fastest ever observed in Type Ia SNe) and at least two distinct homologously-expanding ejecta components: (1) a normal-velocity (12,400 km s$^{-1}$) component consistent with the typical photospheric evolution of Chandrasekhar-mass ejecta; and (2) a high-velocity (23,500 km s$^{-1}$) component visible during the first $\sim$ 3.6 days post-explosion, which locates the component within the outer $<$ 16\% of the ejecta mass. Asymmetric, subsonic explosion processes producing a non-spherical photosphere provide an explanation for the simultaneous presence of the two components, as well as the red excess emission via a slight ${\rm{^{56}Ni}}$ enrichment in the outer $\sim$ 0.5\% of the ejecta mass. Our spectrum from 300 days post-peak advances the constraint against non-degenerate companions and further supports a near-Chandrasekhar-mass explosion origin. Off-center ignited delayed-detonations of Chandrasekhar-mass white dwarfs may be responsible for the observed features of SN 2021aefx in some normal Type Ia SNe. △ Less

Submitted 5 April, 2023; v1 submitted 2 April, 2023; originally announced April 2023.

Comments: Submitted for publication in ApJ. 29 pages, 14 figures, 4 tables

Late-time Hubble Space Telescope Observations of AT 2018cow. II. Evolution of a UV bright Underlying Source 2--4 Yr Post-discovery

Authors: Yuyang Chen, Maria R. Drout, Anthony L. Piro, Charles D. Kilpatrick, Ryan J. Foley, César Rojas-Bravo, M. R. Magee

Abstract: In this second of a two-paper series, we present a detailed analysis of three HST observations taken $\sim$2--4 years post-discovery, examining the evolution of a UV-bright underlying source at the precise position of AT 2018cow. While observations at $\sim$2--3 years post-discovery revealed an exceptionally blue ($L_ν\propto ν^{1.99}$) underlying source with relatively stable optical brightness,… ▽ More In this second of a two-paper series, we present a detailed analysis of three HST observations taken $\sim$2--4 years post-discovery, examining the evolution of a UV-bright underlying source at the precise position of AT 2018cow. While observations at $\sim$2--3 years post-discovery revealed an exceptionally blue ($L_ν\propto ν^{1.99}$) underlying source with relatively stable optical brightness, fading in the NUV was observed at year 4, indicating flattening in the spectrum (to $L_ν\propto ν^{1.64}$). The resulting spectral energy distributions can be described by an extremely hot but small blackbody, and the fading may be intrinsic (cooling) or extrinsic (increased absorption). Considering possible scenarios and explanations, we disfavor significant contributions from stellar sources and dust formation based on the observed color and brightness. By comparing the expected power and the observed luminosity, we rule out interaction with the known radio-producing circumstellar material as well as magnetar spin down with $B\sim10^{15}\,\mathrm{G}$ as possible power sources, though we cannot rule out the possible existence of a denser CSM component (e.g., previously ejected hydrogen envelope) or a magnetar with $B\lesssim10^{14}\,\mathrm{G}$. Finally, we find that a highly-inclined precessing accretion disk can reasonably explain the color, brightness, and evolution of the underlying source. However, a major uncertainty in this scenario is the mass of the central black hole (BH), as both stellar-mass and intermediate-mass BHs face notable challenges that cannot be explained by our simple disk model, and further observations and theoretical works are needed to fully constrain the nature of this underlying source. △ Less

Submitted 28 September, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 24 pages, 12 figures. This is paper II of a two-paper series. v2: updated to the accepted version. Note a mistake in v1 in the caption of Figure 6 is now fixed: removed the word "dereddened" because the plotted observed magnitudes are directly from Table 1 which are not extinction-corrected (the analyses are not affected by this)

Late-time Hubble Space Telescope Observations of AT 2018cow. I. Further Constraints on the Fading Prompt Emission and Thermal Properties 50--60 Days Post-discovery

Authors: Yuyang Chen, Maria R. Drout, Anthony L. Piro, Charles D. Kilpatrick, Ryan J. Foley, César Rojas-Bravo, Kirsty Taggart, Matthew R. Siebert, M. R. Magee

Abstract: The exact nature of the luminous fast blue optical transient AT 2018cow is still debated. In this first of a two-paper series, we present a detailed analysis of three Hubble Space Telescope (HST) observations of AT 2018cow covering $\sim$50--60 days post-discovery in combination with other observations throughout the first two months, and derive significantly improved constraints of the late therm… ▽ More The exact nature of the luminous fast blue optical transient AT 2018cow is still debated. In this first of a two-paper series, we present a detailed analysis of three Hubble Space Telescope (HST) observations of AT 2018cow covering $\sim$50--60 days post-discovery in combination with other observations throughout the first two months, and derive significantly improved constraints of the late thermal properties. By modeling the spectral energy distributions (SEDs), we confirm that the UV-optical emission over 50--60 days was still a smooth blackbody (i.e., optically thick) with a high temperature ($T_{\mathrm{BB}}\sim15000\,\mathrm{K}$) and small radius ($R_{\mathrm{BB}}\lesssim1000\,R_\odot$). Additionally, we report for the first time a break in the bolometric light curve: the thermal luminosity initially declined at a rate of $L_{\mathrm{BB}}\propto t^{-2.40}$, but faded much faster at $t^{-3.06}$ after day 13. Reexamining possible late-time power sources, we disfavor significant contributions from radioactive decay based on the required $^{56}$Ni mass and lack of UV line blanketing in the HST SEDs. We argue that the commonly proposed interaction with circumstellar material may face significant challenges in explaining the late thermal properties, particularly the effects of the optical depth. Alternatively, we find that continuous outflow/wind driven by a central engine can still reasonably explain the combination of a receding photosphere, optically thick and rapidly fading emission, and intermediate-width lines. However, the rapid fading may have further implications on the power output and structure of the system. Our findings may support the hypothesis that AT 2018cow and other "Cow-like transients" are powered mainly by accretion onto a central engine. △ Less

Submitted 28 September, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 22 pages, 10 figures. This is paper I of a two-paper series. v2: updated to the accepted version

Cool, Luminous, and Highly Variable Stars in the Magellanic Clouds. II: Spectroscopic and Environmental Analysis of Thorne-Żytkow Object and Super-AGB Star Candidates

Authors: Anna J. G. O'Grady, Maria R. Drout, B. M. Gaensler, C. S. Kochanek, Kathryn F. Neugent, Carolyn L. Doherty, Joshua S. Speagle, B. J. Shappee, Michael Rauch, Ylva Götberg, Bethany Ludwig, Todd A. Thompson

Abstract: In previous work we identified a population of 38 cool and luminous variable stars in the Magellanic Clouds and examined 11 in detail in order to classify them as either Thorne-Żytkow Objects (TŻOs, red supergiants with a neutron star cores) or super-AGB stars (the most massive stars that will not undergo core collapse). This population includes HV\,2112, a peculiar star previously considered in o… ▽ More In previous work we identified a population of 38 cool and luminous variable stars in the Magellanic Clouds and examined 11 in detail in order to classify them as either Thorne-Żytkow Objects (TŻOs, red supergiants with a neutron star cores) or super-AGB stars (the most massive stars that will not undergo core collapse). This population includes HV\,2112, a peculiar star previously considered in other works to be either a TŻO or high-mass AGB star. Here we continue this investigation, using the kinematic and radio environments and local star formation history of these stars to place constraints on the age of the progenitor systems and the presence of past supernovae. These stars are not associated with regions of recent star formation, and we find no evidence of past supernovae at their locations. Finally, we also assess the presence of heavy elements and lithium in their spectra compared to red supergiants. We find strong absorption in Li and s-process elements compared to RSGs in most of the sample, consistent with super-AGB nucleosynthesis, while HV\,2112 shows additional strong lines associated with TŻO nucleosynthesis. Coupled with our previous mass estimates, the results are consistent with the stars being massive (~4-6.5M$_{\odot}$) or super-AGB (~6.5-12M$_{\odot}$) stars in the thermally pulsing phase, providing crucial observations of the transition between low- and high-mass stellar populations. HV\,2112 is more ambiguous; it could either be a maximally massive sAGB star, or a TŻO if the minimum mass for stability extends down to <13 M$_\odot$. △ Less

Submitted 22 November, 2022; originally announced November 2022.

Comments: 32 pages, 10 figures, 7 tables, spectroscopic data available at https://zenodo.org/record/7058608, accepted to The Astrophysical Journal

The Young Supernova Experiment Data Release 1 (YSE DR1): Light Curves and Photometric Classification of 1975 Supernovae

Authors: P. D. Aleo, K. Malanchev, S. Sharief, D. O. Jones, G. Narayan, R. J. Foley, V. A. Villar, C. R. Angus, V. F. Baldassare, M. J. Bustamante-Rosell, D. Chatterjee, C. Cold, D. A. Coulter, K. W. Davis, S. Dhawan, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, J. Hjorth, M. E. Huber, W. V. Jacobson-Galán, C. D. Kilpatrick, D. Langeroodi , et al. (58 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 21 February, 2023; v1 submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted to ApJS; 64 pages; 35 figures; 10 tables

Spectroscopic r-Process Abundance Retrieval for Kilonovae I: The Inferred Abundance Pattern of Early Emission from GW170817

Authors: Nicholas Vieira, John J. Ruan, Daryl Haggard, Nicole Ford, Maria R. Drout, Rodrigo Fernández, N. R. Badnell

Abstract: Freshly-synthesized r-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features encode insights into the physical conditions of the r-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computation… ▽ More Freshly-synthesized r-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features encode insights into the physical conditions of the r-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computationally challenging. We introduce Spectroscopic r-Process Abundance Retrieval for Kilonovae (SPARK), a modular framework to perform Bayesian inference on kilonova spectra with the goals of inferring elemental abundance patterns and identifying absorption features at early times. SPARK inputs an atomic line list and abundance patterns from reaction network calculations into the TARDIS radiative transfer code. It then performs fast Bayesian inference on observed kilonova spectra by training a Gaussian process surrogate for the approximate posteriors of kilonova ejecta parameters, via active learning. We use the spectrum of GW170817 at 1.4 days to perform the first inference on a kilonova spectrum, and recover a complete abundance pattern. Our inference shows that this ejecta was generated by an r-process with either (1) high electron fraction Y_e ~ 0.35 and high entropy s/k_B ~ 25, or, (2) a more moderate Y_e ~ 0.30 and s/k_B ~ 14. These parameters are consistent with a shocked, polar dynamical component, and a viscously-driven outflow from a remnant accretion disk, respectively. We also recover previous identifications of strontium absorption at ~8000 AA, and tentatively identify yttrium and/or zirconium at < 4500 AA. Our approach will enable computationally-tractable inference on the spectra of future kilonovae discovered through multi-messenger observations. △ Less

Submitted 22 December, 2022; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: Accepted to ApJ; 34 pages, 13 figures

A fast rising tidal disruption event from a candidate intermediate mass black hole

Authors: C. R. Angus, V. F. Baldassare, B. Mockler, R. J. Foley, E. Ramirez-Ruiz, S. I. Raimundo, K. D. French, K. Auchettl, H. Pfister, C. Gall, J. Hjorth, M. R. Drout, K. D. Alexander, G. Dimitriadis, T. Hung, D. O. Jones, A. Rest, M. R. Siebert, K. Taggart, G. Terreran, S. Tinyanont, C. M. Carroll, L. DeMarchi, N. Earl, A. Gagliano , et al. (14 additional authors not shown)

Abstract: 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 emit… ▽ More 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. △ Less

Submitted 5 September, 2022; v1 submitted 31 August, 2022; originally announced September 2022.

Comments: Accepted for publication in Nature Astronomy

A low-cost ultraviolet-to-infrared absolute quantum efficiency characterization system of detectors

Authors: Ajay S. Gill, Mohamed M. Shaaban, Aaron Tohuvavohu, Suresh Sivanandam, Roberto G. Abraham, Seery Chen, Maria R. Drout, Deborah Lokhorst, Christopher D. Matzner, Stefan W. Mochnacki, Calvin B. Netterfield

Abstract: We present a low-cost ultraviolet to infrared absolute quantum efficiency detector characterization system developed using commercial off-the-shelf components. The key components of the experiment include a light source,a regulated power supply, a monochromator, an integrating sphere, and a calibrated photodiode. We provide a step-by-step procedure to construct the photon and quantum efficiency tr… ▽ More We present a low-cost ultraviolet to infrared absolute quantum efficiency detector characterization system developed using commercial off-the-shelf components. The key components of the experiment include a light source,a regulated power supply, a monochromator, an integrating sphere, and a calibrated photodiode. We provide a step-by-step procedure to construct the photon and quantum efficiency transfer curves of imaging sensors. We present results for the GSENSE 2020 BSI CMOS sensor and the Sony IMX 455 BSI CMOS sensor. As a reference for similar characterizations, we provide a list of parts and associated costs along with images of our setup. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: 15 pages, 14 figures, SPIE Astronomical Telescopes and Instrumentation, Conference 12191: X-ray, Optical, and Infrared Detectors for Astronomy X, Paper Number: 12191-39 (Montreal, July 2022)

The origin and evolution of the normal Type Ia SN 2018aoz with infant-phase reddening and excess emission

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Abigail Polin, David J. Sand, Santiago GonzÁlez-GaitÁn, Sang Chul Kim, Youngdae Lee, Hong Soo Park, D. Andrew Howell, Peter E. Nugent, Anthony L. Piro, Peter J. Brown, LluÍs Galbany, Jamison Burke, Daichi Hiramatsu, Griffin Hosseinzadeh, Stefano Valenti, Niloufar Afsariardchi, Jennifer E. Andrews, John Antoniadis, Rachael L. Beaton, K. Azalee Bostroem, Raymond G. Carlberg, S. Bradley Cenko , et al. (18 additional authors not shown)

Abstract: SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is int… ▽ More SN~2018aoz is a Type Ia SN with a $B$-band plateau and excess emission in the infant-phase light curves $\lesssim$ 1 day after first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN~2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show the SN is intermediate between two subtypes of normal Type Ia: Core-Normal and Broad-Line. The excess emission could have contributions from the radioactive decay of surface iron-peak elements as well as ejecta interaction with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on H$α$ and He~I favour a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O~I] and He~I disfavours a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN~2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe~II] and [Ni~II]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1-D models are incompatible with the infant-phase excess emission, $B_{\rm max}-V_{\rm max}$ color, and absence of nebular-phase [Ca~II]. Although the explosion processes of SN~2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear normal after $\sim$ 1 day. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: Submitted for publication in ApJ. 35 pages, 16 figures, 7 tables

Infant-phase reddening by surface Fe-peak elements in a normal Type Ia Supernova

Authors: Yuan Qi Ni, Dae-Sik Moon, Maria R. Drout, Abigail Polin, David J. Sand, Santiago Gonzalez-Gaitan, Sang Chul Kim, Youngdae Lee, Hong Soo Park, D. Andrew Howell, Peter E. Nugent, Anthony L. Piro, Peter J. Brown, Lluis Galbany, Jamison Burke, Daichi Hiramatsu, Griffin Hosseinzadeh, Stefano Valenti, Niloufar Afsariardchi, Jennifer E. Andrews, John Antoniadis, Iair Arcavi, Rachael L. Beaton, K. Azalee Bostroem, Raymond G. Carlberg , et al. (19 additional authors not shown)

Abstract: Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of… ▽ More Type Ia Supernovae are thermonuclear explosions of white dwarf stars. They play a central role in the chemical evolution of the Universe and are an important measure of cosmological distances. However, outstanding questions remain about their origins. Despite extensive efforts to obtain natal information from their earliest signals, observations have thus far failed to identify how the majority of them explode. Here, we present infant-phase detections of SN 2018aoz from a brightness of -10.5 absolute AB magnitudes -- the lowest luminosity early Type Ia signals ever detected -- revealing a hitherto unseen plateau in the $B$-band that results in a rapid redward color evolution between 1.0 and 12.4 hours after the estimated epoch of first light. The missing $B$-band flux is best-explained by line-blanket absorption from Fe-peak elements in the outer 1% of the ejected mass. The observed $B-V$ color evolution of the SN also matches the prediction from an over-density of Fe-peak elements in the same outer 1% of the ejected mass, whereas bluer colors are expected from a purely monotonic distribution of Fe-peak elements. The presence of excess nucleosynthetic material in the extreme outer layers of the ejecta points to enhanced surface nuclear burning or extended sub-sonic mixing processes in some normal Type Ia Supernova explosions. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Nature Astronomy. Main text = 8 pages, 4 figures, 1 table; Full document = 46 pages, with Methods, Supplementary Information, 7 Supplementary figures, 2 Supplementary tables and references. Nat Astron (2022)

Discovery of a Short-Period and Unusually Helium-Deficient Dwarf Nova KSP-OT-201701a by the KMTNet Supernova Program

Authors: Youngdae Lee, Sang Chul Kim, Dae-Sik Moon, Hong Soo Park, Maria R. Drout, Yuan Qi Ni, Hyobin Im

Abstract: We present the first ever discovery of a short-period and unusually helium-deficient dwarf nova KSP-OT-201701a by the Korea Microlensing Telescope Network Supernova Program. The source shows three superoutbursts, each led by a precursor outburst, and several normal outbursts in BVI during the span of ~2.6 years with supercycle and normal cycle lengths of about 360 and 76 days, respectively. Spectr… ▽ More We present the first ever discovery of a short-period and unusually helium-deficient dwarf nova KSP-OT-201701a by the Korea Microlensing Telescope Network Supernova Program. The source shows three superoutbursts, each led by a precursor outburst, and several normal outbursts in BVI during the span of ~2.6 years with supercycle and normal cycle lengths of about 360 and 76 days, respectively. Spectroscopic observations near the end of a superoutburst reveal the presence of strong double-peaked HI emission lines together with weak HeI emission lines. The helium-to-hydrogen intensity ratios measured by HeIλ5876 and Hα lines are 0.10 {\pm} 0.01 at a quiescent phase and 0.26 {\pm} 0.04 at an outburst phase, similar to the ratios found in long-period dwarf novae while significantly lower than those in helium cataclysmic variables (He CVs). Its orbital period of 51.91 {\pm} 2.50 minutes, which is estimated based on time series spectroscopy, is a bit shorter than the superhump period of 56.52 {\pm} 0.19 minutes, as expected from the gravitational interaction between the eccentric disk and the secondary star. We measure its mass ratio to be 0.37^{+0.32}_{-0.21} using the superhump period excess of 0.089 {\pm} 0.053. The short orbital period, which is under the period minimum, the unusual helium deficiency, and the large mass ratio suggest that KSP-OT-201701a is a transition object evolving to a He CV from a long-period dwarf nova with an evolved secondary star. △ Less

Submitted 3 February, 2022; originally announced February 2022.

Comments: 11 pages, 5 figures, published in ApJL

Journal ref: https://iopscience.iop.org/article/10.3847/2041-8213/ac4c41

Searching for Anomalies in the ZTF Catalog of Periodic Variable Stars

Authors: H. S. Chan, V. Ashley Villar, S. H. Cheung, Shirley Ho, Anna J. G. O'Grady, Maria R. Drout, Mathieu Renzo

Abstract: Periodic variables illuminate the physical processes of stars throughout their lifetime. Wide-field surveys continue to increase our discovery rates of periodic variable stars. Automated approaches are essential to identify interesting periodic variable stars for multi-wavelength and spectroscopic follow-up. Here, we present a novel unsupervised machine learning approach to hunt for anomalous peri… ▽ More Periodic variables illuminate the physical processes of stars throughout their lifetime. Wide-field surveys continue to increase our discovery rates of periodic variable stars. Automated approaches are essential to identify interesting periodic variable stars for multi-wavelength and spectroscopic follow-up. Here, we present a novel unsupervised machine learning approach to hunt for anomalous periodic variables using phase-folded light curves presented in the Zwicky Transient Facility Catalogue of Periodic Variable Stars by \citet{Chen_2020}. We use a convolutional variational autoencoder to learn a low dimensional latent representation, and we search for anomalies within this latent dimension via an isolation forest. We identify anomalies with irregular variability. Most of the top anomalies are likely highly variable Red Giants or Asymptotic Giant Branch stars concentrated in the Milky Way galactic disk; a fraction of the identified anomalies are more consistent with Young Stellar Objects. Detailed spectroscopic follow-up observations are encouraged to reveal the nature of these anomalies. △ Less

Submitted 6 December, 2021; originally announced December 2021.

Comments: 26 pages, 17 figures. The full version of Table 4 and Table 5 are available upon request

Revised Stellar Parameters for V471 Tau, A Post-common Envelope Binary in the Hyades

Authors: Philip S. Muirhead, Jason Nordhaus, Maria R. Drout

Abstract: V471 Tau is a post-common-envelope binary consisting of an eclipsing DA white dwarf and a K-type main-sequence star in the Hyades star cluster. We analyzed publicly available photometry and spectroscopy of V471 Tau to revise the stellar and orbital parameters of the system. We used archival K2 photometry, archival Hubble Space Telescope spectroscopy, and published radial-velocity measurements of t… ▽ More V471 Tau is a post-common-envelope binary consisting of an eclipsing DA white dwarf and a K-type main-sequence star in the Hyades star cluster. We analyzed publicly available photometry and spectroscopy of V471 Tau to revise the stellar and orbital parameters of the system. We used archival K2 photometry, archival Hubble Space Telescope spectroscopy, and published radial-velocity measurements of the K-type star. Employing Gaussian processes to fit for rotational modulation of the system flux by the main-sequence star, we recovered the transits of the white dwarf in front of the main-sequence star for the first time. The transits are shallower than would be expected from purely geometric occultations owing to gravitational microlensing during transit, which places an additional constraint on the white-dwarf mass. Our revised mass and radius for the main-sequence star is consistent with single-star evolutionary models given the age and metallicity of the Hyades. However, as noted previously in the literature, the white dwarf is too massive and too hot to be the result of single-star evolution given the age of the Hyades, and may be the product of a merger scenario. We independently estimate the conditions of the system at the time of common envelope that would result in the measured orbital parameters today. △ Less

Submitted 30 November, 2021; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: 11 pages, 7 figures, Accepted to The Astronomical Journal

Final Moments I: Precursor Emission, Envelope Inflation, and Enhanced Mass loss Preceding the Luminous Type II Supernova 2020tlf

Authors: Wynn Jacobson-Galán, Luc Dessart, David Jones, Raffaella Margutti, Deanne Coppejans, Georgios Dimitriadis, Ryan J. Foley, Charles D. Kilpatrick, David J. Matthews, Sofia Rest, Giacomo Terreran, Patrick D. Aleo, Katie Auchettl, Peter K. Blanchard, David A. Coulter, Kyle W. Davis, Thomas de Boer, Lindsay DeMarchi, Maria R. Drout, Nicholas Earl, Alexander Gagliano, Christa Gall, Jens Hjorth, Mark E. Huber, Adaeze L. Ibik , et al. (11 additional authors not shown)

Abstract: We present panchromatic observations and modeling of supernova (SN) 2020tlf, the first normal type II-P/L SN with confirmed precursor emission, as detected by the Young Supernova Experiment transient survey with the Pan-STARRS1 telescope. Pre-explosion emission was detected in $riz-$bands at 130 days prior to SN 2020tlf and persisted at relatively constant flux until first light. Soon after discov… ▽ More We present panchromatic observations and modeling of supernova (SN) 2020tlf, the first normal type II-P/L SN with confirmed precursor emission, as detected by the Young Supernova Experiment transient survey with the Pan-STARRS1 telescope. Pre-explosion emission was detected in $riz-$bands at 130 days prior to SN 2020tlf and persisted at relatively constant flux until first light. Soon after discovery, "flash" spectroscopy of SN 2020tlf revealed prominent narrow symmetric emission lines ($v_w < 300$ km s$^{-1}$) that resulted from the photo-ionization of unshocked circumstellar material (CSM) shedded in progenitor mass loss episodes in the final weeks to months before explosion. Surprisingly, this novel display of pre-SN emission and associated mass loss occurred in a RSG progenitor with ZAMS mass of only 10-12 M$_{\odot}$, as inferred from nebular spectra. Modeling of the light curve and multi-epoch spectra with the non-LTE radiative transfer code CMFGEN and radiation-hydrodynamical (RHD) code HERACLES suggests a dense CSM limited to $r \approx 10^{15}$ cm, and mass loss rate of $10^{-2}$ M$_{\odot}$ yr$^{-1}$. The subsequent luminous light-curve plateau and persistent blue excess indicates an extended progenitor, compatible with a RSG model with $R_{\star} = 1100$ R$_{\odot}$. Inferences from the limits on the shock-powered X-ray and radio luminosity are consistent with these conclusions and suggest a CSM density of $ρ< 2 \times 10^{-16}$ g cm$^{-3}$ for distances of $r \approx 5 \times 10^{15}$ cm, as well as a mass loss rate of $\dot M<1.3 \times 10^{-5}\,\rm{M_{\odot}\,yr^{-1}}$ at larger distances. A promising power source for the observed precursor emission is the ejection of stellar material following energy disposition into the stellar envelope as a result of gravity waves emitted during either neon/oxygen burning or a nuclear flash from silicon combustion. △ Less

Submitted 13 December, 2021; v1 submitted 24 September, 2021; originally announced September 2021.

Comments: 36 pages, 20 figures. Accepted for publication in ApJ

Final Targeting Strategy for the SDSS-IV APOGEE-2S Survey

Authors: Felipe A. Santana, Rachael L. Beaton, Kevin R. Covey, Julia E. O'Connell, Penélope Longa-Peña, Roger Cohen, José G. Fernández-Trincado, Christian R. Hayes, Gail Zasowski, Jennifer S. Sobeck, Steven R. Majewski, S. D. Chojnowski, Nathan De Lee, Ryan J. Oelkers, Guy S. Stringfellow, Andrés Almeida, Borja Anguiano, John Donor, Peter M. Frinchaboy, Sten Hasselquist, Jennifer A. Johnson, Juna A. Kollmeier, David L. Nidever, Adrian. M. Price-Whelan, Alvaro Rojas-Arriagada , et al. (21 additional authors not shown)

Abstract: APOGEE is a high-resolution (R sim 22,000), near-infrared, multi-epoch, spectroscopic survey of the Milky Way. The second generation of the APOGEE project, APOGEE-2, includes an expansion of the survey to the Southern Hemisphere called APOGEE-2S. This expansion enabled APOGEE to perform a fully panoramic mapping of all the main regions of the Milky Way; in particular, by operating in the H-band, A… ▽ More APOGEE is a high-resolution (R sim 22,000), near-infrared, multi-epoch, spectroscopic survey of the Milky Way. The second generation of the APOGEE project, APOGEE-2, includes an expansion of the survey to the Southern Hemisphere called APOGEE-2S. This expansion enabled APOGEE to perform a fully panoramic mapping of all the main regions of the Milky Way; in particular, by operating in the H-band, APOGEE is uniquely able to probe the dust-hidden inner regions of the Milky Way that are best accessed from the Southern Hemisphere. In this paper we present the targeting strategy of APOGEE-2S, with special attention to documenting modifications to the original, previously published plan. The motivation for these changes is explained as well as an assessment of their effectiveness in achieving their intended scientific objective. In anticipation of this being the last paper detailing APOGEE targeting, we present an accounting of all such information complete through the end of the APOGEE-2S project; this includes several main survey programs dedicated to exploration of major stellar populations and regions of the Milky Way, as well as a full list of programs contributing to the APOGEE database through allocations of observing time by the Chilean National Time Allocation Committee (CNTAC) and the Carnegie Institution for Science (CIS). This work was presented along with a companion article, R. Beaton et al. (submitted; AAS29028), presenting the final target selection strategy adopted for APOGEE-2 in the Northern Hemisphere. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Comments: 43 pages; 5 figures; 6 Tables; 1 Appendix; Submitted to Journal and Under Review; Posting to accompany papers using the SDSS-IV/APOGEE-2 Data Release 17 scheduled for December 2021

An Early-Time Optical and Ultraviolet Excess in the type-Ic SN 2020oi

Authors: Alexander Gagliano, Luca Izzo, Charles D. Kilpatrick, Brenna Mockler, Wynn Vincente Jacobson-Galán, Giacomo Terreran, Georgios Dimitriadis, Yossef Zenati, Katie Auchettl, Maria R. Drout, Gautham Narayan, Ryan J. Foley, R. Margutti, Armin Rest, D. O. Jones, Christian Aganze, Patrick D. Aleo, Adam J. Burgasser, D. A. Coulter, Roman Gerasimov, Christa Gall, Jens Hjorth, Chih-Chun Hsu, Eugene A. Magnier, Kaisey S. Mandel , et al. (8 additional authors not shown)

Abstract: We present photometric and spectroscopic observations of Supernova 2020oi (SN 2020oi), a nearby ($\sim$17 Mpc) type-Ic supernova (SN Ic) within the grand-design spiral M100. We undertake a comprehensive analysis to characterize the evolution of SN 2020oi and constrain its progenitor system. We detect flux in excess of the fireball rise model $δt \approx 2.5$ days from the date of explosion in mult… ▽ More We present photometric and spectroscopic observations of Supernova 2020oi (SN 2020oi), a nearby ($\sim$17 Mpc) type-Ic supernova (SN Ic) within the grand-design spiral M100. We undertake a comprehensive analysis to characterize the evolution of SN 2020oi and constrain its progenitor system. We detect flux in excess of the fireball rise model $δt \approx 2.5$ days from the date of explosion in multi-band optical and UV photometry from the Las Cumbres Observatory and the Neil Gehrels Swift Observatory, respectively. The derived SN bolometric luminosity is consistent with an explosion with $M_{\rm ej} = 0.81 \pm 0.03 M_{\odot}$, $E_{k}= 0.79 \pm 0.09 \times 10^{51} \rm{erg} \rm{s}^{-1}$, and $M_{\rm Ni56} = 0.08 \pm 0.02 M_{\odot}$. Inspection of the event's decline reveals the highest $Δm_{15,\rm{bol}}$ reported for a stripped-envelope event to date. Modeling of optical spectra near event peak indicates a partially mixed ejecta comparable in composition to the ejecta observed in SN 1994I, while the earliest spectrum shows signatures of a possible interaction with material of a distinct composition surrounding the SN progenitor. Further, Hubble Space Telescope (HST) pre-explosion imaging reveals a stellar cluster coincident with the event. From the cluster photometry, we derive the mass and age of the SN progenitor using stellar evolution models implemented in the BPASS library. Our results indicate that SN 2020oi occurred in a binary system from a progenitor of mass $M_{\rm ZAMS} \approx 9.5 \pm 1.0 M_{\odot}$, corresponding to an age of $27 \pm 7$ Myr. SN 2020oi is the dimmest SN Ic event to date for which an early-time flux excess has been observed, and the first in which an early excess is unlikely to be associated with shock-cooling. △ Less

Submitted 1 November, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: 49 pages, 21 figures; accepted to ApJ

Rapidly-Declining Hostless Type Ia Supernova KSP-OT-201509b from the KMTNet Supernova Program: Transitional Nature and Constraint on $^{56}$Ni Distribution and Progenitor Type

Authors: Dae-Sik Moon, Yuan Qi Ni, Maria R. Drout, Santiago González-Gaitán, Niloufar Afsariardchi, Hong Soo Park, Youngdae Lee, Sang Chul Kim, John Antoniadis, Dong-Jin Kim, Yongseok Lee

Abstract: We report the early discovery and multi-color ($BVI$) high-cadence light curve analyses of a rapidly-declining sub-Chandrasekhar Type Ia supernova KSP-OT-201509b (= AT2015cx) from the KMTNet Supernova Program. The Phillips parameter and color stretch parameter of KSP-OT-201509b (= AT2015cx) are $Δ$$M_{B,15}$ $\simeq$ 1.62 mag and $s_{BV}$ $\simeq$ 0.54, respectively, at an inferred redshift of 0.0… ▽ More We report the early discovery and multi-color ($BVI$) high-cadence light curve analyses of a rapidly-declining sub-Chandrasekhar Type Ia supernova KSP-OT-201509b (= AT2015cx) from the KMTNet Supernova Program. The Phillips parameter and color stretch parameter of KSP-OT-201509b (= AT2015cx) are $Δ$$M_{B,15}$ $\simeq$ 1.62 mag and $s_{BV}$ $\simeq$ 0.54, respectively, at an inferred redshift of 0.072. These, together with other measured parameters (such as the strength of the secondary $I$-band peak, colors and luminosity), identify the source to be a rapidly-declining Type Ia of transitional nature that is closer to Branch Normal than 91bg-like. Its early light curve evolution and bolometric luminosity are consistent with those of homologously expanding ejecta powered by radioactive decay and a Type Ia SN explosion with 0.32 $\pm$ 0.01 $M_{\odot}$ of synthesized $^{56}$Ni mass, 0.84 $\pm$ 0.12 $M_{\odot}$ of ejecta mass and (0.61 $\pm$ 0.14) $\times$ 10$^{51}$ erg of ejecta kinetic energy. While its $B-V$ and $V-I$ colors evolve largely synchronously with the changes in the $I$-band light curve as found in other supernovae, we also find the presence of an early redward evolution in $V-I$ prior to --10 days since peak. The bolometric light curve of the source is compatible with a stratified $^{56}$Ni distribution extended to shallow layers of the exploding progenitor. Comparisons between the observed light curves and those predicted from ejecta-companion interactions clearly disfavor Roche Lobe-filling companion stars at large separation distances, thus supporting a double-degenerate scenario for its origin. The lack of any apparent host galaxy in our deep stack images reaching a sensitivity limit of $\sim$ 28 $\rm mag\;arcsec^{-2}$ makes KSP-OT-201509b a hostless Type Ia supernova and offers new insights into supernova host galaxy environments. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: 39 pages, 20 figures, To appear in ApJ

FIRST J153350.8+272729: the radio afterglow of a decades-old tidal disruption event

Authors: Vikram Ravi, Hannah Dykaar, Jackson Codd, Ginevra Zaccagnini, Dillon Dong, Maria R. Drout, Bryan M. Gaensler, Gregg Hallinan, Casey Law

Abstract: We present the discovery of the fading radio transient FIRST J153350.8+272729. The source had a maximum observed 5-GHz radio luminosity of $8\times10^{39}$ erg s$^{-1}$ in 1986, but by 2019 had faded by a factor of nearly 400. It is located 0.15 arcsec from the center of a galaxy (SDSS J153350.89+272729) at 147 Mpc, which shows weak Type II Seyfert activity. We show that a tidal disruption event (… ▽ More We present the discovery of the fading radio transient FIRST J153350.8+272729. The source had a maximum observed 5-GHz radio luminosity of $8\times10^{39}$ erg s$^{-1}$ in 1986, but by 2019 had faded by a factor of nearly 400. It is located 0.15 arcsec from the center of a galaxy (SDSS J153350.89+272729) at 147 Mpc, which shows weak Type II Seyfert activity. We show that a tidal disruption event (TDE) is the preferred scenario for FIRST J153350.8+272729, although it could plausibly be interpreted as the afterglow of a long-duration gamma-ray burst. This is only the second TDE candidate to be first discovered at radio wavelengths. Its luminosity fills a gap between the radio afterglows of sub-relativistic TDEs in the local universe, and relativistic TDEs at high redshifts. The unusual properties of FIRST J153350.8+272729 (ongoing nuclear activity in the host galaxy, high radio luminosity) motivate more extensive TDE searches in untargeted radio surveys. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 10 pages, 3 figures, 1 table, submitted to AAS Journals (ApJ)

Gaia EDR3 confirms that Westerlund 1 is closer and older than previously thought

Authors: Mojgan Aghakhanloo, Jeremiah W. Murphy, Nathan Smith, John Parejko, Mariangelly Díaz-Rodríguez, Maria R. Drout, Jose H. Groh, Joseph Guzman, Keivan G. Stassun

Abstract: Using {\it Gaia} Early Data Release 3 (EDR3) parallaxes and Bayesian inference, we infer a parallax of the Westerlund 1 (Wd1) cluster. We find a parallax of $0.34\pm{0.05}$ mas corresponding to a distance of $2.8^{+0.7}_{-0.6}$ kpc. The new {\it Gaia} EDR3 distance is consistent with our previous result using {\it Gaia} DR2 parallaxes. This confirms that Wd1 is less massive and older than previous… ▽ More Using {\it Gaia} Early Data Release 3 (EDR3) parallaxes and Bayesian inference, we infer a parallax of the Westerlund 1 (Wd1) cluster. We find a parallax of $0.34\pm{0.05}$ mas corresponding to a distance of $2.8^{+0.7}_{-0.6}$ kpc. The new {\it Gaia} EDR3 distance is consistent with our previous result using {\it Gaia} DR2 parallaxes. This confirms that Wd1 is less massive and older than previously assumed. Compared to DR2, the EDR3 individual parallax uncertainties for each star decreased by 30\%. However, the aggregate parallax uncertainty for the cluster remained the same. This suggests that the uncertainty is dominated by systematics, which is possibly due to crowding, motions within the cluster, or motions due to binary orbits. △ Less

Submitted 15 January, 2021; originally announced January 2021.

Comments: 3 pages, 1 figure

The Center of Expansion and Age of the Oxygen-rich Supernova Remnant 1E 0102.2-7219

Authors: John Banovetz, Dan Milisavljevic, Niharika Sravan, Robert A. Fesen, Daniel J. Patnaude, Paul P. Plucinsky, William P. Blair, Kathryn E. Weil, Jon A. Morse, Raffaella Margutti, Maria R. Drout

Abstract: We present new proper motion measurements of optically emitting oxygen-rich knots of supernova remnant 1E 0102.2-7219 (E0102), which are used to estimate the remnant's center of expansion and age. Four epochs of high resolution Hubble Space Telescope images spanning 19 yr were retrieved and analyzed. We found a robust center of expansion of alpha=1:04:02.48 and delta=-72:01:53.92 (J2000) with 1-si… ▽ More We present new proper motion measurements of optically emitting oxygen-rich knots of supernova remnant 1E 0102.2-7219 (E0102), which are used to estimate the remnant's center of expansion and age. Four epochs of high resolution Hubble Space Telescope images spanning 19 yr were retrieved and analyzed. We found a robust center of expansion of alpha=1:04:02.48 and delta=-72:01:53.92 (J2000) with 1-sigma uncertainty of 1.77 arcseconds using 45 knots from images obtained with the Advanced Camera for Surveys using the F475W filter in 2003 and 2013 having the highest signal-to-noise ratio. We also estimate an upper limit explosion age of 1738 +/- 175 yr by selecting knots with the highest proper motions, that are assumed to be the least decelerated. We find evidence of an asymmetry in the proper motions of the knots as a function of position angle. We conclude that these asymmetries were most likely caused by interaction between E0102's original supernova blast wave and an inhomogeneous surrounding environment, as opposed to intrinsic explosion asymmetry. The observed non-homologous expansion suggests that the use of a free expansion model inaccurately offsets the center of expansion and leads to an overestimated explosion age. We discuss our findings as they compare to previous age and center of expansion estimates of E0102 and their relevance to a recently identified candidate central compact object. △ Less

Submitted 13 January, 2021; originally announced January 2021.

Comments: 14 pages, 11 figures, revised according to referee comments and resubmitted to ApJ

A Cool and Inflated Progenitor Candidate for the Type Ib Supernova 2019yvr at 2.6 Years Before Explosion

Authors: Charles D. Kilpatrick, Maria R. Drout, Katie Auchettl, Georgios Dimitriadis, Ryan J. Foley, David O. Jones, Lindsay DeMarchi, K. Decker French, Christa Gall, Jens Hjorth, Wynn V. Jacobson-Galan, Raffaella Margutti, Anthony L. Piro, Enrico Ramirez-Ruiz, Armin Rest, Cesar Rojas-Bravo

Abstract: We present Hubble Space Telescope imaging of a pre-explosion counterpart to SN 2019yvr obtained 2.6 years before its explosion as a type Ib supernova (SN Ib). Aligning to a post-explosion Gemini-S/GSAOI image, we demonstrate that there is a single source consistent with being the SN 2019yvr progenitor system, the second SN Ib progenitor candidate after iPTF13bvn. We also analyzed pre-explosion Spi… ▽ More We present Hubble Space Telescope imaging of a pre-explosion counterpart to SN 2019yvr obtained 2.6 years before its explosion as a type Ib supernova (SN Ib). Aligning to a post-explosion Gemini-S/GSAOI image, we demonstrate that there is a single source consistent with being the SN 2019yvr progenitor system, the second SN Ib progenitor candidate after iPTF13bvn. We also analyzed pre-explosion Spitzer/IRAC imaging, but we do not detect any counterparts at the SN location. SN 2019yvr was highly reddened, and comparing its spectra and photometry to those of other, less extinguished SNe Ib we derive $E(B-V)=0.51\substack{+0.27\\-0.16}$ mag for SN 2019yvr. Correcting photometry of the pre-explosion source for dust reddening, we determine that this source is consistent with a $\log(L/L_{\odot}) = 5.3 \pm 0.2$ and $T_{\mathrm{eff}} = 6800\substack{+400\\-200}$ K star. This relatively cool photospheric temperature implies a radius of 320$\substack{+30\\-50} R_{\odot}$, much larger than expectations for SN Ib progenitor stars with trace amounts of hydrogen but in agreement with previously identified SN IIb progenitor systems. The photometry of the system is also consistent with binary star models that undergo common envelope evolution, leading to a primary star hydrogen envelope mass that is mostly depleted but seemingly in conflict with the SN Ib classification of SN 2019yvr. SN 2019yvr had signatures of strong circumstellar interaction in late-time ($>$150 day) spectra and imaging, and so we consider eruptive mass loss and common envelope evolution scenarios that explain the SN Ib spectroscopic class, pre-explosion counterpart, and dense circumstellar material. We also hypothesize that the apparent inflation could be caused by a quasi-photosphere formed in an extended, low-density envelope or circumstellar matter around the primary star. △ Less

Submitted 8 January, 2021; originally announced January 2021.

Comments: 22 pages, 9 figures, submitted to MNRAS

Searches after Gravitational Waves Using ARizona Observatories (SAGUARO): Observations and Analysis from Advanced LIGO/Virgo's Third Observing Run

Authors: K. Paterson, M. J. Lundquist, J. C. Rastinejad, W. Fong, D. J. Sand, J. E. Andrews, R. C. Amaro, O. Eskandari, S. Wyatt, P. N. Daly, H. Bradley, S. Zhou-Wright, S. Valenti, S. Yang, E. Christensen, A. R. Gibbs, F. Shelly, C. Bilinski, L. Chomiuk, A. Corsi, M. R. Drout, R. J. Foley, P. Gabor, P. Garnavich, C. J. Grier , et al. (12 additional authors not shown)

Abstract: 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$σ$ transient detec… ▽ More 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$σ$ 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. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: 44 pages, 21 figures, 2 tables, submitted to ApJ

The Red Supergiant Content of M31 and M33

Authors: Philip Massey, Kathryn F. Neugent, Emily M. Levesque, Maria R. Drout, Stephane Courteau

Abstract: We identify red supergiants (RSGs) in our spiral neighbors M31 and M33 using near-IR (NIR) photometry complete to a luminosity limit of log L/Lo=4.0. Our archival survey data cover 5 deg^2 of M31, and 3 deg^2 for M33, and are likely spatially complete for these massive stars. Gaia is used to remove foreground stars, after which the RSGs can be separated from asymptotic giant branch (AGB) stars in… ▽ More We identify red supergiants (RSGs) in our spiral neighbors M31 and M33 using near-IR (NIR) photometry complete to a luminosity limit of log L/Lo=4.0. Our archival survey data cover 5 deg^2 of M31, and 3 deg^2 for M33, and are likely spatially complete for these massive stars. Gaia is used to remove foreground stars, after which the RSGs can be separated from asymptotic giant branch (AGB) stars in the color-magnitude diagram. The photometry is used to derive effective temperatures and bolometric luminosities via MARCS stellar atmosphere models. The resulting H-R diagrams show superb agreement with the evolutionary tracks of the Geneva evolutionary group. Our census includes 6400 RSGs in M31 and 2850 RSGs in M33 within their Holmberg radii; by contrast, only a few hundred RSGs are known so far in the Milky Way. Our catalog serves as the basis for a study of the RSG binary frequency being published separately, as well as future studies relating to the evolution of massive stars. Here we use the matches between the NIR-selected RSGs and their optical counterparts to show that the apparent similarity in the reddening of OB stars in M31 and M33 is the result of Malmquist bias; the average extinction in M31 is likely higher than that of M33. As expected, the distribution of RSGs follows that of the spiral arms, while the much older AGB population is more uniformly spread across each galaxy's disk. △ Less

Submitted 7 December, 2020; v1 submitted 26 November, 2020; originally announced November 2020.

Comments: Accepted for publication in the Astronomical Journal. Revision is to fix minor typo in Table 2, where M31 and M33 temperature equations had been interchanged

Collapsar R-Process Yields Can Reproduce [Eu/Fe] Abundance Scatter in Metal-Poor Stars

Authors: Kaley Brauer, Alexander P. Ji, Maria R. Drout, Anna Frebel

Abstract: It is unclear if neutron star mergers can explain the observed r-process abundances of metal-poor stars. Collapsars, defined here as rotating massive stars whose collapse results in a rapidly accreting disk around a black hole that can launch jets, are a promising alternative. We find that we can produce a self-consistent model in which a population of collapsars with stochastic europium yields sy… ▽ More It is unclear if neutron star mergers can explain the observed r-process abundances of metal-poor stars. Collapsars, defined here as rotating massive stars whose collapse results in a rapidly accreting disk around a black hole that can launch jets, are a promising alternative. We find that we can produce a self-consistent model in which a population of collapsars with stochastic europium yields synthesizes all of the r-process material in metal-poor ([Fe/H] < -2.5) stars. Our model reproduces the observed scatter and evolution of scatter of [Eu/Fe] abundances. We find that if collapsars are the dominant r-process site for metal-poor stars, r-process synthesis may be linked to supernovae that produce long gamma-ray bursts. Our results also allow for the possibility that core-collapse supernovae beyond those that launch gamma-ray bursts also produce r-process material (e.g., potentially a subset of Type Ic-BL supernovae). Furthermore, we identify collapsar jet properties (isotropic energy, engine luminosity, or engine time) which may trace r-process yield and verify that the amount of r-process yield produced per collapsar in our model (~0.07 Msun) is consistent with other independent estimates. In the future, achieving 0.05 dex precision on distribution scatter or a reliable selection function would further constrain our probe of r-process production. Our model would also hold for another prompt r-process site with a power-law yield, and work is needed to determine if, for example, fast-merging neutron stars can also explain abundance scatter. △ Less

Submitted 19 May, 2021; v1 submitted 29 October, 2020; originally announced October 2020.

Comments: 17 pages, 8 figures. Accepted by ApJ

The Young Supernova Experiment: Survey Goals, Overview, and Operations

Authors: D. O. Jones, R. J. Foley, G. Narayan, J. Hjorth, M. E. Huber, P. D. Aleo, K. D. Alexander, C. R. Angus, K. Auchettl, V. F. Baldassare, S. H. Bruun, K. C. Chambers, D. Chatterjee, D. L. Coppejans, D. A. Coulter, L. DeMarchi, G. Dimitriadis, M. R. Drout, A. Engel, K. D. French, A. Gagliano, C. Gall, T. Hung, L. Izzo, W. V. Jacobson-Galán , et al. (46 additional authors not shown)

Abstract: 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 tele… ▽ More 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. △ Less

Submitted 5 January, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: ApJ, in press; more information at https://yse.ucsc.edu/

The Nickel Mass Distribution of Stripped-Envelope Supernovae: Implications for Additional Power Sources

Authors: Niloufar Afsariardchi, Maria R. Drout, David Khatami, Christopher D. Matzner, Dae-Sik Moon, Yuan Qi Ni

Abstract: We perform a systematic study of the $^{56}$Ni mass ($M_{\rm Ni}$) of 27 stripped envelope supernovae (SESNe) by modeling their light-curve tails, highlighting that use of ``Arnett's rule'' overestimates $M_{\rm Ni}$ for SESN by a factor of $\sim$2. Recently, \citet{Khatami2019} presented a new model relating the peak time ($t_{\rm p}$) and luminosity ($L_{\rm p}$) of a radioactive-powered SN to i… ▽ More We perform a systematic study of the $^{56}$Ni mass ($M_{\rm Ni}$) of 27 stripped envelope supernovae (SESNe) by modeling their light-curve tails, highlighting that use of ``Arnett's rule'' overestimates $M_{\rm Ni}$ for SESN by a factor of $\sim$2. Recently, \citet{Khatami2019} presented a new model relating the peak time ($t_{\rm p}$) and luminosity ($L_{\rm p}$) of a radioactive-powered SN to its $M_{\rm Ni}$ that addresses several limitations of Arnett-like models, but depends on a dimensionless parameter, $β$. Using observed $t_{\rm p}$, $L_{\rm p}$, and tail-measured $M_{\rm Ni}$ values for 27 SESN, we observationally calibrate $β$ for the first time. Despite scatter, we demonstrate that the model of \citet{Khatami2019} with empirically-calibrated $β$ values provides significantly improved measurements of $M_{\rm Ni}$ when only photospheric data is available. However, these observationally-constrained $β$ values are systematically lower than those inferred from numerical simulations, primarily because the observed sample has significantly higher (0.2-0.4 dex) $L_{\rm p}$ for a given $M_{\rm Ni}$. While effects due to composition, mixing, and asymmetry can increase $L_{\rm p}$ current models cannot explain the systematically low $β$ values. However, the discrepancy can be alleviated if $\sim$7--50\% of $L_{\rm p}$ for the observed sample originates from sources other than $^{56}$Ni. Either shock cooling or magnetar spin-down could provide the requisite luminosity. Finally, we find that even with our improved measurements, the $M_{\rm Ni}$ values of SESN are still a factor of $\sim$3 larger than those of hydrogen-rich Type II SN, indicating that these supernovae are inherently different in terms of their progenitor initial mass distributions or explosion mechanisms. △ Less

Submitted 15 June, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 27 pages, 12 figures, accepted for publication in ApJ

Cool, Luminous, and Highly Variable Stars in the Magellanic Clouds from ASAS-SN: Implications for Thorne-Żytkow Objects and Super-Asymptotic Giant Branch Stars

Authors: Anna J. G. O'Grady, Maria R. Drout, B. J. Shappee, Evan B. Bauer, Jim Fuller, C. S. Kochanek, T. Jayasinghe, B. M. Gaensler, K. Z. Stanek, Thomas W. -S. Holoien, J. L. Prieto, Todd A. Thompson

Abstract: Stars with unusual properties can provide a wealth of information about rare stages of stellar evolution and exotic physics. However, determining the true nature of peculiar stars is often difficult. In this work, we conduct a systematic search for cool and luminous stars in the Magellanic Clouds with extreme variability, motivated by the properties of the unusual SMC star and Thorne-Żytkow Object… ▽ More Stars with unusual properties can provide a wealth of information about rare stages of stellar evolution and exotic physics. However, determining the true nature of peculiar stars is often difficult. In this work, we conduct a systematic search for cool and luminous stars in the Magellanic Clouds with extreme variability, motivated by the properties of the unusual SMC star and Thorne-Żytkow Object (TŻO) candidate HV2112. Using light curves from ASAS-SN we identify 38 stars with surface temperatures T $<$ 4800K, luminosities $\log$(L/L$_\odot$) $>$ 4.3, variability periods $>$ 400 days, and variability amplitudes $Δ$V $>$ 2.5 mag. Eleven of these stars possess the distinctive double-peaked light curve morphology of HV2112. We use the pulsation properties and derived occurrence rates for these 12 objects to constrain their nature. From comparisons to stellar populations and models, we find that one star may be a red supergiant with large amplitude pulsations. For the other 11 stars we derive current masses of $\sim$5-10 M$_{\odot}$, below the theoretical minimum mass of $\sim$15 M$_{\odot}$ for TŻOs to be stable, casting doubt on this interpretation. Instead, we find that the temperatures, luminosities, mass-loss rates, and periods of these stars are consistent with predictions for super-Asymptotic Giant Branch (s-AGB) stars that have begun carbon burning but have not reached the superwind phase. We infer lifetimes in this phase of $\sim($1$-$7) $\times$ 10$^{4}$ years, also consistent with an s-AGB interpretation. If confirmed, these objects would represent the first identified population of s-AGB stars, illuminating the transition between low- and high-mass stellar evolution. △ Less

Submitted 14 August, 2020; originally announced August 2020.

Comments: 29 pages, 14 figures, 4 tables, plus Appendix. Accepted for publication in the Astrophysical Journal

SuperRAENN: A Semi-supervised Supernova Photometric Classification Pipeline Trained on Pan-STARRS1 Medium Deep Survey Supernovae

Authors: V. Ashley Villar, Griffin Hosseinzadeh, Edo Berger, Michelle Ntampaka, David O. Jones, Peter Challis, Ryan Chornock, Maria R. Drout, Ryan J. Foley, Robert P. Kirshner, Ragnhild Lunnan, Raffaella Margutti, Dan Milisavljevic, Nathan Sanders, Yen-Chen Pan, Armin Rest, Daniel M. Scolnic, Eugene Magnier, Nigel Metcalfe, Richard Wainscoat, Christopher Waters

Abstract: Automated classification of supernovae (SNe) based on optical photometric light curve information is essential in the upcoming era of wide-field time domain surveys, such as the Legacy Survey of Space and Time (LSST) conducted by the Rubin Observatory. Photometric classification can enable real-time identification of interesting events for extended multi-wavelength follow-up, as well as archival p… ▽ More Automated classification of supernovae (SNe) based on optical photometric light curve information is essential in the upcoming era of wide-field time domain surveys, such as the Legacy Survey of Space and Time (LSST) conducted by the Rubin Observatory. Photometric classification can enable real-time identification of interesting events for extended multi-wavelength follow-up, as well as archival population studies. Here we present the complete sample of 5,243 "SN-like" light curves (in griz) from the Pan-STARRS1 Medium-Deep Survey (PS1-MDS). The PS1-MDS is similar to the planned LSST Wide-Fast-Deep survey in terms of cadence, filters and depth, making this a useful training set for the community. Using this dataset, we train a novel semi-supervised machine learning algorithm to photometrically classify 2,315 new SN-like light curves with host galaxy spectroscopic redshifts. Our algorithm consists of a random forest supervised classification step and a novel unsupervised step in which we introduce a recurrent autoencoder neural network (RAENN). Our final pipeline, dubbed SuperRAENN, has an accuracy of 87% across five SN classes (Type Ia, Ibc, II, IIn, SLSN-I). We find the highest accuracy rates for Type Ia SNe and SLSNe and the lowest for Type Ibc SNe. Our complete spectroscopically- and photometrically-classified samples break down into: 62.0% Type Ia (1839 objects), 19.8% Type II (553 objects), 4.8% Type IIn (136 objects), 11.7% Type Ibc (291 objects), and 1.6% Type I SLSNe (54 objects). Finally, we discuss how this algorithm can be modified for online LSST data streams. △ Less

Submitted 11 August, 2020; originally announced August 2020.

Comments: Submitted to ApJ; Companion paper to Hosseinzadeh et al.; Tables 1 and 2 available on ashleyvillar.com prior to publication

Photometric Classification of 2315 Pan-STARRS1 Supernovae with Superphot

Authors: Griffin Hosseinzadeh, Frederick Dauphin, V. Ashley Villar, Edo Berger, David O. Jones, Peter Challis, Ryan Chornock, Maria R. Drout, Ryan J. Foley, Robert P. Kirshner, Ragnhild Lunnan, Raffaella Margutti, Dan Milisavljevic, Yen-Chen Pan, Armin Rest, Daniel M. Scolnic, Eugene Magnier, Nigel Metcalfe, Richard Wainscoat, Christopher Waters

Abstract: The classification of supernovae (SNe) and its impact on our understanding of the explosion physics and progenitors have traditionally been based on the presence or absence of certain spectral features. However, current and upcoming wide-field time-domain surveys have increased the transient discovery rate far beyond our capacity to obtain even a single spectrum of each new event. We must therefor… ▽ More The classification of supernovae (SNe) and its impact on our understanding of the explosion physics and progenitors have traditionally been based on the presence or absence of certain spectral features. However, current and upcoming wide-field time-domain surveys have increased the transient discovery rate far beyond our capacity to obtain even a single spectrum of each new event. We must therefore rely heavily on photometric classification, connecting SN light curves back to their spectroscopically defined classes. Here we present Superphot, an open-source Python implementation of the machine-learning classification algorithm of Villar et al., and apply it to 2315 previously unclassified transients from the Pan-STARRS1 Medium Deep Survey for which we obtained spectroscopic host-galaxy redshifts. Our classifier achieves an overall accuracy of 82%, with completenesses and purities of >80% for the best classes (SNe Ia and superluminous SNe). For the worst performing SN class (SNe Ibc), the completeness and purity fall to 37% and 21%, respectively. Our classifier provides 1257 newly classified SNe Ia, 521 SNe II, 298 SNe Ibc, 181 SNe IIn, and 58 SLSNe. These are among the largest uniformly observed samples of SNe available in the literature and will enable a wide range of statistical studies of each class. △ Less

Submitted 16 December, 2020; v1 submitted 11 August, 2020; originally announced August 2020.

Comments: Updated to match published version in ApJ. Companion paper to Villar et al. (2020). Full machine-readable tables available as ancillary files at right

The Red Supergiant Binary Fraction of the Large Magellanic Cloud

Authors: Kathryn F. Neugent, Emily M. Levesque, Philip Massey, Nidia I. Morrell, Maria R. Drout

Abstract: The binary fraction of unevolved massive stars is thought to be 70-100% but there are few observational constraints on the binary fraction of the evolved version of a subset of these stars, the red supergiants (RSGs). Here we identify a complete sample of RSGs in the Large Magellanic Cloud (LMC) using new spectroscopic observations and archival UV, IR and broadband optical photometry. We find 4090… ▽ More The binary fraction of unevolved massive stars is thought to be 70-100% but there are few observational constraints on the binary fraction of the evolved version of a subset of these stars, the red supergiants (RSGs). Here we identify a complete sample of RSGs in the Large Magellanic Cloud (LMC) using new spectroscopic observations and archival UV, IR and broadband optical photometry. We find 4090 RSGs with log L/Lo > 3.5 with 1820 of them having log L/Lo > 4, which we believe is our completeness limit. We additionally spectroscopically confirmed 38 new RSG+B star binaries in the LMC, bringing the total known up to 55. We then estimated the binary fraction using a k-nearest neighbors algorithm that classifies stars as single or binary based on photometry with a spectroscopic sample as a training set. We take into account observational biases such as line-of-sight stars and binaries in eclipse while also calculating model-dependent corrections for RSGs with companions that our observations were not designed to detect. Based on our data, we find an initial result of 13.5 +7.56/-6.67% for RSGs with O or B-type companions. Using the Binary Population and Spectral Synthesis (BPASS) models to correct for unobserved systems, this corresponds to a total RSG binary fraction of 19.5 +7.6/-6.7%. This number is in broad agreement with what we would expect given an initial OB binary distribution of 70%, a predicted merger fraction of 20-30% and a binary interaction fraction of 40-50%. △ Less

Submitted 31 July, 2020; originally announced July 2020.

Comments: ApJ, accepted

The $R$-Process Alliance: Fourth Data Release from the Search for $r$-Process-Enhanced Stars in the Galactic Halo

Authors: Erika M. Holmbeck, Terese T. Hansen, Timothy C. Beers, Vinicius M. Placco, Devin D. Whitten, Kaitlin C. Rasmussen, Ian U. Roederer, Rana Ezzeddine, Charli M. Sakari, Anna Frebel, Maria R. Drout, Joshua D. Simon, Ian B. Thompson, Joss Bland-Hawthorn, Brad K. Gibson, Eva K. Grebel, Georges Kordopatis, Andrea Kunder, Jorge Melendez, Julio F. Navarro, Warren A. Reid, George Seabroke, Matthias Steinmetz, Fred Watson, Rosemary F. G. Wyse

Abstract: This compilation is the fourth data release from the $R$-Process Alliance (RPA) search for $r$-process-enhanced stars, and the second release based on "snapshot" high-resolution ($R \sim 30,000$) spectra collected with the du Pont 2.5m Telescope. In this data release, we propose a new delineation between the $r$-I and $r$-II stellar classes at $\mathrm{[Eu/Fe]} = +0.7$, instead of the empirically… ▽ More This compilation is the fourth data release from the $R$-Process Alliance (RPA) search for $r$-process-enhanced stars, and the second release based on "snapshot" high-resolution ($R \sim 30,000$) spectra collected with the du Pont 2.5m Telescope. In this data release, we propose a new delineation between the $r$-I and $r$-II stellar classes at $\mathrm{[Eu/Fe]} = +0.7$, instead of the empirically chosen $\mathrm{[Eu/Fe]} = +1.0$ level previously in use, based on statistical tests of the complete set of RPA data released to date. We also statistically justify the minimum level of [Eu/Fe] for definition of the $r$-I stars, [Eu/Fe] $> +0.3$. Redefining the separation between $r$-I and $r$-II stars will aid in analysis of the possible progenitors of these two classes of stars and whether these signatures arise from separate astrophysical sources at all. Applying this redefinition to previous RPA data, the number of identified $r$-II and $r$-I stars changes to 51 and 121, respectively, from the initial set of data releases published thus far. In this data release, we identify 21 new $r$-II, 111 new $r$-I (plus three re-identified), and 7 new (plus one re-identified) limited-$r$ stars out of a total of 232 target stars, resulting in a total sample of 72 new $r$-II stars, 232 new $r$-I stars, and 42 new limited-$r$ stars identified by the RPA to date. △ Less

Submitted 1 July, 2020; originally announced July 2020.

Comments: 17 pages, 7 figures, 6 tables; accepted to ApJS

Stars stripped in binaries -- the living gravitational wave sources

Authors: Y. Götberg, V. Korol, A. Lamberts, T. Kupfer, K. Breivik, B. Ludwig, M. R. Drout

Abstract: Binary interaction can cause stellar envelopes to be stripped, which significantly reduces the radius of the star. The orbit of a binary composed of a stripped star and a compact object can therefore be so tight that the gravitational radiation the system produces reaches frequencies accessible to the Laser Interferometer Space Antenna (LISA). Two such stripped stars in tight orbits with white dwa… ▽ More Binary interaction can cause stellar envelopes to be stripped, which significantly reduces the radius of the star. The orbit of a binary composed of a stripped star and a compact object can therefore be so tight that the gravitational radiation the system produces reaches frequencies accessible to the Laser Interferometer Space Antenna (LISA). Two such stripped stars in tight orbits with white dwarfs are known so far (ZTF J2130+4420 and CD-30 11223), but many more are expected to exist. These binaries provide important constraints for binary evolution models and may be used as LISA verification sources. We develop a Monte Carlo code that uses detailed evolutionary models to simulate the Galactic population of stripped stars in tight orbits with either neutron star or white dwarf companions. We predict 0-100 stripped star + white dwarf binaries and 0-4 stripped star + neutron star binaries with SNR>5 after 10 years of observations with LISA. More than 90% of these binaries are expected to show large radial velocity shifts of $\gtrsim 200$ km/s, which are spectroscopically detectable. Photometric variability due to tidal deformation of the stripped star is also expected and has been observed in ZTF J2130+4420 and CD-30 11223. In addition, the stripped star + neutron star binaries are expected to be X-ray bright with $L_X \gtrsim 10^{33} - 10^{36}$ erg/s. Our results show that stripped star binaries are promising multi-messenger sources for the upcoming electromagnetic and gravitational wave facilities. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Comments: Submitted to ApJ

A Demonstration of Extremely Low Latency $γ$-ray, X-Ray & UV Follow-Up of a Millisecond Radio Transient

Authors: Aaron Tohuvavohu, Casey J. Law, Jamie A. Kennea, Elizabeth A. K. Adams, Kshitij Aggarwal, Geoffrey Bower, Sarah Burke-Spolaor, Bryan J. Butler, John M. Cannon, S. Bradley Cenko, James DeLaunay, Paul Demorest, Maria R. Drout, Philip A. Evans, Alec S. Hirschauer, T. J. W. Lazio, Justin Linford, Francis E. Marshall, K. McQuinn, Emily Petroff, Evan D. Skillman

Abstract: We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved aro… ▽ More We report results of a novel high-energy follow-up observation of a potential Fast Radio Burst. The radio burst was detected by VLA/realfast and followed-up by the Neil Gehrels Swift Observatory in very low latency utilizing new operational capabilities of Swift (arXiv:2005.01751), with pointed soft X-ray and UV observations beginning at T0+32 minutes, and hard X-ray/gamma-ray event data saved around T0. These observations are $>10$x faster than previous X-ray/UV follow-up of any radio transient to date. No emission is seen coincident with the FRB candidate at T0, with a 0.2s fluence $5σ$ upper limit of $1.35\times10^{-8}$ erg cm$^{-2}$ (14-195 keV) for a SGR 1935+2154-like flare, nor at T0+32 minutes down to $3σ$ upper limits of 22.18 AB mag in UVOT u band, and $3.33\times10^{-13}$ erg cm$^{-2}$ s$^{-1}$ from 0.3-10 keV for the 2 ks observation. The candidate FRB alone is not significant enough to be considered astrophysical, so this note serves as a technical demonstration. These new Swift operational capabilities will allow future FRB detections to be followed up with Swift at even lower latencies than demonstrated here: 15-20 minutes should be regularly achievable, and 5-10 minutes occasionally achievable. We encourage FRB detecting facilities to release alerts in low latency to enable this science. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: Technical note and capability update for the community. We encourage low latency FRB alerts from relevant facilities to enable this science

A Deep CFHT Optical Search for a Counterpart to the Possible Neutron Star -- Black Hole Merger GW190814

Authors: Nicholas Vieira, John J. Ruan, Daryl Haggard, Maria R. Drout, Melania C. Nynka, Hope Boyce, Kristine Spekkens, Samar Safi-Harb, Raymond G. Carlberg, Rodrigo Fernández, Anthony L. Piro, Niloufar Afsariardchi, Dae-Sik Moon

Abstract: We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawai… ▽ More We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star - black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of a NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope. We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multi-band images of the region between 1.7 and 8.7 days post-merger, reaching a 5sigma depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that either the lighter object was tidally disrupted inside of the BH's innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5sigma source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej < 0.015Msun for a 'blue' (kappa = 0.5 cm^2 g^-1) kilonova, and Mej < 0.04Msun for a 'red' (kappa = 5-10 cm^2 g^-1) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational wave events. △ Less

Submitted 13 May, 2020; v1 submitted 20 March, 2020; originally announced March 2020.

Comments: 23 pages, 13 figures including appendices; accepted to ApJ with minor revisions

Thirty Years of Radio Observations of Type Ia SN 1972E and SN 1895B: Constraints on Circumstellar Shells

Authors: Yvette Cendes, Maria R. Drout, Laura Chomiuk, Sumit K. Sarbadhicary

Abstract: We have imaged over 35 years of archival Very Large Array (VLA) observations of the nearby (d$_{\rm{L}}$ $=$ 3.15 Mpc) Type Ia supernovae SN\,1972E and SN\,1895B between 9 and 121 years post-explosion. No radio emission is detected, constraining the 8.5 GHz luminosities of SN\,1972E and SN\,1895B to be L$_{ν,8.5\rm{GHz}}$ $<$ 6.0 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 45 years post-explosion an… ▽ More We have imaged over 35 years of archival Very Large Array (VLA) observations of the nearby (d$_{\rm{L}}$ $=$ 3.15 Mpc) Type Ia supernovae SN\,1972E and SN\,1895B between 9 and 121 years post-explosion. No radio emission is detected, constraining the 8.5 GHz luminosities of SN\,1972E and SN\,1895B to be L$_{ν,8.5\rm{GHz}}$ $<$ 6.0 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 45 years post-explosion and L$_{ν,8.5\rm{GHz}}$ $<$ 8.9 $\times$ 10$^{23}$ erg s$^{-1}$ Hz$^{-1}$ 121 years post-explosion, respectively. These limits imply a clean circumstellar medium (CSM), with $n$ $<$ 0.9 cm$^{-3}$ out to radii of a few $\times$ 10$^{18}$ cm, if the SN blastwave is expanding into uniform density material. Due to the extensive time coverage of our observations, we also constrain the presence of CSM shells surrounding the progenitor of SN\,1972E. We rule out essentially all medium and thick shells with masses of 0.05$-$0.3 M$_\odot$ at radii between $\sim$10$^{17}$ and 10$^{18}$ cm, and thin shells at specific radii with masses down to $\lesssim$0.01 M$_\odot$. These constraints rule out swaths of parameter space for a range of single and double degenerate progenitor scenarios, including recurrent nova, core-degenerate objects, ultra-prompt explosions and white dwarf (WD) mergers with delays of a few hundred years between the onset of merger and explosion. Allowed progenitors include WD-WD systems with a significant ($>$ 10$^{4}$ years) delay from the last episode of common envelope evolution and single degenerate systems undergoing recurrent nova, provided that the recurrence timescale i short and the system has been in the nova phase for $\gtrsim$10$^{4}$ yr, such that a large ($>$ 10$^{18}$ cm) cavity has been evacuated. Future multi-epoch observations of additional intermediate-aged Type Ia SNe will provide a comprehensive view of the large-scale CSM environments around these explosions. △ Less

Submitted 10 January, 2020; originally announced January 2020.

Comments: Accepted for publication in the Astrophysical Journal

The Luminosity Function of Red Supergiants in M31

Authors: Kathryn F. Neugent, Philip Massey, Cyril Georgy, Maria R. Drout, Michael Mommert, Emily M. Levesque, Georges Meynet, Sylvia Ekstrom

Abstract: The mass-loss rates of red supergiant stars (RSGs) are poorly constrained by direct measurements, and yet the subsequent evolution of these stars depends critically on how much mass is lost during the RSG phase. In 2012 the Geneva evolutionary group updated their mass-loss prescription for RSGs with the result that a 20 solar mass star now loses 10x more mass during the RSG phase than in the older… ▽ More The mass-loss rates of red supergiant stars (RSGs) are poorly constrained by direct measurements, and yet the subsequent evolution of these stars depends critically on how much mass is lost during the RSG phase. In 2012 the Geneva evolutionary group updated their mass-loss prescription for RSGs with the result that a 20 solar mass star now loses 10x more mass during the RSG phase than in the older models. Thus, higher mass RSGs evolve back through a second yellow supergiant phase rather than exploding as Type II-P supernovae, in accord with recent observations (the so-called "RSG Problem"). Still, even much larger mass-loss rates during the RSG phase cannot be ruled out by direct measurements of their current dust-production rates, as such mass-loss is episodic. Here we test the models by deriving a luminosity function for RSGs in the nearby spiral galaxy M31 which is sensitive to the total mass loss during the RSG phase. We carefully separate RSGs from asymptotic giant branch stars in the color-magnitude diagram following the recent method exploited by Yang and collaborators in their Small Magellanic Cloud studies. Comparing our resulting luminosity function to that predicted by the evolutionary models shows that the new prescription for RSG mass-loss does an excellent job of matching the observations, and we can readily rule out significantly larger values. △ Less

Submitted 24 November, 2019; originally announced November 2019.

Comments: ApJ, accepted