Showing 1–50 of 69 results for author: Townsley, D M

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

Sensitivity of 3D Convective Urca Simulations to Changes in Urca Reactions

Authors: Brendan Boyd, Alexander Smith Clark, Alan C. Calder, Dean M. Townsley, Michael Zingale

Abstract: A proposed setting for thermonuclear (Type Ia) supernovae is a white dwarf that has gained mass from a companion to the point of carbon ignition in the core. There is a simmering phase in the early stages of burning that involves the formation and growth of a core convection zone. One aspect of this phase is the convective Urca process, a linking of weak nuclear reactions to convection that may al… ▽ More A proposed setting for thermonuclear (Type Ia) supernovae is a white dwarf that has gained mass from a companion to the point of carbon ignition in the core. There is a simmering phase in the early stages of burning that involves the formation and growth of a core convection zone. One aspect of this phase is the convective Urca process, a linking of weak nuclear reactions to convection that may alter the composition and structure of the white dwarf. Convective Urca is not well understood and requires 3D fluid simulations to realistically model. Additionally, the convection is relatively slow (Mach number less than 0.005) so a low-Mach method is needed to make simulating computationally feasible. Using the MAESTROeX low-Mach hydrodynamics code, we investigate recent changes to how the weak reactions are modeled in the convective Urca simulations. We present results that quantify the changes to the reaction rates and their impact on the evolution of the simulation. △ Less

Submitted 25 April, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 9 pages, 7 figures, submitted to the Proceedings of 15th International Conference on Numerical Modeling of Space Plasma Flows (AstroNum)

Journal ref: J. Phys.: Conf. Ser. 2742 (2024)

Challenges Modeling the Low-Luminosity Type Iax Supernovae

Authors: Catherine Feldman, Ellis Eisenberg, Dean M. Townsley, Alan C. Calder

Abstract: Numerical models allow the investigation of phenomena that cannot exist in a laboratory. Computational simulations are therefore essential for advancing our knowledge of astrophysics, however, the very nature of simulation requires making assumptions that can substantially affect their outcome. Here, we present the challenges faced when simulating dim thermonuclear explosions, Type Iax supernovae.… ▽ More Numerical models allow the investigation of phenomena that cannot exist in a laboratory. Computational simulations are therefore essential for advancing our knowledge of astrophysics, however, the very nature of simulation requires making assumptions that can substantially affect their outcome. Here, we present the challenges faced when simulating dim thermonuclear explosions, Type Iax supernovae. This class of dim events produce a slow moving, sparse ejecta that presents challenges for simulation. We investigate the limitations of the equation of state and its applicability to the expanding, cooling ejecta. We also discuss how the "fluff", i.e. the low-density gas on the grid in lieu of vacuum, inhibits the ejecta as it expands. We explore how the final state of the simulation changes as we vary the character of the burning, which influences the outcome of the explosion. These challenges are applicable to a wide range of astrophysical simulations, and are important to discuss and overcome as a community. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: 10 pages, 7 figures, submitted to the Proceedings of 15th International Conference on Numerical Modeling of Space Plasma Flows (AstroNum)

Dimming the Lights: 2D Simulations of Deflagrations of Hybrid C/O/Ne White Dwarfs using FLASH

Authors: Catherine Feldman, Nathanael Gutierrez, Ellis Eisenberg, Donald E. Willcox, Dean M. Townsley, Alan C. Calder

Abstract: The dimmest and most numerous outlier of the Type Ia supernova population, Type Iax events, is increasingly being found in the results of observational campaigns. There is currently no single accepted model to describe these events. This 2D study explores the viability of modeling Type Iax events as a hybrid C/O/Ne white dwarf progenitor undergoing a deflagration using the multi-physics software F… ▽ More The dimmest and most numerous outlier of the Type Ia supernova population, Type Iax events, is increasingly being found in the results of observational campaigns. There is currently no single accepted model to describe these events. This 2D study explores the viability of modeling Type Iax events as a hybrid C/O/Ne white dwarf progenitor undergoing a deflagration using the multi-physics software FLASH. This hybrid was created using the stellar evolution code MESA, and its C-depleted core and mixed structure have demonstrated lower yields than traditional C/O progenitors in previous deflagration-to-detonation studies. To generate a sample, 30 "realizations" of this simulation were performed, the only difference being the shape of the initial matchhead used to start the deflagration. As consistent with earlier work, these realizations produce the familiar hot dense bound remnant surrounded by sparse ejecta. Our results indicate the majority of the star remains unburned (~70%) and bound (>90%). Our realizations produce total ejecta yields on the order of 10$^{-2}$ - 10$^{-1}$ solar masses, ejected $^{56}$Ni yields on the order of 10$^{-4}$ - 10$^{-2}$ solar masses, and ejecta kinetic energies on the order of 10$^{48}$ - 10$^{49}$ ergs. Compared to yields inferred from recent observations of the dimmest Type Iax events - SN 2007qd, SN 2008ha, SN 2010ae, SN 2019gsc, SN 2019muj, SN 2020kyg, and SN 2021fcg - our simulation produces comparable $^{56}$Ni yields, but too-small total yields and kinetic energies. Reignition of the remnant is also seen in some realizations. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: 43 pages, 15 figures, 4 tables. To be published in ApJ

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)

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

Gravity modes on rapidly rotating accreting white dwarfs and their variation after dwarf novae

Authors: Praphull Kumar, Dean M. Townsley

Abstract: Accreting white dwarfs in Cataclysmic variables (CVs) show short-period (tens of minutes) brightness variations that are consistent with non-radial oscillations similar to gravity (g) modes observed in isolated white dwarfs (WDs). GW Librae, a dwarf nova, was the first CV in which non-radial oscillations were observed and continues to be the best studied accreting WD displaying these pulsations. U… ▽ More Accreting white dwarfs in Cataclysmic variables (CVs) show short-period (tens of minutes) brightness variations that are consistent with non-radial oscillations similar to gravity (g) modes observed in isolated white dwarfs (WDs). GW Librae, a dwarf nova, was the first CV in which non-radial oscillations were observed and continues to be the best studied accreting WD displaying these pulsations. Unlike isolated WDs, accreting WDs rotate rapidly, with spin periods comparable to or shorter than typical low-order oscillation periods. Accreting WDs also have a different relationship between their interior temperature and surface temperature. The surface temperature of an accreting WD varies on a months to year timescale between dwarf novae accretion events, allowing study of how this temperature change effects g-mode behavior. Here we show results from adiabatic seismological calculations for accreting WDs, focusing on low-order ($\ell=1$) modes. We demonstrate how g-modes vary in response to temperature changes in the subsurface layers due to a dwarf nova accretion event. These calculations include rotation non-perturbatively, required by the high spin rate. We discuss the thermal history of these accreting WDs, and compare the seismological properties with and without rotation. Comparison of $g$-mode frequencies to observed objects may allow inference of features of the structure of the WD such as mass, surface abundance, accretion history, and more. The variation of mode frequencies during cooling after an outburst provides a novel method of identifying modes. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: 18 pages, 13 figures, Accepted to ApJ

Journal ref: ApJ 951(2023) 122

Flash-X, a multiphysics simulation software instrument

Authors: Anshu Dubey, Klaus Weide, Jared O'Neal, Akash Dhruv, Sean Couch, J. Austin Harris, Tom Klosterman, Rajeev Jain, Johann Rudi, Bronson Messer, Michael Pajkos, Jared Carlson, Ran Chu, Mohamed Wahib, Saurabh Chawdhary, Paul M. Ricker, Dongwook Lee, Katie Antypas, Katherine M. Riley, Christopher Daley, Murali Ganapathy, Francis X. Timmes, Dean M. Townsley, Marcos Vanella, John Bachan , et al. (6 additional authors not shown)

Abstract: Flash-X is a highly composable multiphysics software system that can be used to simulate physical phenomena in several scientific domains. It derives some of its solvers from FLASH, which was first released in 2000. Flash-X has a new framework that relies on abstractions and asynchronous communications for performance portability across a range of increasingly heterogeneous hardware platforms. Fla… ▽ More Flash-X is a highly composable multiphysics software system that can be used to simulate physical phenomena in several scientific domains. It derives some of its solvers from FLASH, which was first released in 2000. Flash-X has a new framework that relies on abstractions and asynchronous communications for performance portability across a range of increasingly heterogeneous hardware platforms. Flash-X is meant primarily for solving Eulerian formulations of applications with compressible and/or incompressible reactive flows. It also has a built-in, versatile Lagrangian framework that can be used in many different ways, including implementing tracers, particle-in-cell simulations, and immersed boundary methods. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 16 pages, 5 Figures, published open access in SoftwareX

