Showing 1–50 of 74 results for author: Shen, K J

Local two-dimensional simulations of the ignition of a helium shell detonation on a white dwarf by an impacting stream

Authors: Nethra Rajavel, Dean M. Townsley, Ken J. Shen

Abstract: The double detonation model is one of the prevalent explosion mechanisms of Type Ia Supernovae (SNe Ia) wherein an outer helium shell detonation triggers a core detonation in the white dwarf (WD). The dynamically driven double degenerate double detonation (D6) is the double detonation of the more massive WD in a binary WD system where the localized impact of the mass transfer stream from the compa… ▽ More The double detonation model is one of the prevalent explosion mechanisms of Type Ia Supernovae (SNe Ia) wherein an outer helium shell detonation triggers a core detonation in the white dwarf (WD). The dynamically driven double degenerate double detonation (D6) is the double detonation of the more massive WD in a binary WD system where the localized impact of the mass transfer stream from the companion sets off the initial helium shell detonation. To have high numerical resolution and control over the stream parameters, we have implemented a study of the local interaction of the stream with the WD surface in 2D. In cases with lower base density of the shell, the stream's impact can cause surface detonation soon after first impact. With higher base densities, after the stream hits the surface, hot material flows around the star and interacts with the incoming stream to produce a denser and narrower impact. Our results therefore show that (1) a directly impacting stream for both a relatively high resolution and for a range of stream parameters can produce a surface detonation, (2) thinner helium shells ignite more promptly via impact, doing so sooner, and (3) there are lower limits on ignition in both shell density and incoming stream speed with lower limits on density being well below those shown by other work to be required for normal appearing SN Ia. This supports stream ignition and therefore the D6 scenario, as a viable mechanism for normal SNe Ia. △ Less

Submitted 20 August, 2024; originally announced August 2024.

Comments: 15 pages, 8 figures, 3 tables, submitted to Apj

Ultraviolet spectroscopy of the supernova Ia hypervelocity runaway white dwarf J0927-6335

Authors: Klaus Werner, Kareem El-Badry, Boris T. Gänsicke, Ken J. Shen

Abstract: The hot white dwarf (WD) J0927-6335 (Gaia DR3 5250394728194220800, effective temperature T$_{\rm eff}$ = 60,000 K, surface gravity log g = 7) was detected as the fastest known Galactic hypervelocity star with a space velocity of $\approx$2800 km s$^{-1}$ and an atmosphere dominated by carbon and oxygen. It is thought to be the surviving WD donor predicted by the "dynamically driven double-degenera… ▽ More The hot white dwarf (WD) J0927-6335 (Gaia DR3 5250394728194220800, effective temperature T$_{\rm eff}$ = 60,000 K, surface gravity log g = 7) was detected as the fastest known Galactic hypervelocity star with a space velocity of $\approx$2800 km s$^{-1}$ and an atmosphere dominated by carbon and oxygen. It is thought to be the surviving WD donor predicted by the "dynamically driven double-degenerate double-detonation" (D$^6$) type Ia supernova formation model. We analysed an ultraviolet spectrum of J0927-6335 obtained recently with the Hubble Space Telescope and found very high abundances of iron and nickel. This could originate in the pollution of the remnant by the SN Ia explosion but it is uncertain to what extent atomic diffusion altered the chemical composition of the accreted material. △ Less

Submitted 15 August, 2024; originally announced August 2024.

Comments: A&A Letters, accepted

Supernova Shocks Cannot Explain the Inflated State of Hypervelocity Runaways from White Dwarf Binaries

Authors: Aakash Bhat, Evan B. Bauer, Rüdiger Pakmor, Ken J. Shen, Ilaria Caiazzo, Abinaya Swaruba Rajamuthukumar, Kareem El-Badry, Wolfgang E. Kerzendorf

Abstract: Recent observations have found a growing number of hypervelocity stars with speeds of $\approx 1500-2500\,$km\,s$^{-1}$ which could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been… ▽ More Recent observations have found a growing number of hypervelocity stars with speeds of $\approx 1500-2500\,$km\,s$^{-1}$ which could have only been produced through thermonuclear supernovae in white dwarf binaries. Most of the observed hypervelocity runaways in this class display a surprising inflated structure: their current radii are roughly an order of magnitude greater than they would have been as white dwarfs filling their Roche lobe. While many simulations exist studying the dynamical phase leading to supernova detonation in these systems, no detailed calculations of the long-term structure of the runaways have yet been performed. We use an existing \textsc{Arepo} hydrodynamical simulation of a supernova in a white dwarf binary as a starting point for the evolution of these stars with the 1 dimensional stellar evolution code MESA. We show that the supernova shock is not enough to inflate the white dwarf over timescales longer than a few thousand years, significantly shorter than the $10^{5-6}$ year lifetimes inferred for observed hypervelocity runaways. Despite experiencing a shock from a supernova less than $\approx 0.02\,R_\odot$ away, our models do not experience significant interior heating, and all contract back to radii around $0.01\,R_\odot$ within about $10^4$\,years. Explaining the observed inflated states requires either an additional source of significant heating or some other physics that is not yet accounted for in the subsequent evolution. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: Submitted to A\&A. 15 pages, 17 figures

Almost All Carbon/Oxygen White Dwarfs Can Support Double Detonations

Authors: Ken J. Shen, Samuel J. Boos, Dean M. Townsley

Abstract: Double detonations of sub-Chandrasekhar-mass white dwarfs (WDs) in unstably mass-transferring double WD binaries have become a leading contender to explain most, if not all, Type Ia supernovae. However, past theoretical studies of the explosion process have assumed relatively ad hoc initial conditions for the helium shells in which the double detonations begin. In this work, we construct realistic… ▽ More Double detonations of sub-Chandrasekhar-mass white dwarfs (WDs) in unstably mass-transferring double WD binaries have become a leading contender to explain most, if not all, Type Ia supernovae. However, past theoretical studies of the explosion process have assumed relatively ad hoc initial conditions for the helium shells in which the double detonations begin. In this work, we construct realistic C/O WDs to use as the starting points for multidimensional double detonation simulations. We supplement these with simplified one-dimensional detonation calculations to gain a physical understanding of the conditions under which shell detonations can propagate successfully. We find that C/O WDs <= 1.0 Msol, which make up the majority of C/O WDs, are born with structures that can support double detonations. More massive C/O WDs require ~1e-3 Msol of accretion before detonations can successfully propagate in their shells, but such accretion may be common in the double WD binaries that host massive WDs. Our findings strongly suggest that if the direct impact accretion stream reaches high enough temperatures and densities during mass transfer from one WD to another, the accreting WD will undergo a double detonation. Furthermore, if the companion is also a C/O WD <= 1.0 Msol, it will undergo its own double detonation when impacted by the ejecta from the first explosion. Exceptions to this outcome may explain the newly discovered class of hypervelocity supernova survivors. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: Submitted

Type Ia Supernovae Can Arise from the Detonations of Both Stars in a Double Degenerate Binary

Authors: Samuel J. Boos, Dean M. Townsley, Ken J. Shen

Abstract: The precise origin of Type Ia supernovae (SNe Ia) is unknown despite their value to numerous areas in astronomy. While it is a long-standing consensus that they arise from an explosion of a carbon/oxygen white dwarf, the exact progenitor configurations and explosion mechanisms that lead to SNe Ia are still debated. One popular theory is the double detonation in which a helium layer, accreted from… ▽ More The precise origin of Type Ia supernovae (SNe Ia) is unknown despite their value to numerous areas in astronomy. While it is a long-standing consensus that they arise from an explosion of a carbon/oxygen white dwarf, the exact progenitor configurations and explosion mechanisms that lead to SNe Ia are still debated. One popular theory is the double detonation in which a helium layer, accreted from a binary companion, detonates on the surface of the primary star, leading to a converging shock-induced detonation of the underlying core. It has recently been seen in simulations that a helium-rich degenerate companion may undergo its own explosion triggered by the impact from the ejecta of the primary star. We show 2D simulations that approximate a white dwarf undergoing a double detonation which triggers the explosion of the degenerate companion, leading to either a triple or quadruple detonation. We also present the first multi-dimensional radiative transfer results from the triple and quadruple detonation scenario. We find that within a range of mass configurations of the degenerate binary, the synthetic light curves and spectra of these events match observations as well as theoretical models of isolated double detonations do. Notably, double and quadruple detonations that are spectrally similar and reach the same peak brightnesses have drastically different ejection masses and produce different amounts of Si- and Fe-group elements. Further understanding of this scenario is needed in order to determine if at least some observed SNe Ia actually originate from two stars exploding. △ Less

Submitted 8 July, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: 22 pages, 21 figures, 3 tables. Accepted for publication in ApJ

Evidence for saturated and disrupted magnetic braking from samples of detached close binaries with M and K dwarfs

Authors: Diogo Belloni, Matthias R. Schreiber, Maxwell Moe, Kareem El-Badry, Ken J. Shen

Abstract: Context. Recent observations of close detached eclipsing M and K dwarf binaries have provided substantial support for magnetic saturation when stars rotate sufficiently fast, leading to a magnetic braking (MB) torque proportional to the spin of the star. Aims. We investigated here how strong MB torques need to be to reproduce the observationally-inferred relative numbers of white dwarf plus M dw… ▽ More Context. Recent observations of close detached eclipsing M and K dwarf binaries have provided substantial support for magnetic saturation when stars rotate sufficiently fast, leading to a magnetic braking (MB) torque proportional to the spin of the star. Aims. We investigated here how strong MB torques need to be to reproduce the observationally-inferred relative numbers of white dwarf plus M dwarf post-common-envelope binaries under the assumption of magnetic saturation. Methods. We carried out binary population simulations with the BSE code adopting empirically-derived inter-correlated main-sequence binary distributions as initial binary populations and compared the simulation outcomes with observations. Results. We found that the dearth of extreme mass ratio binaries in the inter-correlated initial distributions is key to reproduce the large fraction of post-common-envelope binaries hosting low-mass M dwarfs (${\sim0.1-0.2}$ M$_\odot$). In addition, orbital angular momentum loss rates due to MB should be high for M dwarfs with radiative cores and orders of magnitude smaller for fully convective stars to explain the observed dramatic change of the fraction of short-period binaries at the fully convective boundary. Conclusions. We conclude that saturated but disrupted, that is, dropping drastically at the fully convective boundary, MB can explain the observations of both close main-sequence binaries containing M and K dwarfs and post-common-envelope binaries. Whether a similar prescription can explain the spin down rates of single stars and of binaries containing more massive stars needs to be tested. △ Less

Submitted 7 November, 2023; originally announced November 2023.

Comments: Accepted for publication in A&A

Revisiting the classics: On the evolutionary origin of the "Fe II" and "He/N" spectral classes of novae

Authors: E. Aydi, L. Chomiuk, J. Strader, K. V. Sokolovsky, R. E. Williams, D. A. H. Buckley, A. Ederoclite, L. Izzo, R. Kyer, J. D. Linford, A. Kniazev, B. D. Metzger, J. Mikolajewska, P. Molaro, I. Mollina, K. Mukai, U. Munari, M. Orio, T. Panurach, B. J. Shappee, K. J. Shen, J. L. Sokoloski, R. Urquhart, F. M. Walter

Abstract: The optical spectra of novae are characterized by emission lines from the hydrogen Balmer series and either Fe II or He/N, leading to their traditional classification into two spectral classes: "Fe II" and "He/N". For decades, the origins of these spectral features were discussed in the literature in the contexts of different bodies of gas or changes in the opacity of the ejecta, particularly asso… ▽ More The optical spectra of novae are characterized by emission lines from the hydrogen Balmer series and either Fe II or He/N, leading to their traditional classification into two spectral classes: "Fe II" and "He/N". For decades, the origins of these spectral features were discussed in the literature in the contexts of different bodies of gas or changes in the opacity of the ejecta, particularly associated with studies by R. E. Williams and S. N. Shore. Here, we revisit these major studies with dedicated, modern data sets, covering the evolution of several novae from early rise to peak all the way to the nebular phase. Our data confirm previous suggestions in the literature that the "Fe II" and "He/N" spectral classes are phases in the spectroscopic evolution of novae driven primarily by changes in the opacity, ionization, and density of the ejecta, and most if not all novae go through at least three spectroscopic phases as their eruptions evolve: an early He/N (phase 1; observed during the early rise to visible peak and characterized by P Cygni lines of He I, N II, and N III), then an Fe II (phase 2; observed near visible peak and characterized by P Cygni lines of Fe II and O I), and then a later He/N (phase 3; observed during the decline and characterized by emission lines of He I. He II, N II, and N III), before entering the nebular phase. This spectral evolution seems to be ubiquitous across novae, regardless of their speed class; however the duration of each of these phase differs based on the speed class of the nova. △ Less

