next previous
Up: Nucleosynthesis of light

8. Wind composition along TP-AGB phase

In this last section, we present our predictions concerning (i) the evolution of all the surface isotopic ratios along the TP-AGB phase (by including our evolutionary models and extrapolations up to the convective envelope removal) and (ii) the total amount of material ejected at different ages for all the species up to tex2html_wrap_inline13889. We refer to the previous sections for detailed explanations. Thus, in the following, we will just comment these predictions.

8.1. Surface isotopic ratios along TP-AGB

Figs. 24 (click here) to 37 (click here) display, for each computed AGB star, the evolution of the isotopic ratios as a function of the total remaining mass that decreases due to mass loss.

Let us stress that we stop our calculations at the last thermal pulse that occurs before the convective envelope mass is reduced below tex2html_wrap_inline13891. This stage indeed roughly corresponds to the PN ejection that we do not maintain to model here. This explains the quite large remaining total stellar masses in Figs. 24 (click here) to 37 (click here), especially for those computed with the highest mass loss rates.

The first comment concerns the formation of carbon stars. tex2html_wrap_inline13893 AGB stars with Z = 0.02 and 3 and tex2html_wrap_inline13897 AGB stars with Z = 0.005 become carbon stars during their TP-AGB phase. The surface tex2html_wrap_inline13901/tex2html_wrap_inline13903 becomes greater than unity more rapidly for (i) decreasing total mass and/or (ii) decreasing Z. This metallicity dependence is confirmed by the increase of C stars compared to M stars from the galactic bulge to LMC and SMC (Blanco et al. 1978). Our tex2html_wrap_inline13907 (tex2html_wrap_inline13909) AGB stars with Z = 0.02 (0.005) can also become C stars with our lowest mass loss rates. In such cases indeed, in spite of very efficient HBB transforming tex2html_wrap_inline13913 in tex2html_wrap_inline13915 and tex2html_wrap_inline13917, many 3DUP events can still occur with very reduced envelope masses, so that the relatively short inter-pulse phases are not long enough to substantially destroy tex2html_wrap_inline13919. We will see however, together with other isotopic ratios, that these very low mass loss rates seem actually extreme. Let us emphasize that the only super-lithium-rich star (see Sect. 7.2) we found to be also a carbon star is the tex2html_wrap_inline13921 with Z = 0.005 one (with our medium mass loss rate). However, as already mentioned in Sect. 6.2, observations indicate that Z = 0.02 carbon stars are formed easier (i.e. at lower luminosities) than predicted by all the models. Let us also stress the very high tex2html_wrap_inline13927/tex2html_wrap_inline13929 isotopic ratios that can be expected at the surface of tex2html_wrap_inline13931 tex2html_wrap_inline13933 AGB stars, whatever Z. On the other hand, more massive AGB stars experience HBB that progressively leads the tex2html_wrap_inline13937/tex2html_wrap_inline13939 isotopic ratio close to its CN cycle equilibrium value (formation of J stars). The case of tex2html_wrap_inline13941 stars with Z = 0.005 very clearly shows the competition between the 3DUP episodes (increasing tex2html_wrap_inline13945/tex2html_wrap_inline13947) and HBB (decreasing it). For lower mass loss rates, the 3DUP events are more numerous for the same total mass decrease, leading to a more rapid tex2html_wrap_inline13949/tex2html_wrap_inline13951 increase. However, the corresponding longer TP-AGB phase allows the base of the convective envelope to reach higher temperatures that finally destroy tex2html_wrap_inline13953 through the CN cycle. This efficient transformation of tex2html_wrap_inline13955 to tex2html_wrap_inline13957 also reduces the tex2html_wrap_inline13959/tex2html_wrap_inline13961 ratio, i.e. the possibility to form carbon stars. It happens if the temperatures at the base of the convective envelopes become higher than tex2html_wrap_inline13963 K. Such stars have corresponding surface luminosities > tex2html_wrap_inline13967 (tex2html_wrap_inline13969) for Z = 0.02 (0.005), i.e. tex2html_wrap_inline13973 (-6.45). We consequently predict that C stars cannot be formed at higher luminosities than this threshold, due to strong HBB. This is in agreement with Boothroyd et al. (1993) predictions. This also corresponds to the observed threshold for AGB stars of the LMC (Cohen et al. 1981).

