Search

The recently discovered Pa 30 nebula, the putative type Iax supernova remnant associated with the historical supernova of 1181 AD, shows puzzling characteristics that make it unique among known supernova remnants. In particular, Pa 30 exhibits a complex morphology, with a unique radial and filamentary structure, and it hosts a hot stellar remnant at its center, which displays oxygen-dominated, ultra-fast winds. Because of the surviving stellar remnant and the lack of hydrogen and helium in its filaments, it has been suggested that Pa 30 is the product of a failed thermonuclear explosion in a near- or super-Chandrasekhar white dwarf, which created a sub-luminous transient, a rare sub-type of the Ia class of supernovae called type Iax. We here present a detailed study of the 3D structure and velocities of a full radial section of the remnant. The Integral Field Unit (IFU) observations, obtained with the new red channel of the Keck Cosmic Web Imager spectrograph, reveal that the ejecta are consistent with being ballistic, with velocities close to the free-expansion velocity. Additionally, we detect a large cavity inside the supernova remnant and a sharp inner edge to the filamentary structure, which coincides with the outer edge of a bright ring detected in infrared images. Finally, we detect a strong asymmetry in the amount of ejecta along the line of sight, which might hint to an asymmetric explosion. Our analysis provides strong confirmation that the explosion originated from SN 1181., Comment: 13 pages, 8 figures. Accepted for publication in ApJL

Understanding the nature of the luminous 1991T-like supernovae is of great importance to supernova cosmology as they are likely to have been more common in the early universe. In this paper we explore the observational properties of 1991T-like supernovae to study their relationship to other luminous, slow-declining Type~Ia supernovae (SNe Ia). From the spectroscopic and photometric criteria defined in Phillips et al. (1992), we identify 17 1991T-like supernovae from the literature. Combining these objects with ten 1991T-like supernovae from the Carnegie Supernova Project-II, the spectra, light curves, and colors of these events, along with their host galaxy properties, are examined in detail. We conclude that 1991T-like supernovae are closely related in essentially all of their UV, optical, and near-infrared properties -- as well as their host galaxy parameters -- to the slow-declining subset of Branch core-normal supernovae and to the intermediate 1999aa-like events, forming a continuum of luminous SNe Ia. The overriding difference between these three subgroups appears to be the extent to which $^{56}$Ni mixes into the ejecta, producing the pre-maximum spectra dominated by Fe III absorption, the broader UV light curves, and the higher luminosities that characterize the 1991T-like events. Nevertheless, the association of 1991T-like SNe with the rare Type Ia CSM supernovae would seem to run counter to this hypothesis, in which case 1991T-like events may form a separate subclass of SNe Ia, possibly arising from single-degenerate progenitor systems., Comment: Accepted for publication in ApJS

Modern transient surveys now routinely discover flares resulting from tidal disruption events (TDEs) which occur when stars, typically $\sim0.5-2$ M$_{\odot}$, are ripped apart after passing too close to a supermassive black hole. We present three examples of a new class of extreme nuclear transients (ENTs) that we interpret as the tidal disruption of intermediate mass ($\sim3-10$ M$_{\odot}$) stars. Each is coincident with their host-galaxy nucleus and exhibits a smooth ($<10$% excess variability), luminous ($2-7\times10^{45}$ erg s$^{-1}$), and long-lived ($>150$ days) flare. ENTs are extremely rare ($\geq1\times10^{-3}$ Gpc$^{-1}$ yr$^{-1}$) compared to any other known class of transients. They are at least twice as energetic ($0.5-2.5\times 10^{53}$ erg) as any other known transient and these extreme energetics rule out stellar origins., Comment: Submitted to Science

We demonstrate in a proof-of-concept numerical hydrodynamics calculation that the narrow radial filamentary structures seen in Pa 30 could be generated through highly efficient cooling (e.g. via line emission) in the ejecta. Efficient cooling in the ejecta causes a drop of pressure support in Rayleigh-Taylor fingers, leading them to be compressed, and suppressing the growth of Kelvin-Helmholtz instability. Following this result, we make three predictions that could determine whether this is the mechanism responsible for shaping Pa 30: First, we predict very strong emission lines, strong enough to cool a significant fraction of the shock energy in an expansion time. Secondly, we predict that the forward shock should be highly corrugated on small scales, with the shock front closely following the structure of the filaments. Third, we predict that these filaments should be nearly ballistic, with velocities around 90% of the free-expansion velocity ($v \approx 0.9 ~r/t$). These predictions should be falsifiable in follow-up observations of this remnant., Comment: Rejected from both ApJ Letters and MNRAS letters. I'm giving up on this one; I'll just leave the updated version here while we work on the follow-up. Please let me know if you have any comments

