Search

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., Comment: 37 pages, 17 figures

Dust from core-collapse supernovae (CCSNe), specifically Type IIP SNe, has been suggested to be a significant source of the dust observed in high-redshift galaxies. CCSNe eject large amounts of newly formed heavy elements, which can condense into dust grains in the cooling ejecta. However, infrared (IR) observations of typical CCSNe generally measure dust masses that are too small to account for the dust production needed at high redshifts. Type IIn SNe, classified by their dense circumstellar medium (CSM), are also known to exhibit strong IR emission from warm dust, but the dust origin and heating mechanism have generally remained unconstrained because of limited observational capabilities in the mid-IR. Here, we present a JWST/MIRI Medium Resolution Spectrograph (MRS) spectrum of the Type IIn SN 2005ip nearly 17 years post-explosion. The Type IIn SN 2005ip is one of the longest-lasting and most well-studied SNe observed to date. Combined with a Spitzer mid-IR spectrum of SN 2005ip obtained in 2008, this data set provides a rare 15-year baseline, allowing for a unique investigation of the evolution of dust. The JWST spectrum shows a new high-mass dust component ($\gtrsim0.08$ M$_{\odot}$) that is not present in the earlier Spitzer spectrum. Our analysis shows dust likely formed over the past 15 years in the cold, dense shell (CDS), between the forward and reverse shocks. There is also a smaller mass of carbonaceous dust ($\gtrsim0.005$ M$_{\odot}$) in the ejecta. These observations provide new insights into the role of SN dust production, particularly within the CDS, and its potential contribution to the rapid dust enrichment of the early Universe.

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$, Comment: 19 pages, 8 figures, 1 table

We present a comprehensive, configurable open-source framework for estimating the rate of electromagnetic detection of kilonovae (KNe) associated with gravitational wave detections of binary neutron star (BNS) mergers. We simulate the current LIGO-Virgo-KAGRA (LVK) observing run (O4) using up-to-date sensitivity and up-time values as well as the next observing run (O5) using predicted sensitivities. We find the number of discoverable kilonovae during LVK O4 to be ${ 1}_{- 1}^{+ 4}$ or ${ 2 }_{- 2 }^{+ 3 }$, (at 90% confidence) depending on the distribution of NS masses in coalescing binaries, with the number increasing by an order of magnitude during O5 to ${ 19 }_{- 11 }^{+ 24 }$. Regardless of mass model, we predict at most five detectable KNe (at 95% confidence) in O4. We also produce optical and near-infrared light curves that correspond to the physical properties of each merging system. We have collated important information for allocating observing resources and directing search and follow-up observations including distributions of peak magnitudes in several broad bands and timescales for which specific facilities can detect each KN. The framework is easily adaptable, and new simulations can quickly be produced as input information such as merger rates and NS mass distributions are refined. Finally, we compare our suite of simulations to the thus-far completed portion of O4 (as of October 14, 2023), finding a median number of discoverable KNe of 0 and a 95-percentile upper limit of 2, consistent with no detection so far in O4., Comment: 16 pages, 13 figures, submitted to MNRAS

We present mid-infrared (mid-IR) imaging of the Type IIL supernova (SN) 1980K with the James Webb Space Telescope (JWST) more than 40 yr post-explosion. SN 1980K, located in the nearby ($D\approx7$ Mpc) "SN factory" galaxy NGC 6946, was serendipitously captured in JWST/MIRI images taken of the field of SN 2004et in the same galaxy. SN 1980K serves as a promising candidate for studying the transitional phase between young SNe and older SN remnants and also provides a great opportunity to investigate its the close environment. SN 1980K can be identified as a clear and bright point source in all eight MIRI filters from F560W up to F2550W. We fit analytical dust models to the mid-IR spectral energy distribution that reveal a large amount ($M_d \approx 0.002 {M}_{\odot}$) of Si-dominated dust at $T_{dust}\approx 150$ K (accompanied by a hotter dust/gas component), and also computed numerical SED dust models. Radiative transfer modeling of a late-time optical spectrum obtained recently with Keck discloses that an even larger ($\sim 0.24-0.58~{M}_{\odot}$) amount of dust is needed in order for selective extinction to explain the asymmetric line profile shapes observed in SN 1980K. As a conclusion, with JWST, we may see i) pre-existing circumstellar dust heated collisionally (or, partly radiatively), analogous to the equatorial ring of SN 1987A, or ii) the mid-IR component of the presumed newly-formed dust, accompanied by much more colder dust present in the ejecta (as suggested by the late-time the optical spectra)., Comment: 14 pages, 9 figures

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

We present optical, infrared, ultraviolet, and radio observations of SN 2022xkq, an underluminous fast-declining type Ia supernova (SN Ia) in NGC 1784 ($\mathrm{D}\approx31$ Mpc), from $<1$ to 180 days after explosion. The high-cadence observations of SN 2022xkq, a photometrically transitional and spectroscopically 91bg-like SN Ia, cover the first days and weeks following explosion which are critical to distinguishing between explosion scenarios. The early light curve of SN 2022xkq has a red early color and exhibits a flux excess which is more prominent in redder bands; this is the first time such a feature has been seen in a transitional/91bg-like SN Ia. We also present 92 optical and 19 near-infrared (NIR) spectra, beginning 0.4 days after explosion in the optical and 2.6 days after explosion in the NIR. SN 2022xkq exhibits a long-lived C I 1.0693 $\mu$m feature which persists until 5 days post-maximum. We also detect C II $\lambda$6580 in the pre-maximum optical spectra. These lines are evidence for unburnt carbon that is difficult to reconcile with the double detonation of a sub-Chandrasekhar mass white dwarf. No existing explosion model can fully explain the photometric and spectroscopic dataset of SN 2022xkq, but the considerable breadth of the observations is ideal for furthering our understanding of the processes which produce faint SNe Ia., Comment: 38 pages, 16 figures, accepted for publication in ApJ, the figure 15 input models and synthetic spectra are now available at https://zenodo.org/record/8379254

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., Comment: 35 pages, 17 figures, 7 tables, Accepted for Publication in ApJ

Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon-oxygen rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity ($M_{B}=-18.9$ mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak $B$-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O I] $\lambda\lambda 6300,\ 6364$ (${\rm FWHM} \approx 2{,}000$ km s$^{-1}$), strong, broad emission from [Ca II] $\lambda\lambda 7291,\ 7323$ (${\rm FWHM} \approx 7{,}300$ km s$^{-1}$), and a rapid Fe III to Fe II ionization change. Finally, we present the first-ever optical-to-mid-infrared (MIR) nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with $T \approx 500$ K), combined with prominent [O I] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within carbon/oxygen-rich CSM., Comment: 23 pages, 11 figures, submitted to ApJ

