White Dwarf Colors in Low Accretion Rate Binaries

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.Comment: 2 pages, 1 figure, uses CRCKAPB.sty (included); to appear in the `NATO Science Series II - Mathematics, Physics and Chemistry', Kluwer Academic Publisher

Compressional Heating of Accreting White Dwarfs in CV's

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, ~10^{-11}M_sun/yr, applicable to the post-turnaround CV's (the ``TOADS''). Nuclear burning is unstable at these '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 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 CV's in quiescence that are applicable to searches for faint CVs in the field and galactic globular clusters.Comment: 10 pages, 5 figures; uses newpasp.sty (included); to appear in `The Physics of Cataclysmic Variables and Related Objects', ASP Conf. Ser., eds. B.T. Gaensicke, K. Beuermann, K. Reinsc

Seismology of Rapidly Rotating Accreting White Dwarfs

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.Comment: 6 pages, 4 figures, proceedings from 2010 conference "The Physics of Accreting Compact Binaries" Kyoto, Japan. Submitted to proceedings in 201

Pulsational Instabilities in Accreting White Dwarfs

(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.Comment: submitted to ApJL. 3 figure

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

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.Comment: 14 pages, 10 figures. Submitted to ApJ. Comments welcom

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

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.Comment: 20 pages, 12 figures, accepted to the Astrophysical Journa

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

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.Comment: 14 pages, 19 figures, submitted to the Astrophysical Journa

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

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.Comment: Submitted to Ap

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

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.Comment: 28 pages, 16 figures, Accepted to the Astrophysical Journal Supplementa

Thermonuclear (Type Ia) Supernovae and Progenitor Evolution

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.Comment: 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, 201