Supernova remnants (SNRs) exhibit varying degrees of anisotropy, which have been extensively modeled using numerical methods. We implement a technique to measure anisotropies in SNRs by calculating power spectra from their high-resolution images. To test this technique, we develop 3D hydrodynamical models of supernova remnants and generate synthetic x-ray images from them. Power spectra extracted from both the 3D models and the synthetic images exhibit the same dominant angular scale, which separates large scale features from small scale features due to hydrodynamic instabilities. The angular power spectrum at small length scales during relatively early times is too steep to be consistent with Kolmogorov turbulence, but it transitions to Kolmogorov turbulence at late times. As an example of how this technique can be applied to observations, we extract a power spectrum from a \textit{Chandra} observation of Tycho's SNR and compare with our models. Our predicted power spectrum picks out the angular scale of Tycho's fleece-like structures and also agrees with the small-scale power seen in Tycho. We use this to extract an estimate for the density of the circumstellar gas ($n \sim 0.28/\mathrm{cm^3}$), consistent with previous measurements of this density by other means. The power spectrum also provides an estimate of the density profile of the outermost ejecta. Moreover, we observe additional power at large scales which may provide important clues about the explosion mechanism itself., Comment: 13 pages, 8 figures, submitted to ApJ. Comments are most welcome

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 optical and near-infrared observations of SN~2022crv, a stripped envelope supernova in NGC~3054, discovered within 12 hrs of explosion by the Distance Less Than 40 Mpc Survey. We suggest SN~2022crv is a transitional object on the continuum between SNe Ib and SNe IIb. A high-velocity hydrogen feature ($\sim$$-$20,000 -- $-$16,000 $\rm km\,s^{-1}$) was conspicuous in SN~2022crv at early phases, and then quickly disappeared around maximum light. By comparing with hydrodynamic modeling, we find that a hydrogen envelope of $\sim 10^{-3}$ \msun{} can reproduce the behaviour of the hydrogen feature observed in SN~2022crv. The early light curve of SN~2022crv did not show envelope cooling emission, implying that SN~2022crv had a compact progenitor with extremely low amount of hydrogen. The analysis of the nebular spectra shows that SN~2022crv is consistent with the explosion of a He star with a final mass of $\sim$4.5 -- 5.6 \msun{} that has evolved from a $\sim$16 -- 22 \msun{} zero-age main sequence star in a binary system with about 1.0 -- 1.7 \msun{} of oxygen finally synthesized in the core. The high metallicity at the supernova site indicates that the progenitor experienced a strong stellar wind mass loss. In order to retain a small amount of residual hydrogen at such a high metallicity, the initial orbital separation of the binary system is likely larger than $\sim$1000~$\rm R_{\odot}$. The near-infrared spectra of SN~2022crv show a unique absorption feature on the blue side of He I line at $\sim$1.005~$\mu$m. This is the first time that such a feature has been observed in a Type Ib/IIb, and could be due to \ion{Sr}{2}. Further detailed modelling on SN~2022crv can shed light on the progenitor and the origin of the mysterious absorption feature in the near infrared., Comment: accepted for publication in ApJ

We study the coupling of hydrodynamics and reactions in simulations of the double detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator / Strang splitting vs.\ the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and the nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation timestep to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code., Comment: accepted to Astrophysical Journal this version has a new appendix (on shock burning) compared to previous in response to the referee

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 develop a suite of 3D hydrodynamic models of supernova remnants (SNRs) expanding against the circumstellar medium (CSM). We study the Rayleigh-Taylor Instability (RTI) forming at the expansion interface by calculating an angular power spectrum for each of these models. The power spectra of young SNRs is seen to exhibit a dominant angular mode, which is a diagnostic of their ejecta density profile as found by previous studies. The steep scaling of power at smaller modes and the time evolution of the spectra is indicative of absence of a turbulent cascade. Instead, as the time evolution of the spectra suggests, they may be governed by an angular mode dependent net growth rate. We also study the impact of anisotropies in the ejecta as well as in the CSM on the power spectra of velocity and density. We confirm that perturbations in the density field (whether imposed on the ejecta or the CSM) do not influence the anisotropy of the remnant significantly unless they have a very large amplitude and form large-scale coherent structures. In any case, these clumps can only affect structures on large angular scales. The power spectra on small angular scales is completely independent of the initial clumpiness and only governed by the growth and saturation of the Rayleigh-Taylor instability., Comment: 13 pages, 9 figures; accepted for publication in ApJ

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

