Free Access
Issue
A&A
Volume 571, November 2014
Article Number A81
Number of page(s) 8
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201424540
Published online 14 November 2014

© ESO, 2014

1. Introduction

In the last decade, the classical view of Galactic globular clusters (GGCs) as simple (single-age, single initial chemical composition) stellar populations has changed, thanks to a large body of photometric and spectroscopic observations (see e.g. Gratton et al. 2012; Piotto et al. 2012, and references therein for reviews). We now know that probably all GGCs host multiple stellar populations, each one with its own chemical abundance pattern. A typical GGC is populated by a stellar component born with the standard α-enhanced metal distribution typical of the halo field population, plus additional sub-populations showing light element (anti-)correlations, that is, a range of C and O, and sometimes Mg depletions, together with N and Na, and sometimes Al enhancements. These additional sub-populations are also characterised by helium enhancements (see e.g. Milone et al. 2014). Depending on the photometric filters (see e.g. Sbordone et al. 2011; Cassisi et al. 2013), the cluster sub-populations may lie on separate sequences in the colour-magnitude diagrams (CMDs) as a result of the effect of the multiple abundance distributions on isochrones and/or bolometric corrections. This broad picture is nowadays denoted as the globular cluster multiple population phenomenon.

The main scenarios devised to explain the formation of sub-populations within individual clusters (see e.g. Decressin et al. 2007; de Mink et al. 2009; D’Ercole et al. 2011; Conroy & Spergel 2011; Valcarce & Catelan 2011) assume that the chemical patterns are produced by multiple star formation episodes during the early stages of the cluster evolution. After the formation of a first generation (FG) of stars with metal abundance ratios (and He abundance) typical of the halo field population, successive generations (SG stars) originate from matter ejected by pre-existing FG stars (polluters; massive asymptotic giant branch stars or massive stars) diluted with material of FG composition not yet involved in star formation episodes. A different scenario proposed recently by Bastian et al. (2013) assumes instead that stars with SG composition are actually FG low mass stars polluted during their fully convective pre-main sequence phase (pre-MS; hence their chemical composition at the beginning of the MS is essentially uniform, like the multiple star formation scenario) by the ejecta of short-lived massive stars. Throughout the paper we denote as FG stars the cluster sub-population with metal mixture and He abundances typical of the Galactic halo field population, and as SG stars the other populations, irrespective of their origin.

Tracing the distribution of SG stars along the observed cluster CMDs is very important in order to study the evolutionary properties of SG stars, and to shed more light on their formation and the nature of the polluters.

One striking feature of SG stars is, at least in some clusters, the difference between the number ratios of SG-to-FG objects along the asymptotic giant branch (AGB) and the red giant branch (RGB). While stars with very strong CN bands (a signature, together with high-Na and low-O abundances, of SG composition) are abundant along the RGB of any GGC, the majority of AGB stars in these clusters usually seem to be FG objects. This early discovery by Norris et al. (1981) has been confirmed by more recent investigations (Campbell et al. 2010; Gratton et al. 2010). Using spectroscopic measurements of Na, Campbell et al. (2013) have found no SG stars in a sample of AGB stars hosted by NGC 6752. On the other hand, in clusters such as M 5 (Smith & Norris 1993) and possibly 47 Tuc (Campbell et al. 2006), CN-strong/Na-rich AGB stars have been observed, and the discrepancy with the ratios observed along the RGB is probably not large. It is evident that additional observational effort has to be devoted to this issue.

Recently, an attempt to connect the lower incidence of SG stars along the AGB to the multiple population phenomenon has been made by Gratton et al. (2010; but see also the pioneering works by Norris et al. 1981; Norris 1981). These authors found a correlation between the relative frequency of AGB stars and the minimum mass of stars along the horizontal branch (HB) in a given cluster (that is a proxy of the maximum effective temperature along the HB, and hence of the HB morphology), with a further dependence on the cluster metallicity. Clusters with the lowest number ratio of AGB to HB stars (the so-called R2 parameter, Caputo et al. 1989) are those with the bluest HB morphology. In fact, objects with extended HB blue tails like NGC 2808 and NGC 2419 have the lowest R2 values, e.g. R2 ~ 0.06, to be compared with typical values R2 ~ 0.12 for redder morphologies (see e.g. Sandquist 2000), and the very high R2 = 0.18 ± 0.02 for M5 (Sandquist & Bolte 2004).