Finally, note the surface pollution in tex2html_wrap_inline13975 that is more important, compared to that of tex2html_wrap_inline13977, for decreasing total masses. For tex2html_wrap_inline13979 AGB stars, the surface tex2html_wrap_inline13981/tex2html_wrap_inline13983 ratio is quite independent of the mass loss rate. This is due to the inter-pulse duration that is significantly greater than the tex2html_wrap_inline13985 lifetime, so that the tex2html_wrap_inline13987 decay is the dominant factor for the tex2html_wrap_inline13989/tex2html_wrap_inline13991 ratio evolution at the surface of such stars. The quite complicated behavior of the tex2html_wrap_inline13993/tex2html_wrap_inline13995 isotopic ratio at the surface of our most massive models is again due to the competition between the 3DUP that adds tex2html_wrap_inline13997 into the convective envelope and HBB that can partially destroy it. By extrapolation, it is reasonable to expect even lower tex2html_wrap_inline13999/tex2html_wrap_inline14001 ratios in low-mass AGB stars. As mentioned in Forestini et al. (1996), tex2html_wrap_inline14003 could become detectable for < tex2html_wrap_inline14007 evolved AGB stars.

Due to HBB, those stars showing very low surface tex2html_wrap_inline14009/tex2html_wrap_inline14011 ratios are also expected to have very high tex2html_wrap_inline14013/tex2html_wrap_inline14015 ones. It has to be stressed indeed that in our most massive AGB stars, the CN cycle operating inside the convective envelope produces large amounts of tex2html_wrap_inline14017.

The tex2html_wrap_inline14019/tex2html_wrap_inline14021 and tex2html_wrap_inline14023/tex2html_wrap_inline14025 surface isotopic ratios also evolve as expected from our discussions in Sect. 7. Note that quite low tex2html_wrap_inline14027/tex2html_wrap_inline14029 and very high tex2html_wrap_inline14031/tex2html_wrap_inline14033 ratios, which are clear signatures of the operation of the ON cycle at the base of the convective envelope (HBB), appear rather simultaneously with the CN cycle signatures in our most massive objects. Also remark the progressive tex2html_wrap_inline14035/tex2html_wrap_inline14037 decrease in our tex2html_wrap_inline14039 (tex2html_wrap_inline14041) AGB star with Z = 0.02 (0.005). For these objects, the temperature at the base of the convective envelope finally becomes high enough for tex2html_wrap_inline14045 to partially burn.

The tex2html_wrap_inline14047 surface enhancements we predict are displayed in a special form, corresponding to the observation presentations adopted by Jorissen et al. (1992). Compared to their Fig. 8, two conclusions can be stated.

As expected, the tex2html_wrap_inline14077/tex2html_wrap_inline14079 and tex2html_wrap_inline14081/tex2html_wrap_inline14083 surface isotopic ratios are much less sensitive to the 3DUP episodes and HBB. More precisely, only our most massive models with the lowest mass loss rates show large tex2html_wrap_inline14085 depletions due to a very strong HBB. Again, these are however quite extreme situations.