In recent years, there has been ample evidence that Type Ia supernova (SNe Ia) with high Si 2 velocities near peak brightness are distinguished from SNe Ia of lower velocities and may indeed represent a separate progenitor system. These SNe Ia can contaminate the population of normal events used for cosmological analyses, creating unwanted biases in the final analyses. Given that many current and future surveys using SNe Ia as cosmological probes will not have the resources to take a spectrum of all the events, likely only getting host redshifts long after the SNe Ia have faded, we need to turn to methods that could separate these populations based purely on photometry or host properties. Here, we present a study of a sample of well observed, nearby SNe Ia and their hosts to determine if there are significant enough differences between these populations that can be discerned only from the stellar population properties of their hosts. Our results indicate that the global host properties, including star formation, stellar mass, stellar population age, and dust attenuation, of high velocity SNe Ia do not differ significantly from those of lower velocities. However, we do find that high velocity SNe Ia are more concentrated toward the center of their hosts, suggesting that their local environments may indeed differ. Future work requires strengthening photometric probes of high velocity SNe Ia and their local environments to distinguish these events and determine if they originate from a separate progenitor., Comment: 28 pages, 10 figures, 5 tables, submitted

We present photometric and spectroscopic observations and analysis of SN 2021bxu (ATLAS21dov), a low-luminosity, fast-evolving Type IIb supernova (SN). SN 2021bxu is unique, showing a large initial decline in brightness followed by a short plateau phase. With $M_r = -15.93 \pm 0.16\, \mathrm{mag}$ during the plateau, it is at the lower end of the luminosity distribution of stripped-envelope supernovae (SE-SNe) and shows a distinct $\sim$10 day plateau not caused by H- or He-recombination. SN 2021bxu shows line velocities which are at least $\sim1500\,\mathrm{km\,s^{-1}}$ slower than typical SE-SNe. It is photometrically and spectroscopically similar to Type IIb SNe during the photospheric phases of evolution, with similarities to Ca-rich IIb SNe. We find that the bolometric light curve is best described by a composite model of shock interaction between the ejecta and an envelope of extended material, combined with a typical SN IIb powered by the radioactive decay of $^{56}$Ni. The best-fit parameters for SN 2021bxu include a $^{56}$Ni mass of $M_{\mathrm{Ni}} = 0.029^{+0.004}_{-0.005}\,\mathrm{M_{\odot}}$, an ejecta mass of $M_{\mathrm{ej}} = 0.61^{+0.06}_{-0.05}\,\mathrm{M_{\odot}}$, and an ejecta kinetic energy of $K_{\mathrm{ej}} = 8.8^{+1.1}_{-1.0} \times 10^{49}\, \mathrm{erg}$. From the fits to the properties of the extended material of Ca-rich IIb SNe we find a trend of decreasing envelope radius with increasing envelope mass. SN 2021bxu has $M_{\mathrm{Ni}}$ on the low end compared to SE-SNe and Ca-rich SNe in the literature, demonstrating that SN 2021bxu-like events are rare explosions in extreme areas of parameter space. The progenitor of SN 2021bxu is likely a low mass He star with an extended envelope., Comment: 18 pages, 15 figures, accepted to MNRAS

