Showing 1–50 of 92 results for author: Villar, V A

The Type I Superluminous Supernova Catalog I: Light Curve Properties, Models, and Catalog Description

Authors: Sebastian Gomez, Matt Nicholl, Edo Berger, Peter K. Blanchard, V. Ashley Villar, Sofia Rest, Griffin Hosseinzadeh, Aysha Aamer, Yukta Ajay, Wasundara Athukoralalage, David C. Coulter, Tarraneh Eftekhari, Achille Fiore, Noah Franz, Ori Fox, Alexander Gagliano, Daichi Hiramatsu, D. Andrew Howell, Brian Hsu, Mitchell Karmen, Matthew R. Siebert, Réka Könyves-Tóth, Harsh Kumar, Curtis McCully, Craig Pellegrino , et al. (3 additional authors not shown)

Abstract: We present the most comprehensive catalog to date of Type I Superluminous Supernovae (SLSNe), a class of stripped envelope supernovae (SNe) characterized by exceptionally high luminosities. We have compiled a sample of 262 SLSNe reported through 2022 December 31. We verified the spectroscopic classification of each SLSN and collated an exhaustive data set of UV, optical and IR photometry from both… ▽ More We present the most comprehensive catalog to date of Type I Superluminous Supernovae (SLSNe), a class of stripped envelope supernovae (SNe) characterized by exceptionally high luminosities. We have compiled a sample of 262 SLSNe reported through 2022 December 31. We verified the spectroscopic classification of each SLSN and collated an exhaustive data set of UV, optical and IR photometry from both publicly available data and our own FLEET observational follow-up program, totaling over 30,000 photometric detections. Using these data we derive observational parameters such as the peak absolute magnitudes, rise and decline timescales, as well as bolometric luminosities, temperature and photospheric radius evolution for all SLSNe. Additionally, we model all light curves using a hybrid model that includes contributions from both a magnetar central engine and the radioactive decay of $^{56}$Ni. We explore correlations among various physical and observational parameters, and recover the previously found relation between ejecta mass and magnetar spin, as well as the overall progenitor pre-explosion mass distribution with a peak at $\approx 6.5$ M$_\odot$. We find no significant redshift dependence for any parameter, and no evidence for distinct sub-types of SLSNe. We find that $< 3$\% of SLSNe are best fit with a significant contribution from radioactive decay $\gtrsim 50$\%, representing a set of relatively dim and slowly declining SNe. We provide several analytical tools designed to simulate typical SLSN light curves across a broad range of wavelengths and phases, enabling accurate K-corrections, bolometric scaling calculations, and inclusion of SLSNe in survey simulations or future comparison works. The complete catalog, including all of the photometry, models, and derived parameters, is made available as an open-source resource on GitHub. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 59 pages, 22 Figures, Submitted to MNRAS

Anomaly Detection and Approximate Similarity Searches of Transients in Real-time Data Streams

Authors: P. D. Aleo, A. W. Engel, G. Narayan, C. R. Angus, K. Malanchev, K. Auchettl, V. F. Baldassare, A. Berres, T. J. L. de Boer, B. M. Boyd, K. C. Chambers, K. W. Davis, N. Esquivel, D. Farias, R. J. Foley, A. Gagliano, C. Gall, H. Gao, S. Gomez, M. Grayling, C. -C. Lin, E. A. Magnier, K. S. Mandel, T. Matheson, S. I. Raimundo , et al. (5 additional authors not shown)

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

Submitted 1 April, 2024; originally announced April 2024.

Comments: 43 pages (67 pages with Appendix), 15 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

Kilonova Light-Curve Interpolation with Neural Networks

Authors: Yinglei Peng, Marko Ristić, Atul Kedia, Richard O'Shaughnessy, Christopher J. Fontes, Chris L. Fryer, Oleg Korobkin, Matthew R. Mumpower, V. Ashley Villar, Ryan T. Wollaeger

Abstract: Kilonovae are the electromagnetic transients created by the radioactive decay of freshly synthesized elements in the environment surrounding a neutron star merger. To study the fundamental physics in these complex environments, kilonova modeling requires, in part, the use of radiative transfer simulations. The microphysics involved in these simulations results in high computational cost, prompting… ▽ More Kilonovae are the electromagnetic transients created by the radioactive decay of freshly synthesized elements in the environment surrounding a neutron star merger. To study the fundamental physics in these complex environments, kilonova modeling requires, in part, the use of radiative transfer simulations. The microphysics involved in these simulations results in high computational cost, prompting the use of emulators for parameter inference applications. Utilizing a training set of 22248 high-fidelity simulations, we use a neural network to efficiently train on existing radiative transfer simulations and predict light curves for new parameters in a fast and computationally efficient manner. Our neural network can generate millions of new light curves in under a minute. We discuss our emulator's degree of off-sample reliability and parameter inference of the AT2017gfo observational data. Finally, we discuss tension introduced by multi-band inference in the parameter inference results, particularly with regard to the neural network's recovery of viewing angle. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Report number: LA-UR-24-20638

A Physics-Informed Variational Autoencoder for Rapid Galaxy Inference and Anomaly Detection

Authors: Alexander Gagliano, V. Ashley Villar

Abstract: The Vera C. Rubin Observatory is slated to observe nearly 20 billion galaxies during its decade-long Legacy Survey of Space and Time. The rich imaging data it collects will be an invaluable resource for probing galaxy evolution across cosmic time, characterizing the host galaxies of transient phenomena, and identifying novel populations of anomalous systems. To facilitate these studies, we introdu… ▽ More The Vera C. Rubin Observatory is slated to observe nearly 20 billion galaxies during its decade-long Legacy Survey of Space and Time. The rich imaging data it collects will be an invaluable resource for probing galaxy evolution across cosmic time, characterizing the host galaxies of transient phenomena, and identifying novel populations of anomalous systems. To facilitate these studies, we introduce a convolutional variational autoencoder trained to estimate the redshift, stellar mass, and star-formation rates of galaxies from multi-band imaging data. We train and test our physics-informed CVAE on a spectroscopic sample of $\sim$26,000 galaxies within $z<1$ imaged through the Dark Energy Camera Legacy Survey. We show that our model can infer redshift and stellar mass more accurately than the latest image-based self-supervised learning approaches, and is >100x faster than more computationally-intensive SED-fitting techniques. Using a small sample of Green Pea and Red Spiral galaxies reported in the literature, we further demonstrate how this CVAE can be used to rapidly identify rare galaxy populations and interpret what makes them unique. △ Less

Submitted 27 December, 2023; originally announced December 2023.

Comments: Accepted to the "Machine Learning and the Physical Sciences" Workshop at NeurIPS 2023

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

Hierarchical Cross-entropy Loss for Classification of Astrophysical Transients

Authors: V. Ashley Villar, Kaylee de Soto, Alex Gagliano

Abstract: Astrophysical transient phenomena are traditionally classified spectroscopically in a hierarchical taxonomy; however, this graph structure is currently not utilized in neural net-based photometric classifiers for time-domain astrophysics. Instead, independent classifiers are trained for different tiers of classified data, and events are excluded if they fall outside of these well-defined but flat… ▽ More Astrophysical transient phenomena are traditionally classified spectroscopically in a hierarchical taxonomy; however, this graph structure is currently not utilized in neural net-based photometric classifiers for time-domain astrophysics. Instead, independent classifiers are trained for different tiers of classified data, and events are excluded if they fall outside of these well-defined but flat classification schemes. Here, we introduce a weighted hierarchical cross-entropy objective function for classification of astrophysical transients. Our method allows users to directly build and use physics- or observationally-motivated tree-based taxonomies. Our weighted hierarchical cross-entropy loss directly uses this graph to accurately classify all targets into any node of the tree, re-weighting imbalanced classes. We test our novel loss on a set of variable stars and extragalactic transients from the Zwicky Transient Facility, showing that we can achieve similar performance to fine-tuned classifiers with the advantage of notably more flexibility in downstream classification tasks. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: Accepted to Machine Learning and the Physical Sciences Workshop at NeurIPS 2023

Constructing Impactful Machine Learning Research for Astronomy: Best Practices for Researchers and Reviewers

Authors: D. Huppenkothen, M. Ntampaka, M. Ho, M. Fouesneau, B. Nord, J. E. G. Peek, M. Walmsley, J. F. Wu, C. Avestruz, T. Buck, M. Brescia, D. P. Finkbeiner, A. D. Goulding, T. Kacprzak, P. Melchior, M. Pasquato, N. Ramachandra, Y. -S. Ting, G. van de Ven, S. Villar, V. A. Villar, E. Zinger

Abstract: Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best pr… ▽ More Machine learning has rapidly become a tool of choice for the astronomical community. It is being applied across a wide range of wavelengths and problems, from the classification of transients to neural network emulators of cosmological simulations, and is shifting paradigms about how we generate and report scientific results. At the same time, this class of method comes with its own set of best practices, challenges, and drawbacks, which, at present, are often reported on incompletely in the astrophysical literature. With this paper, we aim to provide a primer to the astronomical community, including authors, reviewers, and editors, on how to implement machine learning models and report their results in a way that ensures the accuracy of the results, reproducibility of the findings, and usefulness of the method. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 14 pages, 3 figures; submitted to the Bulletin of the American Astronomical Society

JWST Observations of the Extraordinary GRB 221009A Reveal an Ordinary Supernova Without Signs of $r$-Process Enrichment in a Low-Metallicity Galaxy

Authors: Peter K. Blanchard, V. Ashley Villar, Ryan Chornock, Tanmoy Laskar, Yijia Li, Joel Leja, Justin Pierel, Edo Berger, Raffaella Margutti, Kate D. Alexander, Jennifer Barnes, Yvette Cendes, Tarraneh Eftekhari, Daniel Kasen, Natalie LeBaron, Brian D. Metzger, James Muzerolle Page, Armin Rest, Huei Sears, Daniel M. Siegel, S. Karthik Yadavalli

