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.

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.

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.

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.

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.

Interacting binaries containing white dwarfs can lead to a variety of outcomes that range from powerful thermonuclear explosions, which are important in the chemical evolution of galaxies and as cosmological distance estimators, to strong sources of low frequency gravitational wave radiation, which makes them ideal calibrators for the gravitational low-frequency wave detector LISA mission. However, current theoretical evolution models still fail to explain the observed properties of the known populations of white dwarfs in both interacting and detached binaries. Major limitations are that the existing population models have generally been developed to explain the properties of sub-samples of these systems, occupying small volumes of the vast parameter space, and that the observed samples are severely biased. The overarching goal for the next decade is to assemble a large and homogeneous sample of white dwarf binaries that spans the entire range of evolutionary states, to obtain precise measurements of their physical properties, and to further develop the theory to satisfactorily reproduce the properties of the entire population. While ongoing and future all-sky high- and low-resolution optical spectroscopic surveys allow us to enlarge the sample of these systems, high-resolution ultraviolet spectroscopy is absolutely essential for the characterization of the white dwarfs in these binaries. The Hubble Space Telescope is currently the only facility that provides ultraviolet spectroscopy, and with its foreseeable demise, planning the next ultraviolet mission is of utmost urgency.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Recent Suzaku X-ray spectra of SNR 3C 397 indicate enhanced stable iron-group element abundances of Ni, Mn, Cr, and Fe. Seeking to address key questions about the progenitor and explosion mechanism of 3C 397, we compute nucleosynthetic yields from a suite of multidimensional hydrodynamics models in the near-Chandrasekhar mass, single-degenerate paradigm for supernova Type Ia. Varying the progenitor white dwarf internal structure, composition, ignition, and explosion mechanism, we find the best match to the observed iron-peak elements of 3C 397 are dense (central density $\ge$ 6$\times$10$^{9}$ g cm$^{-3}$), low-carbon white dwarfs that undergo a weak, centrally-ignited deflagration, followed by a subsequent detonation. The amount of $^{56}$Ni produced is consistent with a normal or bright normal supernova Type Ia. A pure deflagration of a centrally-ignited, low central density ($\simeq$ 2$\times$10$^{9}$ g cm$^{-3}$) progenitor white dwarf, frequently considered in the literature, is also found to produce good agreement with 3C 397 nucleosynthetic yields, but leads to a subluminous SN Ia event, in conflict with X-ray linewidth data. Additionally, in contrast to prior work which suggested a large super-solar metallicity for the white dwarf progenitor for SNR 3C 397, we find satisfactory agreement for solar and sub-solar metallicity progenitors. We discuss a range of implications our results have for the single-degenerate channel.

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.

We present HST spectroscopy for 45 cataclysmic variables (CVs), observed with HST/COS and HST/STIS. For 36 CVs, the white dwarf is recognisable through its broad Ly$\alpha$ absorption profile and we measure the white dwarf effective temperatures ($T_{\mathrm{eff}}$) by fitting the HST data assuming $\log\,g=8.35$, which corresponds to the average mass for CV white dwarfs ($\simeq\,0.8\,\mathrm{M}_\odot$). Our results nearly double the number of CV white dwarfs with an accurate temperature measurement. We find that CVs above the period gap have, on average, higher temperatures ($\langle T_{\mathrm{eff}} \rangle \simeq 23\,000\,$K) and exhibit much more scatter compared to those below the gap ($\langle T_{\mathrm{eff}} \rangle \simeq 15\,000\,$K). While this behaviour broadly agrees with theoretical predictions, some discrepancies are present: (i) all our new measurements above the gap are characterised by lower temperatures ($T_{\mathrm{eff}} \simeq 16\,000 - 26\,000\,$K) than predicted by the present day CV population models ($T_{\mathrm{eff}} \simeq 38\,000 - 43\,000\,$K); (ii) our results below the gap are not clustered in the predicted narrow track and exhibit in particular a relatively large spread near the period minimum, which may point to some shortcomings in the CV evolutionary models. Finally, in the standard model of CV evolution, reaching the minimum period, CVs are expected to evolve back towards longer periods with mean accretion rates $\dot{M}\lesssim 2 \times 10^{-11}\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$, corresponding to $T_\mathrm{eff}\lesssim 11\,500\,$K. We do not unambiguously identify any such system in our survey, suggesting that this major component of the predicted CV population still remains elusive to observations.

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.

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.

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.

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.

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\,\AA\ and a Fe feature around 5200\,\AA\@. 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.

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

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.

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.

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.

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.

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.

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.

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$\sigma$ 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.

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.