We present a photometric and spectroscopic analysis of the ultra-luminous and slowly evolving 03fg-like Type Ia SN 2021zny. Our observational campaign starts from $\sim5.3$ hours after explosion (making SN 2021zny one of the earliest observed members of its class), with dense multi-wavelength coverage from a variety of ground- and space-based telescopes, and is concluded with a nebular spectrum $\sim10$ months after peak brightness. SN 2021zny displayed several characteristics of its class, such as the peak brightness ($M_{B}=-19.95$ mag), the slow decline ($\Delta m_{15}(B) = 0.62$ mag), the blue early-time colours, the low ejecta velocities and the presence of significant unburned material above the photosphere. However, a flux excess for the first $\sim1.5$ days after explosion is observed in four photometric bands, making SN 2021zny the third 03fg-like event with this distinct behavior, while its $+313$ d spectrum shows prominent [O I] lines, a very unusual characteristic of thermonuclear SNe. The early flux excess can be explained as the outcome of the interaction of the ejecta with $\sim0.04\:\mathrm{M_{\odot}}$ of H/He-poor circumstellar material at a distance of $\sim10^{12}$ cm, while the low ionization state of the late-time spectrum reveals low abundances of stable iron-peak elements. All our observations are in accordance with a progenitor system of two carbon/oxygen white dwarfs that undergo a merger event, with the disrupted white dwarf ejecting carbon-rich circumstellar material prior to the primary white dwarf detonation., Comment: 19 pages, 16 figures, accepted for publication in MNRAS

1991T-like supernovae are the luminous, slow-declining extreme of the Branch shallow-silicon (SS) subclass of Type Ia supernovae. They are distinguished by extremely weak Ca II H & K and Si II $\lambda6355$ and strong Fe III absorption features in their optical spectra at pre-maximum phases, and have long been suspected to be over-luminous compared to normal Type Ia supernovae. In this paper, the pseudo equivalent width of the Si II $\lambda$6355 absorption obtained at light curve phases from $\leq+10$ days is combined with the morphology of the $i$-band light curve to identify a sample of 1991T-like supernovae in the Carnegie Supernova Project-II. Hubble diagram residuals show that, at optical as well as near-infrared wavelengths, these events are over-luminous by $\sim$0.1-0.5 mag with respect to the less extreme Branch SS (1999aa-like) and Branch core-normal supernovae with similar $B$-band light curve decline rates., Comment: Submitted to ApJ (Aug 30, 2022)

The detonation of a thin ($\lesssim$$0.03\,\mathrm{M_\odot}$) helium shell (He-shell) atop a $\sim$$1\,\mathrm{M_\odot}$ white dwarf (WD) is a promising mechanism to explain normal Type Ia supernovae (SNe Ia), while thicker He-shells and less massive WDs may explain some recently observed peculiar SNe Ia. We present observations of SN 2020jgb, a peculiar SN Ia discovered by the Zwicky Transient Facility (ZTF). Near maximum light, SN 2020jgb is slightly subluminous (ZTF $g$-band absolute magnitude $M_g$ between $-18.2$ and $-18.7$ mag depending on the amount of host galaxy extinction) and shows an unusually red color ($g_\mathrm{ZTF}-r_\mathrm{ZTF}$ between 0.4 and 0.2 mag) due to strong line-blanketing blueward of $\sim$5000 $\r{A}$. These properties resemble those of SN 2018byg, a peculiar SN Ia consistent with a thick He-shell double detonation (DDet) SN. Using detailed radiative transfer models, we show that the optical spectroscopic and photometric evolution of SN 2020jgb are broadly consistent with a $\sim$$0.95\,\mathrm{M_\odot}$ (C/O core + He-shell; up to $\sim$$1.00\,\mathrm{M_\odot}$ depending on the total host extinction) progenitor ignited by a thick ($\sim$$0.13\,\mathrm{M_\odot}$) He-shell. We detect a prominent absorption feature at $\sim$1 $\mu\mathrm{m}$ in the near-infrared (NIR) spectrum of SN 2020jgb, which could originate from unburnt helium in the outermost ejecta. While the sample size is limited, similar 1 $\mu\mathrm{m}$ features have been detected in all the thick He-shell DDet candidates with NIR spectra obtained to date. SN 2020jgb is also the first subluminous, thick He-shell DDet SN discovered in a star-forming galaxy, indisputably showing that He-shell DDet objects occur in both star-forming and passive galaxies, consistent with the normal SN Ia population., Comment: 23 pages, 10 figures. Updated to accepted version (ApJ)

