Free Access
Issue
A&A
Volume 557, September 2013
Article Number L17
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201322422
Published online 19 September 2013

© ESO, 2013

1. Introduction

Long thought to be composed of coeval stars born out of homogeneous material, Galactic globular clusters are now known to host at least two generations of long-lived low-mass stars with different chemical properties (e.g. Gratton et al. 2012,and references therein). First-generation stars have normal sodium abundance similar to that of halo stars. On the other hand, second-generation stars are identified thanks to their sodium overabundances. Both populations are found at the main-sequence turnoff, the subgiant branch, and the red giant branch (RGB) in all Milky Way globular clusters studied so far (e.g. Gratton et al. 2001; Lind et al. 2009; Carretta et al. 2009). At these evolution phases second-generation stars are counted to be about twice as numerous as their first-generation counterparts (Carretta et al. 2010).

However, a recent spectroscopic study (Campbell et al. 2013,hereafter C13) has revealed that the early asymptotic giant branch (AGB) of the metal-poor globular cluster NGC 6752 hosts only stars with normal and slightly enhance sodium abundance, while its RGB is populated by both sodium-normal and sodium-rich stars (Carretta et al. 2007). The result of C13 can be directly related to the deficiency of CN-strong AGB stars first discovered in NGC 6752 (Norris et al. 1981), which was subsequently reported for several other globular clusters (see Gratton et al. 2010b; Campbell et al. 2012,and references therein), although with strong cluster-to-cluster variations. It confirms the suggestion by Gratton et al. (2010b) that the most sodium-rich (and oxygen-poor) stars in globular clusters, which probably have a high helium and nitrogen content, could become AGB-manqué because of their small mass on the horizontal branch (see Sect. 3).

Here we show that the correlation between helium and sodium enhancement in the initial composition of second-generation low-mass stars that was predicted by the so-called fast-rotating massive star polluter model (FRMS) for NGC 6752 (Decressin et al. 2007a) provides a very simple and natural explanation to the different sodium distributions observed on the RGB and AGB in this cluster. Contrary to C13, who interpreted the data as a failure for all second-generation stars to climb the AGB, we conclude that among the second generation, only stars born with initial helium and sodium abundances above a cutoff value evolve directly to the white dwarf phase after central helium-burning, while those with lower helium and sodium content climb the AGB in the same manner as their first generation counterparts. Note that because of the large uncertainties on sodium abundance determination (typical error bars of 0.2 dex), it is currently not possible to observationally distinguish first- and second-generation stars below the sodium cutoff value.

2. Helium-sodium correlation in the framework of the fast rotating massive stars scenario

2.1. General framework

It is now well established that second-generation globular cluster stars formed out of the sodium-rich ashes of hydrogen burning at high temperature that were ejected by more massive first-generation globular cluster stars and were diluted with interstellar gas (Prantzos et al. 2007). Second-generation stars are therefore expected to start their life with a higher helium content than their first-generation counterparts. The impact of this helium enrichment on the extension of the horizontal branch has been extensively studied in the literature (e.g. D’Antona et al. 2002; Gratton et al. 2010a; Dotter et al. 2010), and the connection with missing AGB stars in globular clusters has also been reported (Gratton et al. 2010b).

However, the extent of the helium enrichment for a given sodium distribution strongly differs between the different polluter models. As we discuss here, this has serious implications because stars with different helium contents burn their nuclear fuel at different rates and are therefore expected to have different lifetimes.

The nature of the polluters that have long since disappeared is debated (see e.g. Prantzos & Charbonnel 2006). Possible candidates are fast-rotating massive stars (FRMS) with masses larger than 25 M (Decressin et al. 2007a), and massive (6–10 M) AGB stars (Ventura et al. 2001). In both cases the sodium-rich ejecta are enriched in helium. In the AGB polluter model, helium-enrichement is due to the second dredge-up in intermediate-mass stars (Forestini & Charbonnel 1997; Siess 2007; D’Ercole et al. 2010; Ventura et al. 2013). It occurs before the thermal pulse phase on the AGB, when the competition between third dredge-up and hot-bottom-burning is supposed to enrich the stellar ejecta in sodium. Helium-enrichement of second-generation stars is thus limited and is not directly correlated with the sodium content of the polluter yields. Consequently, the sodium spread observed in the cluster M4, which is quite similar to NGC 6752, is reproduced with limited helium enrichment almost identical for all second-generation stars (D’Ercole et al. 2012). Therefore, and as shown by C13, RGB and AGB second-generation stars are expected to cover the same range in sodium abundances in this case, which is at odds with the observations in NGC 6752. To reconcile the AGB polluter model with the observations, C13 suggested that the mass loss rates during central helium-burning on the horizontal branch and latter phases is enhanced by an ad hoc factor of about 20.