Accreted helium layers on white dwarfs have been highlighted for many decades as a possible site for a detonation triggered by a thermonuclear runaway. In this paper, we find the minimum helium layer thickness that will sustain a steady laterally propagating detonation and show that it depends on the density and composition of the helium layer, specifically C12 and O16. Detonations in these thin helium layers have speeds slower than the Chapman-Jouget (CJ) speed from complete helium burning, v_cj=1.5e9 cm/s. Though gravitationally unbound, the ashes still have unburned helium (~ 80% in the thinnest cases) and only reach up to heavy elements such as Ca40, Ti44, Cr48, and Fe52. It is rare for these thin shells to generate large amounts of Ni56. We also find a new set of solutions that can propagate in even thinner helium layers when O16 is present at a minimum mass fraction of ~0.07. Driven by energy release from alpha captures on O16 and subsequent elements, these slow detonations only create ashes up to Si28 in the outer detonated He shell. We close by discussing how the unbound helium burning ashes may create faint and fast ".Ia" supernovae as well as events with virtually no radioactivity, and speculate on how the slower helium detonation velocities impact the off-center ignition of a carbon detonation that could cause a Type Ia supernova in the double detonation scenario.

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.

We present an ultraviolet spectrum and light curve of the short orbital period cataclysmic variable EZ Lyn obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope 14 months after its dwarf nova outburst, along with ground-based optical photometry. The UV spectrum can be fit with a 13,100K, log g=8 white dwarf using 0.5 solar composition, while fits to the individual lines are consistent with solar abundance for Si and Al, but only 0.3 solar for C. The Discrete Fourier Transforms of the UV and optical light curves at 14 months following outburst show a prominent period at 256 sec. This is the same period reported by Pavlenko in optical data obtained 7 months and one year after outburst, indicating its long term stability over several months, but this period is not evident in the pre-outburst data and is much shorter than the 12.6 min period that was seen in observations obtained during an interval from 8 months to 2.5 years after the 2006 outburst. In some respects, the long and short periods are similar to the behavior seen in GW Lib after its outburst but the detailed explanation for the appearance and disappearance of these periods and their relation to non-radial pulsation modes remain to be explored with theoretical models.

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.

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.

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)

Since the large amplitude 2007 outburst which heated its accreting, pulsating white dwarf, the dwarf nova system GW Librae has been cooling to its quiescent temperature. Our Hubble Space Telescope ultraviolet spectra combined with ground-based optical coverage during the 3rd and 4th year after outburst show that the fluxes and temperatures are still higher than quiescence (T=19,700K and 17,300K vs 16,000K pre-outburst for a log g=8.7 and d=100 pc). The K{wd} of 7.6+/-0.8 km/s determined from the CI1463 absorption line, as well as the gravitational redshift implies a white dwarf mass of 0.79+/-0.08 Msun. The widths of the UV lines imply a white dwarf rotation velocity vsin i of 40 km/s and a spin period of 209 s (for an inclination of 11 deg and a white dwarf radius of 7x10^{8} cm). Light curves produced from the UV spectra in both years show a prominent multiplet near 290 s, with higher amplitude in the UV compared to the optical, and increased amplitude in 2011 vs 2010. As the presence of this set of periods is intermittent in the optical on weekly timescales, it is unclear how this relates to the non-radial pulsations evident during quiescence.

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.

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.

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.

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.