We present a method for analyzing supernova remnants (SNRs) by diagnosing the drivers responsible for structure at different angular scales. First, we perform a suite of hydrodynamic models of the Rayleigh-Taylor instability (RTI) as a supernova collides with its surrounding medium. Using these models we demonstrate how power spectral analysis can be used to attribute which scales in a SNR are driven by RTI and which must be caused by intrinsic asymmetries in the initial explosion. We predict the power spectrum of turbulence driven by RTI and identify a dominant angular mode which represents the largest scale that efficiently grows via RTI. We find that this dominant mode relates to the density scale height in the ejecta, and therefore reveals the density profile of the SN ejecta. If there is significant structure in a SNR on angular scales larger than this mode, then it is likely caused by anisotropies in the explosion. Structure on angular scales smaller than the dominant mode exhibits a steep scaling with wavenumber, possibly too steep to be consistent with a turbulent cascade, and therefore might be determined by the saturation of RTI at different length scales (although systematic 3D studies are needed to investigate this). We also demonstrate, consistent with previous studies, that this power spectrum is independent of the magnitude and length scales of perturbations in the surrounding medium and therefore this diagnostic is unaffected by ``clumpiness" in the CSM., Comment: 9 pages, 5 figures. Published in ApJL

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

A thermonuclear explosion triggered by a helium-shell detonation on a carbon-oxygen white dwarf core has been predicted to have strong UV line blanketing at early times due to the iron-group elements produced during helium-shell burning. We present the photometric and spectroscopic observations of SN 2016dsg, a sub-luminous peculiar Type I SN consistent with a thermonuclear explosion involving a thick He shell. With a redshift of 0.04, the $i$-band peak absolute magnitude is derived to be around -17.5. The object is located far away from its host, an early-type galaxy, suggesting it originated from an old stellar population. The spectra collected after the peak are unusually red, show strong UV line blanketing and weak O I $\lambda$7773 absorption lines, and do not evolve significantly over 30 days. An absorption line around 9700-10500 \AA is detected in the near-infrared spectrum and is likely from the unburnt helium in the ejecta. The spectroscopic evolution is consistent with the thermonuclear explosion models for a sub-Chandrasekhar mass white dwarf with a thick helium shell, while the photometric evolution is not well described by existing models., Comment: accepted for publication in ApJ

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

We present nebular spectra of the Type Ia supernova (SN Ia) SN 2019yvq, which had a bright flash of blue and ultraviolet light after exploding, followed by a rise similar to other SNe Ia. Although SN 2019yvq displayed several other rare characteristics such as persistent high ejecta velocity near peak brightness, it was not especially peculiar and if the early "excess" emission were not observed, it would likely be included in cosmological samples. The excess flux can be explained by several different physical models linked to the details of the progenitor system and explosion mechanism. Each has unique predictions for the optically thin emission at late times. In our nebular spectra, we detect strong [Ca II] $\lambda\lambda$7291, 7324 and Ca NIR triplet emission, consistent with a double-detonation explosion. We do not detect H, He, or [O I] emission, predictions for some single-degenerate progenitor systems and violent white dwarf mergers. The amount of swept-up H or He is <2.8 x 10^-4 and 2.4 x 10^-4 Msun, respectively. Aside from strong Ca emission, the SN 2019yvq nebular spectrum is similar to those of typical SNe Ia with the same light-curve shape. Comparing to double-detonation models, we find that the Ca emission is consistent with a model with a total progenitor mass of 1.15 Msun. However, we note that a lower progenitor mass better explains the early light-curve and peak luminosity. The unique properties of SN 2019yvq suggest that thick He-shell double-detonations only account for $1.1^{+2.1}_{-1.1}\%$ of the total "normal" SN Ia rate. SN 2019yvq is one of the best examples yet that multiple progenitor channels appear necessary to reproduce the full diversity of "normal" SNe Ia., Comment: 16 pages, 6 figures, accepted in ApJL on August 12 2020

(Abridged) Using the Zwicky Transient Facility alert stream, we are conducting a large campaign to spectroscopically classify all transients occurring in galaxies in the Census of the Local Universe (CLU) catalog. The aim of the experiment is to construct a spectroscopically complete, volume-limited sample of transients coincident within 100" of CLU galaxies out to 200 Mpc, and to a depth of 20 mag. We describe the survey design and spectroscopic completeness from the first 16 months of operations. We present results from a systematic search for Calcium rich gap transients in the sample of 22 low luminosity (peak absolute magnitude $M > -17$), hydrogen poor events found in the experiment (out of 754 spectroscopically classified SNe). We report the detection of eight Calcium rich gap transients, and constrain their volumetric rate to be at least $\approx 15\pm5$% of the SN Ia rate. Combining this sample with ten events from the literature, we find a likely continuum of spectroscopic properties ranging from events with SN Ia-like features (Ca-Ia objects) to SN Ib/c-like features (Ca-Ib/c objects) at peak light. Within the Ca-Ib/c events, we find two populations of events distinguished by their red ($g - r \approx 1.5$ mag) or green ($g - r \approx 0.5$ mag) spectral colors at $r$-band peak, wherein redder events show strong line blanketing signatures, slower light curves, weaker He lines and lower [Ca II]/[O I] in the nebular phase. Together, we find that the spectroscopic continuum, volumetric rates and striking old environments are consistent with the explosive burning of He shells on low mass white dwarfs. We posit that Ca-Ia and red Ca-Ib/c objects are consistent with the double detonation of He shells with high He burning efficiency, while green Ca-Ib/c objects could arise from less efficient He burning scenarios such as detonations in low density He shells or He shell deflagrations., Comment: 30 figures, 6 tables. Submitted to ApJ