We present an analysis of ground-based and JWST observations of SN~2022pul, a peculiar "03fg-like" (or "super-Chandrasekhar") Type Ia supernova (SN Ia), in the nebular phase at 338d post explosion. Our combined spectrum continuously covers 0.4--14 $\mu$m and includes the first mid-infrared spectrum of an 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. Strong, broad, centrally peaked [Ne II] line at 12.81 $\mu$m was previously predicted as a hallmark of "violent merger'' SN Ia models, where dynamical interaction between two sub-$M_{ch}$ white dwarfs (WDs) causes disruption of the lower mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to better reproduce the optical iron emission, and add mass in the innermost region ($< 2000$ km s$^{-1}$) to account for the observed narrow [O I]~$\lambda\lambda6300$, 6364 emission. A violent WD-WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SN Ia., Comment: 20 pages, 10 figures, published in ApJ

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., Comment: 31 pages, 12 figures, accepted to MNRAS

We present nebular spectroscopy of SN 2020hvf, a Type Ia supernova (SN Ia) with an early bump in its light curve. SN 2020hvf shares many spectroscopic and photometric similarities to the carbon-rich high-luminosity "03fg-like" SNe Ia. At $>$240 days after peak brightness, we detect unambiguous emission from [Ca II] $\lambda\lambda$7291, 7324 which is never observed in normal-SNe Ia and only seen in peculiar subclasses. SN 2020hvf displays "saw-tooth" emission profiles near 7300 A that cannot be explained with single symmetric velocity components of [Fe II], [Ni II], and [Ca II], indicating an asymmetric explosion. The broad [Ca II] emission is best modeled by two velocity components offset by 1,220 km s$^{-1}$, which could be caused by ejecta associated with each star in the progenitor system, separated by their orbital velocity. For the first time in a SN Ia, we identify narrow (${\rm FWHM} = 180\pm40$ km s$^{-1}$) [Ca II] emission, which we associate with a wind from a surviving, puffed-up companion star. Few published spectra have sufficient resolution and signal-to-noise ratio necessary to detect similar narrow [Ca II] emission, however, we have detected similar line profiles in other 03fg-like SNe Ia. The extremely narrow velocity width of [Ca II] has only otherwise been observed in SNe Iax at late times. Since this event likely had a double-degenerate "super-Chandrasekhar" mass progenitor system, we suggest that a single white dwarf (WD) was fully disrupted and a wind from a surviving companion WD is producing the observed narrow emission. It is unclear if this unique progenitor and explosion scenario can explain the diversity of 03fg-like SNe Ia, potentially indicating that multiple progenitor channels contribute to this subclass., Comment: 18 pages, 10 figures, Published in ApJ

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 $\mu$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 $\mu$m and 4.5 $\mu$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 $\kappa$-mechanism pulsations observed in RSGs, albeit with a much larger amplitude than RSGs such as $\alpha$ Orionis (Betelgeuse)., Comment: 14 pages, 5 figures, submitted to ApJL, comments welcome

We present the optical photometric and spectroscopic analysis of two type Iax SNe 2018cni and 2020kyg. SN 2018cni is a bright type Iax SN (M$_{V,peak}$ = $-$17.81$\pm$0.21 mag) whereas SN 2020kyg (M$_{V,peak}$ = $-$14.52$\pm$0.21 mag) is a faint one. We derive $^{56}$Ni mass of 0.07 and 0.002 M${_\odot}$, ejecta mass of 0.48 and 0.14 M${_\odot}$ for SNe 2018cni and 2020kyg, respectively. A combined study of the bright and faint type Iax SNe in $R/r$- band reveals that the brighter objects tend to have a longer rise time. However, the correlation between the peak luminosity and decline rate shows that bright and faint type Iax SNe exhibit distinct behaviour. Comparison with standard deflagration models suggests that SN 2018cni is consistent with the deflagration of a CO white dwarf whereas the properties of SN 2020kyg can be better explained by the deflagration of a hybrid CONe white dwarf. The spectral features of both the SNe point to the presence of similar chemical species but with different mass fractions. Our spectral modelling indicates stratification at the outer layers and mixed inner ejecta for both the SNe., Comment: 18 pages, 18 figures, Accepted for Publication in The Astrophysical Journal