To connect this result to the cluster multiple populations, we note that moving towards the blue side of the HB (e.g. towards lower masses) when the evolving mass drops below ~0.6 M (the exact value depending on the initial chemical composition) the ratio of HB to AGB timescales increases steadily, until reaching a minimum HB mass (equal to about 0.5 M) below which stars do not evolve to the AGB phase, the so-called AGB-manqué objects (Greggio & Renzini 1990; Dorman et al. 1993). After central-He exhaustion these stars move to the white dwarf cooling sequence. Second generation stars are indeed expected to have a lower mass along the HB because they are typically He-enhanced, hence they originate from RGB stars with a lower mass (at a given age, He-enhanced stars have shorter main sequence lifetimes) and end up with, on average, a lower mass along the HB if the RGB mass loss is approximately the same for FG and SG sub-populations (see e.g. D’Antona et al. 2002; Dalessandro et al. 2011, 2013; Gratton et al. 2011; Marino et al. 2011, 2014, and references therein for a more detailed explanation). One expects, therefore, that the AGB of clusters with a blue HB should lack at least part of the SG component, compared to what is seen along the previous RGB phase.

In this paper we study in detail the case of the AGB population of NGC 6752, in light of the recent results by Campbell et al. (2013), who performed an accurate spectroscopic analysis of both RGB and AGB stars, determining Na abundances. They also stated that Na abundance estimates in cool giants should be more robust than other light elements such as C and N, because they are less affected by molecular band formation and surface gravity estimates, and/or by the possibility of in situ variations due to evolutionary process(es) (see e.g. Norris et al. 1981; Salaris et al. 2002; Denissenkov et al. 2014, and references therein). Campbell et al. (2013) found that all 20 AGB stars in their sample show Na abundances typical of FG stars. To explain the lack of SG AGB stars and the observed low value of R2 (R2 ~ 0.06), they envisaged the following scenario. All FG stars populate the redder side of the observed HB, and all objects hotter than ~11500 K –the so-called Grundahl’s jump (Grundahl et al. 1999), corresponding to the observed effective temperature of the onset of radiative levitation for HB stars– miss the AGB phase, even though standard stellar evolution predicts that at the cluster metallicity only stars evolving on the HB at Teff larger than ~24000 K skip this phase (Dorman et al. 1993; Pietrinferni et al. 2004).

The implication is that objects in the Teff range between ~11 500 K and at least ~24 000 K, must have lost mass very efficiently during the HB phase, to avoid the AGB phase. We note that the effect of an efficient mass loss during the HB evolution has been already explored in previous analyses (see e.g. Yong et al. 2000, and references therein). By employing the Reimers (Reimers 1975) mass loss law, Campbell et al. (2013) estimated a value for the free parameter η ~10 to explain the lack of SG stars along the AGB, an enhancement of a factor of ~20 compared to the standard value used in the calculation of RGB models. Such an enhancement compared to standard values of η is tentatively attributed to the metal enhancement in the envelopes of these stars, caused by radiative levitation (see e.g. Michaud et al. 2008). If the evolution of NGC 6752 HB stars is typical of GGCs with long blue HB tails, this scenario implies that a cluster like M 13, with a very similar HB morphology, [Fe/H], and age, should show the same lack of SG stars along the AGB.

In this analysis we first revisit the HB mass loss scenario discussed by Campbell et al. (2013), considering additional constraints from independent theoretical considerations and observations (Sect. 2). As a second step, we perform synthetic HB modelling, using canonical HB models without mass loss during the core He-burning stage to establish more accurately whether the observed low value of R2 ~ 0.06 does really require very efficient mass loss along the HB. We then study the AGB population in our synthetic samples, to check its expected composition in terms of FG and SG components (Sect. 3). We conclude with a discussion of the results and a comparison with the abundances measured in M 13, a cluster whose properties are very similar to those of NGC 6752 (Sect. 4).

2. Mass loss on the hot HB

In the scenario envisaged by Campbell et al. (2013), all NGC 6752 HB stars hotter than Teff ~ 11 500 K fail to reach the AGB because of enhanced mass loss, possibly associated to the surface metal enhancement caused by radiative levitation. This temperature is assumed to mark the transition between FG and SG stars along the HB, consistent with spectroscopic measurements (Villanova et al. 2009) of nine HB stars in the Teff range between ~ 8500 K and ~ 8800 K, that show how all objects but one have [Na/Fe] ratios typical of FG stars, and the average He abundance is consistent also with the FG component.

thumbnail Fig. 1

Observed distributions of [Na/Fe] abundance ratios in NGC 6752 RGB stars as determined by Campbell et al. (2013) and Carretta et al. (2007) – dashed and solid line, respectively – and in AGB stars (Campbell et al. 2013, dash-dotted line – see text for details).

Open with DEXTER