Abstract: Identifying the astrophysical sites of the $r$-process, one of the primary mechanisms by which heavy elements are formed, is a key goal of modern astrophysics. The discovery of the brightest gamma-ray burst of all time, GRB 221009A, at a relatively nearby redshift, presented the first opportunity to spectroscopically test the idea that $r$-process elements are produced following the collapse of ra… ▽ More Identifying the astrophysical sites of the $r$-process, one of the primary mechanisms by which heavy elements are formed, is a key goal of modern astrophysics. The discovery of the brightest gamma-ray burst of all time, GRB 221009A, at a relatively nearby redshift, presented the first opportunity to spectroscopically test the idea that $r$-process elements are produced following the collapse of rapidly rotating massive stars. Here we present spectroscopic and photometric $\textit{James Webb Space Telescope}$ (JWST) observations of GRB 221009A obtained $+168$ and $+170$ rest-frame days after the initial gamma-ray trigger, and demonstrate they are well-described by a supernova (SN) and power-law afterglow, with no evidence for an additional component from $r$-process emission, and that the SN component strongly resembles the near-infrared spectra of previous SNe, including SN 1998bw. We further find that the SN associated with GRB 221009A is slightly fainter than the expected brightness of SN 1998bw at this phase, concluding that the SN is therefore not an unusual GRB-SN. We infer a nickel mass of $\approx0.09$ M$_{\odot}$, consistent with the lack of an obvious SN detection in the early-time data. We find that the host galaxy of GRB 221009A has a very low metallicity of $\approx0.12$ Z$_{\odot}$ and our resolved host spectrum shows that GRB 221009A occurred in a unique environment in its host characterized by strong H$_2$ emission lines consistent with recent star formation, which may hint at environmental factors being responsible for its extreme energetics. △ Less

Submitted 27 August, 2023; originally announced August 2023.

Comments: 32 pages, 14 Figures, Submitted to Nature Astronomy

SN 2022oqm: A Bright and Multi-peaked Calcium-rich Transient

Authors: S. Karthik Yadavalli, V. Ashley Villar, Luca Izzo, Yossef Zenati, Ryan J. Foley, J. Craig Wheeler, Charlotte R. Angus, Dominik Bánhidi, Katie Auchettl, Barna Imre Bíró, Attila Bódi, Zsófia Bodola, Thomas de Boer, Kenneth C. Chambers, Ryan Chornock, David A. Coulter, István Csányi, Borbála Cseh, Srujan Dandu, Kyle W. Davis, Connor Braden Dickinson, Diego Farias, Joseph Farah, Christa Gall, Hua Gao , et al. (38 additional authors not shown)

Abstract: We present the photometric and spectroscopic evolution of SN 2022oqm, a nearby multi-peaked hydrogen- and helium-weak calcium-rich transient (CaRT). SN 2022oqm was detected 13.1 kpc from its host galaxy, the face-on spiral galaxy NGC 5875. Extensive spectroscopic coverage reveals an early hot (T >= 40,000 K) continuum and carbon features observed $\sim$1~day after discovery, SN Ic-like photospheri… ▽ More We present the photometric and spectroscopic evolution of SN 2022oqm, a nearby multi-peaked hydrogen- and helium-weak calcium-rich transient (CaRT). SN 2022oqm was detected 13.1 kpc from its host galaxy, the face-on spiral galaxy NGC 5875. Extensive spectroscopic coverage reveals an early hot (T >= 40,000 K) continuum and carbon features observed $\sim$1~day after discovery, SN Ic-like photospheric-phase spectra, and strong forbidden calcium emission starting 38 days after discovery. SN 2022oqm has a relatively high peak luminosity (MB = -17 mag) for (CaRTs), making it an outlier in the population. We determine that three power sources are necessary to explain the light curve (LC), with each corresponding to a distinct peak. The first peak is powered by an expanding blackbody with a power law luminosity, suggesting shock cooling by circumstellar material (CSM). Subsequent LC evolution is powered by a double radioactive decay model, consistent with two sources of photons diffusing through optically thick ejecta. From the LC, we derive an ejecta mass and 56Ni mass of ~0.6 solar masses and ~0.09 solar masses. Spectroscopic modeling suggests 0.6 solar masses of ejecta, and with well-mixed Fe-peak elements throughout. We discuss several physical origins for SN 2022oqm and find either a surprisingly massive white dwarf progenitor or a peculiar stripped envelope model could explain SN 2022oqm. A stripped envelope explosion inside a dense, hydrogen- and helium-poor CSM, akin to SNe Icn, but with a large 56Ni mass and small CSM mass could explain SN 2022oqm. Alternatively, helium detonation on an unexpectedly massive white dwarf could also explain SN 2022oqm. △ Less

Submitted 4 April, 2024; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 35 pages, 17 figures, 7 tables, Accepted for Publication in ApJ

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

Roman CCS White Paper: Characterizing Superluminous Supernovae with Roman

Authors: Sebastian Gomez, Kate Alexander, Edo Berger, Peter K. Blanchard, Floor Broekgaarden, Tarraneh Eftekhari, Ori Fox, Kiranjyot Gill, Daichi Hiramatsu, Bhavin Joshi, Mitchell Karmen, Takashi Moriya, Matt Nicholl, Robert Quimby, Eniko Regos, Armin Rest, Benjamin Rose, Melissa Shahbandeh, V. Ashley Villar

Abstract: Type-I Superluminous Supernovae (SLSNe) are an exotic class of core-collapse SN (CCSN) that can be up to 100 times brighter and more slowly-evolving than normal CCSNe. SLSNe represent the end-stages of the most massive stripped stars, and are thought to be powered by the spin-down energy of a millisecond magnetar. Studying them and measuring their physical parameters can help us to better understa… ▽ More Type-I Superluminous Supernovae (SLSNe) are an exotic class of core-collapse SN (CCSN) that can be up to 100 times brighter and more slowly-evolving than normal CCSNe. SLSNe represent the end-stages of the most massive stripped stars, and are thought to be powered by the spin-down energy of a millisecond magnetar. Studying them and measuring their physical parameters can help us to better understand stellar mass-loss, evolution, and explosions. Moreover, thanks to their high luminosities, SLSNe can be seen up to greater distances, allowing us to explore how stellar physics evolves as a function of redshift. The High Latitude Time Domain Survey (HLTDS) will provide us with an exquisite dataset that will discover 100s of SLSNe. Here, we focus on the question of which sets of filters and cadences will allow us to best characterize the physical parameters of these SLSNe. We simulate a set of SLSNe at redshifts ranging from z = 0.1 to z = 5.0, using six different sets of filters, and cadences ranging from 5 to 100 days. We then fit these simulated light curves to attempt to recover the input parameter values for their ejecta mass, ejecta velocity, magnetic field strength, and magnetar spin period. We find that four filters are sufficient to accurately characterize SLSNe at redshifts below $z = 3$, and that cadences faster than 20 days are required to obtain measurements with an uncertainty below 10\%, although a cadence of 70 days is still acceptable under certain conditions. Finally, we find that the nominal survey strategy will not be able to properly characterize the most distant SLSNe at $z = 5$. We find that the addition of 60-day cadence observations for 4 years to the nominal HLTDS survey can greatly improve the prospect of characterizing these most extreme and distant SNe, with only an 8\% increase to the time commitment of the survey. △ Less

Submitted 29 June, 2023; originally announced June 2023.

Comments: Submitted as part of the Nancy Grace Roman Space Telescope's Core Community Surveys call for white papers

The LIGO HET Response (LIGHETR) Project to Discover and Spectroscopically Follow Optical Transients Associated with Neutron Star Mergers

Authors: M. J. Bustamante-Rosell, Greg Zeimann, J. Craig Wheeler, Karl Gebhardt, Aaron Zimmerman, Chris Fryer, Oleg Korobkin, Richard Matzner, V. Ashley Villar, S. Karthik Yadavalli, Kaylee M. de Soto, Matthew Shetrone, Steven Janowiecki, Pawan Kumar, David Pooley, Benjamin P. Thomas, Hsin-Yu Chen, Lifan Wang, Jozsef Vinko, David J. Sand, Ryan Wollaeger, Frederic V. Hessman, Kristen B. McQuinn

Abstract: The LIGO HET Response (LIGHETR) project is an enterprise to follow up optical transients (OT) discovered as gravitational wave merger sources by the LIGO/Virgo collaboration (LVC). Early spectroscopy has the potential to constrain crucial parameters such as the aspect angle. The LIGHETR collaboration also includes the capacity to model the spectroscopic evolution of mergers to facilitate a real-ti… ▽ More The LIGO HET Response (LIGHETR) project is an enterprise to follow up optical transients (OT) discovered as gravitational wave merger sources by the LIGO/Virgo collaboration (LVC). Early spectroscopy has the potential to constrain crucial parameters such as the aspect angle. The LIGHETR collaboration also includes the capacity to model the spectroscopic evolution of mergers to facilitate a real-time direct comparison of models with our data. The principal facility is the Hobby-Eberly Telescope. LIGHETR uses the massively-replicated VIRUS array of spectrographs to search for associated OTs and obtain early blue spectra and in a complementary role, the low-resolution LRS-2 spectrograph is used to obtain spectra of viable candidates as well as a densely-sampled series of spectra of true counterparts. Once an OT is identified, the anticipated cadence of spectra would match or considerably exceed anything achieved for GW170817 = AT2017gfo for which there were no spectra in the first 12 hours and thereafter only roughly once daily. We describe special HET-specific software written to facilitate the program and attempts to determine the flux limits to undetected sources. We also describe our campaign to follow up OT candidates during the third observational campaign of the LIGO and Virgo Scientific Collaborations. We obtained VIRUS spectroscopy of candidate galaxy hosts for 5 LVC gravitational wave events and LRS-2 spectra of one candidate for the OT associated with S190901ap. We identified that candidate, ZTF19abvionh = AT2019pip, as a possible Wolf-Rayet star in an otherwise unrecognized nearby dwarf galaxy. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: 26 pages, 15 figures

Applications of Deep Learning to physics workflows

Authors: Manan Agarwal, Jay Alameda, Jeroen Audenaert, Will Benoit, Damon Beveridge, Meghna Bhattacharya, Chayan Chatterjee, Deep Chatterjee, Andy Chen, Muhammed Saleem Cholayil, Chia-Jui Chou, Sunil Choudhary, Michael Coughlin, Maximilian Dax, Aman Desai, Andrea Di Luca, Javier Mauricio Duarte, Steven Farrell, Yongbin Feng, Pooyan Goodarzi, Ekaterina Govorkova, Matthew Graham, Jonathan Guiang, Alec Gunny, Weichangfeng Guo , et al. (43 additional authors not shown)