Last but not least, we predict substantial tex2html_wrap_inline14087 surface enhancements due to the repetitive dredges-up of the HBS and inter-shell regions. The tex2html_wrap_inline14089/tex2html_wrap_inline14091 isotopic ratio tends to decrease with increasing total mass due to (i) the greater dilution as mixing occurs in a more massive convective envelope and (ii) the thinner HBS and inter-shell regions (in mass). Nevertheless, tex2html_wrap_inline14093/tex2html_wrap_inline14095 could well be increasing again if tex2html_wrap_inline14097 was produced inside the convective envelope itself by strong HBB. Note however that such a very high tex2html_wrap_inline14099 production level is actually expected to be accompanied by very low Mg isotopic ratios, as tex2html_wrap_inline14101 is produced by the MgAl chain in the extreme situations mentioned above. Let us finally mention that Wasserburg et al. (1994) presented computations in which the evolution of the tex2html_wrap_inline14103 surface abundance has been followed self-consistently, by including the effect of neutron irradiation due to an enhanced amount of tex2html_wrap_inline14105 in the inter-shell region (to engender a s-process).

Our results compare very well with those obtained by Boothroyd et al. (1994, 1995) concerning the O isotopic ratios in intermediate-mass stars. They also mention the problem to reproduce the tex2html_wrap_inline14107/tex2html_wrap_inline14109 ratio, especially in low-mass carbon AGB stars. We cannot go further in this comparison without having yet computed the TP-AGB phase of low-mass stars. Boothroyd et al. (1995) already suggested that to conciliate all the data with observations of low-mass stars, one probably has to take into account of slow-particle transports inside radiative zones, what they called ``cool-bottom precessing''. That is what we plan to do (see Sect. 9).

 figure3334
Figure 24:   Surface isotopic ratios as a function of the remaining total stellar mass, that decreases due to the mass loss along the TP-AGB. The solid line refers to the standard mass loss rate (see Sect. 6.3.2 for the corresponding tex2html_wrap_inline14111 parameter values in the Reimers relation). Dotted lines correspond to mass loss rates increased or decreased by a factor of two. The extrapolated thermal pulses (and 3DUP episodes) are indicated by the filled circles, and are of course more numerous for lower mass loss rates. Here are indicated the tex2html_wrap_inline14113/tex2html_wrap_inline14115, tex2html_wrap_inline14117/tex2html_wrap_inline14119, tex2html_wrap_inline14121/tex2html_wrap_inline14123, tex2html_wrap_inline14125/tex2html_wrap_inline14127, tex2html_wrap_inline14129/tex2html_wrap_inline14131 and tex2html_wrap_inline14133/tex2html_wrap_inline14135 ratios, in the case of the tex2html_wrap_inline14137 (Z = 0.02) star. All ratios are in mass fractions, except for the tex2html_wrap_inline14141/tex2html_wrap_inline14143 ratio (given in number abundances)

 figure3367
Figure 25:   tex2html_wrap_inline14145/tex2html_wrap_inline14147, tex2html_wrap_inline14149/tex2html_wrap_inline14151, tex2html_wrap_inline14153/tex2html_wrap_inline14155 ratios as well as [tex2html_wrap_inline14157/tex2html_wrap_inline14159] as a function of tex2html_wrap_inline14161/tex2html_wrap_inline14163, in the case of the tex2html_wrap_inline14165 (Z = 0.02) star. All ratios are in mass fractions and [A/B] means tex2html_wrap_inline14169

 figure3396
Figure 26:   Same as Fig. 24 (click here), in the case of the tex2html_wrap_inline14171 (Z = 0.02) star

 figure3402
Figure 27:   Same as Fig. 25 (click here), in the case of the tex2html_wrap_inline14175 (Z = 0.02) star

 figure3408
Figure 28:   Same as Fig. 24 (click here), in the case of the tex2html_wrap_inline14179 (Z = 0.02) star

 figure3414
Figure 29:   Same as Fig. 25 (click here), in the case of the tex2html_wrap_inline14183 (Z = 0.02) star

 figure3420
Figure 30:   Same as Fig. 24 (click here), in the case of the tex2html_wrap_inline14187 (Z = 0.02) star

 figure3426
Figure 31:   Same as Fig. 25 (click here), in the case of the tex2html_wrap_inline14191 (Z = 0.02) star

 figure3432
