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

PS1-11aop: Probing the Mass Loss History of a Luminous Interacting Supernova Prior to its Final Eruption with Multi-wavelength Observations

Authors: Adaeze L. Ibik, Maria R. Drout, Raffaela Margutti, David Matthews, V. Ashley Villar, Edo Berger, Ryan Chornock, Kate D. Alexander, Tarraneh Eftekhari, Tanmoy Laskar, Ragnhild Lunnan, Ryan J. Foley, David Jones, Dan Milisavljevic, Armin Rest, Daniel Scolnic, Peter K. G. Williams

Abstract: Luminous interacting supernovae are a class of stellar explosions whose progenitors underwent vigorous mass loss in the years prior to core-collapse. While the mechanism by which this material is ejected is still debated, obtaining the full density profile of the circumstellar medium (CSM) could reveal more about this process. Here, we present an extensive multi-wavelength study of PS1-11aop, a lu… ▽ More Luminous interacting supernovae are a class of stellar explosions whose progenitors underwent vigorous mass loss in the years prior to core-collapse. While the mechanism by which this material is ejected is still debated, obtaining the full density profile of the circumstellar medium (CSM) could reveal more about this process. Here, we present an extensive multi-wavelength study of PS1-11aop, a luminous and slowly declining Type IIn SN discovered by the PanSTARRS Medium Deep Survey. PS1-11aop had a peak r-band magnitude of $-$20.5\,mag, a total radiated energy $>$ 8$\times$10$^{50}$\,erg, and it exploded near the center of a star-forming galaxy with super-solar metallicity. We obtained multiple detections at the location of PS1-11aop in the radio and X-ray bands between 4 and 10\,years post-explosion, and if due to the SN, it is one of the most luminous radio supernovae identified to date. Taken together, the multiwavelength properties of PS1-11aop are consistent with a CSM density profile with multiple zones. The early optical emission is consistent with the supernova blastwave interacting with a dense and confined CSM shell which contains multiple solar masses of material that was likely ejected in the final $<$10-100 years prior to the explosion,($\sim$0.05$-$1.0 M$_{\odot}$yr$^{-1}$ at radii of $\lesssim$10$^{16}$\,cm). The radio observations, on the other hand, are consistent with a sparser environment ($\lesssim$2$\times 10^{-3}$ M$_{\odot}$yr$^{-1}$ at radii of $\sim$0.5-1$\times$10$^{17}$\,cm) -- thus probing the history of the progenitor star prior to its final mass loss episode. △ Less

Submitted 19 October, 2024; originally announced October 2024.

Comments: 37 pages, 17 figures

Explaining Non-Merger Gamma-Ray Bursts and Broad-Lined Supernovae with Close Binary Progenitors with Black Hole Central Engine

Authors: Christopher L. Fryer, Eric Burns, Anna Y. Q. Ho, Amy Y. Lien, Daniel A. Perley, Jada L. Vail, V. Ashley Villar

Abstract: For over 25 years, the origin of long-duration gamma-ray bursts (lGRBs) has been linked to the collapse of rotating massive stars. However, we have yet to pinpoint the stellar progenitor powering these transients. Moreover, the dominant engine powering the explosions remains open to debate. Observations of both lGRBs, supernovae associated with these GRBs, such as broad-line (BL) stripped-envelope… ▽ More For over 25 years, the origin of long-duration gamma-ray bursts (lGRBs) has been linked to the collapse of rotating massive stars. However, we have yet to pinpoint the stellar progenitor powering these transients. Moreover, the dominant engine powering the explosions remains open to debate. Observations of both lGRBs, supernovae associated with these GRBs, such as broad-line (BL) stripped-envelope (type Ic) supernovae (hereafter, Ic-BL) supernovae (SNe) and perhaps superluminous SNe, fast blue optical transients, and fast x-ray transients, may provide clues to both engines and progenitors. In this paper, we conduct a detailed study of the tight-binary formation scenario for lGRBs, comparing this scenario to other leading progenitor models. Combining this progenitor scenario with different lGRB engines, we can compare to existing data and make predictions for future observational tests. We find that the combination of the tight-binary progenitor scenario with the black hole accretion disk (BHAD) engine can explain lGRBs, low-luminosity GRBs, ultra-long GRBs, and Ic-BL. We discuss the various progenitor properties required for these different subclasses and note such systems would be future gravitational wave merger sources. We show that the current literature on other progenitor-engine scenarios cannot explain all of these transient classes with a single origin, motivating additional work. We find that the tight-binary progenitor with a magnetar engine is excluded by existing observations. The observations can be used to constrain the properties of stellar evolution, the nature of the GRB and the associated SN engines in lGRBs and Ic-BL. We discuss the future observations needed to constrain our understanding of these rare, but powerful, explosions. △ Less