Abstract: Modern large-scale physics experiments create datasets with sizes and streaming rates that can exceed those from industry leaders such as Google Cloud and Netflix. Fully processing these datasets requires both sufficient compute power and efficient workflows. Recent advances in Machine Learning (ML) and Artificial Intelligence (AI) can either improve or replace existing domain-specific algorithms… ▽ More Modern large-scale physics experiments create datasets with sizes and streaming rates that can exceed those from industry leaders such as Google Cloud and Netflix. Fully processing these datasets requires both sufficient compute power and efficient workflows. Recent advances in Machine Learning (ML) and Artificial Intelligence (AI) can either improve or replace existing domain-specific algorithms to increase workflow efficiency. Not only can these algorithms improve the physics performance of current algorithms, but they can often be executed more quickly, especially when run on coprocessors such as GPUs or FPGAs. In the winter of 2023, MIT hosted the Accelerating Physics with ML at MIT workshop, which brought together researchers from gravitational-wave physics, multi-messenger astrophysics, and particle physics to discuss and share current efforts to integrate ML tools into their workflows. The following white paper highlights examples of algorithms and computing frameworks discussed during this workshop and summarizes the expected computing needs for the immediate future of the involved fields. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: Whitepaper resulting from Accelerating Physics with ML@MIT workshop in Jan/Feb 2023

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

SN 2023ixf in Messier 101: Photo-ionization of Dense, Close-in Circumstellar Material in a Nearby Type II Supernova

Authors: W. V. Jacobson-Galan, L. Dessart, R. Margutti, R. Chornock, R. J. Foley, C. D. Kilpatrick, D. O. Jones, K. Taggart, C. R. Angus, S. Bhattacharjee, L. A. Braff, D. Brethauer, A. J. Burgasser, F. Cao, C. M. Carlile, K. C. Chambers, D. A. Coulter, E. Dominguez-Ruiz, C. B. Dickinson, T. de Boer, A. Gagliano, C. Gall, H. Gao, E. L. Gates, S. Gomez , et al. (43 additional authors not shown)

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

Submitted 21 August, 2023; v1 submitted 7 June, 2023; originally announced June 2023.

Comments: 19 pages, 8 figures. Accepted for publication in ApJL

Multiple Peaks and a Long Precursor in the Type IIn Supernova 2021qqp: An Energetic Explosion in a Complex Circumstellar Environment

Authors: Daichi Hiramatsu, Tatsuya Matsumoto, Edo Berger, Conor Ransome, V. Ashley Villar, Sebastian Gomez, Yvette Cendes, Kishalay De, K. Azalee Bostroem, Joseph Farah, D. Andrew Howell, Curtis McCully, Megan Newsome, Estefania Padilla Gonzalez, Craig Pellegrino, Akihiro Suzuki, Giacomo Terreran

Abstract: We present optical photometry and spectroscopy of the Type IIn supernova (SN) 2021qqp. Its unusual light curve is marked by a long precursor for $\approx300$ days, a rapid increase in brightness for $\approx60$ days, and then a sharp increase of $\approx1.6$ mag in only a few days to a first peak of $M_r \approx -19.5$ mag. The light curve then declines rapidly until it re-brightens to a second di… ▽ More We present optical photometry and spectroscopy of the Type IIn supernova (SN) 2021qqp. Its unusual light curve is marked by a long precursor for $\approx300$ days, a rapid increase in brightness for $\approx60$ days, and then a sharp increase of $\approx1.6$ mag in only a few days to a first peak of $M_r \approx -19.5$ mag. The light curve then declines rapidly until it re-brightens to a second distinct peak of $M_r \approx -17.3$ mag centered at $\approx335$ days after the first peak. The spectra are dominated by Balmer lines with a complex morphology, including a narrow component with a width of $\approx 1300$ km s$^{-1}$ (first peak) and $\approx 2500$ km s$^{-1}$ (second peak) that we associate with the circumstellar medium (CSM) and a P Cygni component with an absorption velocity of $\approx 8500$ km s$^{-1}$ (first peak) and $\approx 5600$ km s$^{-1}$ (second peak) that we associate with the SN-CSM interaction shell. Using the luminosity and velocity evolution, we construct a flexible analytical model, finding two significant mass-loss episodes with peak mass loss rates of $\approx 10$ and $\approx 5\,M_{\odot}$ yr$^{-1}$ about $0.8$ and $2$ yr before explosion, respectively, with a total CSM mass of $\approx 2-4\,M_{\odot}$. We show that the most recent mass-loss episode could explain the precursor for the year preceding the explosion. The SN ejecta mass is constrained to be $\approx 5-30\,M_{\odot}$ for an explosion energy of $\approx (3-10)\times10^{51}$ erg. We discuss eruptive massive stars (luminous blue variable, pulsational pair instability) and an extreme stellar merger with a compact object as possible progenitor channels. △ Less

Submitted 2 April, 2024; v1 submitted 18 May, 2023; originally announced May 2023.

Comments: Updated to match the published article in ApJ, 2024 March 29

Journal ref: The Astrophysical Journal, 964:181 (16pp), 2024 April 1

A Low-Mass Helium Star Progenitor Model for the Type Ibn SN 2020nxt

Authors: Qinan Wang, Anika Goel, Luc Dessart, Ori D. Fox, Melissa Shahbandeh, Sofia Rest, Armin Rest, Jose H. Groh, Andrew Allan, Claes Fransson, Nathan Smith, Griffin Hosseinzadeh, Alexei V. Filippenko, Jennifer Andrews, K. Azalee Bostroem, Thomas G. Brink, Peter Brown, Jamison Burke, Roger Chevalier, Geoffrey C. Clayton, Mi Dai, Kyle W. Davis, Ryan J. Foley, Sebastian Gomez, Chelsea Harris , et al. (33 additional authors not shown)

Abstract: A growing number of supernovae (SNe) are now known to exhibit evidence for significant interaction with a dense, pre-existing, circumstellar medium (CSM). SNe Ibn comprise one such class that can be characterised by both rapidly evolving light curves and persistent narrow He I lines. The origin of such a dense CSM in these systems remains a pressing question, specifically concerning the progenitor… ▽ More A growing number of supernovae (SNe) are now known to exhibit evidence for significant interaction with a dense, pre-existing, circumstellar medium (CSM). SNe Ibn comprise one such class that can be characterised by both rapidly evolving light curves and persistent narrow He I lines. The origin of such a dense CSM in these systems remains a pressing question, specifically concerning the progenitor system and mass-loss mechanism. In this paper, we present multi-wavelength data of the Type Ibn SN 2020nxt, including $HST$/STIS ultraviolet spectra. We fit the data with recently updated CMFGEN models designed to handle configurations for SNe Ibn. The UV coverage yields strong constraints on the energetics and, when combined with the CMFGEN models, offer new insight on potential progenitor systems. We find the most successful model is a $\lesssim4 {\rm M}_\odot$ helium star that lost its $\sim 1\,{\rm M}_\odot$ He-rich envelope in the years preceding core collapse. We also consider viable alternatives, such as a He white dwarf merger. Ultimately, we conclude at least some SNe Ibn do not arise from single, massive ($>30 {\rm M}_\odot$) Wolf-Rayet-like stars. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: 17 pages, 13 figures, 1 table, submitted to MNRAS

Flight of the Bumblebee: the Early Excess Flux of Type Ia Supernova 2023bee revealed by $TESS$, $Swift$ and Young Supernova Experiment Observations

Authors: Qinan Wang, Armin Rest, Georgios Dimitriadis, Ryan Ridden-harper, Matthew R. Siebert, Mark Magee, Charlotte R. Angus, Katie Auchettl, Kyle W. Davis, Ryan J. Foley, Ori D. Fox, Sebastian Gomez, Jacob E. Jencson, David O. Jones, Charles D. Kilpatrick, Justin D. R. Pierel, Anthony L. Piro, Abigail Polin, Collin A. Politsch, César Rojas-bravo, Melissa Shahbandeh, V. Ashley Villar, Yossef Zenati, C. Ashall, Kenneth C. Chambers , et al. (19 additional authors not shown)

Abstract: We present high-cadence ultraviolet through near-infrared observations of the Type Ia supernova (SN Ia) 2023bee in NGC~2708 ($D = 32 \pm 3$ Mpc), finding excess flux in the first days after explosion relative to the expected power-law rise from an expanding fireball. This deviation from typical behavior for SNe Ia is particularly obvious in our 10-minute cadence $TESS$ light curve and $Swift$ UV d… ▽ More We present high-cadence ultraviolet through near-infrared observations of the Type Ia supernova (SN Ia) 2023bee in NGC~2708 ($D = 32 \pm 3$ Mpc), finding excess flux in the first days after explosion relative to the expected power-law rise from an expanding fireball. This deviation from typical behavior for SNe Ia is particularly obvious in our 10-minute cadence $TESS$ light curve and $Swift$ UV data. Compared to a few other normal SNe Ia with detected early excess flux, the excess flux in SN 2023bee is redder in the UV and less luminous. We present optical spectra of SN 2023bee, including two spectra during the period where the flux excess is dominant. At this time, the spectra are similar to those of other SNe Ia but with weaker Si II, C II and Ca II absorption lines, perhaps because the excess flux creates a stronger continuum. We compare the data to several theoretical models that have been proposed to explain the early flux excess in SNe Ia. Interaction with either a nearby companion star or close-in circumstellar material is expected to produce a faster evolution than seen in the data. Radioactive material in the outer layers of the ejecta, either from a double detonation explosion or simply an explosion with a $^{56}$Ni clump near the surface, can not fully reproduce the evolution either, likely due to the sensitivity of early UV observable to the treatment of the outer part of ejecta in simulation. We conclude that no current model can adequately explain the full set of observations. We find that a relatively large fraction of nearby, bright SNe Ia with high-cadence observations have some amount of excess flux within a few days of explosion. Considering potential asymmetric emission, the physical cause of this excess flux may be ubiquitous in normal SNe Ia. △ Less

Submitted 19 November, 2023; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: 21 pages, 12 figures. Accepted by the astrophysical journal

Interpolated kilonova spectra models: necessity for a phenomenological, blue component in the fitting of AT2017gfo spectra

Authors: Marko Ristic, Richard O'Shaughnessy, V. Ashley Villar, Ryan T. Wollaeger, Oleg Korobkin, Chris L. Fryer, Christopher J. Fontes, Atul Kedia