In contrast, the FRMS pollution scenario predicts broad and correlated spreads in both sodium and helium among second-generation low-mass stars (Decressin et al. 2007b). Consequently at a given age, stars with very different helium and sodium abundances are expected to populate the different evolutionary phases. In particular, the observed sodium distribution is expected to be very different between the main sequence and RGB on one hand and the AGB on the other hand, as discussed below.

thumbnail Fig. 1

Correlation predicted by the FRMS model for NGC 6752 for the initial helium and sodium abundances of second-generation stars (Decressin et al. 2007a) that is adopted in our models of low-mass stars.

2.2. Models for second-generation stars

2.2.1. Initial abundances

We computed stellar evolution models with the code STAREVOL (Siess et al. 2000; Lagarde et al. 2012) for low-mass stars with various initial masses and helium abundances with a similar metallicity as NGC 6752 ([Fe/H]  =  −1.56; Harris 1996; Gratton et al. 2005,C13). For our computations we chose a uniform helium mass fraction of 0.248 for first-generation stars (close to the Big Bang abundance deduced from cosmological parameters obtained with Planck; Planck Collaboration 2013), and a sodium abundance [Na/Fe]  =  −0.25 that corresponds to the observed value attributed to NGC 6752 first generation stars. For second-generation stars we adopted the correlation between initial helium and sodium abundances predicted by the FRMS polluter models for NGC 6752 RGB stars (Decressin et al. 2007b) that is shown in Fig. 1. This distribution accounts for the dilution of H-burning ashes ejected by FRMS with pristine gas to reproduce the abundance variations of lithium in NGC 6752 (Pasquini et al. 2005). The initial abundances of carbon, oxygen, nitrogen, magnesium, and aluminium are scaled accordingly. More details on the models will be presented in a forthcoming paper (Chantereau et al., in prep.).

2.2.2. Mass loss and other assumptions

Mass loss during the RGB was accounted for following the Reimers prescription with η = 0.5 (Reimers 1975). Additional models were computed with η = 0.65 for stars with an He abundance greater than 0.263 (see Sect. 3). We neglected atomic diffusion and rotation-induced mixing; this simplification induces an uncertainty on the stellar lifetime of typically 0.2–0.5 Gyr, which is small compared to error bars on the current age estimate of NGC 6752 (13.4 ± 1.1 Gyr; Gratton et al. 2003). The FRMS pollution scenario predicts a maximum delay of 10 Myrs between the formation of the first and second stellar generations (Krause et al. 2013), which is negligible compared to the globular cluster age. Finally, the duration of the AGB phase is typically 0.1% of the total stellar lifetime for the mass and helium domains we explore. These age uncertainties can therefore be neglected in the discussion and we can safely consider that all post-helium-burning stars in NGC 6752 are along or very nearby the 13.4 Gyr isochrone (dashed black line in Fig. 2).

thumbnail Fig. 2

Lifetime and fate of stars versus initial mass and helium content at the metallicity of NGC 6752. The colour code gives ages at the end of central helium burning, or at the RGB tip for stars that fail to ignite helium (dashed area) using the Reimers mass loss rate with η = 0.5. The 13.4 Gyr isochrone is highlighted. The white area corresponds to lifetimes longer than the age of the Universe. The grey line separates stars that end up as helium or carbon-oxygen white dwarfs. The bold black line corresponds to the AGB ascension cutoff. It separates stars that climb the AGB (right-hand side) from the AGB-manqué ones (left).

3. Model predictions

The impact of initial mass and initial helium content on the lifetime and fate of cluster stars is illustrated in Fig. 2. The stars that are currently (i.e. at 13.4 Gyr) either at the main-sequence turnoff or are climbing the RGB have sodium abundances covering the range in [Na/Fe] observed at that phase in the cluster (Carretta et al. 2007,C13) and an initial helium mass fraction between 0.248 and 0.8 (see Sect. 4 for a discussion on the fraction of stars in various bins of helium abundance).