Submitted 14 October, 2024; originally announced October 2024.

Comments: 20 pages, 10 figures, submitted to ApJ

Unveiling the Diversity of Type IIn Supernovae via Systematic Light Curve Modeling

Authors: C. L. Ransome, V. A. Villar

Abstract: Type IIn supernovae (SNeIIn) are a highly heterogeneous subclass of core-collapse supernovae, spectroscopically characterized by signatures of interaction with a dense circumstellar medium (CSM). Here we systematically model the light curves of 142 archival SNeIIn using MOSFiT (the Modular Open Source Fitter for Transients). We find that the observed and inferred properties of SNIIn are diverse, b… ▽ More Type IIn supernovae (SNeIIn) are a highly heterogeneous subclass of core-collapse supernovae, spectroscopically characterized by signatures of interaction with a dense circumstellar medium (CSM). Here we systematically model the light curves of 142 archival SNeIIn using MOSFiT (the Modular Open Source Fitter for Transients). We find that the observed and inferred properties of SNIIn are diverse, but there are some trends. The typical SN CSM is dense ($\sim$10$^{-12}$gcm$^{-3}$) with highly diverse CSM geometry, with a median CSM mass of $\sim$1M$_\odot$. The ejecta are typically massive ($\gtrsim10$M$_\odot$), suggesting massive progenitor systems. We find positive correlations between the CSM mass and the rise and fall times of SNeIIn. Furthermore there are positive correlations between the rise time and fall times and the $r$-band luminosity. We estimate the mass-loss rates of our sample (where spectroscopy is available) and find a high median mass-loss rate of $\sim$10$^{-2}$M$_\odot$yr$^{-1}$, with a range between 10$^{-4}$--1M$_\odot$yr$^{-1}$. These mass-loss rates are most similar to the mass loss from great eruptions of luminous blue variables, consistent with the direct progenitor detections in the literature. We also discuss the role that binary interactions may play, concluding that at least some of our SNeIIn may be from massive binary systems. Finally, we estimate a detection rate of 1.6$\times$10$^5$yr$^{-1}$ in the upcoming Legacy Survey of Space and Time at the Vera C. Rubin Observatory. △ Less

Submitted 16 September, 2024; originally announced September 2024.

Comments: 55 pages, 19 figures, 6 tables

Double "acct": a distinct double-peaked supernova matching pulsational pair-instability models

Authors: C. R. Angus, S. E. Woosley, R. J. Foley, M. Nicholl, V. A. Villar, K. Taggart, M. Pursiainen, P. Ramsden, S. Srivastav, H. F. Stevance, T. Moore, K. Auchettl, W. B. Hoogendam, N. Khetan, S. K. Yadavalli, G. Dimitriadis, A. Gagliano, M. R. Siebert, A. Aamer, T. de Boer, K. C. Chambers, A. Clocchiatti, D. A. Coulter, M. R. Drout, D. Farias , et al. (13 additional authors not shown)

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

Submitted 3 September, 2024; originally announced September 2024.