Abstract: In this work, we present a simple interpolation methodology for spectroscopic time series, based on conventional interpolation techniques (random forests) implemented in widely-available libraries. We demonstrate that our existing library of simulations is sufficient for training, producing interpolated spectra that respond sensitively to varied ejecta parameter, post-merger time, and viewing angl… ▽ More In this work, we present a simple interpolation methodology for spectroscopic time series, based on conventional interpolation techniques (random forests) implemented in widely-available libraries. We demonstrate that our existing library of simulations is sufficient for training, producing interpolated spectra that respond sensitively to varied ejecta parameter, post-merger time, and viewing angle inputs. We compare our interpolated spectra to the AT2017gfo spectral data, and find parameters similar to our previous inferences using broadband light curves. However, the spectral observations have significant systematic short-wavelength residuals relative to our models, which we cannot explain within our existing framework. Similar to previous studies, we argue that an additional blue component is required. We consider a radioactive heating source as a third component characterized by light, slow-moving, lanthanide-free ejecta with $M_{\rm th} = 0.003~M_\odot$, $v_{\rm th} = 0.05$c, and $κ_{\rm th} = 1$ cm$^2$/g. When included as part of our radiative transfer simulations, our choice of third component reprocesses blue photons into lower energies, having the opposite effect and further accentuating the blue-underluminosity disparity in our simulations. As such, we are unable to overcome short-wavelength deficits at later times using an additional radioactive heating component, indicating the need for a more sophisticated modeling treatment. △ Less

Submitted 4 June, 2023; v1 submitted 13 April, 2023; originally announced April 2023.

Comments: 11 pages, 7 figures, submitted to Phys. Rev. Research, presented at April APS session F13.00006

Journal ref: Phys. Rev. Res. 5, 043106 (2023)

SBI++: Flexible, Ultra-fast Likelihood-free Inference Customized for Astronomical Applications

Authors: Bingjie Wang, Joel Leja, V. Ashley Villar, Joshua S. Speagle

Abstract: Flagship near-future surveys targeting $10^8-10^9$ galaxies across cosmic time will soon reveal the processes of galaxy assembly in unprecedented resolution. This creates an immediate computational challenge on effective analyses of the full data-set. With simulation-based inference (SBI), it is possible to attain complex posterior distributions with the accuracy of traditional methods but with a… ▽ More Flagship near-future surveys targeting $10^8-10^9$ galaxies across cosmic time will soon reveal the processes of galaxy assembly in unprecedented resolution. This creates an immediate computational challenge on effective analyses of the full data-set. With simulation-based inference (SBI), it is possible to attain complex posterior distributions with the accuracy of traditional methods but with a $>10^4$ increase in speed. However, it comes with a major limitation. Standard SBI requires the simulated data to have identical characteristics to the observed data, which is often violated in astronomical surveys due to inhomogeneous coverage and/or fluctuating sky and telescope conditions. In this work, we present a complete SBI-based methodology, ``SBI$^{++}$,'' for treating out-of-distribution measurement errors and missing data. We show that out-of-distribution errors can be approximated by using standard SBI evaluations and that missing data can be marginalized over using SBI evaluations over nearby data realizations in the training set. In addition to the validation set, we apply SBI$^{++}$ to galaxies identified in extragalactic images acquired by the James Webb Space Telescope, and show that SBI$^{++}$ can infer photometric redshifts at least as accurately as traditional sampling methods and crucially, better than the original SBI algorithm using training data with a wide range of observational errors. SBI$^{++}$ retains the fast inference speed of $\sim$1 sec for objects in the observational training set distribution, and additionally permits parameter inference outside of the trained noise and data at $\sim$1 min per object. This expanded regime has broad implications for future applications to astronomical surveys. △ Less

Submitted 10 July, 2023; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: Accepted for publication in ApJL; 13 pages, 5 figures. Code and a Jupyter tutorial are made publicly available at https://github.com/wangbingjie/sbi_pp

Journal ref: The Astrophysical Journal Letters, 952, L10 (2023)

Latent Stochastic Differential Equations for Modeling Quasar Variability and Inferring Black Hole Properties

Authors: Joshua Fagin, Ji Won Park, Henry Best, James Hung-Hsu Chan, K. E Saavik Ford, Matthew J. Graham, V. Ashley Villar, Shirley Ho, Matthew O'Dowd

Abstract: Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar's brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described… ▽ More Quasars are bright and unobscured active galactic nuclei (AGN) thought to be powered by the accretion of matter around supermassive black holes at the centers of galaxies. The temporal variability of a quasar's brightness contains valuable information about its physical properties. The UV/optical variability is thought to be a stochastic process, often represented as a damped random walk described by a stochastic differential equation (SDE). Upcoming wide-field telescopes such as the Rubin Observatory Legacy Survey of Space and Time (LSST) are expected to observe tens of millions of AGN in multiple filters over a ten year period, so there is a need for efficient and automated modeling techniques that can handle the large volume of data. Latent SDEs are machine learning models well suited for modeling quasar variability, as they can explicitly capture the underlying stochastic dynamics. In this work, we adapt latent SDEs to jointly reconstruct multivariate quasar light curves and infer their physical properties such as the black hole mass, inclination angle, and temperature slope. Our model is trained on realistic simulations of LSST ten year quasar light curves, and we demonstrate its ability to reconstruct quasar light curves even in the presence of long seasonal gaps and irregular sampling across different bands, outperforming a multioutput Gaussian process regression baseline. Our method has the potential to provide a deeper understanding of the physical properties of quasars and is applicable to a wide range of other multivariate time series with missing data and irregular sampling. △ Less

Submitted 16 April, 2024; v1 submitted 9 April, 2023; originally announced April 2023.

Comments: 22 pages, 12 figures

Journal ref: The Astrophysical Journal, Volume 965, Number 2, April 2024

Neutron star-black hole mergers in next generation gravitational-wave observatories

Authors: Ish Gupta, Ssohrab Borhanian, Arnab Dhani, Debatri Chattopadhyay, Rahul Kashyap, V. Ashley Villar, B. S. Sathyaprakash

Abstract: Observations by the current generation of gravitational-wave detectors have been pivotal in expanding our understanding of the universe. Although tens of exciting compact binary mergers have been observed, neutron star-black hole (NSBH) mergers remained elusive until they were first confidently detected in 2020. The number of NSBH detections is expected to increase with sensitivity improvements of… ▽ More Observations by the current generation of gravitational-wave detectors have been pivotal in expanding our understanding of the universe. Although tens of exciting compact binary mergers have been observed, neutron star-black hole (NSBH) mergers remained elusive until they were first confidently detected in 2020. The number of NSBH detections is expected to increase with sensitivity improvements of the current detectors and the proposed construction of new observatories over the next decade. In this work, we explore the NSBH detection and measurement capabilities of these upgraded detectors and new observatories using the following metrics: network detection efficiency and detection rate as a function of redshift, distributions of the signal-to-noise ratios, the measurement accuracy of intrinsic and extrinsic parameters, the accuracy of sky position measurement, and the number of early-warning alerts that can be sent to facilitate the electromagnetic follow-up. Additionally, we evaluate the prospects of performing multi-messenger observations of NSBH systems by reporting the number of expected kilonova detections with the Vera C. Rubin Observatory and the Nancy Grace Roman Space Telescope. We find that as many as $\mathcal{O}(10)$ kilonovae can be detected by these two telescopes every year, depending on the population of the NSBH systems and the equation of state of neutron stars. △ Less

Submitted 25 July, 2023; v1 submitted 20 January, 2023; originally announced January 2023.

Comments: 30 pages, 15 figures, 15 tables

Journal ref: Phys. Rev. D 107, 124007 (2023)

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

SN 2022ann: A type Icn supernova from a dwarf galaxy that reveals helium in its circumstellar environment

Authors: K. W. Davis, K. Taggart, S. Tinyanont, R. J. Foley, V. A. Villar, L. Izzo, C. R. Angus, M. J. Bustamante-Rosell, D. A. Coulter, N. Earl, D. Farias, J. Hjorth, M. E. Huber, D. O. Jones, P. L. Kelly, C. D. Kilpatrick, D. Langeroodi, H. -Y. Miao, C. M. Pellegrino, E. Ramirez-Ruiz, C. L. Ransome, S. Rest, S. N. Sharief, M. R. Siebert, G. Terreran , et al. (43 additional authors not shown)

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

Submitted 9 November, 2022; originally announced November 2022.

Comments: 21 pages, 13 figures, to be submitted to MNRAS

Amortized Bayesian Inference for Supernovae in the Era of the Vera Rubin Observatory Using Normalizing Flows

Authors: V. Ashley Villar

Abstract: The Vera Rubin Observatory, set to begin observations in mid-2024, will increase our discovery rate of supernovae to well over one million annually. There has been a significant push to develop new methodologies to identify, classify and ultimately understand the millions of supernovae discovered with the Rubin Observatory. Here, we present the first simulation-based inference method using normali… ▽ More The Vera Rubin Observatory, set to begin observations in mid-2024, will increase our discovery rate of supernovae to well over one million annually. There has been a significant push to develop new methodologies to identify, classify and ultimately understand the millions of supernovae discovered with the Rubin Observatory. Here, we present the first simulation-based inference method using normalizing flows, trained to rapidly infer the parameters of toy supernovae model in multivariate, Rubin-like datastreams. We find that our method is well-calibrated compared to traditional inference methodologies (specifically MCMC), requiring only one-ten-thousandth of the CPU hours during test time. △ Less

Submitted 8 November, 2022; originally announced November 2022.

Comments: 5 pages, accepted in the Neurips Machine Learning and the Physical Sciences conference

The First Two Years of FLEET: an Active Search for Superluminous Supernovae

Authors: Sebastian Gomez, Edo Berger, Peter K. Blanchard, Griffin Hosseinzadeh, Matt Nicholl, Daichi Hiramatsu, V. Ashley Villar, Yao Yin

Abstract: In November 2019 we began operating FLEET (Finding Luminous and Exotic Extragalactic Transients), a machine learning algorithm designed to photometrically identify Type I superluminous supernovae (SLSNe) in transient alert streams. Using FLEET, we spectroscopically classified 21 of the 50 SLSNe identified worldwide between November 2019 and January 2022. Based on our original algorithm, we anticip… ▽ More In November 2019 we began operating FLEET (Finding Luminous and Exotic Extragalactic Transients), a machine learning algorithm designed to photometrically identify Type I superluminous supernovae (SLSNe) in transient alert streams. Using FLEET, we spectroscopically classified 21 of the 50 SLSNe identified worldwide between November 2019 and January 2022. Based on our original algorithm, we anticipated that FLEET would achieve a purity of about 50\% for transients with a probability of being a SLSN, \pslsn$>0.5$; the true on-sky purity we obtained is closer to 80\%. Similarly, we anticipated FLEET could reach a completeness of about 30\%, and we indeed measure an upper limit on the completeness of $\approx 33$\%. Here, we present FLEET 2.0, an updated version of FLEET trained on 4,780 transients (almost 3 times more than in FLEET 1.0). FLEET 2.0 has a similar predicted purity to FLEET 1.0, but outperforms FLEET 1.0 in terms of completeness, which is now closer to $\approx 40$\% for transients with \pslsn$>0.5$. Additionally, we explore possible systematics that might arise from the use of FLEET for target selection. We find that the population of SLSNe recovered by FLEET is mostly indistinguishable from the overall SLSN population, in terms of physical and most observational parameters. We provide FLEET as an open source package on GitHub https://github.com/gmzsebastian/FLEET △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 10 pages, 3 figures, submitted to ApJ