Submitted 27 October, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 21 pages, 14 figures, 11 tables, Accepted in MNRAS

SN 2022joj: A Potential Double Detonation with a Thin Helium shell

Authors: E. Padilla Gonzalez, D. A. Howell, G. Terreran, C. McCully, M. Newsome, J. Burke, J. Farah, C. Pellegrino, K. A. Bostroem, G. Hosseinzadeh, J. Pearson, D. J. Sand, M. Shrestha, N. Smith, Y. Dong, N. Meza Retamal, S. Valenti, S. Boos, K. J. Shen, D. Townsley, L. Galbany, L. Piscarreta, R. J. Foley, M. J. Bustamante-Rosell, D. A. Coulter , et al. (12 additional authors not shown)

Abstract: We present photometric and spectroscopic data for SN 2022joj, a nearby peculiar Type Ia supernova (SN Ia) with a fast decline rate ($\rm{Δm_{15,B}=1.4}$ mag). SN 2022joj shows exceedingly red colors, with a value of approximately ${B-V \approx 1.1}$ mag during its initial stages, beginning from $11$ days before maximum brightness. As it evolves the flux shifts towards the blue end of the spectrum,… ▽ More We present photometric and spectroscopic data for SN 2022joj, a nearby peculiar Type Ia supernova (SN Ia) with a fast decline rate ($\rm{Δm_{15,B}=1.4}$ mag). SN 2022joj shows exceedingly red colors, with a value of approximately ${B-V \approx 1.1}$ mag during its initial stages, beginning from $11$ days before maximum brightness. As it evolves the flux shifts towards the blue end of the spectrum, approaching ${B-V \approx 0}$ mag around maximum light. Furthermore, at maximum light and beyond, the photometry is consistent with that of typical SNe Ia. This unusual behavior extends to its spectral characteristics, which initially displayed a red spectrum and later evolved to exhibit greater consistency with typical SNe Ia. We consider two potential explanations for this behavior: double detonation from a helium shell on a sub-Chandrasekhar-mass white dwarf and Chandrasekhar-mass models with a shallow distribution of $\rm{^{56}Ni}$. The shallow nickel models could not reproduce the red colors in the early light curves. Spectroscopically, we find strong agreement between SN 2022joj and double-detonation models with white dwarf masses around 1 $\rm{M_{\odot}}$ and thin He-shell between 0.01 and 0.02 $\rm{M_{\odot}}$. Moreover, the early red colors are explained by line-blanketing absorption from iron-peak elements created by the double detonation scenario in similar mass ranges. However, the nebular spectra composition in SN 2022joj deviates from expectations for double detonation, as we observe strong [Fe III] emission instead of [Ca II] lines as anticipated from double detonation models. More detailed modeling, e.g., including viewing angle effects, is required to test if double detonation models can explain the nebular spectra. △ Less

Submitted 11 August, 2023; originally announced August 2023.

SN 2022joj: A Peculiar Type Ia Supernova Possibly Driven by an Asymmetric Helium-shell Double Detonation

Authors: Chang Liu, Adam A. Miller, Samuel J. Boos, Ken J. Shen, Dean M. Townsley, Steve Schulze, Luke Harvey, Kate Maguire, Joel Johansson, Thomas G. Brink, Umut Burgaz, Georgios Dimitriadis, Alexei V. Filippenko, Saarah Hall, K-Ryan Hinds, Andrew Hoffman, Viraj Karambelkar, Charles D. Kilpatrick, Daniel Perley, Neil Pichay, Huei Sears, Jesper Sollerman, Robert Stein, Jacco H. Terwel, WeiKang Zheng , et al. (6 additional authors not shown)