Comments: 31 pages, 14 figures. Submitted to ApJL, comments welcome

SN 2021foa: The "Flip-Flop" Type IIn / Ibn supernova

Authors: D. Farias, C. Gall, G. Narayan, S. Rest, V. A. Villar, C. R. Angus, K. Auchettl, K. W. Davis, R. Foley, A. Gagliano, J. Hjorth, L. Izzo, C. D. Kilpatrick, H . M. L. Perkins, E. Ramirez-Ruiz, C. L. Ransome, Sarangi. A., R. Yarza, D. A. Coulter, D. O. Jones, N. Khetan, A. Rest, M. R. Siebert, J. J. Swift, K. Taggart , et al. (7 additional authors not shown)

Abstract: We present a comprehensive analysis of the photometric and spectroscopic evolution of SN~2021foa, unique among the class of transitional supernovae for repeatedly changing its spectroscopic appearance from hydrogen-to-helium-to-hydrogen-dominated (IIn-to-Ibn-to-IIn) within 50 days past peak brightness. The spectra exhibit multiple narrow ($\approx$ 300--600~km~s$^{-1}$) absorption lines of hydroge… ▽ More We present a comprehensive analysis of the photometric and spectroscopic evolution of SN~2021foa, unique among the class of transitional supernovae for repeatedly changing its spectroscopic appearance from hydrogen-to-helium-to-hydrogen-dominated (IIn-to-Ibn-to-IIn) within 50 days past peak brightness. The spectra exhibit multiple narrow ($\approx$ 300--600~km~s$^{-1}$) absorption lines of hydrogen, helium, calcium and iron together with broad helium emission lines with a full-width-at-half-maximum (FWHM) of $\sim 6000$~km~s$^{-1}$. For a steady, wind-mass loss regime, light curve modeling results in an ejecta mass of $\sim 8$ M$_{\odot}$ and CSM mass below 1 M$_{\odot}$, and an ejecta velocity consistent with the FWHM of the broad helium lines. We obtain a mass-loss rate of $\approx 2$ M$_{\odot} {\rm yr}^{-1}$. This mass-loss rate is three orders of magnitude larger than derived for normal Type II SNe. We estimate that the bulk of the CSM of SN~2021foa must have been expelled within half a year, about 15 years ago. Our analysis suggests that SN~2021foa had a helium rich ejecta which swept up a dense shell of hydrogen rich CSM shortly after explosion. At about 60 days past peak brightness, the photosphere recedes through the dense ejecta-CSM region, occulting much of the red-shifted emission of the hydrogen and helium lines, which results in observed blue-shift ($\sim -3000$~km~s$^{-1}$). Strong mass loss activity prior to explosion, such as those seen in SN~2009ip-like objects and SN~2021foa as precursor emission, are the likely origin of a complex, multiple-shell CSM close to the progenitor star. △ Less

Submitted 2 September, 2024; originally announced September 2024.

Comments: Submitted to ApJ, 21 figures, 7 tables

Maven: A Multimodal Foundation Model for Supernova Science

Authors: Gemma Zhang, Thomas Helfer, Alexander T. Gagliano, Siddharth Mishra-Sharma, V. Ashley Villar