Journal ref: SoftwareX, Volume 19, 2022, 101168,ISSN 2352-7110

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

The Heating and Pulsations of V386 Serpentis after its 2019 Dwarf Nova Outburst

Authors: Paula Szkody, Patrick Godon, Boris T. Gaensicke, Stella Kafka, Odette F. T. Castillo, Keaton J. Bell, P. B. Cho, Edward M. Sion, Praphull Kumar, Dean M. Townsley, Zach Vanderbosch, Karen I. Winget, Claire J. Olde Loohuis

Abstract: Following the pulsation spectrum of a white dwarf through the heating and cooling involved in a dwarf nova outburst cycle provides a unique view of the changes to convective driving that take place on timescales of months versus millenia for non-accreting white dwarfs. In 2019 January the dwarf nova V386 Ser (one of a small number containing an accreting, pulsating white dwarf), underwent a large… ▽ More Following the pulsation spectrum of a white dwarf through the heating and cooling involved in a dwarf nova outburst cycle provides a unique view of the changes to convective driving that take place on timescales of months versus millenia for non-accreting white dwarfs. In 2019 January the dwarf nova V386 Ser (one of a small number containing an accreting, pulsating white dwarf), underwent a large amplitude outburst. Hubble Space Telescope ultraviolet spectra were obtained 7 and 13 months after outburst along with optical ground-based photometry during this interval and high-speed photometry at 5.5 and 17 months after outburst. The resulting spectral and pulsational analysis shows a cooling of the white dwarf from 21,020 K to 18,750 K (with a gravity log(g) = 8.1) between the two UV observations, along with the presence of strong pulsations evident in both UV and optical at a much shorter period after outburst than at quiescence. The pulsation periods consistently lengthened during the year following outburst, in agreement with pulsation theory. However, it remains to be seen if the behavior at longer times past outburst will mimic the unusual non-monotonic cooling and long periods evident in the similar system GW Lib. △ Less

Submitted 29 April, 2021; originally announced April 2021.

Comments: 25 pages, 14 figures, 5 Tables

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

Type Ia Supernova Explosions from Hybrid Carbon-Oxygen-Neon White Dwarf Progenitors That Have Mixed During Cooling

Authors: Carlyn N. Augustine, Donald E. Willcox, Jared Brooks, Dean M. Townsley, Alan C. Calder

Abstract: The creation of "hybrid" white dwarfs, made of a C-O core within a O-Ne shell has been proposed, and studies indicate that ignition in the C-rich central region makes these viable progenitors for thermonuclear (type Ia) supernovae. Recent work found that the C-O core is mixed with the surrounding O-Ne as the white dwarf cools prior to accretion, which results in lower central C fractions in the ma… ▽ More The creation of "hybrid" white dwarfs, made of a C-O core within a O-Ne shell has been proposed, and studies indicate that ignition in the C-rich central region makes these viable progenitors for thermonuclear (type Ia) supernovae. Recent work found that the C-O core is mixed with the surrounding O-Ne as the white dwarf cools prior to accretion, which results in lower central C fractions in the massive progenitor than previously assumed. To further investigate the efficacy of hybrid white dwarfs as progenitors of thermonuclear supernovae, we performed simulations of thermonuclear supernovae from a new series of hybrid progenitors that include the effects of mixing during cooling. The progenitor white dwarf model was constructed with the one-dimensional stellar evolution code MESA and represented a star evolved through the phase of unstable interior mixing followed by accretion until it reached conditions for the ignition of carbon burning. This MESA model was then mapped to a two-dimensional initial condition for explosions simulated with FLASH. For comparison, similar simulations were performed for a traditional C-O progenitor white dwarf. By comparing the yields of the explosions, we find that, as with earlier studies, the lower C abundance in the hybrid progenitor compared to the traditional C-O progenitor leads to a lower average yield of 56Ni. Although the unmixed hybrid WD showed a similar decrement also in total iron group yield, the mixed case does not and produces a smaller fraction of iron group elements in the form of 56Ni. We attribute this to the higher central density required for ignition and the location, center or off-center, of deflagration ignition. △ Less

Submitted 27 October, 2019; originally announced October 2019.

Comments: 19 pages, 7 figures, accepted by the Astrophysical Journal

Modeling subgrid combustion processes in simulations of thermonuclear supernovae

Authors: Dean M. Townsley, Alan C. Calder, Broxton J. Miles

Abstract: Supernovae of type Ia are thought to arise from the thermonuclear incineration of a carbon-oxygen white dwarf stellar remnant. However, the detailed explosion scenario and stellar evolutionary origin scenario -- or scenarios -- which lead to observed supernovae are still quite uncertain. One of the principal tests of proposed scenarios is comparison with the explosion products inferred, for exampl… ▽ More Supernovae of type Ia are thought to arise from the thermonuclear incineration of a carbon-oxygen white dwarf stellar remnant. However, the detailed explosion scenario and stellar evolutionary origin scenario -- or scenarios -- which lead to observed supernovae are still quite uncertain. One of the principal tests of proposed scenarios is comparison with the explosion products inferred, for example, from the spectrum of the supernovae. Making this comparison requires computation of the combustion dynamics and products through simulation of proposed scenarios. Here we discuss two specific proposed explosion scenarios, the deflagration-detonation transition and the helium shell double detonation, With these two examples in mind, we proceed to discuss challenges to computational modeling of the combustion taking place in these explosions. Both subsonically and supersonically propagating reaction fronts are discussed, called deflagrations and detonations respectively. Several major stages of the combustion occur on length and time scales that are many orders of magnitude smaller than those accessible in simulations of the explosion. Models which attempt to capture this sub-grid behavior and the verification of those models is briefly discussed. △ Less

Submitted 16 August, 2019; originally announced August 2019.

Comments: 7 pages, 3 figures, proceedings from the 13th International Conference on Numerical Modeling of Space Plasma Flows (ASTRONUM-2018)

Journal ref: 2019 J. Phys.: Conf. Ser. 1225, 012004

Thermonuclear (Type Ia) Supernovae and Progenitor Evolution

Authors: Alan C. Calder, Don E. Willcox, Christopher J. DeGrendele, Desmond Shangase, Michael Zingale, Dean M. Townsley

Abstract: Thermonuclear (type Ia) supernovae are bright stellar explosions with the unique property that the light curves can be standardized, allowing them to be used as distance indicators for cosmological studies. Many fundamental questions bout these events remain, however. We provide a critique of our present understanding of these and present results of simulations assuming the single-degenerate proge… ▽ More Thermonuclear (type Ia) supernovae are bright stellar explosions with the unique property that the light curves can be standardized, allowing them to be used as distance indicators for cosmological studies. Many fundamental questions bout these events remain, however. We provide a critique of our present understanding of these and present results of simulations assuming the single-degenerate progenitor model consisting of a white dwarf that has gained mass from a stellar companion. We present results from full three-dimensional simulations of convection with weak reactions comprising the A=23 Urca process in the progenitor white dwarf. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: 10 pages, 4 figures, accepted to the proceedings of ASTRONUM 2018, the 13th International Conference on Numerical Modeling of Space Plasma Flows in Panama City Beach, Florida, USA, on 25-29 June, 2018

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

Evidence for mass accretion driven by spiral shocks onto the white dwarf in SDSS J123813.73-033933.0

Authors: A. F. Pala, B. T. Gänsicke, T. R. Marsh, E. Breedt, J. J. Hermes, J. D. Landstreet, M. R. Schreiber, D. M. Townsley, L. Wang, A. Aungwerojwit, F. -J. Hambsch, B. Monard, G. Myers, P. Nelson, R. Pickard, G. Poyner, D. E. Reichart, R. Stubbings, P. Godon, P. Szkody, D. De Martino, V. S. Dhillon, C. Knigge, S. G. Parsons