Identifying Tidal Disruption Events with an Expansion of the FLEET Machine Learning Algorithm

Authors: Sebastian Gomez, V. Ashley Villar, Edo Berger, Suvi Gezari, Sjoert van Velzen, Matt Nicholl, Peter K. Blanchard, Kate. D. Alexander

Abstract: We present an expansion of FLEET, a machine learning algorithm optimized to select transients that are most likely to be tidal disruption events (TDEs). FLEET is based on a random forest algorithm trained on the light curves and host galaxy information of 4,779 spectroscopically classified transients. For transients with a probability of being a TDE, \ptde$>0.5$, we can successfully recover TDEs w… ▽ More We present an expansion of FLEET, a machine learning algorithm optimized to select transients that are most likely to be tidal disruption events (TDEs). FLEET is based on a random forest algorithm trained on the light curves and host galaxy information of 4,779 spectroscopically classified transients. For transients with a probability of being a TDE, \ptde$>0.5$, we can successfully recover TDEs with a $\approx40$\% completeness and a $\approx30$\% purity when using the first 20 days of photometry, or a similar completeness and $\approx50$\% purity when including 40 days of photometry. We find that the most relevant features for differentiating TDEs from other transients are the normalized host separation, and the light curve $(g-r)$ color during peak. Additionally, we use FLEET to produce a list of the 39 most likely TDE candidates discovered by the Zwicky Transient Facility that remain currently unclassified. We explore the use of FLEET for the Legacy Survey of Space and Time on the Vera C. Rubin Observatory (\textit{Rubin}) and the \textit{Nancy Grace Roman Space Telescope} (\textit{Roman}). We simulate the \textit{Rubin} and \textit{Roman} survey strategies and estimate that $\sim 10^4$ TDEs could be discovered every year by \textit{Rubin}, and $\sim200$ TDEs per year by \textit{Roman}. Finally, we run FLEET on the TDEs in our \textit{Rubin} survey simulation and find that we can recover $\sim 30$\% of those at a redshift $z <0.5$ with \ptde$>0.5$. This translates to $\sim3,000$ TDEs per year that FLEET could uncover from \textit{Rubin}. FLEET is provided as a open source package on GitHub https://github.com/gmzsebastian/FLEET △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 17 pages, 12 figures, submitted to ApJ

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

Photometrically-Classified Superluminous Supernovae from the Pan-STARRS1 Medium Deep Survey: A Case Study for Science with Machine Learning-Based Classification

Authors: Brian Hsu, Griffin Hosseinzadeh, V. Ashley Villar, Edo Berger

Abstract: With the upcoming Vera C.~Rubin Observatory Legacy Survey of Space and Time (LSST), it is expected that only $\sim 0.1\%$ of all transients will be classified spectroscopically. To conduct studies of rare transients, such as Type I superluminous supernovae (SLSNe), we must instead rely on photometric classification. In this vein, here we carry out a pilot study of SLSNe from the Pan-STARRS1 Medium… ▽ More With the upcoming Vera C.~Rubin Observatory Legacy Survey of Space and Time (LSST), it is expected that only $\sim 0.1\%$ of all transients will be classified spectroscopically. To conduct studies of rare transients, such as Type I superluminous supernovae (SLSNe), we must instead rely on photometric classification. In this vein, here we carry out a pilot study of SLSNe from the Pan-STARRS1 Medium-Deep Survey (PS1-MDS) classified photometrically with our SuperRAENN and Superphot algorithms. We first construct a sub-sample of the photometric sample using a list of simple selection metrics designed to minimize contamination and ensure sufficient data quality for modeling. We then fit the multi-band light curves with a magnetar spin-down model using the Modular Open-Source Fitter for Transients (MOSFiT). Comparing the magnetar engine and ejecta parameter distributions of the photometric sample to those of the PS1-MDS spectroscopic sample and a larger literature spectroscopic sample, we find that these samples are overall consistent, but that the photometric sample extends to slower spins and lower ejecta masses, which correspond to lower luminosity events, as expected for photometric selection. While our PS1-MDS photometric sample is still smaller than the overall SLSN spectroscopic sample, our methodology paves the way to an orders-of-magnitude increase in the SLSN sample in the LSST era through photometric selection and study. △ Less

Submitted 17 July, 2023; v1 submitted 20 April, 2022; originally announced April 2022.

Comments: 13 pages, 6 figures, submitted to ApJ

Journal ref: The Astrophysical Journal, Volume 937, Number 1 (2022)

Machine Learning and Cosmology

Authors: Cora Dvorkin, Siddharth Mishra-Sharma, Brian Nord, V. Ashley Villar, Camille Avestruz, Keith Bechtol, Aleksandra Ćiprijanović, Andrew J. Connolly, Lehman H. Garrison, Gautham Narayan, Francisco Villaescusa-Navarro

Abstract: Methods based on machine learning have recently made substantial inroads in many corners of cosmology. Through this process, new computational tools, new perspectives on data collection, model development, analysis, and discovery, as well as new communities and educational pathways have emerged. Despite rapid progress, substantial potential at the intersection of cosmology and machine learning rem… ▽ More Methods based on machine learning have recently made substantial inroads in many corners of cosmology. Through this process, new computational tools, new perspectives on data collection, model development, analysis, and discovery, as well as new communities and educational pathways have emerged. Despite rapid progress, substantial potential at the intersection of cosmology and machine learning remains untapped. In this white paper, we summarize current and ongoing developments relating to the application of machine learning within cosmology and provide a set of recommendations aimed at maximizing the scientific impact of these burgeoning tools over the coming decade through both technical development as well as the fostering of emerging communities. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 32 pages

The Circumstellar Environments of Double-Peaked, Calcium-strong Supernovae 2021gno and 2021inl

Authors: Wynn Jacobson-Galán, Padma Venkatraman, Raffaella Margutti, David Khatami, Giacomo Terreran, Ryan J. Foley, Rodrigo Angulo, Charlotte R. Angus, Katie Auchettl, Peter K. Blanchard, Alexey Bobrick, Joe S. Bright, Cirilla D. Couch, David A. Coulter, Karoli Clever, Kyle W. Davis, Thomas de Boer, Lindsay DeMarchi, Sierra A. Dodd, David O. Jones, Jessica Johnson, Charles D. Kilpatrick, Nandita Khetan, Zhisen Lai, Danial Langeroodi , et al. (20 additional authors not shown)

Abstract: We present panchromatic observations and modeling of calcium-strong supernovae (SNe) 2021gno in the star-forming host galaxy NGC 4165 (D = 30.5 Mpc) and 2021inl in the outskirts of elliptical galaxy NGC 4923 (D = 80 Mpc), both monitored through the Young Supernova Experiment (YSE) transient survey. The multi-color light curves of both SNe show two peaks, the former peak being derived from shock co… ▽ More We present panchromatic observations and modeling of calcium-strong supernovae (SNe) 2021gno in the star-forming host galaxy NGC 4165 (D = 30.5 Mpc) and 2021inl in the outskirts of elliptical galaxy NGC 4923 (D = 80 Mpc), both monitored through the Young Supernova Experiment (YSE) transient survey. The multi-color light curves of both SNe show two peaks, the former peak being derived from shock cooling emission (SCE) and/or shock interaction with circumstellar material (CSM). The primary peak in SN 2021gno is coincident with luminous, rapidly decaying X-ray emission ($L_x = 5 \times 10^{41}$ erg s$^{-1}$) detected by Swift-XRT at $δt = 1$ day after explosion, this observation being the second ever detection of X-rays from a calcium-strong transient. We interpret the X-ray emission from SN 2021gno in the context of shock interaction with dense CSM that extends to $r < 3 \times 10^{14}$ cm. Based on modeling of the SN 2021gno X-ray spectrum, we calculate a CSM mass range of $M_{\rm CSM} = (0.3 - 1.6) \times 10^{-3}$ M$_{\odot}$ and particle densities of $n = (1-4) \times 10^{10}$ cm$^{-3}$. Radio non-detections of SN 2021gno indicate a low-density environment at larger radii ($r > 10^{16}$ cm) and a progenitor mass loss rate of $\dot{M} < 10^{-4}$ M$_{\odot}$ yr$^{-1}$, for $v_w = 500$ km s$^{-1}$. For radiation derived from SCE, modeling of the primary light curve peak in both SNe indicates an extended progenitor envelope mass and radius of $M_e = 0.02 - 0.05$ M$_{\odot}$ and $R_e = 30 - 230$ R$_{\odot}$. The explosion properties of SNe 2021gno and 2021inl suggest progenitor systems containing either a low-mass massive star or a white dwarf (WD), the former being unlikely for either object given the lack of star formation at both explosion sites. Furthermore, the progenitor environments of both SNe are consistent with explosion models for low-mass hybrid He/C/O WD + C/O WD binaries. △ Less

Submitted 7 March, 2022; originally announced March 2022.

Comments: 33 pages, 19 figures. Submitted to ApJ. Comments welcome!

Next Generation Observatories -- Report from the Dawn VI Workshop; October 5-7 2021

Authors: D. H. Shoemaker, Stefan Ballmer, Matteo Barsuglia, E. Berger, Emanuele Berti, Duncan A. Brown, Poonam Chandra, Matthew Evans, Ke Fang, Wen-fai Fong, Andreas Freise, Peter Fritschel, Jenny Greene, C. J. Horowitz, Jeff Kissel, Brian Lantz, Paul D. Lasky, Harald Lueck, M. Coleman Miller, Alexander H. Nitz, David Ottaway, Hiranya V. Peiris, Michele Punturo, D. H. Reitze, Gary H. Sanders , et al. (11 additional authors not shown)

