Luminosities of Carbon-rich Asymptotic Giant Branch stars in the Milky Way

Stars evolving along the Asymptotic Giant Branch can become Carbon-rich in the final part of their evolution. They replenish the inter-stellar medium with nuclear processed material via strong radiative stellar winds. The determination of the luminosity function of these stars, even if far from being conclusive, is extremely important to test the reliability of theoretical models. In particular, strong constraints on the mixing treatment and the mass-loss rate can be derived. We present an updated Luminosity Function of Galactic Carbon Stars obtained from a re-analysis of available data already published in previous papers. Starting from available near- and mid-infrared photometric data, we re-determine the selection criteria. Moreover, we take advantage from updated distance estimates and Period-Luminosity relations and we adopt a new formulation for the computation of Bolometric Corrections. This leads us to collect an improved sample of carbon-rich sources from which we construct an updated Luminosity Function. The Luminosity Function of Galactic Carbon Stars peaks at magnitudes around -4.9, confirming the results obtained in a previous work. Nevertheless, the Luminosity Function presents two symmetrical tails instead of the larger high luminosity tail characterizing the former Luminosity Function. The derived Luminosity Function of Galactic Carbon Stars matches the indications coming from recent theoretical evolutionary Asymptotic Giant Branch models, thus confirming the validity of the choices of mixing treatment and mass-loss history. Moreover, we compare our new Luminosity Function with its counterpart in the Large Magellanic Cloud finding that the two distributions are very similar for dust-enshrouded sources, as expected from stellar evolutionary models. Finally, we derive a new fitting formula aimed to better determine Bolometric Corrections for C-stars.Comment: 7 pages, 5 figures. To be published in Astronomy and Astrophysic

Heavy elements in Globular Clusters: the role of AGB stars

Recent observations of heavy elements in Globular Clusters reveal intriguing deviations from the standard paradigm of the early galactic nucleosynthesis. If the r-process contamination is a common feature of halo stars, s-process enhancements are found in a few Globular Clusters only. We show that the combined pollution of AGB stars with mass ranging between 3 to 6 M$_\odot$ may account for most of the features of the s-process overabundance in M4 and M22. In these stars, the s process is a mixture of two different neutron-capture nucleosynthesis episodes. The first is due to the 13C(a,n)16O reaction and takes place during the interpulse periods. The second is due to the 22Ne(a,n)25Mg reaction and takes place in the convective zones generated by thermal pulses. The production of the heaviest s elements (from Ba to Pb) requires the first neutron burst, while the second produces large overabundances of light s (Sr, Y, Zr). The first mainly operates in the less-massive AGB stars, while the second dominates in the more-massive. From the heavy-s/light-s ratio, we derive that the pollution phase should last for $150\pm 50$ Myr, a period short enough compared to the formation timescale of the Globular Cluster system, but long enough to explain why the s-process pollution is observed in a few cases only. With few exceptions, our theoretical prediction provides a reasonable reproduction of the observed s-process abundances, from Sr to Hf. However, Ce is probably underproduced by our models, while Rb and Pb are overproduced. Possible solutions are discussed.Comment: Accepted by the Ap

On the need of the Light Elements Primary Process (LEPP)

Extant chemical evolution models underestimate the Galactic production of Sr, Y and Zr as well as the Solar System abundances of s-only isotopes with 90<A<130. To solve this problem, an additional (unknown) process has been invoked, the so-called LEPP (Light Element Primary Process). In this paper we investigate possible alternative solutions. Basing on Full Network Stellar evolutionary calculations, we investigate the effects on the Solar System s-only distribution induced by the inclusion of some commonly ignored physical processes (e.g. rotation) or by the variation of the treatment of convective overshoot, mass-loss and the efficiency of nuclear processes. Our main findings are: 1) at the epoch of the formation of the Solar System, our reference model produces super-solar abundances for the whole s-only distribution, even in the range 90<A<130; 2) within errors, the s-only distribution relative to 150Sm is flat; 3) the s-process contribution of the less massive AGB stars (M<1.5 M_SUN) as well as of the more massive ones (M>4.0 M_SUN) are negligible; 4) the inclusion of rotation implies a downward shift of the whole distribution with an higher efficiency for the heavy s-only isotopes, leading to a flatter s-only distribution; 5) different prescriptions on convection or mass-loss produce nearly rigid shifts of the whole distribution. In summary, a variation of the standard paradigm of AGB nucleosynthesis would allow to reconcile models predictions with Solar System s-only abundances. Nonetheless, the LEPP cannot be definitely ruled out, because of the uncertainties still affecting stellar and Galactic chemical evolution models.Comment: Accepted for publication on Ap