Abstract: We present high-time-resolution photometry and phase-resolved spectroscopy of the short-period ($P_\mathrm{orb} = 80.52\,\mathrm{min}$) cataclysmic variable SDSS J123813.73-033933.0, observed with the $\mathit{Hubble}$ $\mathit{Space}$ $\mathit{Telescope}$ $\mathit{(HST)}$, the $\mathit{Kepler/K2}$ mission and the Very Large Telescope (VLT). We also report observations of the first detected super-… ▽ More We present high-time-resolution photometry and phase-resolved spectroscopy of the short-period ($P_\mathrm{orb} = 80.52\,\mathrm{min}$) cataclysmic variable SDSS J123813.73-033933.0, observed with the $\mathit{Hubble}$ $\mathit{Space}$ $\mathit{Telescope}$ $\mathit{(HST)}$, the $\mathit{Kepler/K2}$ mission and the Very Large Telescope (VLT). We also report observations of the first detected super-outburst. SDSS J1238-0339 shows two types of variability: quasi-regular brightenings recurring every $\simeq 8.5$ h during which the system increases in brightness by $\simeq 0.5$mag, and a double hump quasi-sinusoidal modulation at the orbital period. The detailed $\mathit{K2}$ light curve reveals that the amplitude of the double-humps increases during the brightenings and that their phase undergoes a $\simeq 90^{\circ}$ phase shift with respect to the quiescent intervals. The $\mathit{HST}$ data unambiguously demonstrate that these phenomena both arise from the heating and cooling of two relatively large regions on the white dwarf. We suggest that the double-hump modulation is related to spiral shocks in the accretion disc resulting in an enhanced accretion rate heating two localised regions on the white dwarf, with the structure of the shocks fixed in the binary frame explaining the period of the double humps. The physical origin of the 8.5 h brightenings is less clear. However, the correlation between the observed variations of the amplitude and phase of the double-humps with the occurrence of the brightenings is supportive of an origin in thermal instabilities in the accretion disc. △ Less

Submitted 12 December, 2018; v1 submitted 14 November, 2018; originally announced November 2018.

Comments: 26 pages, 22 figures. MNRAS, in press

The Intrinsic Stochasticity of the $^{56}$Ni Distribution of Single-Degenerate Type Ia Supernovae

Authors: Chris Byrohl, Robert T. Fisher, Dean M. Townsley

Abstract: Binary Chandrasekhar-mass white dwarfs accreting mass from non-degenerate stellar companions through the single-degenerate channel have reigned for decades as the leading explanation of Type Ia supernovae. Yet, a comprehensive theoretical explanation has not yet emerged to explain the expected properties of the canonical near-Chandrasekhar-mass white dwarf model. A simmering phase within the conve… ▽ More Binary Chandrasekhar-mass white dwarfs accreting mass from non-degenerate stellar companions through the single-degenerate channel have reigned for decades as the leading explanation of Type Ia supernovae. Yet, a comprehensive theoretical explanation has not yet emerged to explain the expected properties of the canonical near-Chandrasekhar-mass white dwarf model. A simmering phase within the convective core of the white dwarf leads to the ignition of one or more flame bubbles scattered across the core. Consequently, near-Chandrasekhar-mass single-degenerate SNe Ia are inherently stochastic, and are expected to lead to a range of outcomes, from subluminous SN 2002cx-like events, to overluminous SN 1991T-like events. However, all prior simulations of the single-degenerate channel carried through the detonation phase have set the ignition points as free parameters. In this work, for the first time, we place ignition points as predicted by {\it ab initio} models of the convective phase leading up to ignition, and follow through the detonation phase in fully three-dimensional simulations. Single-degenerates in this framework are characteristically overluminous. Using a statistical approach, we determine the $^{56}$Ni mass distribution arising from stochastic ignition. While there is a total spread of $\gtrsim 0.2 M_{\odot}$ for detonating models, the distribution is strongly left-skewed, and with a narrow standard deviation of $\simeq 0.03 M_{\odot}$. Conversely, if single-degenerates are not overluminous but primarily yield normal or failed events, then the models require fine-tuning of the ignition parameters, or otherwise require revised physics or progenitor models. We discuss implications of our findings for the modeling of single-degenerate SNe Ia. △ Less

Submitted 18 October, 2018; originally announced October 2018.

Comments: 14 pages, 10 figures. Submitted to ApJ. Comments welcome

Astrophysical Validation

Authors: A. C. Calder, D. M. Townsley

Abstract: We present examples of validating components of an astrophysical simulation code. Problems of stellar astrophysics are multi-dimensional and involve physics acting on large ranges of length and time scales that are impossible to include in macroscopic models on present computational resources. Simulating these events thus necessitates the development of sub-grid-scale models and the capability to… ▽ More We present examples of validating components of an astrophysical simulation code. Problems of stellar astrophysics are multi-dimensional and involve physics acting on large ranges of length and time scales that are impossible to include in macroscopic models on present computational resources. Simulating these events thus necessitates the development of sub-grid-scale models and the capability to post-process simulations with higher-fidelitymethods.We present an overview of the problem of validating astrophysical models and simulations illustrated with two examples. First, we present a study aimed at validating hydrodynamicswith high energy density laboratory experiments probing shocks and fluid instabilities. Second, we present an effort at validating code modules for use in both macroscopic simulations of astrophysical events and for post processing Lagrangian tracer particles to calculate detailed abundances from thermonuclear reactions occurring during an event. △ Less

Submitted 6 September, 2018; originally announced September 2018.

Comments: 29 pages, 6 figures, To appear in Computer Simulation Validation. Fundamental Concepts, Methodological Frameworks, Philosophical Perspectives. Claus Beisbart and Nicole J. Saam. Eds. (forthcoming). Cham: Springer International Publishing

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

Anatomy of the hyper-runaway star LP 40-365 with Gaia

Authors: R. Raddi, M. A. Hollands, B. T. Gaensicke, D. M. Townsley, J. J. Hermes, N. P. Gentile Fusillo, D. Koester

Abstract: LP 40-365 (aka GD 492) is a nearby low-luminosity hyper-runaway star with an extremely unusual atmospheric composition, which has been proposed as the remnant of a white dwarf that survived a subluminous Type Ia supernova (SN Ia) in a single-degenerate scenario. Adopting the Gaia Data Release (DR2) parallax, 1.58 +/- 0.03 mas, we estimate a radius of 0.18 +/- 0.01 Rsun, confirming LP 40-365 as a s… ▽ More LP 40-365 (aka GD 492) is a nearby low-luminosity hyper-runaway star with an extremely unusual atmospheric composition, which has been proposed as the remnant of a white dwarf that survived a subluminous Type Ia supernova (SN Ia) in a single-degenerate scenario. Adopting the Gaia Data Release (DR2) parallax, 1.58 +/- 0.03 mas, we estimate a radius of 0.18 +/- 0.01 Rsun, confirming LP 40-365 as a subluminous star that is ~ 15 times larger than a typical white dwarf and is compatible with the SN Ia remnant scenario. We present an updated kinematic analysis, making use of the Gaia parallax and proper motion, and confirm that Lp 40-365 is leaving the Milky Way at about 1.5 times the escape velocity of the Solar neighbourhood with a rest-frame velocity of 852 +/- 10 km/s. Integrating the past trajectories of LP 40-365, we confirm it crossed the Galactic disc 5.0 +/- 0.3 Myr ago in the direction of Carina, likely coming from beneath the plane. Finally, we estimate that LP 40-365 was ejected from its progenitor binary with a velocity of at least 600 km/s, which is compatible with theoretical predictions for close binaries containing a white dwarf and a helium-star donor. △ Less

Submitted 8 June, 2018; v1 submitted 25 April, 2018; originally announced April 2018.

Comments: 6 pages, 4 figures, 1 table. Accepted for publication on MNRAS

Further insight on the hypervelocity white dwarf, LP 40-365 (GD 492): a nearby emissary from a single-degenerate Type Ia supernova

Authors: R. Raddi, M. A. Hollands, D. Koester, B. T. Gaensicke, N. P. Gentile-Fusillo, J. J. Hermes, D. M. Townsley

Abstract: The recently discovered hypervelocity white dwarf LP 40-65 (aka GD 492) has been suggested as the outcome of the failed disruption of a white dwarf in a sub-luminous Type Ia supernova (SN Ia). We present new observations confirming GD 492 as a single star with unique spectral features. Our spectroscopic analysis suggests that a helium-dominated atmosphere, with ~ 33 percent neon and 2 percent oxyg… ▽ More The recently discovered hypervelocity white dwarf LP 40-65 (aka GD 492) has been suggested as the outcome of the failed disruption of a white dwarf in a sub-luminous Type Ia supernova (SN Ia). We present new observations confirming GD 492 as a single star with unique spectral features. Our spectroscopic analysis suggests that a helium-dominated atmosphere, with ~ 33 percent neon and 2 percent oxygen by mass, can reproduce most of the observed properties of this highly unusual star. Although our atmospheric model contrasts with the previous analysis in terms of dominant atmospheric species, we confirm that the atmosphere of GD 492 is strongly hydrogen deficient, log(H/He) < -5, and displays traces of eleven other alpha- and iron-group elements (with sulfur, chromium, manganese, and titanium as new detections), indicating nuclear processing of carbon and silicon. We measure a manganese-to-iron ratio seven times larger than Solar. While the observed abundances of GD 492 do not fully match any predicted nuclear yields of a partially-burned supernova remnant, the manganese excess strongly favors a link with a single-degenerate SN Ia event over alternative scenarios. △ Less