Abstract: We present observations of SN 2022joj, a peculiar Type Ia supernova (SN Ia) discovered by the Zwicky Transient Facility (ZTF). SN 2022joj exhibits an unusually red $g_\mathrm{ZTF}-r_\mathrm{ZTF}$ color at early times and a rapid blueward evolution afterward. Around maximum brightness, SN 2022joj shows a high luminosity ($M_{g_\mathrm{ZTF},\mathrm{max}}\simeq-19.7$ mag), a blue broadband color (… ▽ More We present observations of SN 2022joj, a peculiar Type Ia supernova (SN Ia) discovered by the Zwicky Transient Facility (ZTF). SN 2022joj exhibits an unusually red $g_\mathrm{ZTF}-r_\mathrm{ZTF}$ color at early times and a rapid blueward evolution afterward. Around maximum brightness, SN 2022joj shows a high luminosity ($M_{g_\mathrm{ZTF},\mathrm{max}}\simeq-19.7$ mag), a blue broadband color ($g_\mathrm{ZTF}-r_\mathrm{ZTF}\simeq-0.2$ mag), and shallow Si II absorption lines, consistent with those of overluminous, SN 1991T-like events. The maximum-light spectrum also shows prominent absorption around 4200 Å, which resembles the Ti II features in subluminous, SN 1991bg-like events. Despite the blue optical-band colors, SN 2022joj exhibits extremely red ultraviolet minus optical colors at maximum luminosity ($u-v\simeq0.6$ mag and $uvw1 - v\simeq2.5$ mag), suggesting a suppression of flux at $\sim$2500--4000 Å. Strong C II lines are also detected at peak. We show that these unusual spectroscopic properties are broadly consistent with the helium-shell double detonation of a sub-Chandrasekhar mass ($M\simeq1 \mathrm{M_\odot}$) carbon/oxygen (C/O) white dwarf (WD) from a relatively massive helium shell ($M_s\simeq0.04$--$0.1 \mathrm{M_\odot}$), if observed along a line of sight roughly opposite to where the shell initially detonates. None of the existing models could quantitatively explain all the peculiarities observed in SN 2022joj. The low flux ratio of [Ni II] $λ$7378 to [Fe II] $λ$7155 emission in the late-time nebular spectra indicates a low yield of stable Ni isotopes, favoring a sub-Chandrasekhar mass progenitor. The significant blueshift measured in the [Fe II] $λ$7155 line is also consistent with an asymmetric chemical distribution in the ejecta, as is predicted in double-detonation models. △ Less

Submitted 24 November, 2023; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: 24 pages, 11 figures, 6 tables. Updated to accepted version (ApJ)

Journal ref: ApJ 958 178 (2023)

The Q Branch Cooling Anomaly Can Be Explained by Mergers of White Dwarfs and Subgiant Stars

Authors: Ken J. Shen, Simon Blouin, Katelyn Breivik

Abstract: Gaia's exquisite parallax measurements allowed for the discovery and characterization of the Q branch in the Hertzsprung-Russell diagram, where massive C/O white dwarfs (WDs) pause their dimming due to energy released during crystallization. Interestingly, the fraction of old stars on the Q branch is significantly higher than in the population of WDs that will become Q branch stars or that were Q… ▽ More Gaia's exquisite parallax measurements allowed for the discovery and characterization of the Q branch in the Hertzsprung-Russell diagram, where massive C/O white dwarfs (WDs) pause their dimming due to energy released during crystallization. Interestingly, the fraction of old stars on the Q branch is significantly higher than in the population of WDs that will become Q branch stars or that were Q branch stars in the past. From this, Cheng et al. inferred that ~6% of WDs passing through the Q branch experience a much longer cooling delay than that of standard crystallizing WDs. Previous attempts to explain this cooling anomaly have invoked mechanisms involving super-solar initial metallicities. In this paper, we describe a novel scenario in which a standard composition WD merges with a subgiant star. The evolution of the resulting merger remnant leads to the creation of a large amount of 26Mg, which, along with the existing 22Ne, undergoes a distillation process that can release enough energy to explain the Q branch cooling problem without the need for atypical initial abundances. The anomalously high number of old stars on the Q branch may thus be evidence that mass transfer from subgiants to WDs leads to unstable mergers. △ Less

Submitted 13 September, 2023; v1 submitted 8 August, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJL. Added text and a figure to better motivate the initial conditions of the merger remnant evolution. Also amended text regarding the estimated numbers of WD + subgiant mergers

Type Ia Supernova Nucleosynthesis: Metallicity-Dependent Yields

Authors: James D Keegans, Marco Pignatari, Richard J Stancliffe, Claudia Travaglio, Samuel Jones, Brad K Gibson, Dean M Townsley, Broxton J Miles, Ken J Shen, Gareth Few

Abstract: Type Ia supernova explosions (SNIa) are fundamental sources of elements for the chemical evolution of galaxies. They efficiently produce intermediate-mass (with Z between 11 and 20) and iron group elements - for example, about 70% of the solar iron is expected to be made by SNIa. In this work, we calculate complete abundance yields for 39 models of SNIa explosions, based on three progenitors - a 1… ▽ More Type Ia supernova explosions (SNIa) are fundamental sources of elements for the chemical evolution of galaxies. They efficiently produce intermediate-mass (with Z between 11 and 20) and iron group elements - for example, about 70% of the solar iron is expected to be made by SNIa. In this work, we calculate complete abundance yields for 39 models of SNIa explosions, based on three progenitors - a 1.4M deflagration detonation model, a 1.0 double detonation model and a 0.8 M double detonation model - and 13 metallicities, with 22Ne mass fractions of 0, 1x10-7, 1x10-6, 1x10-5, 1x10-4, 1x10-3, 2x10-3, 5x10-3, 1x10-2, 1.4x10-2, 5x10-2, and 0.1 respectively. Nucleosynthesis calculations are done using the NuGrid suite of codes, using a consistent nuclear reaction network between the models. Complete tables with yields and production factors are provided online at Zenodo: Yields. We discuss the main properties of our yields in the light of the present understanding of SNIa nucleosynthesis, depending on different progenitor mass and composition. Finally, we compare our results with a number of relevant models from the literature. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 42 pages, 21 figures. Accepted for publication in ApJS 21-06-23

No plateau observed in late-time near-infrared observations of the underluminous Type Ia supernova 2021qvv

Authors: O. Graur, E. Padilla Gonzalez, J. Burke, M. Deckers, S. W. Jha, L. Galbany, E. Karamenhmetoglu, M. D. Stritzinger, K. Maguire, D. A. Howell, R. Fisher, A. G. Fullard, R. Handberg, D. Hiramatsu, G. Hosseinzadeh, W. E. Kerzendorf, C. McCully, M. Newsome, C. Pellegrino, A. Rest, A. G. Riess, I. R. Seitenzahl, M. M. Shara, K. J. Shen, G. Terreran , et al. (1 additional authors not shown)

Abstract: Near-infrared (NIR) observations of normal Type Ia supernovae (SNe Ia) obtained between 150 to 500 d past maximum light reveal the existence of an extended plateau. Here, we present observations of the underluminous, 1991bg-like SN 2021qvv. Early, ground-based optical and NIR observations show that SN 2021qvv is similar to SN 2006mr, making it one of the dimmest, fastest-evolving 1991bg-like SNe t… ▽ More Near-infrared (NIR) observations of normal Type Ia supernovae (SNe Ia) obtained between 150 to 500 d past maximum light reveal the existence of an extended plateau. Here, we present observations of the underluminous, 1991bg-like SN 2021qvv. Early, ground-based optical and NIR observations show that SN 2021qvv is similar to SN 2006mr, making it one of the dimmest, fastest-evolving 1991bg-like SNe to date. Late-time (170-250 d) Hubble Space Telescope observations of SN 2021qvv reveal no sign of a plateau. An extrapolation of these observations backwards to earlier-phase NIR observations of SN 2006mr suggests the complete absence of a NIR plateau, at least out to 250 d. This absence may be due to a higher ionization state of the ejecta, as predicted by certain sub-Chandrasekhar-mass detonation models, or to the lower temperatures of the ejecta of 1991bg-like SNe, relative to normal SNe Ia, which might preclude their becoming fluorescent and shifting ultraviolet light into the NIR. This suggestion can be tested by acquiring NIR imaging of a sample of 1991bg-like SNe that covers the entire range from slowly-evolving to fast-evolving events ($0.2 \lesssim s_\mathrm{BV} \lesssim 0.6$). A detection of the NIR plateau in slower-evolving, hotter 1991bg-like SNe would provide further evidence that these SNe exist along a continuum with normal SNe Ia. Theoretical progenitor and explosion scenarios would then have to match the observed properties of both SN Ia subtypes. △ Less

Submitted 26 September, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS (14 pages, 9 figures)

The fastest stars in the Galaxy

Authors: Kareem El-Badry, Ken J. Shen, Vedant Chandra, Evan Bauer, Jim Fuller, Jay Strader, Laura Chomiuk, Rohan Naidu, Ilaria Caiazzo, Antonio C. Rodriguez, Pranav Nagarajan, Natsuko Yamaguchi, Zachary P. Vanderbosch, Benjamin R. Roulston, Jan van Roestel, Boris Gänsicke, Jiwon Jesse Han, Kevin B. Burdge, Alexei V. Filippenko, Thomas G. Brink, WeiKang Zheng

Abstract: We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities… ▽ More We report a spectroscopic search for hypervelocity white dwarfs (WDs) that are runaways from Type Ia supernovae (SNe Ia) and related thermonuclear explosions. Candidates are selected from Gaia data with high tangential velocities and blue colors. We find six new runaways, including four stars with radial velocities (RVs) $>1000\,\rm km\,s^{-1}$ and total space velocities $\gtrsim 1300\,\rm km\,s^{-1}$. These are most likely the surviving donors from double-degenerate binaries in which the other WD exploded. The other two objects have lower minimum velocities, $\gtrsim 600\,\rm km\,s^{-1}$, and may have formed through a different mechanism, such as pure deflagration of a WD in a Type Iax supernova. The four fastest stars are hotter and smaller than the previously known "D$^6$ stars," with effective temperatures ranging from $\sim$20,000 to $\sim$130,000 K and radii of $\sim 0.02-0.10\,R_{\odot}$. Three of these have carbon-dominated atmospheres, and one has a helium-dominated atmosphere. Two stars have RVs of $-1694$ and $-2285\rm \,km\,s^{-1}$ -- the fastest systemic stellar RVs ever measured. Their inferred birth velocities, $\sim 2200-2500\,\rm km\,s^{-1}$, imply that both WDs in the progenitor binary had masses $>1.0\,M_{\odot}$. The high observed velocities suggest that a dominant fraction of the observed hypervelocity WD population comes from double-degenerate binaries whose total mass significantly exceeds the Chandrasekhar limit. However, the two nearest and faintest D$^6$ stars have the lowest velocities and masses, suggesting that observational selection effects favor rarer, higher-mass stars. A significant population of fainter low-mass runaways may still await discovery. We infer a birth rate of D$^6$ stars that is consistent with the SN Ia rate. The birth rate is poorly constrained, however, because the luminosities and lifetimes of $\rm D^6$ stars are uncertain. △ Less

Submitted 25 July, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: 26 pages, 17 figures. Accepted to OJA

Catching a nova X-ray/UV flash in the visible? Early spectroscopy of the extremely slow Nova Velorum 2022 (Gaia22alz)

Authors: E. Aydi, L. Chomiuk, J. Mikołajewska, J. Brink, B. D. Metzger, J. Strader, D. A. H. Buckley, E. J. Harvey, T. W. -S. Holoien, L. Izzo, A. Kawash, J. D. Linford, P. Molaro, B. Mollina, P. Mróz, K. Mukai, M. Orio, T. Panurach, P. Senchyna, B. J. Shappee, K. J. Shen, J. L. Sokoloski, K. V. Sokolovsky, R. Urquhart, R. E. Williams

Abstract: We present early spectral observations of the very slow Galactic nova Gaia22alz, over its gradual rise to peak brightness that lasted 180 days. During the first 50 days, when the nova was only 3--4 magnitudes above its normal brightness, the spectra showed narrow (FWHM $\approx$ 400 km s$^{-1}$) emission lines of H Balmer, He I, He II, and C IV, but no P Cygni absorption. A few weeks later, the hi… ▽ More We present early spectral observations of the very slow Galactic nova Gaia22alz, over its gradual rise to peak brightness that lasted 180 days. During the first 50 days, when the nova was only 3--4 magnitudes above its normal brightness, the spectra showed narrow (FWHM $\approx$ 400 km s$^{-1}$) emission lines of H Balmer, He I, He II, and C IV, but no P Cygni absorption. A few weeks later, the high-excitation He II and C IV lines disappeared, and P Cygni profiles of Balmer, He I, and eventually Fe II lines emerged, yielding a spectrum typical of classical novae before peak. We propose that the early spectra of Gaia22alz are produced in the white dwarf's envelope or accretion disk, reprocessing X-ray and ultraviolet emission from the white dwarf after a dramatic increase in the rate of thermonuclear reactions, during a phase known as the ``early X-ray/UV flash''. If true, this would be one of the rare times that the optical signature of the early X-ray/UV flash has been detected. While this phase might last only a few hours in other novae and thus be easily missed, it was possible to detect in Gaia22alz due to its very slow and gradual rise and thanks to the efficiency of new all-sky surveys in detecting transients on their rise. We also consider alternative scenarios that could explain the early spectral features of Gaia22alz and its unusually slow rise. △ Less

Submitted 9 April, 2023; originally announced April 2023.

Comments: 20 pages, 12 figures, 2 tables. Submitted to MNRAS

A Single Binary May Host Recurrent Thermonuclear Supernovae

Authors: Kaela J. Lee, Ken J. Shen

Abstract: The most commonly accepted progenitor system for Type Iax supernovae (SNe Iax) is the partial deflagration of a near-Chandrasekhar-mass white dwarf (WD) accreting from a non-degenerate helium donor star, leaving a bound remnant following the explosion. In this paper, we investigate whether the WD remant can undergo multiple SNe during the system's lifetime. We use Modules for Experiments in Astrop… ▽ More The most commonly accepted progenitor system for Type Iax supernovae (SNe Iax) is the partial deflagration of a near-Chandrasekhar-mass white dwarf (WD) accreting from a non-degenerate helium donor star, leaving a bound remnant following the explosion. In this paper, we investigate whether the WD remant can undergo multiple SNe during the system's lifetime. We use Modules for Experiments in Astrophysics (MESA) to evolve various single-degenerate binaries to determine which could plausibly undergo multiple SNe Iax due to multiple helium accretion phases. We also investigate the possibility for a subsequent Type Ia SN after the formation of a double WD system. Our work concludes that close binaries with relatively high-mass donors produce the highest probability for several thermonuclear SNe. △ Less

Submitted 2 February, 2023; v1 submitted 23 December, 2022; originally announced December 2022.

Comments: 11 pages, 8 figures. Minor updates added in response to community input

StaNdaRT: A repository of standardized test models and outputs for supernova radiative transfer

Authors: Stéphane Blondin, Sergei Blinnikov, Fionntan P. Callan, Christine E. Collins, Luc Dessart, Wesley Even, Andreas Flörs, Andrew G. Fullard, D. John Hillier, Anders Jerkstrand, Daniel Kasen, Boaz Katz, Wolfgang Kerzendorf, Alexandra Kozyreva, Jack O'Brien, Ezequiel A. Pássaro, Nathaniel Roth, Ken J. Shen, Luke Shingles, Stuart A. Sim, Jaladh Singhal, Isaac G. Smith, Elena Sorokina, Victor P. Utrobin, Christian Vogl , et al. (4 additional authors not shown)

Abstract: We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M… ▽ More We present the first results of a comprehensive supernova (SN) radiative-transfer (RT) code-comparison initiative (StaNdaRT), where the emission from the same set of standardized test models is simulated by currently-used RT codes. A total of ten codes have been run on a set of four benchmark ejecta models of Type Ia supernovae. We consider two sub-Chandrasekhar-mass ($M_\mathrm{tot} = 1.0$ M$_\odot$) toy models with analytic density and composition profiles and two Chandrasekhar-mass delayed-detonation models that are outcomes of hydrodynamical simulations. We adopt spherical symmetry for all four models. The results of the different codes, including the light curves, spectra, and the evolution of several physical properties as a function of radius and time, are provided in electronic form in a standard format via a public repository. We also include the detailed test model profiles and several python scripts for accessing and presenting the input and output files. We also provide the code used to generate the toy models studied here. In this paper, we describe in detail the test models, radiative-transfer codes and output formats and provide access to the repository. We present example results of several key diagnostic features. △ Less

Submitted 15 April, 2023; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: Accepted for publication in A&A. 27 pages, 12 figures (v4: updated to match published version). The ejecta models and output files from the simulations are available at https://github.com/sn-rad-trans/data1

Journal ref: A&A 668, A163 (2022)

Searching for a Hypervelocity White Dwarf Companion: A Proper Motion Survey of SN 1006

Authors: Joshua V. Shields, Wolfgang Kerzendorf, Matthew W. Hosek Jr., Ken J. Shen, Armin Rest, Tuan Do, Jessica R. Lu, Andrew G. Fullard, Giovanni Strampelli, Alfredo Zenteno

Abstract: Type Ia Supernovae (SNe Ia) are securely understood to come from the thermonuclear explosion of a white dwarf as a result of binary interaction, but the nature of that binary interaction and the secondary object is uncertain. Recently, a double white dwarf model known as the dynamically driven double-degenerate double-detonation (D6) model has become a promising explanation for these events. One r… ▽ More Type Ia Supernovae (SNe Ia) are securely understood to come from the thermonuclear explosion of a white dwarf as a result of binary interaction, but the nature of that binary interaction and the secondary object is uncertain. Recently, a double white dwarf model known as the dynamically driven double-degenerate double-detonation (D6) model has become a promising explanation for these events. One realization of this scenario predicts that the companion may survive the explosion and reside within the remnant as a fast moving ($V_{peculiar} >1000$ km s$^{-1}$), overluminous ($L > 0.1 L_\odot$) white dwarf. Recently, three objects which appear to have these unusual properties have been discovered in the Gaia survey. We obtained photometric observations of the SN Ia remnant SN 1006 with the Dark Energy Camera over four years to attempt to discover a similar star. We present a deep, high precision astrometric proper motion survey of the interior stellar population of the remnant. We rule out the existence of a high proper motion object consistent with our tested realization of the D6 scenario ($V_{transverse} > 600$ km s$^{-1}$ with $m_r < 21$ corresponding to an intrinsic luminosity of $L > 0.0176 L_\odot$). We conclude that such a star does not exist within the remnant, or is hidden from detection by either strong localized dust or the unlikely possibility of ejection from the binary system near parallel to the line of sight. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: 15 pages, 10 figures

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

Binary Interaction Dominates Mass Ejection in Classical Novae

Authors: Ken J. Shen, Eliot Quataert

Abstract: Recent observations suggest our understanding of mass loss in classical novae is incomplete, motivating a new theoretical examination of the physical processes responsible for nova mass ejection. In this paper, we perform hydrodynamical simulations of classical nova outflows using the stellar evolution code MESA. We find that, when the binary companion is neglected, white dwarfs with masses >= 0.8… ▽ More Recent observations suggest our understanding of mass loss in classical novae is incomplete, motivating a new theoretical examination of the physical processes responsible for nova mass ejection. In this paper, we perform hydrodynamical simulations of classical nova outflows using the stellar evolution code MESA. We find that, when the binary companion is neglected, white dwarfs with masses >= 0.8 Msol successfully launch radiation-pressure-driven optically thick winds that carry away most of the envelope. However, for most of the mass loss phase, these winds are accelerated at radii beyond the white dwarf's Roche radius assuming a typical cataclysmic variable donor. This means that, before a standard optically thick wind can be formed, mass loss will instead be initiated and shaped by the binary interaction. An isotropic optically thick wind is only successfully launched when the acceleration region recedes within the white dwarf's Roche radius, which occurs after most of the envelope has already been ejected. The interaction between these two modes of outflow -- a first phase of slow, binary-driven, equatorially focused mass loss encompassing most of the mass ejection and a second phase consisting of a fast, isotropic, optically thick wind -- is consistent with observations of aspherical ejecta and signatures of multiple outflow components. We also find that isolated lower-mass white dwarfs <= 0.8 Msol do not develop unbound optically thick winds at any stage, making it even more crucial to consider the effects of the binary companion on the resulting outburst. △ Less

Submitted 12 September, 2022; v1 submitted 12 May, 2022; originally announced May 2022.

Comments: Accepted for publication in ApJ. Minor changes to previous version for clarity

Modelling the ionisation state of Type Ia supernovae in the nebular-phase

Authors: Luke J. Shingles, Andreas Flörs, Stuart A. Sim, Christine E. Collins, Friedrich K. Roepke, Ivo R. Seitenzahl, Ken J. Shen

Abstract: The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a me… ▽ More The nebular spectra of Type Ia supernovae ($\gtrapprox$ 100 days after explosion) consist mainly of emission lines from singly- and doubly-ionised Fe-group nuclei. However, theoretical models for many scenarios predict that non-thermal ionisation leads to multiply-ionised species whose recombination photons ionise and deplete Fe$^{+}$ , resulting in negligible [Fe II] emission. We investigate a method to determine the collisional excitation conditions from [Fe II] line ratios independently from the ionisation state and find that it cannot be applied to highly-ionised models due to the influence of recombination cascades on Fe$^{+}$ level populations. When the ionisation state is artificially lowered, the line ratios (and excitation conditions) are too similar to distinguish between explosion scenarios. We investigate changes to the treatment of non-thermal energy deposition as a way to reconcile over-ionised theoretical models with observations and find that a simple work function approximation provides closer agreement with the data for sub-Mch models than a detailed Spencer-Fano treatment with widely-used cross section data. To quantify the magnitude of additional heating processes that would be required to sufficiently reduce ionisation from fast leptons, we artificially boost the rate of energy loss to free electrons. We find that the equivalent of as much as an eight times increase to the plasma loss rate would be needed to reconcile the sub-Mch model with observed spectra. Future studies could distinguish between reductions in the non-thermal ionisation rates and increased recombination rates, such as by clumping. △ Less

Submitted 30 March, 2022; originally announced March 2022.

Comments: 13 pages. Accepted by MNRAS

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

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

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

Submitted 7 March, 2022; originally announced March 2022.

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

X-ray emission from candidate stellar merger remnant TYC 2597-735-1 and its Blue Ring Nebula

Authors: Hans Moritz Günther, Keri Hoadley, Maximilian N. Günther, Brian D. Metzger, P. C. Schneider, Ken J. Shen

Abstract: Tight binary or multiple star systems can interact through mass transfer and follow vastly different evolutionary pathways than single stars. The star TYC 2597-735-1 is a candidate for a recent stellar merger remnant resulting from a coalescence of a low-mass companion with a primary star a few thousand years ago. This violent event is evident in a conical outflow ("Blue Ring Nebula") emitting in… ▽ More Tight binary or multiple star systems can interact through mass transfer and follow vastly different evolutionary pathways than single stars. The star TYC 2597-735-1 is a candidate for a recent stellar merger remnant resulting from a coalescence of a low-mass companion with a primary star a few thousand years ago. This violent event is evident in a conical outflow ("Blue Ring Nebula") emitting in UV light and surrounded by leading shock filaments observed in H$α$ and UV emission. From Chandra data, we report the detection of X-ray emission from the location of TYC 2597-735-1 with a luminosity $\log(L_\mathrm{X}/L_\mathrm{bol})=-5.5$. Together with a previously reported period around 14~days, this indicates ongoing stellar activity and the presence of strong magnetic fields on TYC 2597-735-1. Supported by stellar evolution models of merger remnants, we interpret the inferred stellar magnetic field as dynamo action associated with a newly formed convection zone in the atmosphere of TYC 2597-735-1, though internal shocks at the base of an accretion-powered jet cannot be ruled out. We speculate that this object will evolve into an FK Com type source, i.e. a class of rapidly spinning magnetically active stars for which a merger origin has been proposed but for which no relic accretion or large-scale nebula remains visible. We also detect likely X-ray emission from two small regions close to the outer shock fronts in the Blue Ring Nebula, which may arise from either inhomogenities in the circumstellar medium or in the mass and velocity distribution in the merger-driven outflow. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: accepted by AJ, 15 pages, 5 figures

Impacts of Nuclear Reaction Rate Uncertainties on Sub-Chandrasekhar-Mass White Dwarf Detonations

Authors: Thomas Fitzpatrick, Ken J. Shen

Abstract: The precise nature of Type Ia supernova (SN Ia) progenitors remains a mystery, but the detonation of a sub-Chandrasekhar-mass white dwarf (WD) has become a promising candidate. There is a growing body of work suggesting that the carbon core detonation of a sub-Chandrasekhar-mass WD can be triggered by the detonation of a helium shell accreted from a companion WD, through either inward shock conver… ▽ More The precise nature of Type Ia supernova (SN Ia) progenitors remains a mystery, but the detonation of a sub-Chandrasekhar-mass white dwarf (WD) has become a promising candidate. There is a growing body of work suggesting that the carbon core detonation of a sub-Chandrasekhar-mass WD can be triggered by the detonation of a helium shell accreted from a companion WD, through either inward shock convergence near the center or direct edge-lit detonation. This "double-detonation" SN Ia can be triggered by a small helium shell and is therefore well approximated by the detonation of a bare carbon-oxygen white dwarf (C/O WD). The impacts of uncertainties in experimentally and theoretically determined nuclear reaction rates on nucleosynthesis in the detonations of sub-Chandrasekhar-mass WDs have not yet been fully explored. We investigate the sensitivity of this model to nuclear reaction rate uncertainties to better constrain the nucleosynthetic yields resulting from these phenomena and identify the nuclear reaction rates whose uncertainties have the most significant impacts on nucleosynthesis. We find that the chemical abundances, and in particular those of the iron-group elements, are relatively insensitive to these nuclear reaction rate uncertainties. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Comments: 10 pages, 3 figures, submitted

Masses of White Dwarf Binary Companions to Type Ia Supernovae Measured from Runaway Velocities

Authors: Evan B. Bauer, Vedant Chandra, Ken J. Shen, J. J. Hermes

Abstract: The recently proposed "dynamically driven double-degenerate double-detonation" (D6) scenario posits that Type Ia supernovae (SNe) may occur during dynamically unstable mass transfer between two white dwarfs (WDs) in a binary. This scenario predicts that the donor WD may then survive the explosion and be released as a hypervelocity runaway, opening up the exciting possibility of identifying remnant… ▽ More The recently proposed "dynamically driven double-degenerate double-detonation" (D6) scenario posits that Type Ia supernovae (SNe) may occur during dynamically unstable mass transfer between two white dwarfs (WDs) in a binary. This scenario predicts that the donor WD may then survive the explosion and be released as a hypervelocity runaway, opening up the exciting possibility of identifying remnant stars from D6 SNe and using them to study the physics of detonations that produce Type Ia SNe. Three candidate D6 runaway objects have been identified in Gaia data. The observable runaway velocity of these remnant objects represents their orbital speed at the time of SN detonation. The orbital dynamics and Roche lobe geometry required in the D6 scenario place specific constraints on the radius and mass of the donor WD that becomes the hypervelocity runaway. In this letter, we calculate the radii required for D6 donor WDs as a function of the runaway velocity. Using mass-radius relations for WDs, we then constrain the masses of the donor stars as well. With measured velocities for each of the three D6 candidate objects based on Gaia EDR3, this work provides a new probe of the masses and mass ratios in WD binary systems that produce SN detonations and hypervelocity runaways. △ Less

Submitted 14 December, 2021; v1 submitted 6 December, 2021; originally announced December 2021.

Comments: 8 pages, 4 figures, accepted for publication in ApJL

The SNIa Runaway LP 398-9: Detection of Circumstellar Material and Surface Rotation

Authors: Vedant Chandra, Hsiang-Chih Hwang, Nadia L. Zakamska, Simon Blouin, Andrew Swan, Thomas R. Marsh, Ken J. Shen, Boris T. Gänsicke, J. J. Hermes, Odelia Putterman, Evan B. Bauer, Evan Petrosky, Vikram S. Dhillon, Stuart P. Littlefair, Richard P. Ashley

Abstract: A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back $\approx 10^5$ years to a known supernova r… ▽ More A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back $\approx 10^5$ years to a known supernova remnant. Here we report the discovery of carbon-rich circumstellar material around LP 398-9, revealed by a strong infrared excess and analyzed with follow-up spectroscopy. The circumstellar material is most plausibly composed of inflated layers from the star itself, mechanically and radioactively heated by the past companion's supernova. We also detect a 15.4 hr periodic signal in the UV and optical light curves of LP 398-9, which we interpret as surface rotation. The rotation rate is consistent with theoretical predictions from this supernova mechanism, and the brightness variations could originate from surface inhomogeneity deposited by the supernova itself. Our observations strengthen the case for this double-degenerate SNIa progenitor channel, and motivate the search for more runaway SNIa donors. △ Less

Submitted 13 April, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: Accepted to MNRAS

Multi-Dimensional Radiative Transfer Calculations of Double Detonations of Sub-Chandrasekhar-Mass White Dwarfs

Authors: Ken J. Shen, Samuel J. Boos, Dean M. Townsley, Daniel Kasen

Abstract: Study of the double detonation Type Ia supernova scenario, in which a helium shell detonation triggers a carbon core detonation in a sub-Chandrasekhar-mass white dwarf, has experienced a resurgence in the past decade. New evolutionary scenarios and a better understanding of which nuclear reactions are essential have allowed for successful explosions in white dwarfs with much thinner helium shells… ▽ More Study of the double detonation Type Ia supernova scenario, in which a helium shell detonation triggers a carbon core detonation in a sub-Chandrasekhar-mass white dwarf, has experienced a resurgence in the past decade. New evolutionary scenarios and a better understanding of which nuclear reactions are essential have allowed for successful explosions in white dwarfs with much thinner helium shells than in the original, decades-old incarnation of the double detonation scenario. In this paper, we present the first suite of light curves and spectra from multi-dimensional radiative transfer calculations of thin-shell double detonation models, exploring a range of white dwarf and helium shell masses. We find broad agreement with the observed light curves and spectra of non-peculiar Type Ia supernovae, from subluminous to overluminous subtypes, providing evidence that double detonations of sub-Chandrasekhar-mass white dwarfs produce the bulk of observed Type Ia supernovae. Some discrepancies in spectral velocities and colors persist, but these may be brought into agreement by future calculations that include more accurate initial conditions and radiation transport physics. △ Less

Submitted 31 August, 2021; v1 submitted 27 August, 2021; originally announced August 2021.

Comments: Accepted for publication in ApJ. Same as previous arXiv version

Transients from the Cataclysmic Deaths of Cataclysmic Variables

Authors: Brian D. Metzger, Yossef Zenati, Laura Chomiuk, Ken J. Shen, Jay Strader

Abstract: We explore the observational appearance of the merger of a low-mass star with a white dwarf (WD) binary companion. We are motivated by Schreiber et al. (2016), who found that multiple tensions between the observed properties of cataclysmic variables (CVs) and standard evolution models are resolved if a large fraction of CV binaries merge as a result of unstable mass transfer. Tidal disruption of t… ▽ More We explore the observational appearance of the merger of a low-mass star with a white dwarf (WD) binary companion. We are motivated by Schreiber et al. (2016), who found that multiple tensions between the observed properties of cataclysmic variables (CVs) and standard evolution models are resolved if a large fraction of CV binaries merge as a result of unstable mass transfer. Tidal disruption of the secondary forms a geometrically thick disk around the WD, which subsequently accretes at highly super-Eddington rates. Analytic estimates and numerical hydrodynamical simulations reveal that outflows from the accretion flow unbind a large fraction >~ 90% of the secondary at velocities ~500-1000 km/s within days of the merger. Hydrogen recombination in the expanding ejecta powers optical transient emission lasting about a month with a luminosity > 1e38 erg/s, similar to slow classical novae and luminous red novae from ordinary stellar mergers. Over longer timescales the mass accreted by the WD undergoes hydrogen shell burning, inflating the remnant into a giant of luminosity ~300-5000 L_sun, effective temperature T_eff ~ 3000 K and lifetime ~1e4-1e5 yr. We predict that ~1e3-1e4 Milky Way giants are CV merger products, potentially distinguishable by atypical surface abundances. We explore whether any Galactic historical slow classical novae are masquerading CV mergers by identifying four such post-nova systems with potential giant counterparts for which a CV merger origin cannot be ruled out. We address whether the historical transient CK Vul and its gaseous/dusty nebula resulted from a CV merger. △ Less

Submitted 18 November, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 15 pages, 5 figures, accepted to ApJ after minor edits and errors corrected

Birth of the ELMs: a ZTF survey for evolved cataclysmic variables turning into extremely low-mass white dwarfs

Authors: Kareem El-Badry, Hans-Walter Rix, Eliot Quataert, Thomas Kupfer, Ken J. Shen

Abstract: We present a systematic survey for mass-transferring and recently-detached cataclysmic variables (CVs) with evolved secondaries, which are progenitors of extremely low-mass white dwarfs (ELM WDs), AM CVn systems, and detached ultracompact binaries. We select targets below the main sequence in the Gaia color-magnitude diagram with ZTF light curves showing large-amplitude ellipsoidal variability and… ▽ More We present a systematic survey for mass-transferring and recently-detached cataclysmic variables (CVs) with evolved secondaries, which are progenitors of extremely low-mass white dwarfs (ELM WDs), AM CVn systems, and detached ultracompact binaries. We select targets below the main sequence in the Gaia color-magnitude diagram with ZTF light curves showing large-amplitude ellipsoidal variability and orbital period $P_{\rm orb} < 6$ hr. This yields 51 candidates brighter than G=18, of which we have obtained many-epoch spectra for 21. We confirm all 21 to be completely -- or nearly -- Roche lobe filling close binaries. 13 show evidence of ongoing mass transfer, which has likely just ceased in the other 8. Most of the secondaries are hotter than any previously known CV donors, with temperatures $4700<T_{\rm eff}/{\rm K}<8000$. Remarkably, all secondaries with $T_{\rm eff} \gtrsim 7000\,\rm K$ appear to be detached, while all cooler secondaries are still mass-transferring. This transition likely marks the temperature where magnetic braking becomes inefficient due to loss of the donor's convective envelope. Most of the proto-WD secondaries have masses near $0.15\,M_{\odot}$; their companions have masses near $0.8\,M_{\odot}$. We infer a space density of $\sim 60\,\rm kpc^{-3}$, roughly 80 times lower than that of normal CVs and three times lower than that of ELM WDs. The implied Galactic birth rate, $\mathcal{R}\sim 60\,\rm Myr^{-1}$, is half that of AM CVn binaries. Most systems are well-described by MESA models for CVs in which mass transfer begins only as the donor leaves the main sequence. All are predicted to reach minimum periods $5\lesssim P_{\rm orb}/{\rm min}\lesssim30$ within a Hubble time, where they will become AM CVn binaries or merge. This sample triples the known evolved CV population and offers broad opportunities for improving understanding of the compact binary population. △ Less

Submitted 30 November, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 26 pages, 16 figures, plus appendices. Magnitudes in Table 1 fixed; machine-readable table included in source files. Accepted to MNRAS

8.9-hr Rotation in the Partly Burnt Runaway Stellar Remnant LP 40-365 (GD 492)

Authors: J. J. Hermes, Odelia Putterman, Mark A. Hollands, David J. Wilson, Andrew Swan, Roberto Raddi, Ken J. Shen, Boris T. Gaensicke

Abstract: We report the detection of 8.914-hr variability in both optical and ultraviolet light curves of LP 40-365 (also known as GD 492), the prototype for a class of partly burnt runaway stars that have been ejected from a binary due to a thermonuclear supernova event. We first detected this 1.0% amplitude variation in optical photometry collected by the Transiting Exoplanet Survey Satellite. Re-analysis… ▽ More We report the detection of 8.914-hr variability in both optical and ultraviolet light curves of LP 40-365 (also known as GD 492), the prototype for a class of partly burnt runaway stars that have been ejected from a binary due to a thermonuclear supernova event. We first detected this 1.0% amplitude variation in optical photometry collected by the Transiting Exoplanet Survey Satellite. Re-analysis of observations from the Hubble Space Telescope at the TESS period and ephemeris reveal a 5.8% variation in the ultraviolet of this 9800 K stellar remnant. We propose that this 8.914-hr photometric variation reveals the current surface rotation rate of LP 40-365, and is caused by some kind of surface inhomogeneity rotating in and out of view, though a lack of observed Zeeman splitting puts an upper limit on the magnetic field of <20 kG. We explore ways in which the present rotation period can constrain progenitor scenarios if angular momentum was mostly conserved, which suggests that the survivor LP 40-365 was not the donor star but was most likely the bound remnant of a mostly disrupted white dwarf that underwent advanced burning from an underluminous (Type Iax) supernova. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: 6 pages, 2 figures, accepted for publication in ApJ Letters

LAMOST J0140355+392651: An evolved cataclysmic variable donor transitioning to become an extremely low mass white dwarf

Authors: Kareem El-Badry, Eliot Quataert, Hans-Walter Rix, Daniel R. Weisz, Thomas Kupfer, Ken J. Shen, Maosheng Xiang, Yong Yang, Xiaowei Liu

Abstract: We present LAMOST J0140355+392651 (hereafter J0140), a close ($P_{\rm orb} = 3.81$ hours) binary containing a bloated, low-mass ($M \approx 0.15 M_{\odot}$) proto-white dwarf (WD) and a massive ($M\approx 0.95\,M_{\odot}$) WD companion. The system's optical light curve is dominated by large-amplitude ellipsoidal variability but also exhibits additional scatter, likely driven by pulsations. The pro… ▽ More We present LAMOST J0140355+392651 (hereafter J0140), a close ($P_{\rm orb} = 3.81$ hours) binary containing a bloated, low-mass ($M \approx 0.15 M_{\odot}$) proto-white dwarf (WD) and a massive ($M\approx 0.95\,M_{\odot}$) WD companion. The system's optical light curve is dominated by large-amplitude ellipsoidal variability but also exhibits additional scatter, likely driven by pulsations. The proto-WD is cooler ($T_{\rm eff} = 6800\pm 100$ K) and more puffy ($\log\left[g/\left({\rm cm\,s^{-2}}\right)\right]=4.74\pm0.07$) than any known extremely low mass (ELM) WD, but hotter than any known cataclysmic variable (CV) donor. It either completely or very nearly fills its Roche lobe ($R/R_{\rm Roche\,lobe}=0.99\pm0.01$), suggesting ongoing or recently terminated mass transfer. No dwarf nova-like outbursts have been observed. The spectrum is dominated by the proto-WD but shows tentative hints of H$α$ emission, perhaps due to accretion onto the massive WD. The properties of the system are well-matched by MESA binary evolution models of CVs with donors that underwent significant nuclear evolution before the onset of mass transfer. In these models, the bloated proto-WD is either still losing mass via stable Roche lobe overflow or was doing so until very recently. In either case, it is evolving toward higher temperatures at near-constant luminosity to become an ELM WD. If the system is detached, mass transfer likely ended when the donor became too hot for magnetic braking to remain efficient. Evolutionary models predict that the binary will shrink to $P_{\rm orb}\lesssim 10$ minutes within a few Gyr, when it will either merge or become an AM CVn binary. J0140 provides an observational link between the formation channels of CVs, ELM WDs, detached ultracompact WD binaries, and AM CVn systems. △ Less

Submitted 6 May, 2021; v1 submitted 14 April, 2021; originally announced April 2021.

Comments: 19 pages, 14 figures. Accepted to MNRAS

Non-Local Thermodynamic Equilibrium Radiative Transfer Simulations of Sub-Chandrasekhar-Mass White Dwarf Detonations

Authors: Ken J. Shen, Stéphane Blondin, Daniel Kasen, Luc Dessart, Dean M. Townsley, Samuel Boos, D. John Hillier

Abstract: Type Ia supernovae (SNe Ia) span a range of luminosities and timescales, from rapidly evolving subluminous to slowly evolving overluminous subtypes. Previous theoretical work has, for the most part, been unable to match the entire breadth of observed SNe Ia with one progenitor scenario. Here, for the first time, we apply non-local thermodynamic equilibrium radiative transfer calculations to a rang… ▽ More Type Ia supernovae (SNe Ia) span a range of luminosities and timescales, from rapidly evolving subluminous to slowly evolving overluminous subtypes. Previous theoretical work has, for the most part, been unable to match the entire breadth of observed SNe Ia with one progenitor scenario. Here, for the first time, we apply non-local thermodynamic equilibrium radiative transfer calculations to a range of accurate explosion models of sub-Chandrasekhar-mass white dwarf detonations. The resulting photometry and spectra are in excellent agreement with the range of observed non-peculiar SNe Ia through 15 d after the time of B-band maximum, yielding one of the first examples of a quantitative match to the entire Phillips (1993) relation. The intermediate-mass element velocities inferred from theoretical spectra at maximum light for the more massive white dwarf explosions are higher than those of bright observed SNe Ia, but these and other discrepancies likely stem from the one-dimensional nature of our explosion models and will be improved upon by future non-local thermodynamic equilibrium radiation transport calculations of multi-dimensional sub-Chandrasekhar-mass white dwarf detonations. △ Less

Submitted 7 December, 2023; v1 submitted 16 February, 2021; originally announced February 2021.

Comments: Published in ApJL. This version corrects the previous version, which erroneously contained the R-band light curve of SN 1999dq in the I-band panels of Figures 2 and 6

Multi-Dimensional Parameter Study of Double Detonation Type Ia Supernovae Originating from Thin-Helium-Shell White Dwarfs

Authors: Samuel J. Boos, Dean M. Townsley, Ken J. Shen, Spencer Caldwell, Broxton J. Miles

Abstract: Despite the importance of Type Ia supernovae (SNe Ia) throughout astronomy, the precise progenitor systems and explosion mechanisms that drive SNe Ia are still unknown. An explosion scenario that has gained traction recently is the double detonation in which an accreted shell of He detonates and triggers a secondary detonation in the underlying white dwarf. Our research presents a number of high r… ▽ More Despite the importance of Type Ia supernovae (SNe Ia) throughout astronomy, the precise progenitor systems and explosion mechanisms that drive SNe Ia are still unknown. An explosion scenario that has gained traction recently is the double detonation in which an accreted shell of He detonates and triggers a secondary detonation in the underlying white dwarf. Our research presents a number of high resolution, multi-dimensional, full star simulations of thin-He-shell, sub-Chandrasekhar-mass white dwarf progenitors that undergo a double detonation. This suite of thin-shell progenitors incorporates He shells that are thinner than those in previous multi-dimensional studies. We confirm the viability of the double detonation across a range of He shell parameter space as well as present bulk yields and ejecta profiles for each progenitor. The yields obtained are generally consistent with previous works and indicate the likelihood of producing observables that resemble SNe Ia. The dimensionality of our simulations allow us to examine features of the double detonation more closely, including the details of the off-center secondary ignition and asymmetric ejecta. We find considerable differences in the high-velocity extent of post-detonation products across different lines of sight. The data from this work will be used to generate predicted observables and may further support the viability of the double detonation scenario as a SNe Ia channel as well as show how properties of the progenitor or viewing angle may influence trends in observable characteristics. △ Less

Submitted 3 June, 2021; v1 submitted 28 January, 2021; originally announced January 2021.

Comments: 20 pages, 9 figures, 11 tables, accepted for publication in ApJ

A "blue ring nebula" surrounding a thousands of years old stellar merger

Authors: Keri Hoadley, Christopher Martin, Brian D. Metzger, Mark Seibert, Andrew McWilliam, Ken J. Shen, James D. Neill, Gudmundur Stefansson, Andrew Monson, Bradley E. Schaefer

Abstract: Stellar mergers are a brief-lived but common phase in the evolution of binary star systems. Among the many astrophysical implications of these events include the creation of atypical stars (e.g. magnetic stars, blue stragglers, rapid rotators), interpretation of stellar populations, and formation channels of LIGO-detected compact object mergers. Although stellar mergers are thus commonly invoked p… ▽ More Stellar mergers are a brief-lived but common phase in the evolution of binary star systems. Among the many astrophysical implications of these events include the creation of atypical stars (e.g. magnetic stars, blue stragglers, rapid rotators), interpretation of stellar populations, and formation channels of LIGO-detected compact object mergers. Although stellar mergers are thus commonly invoked phenomena, observations of these events and details of their evolution remain elusive. While a handful of stellar mergers have been directly observed in recent years, the central remnants of these events remain shrouded by an opaque shell of dust and molecules, making it impossible to observe their final state (e.g. as a single merged star or a tighter surviving binary). Here we report observations of an unusual, ring-shaped ultraviolet nebula and the star at its center, TYC 2597-735-1. The nebula shows two opposing fronts, suggesting a bipolar outflow from TYC 2597-735-1. TYC 2597-735-1's spectrum and proximity above the Galactic plane suggest it is an old star, yet it shows abnormally low surface gravity and a detectable long-term luminosity decay, uncharacteristic for its evolutionary stage. TYC 2597-735-1 also exhibits H-alpha emission, radial velocity variations, enhanced ultraviolet radiation, and excess infrared emission, common signposts of dusty circumstellar disks, stellar activity, and accretion. The combined observations, paired with stellar evolution models, suggests TYC 2597-735-1 merged with a lower-mass companion several thousand years ago. TYC 2597-735-1 provides a look at an unobstructed stellar merger found at an evolutionary stage between its dynamic onset and the theorized final equilibrium state, directly inferring how two stars merge into a single star. △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: 25 pages (Main article), 3 Figures and 8 Extended Data Figures, SI (57 pages) included after Extended Data Figures. To appear in Nov 19 publication of Nature. Online release URL (Nov 18): https://www.nature.com/articles/s41586-020-2893-5

New Insights into Classical Novae

Authors: Laura Chomiuk, Brian D. Metzger, Ken J. Shen

Abstract: We survey our understanding of classical novae: non-terminal, thermonuclear eruptions on the surfaces of white dwarfs in binary systems. The recent and unexpected discovery of GeV gamma-rays from Galactic novae has highlighted the complexity of novae and their value as laboratories for studying shocks and particle acceleration. We review half a century of nova literature through this new lens, and… ▽ More We survey our understanding of classical novae: non-terminal, thermonuclear eruptions on the surfaces of white dwarfs in binary systems. The recent and unexpected discovery of GeV gamma-rays from Galactic novae has highlighted the complexity of novae and their value as laboratories for studying shocks and particle acceleration. We review half a century of nova literature through this new lens, and conclude: --The basics of the thermonuclear runaway theory of novae are confirmed by observations. The white dwarf sustains surface nuclear burning for some time after runaway, and until recently, it was commonly believed that radiation from this nuclear burning solely determines the nova's bolometric luminosity. --The processes by which novae eject material from the binary system remain poorly understood. Mass loss from novae is complex (sometimes fluctuating in rate, velocity, and morphology) and often prolonged in time over weeks, months, or years. --The complexity of the mass ejection leads to gamma-ray producing shocks internal to the nova ejecta. When gamma-rays are detected (around optical maximum), the shocks are deeply embedded and the surrounding gas is very dense. --Observations of correlated optical and gamma-ray light curves confirm that the shocks are radiative and contribute significantly to the bolometric luminosity of novae. Novae are therefore the closest and most common "interaction-powered" transients. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: Accepted to Annual Reviews of Astronomy and Astrophysics. 59 pages, 14 figures. Supplementary materials included in this document

Distribution of Si II $λ$6355 Velocities of Type Ia Supernovae and Implications for Asymmetric Explosions

Authors: Keto D. Zhang, WeiKang Zheng, Thomas de Jaeger, Benjamin E. Stahl, Thomas G. Brink, Xuhui Han, Daniel Kasen, Ken J. Shen, Kevin Tang, Alexei V. Filippenko

Abstract: The ejecta velocity is a very important parameter in studying the structure and properties of Type Ia supernovae (SNe Ia). It is also a candidate key parameter in improving the utility of SNe Ia for cosmological distance determinations. Here we study the velocity distribution of a sample of 311 SNe Ia from the kaepora database. The velocities are derived from the Si II $λ$6355 absorption line in o… ▽ More The ejecta velocity is a very important parameter in studying the structure and properties of Type Ia supernovae (SNe Ia). It is also a candidate key parameter in improving the utility of SNe Ia for cosmological distance determinations. Here we study the velocity distribution of a sample of 311 SNe Ia from the kaepora database. The velocities are derived from the Si II $λ$6355 absorption line in optical spectra measured at (or extrapolated to) the time of peak brightness. We statistically show that the observed velocity has a bimodal Gaussian distribution consisting of two groups of SNe Ia: Group I with a lower but narrower scatter ($μ_1 = 11000 \text{km s}^{-1}$, $σ_1 = 700 \text{km s}^{-1}$), and Group II with a higher but broader scatter ($μ_2 = 12300 \text{km s}^{-1}$, $σ_2 = 1800 \text{km s}^{-1}$). The population ratio of Group I to Group II is 201:110 (65%:35%). There is substantial degeneracy between the two groups, but for SNe Ia with velocity $v > 12000 \text{km s}^{-1}$, the distribution is dominated by Group II. The true origin of the two components is unknown, though there could be that naturally there exist two intrinsic velocity distributions as observed. However, we try to use asymmetric geometric models through statistical simulations to reproduce the observed distribution assuming all SNe Ia share the same intrinsic distribution. In the two cases we consider, 35\% of SNe Ia are considered to be asymmetric in Case 1, and all SNe Ia are asymmetric in Case 2. Simulations for both cases can reproduce the observed velocity distribution but require a significantly large portion ($>35\%$) of SNe Ia to be asymmetric. In addition, the Case 1 result is consistent with recent polarization observations that SNe Ia with higher Si II $λ$6355 velocity tend to be more polarized. △ Less

Submitted 11 October, 2020; originally announced October 2020.

Comments: Accepted for publication in MNRAS

SN 2019ehk: A Double-Peaked Ca-rich Transient with Luminous X-ray Emission and Shock-Ionized Spectral Features

Authors: Wynn V. Jacobson-Galán, Raffaella Margutti, Charles D. Kilpatrick, Daichi Hiramatsu, Hagai Perets, David Khatami, Ryan J. Foley, John Raymond, Sung-Chul Yoon, Alexey Bobrick, Yossef Zenati, Lluís Galbany, Jennifer Andrews, Peter J. Brown, Régis Cartier, Deanne L. Coppejans, Georgios Dimitriadis, Matthew Dobson, Aprajita Hajela, D. Andrew Howell, Hanindyo Kuncarayakti, Danny Milisavljevic, Mohammed Rahman, César Rojas-Bravo, David J. Sand , et al. (42 additional authors not shown)

Abstract: We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission… ▽ More We present panchromatic observations and modeling of the Calcium-rich supernova 2019ehk in the star-forming galaxy M100 (d$\approx$16.2 Mpc) starting 10 hours after explosion and continuing for ~300 days. SN 2019ehk shows a double-peaked optical light curve peaking at $t = 3$ and $15$ days. The first peak is coincident with luminous, rapidly decaying $\textit{Swift}$-XRT discovered X-ray emission ($L_x\approx10^{41}~\rm{erg~s^{-1}}$ at 3 days; $L_x \propto t^{-3}$), and a Shane/Kast spectral detection of narrow H$α$ and He II emission lines ($v \approx 500$ km/s) originating from pre-existent circumstellar material. We attribute this phenomenology to radiation from shock interaction with extended, dense material surrounding the progenitor star at $r<10^{15}$ cm and the resulting cooling emission. We calculate a total CSM mass of $\sim$ $7\times10^{-3}$ $\rm{M_{\odot}}$ with particle density $n\approx10^{9}\,\rm{cm^{-3}}$. Radio observations indicate a significantly lower density $n < 10^{4}\,\rm{cm^{-3}}$ at larger radii. The photometric and spectroscopic properties during the second light curve peak are consistent with those of Ca-rich transients (rise-time of $t_r =13.4\pm0.210$ days and a peak B-band magnitude of $M_B =-15.1\pm0.200$ mag). We find that SN 2019ehk synthesized $(3.1\pm0.11)\times10^{-2} ~ \rm{M_{\odot}}$ of ${}^{56}\textrm{Ni}$ and ejected $M_{\rm ej} = (0.72\pm 0.040)~\rm{M_{\odot}}$ total with a kinetic energy $E_{\rm k}=(1.8\pm0.10)\times10^{50}~\rm{erg}$. Finally, deep $\textit{HST}$ pre-explosion imaging at the SN site constrains the parameter space of viable stellar progenitors to massive stars in the lowest mass bin (~10 $\rm{M_{\odot}}$) in binaries that lost most of their He envelope or white dwarfs. The explosion and environment properties of SN 2019ehk further restrict the potential WD progenitor systems to low-mass hybrid HeCO WD + CO WD binaries. △ Less

Submitted 1 July, 2020; v1 submitted 4 May, 2020; originally announced May 2020.

Comments: 51 pages, 27 figures. Accepted for publication in ApJ

Direct evidence for shock-powered optical emission in a nova

Authors: Elias Aydi, Kirill V. Sokolovsky, Laura Chomiuk, Elad Steinberg, Kwan Lok Li, Indrek Vurm, Brian D. Metzger, Jay Strader, Koji Mukai, Ondřej Pejcha, Ken J. Shen, Gregg A. Wade, Rainer Kuschnig, Anthony F. J. Moffat, Herbert Pablo, Andrzej Pigulski, Adam Popowicz, Werner Weiss, Konstanze Zwintz, Luca Izzo, Karen R. Pollard, Gerald Handler, Stuart D. Ryder, Miroslav D. Filipović, Rami Z. E. Alsaberi , et al. (17 additional authors not shown)

Abstract: Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical no… ▽ More Classical novae are thermonuclear explosions that occur on the surfaces of white dwarf stars in interacting binary systems (Bode & Evans 2008). It has long been thought that the luminosity of classical novae is powered by continued nuclear burning on the surface of the white dwarf after the initial runaway (Gallaher & Starrfield 1978). However, recent observations of GeV $γ$-rays from classical novae have hinted that shocks internal to the nova ejecta may dominate the nova emission. Shocks have also been suggested to power the luminosity of events as diverse as stellar mergers (Metzger & Pejcha 2017), supernovae (Moriya et al. 2018), and tidal disruption events (Roth et al. 2016), but observational confirmation has been lacking. Here we report simultaneous space-based optical and $γ$-ray observations of the 2018 nova V906 Carinae (ASASSN-18fv), revealing a remarkable series of distinct correlated flares in both bands. The optical and $γ$-ray flares occur simultaneously, implying a common origin in shocks. During the flares, the nova luminosity doubles, implying that the bulk of the luminosity is shock-powered. Furthermore, we detect concurrent but weak X-ray emission from deeply embedded shocks, confirming that the shock power does not appear in the X-ray band and supporting its emergence at longer wavelengths. Our data, spanning the spectrum from radio to $γ$-ray, provide direct evidence that shocks can power substantial luminosity in classical novae and other optical transients. △ Less

Submitted 12 April, 2020; originally announced April 2020.

Comments: 49 pages including supplementary information. Accepted and published in Nature Astronomy

Journal ref: Nature Astronomy 2020

Manganese Indicates a Transition from Sub- to Near-Chandrasekhar Type Ia Supernovae in Dwarf Galaxies

Authors: Mithi A. C. de los Reyes, Evan N. Kirby, Ivo R. Seitenzahl, Ken J. Shen

Abstract: Manganese (Mn) abundances are sensitive probes of the progenitors of Type Ia supernovae (SNe). In this work, we present a catalog of manganese abundances in dwarf spheroidal satellites of the Milky Way, measured using medium-resolution spectroscopy. Using a simple chemical evolution model, we infer the manganese yield of Type Ia SNe in the Sculptor dwarf spheroidal galaxy (dSph) and compare to the… ▽ More Manganese (Mn) abundances are sensitive probes of the progenitors of Type Ia supernovae (SNe). In this work, we present a catalog of manganese abundances in dwarf spheroidal satellites of the Milky Way, measured using medium-resolution spectroscopy. Using a simple chemical evolution model, we infer the manganese yield of Type Ia SNe in the Sculptor dwarf spheroidal galaxy (dSph) and compare to theoretical yields. The sub-solar yield from Type Ia SNe ($\mathrm{[Mn/Fe]}_{\mathrm{Ia}}=-0.30_{-0.03}^{+0.03}$ at $\mathrm{[Fe/H]}=-1.5$ dex, with negligible dependence on metallicity) implies that sub-Chandrasekhar-mass (sub-$M_{\mathrm{Ch}}$) white dwarf progenitors are the dominant channel of Type Ia SNe at early times in this galaxy, although some fraction ($\gtrsim20\%$) of $M_{\mathrm{Ch}}$ Type Ia or Type Iax SNe are still needed to produce the observed yield. However, this result does not hold in all environments. In particular, we find that dSph galaxies with extended star formation histories (Leo I, Fornax dSphs) appear to have higher [Mn/Fe] at a given metallicity than galaxies with early bursts of star formation (Sculptor dSph), suggesting that $M_{\mathrm{Ch}}$ progenitors may become the dominant channel of Type Ia SNe at later times in a galaxy's chemical evolution. △ Less

Submitted 5 February, 2020; v1 submitted 6 January, 2020; originally announced January 2020.

Comments: 24 pages, 8 figures. Accepted to ApJ

Journal ref: ApJ 891 (2020) 85

The Progenitors of Calcium-Strong Transients

Authors: Ken J. Shen, Eliot Quataert, Rüdiger Pakmor

Abstract: A new class of faint, spectroscopically peculiar transients has emerged in the last decade. We term these events "calcium-strong transients" (CaSTs) because of their atypically high calcium-to-oxygen nebular line ratios. Previous studies have struggled to deduce the identity of their progenitors due to a combination of their extremely extended radial distributions with respect to their host galaxi… ▽ More A new class of faint, spectroscopically peculiar transients has emerged in the last decade. We term these events "calcium-strong transients" (CaSTs) because of their atypically high calcium-to-oxygen nebular line ratios. Previous studies have struggled to deduce the identity of their progenitors due to a combination of their extremely extended radial distributions with respect to their host galaxies and their relatively high rate of occurrence. In this work, we find that the CaST radial distribution is consistent with the radial distribution of two populations of stars: old (ages > 5 Gyr), low-metallicity (Z/Zsol < 0.3) stars and globular clusters. While no obvious progenitor scenario arises from considering old, metal-poor stars, the alternative production site of globular clusters leads us to narrow down the list of possible candidates to three binary scenarios: mergers of helium and oxygen/neon white dwarfs; tidal disruptions of helium white dwarfs by neutron stars; and stable accretion from low-mass helium-burning stars onto white dwarfs. While rare in the field, these binary systems can be formed dynamically at much higher rates in globular clusters. Subsequent binary hardening both increases their interaction rate and ejects them from their parent globular clusters prior to mass transfer contact. Their production in, and ejection from, globular clusters may explain their radial distribution and the absence of globular clusters at their explosion site. This model predicts a currently undiscovered high rate of CaSTs in nuclear star clusters. Alternatively, an undetermined progenitor scenario involving old, low-metallicity stars may instead hold the key to understanding CaSTs. △ Less

Submitted 10 December, 2019; v1 submitted 21 August, 2019; originally announced August 2019.

Comments: Accepted to ApJ. Minor changes and updates to previous version. Main conclusions remain unchanged

Evidence for Sub-Chandrasekhar Type Ia Supernovae from Stellar Abundances in Dwarf Galaxies

Authors: Evan N. Kirby, Justin L. Xie, Rachel Guo, Mithi A. C. de los Reyes, Maria Bergemann, Mikhail Kovalev, Ken J. Shen, Anthony L. Piro, Andrew McWilliam

Abstract: There is no consensus on the progenitors of Type Ia supernovae (SNe Ia) despite their importance for cosmology and chemical evolution. We address this question by using our previously published catalogs of Mg, Si, Ca, Cr, Fe, Co, and Ni abundances in dwarf galaxy satellites of the Milky Way to constrain the mass at which the white dwarf explodes during a typical SN Ia. We fit a simple bi-linear mo… ▽ More There is no consensus on the progenitors of Type Ia supernovae (SNe Ia) despite their importance for cosmology and chemical evolution. We address this question by using our previously published catalogs of Mg, Si, Ca, Cr, Fe, Co, and Ni abundances in dwarf galaxy satellites of the Milky Way to constrain the mass at which the white dwarf explodes during a typical SN Ia. We fit a simple bi-linear model to the evolution of [X/Fe] with [Fe/H], where X represents each of the elements mentioned above. We use the evolution of [Mg/Fe] coupled with theoretical supernova yields to isolate what fraction of the elements originated in SNe Ia. Then, we infer the [X/Fe] yield of SNe Ia for all of the elements except Mg. We compare these observationally inferred yields to recent theoretical predictions for two classes of Chandrasekhar-mass (M_Ch) SN Ia as well as sub-M_Ch SNe Ia. Most of the inferred SN Ia yields are consistent with all of the theoretical models, but [Ni/Fe] is consistent only with sub-M_Ch models. We conclude that the dominant type of SN Ia in ancient dwarf galaxies is the explosion of a sub-M_Ch white dwarf. The Milky Way and dwarf galaxies with extended star formation histories have higher [Ni/Fe] abundances, which could indicate that the dominant class of SN Ia is different for galaxies where star formation lasted for at least several Gyr. △ Less

Submitted 22 July, 2019; v1 submitted 24 June, 2019; originally announced June 2019.

Comments: Accepted to ApJ on 2019 June 22. 21 pages, 9 figures. This is a companion paper to arXiv:1906:08284. Figure 6 presents the main result. Minor proof corrections in v2

Double Detonations with Thin, Modestly Enriched Helium Layers Can Make Normal Type Ia Supernovae

Authors: Dean M. Townsley, Broxton J. Miles, Ken J. Shen, Daniel Kasen

Abstract: It has been proposed that Type Ia supernovae (SNe Ia) that are normal in their spectra and brightness can be explained by a double detonation that ignites first in a helium shell on the surface of the white dwarf (WD). This proposition is supported by the satisfactory match between simulated explosions of sub-Chandrasekhar-mass WDs with no surface He layer and observations of normal SNe Ia. Howeve… ▽ More It has been proposed that Type Ia supernovae (SNe Ia) that are normal in their spectra and brightness can be explained by a double detonation that ignites first in a helium shell on the surface of the white dwarf (WD). This proposition is supported by the satisfactory match between simulated explosions of sub-Chandrasekhar-mass WDs with no surface He layer and observations of normal SNe Ia. However, previous calculations of He-ignited double detonations have required either artificial removal of the He shell ashes or extreme enrichment of the surface He layer in order to obtain normal SNe Ia. Here we demonstrate, for the first time in multi-dimensional full-star simulations, that a thin, modestly enriched He layer will lead to a SN Ia that is normal in its brightness and spectra. This strengthens the case for double detonations as a major contributing channel to the population of normal SNe Ia. △ Less

Submitted 19 July, 2019; v1 submitted 26 March, 2019; originally announced March 2019.

Comments: 11 pages, 4 figures, updated to version published in ApJ letters

Journal ref: 2019 ApJ, 878, L38

Partly burnt runaway stellar remnants from peculiar thermonuclear supernovae

Authors: R. Raddi, M. A. Hollands, D. Koester, J. J. Hermes, B. T. Gaensicke, U. Heber, K. J. Shen, D. M. Townsley, A. F. Pala, J. S. Reding, O. F. Toloza, I. Pelisoli, S. Geier, N. P. Gentile Fusillo, U. Munari, J. Strader

Abstract: We report the discovery of three stars that, along with the prototype LP40-365, form a distinct class of chemically peculiar runaway stars that are the survivors of thermonuclear explosions. Spectroscopy of the four confirmed LP 40-365 stars finds ONe-dominated atmospheres enriched with remarkably similar amounts of nuclear ashes of partial O- and Si-burning. Kinematic evidence is consistent with… ▽ More We report the discovery of three stars that, along with the prototype LP40-365, form a distinct class of chemically peculiar runaway stars that are the survivors of thermonuclear explosions. Spectroscopy of the four confirmed LP 40-365 stars finds ONe-dominated atmospheres enriched with remarkably similar amounts of nuclear ashes of partial O- and Si-burning. Kinematic evidence is consistent with ejection from a binary supernova progenitor; at least two stars have rest-frame velocities indicating they are unbound to the Galaxy. With masses and radii ranging between 0.20-0.28 Msun and 0.16-0.60 Rsun, respectively, we speculate these inflated white dwarfs are the partly burnt remnants of either peculiar Type Iax or electron-capture supernovae. Adopting supernova rates from the literature, we estimate that ~20 LP40-365 stars brighter than 19 mag should be detectable within 2 kpc from the Sun at the end of the Gaia mission. We suggest that as they cool, these stars will evolve in their spectroscopic appearance, and eventually become peculiar O-rich white dwarfs. Finally, we stress that the discovery of new LP40-365 stars will be useful to further constrain their evolution, supplying key boundary conditions to the modelling of explosion mechanisms, supernova rates, and nucleosynthetic yields of peculiar thermonuclear explosions. △ Less

Submitted 7 June, 2019; v1 submitted 13 February, 2019; originally announced February 2019.

Comments: 22 pages, 14 figures, 6 tables. Accepted for publication on MNRAS

Delayed Circumstellar Interaction for Type Ia SN 2015cp Revealed by an HST Ultraviolet Imaging Survey

Authors: M. L. Graham, C. E. Harris, P. E. Nugent, K. Maguire, M. Sullivan, M. Smith, S. Valenti, A. Goobar, O. D. Fox, K. J. Shen, P. L. Kelly, C. McCully, T. G. Brink, A. V. Filippenko

Abstract: The nature and role of the binary companion of carbon-oxygen white dwarf stars that explode as Type Ia supernovae (SNe Ia) are not yet fully understood. Past detections of circumstellar material (CSM) that contain hydrogen for a small number of SN Ia progenitor systems suggest that at least some have a nondegenerate companion. In order to constrain the prevalence, location, and quantity of CSM in… ▽ More The nature and role of the binary companion of carbon-oxygen white dwarf stars that explode as Type Ia supernovae (SNe Ia) are not yet fully understood. Past detections of circumstellar material (CSM) that contain hydrogen for a small number of SN Ia progenitor systems suggest that at least some have a nondegenerate companion. In order to constrain the prevalence, location, and quantity of CSM in SN Ia systems, we performed a near-ultraviolet (NUV) survey with the Hubble Space Telescope (HST) to look for the high-energy signature of SN Ia ejecta interacting with CSM. Our survey revealed that SN 2015cp, a SN 1991T-like overluminous SN Ia, was experiencing late-onset interaction between its ejecta and surrounding CSM at $664$ days after its light-curve peak. We present ground- and space-based follow-up observations of SN 2015cp that reveal optical emission lines of H and Ca, typical signatures of ejecta-CSM interaction. We show how SN 2015cp was likely similar to the well-studied SN Ia-CSM event PTF11kx, making it the second case in which an unambiguously classified SN Ia was observed to interact with a distant shell of CSM that contains hydrogen ($R_{\rm CSM} \gtrsim 10^{16}\ {\rm cm}$). The remainder of our HST NUV images of SNe Ia were nondetections that we use to constrain the occurrence rate of observable late-onset CSM interaction. We apply theoretical models for the emission from ejecta-CSM interaction to our NUV nondetections, and place upper limits on the mass and radial extent of CSM in SN Ia progenitor systems. △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: 29 pages, 10 figures, 6 tables; accepted to ApJ

Don't Blink: Constraining the Circumstellar Environment of the Interacting Type Ia Supernova 2015cp

Authors: C. E. Harris, P. E. Nugent, A. Horesh, J. S. Bright, R. P. Fender, M. L. Graham, K. Maguire, M. Smith, N. Butler, S. Valenti, A. V. Filippenko, O. Fox, A. Goobar, P. L. Kelly, K. J. Shen

Abstract: Despite their cosmological utility, the progenitors of Type Ia supernovae (SNe Ia) are still unknown, with many efforts focused on whether accretion from a nondegenerate companion can grow a carbon-oxygen white dwarf to near the Chandrasekhar mass. The association of SNe Ia resembling SN 1991T ("91T-like") with circumstellar interaction may be evidence for this "single-degenerate" channel. However… ▽ More Despite their cosmological utility, the progenitors of Type Ia supernovae (SNe Ia) are still unknown, with many efforts focused on whether accretion from a nondegenerate companion can grow a carbon-oxygen white dwarf to near the Chandrasekhar mass. The association of SNe Ia resembling SN 1991T ("91T-like") with circumstellar interaction may be evidence for this "single-degenerate" channel. However, the observed circumstellar medium (CSM) in these interacting systems is unlike a stellar wind -- of particular interest, it is sometimes detached from the stellar surface, residing at $\sim 10^{16}~{\rm cm}$. A Hubble Space Telescope (HST) program to discover detached CSM around 91T-like SNe Ia successfully discovered interaction nearly two years after explosion in SN 2015cp (Graham et al., 2018). In this work, we present radio and X-ray follow-up observations of SN 2015cp and analyze them in the framework of Harris, Nugent, & Kasen (2016) to limit the properties of a constant-density CSM shell in this system. Assuming the HST detection was shortly after the shock crossed the CSM, we constrain the total CSM mass in this system to be $< 0.5~{\rm M_\odot}$. This limit is comparable to the CSM mass of supernova PTF11kx, but does not rule out lower masses predicted for recurrent novae. From lessons learned modeling PTF11kx and SN 2015cp, we suggest a strategy for future observations of these events to increase the sample of known interacting SNe Ia. △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: 13 pages, 5 figures, published in The Astrophysical Journal

Journal ref: Harris et al., 2018, The Astrophysical Journal, 868, 21H

A Search For Pulsations in the Optical Light Curve of the Nova ASASSN-17hx

Authors: Lee J. Rosenthal, Ken J. Shen, Gregg Hallinan, Navtej Singh, Laura Chomiuk, Raffaella Margutti, Brian D. Metzger

Abstract: We present high-speed optical observations of the nova ASASSN-17hx, taken both immediately after its discovery and close to its first peak in brightness, to search for seconds -- minutes pulsations associated with the convective eddy turnover timescale within the nova envelope. We do not detect any periodic signal with greater than $3σ$ significance. Through injection and recovery, we rule out per… ▽ More We present high-speed optical observations of the nova ASASSN-17hx, taken both immediately after its discovery and close to its first peak in brightness, to search for seconds -- minutes pulsations associated with the convective eddy turnover timescale within the nova envelope. We do not detect any periodic signal with greater than $3σ$ significance. Through injection and recovery, we rule out periodic signals of fractional amplitude $>7.08\times10^{-4}$ on timescales of 2 seconds and fractional amplitude $>1.06\times10^{-3}$ on timescales of 10 minutes. Additional observations of novae are planned to further constrain ongoing simulations of the launch and propagation of nova winds. △ Less

Submitted 18 October, 2018; v1 submitted 16 October, 2018; originally announced October 2018.

Quantifying How Density Gradients and Front Curvature Affect Carbon Detonation Strength During Type Ia Supernovae

Authors: Broxton J. Miles, Dean M. Townsley, Ken J. Shen, F. X. Timmes, Kevin Moore

Abstract: Accurately reproducing the physics behind the detonations of Type Ia supernovae and the resultant nucleosynthetic yields is important for interpreting observations of spectra and remnants. The scales of the processes involved span orders of magnitudes, making the problem computationally impossible to ever fully resolve in full star simulations in the present and near future. In the lower density r… ▽ More Accurately reproducing the physics behind the detonations of Type Ia supernovae and the resultant nucleosynthetic yields is important for interpreting observations of spectra and remnants. The scales of the processes involved span orders of magnitudes, making the problem computationally impossible to ever fully resolve in full star simulations in the present and near future. In the lower density regions of the star, the curvature of the detonation front will slow the detonation, affecting the production of intermediate mass elements. We find that shock strengthening due to the density gradient present in the outer layers of the progenitor is essential for understanding the nucleosynthesis there, with burning extending well below the density at which a steady-state detonation is extinct. We show that a complete reaction network is not sufficient to obtain physical detonations at high densities and modest resolution due to numerical mixing at the unresolved reaction front. At low densities, below 6$\times$10$^{5}$ g cm$^{-3}$, it is possible to achieve high enough resolution to separate the shock and the reaction region,and the abundance structure predicted by fully resolved quasi-steady-state calculations is obtained. For our best current benchmark yields, we utilize a method in which the unresolved portion of Lagrangian histories are reconstructed based on fully resolved quasi-steady-state detonation calculations. These computations demonstrate that under-resolved simulations agree approximately, $\sim$10\% in post-shock values of temperature, pressure, density, and abundances, with expected detonation structures sufficiently far from the under-resolved region, but that there is still room for some improvement in the treatment of subgrid reactions in the hydrodynamics to before better than 1$\%$ can be achieved at all densities. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: Submitted to ApJ

Three Hypervelocity White Dwarfs in Gaia DR2: Evidence for Dynamically Driven Double-Degenerate Double-Detonation Type Ia Supernovae

Authors: Ken J. Shen, Douglas Boubert, Boris T. Gänsicke, Saurabh W. Jha, Jennifer E. Andrews, Laura Chomiuk, Ryan J. Foley, Morgan Fraser, Mariusz Gromadzki, James Guillochon, Marissa M. Kotze, Kate Maguire, Matthew R. Siebert, Nathan Smith, Jay Strader, Carles Badenes, Wolfgang E. Kerzendorf, Detlev Koester, Markus Kromer, Broxton Miles, Rüdiger Pakmor, Josiah Schwab, Odette Toloza, Silvia Toonen, Dean M. Townsley , et al. (1 additional authors not shown)

Abstract: Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 k… ▽ More Double detonations in double white dwarf (WD) binaries undergoing unstable mass transfer have emerged in recent years as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. One potential outcome of this "dynamically driven double-degenerate double-detonation" (D^6) scenario is that the companion WD survives the explosion and is flung away with a velocity equal to its > 1000 km/s pre-SN orbital velocity. We perform a search for these hypervelocity runaway WDs using Gaia's second data release. In this paper, we discuss seven candidates followed up with ground-based instruments. Three sources are likely to be some of the fastest known stars in the Milky Way, with total Galactocentric velocities between 1000 and 3000 km/s, and are consistent with having previously been companion WDs in pre-SN Ia systems. However, although the radial velocity of one of the stars is > 1000 km/s, the radial velocities of the other two stars are puzzlingly consistent with 0. The combined five-parameter astrometric solutions from Gaia and radial velocities from follow-up spectra yield tentative 6D confirmation of the D^6 scenario. The past position of one of these stars places it within a faint, old SN remnant, further strengthening the interpretation of these candidates as hypervelocity runaways from binary systems that underwent SNe Ia. △ Less

Submitted 18 July, 2018; v1 submitted 30 April, 2018; originally announced April 2018.

Comments: Accepted for publication in ApJ. Minor corrections for clarity. D6 spectra are available as ancillary data files

Energetic eruptions leading to a peculiar hydrogen-rich explosion of a massive star

Authors: Iair Arcavi, D. Andrew Howell, Daniel Kasen, Lars Bildsten, Griffin Hosseinzadeh, Curtis McCully, Zheng Chuen Wong, Sarah Rebekah Katz, Avishay Gal-Yam, Jesper Sollerman, Francesco Taddia, Giorgos Leloudas, Christoffer Fremling, Peter E. Nugent, Assaf Horesh, Kunal Mooley, Clare Rumsey, S. Bradley Cenko, Melissa L. Graham, Daniel A. Perley, Ehud Nakar, Nir J. Shaviv, Omer Bromberg, Ken J. Shen, Eran O. Ofek , et al. (28 additional authors not shown)

Abstract: Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curv… ▽ More Every supernova hitherto observed has been considered to be the terminal explosion of a star. Moreover, all supernovae with absorption lines in their spectra show those lines decreasing in velocity over time, as the ejecta expand and thin, revealing slower moving material that was previously hidden. In addition, every supernova that exhibits the absorption lines of hydrogen has one main light-curve peak, or a plateau in luminosity, lasting approximately 100 days before declining. Here we report observations of iPTF14hls, an event that has spectra identical to a hydrogen-rich core-collapse supernova, but characteristics that differ extensively from those of known supernovae. The light curve has at least five peaks and remains bright for more than 600 days; the absorption lines show little to no decrease in velocity; and the radius of the line-forming region is more than an order of magnitude bigger than the radius of the photosphere derived from the continuum emission. These characteristics are consistent with a shell of several tens of solar masses ejected by the star at supernova-level energies a few hundred days before a terminal explosion. Another possible eruption was recorded at the same position in 1954. Multiple energetic pre-supernova eruptions are expected to occur in stars of 95-130 solar masses, which experience the pulsational pair instability. That model, however, does not account for the continued presence of hydrogen, or the energetics observed here. Another mechanism for the violent ejection of mass in massive stars may be required. △ Less

Submitted 7 November, 2017; originally announced November 2017.

Comments: Published in Nature

The Evolution of the Type Ia Supernova Luminosity Function

Authors: Ken J. Shen, Silvia Toonen, Or Graur

Abstract: Type Ia supernovae (SNe Ia) exhibit a wide diversity of peak luminosities and light curve shapes: the faintest SNe Ia are 10 times less luminous and evolve more rapidly than the brightest SNe Ia. Their differing characteristics also extend to their stellar age distributions, with fainter SNe Ia preferentially occurring in old stellar populations and vice versa. In this Letter, we quantify this SN… ▽ More Type Ia supernovae (SNe Ia) exhibit a wide diversity of peak luminosities and light curve shapes: the faintest SNe Ia are 10 times less luminous and evolve more rapidly than the brightest SNe Ia. Their differing characteristics also extend to their stellar age distributions, with fainter SNe Ia preferentially occurring in old stellar populations and vice versa. In this Letter, we quantify this SN Ia luminosity - stellar age connection using data from the Lick Observatory Supernova Search (LOSS). Our binary population synthesis calculations agree qualitatively with the observed trend in the >1 Gyr-old populations probed by LOSS if the majority of SNe Ia arise from prompt detonations of sub-Chandrasekhar mass white dwarfs (WDs) in double WD systems. Under appropriate assumptions, we show that double WD systems with less massive primaries, which yield fainter SNe Ia, interact and explode at older ages than those with more massive primaries. We find that prompt detonations in double WD systems are capable of reproducing the observed evolution of the SN Ia luminosity function, a constraint that any SN Ia progenitor scenario must confront. △ Less

Submitted 6 December, 2017; v1 submitted 25 October, 2017; originally announced October 2017.

Comments: Accepted for publication in ApJL. Minor changes to previous version for clarity. Data used to construct the observational CDFs in Figure 4 are available in an ancillary file

A Search for a Surviving White Dwarf Companion in SN 1006

Authors: W. E. Kerzendorf, G. Strampelli, K. J. Shen, J. Schwab, R. Pakmor, T. Do, J. Buchner, A. Rest

Abstract: Multiple channels have been proposed to produce Type Ia supernovae, with many scenarios suggesting that the exploding white dwarf accretes from a binary companion pre-explosion. In almost all cases, theory suggests that this companion will survive. However, no such companion has been unambiguously identified in ancient supernova remnants -- possibly falsifying the accretion scenario. Existing surv… ▽ More Multiple channels have been proposed to produce Type Ia supernovae, with many scenarios suggesting that the exploding white dwarf accretes from a binary companion pre-explosion. In almost all cases, theory suggests that this companion will survive. However, no such companion has been unambiguously identified in ancient supernova remnants -- possibly falsifying the accretion scenario. Existing surveys, however, have only looked for stars as faint as $\approx 0.1 L_\odot$ and thus would have missed a surviving white dwarf companion. In this work, we present very deep DECAM imaging $(u, g, r, z)$ of the Type Ia supernova remnant SN 1006 specifically to search for a potential surviving white dwarf companion. We find no object within the inner third of the SN 1006 remnant that is consistent with a relatively young cooling white dwarf. We find that if there is a companion white dwarf, it must have formed long ago and cooled undisturbed for $> 10^8$ yr to be consistent with the redder objects in our sample. We conclude that our findings are consistent with the complete destruction of the secondary (such as in a merger) or an unexpectedly cool and thus very dim surviving companion white dwarf. △ Less

Submitted 19 September, 2017; originally announced September 2017.

Comments: 9 pages, 5 figure, submitted to MNRAS - comments welcome

Sub-Chandrasekhar-mass white dwarf detonations revisited

Authors: Ken J. Shen, Daniel Kasen, Broxton J. Miles, Dean M. Townsley

Abstract: The detonation of a sub-Chandrasekhar-mass white dwarf (WD) has emerged as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. Recent studies have suggested that the rapid transfer of a very small amount of helium from one WD to another is sufficient to ignite a helium shell detonation that subsequently triggers a carbon core detonation, yielding a "dynamically-driven double… ▽ More The detonation of a sub-Chandrasekhar-mass white dwarf (WD) has emerged as one of the most promising Type Ia supernova (SN Ia) progenitor scenarios. Recent studies have suggested that the rapid transfer of a very small amount of helium from one WD to another is sufficient to ignite a helium shell detonation that subsequently triggers a carbon core detonation, yielding a "dynamically-driven double degenerate double detonation" SN Ia. Because the helium shell that surrounds the core explosion is so minimal, this scenario approaches the limiting case of a bare C/O WD detonation. Motivated by discrepancies in previous literature and by a recent need for detailed nucleosynthetic data, we revisit simulations of naked C/O WD detonations in this paper. We disagree to some extent with the nucleosynthetic results of previous work on sub-Chandrasekhar-mass bare C/O WD detonations; e.g., we find that a median-brightness SN Ia is produced by the detonation of a 1.0 Msol WD instead of a more massive and rarer 1.1 Msol WD. The neutron-rich nucleosynthesis in our simulations agrees broadly with some observational constraints, although tensions remain with others. There are also discrepancies related to the velocities of the outer ejecta and light curve shapes, but overall our synthetic light curves and spectra are roughly consistent with observations. We are hopeful that future multi-dimensional simulations will resolve these issues and further bolster the dynamically-driven double degenerate double detonation scenario's potential to explain most SNe Ia. △ Less

Submitted 30 January, 2018; v1 submitted 6 June, 2017; originally announced June 2017.

Comments: Accepted for publication in ApJ. Bulk yields are available in an ancillary data file