Nucleation of small silicon carbide dust clusters in AGB stars

Silicon carbide (SiC) grains are a major dust component in carbon-rich AGB stars. The formation pathways of these grains are, however, not fully understood.\ We calculate ground states and energetically low-lying structures of (SiC)$_n$, $n=1,16$ clusters by means of simulated annealing (SA) and Monte Carlo simulations of seed structures and subsequent quantum-mechanical calculations on the density functional level of theory. We derive the infrared (IR) spectra of these clusters and compare the IR signatures to observational and laboratory data.\ According to energetic considerations, we evaluate the viability of SiC cluster growth at several densities and temperatures, characterising various locations and evolutionary states in circumstellar envelopes.\ We discover new, energetically low-lying structures for Si$_{4}$C$_{4}$, Si$_{5}$C$_{5}$, Si$_{15}$C$_{15}$ and Si$_{16}$C$_{16}$, and new ground states for Si$_{10}$C$_{10}$ and Si$_{15}$C$_{15}$. The clusters with carbon-segregated substructures tend to be more stable by 4-9 eV than their bulk-like isomers with alternating Si-C bonds. However, we find ground states with cage ("bucky"-like) geometries for Si$_{12}$C$_{12}$ and Si$_{16}$C$_{16}$ and low-lying, stable cage structures for n $\ge$ 12. The latter findings indicate thus a regime of clusters sizes that differs from small clusters as well as from large-scale crystals. Thus, and owing to their stability and geometry, the latter clusters may mark a transition from a quantum-confined cluster regime to crystalline, solid bulk-material. The calculated vibrational IR spectra of the ground-state SiC clusters shows significant emission. They include the 10-13 $\mu$m wavelength range and the 11.3 $\mu$m feature inferred from laboratory measurements and observations, respectively, though the overall intensities are rather low.Comment: 16 pages, 25 figures, 3 tables, accepted for publication in Ap

Interplay between pulsation, mass loss, and third dredge-up: More about Miras with and without technetium

