Chemical equilibrium and stable stratification of a multi-component fluid: thermodynamics and application to neutron stars

A general thermodynamic argument shows that multi-component matter in full chemical equilibrium, with uniform entropy per baryon, is generally stably stratified. This is particularly relevant for neutron stars, in which the effects of entropy are negligible compared to those of the equilibrium composition gradient established by weak interactions. It can therefore be asserted that, regardless of the uncertainties in the equation of state of dense matter, neutron stars are stably stratified. This has important, previously discussed consequences for their oscillation modes, magnetic field evolution, and internal angular momentum transport.Comment: AASTeX, 8 pages, including 1 PS figure. Accepted for publication in The Astrophysical Journa

Active Carbon and Oxygen Shell Burning Hydrodynamics

We have simulated 2.5$\times10^3$ s of the late evolution of a $23 \rm M_\odot$ star with full hydrodynamic behavior. We present the first simulations of a multiple-shell burning epoch, including the concurrent evolution and interaction of an oxygen and carbon burning shell. In addition, we have evolved a 3D model of the oxygen burning shell to sufficiently long times (300 s) to begin to assess the adequacy of the 2D approximation. We summarize striking new results: (1) strong interactions occur between active carbon and oxygen burning shells, (2) hydrodynamic wave motions in nonconvective regions, generated at the convective-radiative boundaries, are energetically important in both 2D and 3D with important consequences for compositional mixing, and (3) a spectrum of mixed p- and g-modes are unambiguously identified with corresponding adiabatic waves in these computational domains. We find that 2D convective motions are exaggerated relative to 3D because of vortex instability in 3D. We discuss the implications for supernova progenitor evolution and symmetry breaking in core collapse.Comment: 5 pages, 4 figures in emulateapj format. Accepted for publication in ApJ Letters. High resolution figure version available at http://spinach.as.arizona.ed

The Impact of Hydrodynamic Mixing on Supernova Progenitors

Recent multidimensional hydrodynamic simulations have demonstrated the importance of hydrodynamic motions in the convective boundary and radiative regions of stars to transport of energy, momentum, and composition. The impact of these processes increases with stellar mass. Stellar models which approximate this physics have been tested on several classes of observational problems. In this paper we examine the implications of the improved treatment on supernova progenitors. The improved models predict substantially different interior structures. We present pre-supernova conditions and simple explosion calculations from stellar models with and without the improved mixing treatment at 23 solar masses. The results differ substantially.Comment: 12 pages, 2 figures, accepted for publication in the Astrophysical Journal Letter

Theoretical Studies of Accretion of Matter onto White Dwarfs and the Single Degenerate Scenario for Supernovae of Type Ia

We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown.Comment: 25 Pages, Bulletin of the Astronomical Society of India (BASI) in pres

Evidence for a Mid-Atomic-Number Atmosphere in the Neutron Star 1E1207.4-5209

Recently Sanwal et al. (2002) reported the first clear detection of absorption features in an isolated neutron star, 1E1207.4-5209. Remarkably their spectral modeling demonstrates that the atmosphere cannot be Hydrogen. They speculated that the neutron star atmosphere is indicative of ionized Helium in an ultra-strong (~1.5x10^{14} G) magnetic field. We have applied our recently developed atomic model (Mori & Hailey 2002) for strongly-magnetized neutron star atmospheres to this problem. We find that this model, along with some simp le atomic physics arguments, severely constrains the possible composition of the atmosphere. In particular we find that the absorption features are naturally associated with He-like Oxygen or Neon in a magnetic field of ~10^{12} G, comparable to the magnetic field derived from the spin parameters of the neutron star. This interpretation is consistent with the relative line strengths and widths and is robust. Our model predicts possible substructure in the spectral features, which has now been reported by XMM-Newton (Mereghetti et al. 2002). However we show the Mereghetti et al. claim that the atmosphere is Iron or some comparable high-Z element at ~ 10^{12} G is easily ruled out by the Chandra and XMM-Newton data.Comment: 5 pages, AASTeX, Revised version. Accepted for publication in ApJ Letter

Modeling core collapse supernovae in 2 and 3 dimensions with spectral neutrino transport

