We performed 3-D simulations of proton-rich material entrainment into czw-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-AGB star Sakurais object. Observations of the transient nature and anomalous abu
ndance features are available to validate our method and assumptions, with the aim to apply them to very low metallicity stars in the future. We include nuclear energy feedback from H burning and cover the full $4pi$ geometry of the shell. Runs on $768^3$ and $1536^3$ grids agree well with each other and have been followed for $1500mathrm{min}$ and $1200mathrm{min}$. After a $850mathrm{min}$ long quiescent entrainment phase the simulations enter into a global non-spherical oscillation that is launched and sustained by individual ignition events of H-rich fluid pockets. Fast circumferential flows collide at the antipode and cause the formation and localized ignition of the next H-overabundant pocket. The cycle repeats for more than a dozen times while its amplitude decreases. During the global oscillation the entrainment rate increases temporarily by a factor $approx 100$. Entrained entropy quenches convective motions in the upper layer until the burning of entrained H establishes a separate convection zone. The lower-resolution run hints at the possibility that another global oscillation, perhaps even more violent will follow. The location of the H-burning convection zone agrees with a 1-D model in which the mixing efficiency is calibrated to reproduce the light curve. The simulations have been performed at the NSF Blue Waters supercomputer at NCSA.
We present the first 3-dimensional, fully compressible gas-dynamics simulations in $4pi$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed c
onsequences of the H-ingestion flash in post-AGB stars (Sakurais object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that, for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique simulations with $1024^3$ cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the $1536^3$ simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with veryvstiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of $4.38 pm 1.48 times 10^{-13}M_odot mathrm{/s}$.
The R Coronae Borealis (RCB) stars are hydrogen-deficient, variable stars that are most likely the result of He-CO WD mergers. They display extremely low oxygen isotopic ratios, 16O/18O ~ 1 - 10, 12C/13C>=100, and enhancements up to 2.6dex in F and i
n s-process elements from Zn to La, compared to solar. These abundances provide stringent constraints on the physical processes during and after the double-degenerate merger. As shown before O-isotopic ratios observed in RCB stars cannot result from the dynamic double-degenerate merger phase, and we investigate now the role of the long-term 1D spherical post-merger evolution and nucleosynthesis based on realistic hydrodynamic merger progenitor models. We adopt a model for extra envelope mixing to represent processes driven by rotation originating in the dynamical merger. Comprehensive nucleosynthesis post-processing simulations for these stellar evolution models reproduce, for the first time, the full range of the observed abundances for almost all the elements measured in RCB stars: 16O/18O ratios between 9 and 15, C-isotopic ratios above 100, and ~1.4 - 2.35dex F enhancements, along with enrichments in s-process elements. The nucleosynthesis processes in our models constrain the length and temperature in the dynamic merger shell-of-fire feature as well as the envelope mixing in the post-merger phase. s-process elements originate either in the shell-of-fire merger feature or during the post-merger evolution, but the contribution from the AGB progenitors is negligible. The post-merger envelope mixing must eventually cease ~ 10^6yr after the dynamic merger phase, before the star enters the RCB phase.
[Abbreviated] We have investigated the color-magnitude diagram of Omega Centauri and find that the blue main sequence (bMS) can be reproduced only by models that have a of helium abundance in the range Y=0.35-$0.40. To explain the faint subgiant bran
ch of the reddest stars (MS-a/RG-a sequence), isochrones for the observed metallicity ([Fe/H]approx0.7) appear to require both a high age (~13Gyr) and enhanced CNO abundances ([CNO/Fe]approx0.9$). Y~0.35 must also be assumed in order to counteract the effects of high CNO on turnoff colors, and thereby to obtain a good fit to the relatively blue turnoff of this stellar population. This suggest a short chemical evolution period of time (<1Gyr) for Omega Cen. Our intermediate-mass (super-)AGB models are able to reproduce the high helium abundances, along with [N/Fe]~2 and substantial O depletions if uncertainties in the treatment of convection are fully taken into account. These abundance features distinguish the bMS stars from the dominant [Fe/H] $approx1.7$ population. The most massive super-AGB stellar models (M_zams>=6.8M_sun, M_He,core>=1.245M_sun) predict too large N-enhancements, which limits their role in contributing to the extreme populations. We show quantitatively that highly He- and N-enriched AGB ejecta have particularly efficient cooling properties. Based on these results and on the reconstruction of the orbit of Omega Cen with respect to the Milky Way we propose the galactic plane passage gas purging scenario for the chemical evolution of this cluster. Our model addresses the formation and properties of the bMS population (including their central location in the cluster). We follow our model descriptively through four passage events, which could explain not only some key properties of the bMS, but also of the MS-a/RGB-a and the s-enriched stars.
Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convective
ly mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurais object (V4334 Sagittarii). Asplund etal. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurais object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(alpha, n)O16 reaction.
The He-shell flash convection in AGB stars is the site for the high-temperature component of the s-process in low- and intermediate mass giants, driven by the Ne22 neutron source. [...] The upper convection boundary plays a critical role during the H
-ingestion episode that may lead to neutron-bursts in the most metal-poor AGB stars. We address these problems through global 3-dimensional hydrodynamic simulations including the entire spherical He-shell flash convection zone (as oposed to the 3D box-in-a-star simulations). An important aspect of our current effort is to establish the feasibility of our appoach. We explain why we favour the explicit treatment over the anelastic approximation for this problem. The simulations presented in this paper use a Cartesian grid of 512^3 cells and have been run on four 8-core workstations for four days to simulate ~5000s, which corresponds to almost ten convective turn-over times. The convection layer extends radially at the simulated point in the flash evolution over 7 H_p pressure scale-heights and exceeds the size of the underlying core. Convection is dominated by large convective cells that fill more than an entire octant. [...]
The evolution of the first generations of stars at zero or extremly low metallicity, and especially some crucial properties like the primary N14 production, is charactarized by convective-reactive mixing events that are mostly absent from similar evo
lution phases at solar-like metallicity. These episodes occur when unprocessed H-rich material is mixed accross a convective boundary into C12 rich He-burning material, as for example in He-shell flashes of extremely-low metallicity AGB stars. In this paper we describe the astrophysical context of such convective-reactive events, including the difficulty of current one-dimensional stellar evolution models to correctly simulate these evolutionary phases. We then describe the requirements and current state of modeling convective-reactive processes in the first stars environment. We demonstrate some of the new concepts that we are applying to this problem, i.e. the highly accurate PPB advection scheme in the framework of PPM hydrodynamic simulations of mixing accross a very stiff convective boundary. We show initial results of such simulations that address the first non-reactive step of this problem, which is the entrainment of H at the top boundary of the He-shell flash convection zone.
We review the current state of modeling convective mixing in AGB stars. The focus is on results obtained through multi-dimensional hydrodynamic simulations of AGB convection, both in the envelope and the unstable He-shell. Using two different codes a
nd a wide range of resolutions and modeling assumptions we find that mixing across convective boundaries is significant for He-shell flash convection. We present a preliminary quantitative analysis of this convectively induced extra mixing, based on a sub-set of our simulations. Other non-standard mixing will be discussed briefly.