In late 2014, four images of Supernova (SN) "Refsdal," the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images' discovery, the SN was predicted to reappear within hundreds of days at a new position ~8 arcseconds away in the field. The observed reappearance in late 2015 makes it possible to carry out Refsdal's (1964) original proposal to use a multiply imaged SN to measure the Hubble constant H0, since the time delay between appearances should vary inversely with H0. Moreover, the position, brightness, and timing of the reappearance enable a novel test of the blind predictions of galaxy-cluster models, which are typically constrained only by the positions of multiply imaged galaxies. We have developed a new photometry pipeline that uses DOLPHOT to measure the fluxes of the five images of SN Refsdal from difference images. We apply four separate techniques to perform a blind measurement of the relative time delays and magnification ratios (mu_i/mu_1) between the last image SX and the earlier images S1-S4. We measure the relative time delay of SX-S1 to be 376.0+5.6-5.5 days and the relative magnification to be 0.30+0.05-0.03. This corresponds to a 1.5% precision on the time delay and 17% precision for the magnification ratios, and includes uncertainties due to millilensing and microlensing. In an accompanying paper, we place initial and blind constraints on the value of the Hubble constant., Comment: Published in ApJ. Companion paper presenting H0 constraints published in Science (DOI: 10.1126/science.abh1322)

The gravitationally lensed Supernova Refsdal appeared in multiple images, produced through gravitational lensing by a massive foreground galaxy cluster. After the supernova appeared in 2014, lens models of the galaxy cluster predicted an additional image of the supernova would appear in 2015, which was subsequently observed. We use the time delays between the images to perform a blinded measurement of the expansion rate of the Universe, quantified by the Hubble constant (H0). Using eight cluster lens models, we infer H0 = 64.8 +4.4-4.3 km / s / Mpc, where Mpc is the megaparsec. Using the two models most consistent with the observations, we find H0 = 66.6 +4.1-3.3 km / s / Mpc. The observations are best reproduced by models that assign dark-matter halos to individual galaxies and the overall cluster., Comment: Published in Science on May 11, 2023; this version updated to reflect minor edits to galley proofs. Companion paper presenting time-delay and relative magnification measurements published in ApJ (DOI: 10.3847/1538-4357/ac4ccb)

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., Comment: 17 pages, 13 figures, 1 table, submitted to MNRAS

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., Comment: 21 pages, 12 figures. Accepted by the astrophysical journal

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

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

We present reverberation mapping measurements for the prominent ultraviolet broad emission lines of the active galactic nucleus Mrk817 using 165 spectra obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our ultraviolet observations are accompanied by X-ray, optical, and near-infrared observations as part of the AGN Space Telescope and Optical Reverberation Mapping Program 2 (AGN STORM 2). Using the cross-correlation lag analysis method, we find significant correlated variations in the continuum and emission-line light curves. We measure rest-frame delayed responses between the far-ultraviolet continuum at 1180 A and Ly$\alpha$ $\lambda1215$ A ($10.4_{-1.4}^{+1.6}$ days), N V $\lambda1240$ A ($15.5_{-4.8}^{+1.0}$days), SiIV + OIV] $\lambda1397$ A ($8.2_{-1.4}^{+1.4}$ days), CIV $\lambda1549$ A ($11.8_{-2.8}^{+3.0}$ days), and HeII $\lambda1640$ A ($9.0_{-1.9}^{+4.5}$ days) using segments of the emission-line profile that are unaffected by absorption and blending, which results in sampling different velocity ranges for each line. However, we find that the emission-line responses to continuum variations are more complex than a simple smoothed, shifted, and scaled version of the continuum light curve. We also measure velocity-resolved lags for the Ly$\alpha$, and CIV emission lines. The lag profile in the blue wing of Ly$\alpha$ is consistent with virial motion, with longer lags dominating at lower velocities, and shorter lags at higher velocities. The CIV lag profile shows the signature of a thick rotating disk, with the shortest lags in the wings, local peaks at $\pm$ 1500 $\rm km\,s^{-1}$, and a local minimum at line center. The other emission lines are dominated by broad absorption lines and blending with adjacent emission lines. These require detailed models, and will be presented in future work., Comment: Submitted to ApJ. 25 pages, 8 figures, and 6 tables

Type Iax supernovae (SN Iax) are the largest known class of peculiar white dwarf supernovae, distinct from normal Type Ia supernovae (SN Ia). The unique properties of SN Iax, especially their strong photospheric lines out to extremely late times, allow us to model their optical spectra and derive physical parameters for the long-lasting photosphere. We present an extensive spectral timeseries, including 21 new spectra, of SN Iax 2014dt from +11 to +562 days after maximum light. We are able to reproduce the entire timeseries with a self-consistent, nearly unaltered deflagration explosion model from Fink et al. (2014) using TARDIS, an open-source radiative transfer code (Kerzendorf & Sim 2014; Kerzendorf et al. 2023). We find that the photospheric velocity of SN 2014dt slows its evolution between +64 and +148 days, which closely overlaps the phase when we see SN 2014dt diverge from the normal spectral evolution of SN Ia (+90 to +150 days). The photospheric velocity at these epochs, ~400$-$1000 km s$^{-1}$, may demarcate a boundary within the ejecta below which the physics of SN Iax and normal SN Ia differ. Our results suggest that SN 2014dt is consistent with a weak deflagration explosion model that leaves behind a bound remnant and drives an optically thick, quasi-steady-state wind creating the photospheric lines at late times. The data also suggest that this wind may weaken at epochs past +450 days, perhaps indicating a radioactive power source that has decayed away., Comment: Accepted to ApJ, 22 pages, 8 figures, 3 tables

Supernova (SN) explosions have been sought for decades as a possible source of dust in the Universe, providing the seeds of galaxies, stars, and planetary systems. SN 1987A offers one of the most promising examples of significant SN dust formation, but until the James Webb Space Telescope (JWST), instruments have traditionally lacked the sensitivity at both late times (>1 yr post-explosion) and longer wavelengths (i.e., >10 um) to detect analogous dust reservoirs. Here we present JWST/MIRI observations of two historic Type IIP SNe, 2004et and SN 2017eaw, at nearly 18 and 5 yr post-explosion, respectively. We fit the spectral energy distributions as functions of dust mass and temperature, from which we are able to constrain the dust geometry, origin, and heating mechanism. We place a 90% confidence lower limit on the dust masses for SNe 2004et and 2017eaw of >0.014 and >4e-4 M_sun, respectively. More dust may exist at even colder temperatures or may be obscured by high optical depths. We conclude dust formation in the ejecta to be the most plausible and consistent scenario. The observed dust is radiatively heated to ~100-150 K by ongoing shock interaction with the circumstellar medium. Regardless of the best fit or heating mechanism adopted, the inferred dust mass for SN 2004et is the second highest (next to SN 1987A) inferred dust mass in extragalactic SNe thus far, promoting the prospect of SNe as potential significant sources of dust in the Universe., Comment: 12 pages, 7 figures, submitting to MNRAS