Colors of Type Ia supernovae in the first few days after explosion provide a potential discriminant between different models. In this paper, we present $g-r$ colors of 65 Type Ia supernovae discovered within 5 days from first light by the Zwicky Transient Facility in 2018, a sample that is about three times larger than that in the literature. We find that $g-r$ colors are intrinsically rather homogeneous at early phases, with about half of the dispersion attributable to photometric uncertainties ($\sigma_\mathrm{noise}\sim\sigma_\mathrm{int}\sim$ 0.18 mag). Colors are nearly constant starting from 6 days after first light ($g-r\sim-0.15$ mag), while the time evolution at earlier epochs is characterized by a continuous range of slopes, from events rapidly transitioning from redder to bluer colors (slope of $\sim-0.25$ mag day$^{-1}$) to events with a flatter evolution. The continuum in the slope distribution is in good agreement both with models requiring some amount of $^{56}$Ni mixed in the outermost regions of the ejecta and with "double-detonation" models having thin helium layers ($M_\mathrm{He}=0.01\,M_\odot$) and varying carbon-oxygen core masses. At the same time, six events show evidence for a distinctive "red bump" signature predicted by "double-detonation" models with larger helium masses. We finally identify a significant correlation between the early-time $g-r$ slopes and supernova brightness, with brighter events associated to flatter color evolution (p-value=0.006). The distribution of slopes, however, is consistent with being drawn from a single population, with no evidence for two components as claimed in the literature based on $B-V$ colors., Comment: 17 pages, 9 figures, 1 table; accepted for publication in ApJ; fixed error in time of first light + selection criterion changed; sample increased from 38 to 65 events; conclusions unchanged

We use non-local thermal equilibrium radiative transport modeling to examine observational signatures of sub-Chandrasekhar mass double detonation explosions in the nebular phase. Results range from spectra that look like typical and subluminous Type Ia supernovae (SNe) for higher mass progenitors to spectra that look like Ca-rich transients for lower mass progenitors. This ignition mechanism produces an inherent relationship between emission features and the progenitor mass as the ratio of the nebular [Ca II]/[Fe III] emission lines increases with decreasing white dwarf mass. Examining the [Ca II]/[Fe III] nebular line ratio in a sample of observed SNe we find further evidence for the two distinct classes of SNe Ia identified in Polin et al. by their relationship between Si II velocity and B-band magnitude, both at time of peak brightness. This suggests that SNe Ia arise from more than one progenitor channel, and provides an empirical method for classifying events based on their physical origin. Furthermore, we provide insight to the mysterious origin of Ca-rich transients. Low-mass double detonation models with only a small mass fraction of Ca (1%) produce nebular spectra that cool primarily through forbidden [Ca II] emission., Comment: 14 pages, 9 figures, published by ApJ

We present observations and modeling of SN 2016hnk, a Ca-rich supernova (SN) that is consistent with being the result of a He-shell double-detonation explosion of a C/O white dwarf. We find that SN 2016hnk is intrinsically red relative to typical thermonuclear SNe and has a relatively low peak luminosity ($M_B = -15.4$ mag), setting it apart from low-luminosity Type Ia supernovae (SNe Ia). SN 2016hnk has a fast-rising light curve that is consistent with other Ca-rich transients ($t_r = 15$ d). We determine that SN 2016hnk produced $0.03 \pm 0.01 M_{\odot}$ of ${}^{56}\textrm{Ni}$ and $0.9 \pm 0.3 M_{\odot}$ of ejecta. The photospheric spectra show strong, high-velocity Ca II absorption and significant line blanketing at $\lambda < 5000$ Angstroms, making it distinct from typical (SN 2005E-like) Ca-rich SNe. SN 2016hnk is remarkably similar to SN 2018byg, which was modeled as a He-shell double-detonation explosion. We demonstrate that the spectra and light curves of SN 2016hnk are well modeled by the detonation of a $0.02 \ M_{\odot}$ helium shell on the surface of a $0.85 \ M_{\odot}$ C/O white dwarf. This analysis highlights the second observed case of a He-shell double-detonation and suggests a specific thermonuclear explosion that is physically distinct from SNe that are defined simply by their low luminosities and strong [Ca II] emission., Comment: 33 pages, 18 Figures. Accepted for publication in ApJ