As discussed in the introduction, this enhanced mass loss efficiency is invoked by Campbell et al. (2013) to explain both the observed low value of R2, and the observed [Na/Fe] distribution along the AGB when compared to the RGB counterpart, as can be seen in Fig. 1. In the same figure we also show the [Na/Fe] RGB distribution determined by Carretta et al. (2007) with a larger sample of stars. A Kolmogorov-Smirnov test shows that the two RGB abundance distributions determined by Campbell et al. (2013) and Carretta et al. (2007) are statistically consistent. The three peaks in the histogram of the RGB abundances correspond to [Na/Fe] ~ − 0.05 (FG), ~+0.35 and ~+0.60 (both SG), respectively, and can be nicely associated on a one-to-one basis to the components of the cluster triple MS discovered by Milone et al. (2013). It is very clear that the AGB abundances match the distribution of the lowest peak of the RGB Na distribution, corresponding to the FG composition.

thumbnail Fig. 2

Right panel: Hertzsprung-Russell (H-R) diagrams of stellar models with initial mass equal to 0.59 M, evolved from the ZAHB using the Reimers’ mass loss law with various values of the free parameter η (see labels). Left panel: as the right panel but for a model with initial mass equal to 0.55 M.

Open with DEXTER

Regarding the overall fraction and He abundance of SG stars in the cluster, the photometric analysis by Milone et al. (2013) revealed the presence of a triple MS with a normal cosmological Y ~ 0.245 (the FG component), Y ~ 0.254 and ~0.270.28 respectively, these latter two sequences being identified as the SG component. The number ratios among the three populations, that could also be traced along sub-giant branch and RGB are approximately 30:40:30, which correspond to a ~70% fraction of SG stars. This fraction is consistent with the RGB distribution of Na, shown in Fig. 1. Milone et al. (2013) do not find any evidence of radial gradients of the population ratios out to about 2.5 half-mass radii.

Our first test was to revisit the Campbell et al. (2013) results about the required efficiency of the Reimers mass loss law to deplete the AGB population of NGC 6752. To this purpose, we have calculated α-enhanced ([α/ Fe] = 0.4) HB models with Y = 0.273, Z = 0.001 –corresponding to [Fe/H] ~ − 1.6, consistent with the spectroscopic estimate [Fe/H] = −1.56 reported by Carretta et al. (2010), and He-abundance appropriate for the more He-rich SG population of NGC 6752– for a wide range of masses, with the same code, input physics, metal mixture and procedures of the BaSTI database (Pietrinferni et al. 2006)1. From the BaSTI database we have considered a set of α-enhanced HB stellar models for Y = 0.246, Z = 0.001, the counterpart of the FG cluster composition, which will be used later on in the next section2.

From the calculations with Y = 0.273 we singled out the track with mass equal to 0.59 M, whose Zero Age HB (ZAHB) Teff is ~ 11500 K – corresponding to the Grundahl’s jump – and a hotter one of 0.55 M, with ZAHB Teff ~ 16000 K. For these two masses we computed additional models considering the Reimers mass loss law with several choices of the free parameter η. The Hertzsprung-Russell (H-R) diagrams of the resulting evolutionary tracks are shown in Fig. 2.

The 0.59 M tracks reach the AGB for η increasing up to ~2.0. This model climbs the AGB to a luminosity log (L/L) ≈ 2.7, before moving towards the hot side of the H-R diagram. The track with η = 4 barely reaches the AGB stage, whereas for η ≥ 10.0 the track behaves as an AGB-manqué star. Typical mass loss rates for the tracks with η = 10 are a few times 10-9 M yr-1.

The hotter 0.55 M models behaves as an AGB-manqué star for η ≥ 4. In this case the critical mass loss rates are also a few times 10-9 M yr-1. These threshold values of η are a factor of 1025 higher than the typical values adopted for RGB evolutionary calculations with mass loss (see i.e. Pietrinferni et al. 2004, 2006; Campbell et al. 2013). In fact, employing such high values of η would prevent RGB models from experiencing the He-flash (Castellani & Castellani 1993).

Our results are qualitatively in agreement with those obtained by Campbell et al. (2013), even though they considered a more massive HB model with a higher Y. These mass loss rates are however much larger than constraints from the more physically motivated mass loss law by Vink & Cassisi (2002). These authors computed HB mass-loss rates under the hypothesis that radiation pressure on spectral lines drives a stellar wind during the He-burning stage; their mass loss law3 is expected to be much more appropriate for HB models than the Reimers formula, and is applicable in the range between 12 500 K and 35 000 K. The predicted rates are ~23 orders of magnitude lower than required. Even considering solar surface abundances for the HB models, to mimic the effect of radiative levitation4 we found that the rates predicted by the Vink & Cassisi (2002) equation must be increased by a factor of ~2050 (depending on the mass of the HB model) in order to attain the mass loss efficiency required to force a HB star to miss the AGB stage.