A tight positive correlation between the stellar mass and the gas-phase metallicity of galaxies has been observed at low redshifts. The redshift evolution of this correlation can strongly constrain theories of galaxy evolution. The advent of JWST allows probing the mass-metallicity relation at redshifts far beyond what was previously accessible. Here we report the discovery of two emission-line galaxies at redshifts 8.15 and 8.16 in JWST NIRCam imaging and NIRSpec spectroscopy of targets gravitationally lensed by the cluster RXJ2129.4$+$0005. We measure their metallicities and stellar masses along with nine additional galaxies at $7.2 < z_{\rm spec} < 9.5$ to report the first quantitative statistical inference of the mass-metallicity relation at $z\approx8$. We measure $\sim 0.9$ dex evolution in the normalization of the mass-metallicity relation from $z \approx 8$ to the local Universe; at fixed stellar mass, galaxies are 8 times less metal enriched at $z \approx 8$ compared to the present day. Our inferred normalization is in agreement with the predictions of the FIRE simulations. Our inferred slope of the mass-metallicity relation is similar to or slightly shallower than that predicted by FIRE or observed at lower redshifts. We compare the $z \approx 8$ galaxies to extremely low metallicity analog candidates in the local Universe, finding that they are generally distinct from extreme emission-line galaxies or "green peas" but are similar in strong emission-line ratios and metallicities to "blueberry galaxies". Despite this similarity, at fixed stellar mass, the $z \approx 8$ galaxies have systematically lower metallicities compared to blueberry galaxies., Comment: Published in ApJ

We present observations of a peculiar hydrogen- and helium-poor stripped-envelope (SE) supernova (SN) 2020wnt, primarily in the optical and near-infrared (near-IR). Its peak absolute bolometric magnitude of -20.9 mag and a rise time of 69~days are reminiscent of hydrogen-poor superluminous SNe (SLSNe~I), luminous transients potentially powered by spinning-down magnetars. Before the main peak, there is a brief peak lasting <10 days post-explosion, likely caused by interaction with circumstellar medium (CSM) ejected ~years before the SN explosion. The optical spectra near peak lack a hot continuum and OII absorptions, which are signs of heating from a central engine; they quantitatively resemble those of radioactivity-powered H/He-poor Type Ic SESNe. At ~1 year after peak, nebular spectra reveal a blue pseudo-continuum and narrow OI recombination lines associated with magnetar heating. Radio observations rule out strong CSM interactions as the dominant energy source at +266 days post peak. Near-IR observations at +200-300 day reveal carbon monoxide and dust formation, which causes a dramatic optical light curve dip. Pair-instability explosion models predict slow light curve and spectral features incompatible with observations. SN 2020wnt is best explained as a magnetar-powered core-collapse explosion of a 28 Msun pre-SN star. The explosion kinetic energy is significantly larger than the magnetar energy at peak, effectively concealing the magnetar-heated inner ejecta until well after peak. SN 2020wnt falls into a continuum between normal SNe Ic and SLSNe I and demonstrates that optical spectra at peak alone cannot rule out the presence of a central engine., Comment: Submitted to ApJ

We present early observations and analysis of the double-peaked Type IIb supernova (SN IIb) 2021zby. $TESS$ captured the prominent early shock cooling peak of SN 2021zby within the first $\sim$10 days after explosion with a 30-minute cadence. We present optical and near-infrared spectral series of SN 2021zby, including three spectra during the shock cooling phase. Using a multi-band model fit, we find that the inferred properties of its progenitor are consistent with a red supergiant or yellow supergiant, with an envelope mass of $\sim$0.3-3.0 M$_\odot$ and an envelope radius of $\sim$50-350$ R_\odot$. These inferred progenitor properties are similar to those of other SNe IIb with double-peak feature, such as SNe 1993J, 2011dh, 2016gkg and 2017jgh. This study further validates the importance of the high cadence and early coverage in resolving the shape of the shock cooling light curve, while the multi-band observations, especially UV, is also necessary to fully constrain the progenitor properties., Comment: 12 pages, 5 figures, 2 tables, submitted to ApJL

Once only accessible in nearby galaxies, we can now study individual stars across much of the observable universe aided by galaxy-cluster gravitational lenses. When a star, compact object, or multiple such objects in the foreground galaxy-cluster lens become aligned, they can magnify a background individual star, and the timescale of a magnification peak can limit its size to tens of AU. The number and frequency of microlensing events therefore opens a window into the population of stars and compact objects, as well as high-redshift stars. To assemble the first statistical sample of stars in order to constrain the initial mass function (IMF) of massive stars at redshift z=0.7-1.5, the abundance of primordial black holes in galaxy-cluster dark matter, and the IMF of the stars making up the intracluster light, we are carrying out a 192-orbit program with the Hubble Space Telescope called "Flashlights," which is now two-thirds complete owing to scheduling challenges. We use the ultrawide F200LP and F350LP long-pass WFC3 UVIS filters and conduct two 16-orbit visits separated by one year. Having an identical roll angle during both visits, while difficult to schedule, yields extremely clean subtraction. Here we report the discovery of more than a dozen bright microlensing events, including multiple examples in the famous "Dragon Arc" discovered in the 1980s, as well as the "Spocks" and "Warhol" arcs that have hosted already known supergiants. The ultradeep observer-frame ultraviolet-through-optical imaging is sensitive to hot stars, which will complement deep James Webb Space Telescope infrared imaging. We are also acquiring Large Binocular Telescope LUCI and Keck-I MOSFIRE near-infrared spectra of the highly magnified arcs to constrain their recent star-formation histories.

