We compare in Fig. 3 (click here) the yields predicted by the standard model with those given by RV for stars formed with initial metallicities Z = 0.02 (Y = 0.28) and Z = 0.004 (Y = 0.232), respectively. We take into account the observational fact that HBB takes place in stars more massive than . Accordingly, we use the yields given by RV in case Z = 0.02 with for m > 3.3 and for (i.e. their Tables 3 (click here)e and 3 (click here)a, respectively). Similarly, we use their yields in case Z = 0.004 with for m > 3.25 and for m < 3 (i.e. their tables 3i and 3h). For Z=0.004 RV did not tabulate results for .
Figure 3 (click here) shows that the standard model (pre-AGB evolution as described in Sect. 3) results in yields which differ from the selected yields of RV within a factor 2-3. The standard model predicts larger yields predominantly for AGB stars with in case of carbon and oxygen, and for AGB stars more massive than in case of nitrogen. In this comparison, we have neglected the effect on the resulting yields of differences (up to ) in the initial C, N, and O abundances between the standard model and that of RV.
It was derived in GJ that is needed to fit the observed initial-final mass relation for stars in the Galactic disk (see Weidemann & Koester 1983). However, the selected yields of RV were computed for while for the standard model . Therefore, due to the high mass loss rates of AGB stars much fewer thermal pulses on the AGB occur in the standard model compared to the RV models. The main part of the observed differences between the yields predicted by the standard model and those of RV are probably due to this effect (apart from differences in the detailed description of evolution along the AGB, in particular the efficiency of third dredge-up).
In conclusion, we find that the selected yields given by RV differ from those predicted by the standard model by a factor 2-3. In particular, for high mass AGB stars ( ), the effect of HBB on the nitrogen yields for the selected RV models is much larger than that for the standard model. This suggests that values of the mixing length parameter may be more appropriate for massive AGB stars as we will argue below.
We have presented the yields of intermediate mass AGB stars for appropriate ranges in mass, initial composition, mass loss parameter , and effects of second dredge-up and HBB. We have shown that the yields of such stars are determined by their final stages and are important for the carbon and nitrogen enrichment of the Galactic disk ISM. In particular, AGB stars account probably for more than of the interstellar nitrogen in the disk (depending on the shape of the IMF at low and intermediate mass stars).
From the results presented in GJ and GHJ, we argued that the standard model with provides a reasonable approximation of the yields of intermediate mass AGB stars in the Galactic disk and the Large Magellanic Cloud. These systems have a metallicity that differ by only a factor of 2 (e.g. Russell & Dopita 1992). In galaxies with a substantial lower metallicity, one may expect a lower value of to be more appropriate. We like to emphasize that using a fixed value of does not necessarily mean identical mass loss rates as two stars of the same initial mass evolve differently in the synthetic model due to the explicit metallicity dependence of the recipes used.
Direct observational information on the metallicity dependence of mass loss and element yields in AGB stars is rare. In Groenewegen et al. (1995), the spectral energy distributions and spectra of three long-period variables (one each in the SMC, LMC and Galaxy) with roughly the same period were fitted. From the derived ratios of the dust optical depths in these stars, it was argued that the mass loss rates of AGB stars in the Galaxy, LMC, and SMC are roughly in the ratio of 4:3:1. This corroborates that could be similar for AGB stars in the Galaxy and LMC. Furthermore, this suggests that for AGB stars in low metallicity systems like the SMC (Z 0.004), values of 1-2 may be more appropriate.
In GHJ, we compared the mean abundances in the envelopes of AGB stars predicted by the standard model with the abundances observed in PNe in the Galactic disk. Here we repeat part of this analysis with improved model input and put emphasis on the differences in the description of pre-AGB evolution between the Geneva models and that outlined in GJ. In particular, we consider in more detail the effects of second dredge-up and HBB on the predicted abundances and address the uncertainties involved.
In the model, the abundances within PNe are estimated by averaging the abundances in the ejecta of AGB stars over the final 25000 yr (e.g. Pottasch 1995). We neglect any changes in the ejected shell abundances during the post-AGB phase, e.g. due to a late thermal pulse (Schönberner 1983), which is expected to be a rare event, or due to selective element depletion by dust formation. The latter process may affect the composition both in the wind of an AGB star and during the post-AGB phase (e.g. Bond 1992; van Winckel et al. 1992) but is neglected here for simplicity.
We assume an upper mass limit of 8 for stars that ultimately can become a PN (with final core mass less than 1.2 ) and ignore the possibility that not all our model AGB stars will become PNe. In fact, some of the low-mass AGB stars may evolve so slowly during the post-AGB phase that the material previously collected in the wind is dispersed before the central star has become hot enough to ionize this material. Also, the upper mass limit for AGB stars is matter of debate and may range between 6 and , depending on the critical mass for carbon ignition in an electron degenerate core and on details of the stellar mass-loss scenario (cf. GJ; Vassiliadis & Wood 1993; Hashimoto et al. 1993). Furthermore, we assume a constant value of yr. In reality, the time during which the mass accumulated in a PN has been swept up on the AGB may depend on the mass and initial composition of the progenitor. Nevertheless, we do not expect that these simplifications will alter our qualitative conclusions given below.
The observed PNe abundances are taken from various sources, i.e. mainly from Aller & Cryzack (1983), Zuckerman & Aller (1986), Aller & Keyes (1987), and Kaler et al. (1990). The few halo PNe are excluded as the present comparison concentrates on AGB stars in the Galactic disk. Errors in the observed abundances are typically 0.015 in He/H and about 0.2-0.25 dex in all other ratios considered in Fig. 4 (click here).
Figure 4: Planetary nebulae abundances (by number) predicted by the standard
model with pre-AGB evolution according to the Geneva tracks (solid
curves) and according to the recipes outlined in Sect. 3 (dotted
curves). The latter model without HBB is shown for comparison (dashed
curves). Abundances observed in PNe in the Galactic disk are shown by
open circles (data mainly from Aller & Cryzack (1983),
Zuckerman &
Aller (1986), Aller & Keyes (1987), and
Kaler et al. (1990)). Typical errors
in the observations are indicated at the bottom right corner of each panel
The PNe nowadays observed in the Galactic disk probably originate from AGB stars covering a wide range in initial mass, i.e. . This means that the progenitors of these PNe were formed at galactic ages ranging from about 10-15 Gyr to 50 Myr ago (see e.g.\ Schaller et al. 1992). Therefore, the initial element abundances of these PN progenitors are expected to differ considerably since the enrichment of the Galactic disk ISM over this time interval has been substantial (e.g. Twarog 1980; Edvardsson et al. 1993). When comparing the abundances predicted in the envelopes of final stage AGB stars with those observed in PNe, we take this important effect into account by incorporating a self-consistent model for the chemical evolution of the Galactic disk (van den Hoek et al. 1997; GHJ). Since the metallicity dependent AGB yields and the resulting chemical evolution of the Galactic disk are mutually dependent, an iterative solution method was applied. The adopted star formation history (SFR) and initial mass function (IMF) in this model were derived using observational constraints to the abundance-abundance variations with age of stars in the solar neighbourhood, the metallicity and age distributions of long-living stars as well as constraints to the current space density and formation rate of several post-main-sequence star populations.
Resulting abundance-ratios (by number) in PNe are shown in Fig. 4 (click here) in case of the standard model assuming pre-AGB evolution according to the Geneva tracks (Tables 2 (click here)-20 (click here)). We verified that the resulting abundances are insensitive to the adopted PN lifetime up to 50000 yr. In general, good agreement is found between the observed and predicted PN abundances despite the uncertainties involved. We find that the overall trend of the observations is reproduced well by the standard model independent of the adopted chemical enrichment history of the Galactic disk. However, some discrepancies are present between the standard model (with pre-AGB evolution according to the Geneva tracks) and the observations, in particular at large values of He/H 0.15.
For comparison, we show in Fig. 4 (click here) the PN abundances predicted by the standard model with pre-AGB evolution according to the recipes outlined in Sect. 3. In this case, the enhanced effect of second dredge-up can account for massive AGB stars with He/H up to 0.18 in their envelopes. This suggests that second dredge-up has been relatively important at least for some of the PNe in our sample with He/H 0.15. Alternatively, a substantial fraction of the hydrogen contained in the outer envelope may have turned into helium. Since PNe may evolve from a H and/or He-shell burning AGB star, this will determine the distribution of He/H abundance ratios observed for a given progenitor mass. We have included in Table 31 (click here) the yields of H and He for the standard model with second dredge-up as described by RV (cf. Sect. 3.2) which provides reasonable agreement with the observed PN abundances of He/H up to 0.2, in particular for the more massive PNe.
The effect of HBB on the predicted abundances can be seen in Fig. 4 (click here) by comparison of the standard model with and 1.3 (i.e. no HBB), respectively. Our results indicate that the standard model overestimates the effect of HBB on the resulting N/O abundance ratios in PNe with progenitors mass . We note that the standard model takes into account the maximum effect of HBB as described by RV so that values of the mixing length parameter < 2 in case of RV are probably more appropriate for massive AGB stars. On the other hand, models without HBB are inconsistent with the observed N/O abundances as well as with independent observations discussed in Sect. 4.2. Therefore, the range of N/O abundances observed in the envelopes of post-AGB stars allows for variations in the importance of HBB roughly covering the range from to 0.9 . In case of reduced HBB (i.e. ), the CNO yields of massive stars in Tables 21 (click here)-30 (click here) are more suitable than those given for the standard model.
The procedure to approximate the effect of HBB in a semi-analytical way has been described in the Appendix of GJ. Here the basic parameters were determined by fitting the RV (, ) case for which the effect of HBB is largest. Thus, as the standard model has , possible effects of HBB varying with in particular mass loss were neglected. In fact, the temperature structure of the envelope is expected to change when the number of thermal pulses decreases with increasing values of . This may reduce the amount of HBB occuring in the convective envelope and affect the resulting abundances as observed for PNe with log (N/O) and He/H (cf. Fig. 4 (click here)).
We emphasize that the resulting abundances of PNe do depend strongly on the initial element abundances of their progenitors, i.e. are very sensitive to the detailed chemical enrichment of the Galactic disk. A considerable part of the scatter observed in Fig. 4 (click here) is expected to be caused (in addition to experimental errors) by substantial variations in the initial abundances of PN progenitors due to the inhomogeneous chemical evolution of the Galactic disk ISM (e.g. van den Hoek & de Jong 1997). Furthermore, the progenitors of the PNe nowadays observed in the solar neighbourhood probably have been formed over a large range in galactocentric distance (e.g. Wielen et al. 1996) and thus with a large range in initial metallicity according to the radial abundance gradients in the disk ISM (e.g. Shaver et al. 1983). Therefore, we expect the agreement between the predicted and observed PN abundance-ratios to improve further when averaging over a range in initial composition for a given progenitor mass.
We conclude that the abundance-ratios predicted by the standard model are consistent with the observed abundances in virtually all the PNe in our sample when we allow for plausible variations in strength of second dredge-up and HBB.
The primary application of the stellar yields presented in this paper is probably in galactic chemical evolution studies.
As models with the default parameters for the mass loss on the AGB,
third dredge-up efficiency, and HBB fit many constrains in our
galaxy and the LMC (GJ/GHJ), the corresponding yields (Tables 2 (click here)-20 (click here)) are
probably the most appropriate ones to use. Possible alternatives are
models with less HBB (Tables 21 (click here)-30 (click here)), or using the pre-AGB evolution from the
synthetic model (Table 31 (click here)).
We argued in Sect. 5.2 that the scaling factor of the Reimers
law may be different for low metallicities. To simulate this effect one may
want to use the yields for models with for Z =
0.004 and for Z = 0.001 (Tables 32 (click here)-37 (click here)).
Acknowledgements
We like to thank the referee Georges Meynet for careful reading of the paper and encouraging remarks. It is a pleasure to thank Achim Weiss for his comments on an earlier version of this paper. The research of LBH and MG has been supported under grants 782-372-028 and 782-373-030 by the Netherlands Foundation for Research in Astronomy (ASTRON) which is financially supported by the Netherlands Organisation for Scientific Research (NWO).
Table 2:
Pre-AGB yields for ">, ,
and ">
Table 3:
AGB yields for ">, ,
and ">
Table 4:
Final AGB yields for ">, ,
and ">
Table 5:
Total yields for ">, ,
and ">
Table 6:
Pre-AGB yields for ">,
">, and
Table 7:
AGB yields for ">,
">, and
Table 8:
Final AGB yields for ">,
">, and
Table 9:
Total yields for ">,
">, and
Table 10:
Pre-AGB yields for ">,
">, and
Table 11:
AGB yields for ">, , and
">
Table 12:
Final AGB yields for ">,
">, and
Table 13:
Total yields for ">,
">, and
Table 14:
Pre-AGB yields for ">,
">, and
Table 15: AGB yields for ">,
">, and
Table 16:
Final AGB yields for ">,
">, and
Table 17: Total yields for ">,
">, and
Table 18: Pre-AGB yields for
">,
">, and
Table 19:
AGB yields for ">,
">, and
Table 20:
Final AGB yields for ">, ,
and ">
Table 21:
Total yields for ">,
">, and
Table 22:
AGB yields for ">,
">, and
Table 23:
AGB yields for ">,
">, and
Table 24:
AGB yields for ">,
">, and
Table 25:
AGB yields for ">,
">, and
Table 26:
AGB yields for ">,
">, and
Table 27:
Final AGB yields for ">,
">, and
Table 28:
AGB yields for ">,
">, and
Table 29:
Final AGB yields for ">,
">, and
Table 30: AGB yields for ">,
">, and
Table 31: Final AGB yields for ">,
">, and
">
Table 32: Total AGB yields for H, He for synthetic evolution model
Table 33:
AGB yields for ">,
">, and
Table 34:
Final AGB yields for ">,
">, and
Table 35:
Total yields for ">,
">, and
Table 36:
AGB yields for ">,
">, and
Table 37: Final AGB yields for ">,
">, and
Table 38:
Total yields for ">,
">, and