There are also some independent constraints on the efficiency of mass loss in HB stars, albeit for very hot objects, in a Teff range where models at constant mass already skip the AGB. The constraints arise because mass loss affects the efficiency of atomic diffusion and radiative levitation, and excessive mass loss rates would produce surface abundances inconsistent with observations. Michaud et al. (2011) considered a HB model with Teff ≈ 30 000 K and determined that mass loss rates compatible with the observed abundances for these stars cannot be larger than (3 − 5) × 10-14 M yr-1. Another independent constraint on the mass loss efficiency in hot HB stars come from the evidence that some B subdwarf stars undergo non radial pulsations. A recent analysis (see Hu et al. 2011, and references therein) shows that, if mass loss is efficient in these objects, the rates cannot by larger than ~ 10-15 M/ yr, otherwise the driving mechanism of the pulsations (an opacity-driven mechanism enabled by Fe that accumulate diffusively in the stellar envelope) would not work. Assuming the mass loss scaling (rather than the absolute values) along the HB predicted by Vink & Cassisi (2002) is correct, the expected variation between temperatures of the order of 30 000 K and ~12 000 K is within a factor of 23. In conclusion, the huge mass loss rates required to force all HB stars with a ZAHB location hotter than Teff ≈ 11 500 K to miss the AGB stage do not seem to be substantiated by the previous theoretical and observational considerations. Additional work on both the theoretical and observational side has to be done in order to settle this issue.

thumbnail Fig. 3

Upper panel: CMD of NGC 6752 HB stars from Momany et al. (2002). Lower panel: one synthetic realisation of the cluster HB, with the same number of stars as observed. Open circles denote FG stars, filled triangles SG stars with Y = 0.254, open squares SG stars with Y = 0.273. The vertical lines define the boundaries of the three boxes in which the observed CMD has been partitioned. The location of the reddest vertical line coincides with the location of the Grundahl’s jump (see text for details).

Open with DEXTER

3. Synthetic HB simulations

With the goal to determine whether the observed value of R2 is consistent with a large fraction of HB stars missing the AGB, we have performed an analysis based on synthetic HB calculations (see e.g. Rood 1973; Catelan 1993, for details about this technique), made with the code described by Dalessandro et al. (2013). As an observational counterpart we have considered the UV CMD by Momany et al. (2002), whose HB (cleared of the field contamination) is shown in Fig. 3. The shape of the cluster blue HB in the U-(U-V) CMD is more sensitive to the He-abundance than classical optical CMDs. Concerning the theoretical models, we have employed BaSTI α-enhanced HB tracks for Z = 0.001, with Y = 0.246 and Y = 0.273, described in the previous section. We are assuming that the age of SG and FG stars is the same, equal to 12.5 Gyr. This assumption is appropriate because studies of multiple populations in GGCs have shown that age differences between FG and SG stars is at most of the order of ~100 Myr (see e.g. Cassisi et al. 2008; Marino et al. 2012, and references therein). The precise value of this common age is also not crucial, because a different absolute age would simply cause a variation of the RGB mass loss required to reproduce the observed mass distribution along the HB. The input physics of BaSTI models is detailed in Pietrinferni et al. (2004) and references therein. Crucial to the prediction of HB and AGB evolutionary timescales is the treatment of core mixing during central He-burning. Our models include semiconvection with breathing pulses suppressed (see Pietrinferni et al. 2004, and references therein). We also note that very recently Gabriel et al. (2014) have shown that the inclusion of semiconvection in the calculation of HB models (that extends the He-burning lifetime and the final He-depleted core mass) seems to be an adequate approach to treat core mixing during central He-burning. As far it concerns the breathing pulses (a mixing instability occurring at the boundary of the convective core when central He is substantially depleted), the analysis by Cassisi et al. (2001) shows that their inclusion would smooth out the AGB clump, at odds with the observations, that show very concentrated and prominent AGB clumps in well populated photometries of GGCs (see e.g. M 3 and 47 Tuc photometries by Buonanno et al. 1994; Beccari et al. 2006, respectively, as just two examples). As discussed by Cassisi et al. (2003), the different techniques available to suppress this instability do not affect the predicted value of R2. In general, the inclusion of breathing pulses would decrease the predicted value of R2 by ~30% for models at the red side of the HB.

Once the mixing treatment is fixed, the main source of uncertainty in the predicted R2 values is given by the current uncertainty in the , that affect this parameter at the level of ±0.01.