We present near-infrared (NIR) and optical observations of the Type Ic supernova (SN Ic) SN 2021krf obtained between days 13 and 259 at several ground-based telescopes. The NIR spectrum at day 68 exhibits a rising $K$-band continuum flux density longward of $\sim$ 2.0 $\mu$m, and a late-time optical spectrum at day 259 shows strong [O I] 6300 and 6364 \r{A} emission-line asymmetry, both indicating the presence of dust, likely formed in the SN ejecta. We estimate a carbon-grain dust mass of $\sim$ 2 $\times$ 10$^{-5}$ M$_{\odot}$ and a dust temperature of $\sim$ 900 - 1200 K associated with this rising continuum and suggest the dust has formed in SN ejecta. Utilizing the one-dimensional multigroup radiation hydrodynamics code STELLA, we present two degenerate progenitor solutions for SN 2021krf, characterized by C-O star masses of 3.93 and 5.74 M$_{\odot}$, but with the same best-fit $^{56}$Ni mass of 0.11 M$_{\odot}$ for early times (0-70 days). At late times (70-300 days), optical light curves of SN 2021krf decline substantially more slowly than that expected from $^{56}$Co radioactive decay. Lack of H and He lines in the late-time SN spectrum suggests the absence of significant interaction of the ejecta with the circumstellar medium. We reproduce the entire bolometric light curve with a combination of radioactive decay and an additional powering source in the form of a central engine of a millisecond pulsar with a magnetic field smaller than that of a typical magnetar., Comment: Accepted for publication in ApJ, 27 pages, 21 figures, 6 tables. Previous arXiv submission (arXiv:2211.00205) replaced after acceptance

We present JWST near- and mid-infrared spectroscopic observations of the nearby normal Type Ia supernova SN 2021aefx in the nebular phase at $+255$ days past maximum light. Our Near Infrared Spectrograph (NIRSpec) and Mid Infrared Instrument (MIRI) observations, combined with ground-based optical data from the South African Large Telescope (SALT), constitute the first complete optical $+$ NIR $+$ MIR nebular SN Ia spectrum covering 0.3$-$14 $\mu$m. This spectrum unveils the previously unobserved 2.5$-$5 $\mu$m region, revealing strong nebular iron and stable nickel emission, indicative of high-density burning that can constrain the progenitor mass. The data show a significant improvement in sensitivity and resolution compared to previous Spitzer MIR data. We identify numerous NIR and MIR nebular emission lines from iron-group elements and as well as lines from the intermediate-mass element argon. The argon lines extend to higher velocities than the iron-group elements, suggesting stratified ejecta that are a hallmark of delayed-detonation or double-detonation SN Ia models. We present fits to simple geometric line profiles to features beyond 1.2 $\mu$m and find that most lines are consistent with Gaussian or spherical emission distributions, while the [Ar III] 8.99 $\mu$m line has a distinctively flat-topped profile indicating a thick spherical shell of emission. Using our line profile fits, we investigate the emissivity structure of SN 2021aefx and measure kinematic properties. Continued observations of SN 2021aefx and other SNe Ia with JWST will be transformative to the study of SN Ia composition, ionization structure, density, and temperature, and will provide important constraints on SN Ia progenitor and explosion models., Comment: published in ApJ Letters, 17 pages, 12 figures

Ultraviolet light from early galaxies is thought to have ionized gas in the intergalactic medium. However, there are few observational constraints on this epoch because of the faintness of those galaxies and the redshift of their optical light into the infrared. We report the observation, in JWST imaging, of a distant galaxy that is magnified by gravitational lensing. JWST spectroscopy of the galaxy, at rest-frame optical wavelengths, detects strong nebular emission lines that are attributable to oxygen and hydrogen. The measured redshift is z = 9.51 +- 0.01, corresponding to 510 million years after the Big Bang. The galaxy has a radius of 16.2+4.6-7.2 parsecs, which is substantially more compact than galaxies with equivalent luminosity at z = 6 to 8, leading to a high star formation rate surface density.

We present Young Supernova Experiment $grizy$ photometry of SN 2021hpr, the third Type Ia supernova sibling to explode in the Cepheid calibrator galaxy, NGC 3147. Siblings are useful for improving SN-host distance estimates, and investigating the contributions towards the SN Ia intrinsic scatter (post-standardisation residual scatter in distance estimates). We thus develop a principled Bayesian framework for analyzing SN Ia siblings. At its core is the cosmology-independent relative intrinsic scatter parameter, $\sigma_{Rel}$: the dispersion of siblings distance estimates relative to one another within a galaxy. It quantifies the contribution towards the total intrinsic scatter, $\sigma_0$, from within-galaxy variations about the siblings' common properties. It also affects the combined-distance uncertainty. We present analytic formulae for computing a $\sigma_{Rel}$-posterior from individual siblings distances (estimated using any SN-model). Applying a newly trained BayeSN model, we fit the light curves of each sibling in NGC 3147 individually, to yield consistent distance estimates. However, the wide $\sigma_{Rel}$-posterior means $\sigma_{Rel}\approx\sigma_0$ is not ruled out. We thus combine the distances by marginalizing over $\sigma_{Rel}$ with an informative prior: $\sigma_{Rel}\sim U(0,\sigma_0)$. Simultaneously fitting the trio's light curves improves constraints on distance, and each sibling's individual dust parameters, compared to individual fits. Higher correlation also tightens dust parameter constraints. Therefore, $\sigma_{Rel}$-marginalization yields robust estimates of siblings distances for cosmology, and dust parameters for siblings-host correlation studies. Incorporating NGC 3147's Cepheid-distance yields $H_0=78.4\pm 6.5\,$km/s/Mpc. Our work motivates analyses of homogeneous siblings samples, to constrain $\sigma_{Rel}$, and its SN-model dependence., Comment: Accepted for publication in ApJ. 22 pages, 11 figures, 8 tables. Single-galaxy siblings analysis

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, Comment: 23 pages, Comments are welcome, Submitted to ApJ