Abstract: The workshop Dawn VI: Next Generation Observatories took place online over three days, 5-7 October, 2021. More than 200 physicists and astronomers attended to contribute to, and learn from, a discussion of next-generation ground-based gravitational-wave detectors. The program was centered on the next generation of ground-based gravitational-wave observatories and their synergy with the greater lan… ▽ More The workshop Dawn VI: Next Generation Observatories took place online over three days, 5-7 October, 2021. More than 200 physicists and astronomers attended to contribute to, and learn from, a discussion of next-generation ground-based gravitational-wave detectors. The program was centered on the next generation of ground-based gravitational-wave observatories and their synergy with the greater landscape of scientific observatories of the 2030s. Cosmic Explorer (CE), a concept developed with US National Science Foundation support, was a particular focus; Einstein Telescope (ET), the European next generation concept, is an important complement and partner in forming a network. The concluding summary of the meeting expressed the sentiment that the observational science accessible to CE and ET, also in combination with data from other non-GW observatories, will stimulate a very broad community of analysts and yield insights which are exciting given the access to GWs from the entire universe. The need, and desire, for closer collaboration between ET and CE was expressed; a three-detector network is optimal for delivering much of the science. The science opportunities afforded by CE and ET are broad and compelling, impacting a wide range of disciplines in physics and high energy astrophysics. There was a consensus that CE is a concept that can deliver the promised science. A strong endorsement of Cosmic Explorer, as described in the CE Horizon Study, is a primary outcome of DAWN VI. △ Less

Submitted 20 February, 2022; v1 submitted 23 December, 2021; originally announced December 2021.

Comments: Proceedings of the Dawn VI workshop

A New Classification Model for the ZTF Catalog of Periodic Variable Stars

Authors: Siu-Hei Cheung, V. Ashley Villar, Ho-Sang Chan, Shirley Ho

Abstract: Using the second data release from the Zwicky Transient Facility (ZTF, Bellm et al. 2019), Chen et al. (2020) created a ZTF Catalog of Periodic Variable Stars (ZTF CPVS) of 781, 602 periodic variables stars (PVSs) with 11 class labels. Here, we provide a new classification model of PVSs in the ZTF CPVS using a convolutional variational autoencoder and hierarchical random forest. We cross-match the… ▽ More Using the second data release from the Zwicky Transient Facility (ZTF, Bellm et al. 2019), Chen et al. (2020) created a ZTF Catalog of Periodic Variable Stars (ZTF CPVS) of 781, 602 periodic variables stars (PVSs) with 11 class labels. Here, we provide a new classification model of PVSs in the ZTF CPVS using a convolutional variational autoencoder and hierarchical random forest. We cross-match the sky-coordinate of PVSs in the ZTF CPVS with those presented in the SIMBAD catalog. We identify non-stellar objects that are not previously classified, including extragalactic objects such as Quasi-Stellar Objects, Active Galactic Nuclei, supernovae and planetary nebulae. We then create a new labelled training set with 13 classes in two levels. We obtain a reasonable level of completeness (> 90 %) for certain classes of PVSs, although we have poorer completeness in other classes (~ 40 % in some cases). Our new labels for the ZTF CPVS are available via Zenodo Cheung et al. (2021). △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 4 pages, 1 figure

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

TESSreduce: transient focused TESS data reduction pipeline

Authors: R. Ridden-Harper, A. Rest, R. Hounsell, T. E. Müller-Bravo, Q. Wang, V. A. Villar

Abstract: Since its launch, TESS has provided high cadence observations for objects across the sky. Although high cadence TESS observations provide a unique possibility to study the rapid time evolution of numerous objects, artifacts in the data make it particularly challenging to use in studying transients. Furthermore, the broadband red filter of TESS, makes calibrating it to physical flux units, or magni… ▽ More Since its launch, TESS has provided high cadence observations for objects across the sky. Although high cadence TESS observations provide a unique possibility to study the rapid time evolution of numerous objects, artifacts in the data make it particularly challenging to use in studying transients. Furthermore, the broadband red filter of TESS, makes calibrating it to physical flux units, or magnitudes, challenging. Here we present TESSreduce an open-source, and user-friendly Python package which is built to lower the barrier to entry for transient science with TESS. In a few commands users can produce a reliable TESS light curve, accounting for systematic biases that are present in other models (such as instrument drift and the varied TESS background) and calculate a zeropoint to percent level precision. With this package anyone can use TESS for science, such as studying rapid transients and constraining progenitors of supernovae. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 14 pages, 14 figures

A Convolutional Autoencoder-Based Pipeline for Anomaly Detection and Classification of Periodic Variables

Authors: H. S. Chan, S. H. Cheung, V. Ashley Villar, Shirley Ho

Abstract: The periodic pulsations of stars teach us about their underlying physical process. We present a convolutional autoencoder-based pipeline as an automatic approach to search for out-of-distribution anomalous periodic variables within The Zwicky Transient Facility Catalog of Periodic Variable Stars (ZTF CPVS). We use an isolation forest to rank each periodic variable by its anomaly score. Our overall… ▽ More The periodic pulsations of stars teach us about their underlying physical process. We present a convolutional autoencoder-based pipeline as an automatic approach to search for out-of-distribution anomalous periodic variables within The Zwicky Transient Facility Catalog of Periodic Variable Stars (ZTF CPVS). We use an isolation forest to rank each periodic variable by its anomaly score. Our overall most anomalous events have a unique physical origin: they are mostly highly variable and irregular evolved stars. Multiwavelength data suggest that they are most likely Red Giant or Asymptotic Giant Branch stars concentrated in the Milky Way galactic disk. Furthermore, we show how the learned latent features can be used for the classification of periodic variables through a hierarchical random forest. This novel semi-supervised approach allows astronomers to identify the most anomalous events within a given physical class, significantly increasing the potential for scientific discovery. △ Less

Submitted 27 November, 2021; originally announced November 2021.

Comments: 7 pages, 4 figures

"Super-Kilonovae" from Massive Collapsars as Signatures of Black-Hole Birth in the Pair-instability Mass Gap

Authors: Daniel M. Siegel, Aman Agarwal, Jennifer Barnes, Brian D. Metzger, Mathieu Renzo, V. Ashley Villar

Abstract: The core collapse of rapidly rotating massive ~10 Msun stars ("collapsars"), and resulting formation of hyper-accreting black holes, are a leading model for the central engines of long-duration gamma-ray bursts (GRB) and promising sources of r-process nucleosynthesis. Here, we explore the signatures of collapsars from progenitors with extremely massive helium cores >130 Msun above the pair-instabi… ▽ More The core collapse of rapidly rotating massive ~10 Msun stars ("collapsars"), and resulting formation of hyper-accreting black holes, are a leading model for the central engines of long-duration gamma-ray bursts (GRB) and promising sources of r-process nucleosynthesis. Here, we explore the signatures of collapsars from progenitors with extremely massive helium cores >130 Msun above the pair-instability mass gap. While rapid collapse to a black hole likely precludes a prompt explosion in these systems, we demonstrate that disk outflows can generate a large quantity (up to >50 Msun) of ejecta, comprised of >5-10 Msun in r-process elements and ~0.1-1 Msun of $^{56}$Ni, expanding at velocities ~0.1c. Radioactive heating of the disk-wind ejecta powers an optical/infrared transient, with a characteristic luminosity $\sim 10^{42}$ erg s$^{-1}$ and spectral peak in the near-infrared (due to the high optical/UV opacities of lanthanide elements) similar to kilonovae from neutron star mergers, but with longer durations $\gtrsim$ 1 month. These "super-kilonovae" (superKNe) herald the birth of massive black holes >60 Msun, which, as a result of disk wind mass-loss, can populate the pair-instability mass gap 'from above' and could potentially create the binary components of GW190521. SuperKNe could be discovered via wide-field surveys such as those planned with the Roman Space Telescope or via late-time infrared follow-up observations of extremely energetic GRBs. Gravitational waves of frequency ~0.1-50 Hz from non-axisymmetric instabilities in self-gravitating massive collapsar disks are potentially detectable by proposed third-generation intermediate and high-frequency observatories at distances up to hundreds of Mpc; in contrast to the "chirp" from binary mergers, the collapsar gravitational-wave signal decreases in frequency as the disk radius grows ("sad trombone"). △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 36 pages, 21 figures, 2 tables

Target of Opportunity Observations of Gravitational Wave Events with Vera C. Rubin Observatory

Authors: Igor Andreoni, Raffaella Margutti, Om Sharan Salafia, B. Parazin, V. Ashley Villar, Michael W. Coughlin, Peter Yoachim, Kris Mortensen, Daniel Brethauer, S. J. Smartt, Mansi M. Kasliwal, Kate D. Alexander, Shreya Anand, E. Berger, Maria Grazia Bernardini, Federica B. Bianco, Peter K. Blanchard, Joshua S. Bloom, Enzo Brocato, Mattia Bulla, Regis Cartier, S. Bradley Cenko, Ryan Chornock, Christopher M. Copperwheat, Alessandra Corsi , et al. (30 additional authors not shown)

Abstract: The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exp… ▽ More The discovery of the electromagnetic counterpart to the binary neutron star merger GW170817 has opened the era of gravitational-wave multi-messenger astronomy. Rapid identification of the optical/infrared kilonova enabled a precise localization of the source, which paved the way to deep multi-wavelength follow-up and its myriad of related science results. Fully exploiting this new territory of exploration requires the acquisition of electromagnetic data from samples of neutron star mergers and other gravitational wave sources. After GW170817, the frontier is now to map the diversity of kilonova properties and provide more stringent constraints on the Hubble constant, and enable new tests of fundamental physics. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) can play a key role in this field in the 2020s, when an improved network of gravitational-wave detectors is expected to reach a sensitivity that will enable the discovery of a high rate of merger events involving neutron stars (about tens per year) out to distances of several hundred Mpc. We design comprehensive target-of-opportunity observing strategies for follow-up of gravitational-wave triggers that will make the Rubin Observatory the premier instrument for discovery and early characterization of neutron star and other compact object mergers, and yet unknown classes of gravitational wave events. △ Less

Submitted 20 April, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: Accepted for publication in ApJS. arXiv admin note: text overlap with arXiv:1812.04051

Extragalactic Millimeter Transients in the Era of Next Generation CMB Surveys

Authors: T. Eftekhari, E. Berger, B. D. Metzger, T. Laskar, V. A. Villar, K. D. Alexander, G. P. Holder, J. D. Vieira, N. Whitehorn, P. K. G. Williams