The detonation of a helium shell on a white dwarf has been proposed as a possible explosion triggering mechanism for Type Ia supernovae. Here, we report ZTF 18aaqeasu (SN 2018byg/ATLAS 18pqq), a peculiar Type I supernova, consistent with being a helium-shell double-detonation. With a rise time of $\approx 18$ days from explosion, the transient reached a peak absolute magnitude of $M_R \approx -18.2$ mag, exhibiting a light curve akin to sub-luminous SN 1991bg-like Type Ia supernovae, albeit with an unusually steep increase in brightness within a week from explosion. Spectra taken near peak light exhibit prominent Si absorption features together with an unusually red color ($g-r \approx 2$ mag) arising from nearly complete line blanketing of flux blue-wards of 5000 \AA. This behavior is unlike any previously observed thermonuclear transient. Nebular phase spectra taken at and after $\approx 30$ days from peak light reveal evidence of a thermonuclear detonation event dominated by Fe-group nucleosynthesis. We show that the peculiar properties of ZTF 18aaqeasu are consistent with the detonation of a massive ($\approx 0.15$ M$_\odot$) helium shell on a sub-Chandrasekhar mass ($\approx 0.75$ M$_\odot$) white dwarf after including mixing of $\approx 0.2$ M$_\odot$ of material in the outer ejecta. These observations provide evidence of a likely rare class of thermonuclear supernovae arising from detonations of massive helium shells., Comment: 10 pages, 6 figures. Submitted to ApJL

We present a numerical parameter survey of sub-Chandrasekhar mass white dwarf (WD) explosions. Carbon-oxygen WDs accreting a helium shell have the potential to explode in the sub-Chandrasekhar mass regime. Previous studies have shown how the ignition of a helium shell can either directly ignite the WD at the core-shell interface or propagate a shock wave into the the core causing a central ignition. We examine the explosions of WDs from 0.6 - 1.2 M$_\odot$ with helium shells of 0.01, 0.05 and 0.08 M$_\odot$. Distinct observational signatures of sub-Chandrasekhar mass WD explosions are predicted for two categories of shell size. Thicker-shell models show an early time flux excess, which is caused by the presence of radioactive material in the ashes of the helium shell, and red colors due to these ashes creating significant line blanketing in the UV through the blue portion of the spectrum. Thin shell models reproduce several typical Type Ia supernova signatures. We identify a relationship between Si II velocity and luminosity which, for the first time, identifies a sub-class of observed supernovae that are consistent with these models. This sub-class is further delineated by the absence of carbon in their atmospheres. We suggest that the proposed difference in the ratio of selective to total extinction between the high velocity and normal velocity Type Ia supernovae is not due to differences in the properties of the dust around these events, but is rather an artifact of applying a single extinction correction to two intrinsically different populations of supernovae., Comment: 13 pages, 12 figures, Accepted to ApJ 30 Decemeber 2018

The steadily improving sensitivity of pulsar timing arrays (PTAs) suggests that gravitational waves (GWs) from supermassive black hole binary (SMBHB) systems in the nearby universe will be de- tectable sometime during the next decade. Currently, PTAs assume an equal probability of detection from every sky position, but as evidence grows for a non-isotropic distribution of sources, is there a most likely sky position for a detectable single source of GWs? In this paper, a collection of galactic catalogs is used to calculate various metrics related to the detectability of a single GW source resolv- able above a GW background, assuming that every galaxy has the same probability of containing a SMBHB. Our analyses of these data reveal small probabilities that one of these sources is currently in the PTA band, but as sensitivity is improved regions of consistent probability density are found in predictable locations, specifically around local galaxy clusters., Comment: 9 pages, 9 figures, accepted for submission in ApJ