Figure 32:   Same as Fig. 24 (click here), in the case of the tex2html_wrap_inline14195 (Z = 0.005) star

 figure3438
Figure 33:   Same as Fig. 25 (click here), in the case of the tex2html_wrap_inline14199 (Z = 0.005) star

 figure3444
Figure 34:   Same as Fig. 24 (click here), in the case of the tex2html_wrap_inline14203 (Z = 0.005) star

 figure3450
Figure 35:   Same as Fig. 25 (click here), in the case of the tex2html_wrap_inline14207 (Z = 0.005) star

 figure3456
Figure 36:   Same as Fig. 24 (click here), in the case of the tex2html_wrap_inline14211 (Z = 0.005) star

 figure3462
Figure 37:   Same as Fig. 25 (click here), in the case of the tex2html_wrap_inline14215 (Z = 0.005) star

8.2. Other comparisons with observations

Evolutionary models of TP-AGB stars can be strongly constrained by the confrontation between predicted and observed surface isotopic ratios, especially if different ratios are observed for the same stars. For most of the observed AGB stars, only the C and O isotopic ratios have been determined (see e.g. Harris & Lambert 1984; Harris et al. 1985; Lambert et al. 1986; Harris et al. 1987; Harris & Lambert 1987; Kahane et al. 1992). Globally, we found agreements. General features come out.

Presently, there is only one star for which all the isotopic ratios have been determined. It is Cw Leo. This C star is one of the nearest evolved AGB star ever observed. Its distance is estimated to tex2html_wrap_inline14363 pc (see e.g. Claussen et al. 1987; Keady et al. 1988; Griffin 1990). This corresponds to a surface luminosity tex2html_wrap_inline14365 tex2html_wrap_inline14367 and a present mass loss rate of roughly tex2html_wrap_inline14369 tex2html_wrap_inline14371tex2html_wrap_inline14373. Such a high wind makes it quite optically thick. A very extended circumstellar envelope has been detected to surround Cw Leo, IRC+10216. IRC+10216 shows evidences for a bi-polar structure (see e.g. Guélin et al. 1993), attesting that Cw Leo is very evolved and close to the end of its TP-AGB phase. Consequently, its present convective envelope mass is probably very reduced (tex2html_wrap_inline14375 tex2html_wrap_inline14377 typically).

The circumstellar envelope of Cw Leo allowed very accurate determinations of the isotopic ratios (see Kahane et al. 1988 and 1992; Guélin et al. 1995; Forestini et al. 1996): tex2html_wrap_inline14379/tex2html_wrap_inline14381 tex2html_wrap_inline14383, tex2html_wrap_inline14385/tex2html_wrap_inline14387 > 1400, tex2html_wrap_inline14391/tex2html_wrap_inline14393 > 5300, tex2html_wrap_inline14397/tex2html_wrap_inline14399 tex2html_wrap_inline14401, tex2html_wrap_inline14403/tex2html_wrap_inline14405 tex2html_wrap_inline14407, tex2html_wrap_inline14409/tex2html_wrap_inline14411 tex2html_wrap_inline14413, tex2html_wrap_inline14415/tex2html_wrap_inline14417 tex2html_wrap_inline14419 and tex2html_wrap_inline14421/tex2html_wrap_inline14423 tex2html_wrap_inline14425 (to be confirmed). Furthermore, the Si isotopic ratios appear to have interstellar medium values, i.e. unaltered by the star. As derived in Guélin et al. (1995) and Forestini et al. (1996), almost all these isotopic ratios seem to be consistent (within the error bars), following our predictions, with an AGB star that had a main sequence total mass between 4 and tex2html_wrap_inline14427 (assuming Z = 0.02; see Figs. 26 (click here) and 29 (click here)). However, tex2html_wrap_inline14431/tex2html_wrap_inline14433 is again not well explained. As various surface isotopic ratios (especially tex2html_wrap_inline14435/tex2html_wrap_inline14437 and tex2html_wrap_inline14439/tex2html_wrap_inline14441) are significantly changing between 4 and tex2html_wrap_inline14443 AGB models, we are presently computing the TP-AGB phase of a tex2html_wrap_inline14445 star with Z = 0.02 in order to improve this last prediction.