Abstract: The next generation of wide-field cosmic microwave background (CMB) surveys are uniquely poised to open a new window for time-domain astronomy in the millimeter band. Here we explore the discovery phase space for extragalactic transients with near-term and future CMB experiments to characterize the expected population. We use existing millimeter-band light curves of known transients (gamma-ray bur… ▽ More The next generation of wide-field cosmic microwave background (CMB) surveys are uniquely poised to open a new window for time-domain astronomy in the millimeter band. Here we explore the discovery phase space for extragalactic transients with near-term and future CMB experiments to characterize the expected population. We use existing millimeter-band light curves of known transients (gamma-ray bursts, tidal disruption events, fast blue optical transients, neutron star mergers) and theoretical models, in conjunction with known and estimated volumetric rates. Using Monte Carlo simulations of various CMB survey designs (area, cadence, depth, duration) we estimate the detection rates and the resulting light curve characteristics. We find that existing and near-term surveys will find tens to hundreds of long-duration gamma-ray bursts (LGRBs), driven primarily by detections of the reverse shock emission, and including off-axis LGRBs. Next-generation experiments (CMB-S4, CMB-HD) will find tens of fast blue optical transients (FBOTs) in the nearby universe and will detect a few tidal disruption events. CMB-HD will additionally detect a small number of short gamma-ray bursts, where these will be discovered within the detection volume of next generation gravitational wave experiments like Cosmic Explorer. △ Less

Submitted 11 September, 2022; v1 submitted 11 October, 2021; originally announced October 2021.

Comments: 24 pages, 6 figures; accepted to ApJ

SN 2018agk: A Prototypical Type Ia Supernova with a Smooth Power-law Rise in Kepler (K2)

Authors: Qinan Wang, Armin Rest, Yossef Zenati, Ryan Ridden-Harper, Georgios Dimitriadis, Gautham Narayan, V. Ashley Villar, Mark R. Magee, Ryan J. Foley, Edward J. Shaya, Peter Garnavich, Lifan Wang, Lei Hu, Attila Bodi, Patrick Armstrong, Katie Auchettl, Thomas Barclay, Geert Barentsen, Zsófia Bognár, Joseph Brimacombe, Joanna Bulger, Jamison Burke, Peter Challis, Kenneth Chambers, David A. Coulter , et al. (51 additional authors not shown)

Abstract: We present the 30-min cadence Kepler/K2 light curve of the Type Ia supernova (SN Ia) SN 2018agk, covering approximately one week before explosion, the full rise phase and the decline until 40 days after peak. We additionally present ground-based observations in multiple bands within the same time range, including the 1-day cadence DECam observations within the first $\sim$5 days after the first li… ▽ More We present the 30-min cadence Kepler/K2 light curve of the Type Ia supernova (SN Ia) SN 2018agk, covering approximately one week before explosion, the full rise phase and the decline until 40 days after peak. We additionally present ground-based observations in multiple bands within the same time range, including the 1-day cadence DECam observations within the first $\sim$5 days after the first light. The Kepler early light curve is fully consistent with a single power-law rise, without evidence of any bump feature. We compare SN 2018agk with a sample of other SNe~Ia without early excess flux from the literature. We find that SNe Ia without excess flux have slowly-evolving early colors in a narrow range ($g-i\approx -0.20\pm0.20$ mag) within the first $\sim 10$ days. On the other hand, among SNe Ia detected with excess, SN 2017cbv and SN 2018oh tend to be bluer, while iPTF16abc's evolution is similar to normal SNe Ia without excess in $g-i$. We further compare the Kepler light curve of SN 2018agk with companion-interaction models, and rule out the existence of a typical non-degenerate companion undergoing Roche-lobe overflow at viewing angles smaller than $45^{\circ}$. △ Less

Submitted 28 December, 2021; v1 submitted 31 August, 2021; originally announced August 2021.

Comments: 20 pages, 14 figures, 5 tables. Published in ApJ

Journal ref: The Astrophysical Journal, 2021, Volume 923, Number 2

Evidence for X-ray Emission in Excess to the Jet Afterglow Decay 3.5 yrs After the Binary Neutron Star Merger GW 170817: A New Emission Component

Authors: A. Hajela, R. Margutti, J. S. Bright, K. D. Alexander, B. D. Metzger, V. Nedora, A. Kathirgamaraju, B. Margalit, D. Radice, C. Guidorzi, E. Berger, A. MacFadyen, D. Giannios, R. Chornock, I. Heywood, L. Sironi, O. Gottlieb, D. Coppejans, T. Laskar, Y. Cendes, R. Barniol Duran, T. Eftekhari, W. Fong, A. McDowell, M. Nicholl , et al. (12 additional authors not shown)

Abstract: For the first $\sim3$ years after the binary neutron star merger event GW 170817 the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n $< 0.01\,\rm{cm^{-3}}$. We report on observational evidence for an excess of X-ray emission at $δt>900$ days after the merger. With… ▽ More For the first $\sim3$ years after the binary neutron star merger event GW 170817 the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n $< 0.01\,\rm{cm^{-3}}$. We report on observational evidence for an excess of X-ray emission at $δt>900$ days after the merger. With $L_x\approx5\times 10^{38}\,\rm{erg\,s^{-1}}$ at 1234 days, the recently detected X-ray emission represents a $\ge 3.2\,σ$ (Gaussian equivalent) deviation from the universal post jet-break model that best fits the multi-wavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance $\ge 3.1\,σ$, depending on the fireball collimation, with the most realistic models showing excesses at the level of $\ge 3.7\,σ$. A lack of detectable 3 GHz radio emission suggests a harder broad-band spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e. a kilonova afterglow. In this context, we present a set of ab-initio numerical-relativity BNS merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented. △ Less

Submitted 5 March, 2022; v1 submitted 5 April, 2021; originally announced April 2021.

Comments: 30 pages, 13 figures, Accepted for publication in ApJL

A Deep Learning Approach for Active Anomaly Detection of Extragalactic Transients

Authors: V. Ashley Villar, Miles Cranmer, Edo Berger, Gabriella Contardo, Shirley Ho, Griffin Hosseinzadeh, Joshua Yao-Yu Lin

Abstract: There is a shortage of multi-wavelength and spectroscopic followup capabilities given the number of transient and variable astrophysical events discovered through wide-field, optical surveys such as the upcoming Vera C. Rubin Observatory. From the haystack of potential science targets, astronomers must allocate scarce resources to study a selection of needles in real time. Here we present a variat… ▽ More There is a shortage of multi-wavelength and spectroscopic followup capabilities given the number of transient and variable astrophysical events discovered through wide-field, optical surveys such as the upcoming Vera C. Rubin Observatory. From the haystack of potential science targets, astronomers must allocate scarce resources to study a selection of needles in real time. Here we present a variational recurrent autoencoder neural network to encode simulated Rubin Observatory extragalactic transient events using 1% of the PLAsTiCC dataset to train the autoencoder. Our unsupervised method uniquely works with unlabeled, real time, multivariate and aperiodic data. We rank 1,129,184 events based on an anomaly score estimated using an isolation forest. We find that our pipeline successfully ranks rarer classes of transients as more anomalous. Using simple cuts in anomaly score and uncertainty, we identify a pure (~95% pure) sample of rare transients (i.e., transients other than Type Ia, Type II and Type Ibc supernovae) including superluminous and pair-instability supernovae. Finally, our algorithm is able to identify these transients as anomalous well before peak, enabling real-time follow up studies in the era of the Rubin Observatory. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: 13 pages,9 figures, submitted to AAS Journals. Comments welcome!

The Luminous and Double-Peaked Type Ic Supernova 2019stc: Evidence for Multiple Energy Sources

Authors: Sebastian Gomez, Edo Berger, Griffin Hosseinzadeh, Peter K. Blanchard, Matt Nicholl, V. Ashley Villar

Abstract: We present optical photometry and spectroscopy of SN\,2019stc (=ZTF19acbonaa), an unusual Type Ic supernova (SN Ic) at a redshift of $z=0.117$. SN\,2019stc exhibits a broad double-peaked light curve, with the first peak having an absolute magnitude of $M_r=-20.0$ mag, and the second peak, about 80 rest-frame days later, $M_r=-19.2$ mag. The total radiated energy is large,… ▽ More We present optical photometry and spectroscopy of SN\,2019stc (=ZTF19acbonaa), an unusual Type Ic supernova (SN Ic) at a redshift of $z=0.117$. SN\,2019stc exhibits a broad double-peaked light curve, with the first peak having an absolute magnitude of $M_r=-20.0$ mag, and the second peak, about 80 rest-frame days later, $M_r=-19.2$ mag. The total radiated energy is large, $E_{\rm rad}\approx 2.5\times 10^{50}$ erg. Despite its large luminosity, approaching those of Type I superluminous supernovae (SLSNe), SN\,2019stc exhibits a typical SN Ic spectrum, bridging the gap between SLSNe and SNe Ic. The spectra indicate the presence of Fe-peak elements, but modeling of the first light curve peak with radioactive heating alone leads to an unusually high nickel mass fraction of $f_{\rm Ni}\approx 31\%$ ($M_{\rm Ni}\approx 3.2$ M$_\odot$). Instead, if we model the first peak with a combined magnetar spin-down and radioactive heating model we find a better match with $M_{\rm ej}\approx 4$ M$_\odot$, a magnetar spin period of $P_{\rm spin}\approx 7.2$ ms and magnetic field of $B\approx 10^{14}$ G, and $f_{\rm Ni}\lesssim 0.2$ (consistent with SNe Ic). The prominent second peak cannot be naturally accommodated with radioactive heating or magnetar spin-down, but instead can be explained as circumstellar interaction with $\approx 0.7$ $M_\odot$ of hydrogen-free material located $\approx 400$ AU from the progenitor. Including the remnant mass leads to a CO core mass prior to explosion of $\approx 6.5$ M$_\odot$. The host galaxy has a metallicity of $\approx 0.26$ Z$_\odot$, low for SNe Ic but consistent with SLSNe. Overall, we find that SN\,2019stc is a transition object between normal SNe Ic and SLSNe. △ Less

Submitted 29 July, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: 14 pages, 13 figures, Accepted to ApJ

Journal ref: The Astrophysical Journal, 2021 June 4, Volume 913, Number 2

A Late-Time Galaxy-Targeted Search for the Radio Counterpart of GW190814

Authors: K. D. Alexander, G. Schroeder, K. Paterson, W. Fong, P. Cowperthwaite, S. Gomez, B. Margalit, R. Margutti, E. Berger, P. Blanchard, R. Chornock, T. Eftekhari, T. Laskar, B. D. Metzger, M. Nicholl, V. A. Villar, P. K. G. Williams

Abstract: GW190814 was a compact object binary coalescence detected in gravitational waves by Advanced LIGO and Advanced Virgo that garnered exceptional community interest due to its excellent localization and the uncertain nature of the binary's lighter-mass component (either the heaviest known neutron star, or the lightest known black hole). Despite extensive follow up observations, no electromagnetic cou… ▽ More GW190814 was a compact object binary coalescence detected in gravitational waves by Advanced LIGO and Advanced Virgo that garnered exceptional community interest due to its excellent localization and the uncertain nature of the binary's lighter-mass component (either the heaviest known neutron star, or the lightest known black hole). Despite extensive follow up observations, no electromagnetic counterpart has been identified. Here we present new radio observations of 75 galaxies within the localization volume at $Δt\approx 35-266$ days post-merger. Our observations cover $\sim32$% of the total stellar luminosity in the final localization volume and extend to later timescales than previously-reported searches, allowing us to place the deepest constraints to date on the existence of a radio afterglow from a highly off-axis relativistic jet launched during the merger (assuming that the merger occurred within the observed area). For a viewing angle of $\sim46^{\circ}$ (the best-fit binary inclination derived from the gravitational wave signal) and assumed electron and magnetic field energy fractions of $ε_e=0.1$ and $ε_B=0.01$, we can rule out a typical short gamma-ray burst-like Gaussian jet with isotropic-equivalent kinetic energy $2\times10^{51}$ erg propagating into a constant density medium $n\gtrsim0.01$ cm$^{-3}$. These are the first limits resulting from a galaxy-targeted search for a radio counterpart to a gravitational wave event, and we discuss the challenges, and possible advantages, of applying similar search strategies to future events using current and upcoming radio facilities. △ Less

Submitted 16 September, 2021; v1 submitted 17 February, 2021; originally announced February 2021.

Comments: Accepted to ApJ; text and figures updated to match accepted version. 20 pages, 6 figures, 1 table

Detection and Parameter Estimation of Gravitational Waves from Binary Neutron-Star Mergers in Real LIGO Data using Deep Learning

Authors: Plamen G. Krastev, Kiranjyot Gill, V. Ashley Villar, Edo Berger

Abstract: One of the key challenges of real-time detection and parameter estimation of gravitational waves from compact binary mergers is the computational cost of conventional matched-filtering and Bayesian inference approaches. In particular, the application of these methods to the full signal parameter space available to the gravitational-wave detectors, and/or real-time parameter estimation is computati… ▽ More One of the key challenges of real-time detection and parameter estimation of gravitational waves from compact binary mergers is the computational cost of conventional matched-filtering and Bayesian inference approaches. In particular, the application of these methods to the full signal parameter space available to the gravitational-wave detectors, and/or real-time parameter estimation is computationally prohibitive. On the other hand, rapid detection and inference are critical for prompt follow-up of the electromagnetic and astro-particle counterparts accompanying important transients, such as binary neutron-star and black-hole neutron-star mergers. Training deep neural networks to identify specific signals and learn a computationally efficient representation of the mapping between gravitational-wave signals and their parameters allows both detection and inference to be done quickly and reliably, with high sensitivity and accuracy. In this work we apply a deep-learning approach to rapidly identify and characterize transient gravitational-wave signals from binary neutron-star mergers in real LIGO data. We show for the first time that artificial neural networks can promptly detect and characterize binary neutron star gravitational-wave signals in real LIGO data, and distinguish them from noise and signals from coalescing black-hole binaries. We illustrate this key result by demonstrating that our deep-learning framework classifies correctly all gravitational-wave events from the Gravitational-Wave Transient Catalog, GWTC-1 [Phys. Rev. X 9 (2019), 031040]. These results emphasize the importance of using realistic gravitational-wave detector data in machine learning approaches, and represent a step towards achieving real-time detection and inference of gravitational waves. △ Less

Submitted 28 July, 2021; v1 submitted 23 December, 2020; originally announced December 2020.

Comments: 10 pages, 8 figures, 1 table. Updated grant information. Published in Physics Letters B

Journal ref: Phys. Lett. B 815, 136161 (2021)

Anomaly Detection for Multivariate Time Series of Exotic Supernovae

Authors: V. Ashley Villar, Miles Cranmer, Gabriella Contardo, Shirley Ho, Joshua Yao-Yu Lin

Abstract: Supernovae mark the explosive deaths of stars and enrich the cosmos with heavy elements. Future telescopes will discover thousands of new supernovae nightly, creating a need to flag astrophysically interesting events rapidly for followup study. Ideally, such an anomaly detection pipeline would be independent of our current knowledge and be sensitive to unexpected phenomena. Here we present an unsu… ▽ More Supernovae mark the explosive deaths of stars and enrich the cosmos with heavy elements. Future telescopes will discover thousands of new supernovae nightly, creating a need to flag astrophysically interesting events rapidly for followup study. Ideally, such an anomaly detection pipeline would be independent of our current knowledge and be sensitive to unexpected phenomena. Here we present an unsupervised method to search for anomalous time series in real time for transient, multivariate, and aperiodic signals. We use a RNN-based variational autoencoder to encode supernova time series and an isolation forest to search for anomalous events in the learned encoded space. We apply this method to a simulated dataset of 12,159 supernovae, successfully discovering anomalous supernovae and objects with catastrophically incorrect redshift measurements. This work is the first anomaly detection pipeline for supernovae which works with online datastreams. △ Less

Submitted 21 October, 2020; originally announced October 2020.

Comments: 6 pages, 2 figures, written for non-astronomers, submitted to the NeurIPS workshop Machine Learning and the Physical Sciences. Comments welcome!!

Late-Time Radio and Millimeter Observations of Superluminous Supernovae and Long Gamma Ray Bursts: Implications for Obscured Star Formation, Central Engines, and Fast Radio Bursts

Authors: T. Eftekhari, B. Margalit, C. M. B. Omand, E. Berger, P. K. Blanchard, P. Demorest, B. D. Metzger, K. Murase, M. Nicholl, V. A. Villar, P. K. G. Williams, K. D. Alexander, S. Chatterjee, D. L. Coppejans, J. M. Cordes, S. Gomez, G. Hosseinzadeh, B. Hsu, K. Kashiyama, R. Margutti, Y. Yin

Abstract: We present the largest and deepest late-time radio and millimeter survey to date of superluminous supernovae (SLSNe) and long duration gamma-ray bursts (LGRBs) to search for associated non-thermal synchrotron emission. Using the Karl G. Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA), we observed 43 sources at 6 and 100 GHz on a timescale of $\sim 1 - 19$… ▽ More We present the largest and deepest late-time radio and millimeter survey to date of superluminous supernovae (SLSNe) and long duration gamma-ray bursts (LGRBs) to search for associated non-thermal synchrotron emission. Using the Karl G. Jansky Very Large Array (VLA) and the Atacama Large Millimeter/submillimeter Array (ALMA), we observed 43 sources at 6 and 100 GHz on a timescale of $\sim 1 - 19$ yr post-explosion. We do not detect radio/mm emission from any of the sources, with the exception of a 6 GHz detection of PTF10hgi (Eftekhari et al. 2019), as well as the detection of 6 GHz emission near the location of the SLSN PTF12dam, which we associate with its host galaxy. We use our data to place constraints on central engine emission due to magnetar wind nebulae and off-axis relativistic jets. We also explore non-relativistic emission from the SN ejecta, and place constraints on obscured star formation in the host galaxies. In addition, we conduct a search for fast radio bursts (FRBs) from some of the sources using VLA Phased-Array observations; no FRBs are detected to a limit of $16$ mJy ($7σ$; 10 ms duration) in about 40 min on source per event. A comparison to theoretical models suggests that continued radio monitoring may lead to detections of persistent radio emission on timescales of $\gtrsim {\rm decade}$. △ Less