Submitted 20 March, 2018; originally announced March 2018.

Comments: 12 pages, 7 figures, 4 tables. Accepted for publication in the Astrophysical Journal

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

Nucleosynthesis in thermonuclear supernovae

Authors: Ivo R. Seitenzahl, Dean M. Townsley

Abstract: The explosion energy of thermonuclear (Type Ia) supernovae is derived from the difference in nuclear binding energy liberated in the explosive fusion of light 'fuel' nuclei, predominantly carbon and oxygen, into more tightly bound nuclear 'ash' dominated by iron and silicon group elements. The very same explosive thermonuclear fusion event is also one of the major processes contributing to the nuc… ▽ More The explosion energy of thermonuclear (Type Ia) supernovae is derived from the difference in nuclear binding energy liberated in the explosive fusion of light 'fuel' nuclei, predominantly carbon and oxygen, into more tightly bound nuclear 'ash' dominated by iron and silicon group elements. The very same explosive thermonuclear fusion event is also one of the major processes contributing to the nucleosynthesis of the heavy elements, in particular the iron-group elements. For example, most of the iron and manganese in the sun and its planetary system were produced in thermonuclear supernovae. Here, we review the physics of explosive thermonuclear burning in carbon-oxygen white dwarf material and the methodologies utilized in calculating predicted nucleosynthesis from hydrodynamic explosion models. While the dominant explosion scenario remains unclear, many aspects of the nuclear combustion and nucleosynthesis are common to all models and must occur in some form in order to produce the observed yields. We summarize the predicted nucleosynthetic yields for existing explosion models, placing particular emphasis on characteristic differences in the nucleosynthetic signatures of the different suggested scenarios leading to Type Ia supernovae. Following this, we discuss how these signatures compare with observations of several individual supernovae, remnants, and the composition of material in our galaxy and galaxy clusters. △ Less

Submitted 2 April, 2017; originally announced April 2017.

Comments: Author version of a chapter for the 'Handbook of Supernovae', edited by A. Alsabti and P. Murdin, Springer. 22 pages, 2 figures

Observational evidence for high neutronization in supernova remnants: implications for Type Ia supernova progenitors

Authors: Héctor Martínez-Rodríguez, Carles Badenes, Hiroya Yamaguchi, Eduardo Bravo, F. X. Timmes, Broxton J. Miles, Dean M. Townsley, Anthony L. Piro, Hideyuki Mori, Brett Andrews, Sangwook Park

Abstract: The physical process whereby a carbon--oxygen white dwarf explodes as a Type Ia supernova (SN Ia) remains highly uncertain. The degree of neutronization in SN Ia ejecta holds clues to this process because it depends on the mass and the metallicity of the stellar progenitor, and on the thermodynamic history prior to the explosion. We report on a new method to determine ejecta neutronization using C… ▽ More The physical process whereby a carbon--oxygen white dwarf explodes as a Type Ia supernova (SN Ia) remains highly uncertain. The degree of neutronization in SN Ia ejecta holds clues to this process because it depends on the mass and the metallicity of the stellar progenitor, and on the thermodynamic history prior to the explosion. We report on a new method to determine ejecta neutronization using Ca and S lines in the X-ray spectra of Type Ia supernova remnants (SNRs). Applying this method to \textit{Suzaku} data of Tycho, Kepler, 3C 397 and G337.2$-$0.7 in the Milky Way, and N103B in the Large Magellanic Cloud, we find that the neutronization of the ejecta in N103B is comparable to that of Tycho and Kepler, which suggests that progenitor metallicity is not the only source of neutronization in SNe Ia. We then use a grid of SN Ia explosion models to infer the metallicities of the stellar progenitors of our SNRs. The implied metallicities of 3C 397, G337.2$-$0.7, and N103B are major outliers compared to the local stellar metallicity distribution functions, indicating that progenitor metallicity can be ruled out as the origin of neutronization for these SNRs. Although the relationship between ejecta neutronization and equivalent progenitor metallicity is subject to uncertainties stemming from the $^{12}$C$\,$+$^{16}$O reaction rate, which affects the Ca/S mass ratio, our main results are not sensitive to these details. △ Less

Submitted 28 June, 2017; v1 submitted 24 January, 2017; originally announced January 2017.

Comments: 12 pages, 7 figures, accepted for publication in the ApJ

Journal ref: The Astrophysical Journal, 843, 35 (2017)

Cosmic Chandlery with Thermonuclear Supernovae

Authors: Alan C. Calder, Brendan K. Krueger, Aaron P. Jackson, Don E. Willcox, Broxton J. Miles, Dean M. Townsley

Abstract: Thermonuclear (Type Ia) supernovae are bright stellar explosions, the light curves of which can be calibrated to allow for use as "standard candles" for measuring cosmological distances. Contemporary research investigates how the brightness of an event may be influenced by properties of the progenitor system that follow from properties of the host galaxy such as composition and age. The goals are… ▽ More Thermonuclear (Type Ia) supernovae are bright stellar explosions, the light curves of which can be calibrated to allow for use as "standard candles" for measuring cosmological distances. Contemporary research investigates how the brightness of an event may be influenced by properties of the progenitor system that follow from properties of the host galaxy such as composition and age. The goals are to better understand systematic effects and to assess the intrinsic scatter in the brightness, thereby reducing uncertainties in cosmological studies. We present the results from ensembles of simulations in the single-degenerate paradigm addressing the influence of age and metallicity on the brightness of an event and compare our results to observed variations of brightness that correlate with properties of the host galaxy. We also present results from "hybrid" progenitor models that incorporate recent advances in stellar evolution. △ Less

Submitted 6 December, 2016; originally announced December 2016.

Comments: 14 pages, accepted to the proceedings of Astronum2016 - 11th International Conference on Numerical Modeling of Space Plasma Flows, June 6-10, 2016

GW Librae: Still Hot Eight Years Post-Outburst

Authors: Paula Szkody, Anjum S. Mukadam, Boris T. Gaensicke, Paul Chote, Peter Nelson, Gordon Myers, Odette Toloza, Elizabeth O. Waagen, Edward M. Sion, Denis J. Sullivan, Dean M. Townsley

Abstract: We report continued Hubble Space Telescope (HST) ultraviolet spectra and ground-based optical photometry and spectroscopy of GW Librae eight years after its largest known dwarf nova outburst in 2007. This represents the longest cooling timescale measured for any dwarf nova. The spectra reveal that the white dwarf still remains about 3000 K hotter than its quiescent value. Both ultraviolet and opti… ▽ More We report continued Hubble Space Telescope (HST) ultraviolet spectra and ground-based optical photometry and spectroscopy of GW Librae eight years after its largest known dwarf nova outburst in 2007. This represents the longest cooling timescale measured for any dwarf nova. The spectra reveal that the white dwarf still remains about 3000 K hotter than its quiescent value. Both ultraviolet and optical light curves show a short period of 364-373 s, similar to one of the non-radial pulsation periods present for years prior to the outburst, and with a similar large UV/optical amplitude ratio. A large modulation at a period of 2 h (also similar to that observed prior to outburst) is present in the optical data preceding and during the HST observations, but the satellite observation intervals did not cover the peaks of the optical modulation so it is not possible to determine its corresponding UV amplitude. The similarity of the short and long periods to quiescent values implies the pulsating, fast spinning white dwarf in GW Lib may finally be nearing its quiescent configuration. △ Less

Submitted 2 June, 2016; originally announced June 2016.

Comments: 6 figures, accepted in AJ

A Tracer Method for Computing Type Ia Supernova Yields: Burning Model Calibration, Reconstruction of Thickened Flames, and Verification for Planar Detonations

Authors: Dean M. Townsley, Broxton J. Miles, F. X. Timmes, Alan C. Calder, Edward F. Brown

