Evolution and nucleosynthesis of asymptotic giant branch stellar models of low metallicity

We present stellar evolutionary tracks and nucleosynthetic predictions for a grid of stellar models of low- and intermediate-mass asymptotic giant branch (AGB) stars at $Z=0.001$ ([Fe/H]$=-1.2$). The models cover an initial mass range from 1 M$_{\odot}$ to 7 M$_{\odot}$. Final surface abundances and stellar yields are calculated for all elements from hydrogen to bismuth as well as isotopes up to the iron group. We present the first study of neutron-capture nucleosynthesis in intermediate-mass AGB models, including a super-AGB model, of [Fe/H] = $-1.2$. We examine in detail a low-mass AGB model of 2 M$_{\odot}$ where the $^{13}$C($\alpha$,$n$)$^{16}$O reaction is the main source of neutrons. We also examine an intermediate-mass AGB model of 5 M$_{\odot}$ where intershell temperatures are high enough to activate the $^{22}$Ne neutron source, which produces high neutron densities up to $\sim 10^{14}$ n cm$^{-3}$. Hot bottom burning is activated in models with $M \geq 3$ M$_{\odot}$. With the 3 M$_{\odot}$ model we investigate the effect of varying the extent in mass of the region where protons are mixed from the envelope into the intershell at the deepest extent of each third dredge-up. We compare the results of the low-mass models to three post-AGB stars with a metallicity of [Fe/H] $\simeq -1.2$. The composition is a good match to the predicted neutron-capture abundances except for Pb and we confirm that the observed Pb abundances are lower than what is calculated by AGB models.Comment: 23 pages, 18 figures, 13 tables, accepted for publication in Ap

Carbon-enhanced metal-poor stars: a window on AGB nucleosynthesis and binary evolution. II. Statistical analysis of a sample of 67 CEMP-ss stars

Many observed CEMP stars are found in binary systems and show enhanced abundances of $s$-elements. The origin of the chemical abundances of these CEMP-$s$ stars is believed to be accretion in the past of enriched material from a primary star in the AGB phase. We investigate the mechanism of mass transfer and the process of nucleosynthesis in low-metallicity AGB stars by modelling the binary systems in which the observed CEMP-$s$ stars were formed. For this purpose we compare a sample of $67$ CEMP-$s$ stars with a grid of binary stars generated by our binary evolution and nucleosynthesis model. We classify our sample CEMP-$s$ stars in three groups based on the observed abundance of europium. In CEMP$-s/r$ stars the europium-to-iron ratio is more than ten times higher than in the Sun, whereas it is lower than this threshold in CEMP$-s/nr$ stars. No measurement of europium is currently available for CEMP-$s/ur$ stars. On average our models reproduce well the abundances observed in CEMP-$s/nr$ stars, whereas in CEMP-$s/r$ stars and CEMP-$s/ur$ stars the abundances of the light-$s$ elements are systematically overpredicted by our models and in CEMP-$s/r$ stars the abundances of the heavy-$s$ elements are underestimated. In all stars our modelled abundances of sodium overestimate the observations. This discrepancy is reduced only in models that underestimate the abundances of most of the $s$-elements. Furthermore, the abundance of lead is underpredicted in most of our model stars. These results point to the limitations of our AGB nucleosynthesis model, particularly in the predictions of the element-to-element ratios. Finally, in our models CEMP-$s$ stars are typically formed in wide systems with periods above 10000 days, while most of the observed CEMP-$s$ stars are found in relatively close orbits with periods below 5000 days.Comment: 23 pages, 8 figures, accepted for publication on Astronomy & Astrophysic

Reaction rate uncertainties and the operation of the NeNa and MgAl chains during HBB in intermediate-mass AGB stars

We test the effect of proton-capture reaction rate uncertainties on the abundances of the Ne, Na, Mg and Al isotopes processed by the NeNa and MgAl chains during hot bottom burning (HBB) in asymptotic giant branch (AGB) stars of intermediate mass between 4 and 6 solar masses and metallicities between Z=0.0001 and 0.02. We provide uncertainty ranges for the AGB stellar yields, for inclusion in galactic chemical evolution models, and indicate which reaction rates are most important and should be better determined. We use a fast synthetic algorithm based on detailed AGB models. We run a large number of stellar models, varying one reaction per time for a very fine grid of values, as well as all reactions simultaneously. We show that there are uncertainties in the yields of all the Ne, Na, Mg and Al isotopes due to uncertain proton-capture reaction rates. The most uncertain yields are those of 26Al and 23Na (variations of two orders of magnitude), 24Mg and 27Al (variations of more than one order of magnitude), 20Ne and 22Ne (variations between factors 2 and 7). In order to obtain more reliable Ne, Na, Mg and Al yields from IM-AGB stars the rates that require more accurate determination are: 22Ne(p,g)23Na, 23Na(p,g)24Mg, 25Mg(p,g)26Al, 26Mg(p,g)27Al and 26Al(p,g)27Si. Detailed galactic chemical evolution models should be constructed to address the impact of our uncertainty ranges on the observational constraints related to HBB nucleosynthesis, such as globular cluster chemical anomalies.Comment: accepted for publication on Astronomy & Astrophysic