The overwhelming evidence that the core collapse supernova mechanism is inherently multidimensional, the complexity of the physical processes involved, and the increasing evidence from simulations that the explosion is marginal presents great computational challenges for the realistic modeling of this event, particularly in 3 spatial dimensions. We have developed a code which is scalable to computations in 3 dimensions which couples PPM Lagrangian with remap hydrodynamics [1], multigroup, flux-limited diffusion neutrino transport [2], with many improvements), and a nuclear network [3]. The neutrino transport is performed in a ray-by-ray plus approximation wherein all the lateral effects of neutrinos are included (e.g., pressure, velocity corrections, advection) except the transport. A moving radial grid option permits the evolution to be carried out from initial core collapse with only modest demands on the number of radial zones. The inner part of the core is evolved after collapse along with the rest of the core and mantle by subcycling the lateral evolution near the center as demanded by the small Courant times. We present results of 2-D simulations of a symmetric and an asymmetric collapse of both a 15 and an 11 M progenitor. In each of these simulations we have discovered that once the oxygen rich material reaches the shock there is a synergistic interplay between the reduced ram pressure, the energy released by the burning of the shock heated oxygen rich material, and the neutrino energy deposition which leads to a revival of the shock and an explosion.Comment: 10 pages, 3 figure

Observational Tests and Predictive Stellar Evolution II: Non-standard Models

We examine contributions of second order physical processes to results of stellar evolution calculations amenable to direct observational testing. In the first paper in the series (Young et al. 2001) we established baseline results using only physics which are common to modern stellar evolution codes. In the current paper we establish how much of the discrepancy between observations and baseline models is due to particular elements of new physics. We then consider the impact of the observational uncertainties on the maximum predictive accuracy achievable by a stellar evolution code. The sun is an optimal case because of the precise and abundant observations and the relative simplicity of the underlying stellar physics. The Standard Model is capable of matching the structure of the sun as determined by helioseismology and gross surface observables to better than a percent. Given an initial mass and surface composition within the observational errors, and no additional constraints for which the models can be optimized, it is not possible to predict the sun's current state to better than ~7%. Convectively induced mixing in radiative regions, seen in multidimensional hydrodynamic simulations, dramatically improves the predictions for radii, luminosity, and apsidal motions of eclipsing binaries while simultaneously maintaining consistency with observed light element depletion and turnoff ages in young clusters (Young et al. 2003). Systematic errors in core size for models of massive binaries disappear with more complete mixing physics, and acceptable fits are achieved for all of the binaries without calibration of free parameters. The lack of accurate abundance determinations for binaries is now the main obstacle to improving stellar models using this type of test.Comment: 33 pages, 8 figures, accepted for publication in the Astrophysical Journa

Fastener starter tool

A fastener starter tool includes a number of spring retention fingers for retaining a small part, or combination of parts. The tool has an inner housing, which holds the spring retention fingers, a hand grip, and an outer housing configured to slide over the inner housing and the spring retention fingers toward and away from the hand grip, exposing and opening, or respectively, covering and closing, the spring retention fingers. By sliding the outer housing toward (away from) the hand grip, a part can be released from (retained by) the tool. The tool may include replaceable inserts, for retaining parts, such as screws, and configured to limit the torque applied to the part, to prevent cross threading. The inner housing has means to transfer torque from the hand grip to the insert. The tool may include replaceable bits, the inner housing having means for transferring torque to the replaceable bit

3D hydrodynamic simulations of carbon burning in massive stars

We present the first detailed 3D hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon-burning shell within a 15 M⊙ 1D stellar evolution model. We consider models with 1283, 2563, 5123, and 10243 zones. The turbulent flow properties of these carbon-burning simulations are very similar to the oxygen-burning case. We performed a mean field analysis of the kinetic energy budgets within the Reynolds-averaged Navier–Stokes framework. For the upper convective boundary region, we find that the numerical dissipation is insensitive to resolution for linear mesh resolutions above 512 grid points. For the stiffer, more stratified lower boundary, our highest resolution model still shows signs of decreasing sub-grid dissipation suggesting it is not yet numerically converged. We find that the widths of the upper and lower boundaries are roughly 30 per cent and 10 per cent of the local pressure scaleheights, respectively. The shape of the boundaries is significantly different from those used in stellar evolution models. As in past oxygen-shell-burning simulations, we observe entrainment at both boundaries in our carbon-shell-burning simulations. In the large Péclet number regime found in the advanced phases, the entrainment rate is roughly inversely proportional to the bulk Richardson number, RiB (∝RiB−α, 0.5 ≲ α ≲ 1.0). We thus suggest the use of RiB as a means to take into account the results of 3D hydrodynamics simulations in new 1D prescriptions of convective boundary mixing