Submitted 22 December, 2021; v1 submitted 13 October, 2020; originally announced October 2020.

Comments: 30 pages; 12 figures; accepted to ApJ

FLEET: A Redshift-Agnostic Machine Learning Pipeline to Rapidly Identify Hydrogen-Poor Superluminous Supernovae

Authors: Sebastian Gomez, Edo Berger, Peter K. Blanchard, Griffin Hosseinzadeh, Matt Nicholl, V. Ashley Villar, Yao Yin

Abstract: Over the past decade wide-field optical time-domain surveys have increased the discovery rate of transients to the point that $\lesssim 10\%$ are being spectroscopically classified. Despite this, these surveys have enabled the discovery of new and rare types of transients, most notably the class of hydrogen-poor superluminous supernovae (SLSN-I), with about 150 events confirmed to date. Here we pr… ▽ More Over the past decade wide-field optical time-domain surveys have increased the discovery rate of transients to the point that $\lesssim 10\%$ are being spectroscopically classified. Despite this, these surveys have enabled the discovery of new and rare types of transients, most notably the class of hydrogen-poor superluminous supernovae (SLSN-I), with about 150 events confirmed to date. Here we present a machine-learning classification algorithm targeted at rapid identification of a pure sample of SLSN-I to enable spectroscopic and multi-wavelength follow-up. This algorithm is part of the FLEET (Finding Luminous and Exotic Extragalactic Transients) observational strategy. It utilizes both light curve and contextual information, but without the need for a redshift, to assign each newly-discovered transient a probability of being a SLSN-I. This classifier can achieve a maximum purity of about 85\% (with 20\% completeness) when observing a selection of SLSN-I candidates. Additionally, we present two alternative classifiers that use either redshifts or complete light curves and can achieve an even higher purity and completeness. At the current discovery rate, the FLEET algorithm can provide about $20$ SLSN-I candidates per year for spectroscopic follow-up with 85\% purity; with the Legacy Survey of Space and Time we anticipate this will rise to more than $\sim 10^3$ events per year. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: 17 pages, 12 figures, submitted to ApJ

Extended Self-similar Solution for Circumstellar Material-Supernova Ejecta Interaction

Authors: Brighten Jiang, Shuai Jiang, V. Ashley Villar

Abstract: In this note, we present a detailed self-similar solution to the interaction of a uniformly expanding gas and a stationary ambient medium, with an application to supernovae interacting with preexisting circumstellar media (Type IIn SNe). We implement the generalized solution into the Modular Open Source Fitter for Transients (MOSFiT), an open-source Python package for fitting extragalactic transie… ▽ More In this note, we present a detailed self-similar solution to the interaction of a uniformly expanding gas and a stationary ambient medium, with an application to supernovae interacting with preexisting circumstellar media (Type IIn SNe). We implement the generalized solution into the Modular Open Source Fitter for Transients (MOSFiT), an open-source Python package for fitting extragalactic transient light curves. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: published in RNAAS; implemented in MOSFiT

Journal ref: Res. Notes AAS 4 16 (2020)