J. W. Truran

CNO ABUNDANCES AND HYDRODYNAMIC MODELS OF THE NOVA OUTBURST.

We have used a fully implicit, Lagrangian, hydrodynamic computer code incorporating a nuclear reaction network to follow thermonuclear runaways in the hydrogen-rich envelopes of white dwarfs in order to produce a nova outburst. Because of the short time-scales and the high nuclear burning rates produced in our models, the nuclear reactions are far out of equilibrium and the beta-plus unstable nuclei become the most abundant nuclei in the envelope except for hydrogen and helium. Our models have ejected 1.00017 solar mass with kinetic energies of 8 times 10 to the 44-th power ergs, a value that agrees quite closely with the observed values for novae.

Three-Dimensional Simulations of the Deflagration Phase of the Gravitationally Confined Detonation Model of Type Ia Supernovae

We report the results of a series of three-dimensional (3D) 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 the breakout on the stellar surface. We find that detonation conditions are robustly reached in our 3D 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 56Ni 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.

Charged-particle and neutron-capture processes in the high-entropy wind of core-collapse supernovae

The astrophysical site of the r-process is still uncertain, and a full exploration of the systematics of this process in terms of its dependence on nuclear properties from stability to the neutron drip-line within realistic stellar environments has still to be undertaken. Sufficiently high neutron-to-seed ratios can only be obtained either in very neutron-rich low-entropy environments or moderately neutron-rich high-entropy environments, related to neutron star mergers (or jets of neutron star matter) and the high-entropy wind of core-collapse supernova explosions. As chemical evolution models seem to disfavor neutron star mergers, we focus here on high-entropy environments characterized by entropy S, electron abundance Y-e, and expansion velocity V-exp. We investigate the termination point of charged-particle reactions, and we define a maximum entropy S-final for a given V-exp and Y-e, beyond which the seed production of heavy elements fails due to the very small matter density. We then investigate whether an r-process subsequent to the charged-particle freeze-out can in principle be understood on the basis of the classical approach, which assumes a chemical equilibrium between neutron captures and photodisintegrations, possibly followed by a beta-flow equilibrium. In particular, we illustrate how long such a chemical equilibrium approximation holds, how the freeze-out from such conditions affects the abundance pattern, and which role the late capture of neutrons originating from beta-delayed neutron emission can play. Furthermore, we analyze the impact of nuclear properties from different theoretical mass models on the final abundances after these late freeze-out phases and beta-decays back to stability. As only a superposition of astrophysical conditions can provide a good fit to the solar r-abundances, the question remains how such superpositions are attained, resulting in the apparently robust r-process pattern observed in low metallicity stars.

Morphology of Rising Hydrodynamic and Magnetohydrodynamic Bubbles from Numerical Simulations

Recent Chandra and XMM-Newton observations of galaxy cluster cooling flows have revealed X-ray emission voids of up to 30 kpc in size that have been identified with buoyant, magnetized bubbles. Motivated by these observations, we have investigated the behavior of rising bubbles in stratified atmospheres using the FLASH adaptive-mesh simulation code. We present results from two-dimensional simulations with and without the effects of magnetic fields and with varying bubble sizes and background stratifications. We find purely hydrodynamic bubbles to be unstable; a dynamically important magnetic field is required to maintain a bubbles integrity. This suggests that, even absent thermal conduction, for bubbles to be persistent enough to be regularly observed, they must be supported in large part by magnetic fields. Thermal conduction unmitigated by magnetic fields can dissipate the bubbles even faster. We also observe that the bubbles leave a tail as they rise; the structure of these tails can indicate the history of the dynamics of the rising bubble.

From First Stars to the Spite Plateau: A Possible Reconciliation of Halo Stars Observations with Predictions from Big Bang Nucleosynthesis