NASA will launch the Nancy Grace Roman Space Telescope (Roman) in the second half of this decade, which will allow for a generation-defining measurement of dark energy through multiple probes, including Type Ia supernovae (SNe Ia). To improve decisions on survey strategy, we have created the first simulations of realistic Roman images that include artificial SNe Ia injected as point sources in the images. Our analysis combines work done on Roman simulations for weak gravitational lensing studies as well as catalog-level simulations of SN samples. We have created a time series of images over two years containing $\sim$ 1,050 SNe Ia, covering a 1 square degree subarea of a planned 5 square degree deep survey. We have released these images publicly for community use along with input catalogs of all injected sources. We create secondary products from these images by generating coadded images and demonstrating recovery of transient sources using image subtraction. We perform first-use analyses on these images in order to measure galaxy-detection efficiency, point source-detection efficiency, and host-galaxy association biases. The simulated images can be found here: https://roman.ipac.caltech.edu/sims/SN_Survey_Image_sim.html., Comment: 12 pages, 12 figures. Submitted to MNRAS. For simulated images see https://roman.ipac.caltech.edu/sims/SN_Survey_Image_sim.html

Few published ultraviolet (UV) spectra exist for stripped-envelope supernovae, and none to date for broad-lined Type Ic supernovae (SN Ic-bl). These objects have extremely high ejecta velocities and are the only supernova type directly linked to gamma-ray bursts (GRBs). Here we present two epochs of HST/STIS spectra of the SN Ic-bl 2014ad, the first UV spectra for this class. We supplement this with 26 new epochs of ground-based optical spectra, augmenting a rich spectral time series. The UV spectra do not show strong features and are consistent with broadened versions of other SN Ic spectra observed in the UV. We measure Fe II 5169 Angstrom velocities and show that SN 2014ad has even higher ejecta velocities than most SNe Ic both with and without observed GRBs. We construct models of the SN 2014ad UV+optical spectra using TARDIS, a 1D Monte-Carlo radiative-transfer spectral synthesis code. The models fit the data well at multiple epochs in the optical but underestimate the flux in the UV, likely due to simplifying assumptions. We find that high densities at high velocities are needed to reproduce the spectra, with $\sim$3 M$_\odot$ of material at $v >$ 22,000 km s$^{-1}$, assuming spherical symmetry. Our nebular line fits suggest a steep density profile at low velocities. Together, these results imply a higher total ejecta mass than estimated from previous light curve analysis and expected from theory. This may be reconciled by a flattening of the density profile at low velocity and extra emission near the center of the ejecta., Comment: 26 pages, 14 figures, accepted in ApJ

We present high-cadence optical, ultraviolet (UV), and near-infrared data of the nearby ($D\approx23$ Mpc) Type II supernova (SN) 2021yja. Many Type II SNe show signs of interaction with circumstellar material (CSM) during the first few days after explosion, implying that their red supergiant (RSG) progenitors experience episodic or eruptive mass loss. However, because it is difficult to discover SNe early, the diversity of CSM configurations in RSGs has not been fully mapped. SN 2021yja, first detected within ${\approx}5.4$ hours of explosion, shows some signatures of CSM interaction (high UV luminosity, radio and x-ray emission) but without the narrow emission lines or early light curve peak that can accompany CSM. Here we analyze the densely sampled early light curve and spectral series of this nearby SN to infer the properties of its progenitor and CSM. We find that the most likely progenitor was an RSG with an extended envelope, encompassed by low-density CSM. We also present archival Hubble Space Telescope imaging of the host galaxy of SN 2021yja, which allows us to place a stringent upper limit of ${\lesssim}9\ M_\odot$ on the progenitor mass. However, this is in tension with some aspects of the SN evolution, which point to a more massive progenitor. Our analysis highlights the need to consider progenitor structure when making inferences about CSM properties, and that a comprehensive view of CSM tracers should be made to give a fuller view of the last years of RSG evolution., Comment: updated to match accepted version

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 $\delta 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., Comment: 33 pages, 19 figures. Submitted to ApJ. Comments welcome!

Located on the campus of the Thacher School in Southern California, the Thacher Observatory has a legacy of astronomy research and education that dates back to the late 1950's. In 2016, the observatory was fully renovated with upgrades including a new 0.7-m telescope, a research grade camera, and a slit dome with full automation capabilities. The low-elevation site is bordered by the Los Padres National Forest and therefore affords dark to very dark skies allowing for accurate and precise photometric observations. We present a characterization of the site including sky brightness, weather, and seeing, and we demonstrate the on-sky performance of the facility. Our primary research programs are based around our multi-band photometric capabilities and include photometric monitoring of variable sources, a nearby supernova search and followup program, a quick response transient followup effort, and exoplanet and eclipsing binary light curves. Select results from these programs are included in this work which highlight the broad range of science available to an automated observatory with a moderately sized telescope., Comment: 20 pages, 18 figures, PASP accepted