Carbon-enhanced metal-poor stars: a window on AGB nucleosynthesis and binary evolution. I. Detailed analysis of 15 binary stars with known orbital periods

AGB stars are responsible for producing a variety of elements, including carbon, nitrogen, and the heavy elements produced in the slow neutron-capture process ($s$-elements). There are many uncertainties involved in modelling the evolution and nucleosynthesis of AGB stars, and this is especially the case at low metallicity, where most of the stars with high enough masses to enter the AGB have evolved to become white dwarfs and can no longer be observed. The stellar population in the Galactic halo is of low mass ($\lesssim 0.85M_{\odot}$) and only a few observed stars have evolved beyond the first giant branch. However, we have evidence that low-metallicity AGB stars in binary systems have interacted with their low-mass secondary companions in the past. The aim of this work is to investigate AGB nucleosynthesis at low metallicity by studying the surface abundances of chemically peculiar very metal-poor stars of the halo observed in binary systems. To this end we select a sample of 15 carbon- and $s$-element-enhanced metal-poor (CEMP-$s$) halo stars that are found in binary systems with measured orbital periods. With our model of binary evolution and AGB nucleosynthesis, we determine the binary configuration that best reproduces, at the same time, the observed orbital period and surface abundances of each star of the sample. The observed periods provide tight constraints on our model of wind mass transfer in binary stars, while the comparison with the observed abundances tests our model of AGB nucleosynthesis.Comment: 18 pages, 20 figures, accepted for publication on A&

Modelling the observed properties of carbon-enhanced metal-poor stars using binary population synthesis

The stellar population in the Galactic halo is characterised by a large fraction of CEMP stars. Most CEMP stars are enriched in $s$-elements (CEMP-$s$ stars), and some of these are also enriched in $r$-elements (CEMP-$s/r$ stars). One formation scenario proposed for CEMP stars invokes wind mass transfer in the past from a TP-AGB primary star to a less massive companion star which is presently observed. We generate low-metallicity populations of binary stars to reproduce the observed CEMP-star fraction. In addition, we aim to constrain our wind mass-transfer model and investigate under which conditions our synthetic populations reproduce observed abundance distributions. We compare the CEMP fractions and the abundance distributions determined from our synthetic populations with observations. Several physical parameters of the binary stellar population of the halo are uncertain, e.g. the initial mass function, the mass-ratio and orbital-period distributions, and the binary fraction. We vary the assumptions in our model about these parameters, as well as the wind mass-transfer process, and study the consequent variations of our synthetic CEMP population. The CEMP fractions calculated in our synthetic populations vary between 7% and 17%, a range consistent with the CEMP fractions among very metal-poor stars recently derived from the SDSS/SEGUE data sample. The results of our comparison between the modelled and observed abundance distributions are different for CEMP-$s/r$ stars and for CEMP-$s$ stars. For the latter, our simulations qualitatively reproduce the observed distributions of C, Na, Sr, Ba, Eu, and Pb. Contrarily, for CEMP-$s/r$ stars our model cannot reproduce the large abundances of neutron-rich elements such as Ba, Eu, and Pb. This result is consistent with previous studies, and suggests that CEMP-$s/r$ stars experienced a different nucleosynthesis history to CEMP-$s$ stars.Comment: 17 pages, 11 figures, accepted for publication on Astronomy and Astrophysic

Evolution and nucleosynthesis of helium-rich asymptotic giant branch models