Abstract: A common setting in astronomy is the availability of a small number of high-quality observations, and larger amounts of either lower-quality observations or synthetic data from simplified models. Time-domain astrophysics is a canonical example of this imbalance, with the number of supernovae observed photometrically outpacing the number observed spectroscopically by multiple orders of magnitude. A… ▽ More A common setting in astronomy is the availability of a small number of high-quality observations, and larger amounts of either lower-quality observations or synthetic data from simplified models. Time-domain astrophysics is a canonical example of this imbalance, with the number of supernovae observed photometrically outpacing the number observed spectroscopically by multiple orders of magnitude. At the same time, no data-driven models exist to understand these photometric and spectroscopic observables in a common context. Contrastive learning objectives, which have grown in popularity for aligning distinct data modalities in a shared embedding space, provide a potential solution to extract information from these modalities. We present Maven, the first foundation model for supernova science. To construct Maven, we first pre-train our model to align photometry and spectroscopy from 0.5M synthetic supernovae using a constrastive objective. We then fine-tune the model on 4,702 observed supernovae from the Zwicky Transient Facility. Maven reaches state-of-the-art performance on both classification and redshift estimation, despite the embeddings not being explicitly optimized for these tasks. Through ablation studies, we show that pre-training with synthetic data improves overall performance. In the upcoming era of the Vera C. Rubin Observatory, Maven serves as a Rosetta Stone for leveraging large, unlabeled and multimodal time-domain datasets. △ Less

Submitted 29 August, 2024; originally announced August 2024.

Comments: code: https://github.com/ThomasHelfer/multimodal-supernovae data: https://huggingface.co/datasets/thelfer/multimodal_supernovae

Constraining Dust Formation in the Superluminous Supernova 2017gci with JWST Observations

Authors: Sebastian Gomez, Tea Temim, Ori Fox, V. Ashley Villar, Melissa Shahbandeh, Chris Ashall, Jacob E. Jencson, Danial Langeroodi, Ilse De Looze, Dan Milisavljevic, Justin Pierel, Armin Rest, Tamás Szalai, Samaporn Tinyanont

Abstract: We present JWST/MIRI observations of the Type I superluminous supernova (SLSN) 2017gci taken over 2000 rest-frame days after the supernova (SN) exploded, which represent the latest phase images taken of any known SLSN. We find that archival \WISE detections of SN\,2017gci taken 70 to 200 days after explosion are most likely explained by an IR dust echo from a $\sim 3 \times 10^{-4}$ M$_\odot$ shel… ▽ More We present JWST/MIRI observations of the Type I superluminous supernova (SLSN) 2017gci taken over 2000 rest-frame days after the supernova (SN) exploded, which represent the latest phase images taken of any known SLSN. We find that archival \WISE detections of SN\,2017gci taken 70 to 200 days after explosion are most likely explained by an IR dust echo from a $\sim 3 \times 10^{-4}$ M$_\odot$ shell of pre-existing dust, as opposed to freshly-formed dust. New JWST observations reveal IR emission in the field of SN\,2017gci, which we determine is most likely dominated by the host galaxy of the SN, based on the expected flux of the galaxy and the measurable separation between said emission and the location of the SN. Based on models for IR emission of carbonate dust, we place a $3σ$ upper limit of $0.83$ M$_\odot$ of dust formed in SN\,2017gci, with a lowest $1σ$ limit of $0.44$ M$_\odot$. Infrared (IR) detections of other SLSNe have suggested that SLSNe could be among the most efficient dust producers in the universe. Our results suggest that SLSNe do not necessarily form more dust than other types of SNe, but instead might have a more accelerated dust formation process. More IR observations of a larger sample of SLSNe will be required to determine how efficient dust production is in SLSNe. △ Less

Submitted 27 August, 2024; originally announced August 2024.

Comments: 11 pages, 5 figures, submitted to ApJ

Finding the Fuse: Prospects for the Detection and Characterization of Hydrogen-Rich Core-Collapse Supernova Precursor Emission with the LSST

Authors: A. Gagliano, E. Berger, V. A. Villar, D. Hiramatsu, R. Kessler, T. Matsumoto, A. Gilkis, E. Laplace