Abstract: We refine our previously introduced parameterized model for explosive carbon-oxygen fusion during thermonuclear supernovae (SN Ia) by adding corrections to post-processing of recorded Lagrangian fluid element histories to obtain more accurate isotopic yields. Deflagration and detonation products are verified for propagation in a uniform density medium. A new method is introduced for reconstructing… ▽ More We refine our previously introduced parameterized model for explosive carbon-oxygen fusion during thermonuclear supernovae (SN Ia) by adding corrections to post-processing of recorded Lagrangian fluid element histories to obtain more accurate isotopic yields. Deflagration and detonation products are verified for propagation in a uniform density medium. A new method is introduced for reconstructing the temperature-density history within the artificially thick model deflagration front. We obtain better than 5\% consistency between the electron capture computed by the burning model and yields from post-processing. For detonations, we compare to a benchmark calculation of the structure of driven steady-state planar detonations performed with a large nuclear reaction network and error-controlled integration. We verify that, for steady-state planar detonations down to a density of 5x10^6 g/cc, our post processing matches the major abundances in the benchmark solution typically to better than 10% for times greater than 0.01 s after the shock front passage. As a test case to demonstrate the method, presented here with post-processing for the first time, we perform a two dimensional simulation of a SN Ia in the Chandrasekhar-mass deflagration-detonation transition (DDT) scenario. We find that reconstruction of deflagration tracks leads to slightly more complete silicon burning than without reconstruction. The resulting abundance structure of the ejecta is consistent with inferences from spectroscopic studies of observed SNe Ia. We confirm the absence of a central region of stable Fe-group material for the multi-dimensional DDT scenario. Detailed isotopic yields are tabulated and only change modestly when using deflagration reconstruction. △ Less

Submitted 16 May, 2016; originally announced May 2016.

Comments: 28 pages, 16 figures, Accepted to the Astrophysical Journal Supplemental

GW Librae: A unique laboratory for pulsations in an accreting white dwarf

Authors: O. Toloza, B. T. Gaensicke, J. J. Hermes, D. M. Townsley, M. R. Schreiber, P. Szkody, A. Pala, K. Beuermann, L. Bildsten, E. Breedt, M. Cook, P. Godon, A. A. Henden, I. Hubeny, C. Knigge, K. S. Long, T. R. Marsh, D. de Martino, A. S. Mukadam, G. Myers, P. Nelson, A. Oksanen, J. Patterson, E. M. Sion, M. Zorotovic

Abstract: Non-radial pulsations have been identified in a number of accreting white dwarfs in cataclysmic variables. These stars offer insight into the excitation of pulsation modes in atmospheres with mixed compositions of hydrogen, helium, and metals, and the response of these modes to changes in the white dwarf temperature. Among all pulsating cataclysmic variable white dwarfs, GW Librae stands out by ha… ▽ More Non-radial pulsations have been identified in a number of accreting white dwarfs in cataclysmic variables. These stars offer insight into the excitation of pulsation modes in atmospheres with mixed compositions of hydrogen, helium, and metals, and the response of these modes to changes in the white dwarf temperature. Among all pulsating cataclysmic variable white dwarfs, GW Librae stands out by having a well-established observational record of three independent pulsation modes that disappeared when the white dwarf temperature rose dramatically following its 2007 accretion outburst. Our analysis of HST ultraviolet spectroscopy taken in 2002, 2010 and 2011, showed that pulsations produce variations in the white dwarf effective temperature as predicted by theory. Additionally in May~2013, we obtained new HST/COS ultraviolet observations that displayed unexpected behaviour: besides showing variability at ~275s, which is close to the post-outburst pulsations detected with HST in 2010 and 2011, the white dwarf exhibits high-amplitude variability on a ~4.4h time-scale. We demonstrate that this variability is produced by an increase of the temperature of a region on white dwarf covering up to ~30 per cent of the visible white dwarf surface. We argue against a short-lived accretion episode as the explanation of such heating, and discuss this event in the context of non-radial pulsations on a rapidly rotating star △ Less

Submitted 7 April, 2016; originally announced April 2016.

Comments: 11 pages, 8 figures, 2 tables, accpted in MNRAS

White Paper on Nuclear Astrophysics

Authors: Almudena Arcones, Dan W. Bardayan, Timothy C. Beers, Lee A. Berstein, Jeffrey C. Blackmon, Bronson Messer, B. Alex Brown, Edward F. Brown, Carl R. Brune, Art E. Champagne, Alessandro Chieffi, Aaron J. Couture, Pawel Danielewicz, Roland Diehl, Mounib El-Eid, Jutta Escher, Brian D. Fields, Carla Fröhlich, Falk Herwig, William Raphael Hix, Christian Iliadis, William G. Lynch, Gail C. McLaughlin, Bradley S. Meyer, Anthony Mezzacappa , et al. (18 additional authors not shown)

Abstract: This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation… ▽ More This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade. △ Less

Submitted 24 March, 2016; v1 submitted 4 March, 2016; originally announced March 2016.

Comments: 131 pages, community white paper based on 2012 JINA Town Meeting in Detroit, MI, and 2014 APS Town Meeting in College Station, TX, corrected author names

Type Ia Supernova Explosions from Hybrid Carbon-Oxygen-Neon White Dwarf Progenitors

Authors: Donald E. Willcox, Dean M. Townsley, Alan C. Calder, Pavel A. Denissenkov, Falk Herwig

Abstract: Motivated by recent results in stellar evolution that predict the existence of hybrid white dwarf (WD) stars with a C-O core inside an O-Ne shell, we simulate thermonuclear (Type Ia) supernovae from these hybrid progenitors. We use the FLASH code to perform multidimensional simulations in the deflagration to detonation transition (DDT) explosion paradigm. Our hybrid progenitor models were produced… ▽ More Motivated by recent results in stellar evolution that predict the existence of hybrid white dwarf (WD) stars with a C-O core inside an O-Ne shell, we simulate thermonuclear (Type Ia) supernovae from these hybrid progenitors. We use the FLASH code to perform multidimensional simulations in the deflagration to detonation transition (DDT) explosion paradigm. Our hybrid progenitor models were produced with the MESA stellar evolution code and include the effects of the Urca process, and we map the progenitor model to the FLASH grid. We performed a suite of DDT simulations over a range of ignition conditions consistent with the progenitor's thermal and convective structure assuming multiple ignition points. To compare the results from these hybrid WD stars to previous results from C-O white dwarfs, we construct a set of C-O WD models with similar properties and similarly simulate a suite of explosions. We find that despite significant variability within each suite, trends distinguishing the explosions are apparent in their $^{56}$Ni yields and the kinetic properties of the ejecta. We comment on the feasibility of these explosions as the source of some classes of observed subluminous events. △ Less

Submitted 19 February, 2016; originally announced February 2016.

Comments: 14 pages, 19 figures, submitted to the Astrophysical Journal

Seismology of Rapidly Rotating Accreting White Dwarfs

Authors: Dean M. Townsley, Phil Arras, Lars Bildsten

Abstract: A number of White Dwarfs (WDs) in cataclysmic binaries have shown brightness variations consistent with non-radial oscillations as observed in isolated WDs. A few objects have been well-characterized with photometric campaigns in the hopes of gleaning information about the mass, spin, and possibly internal structural characteristics. The novel aspect of this work is the possiblity to measure or co… ▽ More A number of White Dwarfs (WDs) in cataclysmic binaries have shown brightness variations consistent with non-radial oscillations as observed in isolated WDs. A few objects have been well-characterized with photometric campaigns in the hopes of gleaning information about the mass, spin, and possibly internal structural characteristics. The novel aspect of this work is the possiblity to measure or constrain the interior structure and spin rate of WDs which have spent gigayears accreting material from their companion, undergoing thousands of nova outbursts in the process. In addition, variations in the surface temperature affect the site of mode driving, and provide unique and challenging tests for mode driving theories previously applied to isolated WD's. Having undergone long-term accretion, these WDs are expected to have been spun up. Spin periods in the range 60-100 seconds have been measured by other means for two objects, GW Lib and V455 And. Compared to typical mode frequencies, the spin frequency may be similar or higher, and the Coriolis force can no longer be treated as a small perturbation on the fluid motions. We present the results of a non-perturbative calculation of the normal modes of these WDs, using interior thermal structures appropriate to accreting systems. This includes a discussion of the surface brightness distributions, which are strongly modified from the non-rotating case. Using the measured spin period of approximately 100 seconds, we show that the observed pulsations from GW Lib are consistent with the three lowest azimuthal order rotationally modified modes that have the highest frequency in the stellar frame. The high frequencies are needed for the convective driving, but are then apparently shifted to lower frequencies by a combination of their pattern motion and the WD rotation. △ Less

Submitted 8 January, 2016; originally announced January 2016.

Comments: 6 pages, 4 figures, proceedings from 2010 conference "The Physics of Accreting Compact Binaries" Kyoto, Japan. Submitted to proceedings in 2011

