Deep Mixing of He-3: Reconciling Big Bang and Stellar Nucleosynthesis

Low-mass stars, ~1-2 solar masses, near the Main Sequence are efficient at producing He-3, which they mix into the convective envelope on the giant branch and should distribute into the Galaxy by way of envelope loss. This process is so efficient that it is difficult to reconcile the low observed cosmic abundance of He-3 with the predictions of both stellar and Big Bang nucleosynthesis. In this paper we find, by modeling a red giant with a fully three-dimensional hydrodynamic code and a full nucleosynthetic network, that mixing arises in the supposedly stable and radiative zone between the hydrogen-burning shell and the base of the convective envelope. This mixing is due to Rayleigh-Taylor instability within a zone just above the hydrogen-burning shell, where a nuclear reaction lowers the mean molecular weight slightly. Thus we are able to remove the threat that He-3 production in low-mass stars poses to the Big Bang nucleosynthesis of He-3.Comment: Accepted by Science, and available from Science Express onlin

Compulsory Deep Mixing of 3He and CNO Isotopes in the Envelopes of low-mass Red Giants

Three-dimensional stellar modeling has enabled us to identify a deep-mixing mechanism that must operate in all low mass giants. This mixing process is not optional, and is driven by a molecular weight inversion created by the 3He(3He,2p)4He reaction. In this paper we characterize the behavior of this mixing, and study its impact on the envelope abundances. It not only eliminates the problem of 3He overproduction, reconciling stellar and big bang nucleosynthesis with observations, but solves the discrepancy between observed and calculated CNO isotope ratios in low mass giants, a problem of more than 3 decades' standing. This mixing mechanism, which we call `$\delta\mu$-mixing', operates rapidly (relative to the nuclear timescale of overall evolution, ~ 10^8 yrs) once the hydrogen burning shell approaches the material homogenized by the surface convection zone. In agreement with observations, Pop I stars between 0.8 and 2.0\Msun develop 12C/13C ratios of 14.5 +/- 1.5, while Pop II stars process the carbon to ratios of 4.0 +/- 0.5. In stars less than 1.25\Msun, this mechanism also destroys 90% to 95% of the 3He produced on the main sequence.Comment: Final accepted version (submitted to Astrophys J in Jan 2007...

Low and intermediate-mass close binary evolution and the initial - final mass relation

Using Eggleton's stellar evolution code, we carry out 150 runs of Pop I binary evolution calculations, with the initial primary mass between 1 and 8 solar masses the initial mass ratio between 1.1 and 4, and the onset of Roche lobe overflow (RLOF) at an early, middle, or late Hertzsprung-gap stage. We assume that RLOF is conservative in the calculations, and find that the remnant mass of the primary may change by more than 40 per cent over the range of initial mass ratio or orbital period, for a given primary mass. This is contrary to the often-held belief that the remnant mass depends only on the progenitor mass if mass transfer begins in the Hertzsprung gap. We fit a formula, with an error less than 3.6 per cent, for the remnant (white dwarf) mass as a function of the initial mass of the primary, the initial mass ratio, and the radius of the primary at the onset of RLOF. We also find that a carbon-oxygen white dwarf with mass as low as 0.33 solar masses may be formed if the primary's initial mass is around 2.5 solar masses.Comment: 7 pages for main text, 11 pages for appendix (table A1), 12 figure

The Destruction of 3He by Rayleigh-Taylor Instability on the First Giant Branch

Low-mass stars, ~1-2 solar masses, near the Main Sequence are efficient at producing 3He, which they mix into the convective envelope on the giant branch and distribute into the Galaxy by way of envelope loss. This process is so efficient that it is difficult to reconcile the observed cosmic abundance of 3He with the predictions of Big Bang nucleosynthesis. In this paper we find, by modeling a red giant with a fully three-dimensional hydrodynamic code and a full nucleosynthetic network, that mixing arises in the supposedly stable and radiative zone between the hydrogen-burning shell and the base of the convective envelope. This mixing is due to Rayleigh-Taylor instability within a zone just above the hydrogen-burning shell. In this zone the burning of the 3He left behind by the retreating convective envelope is predominantly by the reaction 3He + 3He -> 4He + 2p, a reaction which, untypically for stellar nuclear reactions, {\it lowers} the mean molecular weight, leading to a local minimum. This local minimum leads to Rayleigh-Taylor instability, and turbulent motion is generated which will continue ultimately up into the normal convective envelope. Consequently material from the envelope is dragged down sufficiently close to the burning shell that the 3He in it is progressively destroyed. Thus we are able to remove the threat that 3He production in low-mass stars poses to the Big Bang nucleosynthesis of 3He. Some slow mixing mechanism has long been suspected, that connects the convective envelope of a red giant to the burning shell. It appears to be necessary to account for progressive changes in the 12C/13C and 14N/12C ratios on the First Giant Branch. We suggest that these phenomena are also due to the Rayleigh-Taylor-unstable character of the 3He-burning region.Comment: Paper presented at IAU-GA, Prague, August 200

A critical test of stellar evolution and convective core overshooting by means of zeta-Aur systems

No description supplie

Evolution in Binary and Triple Stars, with an application to SS Lac

We present equations governing the way in which both the orbit and the intrinsic spins of stars in a close binary should evolve subject to a number of perturbing forces, including the effect of a third body in a possibly inclined wider orbit. We illustrate the solutions in some binary-star and triple-star situations: tidal friction in a wide but eccentric orbit of a radio pulsar about a B star, the Darwin and eccentricity instabilities in a more massive but shorter-period massive X-ray binary, and the interaction of tidal friction with Kozai cycles in a triple such as Algol (beta-Per), at an early stage in that star's life when all 3 components were ZAMS stars. We also attempt to model in some detail the interesting triple system SS Lac, which stopped eclipsing in about 1950. We find that our model of SS Lac is quite constrained by the relatively good observational data of this system, and leads to a specific inclination (29 deg) of the outer orbit relative to the inner orbit at epoch zero (1912). Although the intrinsic spins of the stars have little effect on the orbit, the converse is not true: the spin axes can vary their orientation relative to the close binary by up to 120 deg on a timescale of about a century.Comment: 30 pages, 6 figure

N-body simulations of stars escaping from the Orion nebula

We study the dynamical interaction in which the two single runaway stars AE Aurigae and mu Columbae and the binary iota Orionis acquired their unusually high space velocity. The two single runaways move in almost opposite directions with a velocity greater than 100 km/s away from the Trapezium cluster. The star iota Ori is an eccentric (e=0.8) binary moving with a velocity of about 10 km/s at almost right angles with respect to the two single stars. The kinematic properties of the system suggest that a strong dynamical encounter occurred in the Trapezium cluster about 2.5 Myr ago. Curiously enough, the two binary components have similar spectral type but very different masses, indicating that their ages must be quite different. This observation leads to the hypothesis that an exchange interaction occurred in which an older star was swapped into the original iota Orionis binary. We test this hypothesis by a combination of numerical and theoretical techniques, using N-body simulations to constrain the dynamical encounter, binary evolution calculations to constrain the high orbital eccentricity of iota Orionis and stellar evolution calculations to constrain the age discrepancy of the two binary components. We find that an encounter between two low eccentricity (0.4<e<0.6) binaries with comparable binding energy, leading to an exchange and the ionization of the wider binary, provides a reasonable solution to this problem.Comment: 14 pages, 13 figures, 5 tables, accepted for publication in MNRA

Envelope Ejection: an Alternative Process for some Early Case B Binaries

We discuss the evolution of binaries with moderately high masses (about 10 - 30 solar masses), and with periods of about 3 - 300d, corresponding mostly to early Case B. These are usually thought to evolve either by reasonably conservative Roche-lobe overflow, if the initial mass ratio is fairly mild, or else by highly non-conservative common-envelope evolution, with spiral-in to short periods (hours, typically), if the initial mass ratio is rather extreme. We discuss here a handful of binaries from part of this period range (about 50 - 250d), which appear to have followed a different path: we argue that they must have lost a large proportion of initial mass (about 70 - 80%), but without shortening their periods at all. We suggest that their behaviour may be due to the fact that stars of such masses, when evolved also to rather large radii, are not far from the Humphreys-Davidson limit where single stars lose their envelopes spontaneously in P Cygni winds, and so have envelopes which are only lightly bound to the core. These envelopes therefore may be relatively easily dissipated by the perturbing effect of a companion. In addition, some or all of the stars considered here may have been close to the Cepheid instability strip when they filled their Roche lobes. One or other, or both, of high luminosity and Cepheid instability, in combination with an appropriately close binary companion, may be implicated