We follow-up on a previous finding that AGB Mira variables containing the 3DUP indicator technetium (Tc) in their atmosphere form a different sequence of K-[22] colour as a function of pulsation period than Miras without Tc. A near- to mid-infrared colour such as K-[22] is a good probe for the dust mass-loss rate of the stars. Contrary to what might be expected, Tc-poor Miras show redder K-[22] colours (i.e. higher dust mass-loss rates) than Tc-rich Miras at a given period. Here, the previous sample is extended and the analysis is expanded towards other colours and dust spectra. The most important aim is to investigate if the same two sequences can be revealed in the gas mass-loss rate. We analysed new optical spectra and expanded the sample by including more stars from the literature. Near- and mid-IR photometry and ISO dust spectra of our stars were investigated. Literature data of gas mass-loss rates of Miras and semi-regular variables were collected and analysed. Our results show that Tc-poor Miras are redder than Tc-rich Miras in a broad range of the mid-IR, suggesting that the previous finding based on the K-[22] colour is not due to a specific dust feature in the 22 micron band. We establish a linear relation between K-[22] and the gas mass-loss rate. We also find that the 13 micron feature disappears above K-[22]~2.17 mag, corresponding to \dot{M}_{\rm g}\sim2.6\times10^{-7}M_{\sun}yr^{-1}. No similar sequences of Tc-poor and Tc-rich Miras in the gas mass-loss rate vs. period diagram are found, most probably owing to limitations in the available data. Different hypotheses to explain the observation of two sequences in the P vs. K-[22] explain the observation of two sequences in the P vs. K-[22 diagram are discussed and tested, but so far none of them convincingly explains the observations. Nevertheless, we might have found an hitherto unknown but potentially important process influencing mass loss on the TP-AGB.Comment: 16 pages, 15 figures, 2 online tables, accepted for publication in A&

The chemical composition of carbon stars: The R-type stars

The aim of this work is to shed some light on the problem of the formation of carbon stars of R-type from a detailed study of their chemical composition. We use high-resolution and high signal-to-noise optical spectra of 23 R-type stars selected from the Hipparcos catalogue. The chemical analysis is made using spectral synthesis in LTE and state-of-the-art carbon-rich spherical model atmospheres. We derive their CNO content (including the carbon isotopic ratio), average metallicity, lithium, and light (Sr, Y, Zr) and heavy (Ba, La, Nd, Sm) s-element abundances. The observed properties of the stars (galactic distribution, kinematics, binarity, photometry and luminosity) are also discussed. Our analysis shows that late-R stars are carbon stars with identical chemical and observational characteristics than the normal (N-type) AGB carbon stars. We confirm the results of the sole previous abundance analysis of early-R stars by Dominy (1984, ApJS, 55, 27), namely: they are carbon stars with near solar metallicity showing enhanced nitrogen, low carbon isotopic ratios and no s-element enhancements. In addition, we have found that early-R stars have Li abundances larger than expected for post RGB tip giants. We also find that a significant number (aprox. 40 %) of the early-R stars in our sample are wrongly classified, being probably classical CH stars and normal K giants. In consequence, we suggest that the number of true R stars is considerably lower than previously believed. We briefly discuss the different scenarios proposed for the formation of early-R stars. The mixing of carbon during an anomalous He-flash is favoured, although no physical mechanism able to trigger that mixing has been found yet. The origin of these stars still remains a mystery.Comment: 15 pages, 8 figures, accepted for publication in Astronomy and Astrophysic

Magnetic-buoyancy Induced Mixing in AGB Stars: Presolar SiC Grains

Isotope ratios can be measured in presolar SiC grains from ancient Asymptotic Giant Branch (AGB) stars at permil-level (0.1\%) precision. Such precise grain data permit derivation of more stringent constraints and calibrations on mixing efficiency in AGB models than traditional spectroscopic observations. In this paper we compare SiC heavy-element isotope ratios to a new series of FRUITY models that include the effects of mixing triggered by magnetic fields. Based on 2D and 3D simulations available in the literature, we propose a new formulation, upon which the general features of mixing induced by magnetic fields can be derived. The efficiency of such a mixing, on the other hand, relies on physical quantities whose values are poorly constrained. We present here our calibration by comparing our model results with the heavy-element isotope data of presolar SiC grains from AGB stars. We demonstrate that the isotopic compositions of all measured elements (Ni, Sr, Zr, Mo, Ba) can be simultaneously fitted by adopting a single magnetic field configuration in our new FRUITY models.Comment: 9 pages, 4 figures, 1 table. Accepted for publication on ApJ

New constraints on the major neutron source in low-mass AGB stars

We compare updated Torino postprocessing asymptotic giant branch (AGB) nucleosynthesis model calculations with isotopic compositions of mainstream SiC dust grains from low-mass AGB stars. Based on the data-model comparison, we provide new constraints on the major neutron source, 13C({\alpha},n)16O in the He-intershell, for the s-process. We show that the literature Ni, Sr, and Ba grain data can only be consistently explained by the Torino model calculations that adopt the recently proposed magnetic-buoyancy-induced 13C-pocket. This observation provides strong support to the suggestion of deep mixing of H into the He-intershell at low 13C concentrations as a result of efficient transport of H through magnetic tubes.Comment: ApJ, accepte

The s-Process Nucleosynthesis in Low Mass Stars : Impact of the Uncertainties in the Nuclear Physics Determined by Monte Carlo Variations

© Springer Nature Switzerland AG 2019We investigated the impact of uncertainties in neutron-capture and weak reactions (on heavy elements) on the s-process nucleosynthesis in low-mass stars using a Monte-Carlo based approach. We performed extensive nuclear reaction network calculations that include newly evaluated temperature-dependent upper and lower limits for the individual reaction rates. Our sophisticated approach is able to evaluate the reactions that impact more significantly the final abundances. We found that -decay rate uncertainties affect typically nuclides near s-process branchings, whereas most of the uncertainty in the final abundances is caused by uncertainties in neutron capture rates, either directly producing or destroying the nuclide of interest. Combined total nuclear uncertainties due to reactions on heavy elements are approximately 50%