Nonradial pulsations in the primary white dwarfs of cataclysmic variables can now potentially allow us to explore the stellar interior of these accretors using stellar seismology. In this context, we conducted a multi-site campaign on the accreting pulsator SDSS J161033.64-010223.3 (V386 Ser) using seven observatories located around the world in May 2007 over a duration of 11 days. We report the best fit periodicities here, which were also previously observed in 2004, suggesting their underlying stability. Although we did not uncover a sufficient number of independent pulsation modes for a unique seismological fit, our campaign revealed that the dominant pulsation mode at 609s is an evenly spaced triplet. The even nature of the triplet is suggestive of rotational splitting, implying an enigmatic rotation period of about 4.8 days. There are two viable alternatives assuming the triplet is real: either the period of 4.8 days is representative of the rotation period of the entire star with implications for the angular momentum evolution of these systems, or it is perhaps an indication of differential rotation with a fast rotating exterior and slow rotation deeper in the star. Investigating the possibility that a changing period could mimic a triplet suggests that this scenario is improbable, but not impossible. Using time-series spectra acquired in May 2009, we determine the orbital period of SDSS J161033.64-010223.3 to be 83.8 +/- 2.9 min. Three of the observed photometric frequencies from our May 2007 campaign appear to be linear combinations of the 609s pulsation mode with the first harmonic of the orbital period at 41.5min. This is the first discovery of a linear combination between nonradial pulsation and orbital motion for a variable white dwarf.

We present a theoretical framework for formal study of systematic effects in Supernovae Type Ia (SN Ia) that utilizes 2-d simulations to implement a form of the deflagration-detonation transition (DDT) explosion scenario. The framework is developed from a randomized initial condition that leads to a sample of simulated SN Ia whose Ni56 masses have a similar average and range to those observed, and have many other modestly realistic features such as the velocity extent of intermediate mass elements. The intended purpose is to enable statistically well-defined studies of both physical and theoretical parameters of the SN Ia explosion simulation. We present here a thorough description of the outcome of the SN Ia explosions produced by our current simulations. A first application of this framework is utilized to study the dependence of the SN Ia on the Ne22 content, which is known to be directly influenced by the progenitor stellar population's metallicity. Our study is very specifically tailored to measure how the Ne22 content influences the competition between the rise of plumes of burned material and the expansion of the star before these plumes reach DDT conditions. This competition controls the amount of material in nuclear statistical equilibrium (NSE) and therefore Ni56 produced by setting the density at which nucleosynthesis takes place during the detonation phase of the explosion. Although the outcome following from any particular ignition condition can change dramatically with Ne22 content, with a sample of 20 ignition conditions we find that the systematic change in the expansion of the star prior to detonation is not large enough to compete with the dependence on initial neutron excess discussed by Timmes, Brown & Truran (2003). (Abridged)

Some explosion models for Type Ia supernovae (SN Ia), such as the gravitationally confined detonation (GCD) or the double detonation sub-Chandrasekhar (DDSC) models, rely on the spontaneous initiation of a detonation in the degenerate C/O material of a white dwarf. The length scales pertinent to the initiation of the detonation are notoriously unresolved in multi-dimensional stellar simulations, prompting the use of results of 1D simulations at higher resolution, such as the ones performed for this work, as guidelines for deciding whether or not conditions reached in the higher dimensional full star simulations successfully would lead to the onset of a detonation. Spontaneous initiation relies on the existence of a suitable gradient in self-ignition (induction) times of the fuel, which we set up with a spatially localized non-uniformity of temperature -- a hot spot. We determine the critical (smallest) sizes of such hot spots that still marginally result in a detonation in white dwarf matter by integrating the reactive Euler equations with the hydrodynamics code FLASH. We quantify the dependences of the critical sizes of such hot spots on composition, background temperature, peak temperature, geometry, and functional form of the temperature disturbance, many of which were hitherto largely unexplored in the literature. We discuss the implications of our results in the context of modeling of SNe Ia.

We review the most decisive currently available measurements of the surface effective temperatures, Teff, of white dwarf (WD) primaries in cataclysmic variables (CVs) during accretion quiescence, and use these as a diagnostic for their time averaged accretion rate, <Mdot>. Using time-dependent calculations of the WD envelope, we investigate the sensitivity of the quiescent Teff to long term variations in the accretion rate. We find that the quiescent Teff provides one of the best available tests of predictions for the angular momentum loss and resultant mass transfer rates which govern the evolution of CVs. While gravitational radiation is sufficient to explain the <Mdot> of strongly magnetic CVs at all Porb, faster angular momentum loss is required by the temperatures of dwarf nova primaries (non-magnetic systems). This provides evidence that a normal stellar magnetic field structure near the secondary is essential for the enhanced braking mechanism to work, supporting the well-known stellar wind braking hypothesis. The contrast in <Mdot> is most prominent for orbital periods Porb > 3 hours, above the period gap, but a modest enhancement is also present at shorter Porb. The averaging time which <Mdot> reflects is as much as 10^5 years for low-<Mdot> systems and as little as 10^3 years for high-<Mdot> systems. We discuss the security of conclusions drawn about the CV population in light of these time scales and our necessarily incomplete sample of systems. Measurements for non-magnetic systems above the period gap fall below predictions from traditional stellar wind braking prescriptions, but above more recent predictions with somewhat weaker angular momentum loss. We also discuss the apparently high Teff's found in the VY Scl stars. (abridged)

We study the gravitationally confined detonation (GCD) model of Type Ia supernovae through the detonation phase and into homologous expansion. In the GCD model, a detonation is triggered by the surface flow due to single point, off-center flame ignition in carbon-oxygen white dwarfs. The simulations are unique in terms of the degree to which non-idealized physics is used to treat the reactive flow, including weak reaction rates and a time dependent treatment of material in nuclear statistical equilibrium (NSE). Careful attention is paid to accurately calculating the final composition of material which is burned to NSE and frozen out in the rapid expansion following the passage of a detonation wave over the high density core of the white dwarf; and an efficient method for nucleosynthesis post-processing is developed which obviates the need for costly network calculations along tracer particle thermodynamic trajectories. Observational diagnostics are presented for the explosion models, including abundance stratifications and integrated yields. We find that for all of the ignition conditions studied here, a self regulating process comprised of neutronization and stellar expansion results in final \iso{Ni}{56} masses of $\sim$1.1\msun. But, more energetic models result in larger total NSE and stable Fe peak yields. The total yield of intermediate mass elements is $\sim0.1$\msun and the explosion energies are all around 1.5$\times10^{51}$ ergs. The explosion models are briefly compared to the inferred properties of recent Type Ia supernova observations. The potential for surface detonation models to produce lower luminosity (lower \iso{Ni}{56} mass) supernovae is discussed.

We present time-series VLT spectroscopy and NTT photometry of the cataclysmic variable SDSS J220553.98+115553.7, which contains a pulsating white dwarf. We determine a spectroscopic orbital period of 82.825 +/- 0.089 min from velocity measurements of the Halpha emission line. A period analysis of the light curves reveals a dominant periodicity at 44.779 +/- 0.038 min which is not related to the spectroscopic period. However, the light curves do not exhibit a variation at any frequency which is attributable to GW Lib type pulsations, to a detection limit of 5 mmag. This non-detection is in contrast to previous studies which have found three pulsation frequencies with amplitudes of 9 to 11 mmag at optical wavelengths. Destructive interference and changes to the thermal properties of the driving layer in direct response to accretion can be ruled out as causes of the non-detection. Alternatively, it is feasible that the object has cooled out of the instability strip since a previous (unobserved) dwarf nova superoutburst. This would be the first time this behaviour has been seen in a cataclysmic variable pulsator. Another possibility is that changes in the surface characteristics, possibly induced by accretion phenomena, have modified the surface visibility of the pulsation modes. Further observations, particularly improved constraints on the timescale for changes in the mode spectrum, are needed to distinguish among possible explanations.