We examine now the predictions for the abundance distribution on the AGB. Importantly, among the stars that burn helium at their centre, those with initial helium abundance above 0.31 (in mass fraction) and initial masses below 0.735 M conclude their life without climbing the AGB, as can be seen in Fig. 2. These so-called AGB-manqué (failed-AGB, Greggio & Renzini 1990) lie along the post-helium-burning 13.4 Gyr isochrone in the domain between the grey and black lines in Fig. 2. They undergo central helium-burning at very high effective temperatures on the horizontal branch (Fig. 3) on the left side of the Grundahl jump. At that phase the stellar envelope surrounding the nuclear active core is extremely thin. As a consequence, once the helium fuel is exhausted at the centre, the energy released by the contraction of the core and by hydrogen shell-burning mainly diffuses across the thin, transparent helium-rich external layers and fails to swell them up. Therefore the stars do not develop an extensive outer convection zone and do not move towards the AGB along the Hayashi track. Instead, they first evolve at a constant effective temperature towards higher luminosity before moving towards the white dwarf region (Fig. 3).

thumbnail Fig. 3

Theoretical evolution path in the Hertzsprung-Russell diagram for stars in NGC 6752. Along individual tracks the solid line corresponds to the central helium-burning phase, and the dotted line to the more advanced phases (AGB or AGB-manqué). All models are computed with a value of η = 0.5 in the Reimers mass loss prescription, and we also show the prediction with η = 0.65 for the 0.79 M, Y = 0.263 model (long dashed). They correspond to a first-generation star (black) and to several second-generation stars (colours) that terminate helium burning at approximately 13.4 Gyr. The corresponding initial masses and helium abundances are given in the graph. The 0.82 M star born with pristine helium as well as second-generation stars born with modest helium abundances evolve along the horizontal branch at a relatively low temperature and then all converge towards the AGB. Instead second-generation stars born with higher helium abundances evolve at much higher effective temperatures on the horizontal branch and become AGB-manqué.

We note that the AGB ascension cutoff values for the initial helium abundance and stellar mass are slightly sensitive to the mass loss rate adopted along the RGB. With the assumed FRMS helium-sodium distribution and using a value of 0.5 for the mass loss rate parameter η, the theoretical AGB ascension cutoff is obtained for a helium mass fraction of 0.31 which corresponds to an initial [Na/Fe] value of 0.4 dex. Thus the highest sodium abundance for second-generation stars that move towards the AGB is 0.4 dex in that case. This is higher than the upper envelope value of 0.18 dex observed for AGB stars in NGC 6752. On the other hand, using a value of 0.65 (instead of 0.5) is sufficient to shift the theoretical AGB ascension cutoff along the 13.4 Gyr isochrone down to the observed [Na/Fe] cutoff value of 0.18. This can be seen in Fig. 3, where we show the evolution paths for both values of η for the star that lies at the observational [Na/Fe] cutoff of 0.18 dex on the 13.4 Gyr isochrone (i.e., initial mass and helium mass fraction of 0.79 M and 0.263). Depending on whether η is 0.5 or 0.65, the model either moves towards the AGB or becomes an AGB-manqué. Therefore we found that within the FRMS framework, a modest increase in the mass loss along the RGB for the most helium-rich stars, compatible with modified Reimers-type prescriptions (Schröder & Cuntz 2007) that propose a stronger dependency of the mass-loss rate with stellar effective temperature than the classical Reimers formula, is sufficient to achieve a perfect agreement between the theoretical and observational sodium cutoff on the AGB. This is different from the attempt of C13 to fit the sodium data by assuming an increased mass loss rate by a factor of 20 on the horizontal branch for all second generation stars.

It is interesting to estimate the fraction of horizontal branch stars that evolve along the AGB depending on the value of η. Typically, for η = 0.5 (a limiting value of He = 0.31), models show that ~85% climb the AGB, while a modest increase from η = 0.5 to 0.65 reduces this fraction to 25%, which agrees reasonably well with the 30% given by C13.

4. Discussion and conclusions