We present constraints on cosmological parameters from the Pantheon+ analysis of 1701 light curves of 1550 distinct Type Ia supernovae (SNe Ia) ranging in redshift from $z=0.001$ to 2.26. This work features an increased sample size, increased redshift span, and improved treatment of systematic uncertainties in comparison to the original Pantheon analysis and results in a factor of two improvement in cosmological constraining power. For a Flat$\Lambda$CDM model, we find $\Omega_M=0.334\pm0.018$ from SNe Ia alone. For a Flat$w_0$CDM model, we measure $w_0=-0.90\pm0.14$ from SNe Ia alone, H$_0=73.5\pm1.1$ km s$^{-1}$ Mpc$^{-1}$ when including the Cepheid host distances and covariance (SH0ES), and $w_0=-0.978^{+0.024}_{-0.031}$ when combining the SN likelihood with constraints from the cosmic microwave background (CMB) and baryon acoustic oscillations (BAO); both $w_0$ values are consistent with a cosmological constant. We also present the most precise measurements to date on the evolution of dark energy in a Flat$w_0w_a$CDM universe, and measure $w_a=-0.1^{+0.9}_{-2.0}$ from Pantheon+ alone, H$_0=73.3\pm1.1$ km s$^{-1}$ Mpc$^{-1}$ when including SH0ES, and $w_a=-0.65^{+0.28}_{-0.32}$ when combining Pantheon+ with CMB and BAO data. Finally, we find that systematic uncertainties in the use of SNe Ia along the distance ladder comprise less than one third of the total uncertainty in the measurement of H$_0$ and cannot explain the present "Hubble tension" between local measurements and early-Universe predictions from the cosmological model., Comment: 34 Pages, 16 Figures, 7 Tables. Published in ApJ. Comments welcome. Papers and data release here: https://pantheonplussh0es.github.io

Seeing pristine material from the donor star in a Type Ia supernova (SN Ia) explosion can reveal the nature of the binary system. In this paper, we present photometric and spectroscopic observations of SN 2020esm, one of the best-studied SNe of the class of "super-Chandrasekhar" SNe Ia (SC SNe Ia), with data obtained $-12$ to +360 days relative to peak brightness, obtained from a variety of ground- and space-based telescopes. Initially misclassified as a Type II supernova, SN 2020esm peaked at $M_{B} = -19.9$ mag, declined slowly ($\Delta m_{15}(B) = 0.92$ mag), and had particularly blue UV and optical colors at early times. Photometrically and spectroscopically, SN 2020esm evolved similarly to other SC SNe Ia, showing the usual low ejecta velocities, weak intermediate mass elements (IMEs), and the enhanced fading at late times, but its early spectra are unique. Our first few spectra (corresponding to a phase of $\gtrsim$10~days before peak) reveal a nearly-pure carbon/oxygen atmosphere during the first days after explosion. This composition can only be produced by pristine material, relatively unaffected by nuclear burning. The lack of H and He may further indicate that SN 2020esm is the outcome of the merger of two carbon/oxygen white dwarfs (WDs). Modeling its bolometric light curve, we find a $^{56}$Ni mass of $1.23^{+0.14}_{-0.14}$ M$_{\odot}$ and an ejecta mass of $1.75^{+0.32}_{-0.20}$ M$_{\odot}$, in excess of the Chandrasekhar mass. Finally, we discuss possible progenitor systems and explosion mechanisms of SN 2020esm and, in general, the SC SNe Ia class., Comment: 21 pages, 15 figures, accepted for publication in APJ

Here we present 1701 light curves of 1550 spectroscopically confirmed Type Ia supernovae (SNe Ia) that will be used to infer cosmological parameters as part of the Pantheon+ SN analysis and the SH0ES (Supernovae and H0 for the Equation of State of dark energy) distance-ladder analysis. This effort is one part of a series of works that perform an extensive review of redshifts, peculiar velocities, photometric calibration, and intrinsic-scatter models of SNe Ia. The total number of light curves, which are compiled across 18 different surveys, is a significant increase from the first Pantheon analysis (1048 SNe), particularly at low redshift ($z$). Furthermore, unlike in the Pantheon analysis, we include light curves for SNe with $z<0.01$ such that SN systematic covariance can be included in a joint measurement of the Hubble constant (H$_0$) and the dark energy equation-of-state parameter ($w$). We use the large sample to compare properties of 151 SNe Ia observed by multiple surveys and 12 pairs/triplets of "SN siblings" - SNe found in the same host galaxy. Distance measurements, application of bias corrections, and inference of cosmological parameters are discussed in the companion paper by Brout et al. (2022b), and the determination of H$_0$ is discussed by Riess et al. (2022). These analyses will measure w with $\sim3\%$ precision and H$_0$ with 1 km/s/Mpc precision., Comment: Submitted to ApJL. Comments welcome. Papers and data release here: https://github.com/PantheonPlusSH0ES/PantheonPlusSH0ES.github.io

We present Hubble Space Telescope (HST) observations of the type IIb supernova (SN) 2016gkg at 652, 1698, and 1795 days with the Advanced Camera for Surveys (ACS) and Wide Field Camera 3 (WFC3). Comparing to pre-explosion imaging from 2001 obtained with the Wide Field Planetary Camera 2, we demonstrate that SN 2016gkg is now fainter than its candidate counterpart in the latest WFC3 imaging, implying that the counterpart has disappeared and confirming that it was the SN progenitor star. We show the latest light curve and Keck spectra of SN 2016gkg, which implies that SN 2016gkg is declining more slowly than the expected rate for ${}^{56}$Co decay during its nebular phase. We find that this emission is too luminous to be powered by other radioisotopes, thus we infer that SN 2016gkg is entering a new phase in its evolution where it is powered primarily by interaction with circumstellar matter. Finally, we re-analyze the progenitor star spectral energy distribution and late-time limits in the context of binary evolution models and including emission from a potential companion star and find that all companion stars would be fainter than our limiting magnitudes., Comment: 10 pages, 5 figures

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