We report the results of a series of three-dimensional (3-D) simulations of the deflagration phase of the gravitationally confined detonation mechanism for Type Ia supernovae. In this mechanism, ignition occurs at one or several off-center points, resulting in a burning bubble of hot ash that rises rapidly, breaks through the surface of the star, and collides at a point opposite breakout on the stellar surface. We find that detonation conditions are robustly reached in our 3-D simulations for a range of initial conditions and resolutions. Detonation conditions are achieved as the result of an inwardly-directed jet that is produced by the compression of unburnt surface material when the surface flow collides with itself. A high-velocity outwardly-directed jet is also produced. The initial conditions explored in this paper lead to conditions at detonation that can be expected to produce large amounts of $^{56}$Ni and small amounts of intermediate mass elements. These particular simulations are therefore relevant only to high luminosity Type Ia supernovae. Recent observations of Type Ia supernovae imply a compositional structure that is qualitatively consistent with that expected from these simulations.

We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova, and apply it to study the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified 3-stage burning model and a non-static ash state are integrated with an artificially thickened advection-diffusion-reaction (ADR) flame front in order to provide an accurate but highly efficient representation of the energy release and electron capture in and after the unresolvable flame. We demonstrate that both our ADR and energy release methods do not generate significant acoustic noise, as has been a problem with previous ADR-based schemes. We proceed to model aspects of the deflagration, particularly the role of buoyancy of the hot ash, and find that our methods are reasonably well-behaved with respect to numerical resolution. We show that if a detonation occurs in material swept up by the material ejected by the first rising bubble but gravitationally confined to the white dwarf (WD) surface (the GCD paradigm), the density structure of the WD at detonation is systematically correlated with the distance of the deflagration ignition point from the center of the star. Coupled to a suitably stochastic ignition process, this correlation may provide a plausible explanation for the variety of nickel masses seen in Type Ia Supernovae.

We develop and calibrate a realistic model flame for hydrodynamical simulations of deflagrations in white dwarf (Type Ia) supernovae. Our flame model builds on the advection-diffusion-reaction model of Khokhlov and includes electron screening and Coulomb corrections to the equation of state in a self-consistent way. We calibrate this model flame--its energetics and timescales for energy release and neutronization--with self-heating reaction network calculations that include both these Coulomb effects and up-to-date weak interactions. The burned material evolves post-flame due to both weak interactions and hydrodynamic changes in density and temperature. We develop a scheme to follow the evolution, including neutronization, of the NSE state subsequent to the passage of the flame front. As a result, our model flame is suitable for deflagration simulations over a wide range of initial central densities and can track the temperature and electron fraction of the burned material through the explosion and into the expansion of the ejecta.

(Abridged) The Cataclysmic Variable (CV) population harbors a diverse range of donor stars and accreting white dwarfs (WDs). A range of WD masses is expected, from low mass Helium core WDs, to massive WDs which have previously accreted at rates high enough for Hydrogen to burn steadily. Furthermore, a wide range of Helium enrichment is expected in the accreted material depending on the degree to which the donor star is evolved. We investigate the impact of this diversity on the range of effective temperatures ($T_{\rm eff}$) for which g-modes are unstable. The critical $T_{\rm eff}$ below which modes are unstable ("blue edge") depends on both surface gravity, $g$, and He abundance, $Y$. The Hydrogen/first Helium ionization instability strip is more sensitive to $g$ than $Y$. We find that (for solar composition envelopes), relative to a fiducial WD mass $0.6 M_\odot$, the blue edge for a $0.4 M_\odot$ He core WD shifts downward by $\approx 1000 {\rm K}$, while that for a massive $\approx 1.2 M_\odot$ WD shifts upward by $\approx 2000 {\rm K}$. The second Helium ionization instability strip exhibits strong dependences on both $g$ and $Y$. Surprisingly, increasing $Y$ by only 10% relative to solar creates an instability strip near $15,000 {\rm K}$. Hence CV's below the period gap with evolved donor stars of $Y\ga 0.4$ may have an "intermediate" instability strip well outside of the isolated DA and DB variables. This "intermediate" instability strip also occurs for low mass He WD with solar composition envelopes. The lack of pulsations in CV's with $T_{\rm eff}$ in the pure Hydrogen ZZ Ceti instability strip is also easily explained.