Detailed analyses of the uncertainties affecting He-burning, low-mass stellar models can be found in Cassisi et al. (1998, 1999, 2001, 2003, 2007) and reference therein.

As described before, we considered HB models with two values of Y; synthetic stars with intermediate values of the initial He mass fraction were obtained by interpolation during the synthetic HB calculations. The BaSTI HB tracks also include the AGB evolution until the first thermal pulse (TP); given that in low mass stars the TP lifetime is largely negligible compared to the AGB timescales before the first TP (see e.g. the recent AGB calculations by Weiss & Ferguson 2009, for the case of 1.0 M models at various metallicities), our simulations can predict consistently the value of R2 for a given HB morphology. Second generation stars with enhanced He are represented by models with standard FG chemical composition and the appropriate He abundance.

To account for the effect of radiative levitation in the stellar atmospheres, which affects stars hotter than the Grundahl’s jump, we have applied to the bolometric luminosities predicted by the models, bolometric corrections to the U and V photometric bands for [Fe/H] = 0.0 (and scaled-solar mixture) when Teff is above ~12 000 K5. This is a crude approximation, because of the lack of extended grids of both stellar evolution and atmosphere models with a wide range of chemical compositions, that include consistently the effect of radiative levitation during the whole HB evolution. However, the use of this approximation in synthetic HB modelling has allowed Dalessandro et al. (2011) to recover the different magnitude levels of NGC 2808 sub-populations with the three He-abundances inferred from the multimodal main sequence, in optical-ultraviolet CMDs of the HB. In addition, the Teff distribution of NGC 2808 extreme HB stars determined from the simulations was found to be in good agreement with the spectroscopic estimates by Moehler et al. (2004) As shown by Fig. 3, we divided the observed CMD of the cluster HB into three boxes, with −0.3 < (UV) ≤ 0.4 (red box) , −0.95 < (UV) ≤ − 0.3 (intermediate box), and (UV) ≤ − 0.95 (blue box). The blue edge of the red box corresponds approximately to the colour of the Grundahl’s jump, while the blue edge of the intermediate box is chosen to have approximately 30:40:30 number ratios of HB stars in the three regions (red:intermediate:blue), the same as the estimated ratios of the three cluster sub-populations (in order of increasing Y) along MS and RGB, that we use as input.

In our modelling we assumed that stars belonging to the three sub-populations are fed at a continuous rate onto the HB (see e.g. Rood 1973), and that the FG component is located in the red box, consistent with spectroscopic observations. The slope of the observed HB in the intermediate box is sensitive to Y, and is consistent with the location of synthetic stars from the SG component with Y = 0.254. As a consequence, the most He-enriched stars must be largely confined to the blue box. With these constraints on the simulations, the relative number of stars in the different boxes, plus the magnitude and colour location of the HB were reproduced with E(BV) = 0.03 (consistent with E(BV) = 0.04 ± 0.005 and systematic error of ±0.011, as estimated by Gratton et al. 2003), (mM)0 = 13.12 (consistent i.e. with the cluster white dwarf distance by Renzini et al. 1996) and the mass distribution shown in Fig. 4. The simulations also include a Gaussian random photometric error whose mean value for the various filters was taken from Momany et al. (2002) photometry.

Figure 3 shows one realisation of the HB, with the number of synthetic stars (275) equal to the sample of stars in the Momany et al. (2002) photometry. As discussed before, empirical and theoretical constraints dictate that the three components should mainly populate separate regions, and in this synthetic CMD only evolved SG stars cross the boundaries between boxes, causing just a minor contamination of the neighbouring population.

Assuming a common age of 12.5 Gyr, the mass distribution of Fig. 4 corresponds to the following total RGB mass loss (ΔM) for the FG and SG populations:

  • 1.

    FG – ΔM ranging from 0.180 to 0.190 M with a flat probability distribution;

  • 2.

    SG, Y = 0.254ΔM = (0.220 ± 0.014) M with a Gaussian probability distribution;

  • 3.

    SG, Y = 0.273ΔM = (0.282 ± 0.002) M Gaussian for 80% of stars, and ΔM- 0.268±0.006 for 20% of stars.

Small variations of this age and/or possible small age differences between FG and SG stars will affect the ΔM values necessary to reproduce the mass distribution of Fig. 4. For example, an age change of 0.5 Gyr at around 12.5 Gyr changes the evolving RGB mass by ~0.01 M.

thumbnail Fig. 4

Mass distribution along the synthetic HB of Fig. 3. The meaning of the symbols is as in Fig. 3.

Open with DEXTER