There is now strong evidence that some stars have been born with He mass fractions as high as $Y \approx 0.40$ (e.g., in $\omega$ Centauri). However, the advanced evolution, chemical yields, and final fates of He-rich stars are largely unexplored. We investigate the consequences of He-enhancement on the evolution and nucleosynthesis of intermediate-mass asymptotic giant branch (AGB) models of 3, 4, 5, and 6 M$_\odot$ with a metallicity of $Z = 0.0006$ ([Fe/H] $\approx -1.4$). We compare models with He-enhanced compositions ($Y=0.30, 0.35, 0.40$) to those with primordial He ($Y=0.24$). We find that the minimum initial mass for C burning and super-AGB stars with CO(Ne) or ONe cores decreases from above our highest mass of 6 M$_\odot$ to $\sim$ 4-5 M$_\odot$ with $Y=0.40$. We also model the production of trans-Fe elements via the slow neutron-capture process (s-process). He-enhancement substantially reduces the third dredge-up efficiency and the stellar yields of s-process elements (e.g., 90% less Ba for 6 M$_\odot$, $Y=0.40$). An exception occurs for 3 M$_\odot$, where the near-doubling in the number of thermal pulses with $Y=0.40$ leads to $\sim$ 50% higher yields of Ba-peak elements and Pb if the $^{13}$C neutron source is included. However, the thinner intershell and increased temperatures at the base of the convective envelope with $Y=0.40$ probably inhibit the $^{13}$C neutron source at this mass. Future chemical evolution models with our yields might explain the evolution of s-process elements among He-rich stars in $\omega$ Centauri.Comment: 21 pages, 16 figures, accepted for publication by MNRAS. Stellar yields included as online data table

On the asymptotic giant branch star origin of peculiar spinel grain OC2

Microscopic presolar grains extracted from primitive meteorites have extremely anomalous isotopic compositions revealing the stellar origin of these grains. The composition of presolar spinel grain OC2 is different from that of all other presolar spinel grains. Large excesses of the heavy Mg isotopes are present and thus an origin from an intermediate-mass (IM) asymptotic giant branch (AGB) star was previously proposed for this grain. We discuss the isotopic compositions of presolar spinel grain OC2 and compare them to theoretical predictions. We show that the isotopic composition of O, Mg and Al in OC2 could be the signature of an AGB star of IM and metallicity close to solar experiencing hot bottom burning, or of an AGB star of low mass (LM) and low metallicity suffering very efficient cool bottom processing. Large measurement uncertainty in the Fe isotopic composition prevents us from discriminating which model better represents the parent star of OC2. However, the Cr isotopic composition of the grain favors an origin in an IM-AGB star of metallicity close to solar. Our IM-AGB models produce a self-consistent solution to match the composition of OC2 within the uncertainties related to reaction rates. Within this solution we predict that the 16O(p,g)17F and the 17O(p,a)14N reaction rates should be close to their lower and upper limits, respectively. By finding more grains like OC2 and by precisely measuring their Fe and Cr isotopic compositions, it may be possible in the future to derive constraints on massive AGB models from the study of presolar grains.Comment: 10 pages, 8 figures, accepted for publication on Astronomy & Astrophysic

The Chemical Evolution of Helium in Globular Clusters: Implications for the Self-Pollution Scenario

We investigate the suggestion that there are stellar populations in some globular clusters with enhanced helium (Y from 0.28 to 0.40) compared to the primordial value. We assume that a previous generation of massive Asymptotic Giant Branch (AGB) stars have polluted the cluster. Two independent sets of AGB yields are used to follow the evolution of helium and CNO using a Salpeter initial mass function (IMF) and two top-heavy IMFs. In no case are we able to produce the postulated large Y ~ 0.35 without violating the observational constraint that the CNO content is nearly constant.Comment: accepted for publication in Ap

Partial mixing and the formation of 13C pockets in AGB stars: effects on the s-process elements

The production of the elements heavier than iron via slow neutron captures (the s process) is a main feature of the contribution of asymptotic giant branch (AGB) stars of low mass (< 5 Msun) to the chemistry of the cosmos. However, our understanding of the main neutron source, the 13C(alpha,n)16O reaction, is still incomplete. It is commonly assumed that in AGB stars mixing beyond convective borders drives the formation of 13C pockets. However, there is no agreement on the nature of such mixing and free parameters are present. By means of a parametric model we investigate the impact of different mixing functions on the final s-process abundances in low-mass AGB models. Typically, changing the shape of the mixing function or the mass extent of the region affected by the mixing produce the same results. Variations in the relative abundance distribution of the three s-process peaks (Sr, Ba, and Pb) are generally within +/-0.2 dex, similar to the observational error bars. We conclude that other stellar uncertainties - the effect of rotation and of overshoot into the C-O core - play a more important role than the details of the mixing function. The exception is at low metallicity, where the Pb abundance is significantly affected. In relation to the composition observed in stardust SiC grains from AGB stars, the models are relatively close to the data only when assuming the most extreme variation in the mixing profile.Comment: 17 pages, 8 figures, 6 tables, accepted for publications on Monthly Notices of the Royal Astronomical Societ