Detailed spectropolarimetric studies may hold the key to probing the explosion mechanisms and the progenitor scenarios of Type Ia supernovae (SNe Ia). We present multi-epoch spectropolarimetry and imaging polarimetry of SN 2019ein, an SN Ia showing high expansion velocities at early phases. The spectropolarimetry sequence spans from $\sim -11$ to $+$10 days relative to peak brightness in the $B$-band. We find that the level of the continuum polarization of SN 2019ein, after subtracting estimated interstellar polarization, is in the range $0.0-0.3\%$, typical for SNe Ia. The polarization position angle remains roughly constant before and after the SN light-curve peak, implying that the inner regions share the same axisymmetry as the outer layers. We observe high polarization ($\sim 1\%$) across both the Si II $\lambda6355$ and Ca II near-infrared triplet features. These two lines also display complex polarization modulations. The spectropolarimetric properties of SN 2019ein rule out a significant departure from spherical symmetry of the ejecta for up to a month after the explosion. These observations disfavour merger-induced and double-detonation models for SN 2019ein. The imaging polarimetry shows weak evidence for a modest increase in polarization after $\sim 20$ days since the $B$-band maximum. If this rise is real and is observed in other SNe Ia at similar phases, we may have seen, for the first time, an aspherical interior similar to what has been previously observed for SNe IIP. Future polarization observations of SNe Ia extending to post-peak epochs will help to examine the inner structure of the explosion., Comment: Accepted by MNRAS. 13 pages, 7 figures. Updated Figure 6

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

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

We test the adequacy of ultraviolet (UV) spectra for characterizing the outer structure of Type Ia supernova (SN) ejecta. For this purpose, we perform spectroscopic analysis for ASASSN-14lp, a normal SN Ia showing low continuum in the mid-UV regime. To explain the strong UV suppression, two possible origins have been investigated by mapping the chemical profiles over a significant part of their ejecta. We fit the spectral time series with mid-UV coverage obtained before and around maximum light by HST, supplemented with ground-based optical observations for the earliest epochs. The synthetic spectra are calculated with the one dimensional MC radiative-transfer code TARDIS from self-consistent ejecta models. Among several physical parameters, we constrain the abundance profiles of nine chemical elements. We find that a distribution of $^{56}$Ni (and other iron-group elements) that extends toward the highest velocities reproduces the observed UV flux well. The presence of radioactive material in the outer layers of the ejecta, if confirmed, implies strong constraints on the possible explosion scenarios. We investigate the impact of the inferred $^{56}$Ni distribution on the early light curves with the radiative transfer code TURTLS, and confront the results with the observed light curves of ASASSN-14lp. The inferred abundances are not in conflict with the observed photometry. We also test whether the UV suppression can be reproduced if the radiation at the photosphere is significantly lower in the UV regime than the pure Planck function. In this case, solar metallicity might be sufficient enough at the highest velocities to reproduce the UV suppression., Comment: 18 pages, 16 figures, accepted for publication in MNRAS

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for ''red'' kilonovae and rule out ''blue'' kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources., Comment: 86 pages, 9 figures

Type Iax supernovae represent the largest class of peculiar white-dwarf supernovae. The type Iax SN~2012Z in NGC 1309 is the only white dwarf supernova with a detected progenitor system in pre-explosion observations. Deep \textit{Hubble Space Telescope} images taken before SN~2012Z show a luminous, blue source that we have interpreted as a helium-star companion (donor) to the exploding white dwarf. We present here late-time \textit{HST} observations taken $\sim$1400 days after the explosion to test this model. We find the SN light curve can empirically be fit by an exponential decay model in magnitude units. The fitted asymptotic brightness is within $10\%$ of our latest measurements and approximately twice the brightness of the pre-explosion source. The decline of the light curve is too slow to be powered by $^{56}$Co or $^{57}$Co decay: if radioactive decay is the dominate power source, it must be from longer half-life species like $^{55}$Fe. Interaction with circumstellar material may contribute to the light curve, as may shock heating of the companion star. Companion-star models underpredict the observed flux in the optical, producing most of their flux in the UV at these epochs. A radioactively-heated bound remnant, left after only a partial disruption of the white dwarf, is also capable of producing the observed excess late-time flux. Our analysis suggests that the total ejecta + remnant mass is consistent with the Chandrasekhar mass for a range of type Iax supernovae., Comment: 20 pages, 13 figures, accepted to ApJ

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

We present X-ray and radio observations of what may be the closest type Iax supernova (SN) to date, SN 2014dt (d=12.3-19.3 Mpc) and provide tight constraints on the radio and X-ray emission. We infer a specific radio luminosity of < (1.0-2.4)E25 erg/s/Hz at a frequency of 7.5 GHz and a X-ray luminosity < 1.4E38 erg/s (0.3-10 keV) at ~38-48 days post-explosion. We interpret these limits in the context of Inverse Compton (IC) emission and synchrotron emission from a population of electrons accelerated at the forward shock of the explosion in a power-law distribution $N_e(\gamma_e)\propto \gamma_e^{-p}$ with p=3. Our analysis constrains the progenitor system mass-loss rate to be smaller than 5E-6 solar masses per year at distances where r <= 1E16 cm for an assumed wind velocity v=100 km/s, and a fraction of post-shock energy into magnetic fields and relativistic electrons of epsilon_B=0.01 and epsilon_e=0.1, respectively. This result rules out some of the parameter space of symbiotic giant star companions, and it is consistent with the low mass-loss rates expected from He-star companions. Our calculations also show that the improved sensitivity of the next generation Very Large Array (ngVLA) is needed to probe the very low-density media characteristic of He stars that are the leading model for binary stellar companions of white dwarfs giving origin to type Iax SNe., Comment: Submitted to MNRAS, 9 pages and 7 figures