We calculate the heating and cooling of the accreting white dwarf (WD) in the ultracompact AM Canum Venaticorum (AM CVn) binaries and show that the WD can contribute significantly to their optical and ultraviolet emission. We estimate the WD's effective temperature, Teff, using the optical continuum for a number of observed binaries, and show that it agrees well with our theoretical calculations. Driven by gravitational radiation losses, the time averaged accretion rate, <Mdot>, decreases monotonically with increasing Porb, covering six orders of magnitude. If the short period (Porb<10 min) systems accrete at a rate consistent with gravitational radiation via direct impact, we predict their unpulsed optical/UV light to be that of the Teff>50,000$ K accreting WD. At longer Porb we calculate the Teff and absolute visual magnitude, M_V, that the accreting WD will have during low accretion states, and find that the WD naturally crosses the pulsational instability strip. Discovery and study of pulsations could allow for the measurement of the accumulated helium mass on the accreting WD, as well as its rotation rate. Accretion heats the WD core, but for Porb>40 minutes, the WD's Teff is set by its cooling as <Mdot> plummets. For the two long period AM CVn binaries with measured parallaxes, GP Com and CE 315, we show that the optical broadband colors and intensity are that expected from a pure helium atmosphere WD. This confirms that the WD brightness sets the minimum light in wide AM CVn binaries, allowing for meaningful constraints on their population density from deep optical searches, both in the field and in Globular Clusters.

Classical Novae (CNe) are the brightest manifestation of mass transfer onto a white dwarf in a cataclysmic variable (CV). As such, they are probes of the mass transfer rate, Mdot, and WD mass, Mwd, in these interacting binaries. Our calculations of the dependence of the CN ignition mass, Mign, on Mdot and Mwd yields the recurrence times of these explosions. We show that the observed CNe orbital period distribution is consistent with the interrupted magnetic braking evolutionary scenario, where at orbital periods Porb > 3 hr mass transfer is driven by angular momentum loss via a wind from the companion star and at Porb < 3 hr by gravitational radiation. About 50% of CNe occur in binaries accreting at Mdot ~= 10^{-9} Msun/yr with Porb = 3-4 hr, with the remaining 50% split evenly between Porb longer (higher Mdot) and shorter (lower Mdot) than this. This resolution of the relative contribution to the CN rate from different CVs tells us that 3(9)x10^5 CVs with WD mass 1.0(0.6)Msun are needed to produce one CN per year. Using the K-band specific CN rate measured in external galaxies, we find a CV birthrate of 2(4)x10^{-4}/yr per 10^{10}Lsun,K, very similar to the luminosity specific Type Ia supernova rate in elliptical galaxies. Likewise, we predict that there should be 60-180 CVs for every 10^6Lsun,K in an old stellar population, similar to the number of X-ray identified CVs in the globular cluster 47 Tuc, showing no overabundance relative to the field. Using a two-component steady state model of CV evolution we show that the fraction of CVs which are magnetic (22%) implies a birthrate of 8% relative to non-magnetic CVs, similar to the fraction of strongly magnetic field WDs. (abridged)

