next previous
Up: New theoretical yields

5. Discussion

5.1. Comparison with RV

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 tex2html_wrap_inline3841 tex2html_wrap_inline3843. Accordingly, we use the yields given by RV in case Z = 0.02 with tex2html_wrap_inline3847 for m > 3.3 tex2html_wrap_inline3851 and tex2html_wrap_inline3853 for tex2html_wrap_inline3855 tex2html_wrap_inline3857 (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 tex2html_wrap_inline3861 for m > 3.25 tex2html_wrap_inline3865 and tex2html_wrap_inline3867 for m < 3 tex2html_wrap_inline3871 (i.e. their tables 3i and 3h). For Z=0.004 RV did not tabulate results for tex2html_wrap_inline3875.

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 tex2html_wrap_inline3879 tex2html_wrap_inline3881 in case of carbon and oxygen, and for AGB stars more massive than tex2html_wrap_inline3883 tex2html_wrap_inline3885 in case of nitrogen. In this comparison, we have neglected the effect on the resulting yields of differences (up to tex2html_wrap_inline3887) in the initial C, N, and O abundances between the standard model and that of RV.

It was derived in GJ that tex2html_wrap_inline3889 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 tex2html_wrap_inline3891 while for the standard model tex2html_wrap_inline3893. 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 (tex2html_wrap_inline3897 tex2html_wrap_inline3899), 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 tex2html_wrap_inline3901 may be more appropriate for massive AGB stars as we will argue below.

5.2. AGB stars and the enrichment of the galactic ISM

We have presented the yields of intermediate mass AGB stars for appropriate ranges in mass, initial composition, mass loss parameter tex2html_wrap_inline3903, 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 tex2html_wrap_inline3905 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 tex2html_wrap_inline3907 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 tex2html_wrap_inline3909 to be more appropriate. We like to emphasize that using a fixed value of tex2html_wrap_inline3911 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 tex2html_wrap_inline3913 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 tex2html_wrap_inline3915 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 tex2html_wrap_inline3919 0.004), values of tex2html_wrap_inline3921 tex2html_wrap_inline3923 1-2 may be more appropriate.