Abstract: Enhanced emission in the months to years preceding explosion has been detected for several core-collapse supernovae (SNe). Though the physical mechanisms driving the emission remain hotly debated, the light curves of detected events show long-lived ($\geq$50 days), plateau-like behavior, suggesting hydrogen recombination may significantly contribute to the total energy budget. The Vera C. Rubin Ob… ▽ More Enhanced emission in the months to years preceding explosion has been detected for several core-collapse supernovae (SNe). Though the physical mechanisms driving the emission remain hotly debated, the light curves of detected events show long-lived ($\geq$50 days), plateau-like behavior, suggesting hydrogen recombination may significantly contribute to the total energy budget. The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will provide a decade-long photometric baseline to search for this emission, both in binned pre-explosion observations after an SN is detected and in single-visit observations prior to the SN explosion. In anticipation of these searches, we simulate a range of eruptive precursor models to core-collapse SNe and forecast the discovery rates of these phenomena in LSST data. We find a detection rate of ~40-130 yr$^{-1}$ for SN IIP/IIL precursors and ~110 yr$^{-1}$ for SN IIn precursors in single-epoch photometry. Considering the first three years of observations with the effects of rolling and observing triplets included, this number grows to a total of 150-400 in binned photometry, with the highest number recovered when binning in 100-day bins for 2020tlf-like precursors and in 20-day bins for other recombination-driven models from the literature. We quantify the impact of using templates contaminated by residual light (from either long-lived or separate precursor emission) on these detection rates, and explore strategies for estimating baseline flux to mitigate these issues. Spectroscopic follow-up of the eruptions preceding core-collapse SNe and detected with LSST will offer important clues to the underlying drivers of terminal-stage mass loss in massive stars. △ Less

Submitted 23 August, 2024; originally announced August 2024.

Comments: 20 pages, 12 figures. Comments welcome

One Year of SN 2023ixf: Breaking Through the Degenerate Parameter Space in Light-Curve Models with Pulsating Progenitors

Authors: Brian Hsu, Nathan Smith, Jared A. Goldberg, K. Azalee Bostroem, Griffin Hosseinzadeh, David J. Sand, Jeniveve Pearson, Daichi Hiramatsu, Jennifer E. Andrews, Emma R. Beasor, Yize Dong, Joseph Farah, LluÍs Galbany, Sebastian Gomez, Estefania Padilla Gonzalez, Claudia P. Gutiérrez, D. Andrew Howell, Réka Könyves-Tóth, Curtis McCully, Megan Newsome, Manisha Shrestha, Giacomo Terreran, V. Ashley Villar, Xiaofeng Wang

Abstract: We present and analyze the extensive optical broadband photometry of the Type II SN 2023ixf up to one year after explosion. We find that, when compared to two pre-existing model grids, the pseudo-bolometric light curve is consistent with drastically different combinations of progenitor and explosion properties. This may be an effect of known degeneracies in Type IIP light-curve models. We independ… ▽ More We present and analyze the extensive optical broadband photometry of the Type II SN 2023ixf up to one year after explosion. We find that, when compared to two pre-existing model grids, the pseudo-bolometric light curve is consistent with drastically different combinations of progenitor and explosion properties. This may be an effect of known degeneracies in Type IIP light-curve models. We independently compute a large grid of ${\tt MESA+STELLA}$ single-star progenitor and light-curve models with various zero-age main-sequence masses, mass-loss efficiencies, and convective efficiencies. Using the observed progenitor variability as an additional constraint, we select stellar models consistent with the pulsation period and explode them according to previously established scaling laws to match plateau properties. Our hydrodynamic modeling indicates that SN 2023ixf is most consistent with a moderate-energy ($E_{\rm exp}\approx7\times10^{50}$ erg) explosion of an initially high-mass red supergiant progenitor ($\gtrsim 17\ M_{\odot}$) that lost a significant amount of mass in its prior evolution, leaving a low-mass hydrogen envelope ($\lesssim 3\ M_{\odot}$) at the time of explosion, with a radius $\gtrsim 950\ R_{\odot}$ and a synthesized $^{56}$Ni mass of $0.07\ M_{\odot}$. We posit that previous mass transfer in a binary system may have stripped the envelope of SN 2023ixf's progenitor. The analysis method with pulsation period presented in this work offers a way to break degeneracies in light-curve modeling in the future, particularly with the upcoming Vera C.~Rubin Observatory Legacy Survey of Space and Time, when a record of progenitor variability will be more common. △ Less