Pulsation modes have recently been observed in a handful of white dwarf (WD) primaries of cataclysmic variables, allowing an interesting new probe into the structure of accreting WD's. We briefly discuss the seismology of these objects, how stellar properties may be inferred from the observed mode frequencies, and mode driving mechanisms. For one pulsator, GW Lib, we have shown that a WD mass $M=1.05 M_\odot$ and accreted envelope mass $M_{\rm env}=0.4\times 10^{-4} M_\odot$ give the best match to the observed pulsation periods. A first exploration of mode driving favors $T_{\rm eff} = 14 000 K$ and a massive WD, but more work is necessary.

We present a first analysis of the g-mode oscillation spectrum for the white dwarf (WD) primary of GW Lib, a faint cataclysmic variable (CV). Stable periodicities have been observed from this WD for a number of years, but their interpretation as stellar pulsations has been hampered by a lack of theoretical models appropriate to an accreting WD. Using the results of Townsley and Bildsten, we construct accreting models for the observed effective temperature and approximate mass of the WD in GW Lib. We compute g-mode frequencies for a range of accreted layer masses, Macc, and long term accretion rates, <Mdot>. If we assume that the observed oscillations are from l=1 g-modes, then the observed periods are matched when M ~= 1.02 Msun, Macc ~= 0.31 x 10^-4 Msun and <Mdot> ~= 7.3 x 10^-11 Msun/yr. Much more sensitive observations are needed to discover more modes, after which we will be able to more accurately measure these parameters and constrain or measure the WD's rotation rate.

Our previous theoretical study of the impact of an accreting envelope on the thermal state of an underlying white dwarf (WD) has yielded equilibrium core temperatures, classical nova ignition masses and thermal luminosities for WDs accreting at time averaged rates of <Mdot> = 10^-11 - 10^-8 Msun/yr. These <Mdot>'s are appropriate to WDs in cataclysmic variables (CVs) of P_orb <~ 7 hr, many of which accrete sporadically as Dwarf Novae. Approximately thirty nonmagnetic Dwarf Novae have been observed in quiescence, when the accretion rate is low enough for spectral detection of the WD photosphere, and a measurement of T_eff. We use our theoretical work to translate the measured T_eff's into local time-averaged accretion rates, confirming the factor of ten drop in <Mdot> predicted for CV's as they transit the period gap. For DN below the period gap, we show that if <Mdot> is that given by gravitational radiation losses alone, then the WD masses are > 0.8 Msun. An alternative conclusion is that the masses are closer to 0.6 Msun and <Mdot> is 3-4 times larger than that expected from gravitational radiation losses. In either case, it is very plausible that a subset of CVs with P_orb < 2 hours will have T_eff's low enough for them to become non-radial pulsators, as discovered by van Zyl and collaborators in GW Lib.

White dwarfs experience a thermal renaissance when they receive mass from a stellar companion in a binary. For accretion rates < 10^-8 Msun/yr, the freshly accumulated hydrogen/helium envelope ignites in a thermally unstable manner that results in a classical novae (CN) outburst and ejection of material. We have undertaken a theoretical study of the impact of the accumulating envelope on the thermal state of the underlying white dwarf (WD). This has allowed us to find the equilibrium WD core temperatures (T_c), the classical nova ignition masses (M_ign) and the thermal luminosities for WDs accreting at rates of 10^-11 - 10^-8 Msun/yr. These accretion rates are most appropriate to WDs in cataclysmic variables (CVs) of P_orb <~ 7 hr, many of which accrete sporadically as dwarf novae. We have included ^3He in the accreted material at levels appropriate for CVs and find that it significantly modifies the CN ignition mass. We compare our results with several others from the CN literature and find that the inclusion of ^3He leads to lower M_ign for <Mdot> >~ 10^-10 Msun/yr, and that for <Mdot> below this the particular author's assumption concerning T_c, which we calculate consistently, is a determining factor. Initial comparisons of our CN ignition masses with measured ejected masses find reasonable agreement and point to ejection of material comparable to that accreted.