We showed that the correlation between the helium and sodium abundances expected from the FRMS models by Decressin et al. (2007a) for NGC 6752 provides a reasonable solution for the observed lack of sodium-rich stars along the early-AGB in this globular cluster. The FRMS scenario predicts a high initial dispersion of helium abundances, and one may wonder whether this requirement agrees with the distributions of stars in the observed colour magnitude diagram. For NGC 6752 a recent study by Milone et al. (2013,hereafter M13) obtained that the observed colour dispersions at a given magnitude on the main sequence and on the RGB can be well fitted by a helium dispersion of only 0.03–0.04 in mass fraction, much lower than the helium dispersion predicted in our model. However, a higher helium dispersion than the one inferred by M13 for NGC 6752 can probably not be ruled out for the following reasons: first, the chance of observing stars of a given initial helium abundance at a given magnitude and age decreases strongly for high helium values. Typically, starting from the initial distribution in helium given by Decressin et al. (2007a), which reproduces the observed ratio 1:2 between the first and second generation of stars, we have the following distribution of the number of stars as function of the helium abundance: 82.4% with a helium lower than 0.275, 12.6% with a helium between 0.275 and 0.4, and 5% with a helium between 0.4 and 0.7. These numbers are given for an age of 13.4 Gyr and are estimated at a position corresponding to two magnitudes below the turn off. Second, the dispersion obtained by M13 was deduced from shifts in colours of averaged positions of various main sequence and RGB lines, but there are stars that are bluer than the bluest averaged lines used by M13. This may indicate that there are stars with higher helium abundances than inferred by these authors. Third, the link between colours and helium and/or other element abundances is by far not straightforward, since it requires computating model atmospheres and synthetic spectra that probably have their own uncertainties. Moreover, M13 explored only a relatively narrow set of parameters, and from their study one cannot exclude that higher helium abundances might be compatible with the observations. Thus, until comparisons become available for various atmosphere models and various codes that compute synthetic spectra for a large range of chemical compositions, we cannot discard the possibility that stars with higher helium abundances than those inferred by M13 are present.

According to C13, all stars bluer than the Grundahl jump in NGC 6752 (i.e. bluer than 4.1–4.3 in log Teff) do not ascend the AGB. From Fig. 3 we see that the limit in effective temperature between stars evolving into the AGB phase and those skipping this phase is between 4.1 (for η equal to 0.65 instead of 0.5 for helium-rich stars) and 4.3 (for models with a constant value of η). Therefore our models are quite consistent with the limit in effective temperature suggested by C13. However, as noted above, this limit depends on the calibration process that allows translating colours into effective temperatures and thus may present some uncertainties depending on which atmosphere and synthetic spectrum model is used.

Although the observed cutoff value [Na/Fe]  ~ 0.18 dex is somewhat uncertain (compare with the typical error bars of 0.2 dex usually quoted in the literature for sodium abundance determination), and in view on the uncertainties on the mass loss rate for evolved low-mass stars, we consider the agreement between the observations and the predictions of the current models very satisfactory.

In conclusion, our results, obtained without fine-tuning the initial helium and sodium distribution or the physics of the stellar models, satisfactorily reconcile the different sodium patterns along the RGB and AGB of NGC 6752, and naturally explain the absence of AGB stars with high sodium abundances in this globular cluster. This indicates that the FRMS represents a reasonable solution for the origin of the chemical anomalies in globular clusters. As outlined above, some additional efforts are needed to unambiguously constrain the actual dispersion in helium abundances within these systems.

Acknowledgments

We acknowledge support from the Swiss National Science Foundation, ESF-Eurogenesis, and the French Centre National de la Recherche Scientifique. We are grateful to the referee, A.Weiss, for his constructive comments and to A. Maeder for careful reading of the manuscript.