On Measuring the Metallicity of a Type Ia Supernova's Progenitor

Authors: Broxton J. Miles, Daniel R. van Rossum, Dean M. Townsley, F. X. Timmes, Aaron P. Jackson, Alan C. Calder, Edward F. Brown

Abstract: In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single de… ▽ More In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition scenario with different $^{56}$Ni yields and the W7 simulation to control for differences between explosion models and total yields. While the overall yields of intermediate mass elements (16 $<$ A $\leq$ 40) differ between the three cases, trends in the yields are similar. With increasing metallicity, $^{28}$Si yields remain nearly constant, $^{40}$Ca yields decline, and Ti and $^{54}$Fe yields increase. In the synthetic spectra, we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: a Ti feature around 4200\,Å and a Fe feature around 5200\,Å\@. In all three simulations, their pseudo equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two features may allow differentiation among progenitor metallicities of observed \sneia and potentially help reduce the intrinsic Hubble scatter. △ Less

Submitted 3 May, 2016; v1 submitted 24 August, 2015; originally announced August 2015.

Comments: 20 pages, 18 figures, accepted by the Astrophysical Journal

Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions

Authors: Bill Paxton, Pablo Marchant, Josiah Schwab, Evan B. Bauer, Lars Bildsten, Matteo Cantiello, Luc Dessart, R. Farmer, H. Hu, N. Langer, R. H. D. Townsend, Dean M. Townsley, F. X. Timmes

Abstract: We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of… ▽ More We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development. △ Less

Submitted 8 January, 2017; v1 submitted 9 June, 2015; originally announced June 2015.

Comments: Published in the ApJ Supplement Series; updated to account for published erratum

Journal ref: ApJS, 220, 15 (2015)

Heavy metals in a light white dwarf: Abundances of the metal-rich, extremely low-mass GALEX J1717+6757

Authors: J. J. Hermes, B. T. Gaensicke, D. Koester, M. C. P. Bours, D. M. Townsley, J. Farihi, T. R. Marsh, Stuart Littlefair, V. S. Dhillon, A. Gianninas, E. Breedt, R. Raddi

Abstract: Using the Hubble Space Telescope, we detail the first abundance analysis enabled by far-ultraviolet spectroscopy of a low-mass (~0.19 Msun) white dwarf (WD), GALEX J1717+6757, which is in a 5.9-hr binary with a fainter, more-massive companion. We see absorption from nine metals, including roughly solar abundances of Ca, Fe, Ti, and P. We detect a significantly sub-solar abundance of C, and put upp… ▽ More Using the Hubble Space Telescope, we detail the first abundance analysis enabled by far-ultraviolet spectroscopy of a low-mass (~0.19 Msun) white dwarf (WD), GALEX J1717+6757, which is in a 5.9-hr binary with a fainter, more-massive companion. We see absorption from nine metals, including roughly solar abundances of Ca, Fe, Ti, and P. We detect a significantly sub-solar abundance of C, and put upper limits on N and O that are also markedly sub-solar. Updated diffusion calculations indicate that all metals should settle out of the atmosphere of this 14,900 K, log(g) = 5.67 WD in the absence of radiative forces in less than 20 yr, orders of magnitude faster than the cooling age of hundreds of Myr. We demonstrate that ongoing accretion of rocky material that is often the cause of atmospheric metals in isolated, more massive WDs is unlikely to explain the observed abundances in GALEX J1717+6757. Using new radiative levitation calculations, we determine that radiative forces can counteract diffusion and support many but not all of the elements present in the atmosphere of this WD; radiative levitation cannot, on its own, explain all of the observed abundance patterns, and additional mechanisms such as rotational mixing may be required. Finally, we detect both primary and secondary eclipses using ULTRACAM high-speed photometry, which we use to constrain the low-mass WD radius and rotation rate as well as update the ephemeris from the discovery observations of this WD+WD binary. △ Less

Submitted 29 July, 2014; originally announced July 2014.

Comments: 10 pages, 5 figures, accepted for publication in MNRAS

Power-Law Wrinkling Turbulence-Flame Interaction Model for Astrophysical Flames

Authors: Aaron P. Jackson, Dean M. Townsley, Alan C. Calder

Abstract: We extend a model for turbulence-flame interactions (TFI) to consider astrophysical flames with a particular focus on combustion in type Ia supernovae. The inertial range of the turbulent cascade is nearly always under-resolved in simulations of astrophysical flows, requiring the use of a model in order to quantify the effects of subgrid-scale wrinkling of the flame surface. We provide implementat… ▽ More We extend a model for turbulence-flame interactions (TFI) to consider astrophysical flames with a particular focus on combustion in type Ia supernovae. The inertial range of the turbulent cascade is nearly always under-resolved in simulations of astrophysical flows, requiring the use of a model in order to quantify the effects of subgrid-scale wrinkling of the flame surface. We provide implementation details to extend a well-tested TFI model to low-Prandtl number flames for use in the compressible hydrodynamics code FLASH. A local, instantaneous measure of the turbulent velocity is calibrated for FLASH and verification tests are performed. Particular care is taken to consider the relation between the subgrid rms turbulent velocity and the turbulent flame speed, especially for high-intensity turbulence where the turbulent flame speed is not expected to scale with the turbulent velocity. Finally, we explore the impact of different TFI models in full-star, three-dimensional simulations of type Ia supernovae. △ Less

Submitted 18 February, 2014; originally announced February 2014.

Comments: 20 pages, 12 figures, accepted to the Astrophysical Journal

KIC11911480: the second ZZ Ceti in the $Kepler$ field

Authors: S. Greiss, B. T. Gaensicke, J. J. Hermes, D. Steeghs, D. Koester, G. Ramsay, T. Barclay, D. M. Townsley

Abstract: We report the discovery of the second pulsating hydrogen-rich (DA) white dwarf in the $Kepler$ field, KIC11911480. It was selected from the $Kepler$-INT Survey (KIS) on the basis of its colours and its variable nature was confirmed using ground-based time-series photometry. An atmosphere model fit to an intermediate-resolution spectrum of KIC11911480 places this DA white dwarf close to the blue ed… ▽ More We report the discovery of the second pulsating hydrogen-rich (DA) white dwarf in the $Kepler$ field, KIC11911480. It was selected from the $Kepler$-INT Survey (KIS) on the basis of its colours and its variable nature was confirmed using ground-based time-series photometry. An atmosphere model fit to an intermediate-resolution spectrum of KIC11911480 places this DA white dwarf close to the blue edge of the empirical boundaries of the ZZ Ceti instability strip: $T_\mathrm{eff} = 12\,160 \pm 250$ K and $\log{g} = 7.94 \pm 0.10 $. Assuming a mass-radius relation and cooling models for DA white dwarfs, the atmospheric parameters yield: M$_{\rm WD}$ = 0.57 $\pm$ 0.06 M$_\odot$. We also obtained two quarters (Q12 and Q16) of nearly uninterrupted short-cadence $Kepler$ data on this star. We detect a total of six independent pulsation modes with a $\geq$ 3$σ$ confidence in its amplitude power spectrum. These pulsations have periods ranging between 172.9 s and 324.5 s, typical of the hotter ZZ Ceti stars. Our preliminary asteroseismic study suggest that KIC11911480 has a rotation rate of 3.5$\pm$0.5 days. △ Less

Submitted 16 December, 2013; originally announced December 2013.

Comments: 8 pages, 5 figures, accepted for publication in MNRAS

Hubble Space Telescope and Ground-Based Observations of V455 Andromedae Post-Outburst

Authors: Paula Szkody, Anjum S. Mukadam, Boris T. Gaensicke, Arne Henden, Edward M. Sion, Dean M. Townsley, Damian Christian, Ross E. Falcon, Stylianos Pyrzas, Justin Brown, Kelsey Funkhouser