Since the pioneering observations of Spite & Spite in 1982, the constant lithium abundance of metal-poor ([Fe/H] < -1.3) halo stars near the turnoff has been attributed to a cosmological origin. Closer analysis, however, revealed that the observed abundance lies at Δ7Li ~ 0.4 dex below the predictions of big bang nucleosynthesis (BBN). The measurements of deuterium abundances along the lines of sight toward quasars, and the recent data from the Wilkinson Microwave Anisotropy Probe (WMAP), have independently confirmed this gap. We suggest here that part of the discrepancy (from 0.2 to 0.3 dex) is explained by a first generation of stars that efficiently depleted lithium. Assuming that the models for lithium evolution in halo turnoff stars, as well as the Δ7Li, estimates are correct, we infer that between one-third and one-half of the baryonic matter of the early halo (i.e., ~109 M☉) was processed through Population III stars. This new paradigm proposes a very economical solution to the lingering difficulty of understanding the properties of the Spite plateau and its lack of star-to-star scatter down to [Fe/H] = -2.5. It is moreover in agreement both with the absence of lithium in the most iron-poor turnoff star currently known (HE 1327-2326) and also with new trends of the plateau suggesting its low-metallicity edge may be reached around [Fe/H] = -2.5. We discuss the role of turbulent mixing associated with enhanced supernovae explosions in the early interstellar medium in this picture. We suggest how it may explain the small scatter and also other recent observational features of the lithium plateau. Finally, we show that other chemical properties of the extremely metal-poor stars (such as carbon enrichment) are also in agreement with significant Population III processing in the halo, provided these models include mass loss and rotationally induced mixing.

Flame Evolution During Type Ia Supernovae and the Deflagration Phase in the Gravitationally Confined Detonation Scenario

We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova and apply it in a study of the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified three-stage burning model and a nonstatic 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 neither our ADR nor our energy release methods 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.

A Chandra low energy transmission grating spectrometer observation of V4743 Sagittarii : a supersoft X-ray source and a violently variable light curve

V4743 Sagittarii (Nova Sgr 2002 No. 3) was discovered on 2002 September 20. We obtained a 5 ks ACIS-S spectrum in 2002 November and found that the nova was faint in X-rays. We then obtained a 25 ks Chandra Low Energy Transmission Grating Spectrometer (LETGS) observation on 2003 March 19. By this time, it had evolved into the supersoft X-ray phase exhibiting a continuous spectrum with deep absorption features. The light curve from the observation showed large-amplitude oscillations with a period of 1325 s (22 minutes) followed by a decline in the total count rate after ~13 ks of observations. The count rate dropped from ~40 counts s-1 to practically zero within ~6 ks and stayed low for the rest of the observation (~6 ks). The spectral hardness ratio changed from maxima to minima in correlation with the oscillations and then became significantly softer during the decay. Strong H-like and He-like lines of oxygen, nitrogen, and carbon were found in absorption during the bright phase, indicating temperatures between 1 and 2 MK, but they were shifted in wavelength corresponding to a Doppler velocity of -2400 km s-1. The spectrum obtained after the decline in count rate showed emission lines of C VI, N VI, and N VII, suggesting that we were seeing expanding gas ejected during the outburst, probably originating from CNO-cycled material. An XMM-Newton Target of Opportunity observation, obtained on 2002 April 4 and a later LETGS observation from 2003 July 18 also showed oscillations, but with smaller amplitudes.

Capturing the Fire: Flame Energetics and Neutronization for Type Ia Supernova Simulations

We develop and calibrate a realistic model flame for hydrodynamic 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 postflame 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.

Nucleosynthesis Modes in The High-Entropy Wind of Type II Supernovae: Comparison of Calculations With Halo-Star Observations

While the high-entropy wind (HEW) of Type II supernovae remains one of the more promising sites for the rapid neutron-capture (r-) process, hydrodynamic simulations have yet to reproduce the astrophysical conditions under which the latter occurs. We have performed large-scale network calculations within an extended parameter range of the HEW, seeking to identify or to constrain the necessary conditions for a full reproduction of all r-process residuals N r,☉ = N ☉–N s,☉ by comparing the results with recent astronomical observations. A superposition of weighted entropy trajectories results in an excellent reproduction of the overall N r,☉ pattern beyond Sn. For the lighter elements, from the Fe group via Sr-Y-Zr to Ag, our HEW calculations indicate a transition from the need for clearly different sources (conditions/sites) to a possible co-production with r-process elements, provided a range of entropies are contributing. This explains recent halo-star observations of a clear noncorrelation of Zn and Ge and a weak correlation of Sr-Zr with heavier r-process elements. Moreover, new observational data on Ru and Pd also seem to confirm a partial correlation with Sr as well as the main r-process elements (e.g., Eu).

On Li-7 production in nova explosions

Calculations of Li-7 production occurring as a concomitant of thermonuclear runaways in hydrogen envelopes of white dwarfs are reported. It is found that sufficient Li-7 can be produced in models displaying fast-nova-like features to suggest that the corresponding objects represent significant contributors to the enrichment of galactic matter. The sensitivities of these results to various assumptions and uncertainties are discussed.