Submitted 14 August, 2024; originally announced August 2024.

Comments: 18 pages, 7 figures, submitted to ApJ. Comments welcome

AI methods for approximate compiling of unitaries

Authors: David Kremer, Victor Villar, Sanjay Vishwakarma, Ismael Faro, Juan Cruz-Benito

Abstract: This paper explores artificial intelligence (AI) methods for the approximate compiling of unitaries, focusing on the use of fixed two-qubit gates and arbitrary single-qubit rotations typical in superconducting hardware. Our approach involves three main stages: identifying an initial template that approximates the target unitary, predicting initial parameters for this template, and refining these p… ▽ More This paper explores artificial intelligence (AI) methods for the approximate compiling of unitaries, focusing on the use of fixed two-qubit gates and arbitrary single-qubit rotations typical in superconducting hardware. Our approach involves three main stages: identifying an initial template that approximates the target unitary, predicting initial parameters for this template, and refining these parameters to maximize the fidelity of the circuit. We propose AI-driven approaches for the first two stages, with a deep learning model that suggests initial templates and an autoencoder-like model that suggests parameter values, which are refined through gradient descent to achieve the desired fidelity. We demonstrate the method on 2 and 3-qubit unitaries, showcasing promising improvements over exhaustive search and random parameter initialization. The results highlight the potential of AI to enhance the transpiling process, supporting more efficient quantum computations on current and future quantum hardware. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Off-axis Hartmann wavefront sensing for the GMT-Consortium Large Earth Finder (G-CLEF) red camera optics

Authors: Matthew C. H. Leung, Colby A. Jurgenson, Andrew Szentgyorgyi, Brian McLeod, Cem Onyuksel, Joseph Zajac, David Charbonneau, William Podgorski, Abigail Unger, Mark Mueller, Matthew Smith, Daniel Baldwin, V. Ashley Villar

Abstract: The Hartmann test is a method used to measure the wavefront error in a focal optical system, wherein a mask with a pattern of small holes is placed at the system's aperture stop. By taking an image at a defocused plane, the differences between the ideal and real positions of the reimaged holes (called the transverse ray aberrations) can be measured, which can then be used to estimate the wavefront… ▽ More The Hartmann test is a method used to measure the wavefront error in a focal optical system, wherein a mask with a pattern of small holes is placed at the system's aperture stop. By taking an image at a defocused plane, the differences between the ideal and real positions of the reimaged holes (called the transverse ray aberrations) can be measured, which can then be used to estimate the wavefront error. However, the Hartmann test is usually used with an on-axis field. In this paper, we present a wavefront sensing method which generalizes the classical Hartmann test for off-axis field angles and arbitrary reference wavefronts. Our method involves taking images at two defocused planes, and then using the real reimaged hole positions on both planes to estimate the trajectories of rays from the system's exit pupil, at which the reference wavefront is situated. We then propagate the rays forward from the reference wavefront to one of the two defocused planes, in order to find the ideal reimaged hole positions, from which we can compute transverse ray aberrations. We derive and solve a pair of nonlinear partial differential equations relating transverse ray aberrations to wavefront error, using Zernike decomposition and nonlinear least squares. Our method has been verified on simulated data from the 7-lens f/2.25 red camera system of the GMT-Consortium Large Earth Finder (G-CLEF), a high resolution optical echelle spectrograph which will be a first light instrument for the Giant Magellan Telescope (GMT). △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: 18 pages, 23 figures, SPIE Astronomical Telescopes + Instrumentation 2024

Journal ref: Proc. SPIE 13096, Ground-based and Airborne Instrumentation for Astronomy X, 130964M (18 July 2024)

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