References

  1. Campbell, S. W., Yong, D., Wylie-de Boer, E. C., et al. 2012, ASPCS, 458, 205 [NASA ADS] [Google Scholar]
  2. Campbell, S. W., D’Orazi, V., Yong, D., et al. 2013, Nature, 498, 198 (C13) [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  3. Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2007, A&A, 464, 967 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Carretta, E., Bragaglia, A., Gratton, R., & Lucatello, S. 2009, A&A, 505, 139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2010, A&A, 516, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Chantereau, W., Charbonnel, C., Decressin, T., & Meynet, G. 2013, Mem. Soc. Astron. It., 84, 200 [NASA ADS] [Google Scholar]
  7. D’Antona, F., Caloi, V., Montalbán, J., Ventura, P., & Gratton, R. 2002, A&A, 395, 69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Decressin, T., Meynet, G., Charbonnel, C., Prantzos, N., & Ekström, S. 2007a, A&A, 464, 1029 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Decressin, T., Charbonnel, C., & Meynet, G. 2007b, A&A, 475, 859 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. D’Ercole, A., D’Antona, F., Ventura, P., Vesperini, E., & McMillan, S. L. W. 2010, MNRAS, 407, 854 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  11. D’Ercole, A., D’Antona, F., Carini, R., Vesperini, E., & Ventura, P. 2012, MNRAS, 423, 1521 [NASA ADS] [CrossRef] [Google Scholar]
  12. Dotter, A., Sarajedini, A., Anderson, J., et al. 2010, ApJ, 708, 698 [NASA ADS] [CrossRef] [Google Scholar]
  13. Forestini, M., & Charbonnel, C. 1997, A&A, 123, 241 [Google Scholar]
  14. Gratton, R. G., Bonifacio, P., Bragaglia, A., et al. 2001, A&A, 369, 87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Gratton, R. G., Bragaglia, A., Carretta, E., et al. 2003, A&A, 408, 529 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Gratton, R. G., Bragaglia, A., Carretta, E., et al. 2005, A&A, 440, 901 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Gratton, R. G., Carretta, E., Bragaglia, A., Lucatello, S., & D’Orazi, V. 2010a, A&A, 517, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Gratton, R. G., D’Orazi, V., Bragaglia, A., Carretta, E., & Lucatello, S. 2010b, A&A, 522, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Gratton, R. G., Carretta, E., & Bragaglia, A. 2012, A&ARv, 20, 50 [CrossRef] [Google Scholar]
  20. Greggio, L., & Renzini, A. 1990, ApJ, 364, 35 [Google Scholar]
  21. Harris, W. E. 1996, AJ, 112, 1487 [NASA ADS] [CrossRef] [Google Scholar]
  22. Krause, M., Charbonnel, C., Decressin, T., Meynet, G., & Prantzos, N. 2013, A&A, 552, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Lagarde, N., Decressin, T., Charbonnel, C., et al. 2012, A&A, 543, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Lind, K., Primas, F., Charbonnel, C., Grundahl, F., & Asplund, M. 2009, A&A, 403, 545 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Milone, A. P., Marino, A. F., Piotto, G., et al. 2013, ApJ, 767, 120 (M13) [NASA ADS] [CrossRef] [Google Scholar]
  26. Norris, J., Cottrell, P. L., Freeman, K. C., & Da Costa, G. S. 1981, ApJ, 244, 205 [NASA ADS] [CrossRef] [Google Scholar]
  27. Paresce, F., & De Marchi, G. 2000, ApJ, 534, 870 [NASA ADS] [CrossRef] [Google Scholar]
  28. Pasquini, L., Bonifacio, P., Molaro, P., et al. 2005, A&A, 441, 549 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Planck Collaboration 2013, A&A, submitted [arXiv:1303.5076] [Google Scholar]
  30. Prantzos, N., & Charbonnel, C. 2006, A&A, 458, 135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Prantzos, N., Charbonnel, C., & Iliadis, C. 2007, A&A, 470, 179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Reimers, D. 1975, Mem. Soc. Roy. Sci. Liège, 8, 369 [Google Scholar]
  33. Siess, L. 2007, A&A, 476, 893 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 [NASA ADS] [Google Scholar]
  35. Schröder, K.-P., & Cuntz, M. 2007, A&A, 465, 593 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Ventura, P., D’Antona, F., Mazzitelli, I., & Gratton, R. 2001, ApJ, 550, L65 [NASA ADS] [CrossRef] [Google Scholar]
  37. Ventura, P., Di Criscienzo, M., Carini, R., & D’Antona, F. 2013, MNRAS, 431, 3642 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1

Correlation predicted by the FRMS model for NGC 6752 for the initial helium and sodium abundances of second-generation stars (Decressin et al. 2007a) that is adopted in our models of low-mass stars.

In the text
thumbnail Fig. 2

Lifetime and fate of stars versus initial mass and helium content at the metallicity of NGC 6752. The colour code gives ages at the end of central helium burning, or at the RGB tip for stars that fail to ignite helium (dashed area) using the Reimers mass loss rate with η = 0.5. The 13.4 Gyr isochrone is highlighted. The white area corresponds to lifetimes longer than the age of the Universe. The grey line separates stars that end up as helium or carbon-oxygen white dwarfs. The bold black line corresponds to the AGB ascension cutoff. It separates stars that climb the AGB (right-hand side) from the AGB-manqué ones (left).

In the text
thumbnail Fig. 3

Theoretical evolution path in the Hertzsprung-Russell diagram for stars in NGC 6752. Along individual tracks the solid line corresponds to the central helium-burning phase, and the dotted line to the more advanced phases (AGB or AGB-manqué). All models are computed with a value of η = 0.5 in the Reimers mass loss prescription, and we also show the prediction with η = 0.65 for the 0.79 M, Y = 0.263 model (long dashed). They correspond to a first-generation star (black) and to several second-generation stars (colours) that terminate helium burning at approximately 13.4 Gyr. The corresponding initial masses and helium abundances are given in the graph. The 0.82 M star born with pristine helium as well as second-generation stars born with modest helium abundances evolve along the horizontal branch at a relatively low temperature and then all converge towards the AGB. Instead second-generation stars born with higher helium abundances evolve at much higher effective temperatures on the horizontal branch and become AGB-manqué.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.