Our recent theoretical work (Townsley and Bildsten 2002) on the thermal state of white dwarfs (WDs) in low mass transfer rate binaries allows us to predict the broadband colors of the binary from those of the WD and companion when the disk is dim. The results based on standard CV evolution are presented h ere. These will aid the discovery of such objects in field surveys and proper-motion selected globular cluster surveys with HST; especially for the largely unexplored post period minim um Cataclysmic Variables (CVs) with the lowest accretion rates and degenerate companions. We have also calculated the fraction of time that the WD resides in the ZZ Ceti instability strip thus clarifying that we expect many accreting WDs to exhibit non-radial oscillations. The study of these will provide new insights into the rotational and thermal structure of an actively accreting WD.

We explain the physics of compressional heating of the deep interior of an accreting white dwarf (WD) at accretion rates low enough so that the accumulated hydrogen burns unstably and initiates a classical nova (CN). In this limit, the WD core temperature (T_c) reaches an equilibrium value (T_c,eq) after accreting an amount of mass much less than the WD's mass. Once this equilibrium is reached, the compressional heating from within the envelope exits the surface. This equilibrium yields useful relations between the WD surface temperature, accretion rate and mass that can be employed to measure accretion rates from observed WD effective temperatures, thus testing binary evolution models for cataclysmic variables.

In recent years several Dwarf Novae (DN) systems have been observed in quiescence, when the accretion rate is low and the WD photosphere can be directly detected. The WDs are observed to cool after the DN outburst from high effective temperatures to lower effective temperatures (T_eff) thought to be indicative of the thermal state of the deep interior of the WD. Sion has argued that the most likely energy source for this quiescent luminosity is the gravitational compression of the WD interior, which rejuvenates an otherwise cold WD into a much hotter state. We are undertaking a theoretical study of the compressional heating of WD's, extending down to the very low time averaged accretion rates, <Mdot>~10^{-11}M_sun/yr, applicable to the post-turnaround CV's (the ``TOADS''). Nuclear burning is unstable at these <Mdot>'s, so we have incorporated the recurrent heating and cooling of the WD throughout the classical novae limit cycle. In addition to self-consistently finding the range of T_eff as a function of <Mdot> during the cycle, we also self-consistently find the ignition masses. Comparing these theoretical masses to the observed ejected masses will tell us whether the WD mass in CV's is secularly increasing or decreasing. We close by comparing our results to the accumulated observations of quiescent DN and making predictions for the colors of low <Mdot> CV's in quiescence that are applicable to searches for faint CVs in the field and galactic globular clusters.

Current evolutionary models imply that most Cataclysmic Variables (CVs) have P_orb<2 hours and are Dwarf Nova (DN) systems that are quiescent most of the time. Observations of nearby quiescent DN find that the UV spectrum is dominated by the hot white dwarf (WD), indicating that it provides a significant fraction of the optical light in addition to the quiescent disk and main sequence companion. Hence, identifying a faint, quiescent CV in either the field or a globular cluster (GC) from broadband colors depends on our ability to predict the WD contribution in quiescence. We are undertaking a theoretical study of the compressional heating of WDs, extending down to very low time averaged accretion rates, <Mdot>~10^{-11}M_sun/yr, which allows us to self-consistently find the T_eff} of the WD. We demonstrate here that most of the compressional heating occurs in the freshly accreted envelope and that the WD core temperature reaches a fixed value on a timescale less than typical evolutionary times. Since nuclear burning is unstable at these <Mdot>'s, we have incorporated the recurrent heating and cooling of the WD core throughout the classical novae limit cycle in order to find the T_eff-<Mdot> relations. Comparing to observations of field DN confirms the <Mdot>-P_orb relation of disrupted magnetic braking. We also predict broad-band colors of a quiescent CV as a function of <Mdot> and companion mass and show that this leads to the identification of what may be many CVs in deep HST images of GCs.