Practical and efficient quantum circuit synthesis and transpiling with Reinforcement Learning

Authors: David Kremer, Victor Villar, Hanhee Paik, Ivan Duran, Ismael Faro, Juan Cruz-Benito

Abstract: This paper demonstrates the integration of Reinforcement Learning (RL) into quantum transpiling workflows, significantly enhancing the synthesis and routing of quantum circuits. By employing RL, we achieve near-optimal synthesis of Linear Function, Clifford, and Permutation circuits, up to 9, 11 and 65 qubits respectively, while being compatible with native device instruction sets and connectivity… ▽ More This paper demonstrates the integration of Reinforcement Learning (RL) into quantum transpiling workflows, significantly enhancing the synthesis and routing of quantum circuits. By employing RL, we achieve near-optimal synthesis of Linear Function, Clifford, and Permutation circuits, up to 9, 11 and 65 qubits respectively, while being compatible with native device instruction sets and connectivity constraints, and orders of magnitude faster than optimization methods such as SAT solvers. We also achieve significant reductions in two-qubit gate depth and count for circuit routing up to 133 qubits with respect to other routing heuristics such as SABRE. We find the method to be efficient enough to be useful in practice in typical quantum transpiling pipelines. Our results set the stage for further AI-powered enhancements of quantum computing workflows. △ Less