Finally, another very interesting way to constraint stellar evolution models of AGB stars is to compare the predicted surface isotopic ratios with those determined in primitive meteorites (mainly through grains included inside carbon rich chondrites). Most of the tex2html_wrap_inline14449/tex2html_wrap_inline14451, tex2html_wrap_inline14453/tex2html_wrap_inline14455 and tex2html_wrap_inline14457/tex2html_wrap_inline14459 isotopic anomalies measured in grains (see e.g. Ott 1991; Zinner et al. 1991a) are consistent with those coming from the surface of evolved TP-AGB stars of intermediate-mass. This is of course in favor of the presence, when the solar system has been formed, of one or a few AGB stars in its surroundings. Such grains are indeed well known to be formed in the cool atmosphere of evolved AGB stars (i.e. when mass loss rates become rather high). Concerning the measured O isotopic ratios, the situation is not so clear as many grains have ratios significantly out of the corresponding predicted ranges. Some of these grains at least can have been formed in low-mass AGB stars. Again, as mentioned in Sect. 8.1, Boothroyd et al. (1995) recently suggested that in such stars, cool-bottom processing can help to reduce these discrepancies. On the contrary, some other grains perfectly agree with our predictions. This is for example the case of a tex2html_wrap_inline14461 tex2html_wrap_inline14463 grain of the Bishunpur LL3.1 chondrite by Huss et al. (1994). The derived ratios, tex2html_wrap_inline14465/tex2html_wrap_inline14467 tex2html_wrap_inline14469, tex2html_wrap_inline14471/tex2html_wrap_inline14473 tex2html_wrap_inline14475 and tex2html_wrap_inline14477/tex2html_wrap_inline14479 tex2html_wrap_inline14481, are in perfect agreement with ratios corresponding to a 3 - tex2html_wrap_inline14485 AGB star of nearly solar metallicity. Finally, other SiC tex2html_wrap_inline14487-rich grains raise problems, like that discovered by Zinner et al. (1991b), exhibiting a very high tex2html_wrap_inline14489/tex2html_wrap_inline14491 ratio of tex2html_wrap_inline14493. It seems hard to explain such a value from our models. However, this grain also exhibits an anomalously large tex2html_wrap_inline14495/tex2html_wrap_inline14497 isotopic ratio that seems to exclude its AGB star origin.

8.3. Yields

Tables 8 (click here) to 10 (click here) give the total mass ejected by the wind (in tex2html_wrap_inline14513) for various elements (i) at the top of the first RGB, (ii) at the beginning of the TP-AGB phase and (iii) at the end of the TP-AGB phase (thus including our extrapolated nucleosynthesis calculations). As such data mainly concerns the chemical evolution of galaxies, we also present the net yield (positive or negative) for each element. The data of Tables 8 (click here) to 10 (click here) correspond to the medium mass loss rates (see Sect. 8.1), i.e. the solid lines of Figs. 24 (click here) to 37 (click here). Let us briefly comment.

  Table 8:   Total mass ejected (in tex2html_wrap_inline14515) in the wind for some nuclides at the RGB top, at the end of the E-AGB phase and at the AGB tip, just before the ejection of the planetary nebula. Also given is the net contribution for these elements, i.e. the total final mass ejected reduced by the amount of matter that would have been ejected with the initial composition. Here are presented elements from He to N

  Table 9:   Same as Table 8 for elements from O to Na

  Table 10:   Same as Table 8 for elements from Mg to Si

8.3.1. He