Abstract: Hubble Space Telescope spectra obtained in 2010 and 2011, three and four years after the large amplitude dwarf nova outburst of V455 And, were combined with optical photometry and spectra to study the cooling of the white dwarf, its spin, and possible pulsation periods after the outburst. The modeling of the ultraviolet (UV) spectra show that the white dwarf temperature remains ~600 K hotter than… ▽ More Hubble Space Telescope spectra obtained in 2010 and 2011, three and four years after the large amplitude dwarf nova outburst of V455 And, were combined with optical photometry and spectra to study the cooling of the white dwarf, its spin, and possible pulsation periods after the outburst. The modeling of the ultraviolet (UV) spectra show that the white dwarf temperature remains ~600 K hotter than its quiescent value at three years post outburst, and still a few hundred degrees hotter at four years post outburst. The white dwarf spin at 67.6 s and its second harmonic at 33.8 s are visible in the optical within a month of outburst and are obvious in the later UV observations in the shortest wavelength continuum and the UV emission lines, indicating an origin in high temperature regions near the accretion curtains. The UV light curves folded on the spin period show a double-humped modulation consistent with two-pole accretion. The optical photometry two years after outburst shows a group of frequencies present at shorter periods (250-263 s) than the periods ascribed to pulsation at quiescence, and these gradually shift toward the quiescent frequencies (300-360 s) as time progresses past outburst. The most surprising result is that the frequencies near this period in the UV data are only prominent in the emission lines, not the UV continuum, implying an origin away from the white dwarf photosphere. Thus, the connection of this group of periods with non-radial pulsations of the white dwarf remains elusive. △ Less

Submitted 4 September, 2013; originally announced September 2013.

Comments: 18 pages, 14 figs; accepted for ApJ

On Silicon Group Elements Ejected by Supernovae Type Ia

Authors: Soma De, F. X. Timmes, Themis Athanassiadou, David A. Chamulak, Wendy Hawley, Dennis Jack, Edward F. Brown, Alan C. Calder, Dean M. Townsley

Abstract: There is compelling evidence that the peak brightness of a Type Ia supernova is affected by the electron fraction Ye at the time of the explosion. The electron fraction is set by the aboriginal composition of the white dwarf and the reactions that occur during the pre explosive convective burning. To date, determining the makeup of the white dwarf progenitor has relied on indirect proxies, such as… ▽ More There is compelling evidence that the peak brightness of a Type Ia supernova is affected by the electron fraction Ye at the time of the explosion. The electron fraction is set by the aboriginal composition of the white dwarf and the reactions that occur during the pre explosive convective burning. To date, determining the makeup of the white dwarf progenitor has relied on indirect proxies, such as the average metallicity of the host stellar population. In this paper, we present analytical calculations supporting the idea that the electron fraction of the progenitor systematically influences the nucleosynthesis of silicon group ejecta in Type Ia supernovae. In particular, we suggest the abundances generated in quasi nuclear statistical equilibrium are preserved during the subsequent freezeout. This allows one to potential recovery of Ye at explosion from the abundances recovered from an observed spectra. We show that measurement of 28Si, 32S, 40Ca, and 54Fe abundances can be used to construct Ye in the silicon rich regions of the supernovae. If these four abundances are determined exactly, they are sufficient to recover Ye to 6 percent. This is because these isotopes dominate the composition of silicon-rich material and iron rich material in quasi nuclear statistical equilibrium. Analytical analysis shows that the 28Si abundance is insensitive to Ye, the 32S abundance has a nearly linear trend with Ye, and the 40Ca abundance has a nearly quadratic trend with Ye. We verify these trends with post-processing of 1D models and show that these trends are reflected in model synthetic spectra. △ Less

Submitted 8 May, 2014; v1 submitted 17 April, 2013; originally announced April 2013.

Comments: Submitted to the ApJ

The influence of chemical composition on models of Type Ia supernovae

Authors: Alan C. Calder, Brendan K. Krueger, Aaron P. Jackson, Dean M. Townsley

Abstract: Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite the highly successful use of these events in this capacity, many fundamental questions remain. Contemporary research investigates how properties of the progenitor system that follow from the host galaxy such as composition a… ▽ More Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite the highly successful use of these events in this capacity, many fundamental questions remain. Contemporary research investigates how properties of the progenitor system that follow from the host galaxy such as composition and age influence the brightness of an event with the goal of better understanding and assessing the intrinsic scatter in the brightness. We provide an overview of these supernovae and proposed progenitor systems, all of which involve one or more compact stars known as white dwarfs. We describe contemporary research investigating how the composition and structure of the progenitor white dwarf systematically influences the explosion outcome assuming the progenitor is a single white dwarf that has gained mass from a companion. We present results illustrating some of these systematic effects from our research. △ Less

Submitted 9 March, 2013; originally announced March 2013.

Comments: 22 pages, 3 figures, Frontiers of Physics, in press

Evaluating Systematic Dependencies of Type Ia Supernovae: The Influence of Central Density

Authors: Brendan K. Krueger, Aaron P. Jackson, Alan C. Calder, Dean M. Townsley, Edward F. Brown, Francis X. Timmes

Abstract: We present a study exploring a systematic effect on the brightness of type Ia supernovae using numerical models that assume the single-degenerate paradigm. Our investigation varied the central density of the progenitor white dwarf at flame ignition, and considered its impact on the explosion yield, particularly the production and distribution of radioactive Ni-56, which powers the light curve. We… ▽ More We present a study exploring a systematic effect on the brightness of type Ia supernovae using numerical models that assume the single-degenerate paradigm. Our investigation varied the central density of the progenitor white dwarf at flame ignition, and considered its impact on the explosion yield, particularly the production and distribution of radioactive Ni-56, which powers the light curve. We performed a suite of two-dimensional simulations with randomized initial conditions, allowing us to characterize the statistical trends that we present. The simulations indicate that production of Fe-group material is statistically independent of progenitor central density, but the mass of stable Fe-group isotopes is tightly correlated with central density, with a decrease in the production of Ni-56 at higher central densities. These results imply progenitors with higher central densities produce dimmer events. We provide details of the post-explosion distribution of Ni-56 in the models, including the lack of a consistent centrally-located deficit of Ni-56, which may be compared to observed remnants. By performing a self-consistent extrapolation of our model yields and considering the main-sequence lifetime of the progenitor star and the elapsed time between the formation of the white dwarf and the onset of accretion, we develop a brightness-age relation that improves our prediction of the expected trend for single degenerates and we compare this relation with observations. △ Less

Submitted 9 August, 2012; originally announced August 2012.

Comments: 25 pages, 18 figures, accepted to ApJ on July 30, 2012

Journal ref: Astrophysical Journal 757 (2012) 175

Laterally Propagating Detonations in Thin Helium Layers on Accreting White Dwarfs

Authors: Dean M. Townsley, Kevin Moore, Lars Bildsten

Abstract: Theoretical work has shown that intermediate mass (0.01Msun<M_He<0.1Msun) Helium shells will unstably ignite on the accreting white dwarf (WD) in an AM CVn binary. For more massive (M>0.8Msun) WDs, these helium shells can be dense enough (5x10^5 g/cc) that the convectively burning region runs away on a timescale comparable to the sound travel time across the shell; raising the possibility for an e… ▽ More Theoretical work has shown that intermediate mass (0.01Msun<M_He<0.1Msun) Helium shells will unstably ignite on the accreting white dwarf (WD) in an AM CVn binary. For more massive (M>0.8Msun) WDs, these helium shells can be dense enough (5x10^5 g/cc) that the convectively burning region runs away on a timescale comparable to the sound travel time across the shell; raising the possibility for an explosive outcome. The nature of the explosion (i.e. deflagration or detonation) remains ambiguous. In the case of detonation, this causes a laterally propagating front whose properties in these geometrically thin and low density shells we begin to study here. Our calculations show that the radial expansion time of <0.1 s leads to incomplete helium burning, in agreement with recent work by Sim and collaborators, but that the nuclear energy released is still adequate to realize a self-sustaining detonation propagating laterally at slower than the Chapman-Jouguet speed. Our simulations resolve the subsonic region behind the front and are consistent with a direct computation of the reaction structure from the shock strength. The ashes are typically He rich, and consist of predominantly Ti-44, Cr-48, along with a small amount of Fe-52, with very little Ni-56 and with significant Ca-40 in carbon-enriched layers. If this helium detonation results in a Type Ia Supernova, its spectral signatures would appear for the first few days after explosion. (abridged) △ Less

Submitted 29 May, 2012; originally announced May 2012.

Comments: 7 pages, 5 figure, accepted to the Astrophysical Journal

On Simulating Type Ia Supernovae

Authors: Alan C. Calder, Brendan K. Krueger, Aaron P. Jackson, Dean M. Townsley, Edward F. Brown, Francis X. Timmes

Abstract: Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite their highly successful use in this capacity, the progenitors of these events are incompletely understood. We describe simulating type Ia supernovae in the paradigm of a thermonuclear runaway occurring in a massive white dwa… ▽ More Type Ia supernovae are bright stellar explosions distinguished by standardizable light curves that allow for their use as distance indicators for cosmological studies. Despite their highly successful use in this capacity, the progenitors of these events are incompletely understood. We describe simulating type Ia supernovae in the paradigm of a thermonuclear runaway occurring in a massive white dwarf star. We describe the multi-scale physical processes that realistic models must incorporate and the numerical models for these that we employ. In particular, we describe a flame-capturing scheme that addresses the problem of turbulent thermonuclear combustion on unresolved scales. We present the results of our study of the systematics of type Ia supernovae including trends in brightness following from properties of the host galaxy that agree with observations. We also present performance results from simulations on leadership-class architectures. △ Less

Submitted 4 May, 2012; originally announced May 2012.

Comments: 13 pages, 3 figures, accepted to proceedings of the Conference on Computational Physics, Oct. 30 - Nov. 3, 2011

The Detonation Mechanism of the Pulsationally-Assisted Gravitationally-Confined Detonation Model of Type Ia Supernovae

Authors: G. C. Jordan IV, C. Graziani, R. T. Fisher, D. M. Townsley, C. Meakin, K. Weide, L. B. Reid, J. Norris, R. Hudson, D. Q. Lamb

Abstract: We describe the detonation mechanism comprising the "Pulsationally Assisted" Gravitationally Confined Detonation (GCD) model of Type Ia supernovae SNe Ia. This model is analogous to the previous GCD model reported in Jordan et al.(2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phas… ▽ More We describe the detonation mechanism comprising the "Pulsationally Assisted" Gravitationally Confined Detonation (GCD) model of Type Ia supernovae SNe Ia. This model is analogous to the previous GCD model reported in Jordan et al.(2008); however, the chosen initial conditions produce a substantively different detonation mechanism, resulting from a larger energy release during the deflagration phase. The resulting final kinetic energy and nickel-56 yields conform better to observational values than is the case for the "classical" GCD models. In the present class of models, the ignition of a deflagration phase leads to a rising, burning plume of ash. The ash breaks out of the surface of the white dwarf, flows laterally around the star, and converges on the collision region at the antipodal point from where it broke out. The amount of energy released during the deflagration phase is enough to cause the star to rapidly expand, so that when the ash reaches the antipodal point, the surface density is too low to initiate a detonation. Instead, as the ash flows into the collision region (while mixing with surface fuel), the star reaches its maximally expanded state and then contracts. The stellar contraction acts to increase the density of the star, including the density in the collision region. This both raises the temperature and density of the fuel-ash mixture in the collision region and ultimately leads to thermodynamic conditions that are necessary for the Zel'dovich gradient mechanism to produce a detonation. We demonstrate feasibility of this scenario with three 3-dimensional (3D), full star simulations of this model using the FLASH code. We characterized the simulations by the energy released during the deflagration phase, which ranged from 38% to 78% of the white dwarf's binding energy. We show that the necessary conditions for detonation are achieved in all three of the models. △ Less

Submitted 17 September, 2012; v1 submitted 17 February, 2012; originally announced February 2012.

Comments: 22 pages, 8 figures; ApJ

Evaluating Systematic Dependencies of Type Ia Supernovae

Authors: Alan C. Calder, Brendan K. Krueger, Aaron P. Jackson, Dean M. Townsley, Francis X. Timmes, Edward F. Brown, David A. Chamulak

Abstract: Type Ia supernovae are bright stellar explosions thought to occur when a thermonuclear runaway consumes roughly a solar mass of degenerate stellar material. These events produce and disseminate iron-peak elements, and properties of their light curves allow for standardization and subsequent use as cosmological distance indicators. The explosion mechanism of these events remains, however, only part… ▽ More Type Ia supernovae are bright stellar explosions thought to occur when a thermonuclear runaway consumes roughly a solar mass of degenerate stellar material. These events produce and disseminate iron-peak elements, and properties of their light curves allow for standardization and subsequent use as cosmological distance indicators. The explosion mechanism of these events remains, however, only partially understood. Many models posit the explosion beginning with a deflagration born near the center of a white dwarf that has gained mass from a stellar companion. In order to match observations, models of this single-degenerate scenario typically invoke a subsequent transition of the (subsonic) deflagration to a (supersonic) detonation that rapidly consumes the star. We present an investigation into the systematics of thermonuclear supernovae assuming this paradigm. We utilize a statistical framework for a controlled study of two-dimensional simulations of these events from randomized initial conditions. We investigate the effect of the composition and thermal history of the progenitor on the radioactive yield, and thus brightness, of an event. Our results offer an explanation for some observed trends of mean brightness with properties of the host galaxy. △ Less

Submitted 17 November, 2010; originally announced November 2010.

Comments: 14 pages, to appear in the Proceedings of the SciDAC 2010 meeting

Evaluating Systematic Dependencies of Type Ia Supernovae: The Influence of Deflagration to Detonation Density

Authors: Aaron P. Jackson, Alan C. Calder, Dean M. Townsley, David A. Chamulak, Edward F. Brown, F. X. Timmes

Abstract: We explore the effects of the deflagration to detonation transition (DDT) density on the production of Ni-56 in thermonuclear supernova explosions (type Ia supernovae). Within the DDT paradigm, the transition density sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear statistical equilibrium (NSE) material produced. We employ a theoretica… ▽ More We explore the effects of the deflagration to detonation transition (DDT) density on the production of Ni-56 in thermonuclear supernova explosions (type Ia supernovae). Within the DDT paradigm, the transition density sets the amount of expansion during the deflagration phase of the explosion and therefore the amount of nuclear statistical equilibrium (NSE) material produced. We employ a theoretical framework for a well-controlled statistical study of two-dimensional simulations of thermonuclear supernovae with randomized initial conditions that can, with a particular choice of transition density, produce a similar average and range of Ni-56 masses to those inferred from observations. Within this framework, we utilize a more realistic "simmered" white dwarf progenitor model with a flame model and energetics scheme to calculate the amount of Ni-56 and NSE material synthesized for a suite of simulated explosions in which the transition density is varied in the range 1-3x10^7 g/cc. We find a quadratic dependence of the NSE yield on the log of the transition density, which is determined by the competition between plume rise and stellar expansion. By considering the effect of metallicity on the transition density, we find the NSE yield decreases by 0.055 +/- 0.004 solar masses for a 1 solar metallicity increase evaluated about solar metallicity. For the same change in metallicity, this result translates to a 0.067 +/- 0.004 solar mass decrease in the Ni-56 yield, slightly stronger than that due to the variation in electron fraction from the initial composition. Observations testing the dependence of the yield on metallicity remain somewhat ambiguous, but the dependence we find is comparable to that inferred from some studies. △ Less

Submitted 7 July, 2010; originally announced July 2010.

Comments: 15 pages, 13 figures, accepted to ApJ on July 6, 2010

Journal ref: Astrophys.J. 720 (2010) 99-113

On variations of the brightness of type Ia supernovae with the age of the host stellar population

Authors: Brendan K. Krueger, Aaron P. Jackson, Dean M. Townsley, Alan C. Calder, Edward F. Brown, F. X. Timmes

Abstract: Recent observational studies of type Ia supernovae (SNeIa) suggest correlations between the peak brightness of an event and the age of the progenitor stellar population. This trend likely follows from properties of the progenitor white dwarf (WD), such as central density, that follow from properties of the host stellar population. We present a statistically well-controlled, systematic study utiliz… ▽ More Recent observational studies of type Ia supernovae (SNeIa) suggest correlations between the peak brightness of an event and the age of the progenitor stellar population. This trend likely follows from properties of the progenitor white dwarf (WD), such as central density, that follow from properties of the host stellar population. We present a statistically well-controlled, systematic study utilizing a suite of multi-dimensional SNeIa simulations investigating the influence of central density of the progenitor WD on the production of Fe-group material, particularly radioactive Ni-56, which powers the light curve. We find that on average, as the progenitor's central density increases, production of Fe-group material does not change but production of Ni-56 decreases. We attribute this result to a higher rate of neutronization at higher density. The central density of the progenitor is determined by the mass of the WD and the cooling time prior to the onset of mass transfer from the companion, as well as the subsequent accretion heating and neutrino losses. The dependence of this density on cooling time, combined with the result of our central density study, offers an explanation for the observed age-luminosity correlation: a longer cooling time raises the central density at ignition thereby producing less Ni-56 and thus a dimmer event. While our ensemble of results demonstrates a significant trend, we find considerable variation between realizations, indicating the necessity for averaging over an ensemble of simulations to demonstrate a statistically significant result. △ Less

Submitted 6 July, 2010; originally announced July 2010.

Comments: 5 pages, 4 figures, 1 table, accepted to ApJL

Journal ref: Astrophysical Journal Letters 719 (2010) L5-L9