Submitted 21 May, 2024; originally announced May 2024.

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, D. O. Jones, C. -C. Lin, E. A. Magnier, K. S. Mandel, T. Matheson , et al. (7 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 24 July, 2024; v1 submitted 1 April, 2024; originally announced April 2024.

Comments: 44 pages (68 pages with Appendix), 15 figures, accepted 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 (composed of 412 unique ejecta parameter combinations evaluated at 54 viewing angles), 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 16 September, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: 16 pages, 11 figures

Report number: LA-UR-24-20638

Journal ref: Phys. Rev. Res. 6, 033078 (2024)

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

Gamma-ray Transient Network Science Analysis Group Report

Authors: Eric Burns, Michael Coughlin, Kendall Ackley, Igor Andreoni, Marie-Anne Bizouard, Floor Broekgaarden, Nelson L. Christensen, Filippo D'Ammando, James DeLaunay, Henrike Fleischhack, Raymond Frey, Chris L. Fryer, Adam Goldstein, Bruce Grossan, Rachel Hamburg, Dieter H. Hartmann, Anna Y. Q. Ho, Eric J. Howell, C. Michelle Hui, Leah Jenks, Alyson Joens, Stephen Lesage, Andrew J. Levan, Amy Lien, Athina Meli , et al. (12 additional authors not shown)

Abstract: The Interplanetary Network (IPN) is a detection, localization and alert system that utilizes the arrival time of transient signals in gamma-ray detectors on spacecraft separated by planetary baselines to geometrically locate the origin of these transients. Due to the changing astrophysical landscape and the new emphasis on time domain and multi-messenger astrophysics (TDAMM) from the Pathways to D… ▽ More The Interplanetary Network (IPN) is a detection, localization and alert system that utilizes the arrival time of transient signals in gamma-ray detectors on spacecraft separated by planetary baselines to geometrically locate the origin of these transients. Due to the changing astrophysical landscape and the new emphasis on time domain and multi-messenger astrophysics (TDAMM) from the Pathways to Discovery in Astronomy and Astrophysics for the 2020s, this Gamma-ray Transient Network Science Analysis Group was tasked to understand the role of the IPN and high-energy monitors in this new era. The charge includes describing the science made possible with these facilities, tracing the corresponding requirements and capabilities, and highlighting where improved operations of existing instruments and the IPN would enhance TDAMM science. While this study considers the full multiwavelength and multimessenger context, the findings are specific to space-based high-energy monitors. These facilities are important both for full characterization of these transients as well as facilitating follow-up observations through discovery and localization. The full document reports a brief history of this field, followed by our detailed analyses and findings in some 68 pages, providing a holistic overview of the role of the IPN and high-energy monitors in the coming decades. △ Less

Submitted 5 October, 2023; v1 submitted 8 August, 2023; originally announced August 2023.

Comments: Terms of Reference and additional information on the Science Analysis Group are available at https://pcos.gsfc.nasa.gov/sags/gtn-sag.php

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

The optical light curve of GRB 221009A: the afterglow and the emerging supernova

Authors: M. D. Fulton, S. J. Smartt, L. Rhodes, M. E. Huber, A. V. Villar, T. Moore, S. Srivastav, A. S. B. Schultz, K. C. Chambers, L. Izzo, J. Hjorth, T. -W. Chen, M. Nicholl, R. J. Foley, A. Rest, K. W. Smith, D. R. Young, S. A. Sim, J. Bright, Y. Zenati, T. de Boer, J. Bulger, J. Fairlamb, H. Gao, C. -C. Lin , et al. (24 additional authors not shown)

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

Submitted 23 March, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: Submitted to ApJL on 23rd January 2023, for consideration for publication in the special issue on GRB 221009A. Accepted on 2nd March 2023. The results of this paper are under press embargo until 28th March 2023. 15 pages, 7 figures, 1 table

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

Evidence for Extended Hydrogen-Poor CSM in the Three-Peaked Light Curve of Stripped Envelope Ib Supernova

Authors: Yossef Zenati, Qinan Wang, Alexey Bobrick, Lindsay DeMarchi, Hila Glanz, Mor Rozner, Armin Rest, Brian D. Metzger, Raffaella Margutti, Sebastian Gomez, Nathan Smith, Silvia Toonen, Joe S. Bright, Colin Norman, Ryan J. Foley, Alexander Gagliano, Julian H. Krolik, Stephen J. Smartt, Ashley V. Villar, Gautham Narayan, Ori Fox, Katie Auchettl, Daniel Brethauer, Alejandro Clocchiatti, Sophie V. Coelln , et al. (18 additional authors not shown)

Abstract: We present multi-band ATLAS photometry for SN 2019tsf, a stripped-envelope Type Ib supernova (SESN). The SN shows a triple-peaked light curve and a late (re-)brightening, making it unique among stripped-envelope systems. The re-brightening observations represent the latest photometric measurements of a multi-peaked Type Ib SN to date. As late-time photometry and spectroscopy suggest no hydrogen, t… ▽ More We present multi-band ATLAS photometry for SN 2019tsf, a stripped-envelope Type Ib supernova (SESN). The SN shows a triple-peaked light curve and a late (re-)brightening, making it unique among stripped-envelope systems. The re-brightening observations represent the latest photometric measurements of a multi-peaked Type Ib SN to date. As late-time photometry and spectroscopy suggest no hydrogen, the potential circumstellar material (CSM) must be H-poor. Moreover, late (>150 days) spectra show no signs of narrow emission lines, further disfavouring CSM interaction. On the contrary, an extended CSM structure is seen through a follow-up radio campaign with Karl G. Jansky Very Large Array (VLA), indicating a source of bright optically thick radio emission at late times, which is highly unusual among H-poor SESNe. We attribute this phenomenology to an interaction of the supernova ejecta with spherically-asymmetric CSM, potentially disk-like, and we present several models that can potentially explain the origin of this rare Type Ib supernova. The warped disc model paints a novel picture, where the tertiary companion perturbs the progenitors CSM, that can explain the multi-peaked light curves of SNe, and here we apply it to SN 2019tsf. This SN 2019tsf is likely a member of a new sub-class of Type Ib SNe and among the recently discovered class of SNe that undergo mass transfer at the moment of explosion △ Less

Submitted 23 July, 2022; v1 submitted 14 July, 2022; originally announced July 2022.

Comments: 23 pages, Comments are welcome, Submitted to ApJ

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