It is obvious from Fig. 3 that an exact match of the observed colour (hence mass) distribution of HB stars within each box would require (especially in the intermediate and blue box) more complex mass loss laws. There is no a priori reason why a simple flat or Gaussian probability distribution should represent the correct mathematical form of the total RGB mass loss, but on the other hand they provide an approximation that is sufficient for most purposes. The value of the parameter R2 we obtained from our simulations is basically unaffected (changes of less than 0.01) by modifying the mass distributions within each box, i.e. slightly varying the mean value and/or the spread of the Gaussian distributions, or switching from Gaussian to flat distribution, while keeping the number of objects within each box unchanged.

Simulations with a large total number of HB stars (about 50 times the observed number in our adopted photometry) provide R2 = 0.07, in very good agreement with the observed R2 = 0.06 ± 0.02 (see data in Table 1 of Sandquist 2000). When considering multiple realisations of the HB with the same number of stars in Momany et al. (2002), we obtain statistical variations of R2 – as a result of the small data set – of about ±0.01 around R2 = 0.07. Coincidentally, for the number of HB stars in our adopted CMD, one expects a number of AGB stars approximately equal to the number sampled spectroscopically by Campbell et al. (2013).

This agreement with our models at constant mass does not coincide with the possibility that a large number of stars miss the AGB phase because of mass loss during the core He-burning stage, apart from the bluest object that naturally skip the AGB. In the hypothesis that all SG stars along the HB miss the AGB phase, our simulations predict a value R2 = 0.03 that is too low. The agreement with the observed value of R2 comes instead from the decreasing ratio between AGB and HB lifetimes when moving towards the blue side of the HB. This is shown clearly by Fig. 5, which shows the ratio of AGB to HB lifetimes for the models used in our simulations. We note the steep drop below ~0.560.58 M; for the sake of comparison the upper mass limit in the synthetic CMD is ~0.6 M, see Fig. 4. In this region the lifetime ratio is independent of the composition. Helium-rich models have slightly higher ratios above this mass threshold, because of a more efficient H-burning shell during the HB phase.

thumbnail Fig. 5

Ratio of AGB to HB lifetimes as a function of the total mass for the models used in our simulations (see text for details).

Open with DEXTER

As a test, we have calculated a well-populated synthetic HB with the parameters listed above, considering just the FG component. We obtained in this case R2 = 0.11. The same test considering only the SG population with Y = 0.254 provided R2 = 0.09, while the case with just the SG population with Y = 0.273 gave R2 = 0.02. The progression towards decreasing values of R2 when considering increasingly bluer morphologies is clear. It is important to note that all stars belonging to the SG population with Y = 0.254 are expected to evolve to the AGB. All stars predicted to miss the AGB belong to the SG population with Y = 0.273; however, a very small fraction of this population is also expected to evolve to the AGB.

By examining the AGB population in our simulations we found that about 50% of the stars belong to the FG, and about 50% to the SG sub-population with Y = 0.254. From multiple synthetic realisations of the HB we determined a probability of only ~60% to find a single AGB star belonging to the SG with Y = 0.273 within a sample of 20 AGB objects, i.e. the size of Campbell et al. (2013) AGB sample.

We close this section by noting that the recent study by Charbonnel et al. (2013) also considered the enhanced initial He in SG stars as the likely explanation for the results by Campbell et al. (2013). The authors, however, did not perform any detailed synthetic HB simulations and employed a somewhat larger He-abundance range in SG stars compared to the current observational constraints (see also their discussion on this subject).

3.1. Consistency checks

To establish whether our theoretical predictions of R2 are biased by possible systematic errors in the evolutionary lifetimes predicted by the HB tracks hence the agreement with NGC 6752 is just fortuitous, we have performed additional tests. The most important one was to check whether our simulations are able to reproduce the value of R2 also for a cluster with a redder morphology, for which one expects all HB stars to end up evolving along the AGB. To this purpose we examined the detailed simulations of the HB of M 3, a cluster with metallicity very similar to NGC 6752 but with essentially all HB stars cooler than the Grundahl’s jump, presented in Dalessandro et al. (2013), obtained with the same evolutionary tracks and code employed here. The work by Dalessandro et al. (2013) makes use of far-ultraviolet(UV)/optical colour–magnitude diagrams for clusters with extended HBs, that enable to estimate the initial He abundance along the HB blue-tails, because in UV filters the HB becomes almost horizontal at high effective temperatures, and He variations cause variations of the UV magnitude levels (see Dalessandro et al. 2011, 2013, for details).

The analysis by Cho et al. (2005) determined an observed value R2 = 0.12 ± 0.03 for M 3 (we estimated the error bar on R2 using Poisson statistics). The simulations by Dalessandro et al. (2013) match the observed stellar distribution and magnitude levels along the cluster HB with a range of initial Y abundances equal to ~0.02, starting from a minimum Y = 0.246. From these simulations we derived R2 = 0.11, in agreement with the observed value.

We then considered another cluster with an extended blue HB and basically the same metallicity and age as NGC 6752, namely M 13 (see e.g. Dalessandro et al. 2013, and references therein). The HB ratio, HBR = (BR)/(B + V + R) –where B, V, and R are the number of stars at the blue side, within, and at the red side of the RR Lyrae instability strip, respectively– is 0.97 for M 13 and 1.0 for NGC 6752 (see e.g. Carretta et al. 2007).

The study by Cho et al. (2005) derived R2 = 0.07 ± 0.03, consistent with the result for NGC 6752. Our detailed simulations for the HB of M 13 have been discussed by Dalessandro et al. (2013), who estimated a range of initial Y from 0.246 up to ~0.31. This higher range of Y values compared to NGC 6752 is consistent with the fact that M 13 shows a wider range of Na abundances. Our HB simulations provide R2 = 0.06 (as for NGC 6752) again in very good agreement with the observations.

The study by Cho et al. (2005) provides another very interesting test for the evolutionary timescales of the models. They determined R2 = 0.16 ± 0.07 (we calculated the error bar from Poisson statistics) when considering all AGB stars but only the reddest component along the HB, e.g. all HB stars located at Teff roughly lower than the Grundahl’s jump location. On the simulation side, this reddest HB component corresponds to the synthetic HB stars with Y below 0.26 in the analysis by Dalessandro et al. (2013). Restricting the synthetic HB population to this component, while keeping the full AGB population, the simulations provide R2 = 0.15, again in agreement with observations. It is worth noticing that, within this test, if synthetic stars with Y higher than 0.26 (all HB stars hotter than the Grundahl’s jump) do not contribute to the AGB population, the simulations would give R2 = 0.11, much lower than the value observed. Unfortunately the large statistical error on the observed R2 of this test only provides an additional qualitative indication that HB models at this metallicity predict R2 values consistent with observations.

4. Discussion and conclusions

In this paper we have explored the scenario proposed by Campbell et al. (2013) to explain the [Na/Fe] abundances measured along the AGB of NGC 6752, and the low value of the R2 parameter, that envisages strong mass loss for all SG HB stars (located at Teff hotter than the temperature for the onset of radiative levitation), that prevents their evolution to the AGB. We have found that the required mass loss rates are of the order of 10-9 M yr-1, much higher than current empirical and theoretical constraints.

Our synthetic modelling – with our adopted stellar models – of the observed HB has provided a value of R2 consistent with the observations. Assuming that SG stars do not reach the AGB phase would produce a value of R2 that is too low compared to the observed one. Consistent synthetic HB analyses of M 3 (a cluster with the same metallicity and age, but a much redder HB morphology) and M 13 (a cluster with almost the same HB morphology of NGC 6752) provide R2 values again in agreement with observations. As a conclusion, our HB models correctly predict R2, at least at the metallicity of NGC 6752, and no strong mass loss during the core He-burning stage is required to reproduce the observed cluster R2.

The HB simulations also predict that for the size of NGC 6752 spectroscopic AGB sample, about 50% of the population is made of FG stars, and the rest belongs to the SG population with the lower Y mass fraction. The theoretical prediction that SG stars populate the AGB is consistent with the spectroscopic analysis by Villanova et al. (2009), who found in their small sample of HB stars cooler than the Grundahl’s jump, one object with high Na, brighter than the bulk of the HB stars. This object can be understood to be an evolved HB SG star that is moving towards the AGB. The presence of some evolved SG stars at Teff below the Grundahl’s jump is expected from the HB simulations (see Fig. 3).

On the other hand, only FG stars are observed in the AGB sample analysed spectroscopically (see Fig. 1). This is hard to explain, given the indications presented in this paper against the hypothesised enhanced mass loss. The situation becomes even more puzzling when considering the results for AGB stars in M 13. The spectroscopic analysis of a sample of its AGB and RGB stars by Johnson & Pilachowski (2012) shows that the AGB stars span almost the whole range of [Na/Fe] values covered by RGB stars. There is no indication that only FG stars inhabit the AGB. Given the close similarity of M 13 and NGC 6752 also in terms of the HB morphology, it is difficult to find a mechanism to deplete SG AGB stars that works in NGC 6752 but not in M 13.

We conclude by mentioning a final point that may be relevant to this issue. Campbell et al. (2013) measured Na abundances in RGB and AGB stars, and derived [Na/Fe] ratios by considering a uniform [Fe/H] value for all stars, as determined in the literature on RGB stars. In the case of M 13, Johnson & Pilachowski (2012) determined both Na and Fe abundances for their AGB stars. What would be important to consider is the effect of the interplay between radiative levitation and mass loss. For example, Michaud et al. (2011) show that both Na and Fe are strongly enhanced in the outer layers of a 0.59 M (initial Z = 0.0001) HB model at 14 000 K. If mass loss is efficient even only at the low level required to preserve the efficiency of radiative levitation, some amount of Na and Fe is certainly lost from the star. When convection homogenises the remaining envelope on the AGB, the surface abundances of these two elements will be depleted to some degree (maybe negligible) compared to the RGB abundances. If both Na and Fe are depleted by roughly the same amount, the measured [Na/Fe] on the AGB should reflect the original ratios along the RGB. If, however, only Na is measured and Fe is assumed to be equal to the RGB values, one might in principle underestimate the AGB [Na/Fe] ratio for objects affected by radiative levitation along the HB (as is the case of SG stars in NGC 6752 and M 13).

At the moment there are no available calculations of the full evolution from the HB to the AGB, self-consistently accounting for diffusion and radiative levitation, and mass loss. It would be certainly important to attempt these computations, with the aim of establishing whether the interplay between radiative levitation and mass loss plays a major role in explaining the puzzling chemical composition of AGB stars in NGC 6752 and the difference with M 13.


1

The whole BaSTI library is available at: http://www.oa-teramo.inaf.it/BASTI

2

Although SG stars appear to be O-depleted compared to the FG counterpart, as long as the C + N + O sum is kept constant, the evolutionary and structural properties of HB stars are barely affected (Pietrinferni et al. 2009). Given that SG stars in NGC 6752 do not show any evidence of CNO enhancement (Carretta et al. 2005; Yong et al. 2013), our assumption of adopting the same αenhanced mixture for stellar models representing both FG and SG stars is justified.

3

log  VC02 = −11.70 + 1.07 log (Teff/ 20000) + 2.13 (log (L/L) − 1.5) − 1.09 log (2M/M) + 0.97 log (Z/Z), where Teff is in Kelvin, and the other symbols have their usual meaning.

4

Spectroscopic analyses (see e.g. Behr 2003; Pace et al. 2006) have shown that iron and other heavy elements such as Ti and Cr can reach values around solar or even higher, as a consequence of radiative levitation in the atmosphere.

5

Only the morphology of the evolutionary tracks is affected by the choice of the bolometric corrections, as shown in Fig. 7 of Dalessandro et al. (2011), not the evolutionary timescales.

Acknowledgments

We warmly thank R. Gratton for the very useful suggestions and discussions about uncertainties of light element spectroscopy in globular cluster cool giants. S.C. acknowledges financial support from PRIN-INAF 2011, Multiple Populations in Globular Clusters: their role in the Galaxy assembly, (PI: E. Carretta), and from PRIN MIUR 2010-2011, project The Chemical and Dynamical Evolution of the Milky Way and Local Group Galaxies, prot. 2010LY5N2T (PI: F. Matteucci).

References

All Figures

thumbnail Fig. 1

Observed distributions of [Na/Fe] abundance ratios in NGC 6752 RGB stars as determined by Campbell et al. (2013) and Carretta et al. (2007) – dashed and solid line, respectively – and in AGB stars (Campbell et al. 2013, dash-dotted line – see text for details).

Open with DEXTER
In the text
thumbnail Fig. 2

Right panel: Hertzsprung-Russell (H-R) diagrams of stellar models with initial mass equal to 0.59 M, evolved from the ZAHB using the Reimers’ mass loss law with various values of the free parameter η (see labels). Left panel: as the right panel but for a model with initial mass equal to 0.55 M.

Open with DEXTER
In the text
thumbnail Fig. 3

Upper panel: CMD of NGC 6752 HB stars from Momany et al. (2002). Lower panel: one synthetic realisation of the cluster HB, with the same number of stars as observed. Open circles denote FG stars, filled triangles SG stars with Y = 0.254, open squares SG stars with Y = 0.273. The vertical lines define the boundaries of the three boxes in which the observed CMD has been partitioned. The location of the reddest vertical line coincides with the location of the Grundahl’s jump (see text for details).

Open with DEXTER
In the text
thumbnail Fig. 4

Mass distribution along the synthetic HB of Fig. 3. The meaning of the symbols is as in Fig. 3.

Open with DEXTER
In the text
thumbnail Fig. 5

Ratio of AGB to HB lifetimes as a function of the total mass for the models used in our simulations (see text for details).

Open with DEXTER
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.