5.3. Comparison with PN abundances observed in the Galactic disk

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 tex2html_wrap_inline3927 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 tex2html_wrap_inline3929 for stars that ultimately can become a PN (with final core mass less than tex2html_wrap_inline39311.2 tex2html_wrap_inline3933) 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 tex2html_wrap_inline3935 tex2html_wrap_inline3937, 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 tex2html_wrap_inline3939 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).

  figure719
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. tex2html_wrap_inline3943 tex2html_wrap_inline3945. 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 tex2html_wrap_inline3949 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 tex2html_wrap_inline39510.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 tex2html_wrap_inline39530.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 tex2html_wrap_inline39550.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 tex2html_wrap_inline39570.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 tex2html_wrap_inline3959 and 1.3 tex2html_wrap_inline3961 (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 tex2html_wrap_inline3963 tex2html_wrap_inline3965. 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 tex2html_wrap_inline3967 < 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 tex2html_wrap_inline3971 to 0.9 tex2html_wrap_inline3973. In case of reduced HBB (i.e. tex2html_wrap_inline3975 tex2html_wrap_inline3977), 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 (tex2html_wrap_inline3979, tex2html_wrap_inline3981) case for which the effect of HBB is largest. Thus, as the standard model has tex2html_wrap_inline3983, 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 tex2html_wrap_inline3985. 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) tex2html_wrap_inline3987 and He/H tex2html_wrap_inline3989 (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.

5.4. Final remarks

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 tex2html_wrap_inline3993 of the Reimers law may be different for low metallicities. To simulate this effect one may want to use the yields for models with tex2html_wrap_inline3995 for Z = 0.004 and tex2html_wrap_inline3999 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 tex2html_wrap_inline4003">, tex2html_wrap_inline4005, and tex2html_wrap_inline4007">

   Table 3: AGB yields for tex2html_wrap_inline4009">, tex2html_wrap_inline4011, and tex2html_wrap_inline4013">

    Table 4: Final AGB yields for tex2html_wrap_inline4015">, tex2html_wrap_inline4017, and tex2html_wrap_inline4019">

    Table 5: Total yields for tex2html_wrap_inline4021">, tex2html_wrap_inline4023, and tex2html_wrap_inline4025">

   Table 6: Pre-AGB yields for tex2html_wrap_inline4027">, tex2html_wrap_inline4029">, and tex2html_wrap_inline4031

   Table 7: AGB yields for tex2html_wrap_inline4033">, tex2html_wrap_inline4035">, and tex2html_wrap_inline4037

   Table 8: Final AGB yields for tex2html_wrap_inline4039">, tex2html_wrap_inline4041">, and tex2html_wrap_inline4043

   Table 9: Total yields for tex2html_wrap_inline4045">, tex2html_wrap_inline4047">, and tex2html_wrap_inline4049

   Table 10: Pre-AGB yields for tex2html_wrap_inline4051">, tex2html_wrap_inline4053">, and tex2html_wrap_inline4055

   Table 11: AGB yields for tex2html_wrap_inline4057">, tex2html_wrap_inline4059, and tex2html_wrap_inline4061">

   Table 12: Final AGB yields for tex2html_wrap_inline4063">, tex2html_wrap_inline4065">, and tex2html_wrap_inline4067

   Table 13: Total yields for tex2html_wrap_inline4069">, tex2html_wrap_inline4071">, and tex2html_wrap_inline4073

   Table 14: Pre-AGB yields for tex2html_wrap_inline4075">, tex2html_wrap_inline4077">, and tex2html_wrap_inline4079

   Table 15: AGB yields for tex2html_wrap_inline4081">, tex2html_wrap_inline4083">, and tex2html_wrap_inline4085

   Table 16: Final AGB yields for tex2html_wrap_inline4087">, tex2html_wrap_inline4089">, and tex2html_wrap_inline4091

   Table 17: Total yields for tex2html_wrap_inline4093">, tex2html_wrap_inline4095">, and tex2html_wrap_inline4097

   Table 18: Pre-AGB yields for tex2html_wrap_inline4099">, tex2html_wrap_inline4101">, and tex2html_wrap_inline4103

   Table 19: AGB yields for tex2html_wrap_inline4099">, tex2html_wrap_inline4101">, and tex2html_wrap_inline4103

   Table 20: Final AGB yields for tex2html_wrap_inline4105">, tex2html_wrap_inline4107, and tex2html_wrap_inline4109">

   Table 21: Total yields for tex2html_wrap_inline4111">, tex2html_wrap_inline4113">, and tex2html_wrap_inline4115

   Table 22: AGB yields for tex2html_wrap_inline4117">, tex2html_wrap_inline4119">, and tex2html_wrap_inline4121

   Table 23: AGB yields for tex2html_wrap_inline4123">, tex2html_wrap_inline4125">, and tex2html_wrap_inline4127

   Table 24: AGB yields for tex2html_wrap_inline4129">, tex2html_wrap_inline4131">, and tex2html_wrap_inline4133

   Table 25: AGB yields for tex2html_wrap_inline4135">, tex2html_wrap_inline4137">, and tex2html_wrap_inline4139

   Table 26: AGB yields for tex2html_wrap_inline4141">, tex2html_wrap_inline4143">, and tex2html_wrap_inline4145

   Table 27: Final AGB yields for tex2html_wrap_inline4147">, tex2html_wrap_inline4149">, and tex2html_wrap_inline4151

   Table 28: AGB yields for tex2html_wrap_inline4153">, tex2html_wrap_inline4155">, and tex2html_wrap_inline4157

   Table 29: Final AGB yields for tex2html_wrap_inline4159">, tex2html_wrap_inline4161">, and tex2html_wrap_inline4163

   Table 30: AGB yields for tex2html_wrap_inline4165">, tex2html_wrap_inline4167">, and tex2html_wrap_inline4169

   Table 31: Final AGB yields for tex2html_wrap_inline4171">, tex2html_wrap_inline4173">, and tex2html_wrap_inline4175">

   Table 32: Total AGB yields for H, He for synthetic evolution model

   Table 33: AGB yields for tex2html_wrap_inline4177">, tex2html_wrap_inline4179">, and tex2html_wrap_inline4181

   Table 34: Final AGB yields for tex2html_wrap_inline4183">, tex2html_wrap_inline4185">, and tex2html_wrap_inline4187

   Table 35: Total yields for tex2html_wrap_inline4189">, tex2html_wrap_inline4191">, and tex2html_wrap_inline4193

   Table 36: AGB yields for tex2html_wrap_inline4195">, tex2html_wrap_inline4197">, and tex2html_wrap_inline4199

   Table 37: Final AGB yields for tex2html_wrap_inline4201">, tex2html_wrap_inline4203">, and tex2html_wrap_inline4205

   Table 38: Total yields for tex2html_wrap_inline4207">, tex2html_wrap_inline4209">, and tex2html_wrap_inline4211


next previous
Up: New theoretical yields

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