tex2html_wrap_inline15929 is produced during the pre-main sequence and the main sequence phases. For the stars with masses lower than tex2html_wrap_inline15931 (tex2html_wrap_inline15933) with Z = 0.02 (0.005), tex2html_wrap_inline15937 survives during the following phases of evolution and is injected in the ISM by stellar wind, mainly during the TP-AGB ultimate phase. These objects thus contribute, while more modestly than low-mass (< tex2html_wrap_inline15941) stars to the galactic enrichment of tex2html_wrap_inline15943. In more massive stars however, tex2html_wrap_inline15945 is progressively burnt at the benefits of tex2html_wrap_inline15947 and tex2html_wrap_inline15949 by HBB. At the end of their TP-AGB phase, tex2html_wrap_inline15951 can even be completely destroyed, so that the corresponding net yields are negative.

Intermediate-mass stars contribute to the galactic enrichment of tex2html_wrap_inline15953. This pollution mainly occurs at the end of the TP-AGB phase and is higher for the highest stellar masses and lower metallicity models. Indeed, tex2html_wrap_inline15955 is brought to the surface during the successive 3DUP episodes. Moreover, in our most massive stars, HBB further contributes to a non-negligible production of this element inside the convective envelope.

8.3.2. tex2html_wrap_inline15963

TP-AGB stars appear to be major contributors to the tex2html_wrap_inline15965 galactic enrichment. The very nuclearly fragile tex2html_wrap_inline15967 is destroyed at very low temperatures in stellar interiors. As was discussed in Sect. 7.2, it can however be produced at the expense of tex2html_wrap_inline15969 when HBB is effective. As a result, the global yield of tex2html_wrap_inline15971 is positive in our 5 and tex2html_wrap_inline15973 stars with Z = 0.02, whereas it is negative for the other stars. Note that the tex2html_wrap_inline15977 production in our tex2html_wrap_inline15979 AGB star with Z = 0.005 is not large enough to compensate its depletion during the previous stellar evolution stages.

8.3.3. CNO elements and tex2html_wrap_inline15987

As already mentioned (see Sects. 6.2 and 8.1) 3DUP events significantly pollute the convective envelope in tex2html_wrap_inline15989. However, if HBB occurs, CNO-burning leads to primary production of tex2html_wrap_inline15991 and tex2html_wrap_inline15993 at the expense of tex2html_wrap_inline15995 (and possibly tex2html_wrap_inline15997). The net yields of tex2html_wrap_inline15999 thus decrease for increasing stellar masses and even become negative for our 5 and tex2html_wrap_inline16001 (tex2html_wrap_inline16003) stars with Z = 0.02 (0.005). Owing to that, the net yields of tex2html_wrap_inline16007 and tex2html_wrap_inline16009 rapidly increase with initial stellar mass and metallicity.

The more massive an AGB star, the less efficient the tex2html_wrap_inline16011 production (see Sects. 7.3 and 8.1). Due to HBB, our most massive models even present negative net yields for that element. As explained in Sect. 8.1, low-mass AGB stars are probably the most efficient tex2html_wrap_inline16013 producers in galaxies.

8.3.4. Ne and Na

As a net result of the TP-AGB phase, intermediate-mass stars somewhat contribute to the ISM enrichment in tex2html_wrap_inline16019 and, to a still lower extent, tex2html_wrap_inline16021, to the expense of tex2html_wrap_inline16023.

8.3.5. Mg and Al

Massive AGB stars (higher than tex2html_wrap_inline16025) which experience strong HBB are the most efficient sources of tex2html_wrap_inline16027 among the AGB stars. They even significantly destroy tex2html_wrap_inline16029. In addition, tex2html_wrap_inline16031 is also substantially produced by all our models, except the tex2html_wrap_inline16033 star with Z = 0.005. We also note a non-negligible tex2html_wrap_inline16037 production inside the same massive AGB stars.


next previous
Up: Nucleosynthesis of light

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr