Free Access
Volume 534, October 2011
Article Number A123
Number of page(s) 10
Section Galactic structure, stellar clusters and populations
Published online 19 October 2011

© ESO, 2011

1. Introduction

In the recent years, it has become clear that the formation of the most massive stellar clusters, the globular clusters (GCs), is a complex phenomenon, involving several episodes of star formation, with the latest generations of stars forming from material polluted by the ejecta of a fraction of the earliest (primordial) population, as first suggested by Norris et al. (1981; see also Gratton et al. 2001, 2004; Carretta et al. 2010). Evidence for the presence of these different populations stems from their chemical composition, in particular the abundances of the light elements involved in p-capture processes (C, N, O, Na, Mg, Al, Si: Osborn 1971; see also Cottrell & Da Costa 1981; Kraft 1994; Denisenkov & Denisenkova 1989; Carretta et al. 2009a,b, and references therein), as well as from the splitting of sequences in the colour magnitude diagram (see e.g., Piotto 2008). While in the recent years some spectacular splitting of the main sequence (MS) have been found (see e.g., Piotto et al. 2007), the most dramatic evidence concerns the horizontal branch (HB) and was found already several decades ago (Harris 1974). This splitting of the HB has long been related to the second parameter problem (Sandage & Wildey 1967; van den Bergh 1967) but it was only quite recently that a plausible cause has been identified in variations of the He content related to the multiple populations (D’Antona & Caloi 2004; but see Norris et al. 1981, for a very similar early suggestion). As pointed out by D’Antona et al. (2002), variations in the He content among different stellar generations in GCs are expected to be present, correlated to the variation of the abundances of the p-capture process elements. These He abundance variations are invoked to explain the multiple MSs seen in mono-metallic GCs like NGC 2808. However, because He-rich stars are expected to burn H at a faster rate than the He-poorer ones while on the MS, their progeny on the HB is expected to be less massive, then bluer, than the progeny of He-normal stars. Carretta et al. (2007) and Gratton et al. (2010) discussed the general correlation between the abundances of p − capture elements and the colour of HB stars, and found indeed that a close relation exists between the extension of the Na-O and Mg-Al anticorrelations, and that of the HB, precisely as expected in this interpretative scenario. These studies also showed that the main parameter driving the whole phenomenon of multiple populations is the cluster mass (confirming previous similar results obtained by Recio-Blanco et al. 2006; Carretta et al. 2006, 2009a), and that the second parameter can be explained almost entirely by a combination of variations of the age and mass of GCs.

While this is extremely intriguing, a definitive confirmation of this scenario with direct determination of the chemical composition of HB stars is required. The surface composition of the hottest stars on the HB, those with temperatures  >11 500 K (the location of the so-called u-jump: Grundahl et al. 1999), are known to be heavily influenced by sedimentation caused by diffusion and by levitation through radiative pressure effects (see, e.g., Behr et al. 1999, 2003; see also Grundahl et al. 1999; Caloi 1999). However, the same studies revealed that abundances for cooler stars can be reliably determined, although the low S/N and wavelength range of the spectra available at the epoch prevented an accurate determination of the abundances of those elements of interest here (O, Na). Behr et al. (2003) first obtained He abundances for cool BHB stars but the large error bars did not allow the authors to use them as diagnostics for the multiple population scenario. Very recently, a few investigations tried to measure the abundances of O and Na, and even of He, in such stars with sufficiently good accuracy. The pioneering study was the analysis of five blue HB stars in NGC 6752 by Villanova et al. (2009). This cluster has a very blue HB, and the observed stars (with temperatures  >11 500 K) are among the coolest HB stars. They behave as expected, most of them being Na-poor and O-rich, save for one star, which is more luminous than the others, and likely is an evolved object that probably started its HB evolution at a much higher temperature. Two other studies (Marino et al. 2011; and Villanova et al. 2011) focused on M 4, the nearest GC. This cluster has a much redder HB than NGC 6752, so that stars on both the red and blue sides of the RR Lyrae instability strip could be studied. Again, as expected, the red HB stars were found to be Na-poor and O-rich, and the blue ones Na-rich and O-poor.

In this paper we present a similar investigation of the HB of NGC 2808. NGC 2808 is among the most massive Galactic GCs, and has become the most classical example of multiple populations. For a long time it was known to have a very peculiar HB, with a multimodal distribution of colours and masses (Harris 1974; see Bedin et al. 2000; Momany et al. 2004; and Dalessandro et al. 2011, for recent accurate photometry). D’Antona et al. (2004, 2005) observed a widening of the MS, later found by Piotto et al. (2007) to consist of three distinct sequences, which can be explained with different values of the He content. D’Antona et al. (2004) correlated the different MSs with the different populations of HB stars. Low S/N, high-resolution spectra along the BHB of NGC 2808 were obtained by Recio-Blanco et al. (2004) and were used to derive radial and rotational velocities. Pace et al. (2006) obtained intermediate-resolution spectra for a number of HB stars, but they were only able to derive the Fe abundances from these spectra. They found that stars cooler than the u-jump at  ~12 000 K provide values similar to those obtained for red giants, while metallicity increases with temperature for stars warmer than this limit. Carretta et al. (2004, 2006) studied the Na-O anticorrelation in several hundred red giants, and found an extended anticorrelation, with different peaks in the distribution. Bragaglia et al. (2010a) used this same observational material to derive evidence of the variation of the He content along the RGB, with results that on the whole agree with those obtained from the colours of MS and HB. However, while very suggestive, the connection between the chemical anomalies along the RGB and the multiple MS and HB could not be definitively shown by investigations that used very different techniques on stars in different evolutionary phases. To provide this definitive evidence, about one year ago, Bragaglia et al. (2010b) obtained spectra for one star on the blue MS, and another on the red one. They found that as expected, the blue MS star is Na- and Al-rich and Mg-poor, while the red one is Na- and Al-poor and Mg-rich. In this paper, we intend to complete the observational pattern by determining the abundances of O and Na for several tens of HB stars, using spectra acquired with FLAMES+GIRAFFE at VLT. We noticed immediately that while we possibly detected the He I line at 5876 Å in the spectra of several of the BHB stars, stars are so faint that only prohibitively long exposure times would have ensured the S/N required for accurate measure of the He abundance, as needed to distinguish primordial abundances from those expected for second generation stars in the multiple population scenarios1. Therefore we limited our analysis to Na and O abundances, for which the S/N required to provide sensible results is much lower.

The structure of the paper is the following: in Sect. 2 we present the observational data; in Sect. 3 we explain our analysis methods; in Sect. 4 we discuss the Na-O anticorrelation; conclusions are drawn in Sect. 5.

2. Observation

We used the high multiplex capability of GIRAFFE at VLT (Pasquini et al. 2004) to acquire spectra for 49 stars along the red HB (RHB) and 36 on the blue HB (BHB) of NGC 2808. This was the maximum number of fibres that we were able to allocate to stars in these evolutionary phases given the practical limitations caused by the finite size of the fibre heads2. In addition, 34 fibres were used to point faint stars along the RGB (below the RGB bump), the remaining ones being used to acquire sky spectra. The median spectra from these last fibres were subtracted from those used for the stars. This was of particular relevance here, because the observed stars are typically very faint. Two spectral configurations were used, HR12 and HR19, providing high-resolution spectra including the strongest features of O I (the IR triplet at 7771-74 Å) and Na I (the resonance D doublet at 5890-96 Å, as well as the subordinate strong doublet at 8183-94 Å) accessible from ground. The O I triplet and the D lines of Na I are the only lines of these two elements that can be used to determine O and Na abundances without a prohibitively observing time for NGC 2808. A few lines of Mg, Al, Si, Ca, Fe, and Ba were also included in the selected observing ranges.

Our programme was executed in service mode, and not all planned observations were actually carried out, so that in total we had 1.5 h (2 × 45 min exposures) and 6.15 h (7 × 53 min) of observation with the gratings HR12 and HR19, respectively. The S/N of the summed spectra is typically  ~20 for the RHB stars. For the BHB stars, spectra with grating HR12 have S/N ~ 15 and those with grating HR19 S/N ~ 10. The spectra were reduced using the ESO FLAMES GIRAFFE pipeline version 2.8.7. Sky subtraction, combination of individual exposures for each star, translation to rest-frame and continuum tracing were performed within IRAF3. Telluric lines were removed from the longest wavelength spectra by subtracting the average spectrum of those BHB stars for which the observations have S/N > 8. This combined spectrum has a S/N ~ 50 (much higher than the S/N of the individual programme stars), and was obviously taken with the same airmass than the programme star, so that the excision of the telluric lines turned out to be excellent. Examples of spectra are shown in Fig. 1.

thumbnail Fig. 1

Portions of the spectra of the RHB star #8342 (left column) and of the BHB star #7189 (right column). From top to bottom: the regions including the O I triplet at 7771-74 Å, the Na D lines at 5890-96 Å, and the Na doublet at 8183-94 Å. The spectra in the region of the O I triplet have been smoothed by a Gaussian with FWHM = 0.3 Å. The Na D interstellar lines are clearly visible in the median panels.

thumbnail Fig. 2

Colour magnitude diagram of the inner region of NGC 2808. Filled circles are confirmed cluster members, while open circles (13 and 14 red HB and red giants respectively) are non cluster members. Different colours are used for stars on the RGB (green), RHB (red), and BHB (blue). Stars not observed are shown as dots.

Table 1

Basic data for program stars.

Figure 2 shows the location of the programme stars on the colour magnitude diagram of NGC 2808. Our NGC 2808 ground-based photometric catalogue (see Momany et al. 2004) consists of UBV observations obtained at the Wide-Field Imager (WFI) mounted on the 2.2 m ESO-MPI (La Silla, Chile). The WFI catalogue has a total field of view of 34′ ×    33′. Photometric data for the programme stars are listed in Table 1. The K magnitudes are from the 2MASS point source catalogue (Skrutskie et al. 2006). Note that no information on membership of the programme stars to the cluster was available prior to the observations. The stars selected for observations lie close to the mean loci of the colour magnitude diagram of NGC 2808. Virtually no field contamination is expected for the BHB, while some field interlopers may be present in our RHB sample. We then determined their membership from the radial velocities (also listed in Table 1), because the fairly high radial velocity of NGC 2808 ( + 101.6 ± 0.7 kms-1, Harris 1996) ensured that very few field stars with similar radial velocities were to be expected. Thirty-seven RHB stars turned out to be cluster members based on this criterion. A close inspection revealed that most of the cluster members are grouped in a very narrow region of the colour magnitude diagram, while the field stars are typically slightly bluer and fainter than the real RHB stars. After the abundance analysis, we excluded another star (#36742) because its Fe abundance is much larger than typical for stars in NGC 2808. Note that the radial velocity and the K magnitude of this star only marginally agree with those of the NGC 2808 RHB members. A similar analysis for the RGB stars indicated that 22 stars are cluster members. The average radial velocity of the 36 RHB stars is  +101.0 ± 0.9 kms-1, in agreement with the value listed by Harris (1996). The rms scatter is 5.7 kms-1. The average radial velocity of the lower RGB (21 stars) is  + 99.9 ± 3.0 kms-1, with an rms scatter of 9.2 kms-1. The difference between the rms scatter obtained for RHB and RGB stars is not significant.

Not unexpectedly, most of the BHB stars of NGC 2808 turned out to be hotter than 11 500 K, with only six stars cooler than this limit. In the following, we will present our result for these last stars and for the RHB member stars. Results for the other stars will be presented elsewhere.

We also measured radial velocities for the BHB stars cooler than the Grundahl u-jump, though with fairly large errors because they are based on the O I triplet lines alone. Three of these stars are in common with the study of Recio-Blanco et al. (2004). These are 48 003 = 2333, 10 449 = 3949, and 37 289 = 3841 (the first number is our identification, the second that from Recio-Blanco et al. 2004); however, we have radial velocity only for two of them. The comparison is fairly good, with a difference of 2 ± 3 kms-1.

Various authors found that horizontal branch stars rotate with velocities up to a few tens of kms-1 for BHB stars (Peterson et al. 1995; Behr et al. 2000a,b) and slightly less for RHB ones (Carney et al. 2008). No really fast rotator was found among the 16 BHB stars observed by Recio-Blanco et al. (2004) in NGC 2808, where the maximum rotation velocity was 13 ± 4 kms-1 for their star 3841, which is our star 37 289. We have checked our sample for fast rotators. To this aim, we considered the FWHM of the fit to the cross-correlation peak obtained using as a template a spectrum of star of similar atmospheric parameter taken with the same instrument. To derive vsini from the FWHM the contribution caused by the instrumental profile, macro- and micro-turbulence, and thermal Doppler broadening should be taken into account, both for the stars and the template, and the relation should be calibrated (see Lucatello & Gratton 2003). However, the FWHM of the cross correlation of spectra of stars with similar atmospheric parameters with the same template gives a qualitative indication of the rotational rate. To this aim we used the spectra taken with grating HR12, because in HR19 the number of stellar lines is much smaller and the spectra are dominated by telluric lines.

Among red HB stars we found no evidence of a considerable variation in the measured FWHMs, which range from  ~ 28 to  ~ 40 kms-1. While these values cannot exclude the presence of stars rotating as fast as vsini ~ 5–7 kms-1 (typical values for field RHB stars observed by Carney et al. 2008), just the accounting of the instrumental profile of template and object stars at the spectral resolution used (R ≃ 18   700) leaves no room for high rotational velocities. We cannot draw any conclusions concerning blue HB stars: the paucity of lines combined with the quite low signal-to-noise ratio makes the estimates of FWHM fairly uncertain as well as highly affected by noise, residuals of telluric subtraction, etc.

3. Analysis

3.1. Atmospheric parameters

The programme stars cover a wide range in effective temperature. It is therefore difficult to mantain a uniform criterium to derive atmospheric parameters. In addition, the low S/N and the limited spectral coverage of the spectra prevent a full spectroscopic analysis. The adopted procedure was the following.

For the RHB and RGB stars, effective temperatures were derived from the B − V and V − K colours, using the calibration of Alonso et al. (1999, with the erratum of Alonso et al. 2001). The colours were dereddened using the E(B − V) values from the updated on-line version of the Harris (1996) catalogue and the E(V − K)/E(B − V) value from Cardelli et al. (1989). The calibrations require input values for the metallicity [A/H]. We adopted the value obtained by Carretta et al. (2006). We assign weight 4 to the B − V colours, and 1 to the V − K ones, because the programme stars are very faint for the 2MASS observations. For the blue HB stars we used the (B − V) − Teff calibration by Kurucz4.

Errors are mainly caused by uncertainties in the colours. For the RHB stars, an error of 0.01 mag in B − V causes an error of 25 K in the effective temperature; a slightly smaller error of 20 K is obtained for red giants. A similar photometric error for the BHB stars causes an error of  ~ 300 K. We assumed errors of 50 and 300 K as representative values for the internal errors in the temperatures. Systematic errors caused by scale errors or incorrect parameters for the cluster are likely larger. However, they are not really important in our discussion.

The surface gravities were obtained from the masses, luminosities, and effective temperatures. For the masses we adopted values of 0.9, 0.7, and 0.5 M for stars on the RGB, RHB, and BHB (see Gratton et al. 2010, for a discussion of adequate values for the different sequences). Note that these values are not very critical, because a variation of 0.1 M changes gravities by 0.07 and 0.10 dex for RHB and BHB stars. The effects on [Na/Fe] are of 0.11 and 0.08 dex, and those on [O/Fe] 0.02 and 0.13 dex, for BHB and RHB respectively. The bolometric corrections needed in these derivations were obtained using calibrations consistent with those used for the effective temperatures (Alonso et al. 1999, for the red giant and RHB stars, and Kurucz for the BHB stars). The distance moduli used were taken from the Harris catalogue.

Errors in gravities are very small. The assumption about masses is likely correct within 10% (0.04 dex error in the gravities), while that on the effective temperature and luminosity causes errors in gravities not larger than  ~2% for the RHB stars and red giants, and  ~12% for the BHB stars. The error in gravities is then not larger than 0.05 dex for the cool stars, while it may be 0.10 dex for the warm ones.

The same metal abundance of [A/H] = –1.14 and microturbulence velocity of 1.8 kms-1 were adopted for all RHB stars. The metal abundance is the average value determined in Carretta et al. (2006), and agrees well with the average [Fe/H] =  − 1.18 ± 0.01 (rms of 0.07 dex) we determined from the 36 RHB stars. The adopted microturbulence velocity sets to zero the average of the trends of Fe abundances with expected line strength (see Magain 1984, for a justification of this approach). We prefer to adopt this average value rather than values appropriate for each star because we were able to measure equivalent widths (EW) only for a limited number of Fe lines (in the range 6–16) for each star, with fairly large errors (~13 mÅ, evaluated from the rms of measures for the same line across the RHB sample), and over a limited range of EWs. The uncertainties in the individual slopes, which agree with the observed rms scatter, would translate into a large error bar of  ± 0.5 kms-1. This is likely more than the expected star-to-star scatter, which we estimated to be  ~0.3 kms-1 in analogy with the RGB stars studied by Carretta et al. (2009a). We will adopt this last value as the error bar of our microturbulent velocities. We were unable to derive microturbulent velocities for BHB stars because too few lines were measured. The value we adopted (2 kms-1) is in the middle of the range usually found in previous analysis of BHB stars (Lambert et al. 1992; Behr et al. 1999, 2000b; Kinman et al. 2000; Fabbian et al. 2005; Villanova et al. 2009; Marino et al. 2011).

While a single Fe II line was measured in our spectra (at 5991.38 Å), abundances derived from this line agree well with those obtained from the Fe I lines: on average we obtained [Fe II/H] =  − 1.13 ± 0.03. This supports the choice of the atmospheric parameters adopted throughout this paper.

Table 2

Atmospheric parameters and abundances.

Table 2 lists the effective temperatures Teff and surface gravities log g we used in the analysis of the programme stars, as well as the abundances we obtained from our analysis for lines of Fe I, Fe II, O I, Na I, Mg I, Si I, Ca I, and Ba II. Abundances were estimated from equivalent widths. The oscillator strengths for Fe lines were generally obtained from the VALD database (Kupka et al. 2000). The same source was used for the Mg line at 8213.04 Å, for four lines of Si (5948.55, 7918.38, 7932.35 and 7944.00 Å), for the Ca line at 5857.46 Å, and for the Ba II line at 5853.69 Å. Note that not all stars were observed with both gratings. As a consequence, we were not able to measure the abundances of all elements in each star. For this reason we have no O abundance for star #10449, for instance.

For O I and Na I lines, the abundances include corrections due to departures from LTE, which are large for the transition and stars of interest here. The non-LTE corrections used are estimated from standard statistical equilibrium calculations. They include integration of the radiation transfer equation throughout the stellar atmosphere, with corrections to the populations of the individual levels and to the source function computed at all optical depth. There are uncertainties, related to the approximation present and to the poor knowledge of the relevant cross sections. In particular the treatment of collisions with H I atoms is critical and poorly known. Because different recipes are used, various non-LTE estimates produce somewhat different results, although there is generally quite a good consensus on the sign and order of magnitude of the corrections, in particular for quite simple atoms like O I and Na I. For the RHB stars, we used the corrections listed by Gratton et al. (1999), whose cross sections for collisions with H I atoms were tuned to reproduce observations of RR Lyrae stars. The corrections are  ~0.5 dex for O I and  ~0.2 dex for Na, in both cases decreasing the LTE abundances. Very recently, new Na non-LTE corrections were published by Lind et al. (2011), making use of a more extended model atom and moreover updated data about cross sections. For the lines used and the combination of atomic parameters appropriate for the programme stars, these corrections are about 0.3 dex larger than those adopted here and this difference is quite uniform. Were we to adopt these corrections, the Na abundances would be about 0.3 dex smaller. However, we prefer to keep here the older correction by Gratton et al. (1999) for uniformity with the red giants studied by Carretta et al. (2006). However, even these corrections were not computed for stars as hot as our BHB stars. In this case, we used the results of the statistical equilibrium calculations by Takeda (1997) for O I (with corrections to O abundances in the range 1 − 1.5 dex, and by Mashonkina et al. (2000) for Na I (in this case the corrections are  ~0.5 dex). We alert the reader that there might be some offsets between these different sets of corrections. These offsets are difficult to quantify accurately. If we compare non-LTE abundance corrections for O triplet lines by Takeda (1997) with those by Gratton et al. (1999) for stars with Teff = 7000 K, however, where both are available, we find that the first ones are smaller by  ~0.1 dex (~0.5 vs.  ~0.6 dex). A similar comparison for the Na non-LTE corrections of Mashonkina et al. (2000) shows that these are smaller by  ~0.2 dex than those by Gratton et al. (1999) (~0.4 vs.  ~0.6 dex). We conclude that uncertainties in the non-LTE corrections may be likely as large as a few tenths of a dex (see also comparisons in Lind et al. 2011, which are appropriate to cooler stars, however). Hence, while these uncertainties are not negligible, it is difficult to conceive that they may cause the main results of this paper.

The doublet of Al I at 7835-36 Å is actually in the observed spectral range. We carefully examined the relevant region in the spectra of all RHB stars, but we did not find any evidence for it. The line is also not detectable in the spectrum we obtained by summing all 36 HB stars. From this last spectrum (which has a S/N > 100 at the relevant wavelength) we derived an upper limit of 10 mÅ to the equivalent width for each component of the doublet. This yields an upper limit of  [Al/Fe]  < 0.15, which is similar to the Al abundances obtained for the most Al-poor RGB stars in NGC 2808 (Carretta et al. 2009b). No Al detection is possible in the low S/N spectra of BHB stars.

3.2. Error analysis

Table 3

Sensitivity and error analysis.

Error analysis was made in the usual way, by repeating the abundance derivation by modifying a single parameter each time. Relevant data are given in Table 3. The last column gives an estimate of the total internal errors, estimated using the sensitivities listed above, as well as the errors in the individual parameters given in lines 2 and 3 for blue and red HB stars, respectively. In general, the main source of uncertainty is the measure of the equivalent widths, which is not surprising considering the fairly low S/N of our spectra. In a few cases (Fe, O, and Na) a significant contribution comes from errors in the effective temperatures and in the microturbulent velocities.

For most elements (Fe, Mg, Si, Ca, and Ba) the dispersion of individual values (for RHB stars alone) agrees fairly well with these estimates of the internal errors. The case of Na and O will be discussed in the next section.

4. The Na-O anticorrelation along the HB of NGC 2808

Figure 3 shows the Na-O anticorrelation we obtained for the HB stars of NGC 2808. Different symbols are used for blue and red HB stars. For comparison we also plotted the Na-O anticorrelation for red giants by Carretta et al. (2006). We recall that the extreme blue HB stars, hotter than 11 500 K, are not considered here, because their surface abundance is not related in a simple way to their original composition. The moderately BHB stars for which we may provide sensible Na and O abundances are all O-poor and Na-rich, with average values of [O/Fe] =  −0.59 ± 0.12 and [Na/Fe] =  + 0.44 ± 0.06. The rms scatter of individual values of 0.24 dex and 0.15 dex agree with the internal errors of Table 3, so that we cannot conclude about any internal scatter in this small sample of six stars, even if there is a hint that the stars might be distributed into two groups. We also notice that systematic errors caused by, e.g. the corrections for non-LTE effects (which are rather large for the O lines) are certainly much larger than the error bars given above for the average abundances, which stems from the dispersion of individual values alone. However, the comparison of the Na and O abundances for these stars with those obtained for red giants suggests that they descend from the moderately O-poor red giants.

On the other hand, all RHB stars are much richer in O than the BHB stars, with positive values of [O/Fe]. These stars span a fairly large range of Na abundances, and the rms scatter of [Na/Fe] values of 0.18 ± 0.03 dex appears to be significantly larger than the observational error of 0.11 dex. While not extremely significant, there is a clear hint for a Na-O anticorrelation. The correlation coefficient is r =  −0.34 over 31 stars, which has less than 5% chance of being a random effect. In addition, we notice that there is a very significant correlation between the Na abundances obtained from the resonance D-doublet and the subordinate one at 8183-94 Å. This suggests that the red HB stars have not a single value of Na and O abundances, with only observational scatter; instead they seem to include at least two different populations. This means that only a fraction of the RHB stars of NGC 2808 descends from the primordial population (see Carretta et al. 2010, for a definition of primordial population). The comparison with red giants shows that there is a very good correspondence between the RHB stars and the O-rich sequence in NGC 2808; we incidentally note here that also O-rich red giants display a distinct spread in Na abundances, so that many of them are classified as belonging to the intermediate population according to the definition by Carretta et al. (2009a).

thumbnail Fig. 3

Left panel: Na-O anticorrelation for HB stars of NGC 2808. Filled circles are BHB stars; open circles are RHB stars. Right panel: the same, but for red giants (from Carretta et al. 2006).

While we cannot determine directly the He abundances for the various populations on the HB of NGC 2808, we may propose values that are compatible with their colours and magnitudes. This matter has been considered in depth in previous investigations. For instance, D’Antona et al. (2005) considered three populations with Y = 0.24, Y = 0.26–0.29, and Y = 0.40, with respectively 50, 30, and 20% of the stars, and found that this mix is able to reproduce both the main sequence and the HB of NGC 2808. The first group corresponds to the RHB, the second to the BHB, and the third to extreme BHB (not sampled in this paper). The same proportions are found for P, I, and E stars along the RGB by Carretta et al. (2010). This suggests the identification of the He-poor stars and RHB with the P-stars, of the He-intermediate and BHB stars with the I-stars, and of the He-rich and extreme BHB ones with the E-stars. As expected from this scenario, we indeed found a clear correlation between the Na/O abundance ratio and colour and temperature along the HB in NGC 2808 (see Fig. 4). This is obvious when comparing the O-poor, Na-rich BHB with the O-rich RHB stars. While we were unable to analyse the warmer stars along the HB of NGC 2808, we think this is a strong circumstantial confirmation of the scenario where the location of stars along the HB is determined also by their helium content. However, because there is a hint for a trend of Na and O abundances with colour also among the red HB stars alone, we suggest that there are more than three stellar populations in NGC 2808, and that the He-poor group itself divides into further groups. In order to produce the distribution in colours of the stars, the helium abundance differences among these groups should be small ( ≤ 0.01). However, the presence of significant variations in Na content has implications in a scenario for the formation of this cluster.

thumbnail Fig. 4

Left: run of the [Na/O] abundance ratio with B − V colour along the HB of NGC 2808. Filled circles and open circles are for blue and red HB stars, respectively. Right: the same, but with Teff rather than colour.

5. Conclusions

We have obtained spectra for several tens stars on the HB of NGC 2808, which is the prototype mono-metallic, multiple stellar population cluster. Several red HB stars were observed, and six on the blue HB cooler than the gap at 11 500 K studied by Grundahl, Caloi, and Behr. We found that as expected, the blue HB stars are O-poor and Na-rich, and the red HB ones are O-rich and Na-poor. This strongly supports the identification of the BHB as the progeny of the intermediate MS, and of the RHB with the red MS, made by D’Antona et al. (2005, no reliable abundances could be obtained for the extremely blue HB stars, which should be the progeny of the blue MS stars). As additional information, we found that some Na-O anticorrelation can be found also among the RHB stars alone, suggesting that this sequence, which includes almost half of the HB stars in NGC 2808 itself, is made of at least two distinct groups of stars. Hence there are more than three stellar populations in NGC 2808.


In practice, among stars with temperature  <11 500 K, the He I line is detected only for star #32 324, and even in that case with such an error on the equivalent width that the estimate on He abundance is not useful. For the coolest stars the He line is not detectable in our noisy spectra even when they were summed together.


All stars were chosen to be free from any companion closer than 2 arcsec and brighter than V + 2 mag, where V is the target magnitude.


IRAF is distributed by the National Optical Astronomical Observatory, which are operated by the Association of Universities for Research in Astronomy, under contract with the National Science Foundation.




This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research has made use of the NASA’s Astrophysical Data System. This research has been funded by PRIN INAF “Formation and Early Evolution of Massive Star Clusters”. We thank an anonymous referee for his/her useful suggestions.


  1. Alonso, A., Arribas, S., & Martinez-Roger, C. 1999, A&AS, 140, 261 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Alonso, A., Arribas, S., & Martinez-Roger, C. 2001, A&A, 376, 1039 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Bedin, L. R., Piotto, G., Zoccali, M., et al. 2000, A&A, 363, 159 [NASA ADS] [Google Scholar]
  4. Behr, B. B., Cohen, J. G., McCarthy, J. K., & Djorgowski, S. G. 1999, ApJ, 517, L135 [NASA ADS] [CrossRef] [Google Scholar]
  5. Behr, B. C., Djorgovski, S. G., Cohen, J. G., et al. 2000a, ApJ, 528, 849 [NASA ADS] [CrossRef] [Google Scholar]
  6. Behr, B. C., Cohen, J. G., & McCarthy, J. K. 2000b, ApJ, 531, L37 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  7. Bragaglia, A., Carretta, E., Gratton, R. G., et al. 2010a, A&A, 519, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bragaglia, A., Carretta, E., Gratton, R. G., et al. 2010b, ApJ, 720, L41 [NASA ADS] [CrossRef] [Google Scholar]
  9. Caloi, V. 1999, A&A, 343, 904 [NASA ADS] [Google Scholar]
  10. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
  11. Carney, B. W., Gray, D. F., Yong, D., et al. 2008, AJ, 135, 892 [NASA ADS] [CrossRef] [Google Scholar]
  12. Carretta, E. 2006, AJ, 131, 1766 [NASA ADS] [CrossRef] [Google Scholar]
  13. Carretta, E., Bragaglia, A., & Cacciari, C. 2004, ApJ, 610, L25 [NASA ADS] [CrossRef] [Google Scholar]
  14. Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2006, A&A, 450, 523 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Carretta, E., Recio-Blanco, A., Gratton, R. G., Piotto, G., & Bragaglia, A. 2007, ApJ, 671, L125 [NASA ADS] [CrossRef] [Google Scholar]
  16. Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2009a, A&A, 505, 117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2009b, A&A, 505, 139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2010, A&A, 516, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Cottrell, P. L., & Da Costa, G. S. 1981, ApJ, 245, L79 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dalessandro, E., Salaris, M., Ferraro, F. R., et al. 2011, MNRAS, 410, 694 [NASA ADS] [CrossRef] [Google Scholar]
  21. D’Antona, F., & Caloi, V. 2004, ApJ, 611, 871 [NASA ADS] [CrossRef] [Google Scholar]
  22. D’Antona, F., Caloi, V., Montalban, J., Ventura, P., & Gratton, R. 2002, A&A, 395, 69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. DAntona, F., Bellazzini, M., Caloi, V., et al. 2005, ApJ, 631, 868 [NASA ADS] [CrossRef] [Google Scholar]
  24. Denisenkov, P. A., & Denisenkova, S. N. 1989, A. Tsir., 1538, 11 [NASA ADS] [Google Scholar]
  25. Fabbian, D., Recio-Blanco, A., Gratton, R. G., & Piotto, G. 2005, A&A, 434, 235 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Gratton, R. G., Carretta, E., Eriksson, K., & Gustafsson, B. 1999, A&A, 350, 955 [NASA ADS] [Google Scholar]
  27. Gratton, R. G., Bonifacio, P., & Bragaglia, A. 2001, A&A, 369, 87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Gratton, R., Sneden, C., & Carretta, E. 2004, ARA&A, 42, 385 [NASA ADS] [CrossRef] [Google Scholar]
  29. Gratton, R. G., Carretta, E., Bragaglia, A., Lucatello, S., & D’Orazi, V. 2010, A&A, 517, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Grundahl, F., Catelan, M., Landsman, W. B., Stetson, P. B., & Andersen, M. I. 1999, ApJ, 524, 242 [NASA ADS] [CrossRef] [Google Scholar]
  31. Harris, W. E. 1974, ApJ, 192, L161 [Google Scholar]
  32. Harris, W. E. 1996, AJ, 112, 1487 [NASA ADS] [CrossRef] [Google Scholar]
  33. Kinman, T., Castelli, F., Cacciari, C., et al. 2000, A&A, 364, 102 [NASA ADS] [Google Scholar]
  34. Kraft, R. P. 1994, PASP, 106, 553 [NASA ADS] [CrossRef] [Google Scholar]
  35. Kupka, F., Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C., & Weiss, W. W. 2000, Balt. Astron., 9, 590 [NASA ADS] [Google Scholar]
  36. Lambert, D. L., McWilliam, A., & Smith, V. V. 1992, ApJ, 386, L685 [NASA ADS] [CrossRef] [Google Scholar]
  37. Lind, K., Asplund, M., Barklem, P. S., & Belyaev, A. K. 2011, A&A, 528, A103 [Google Scholar]
  38. Lucatello, S., & Gratton, R. G. 2003, A&A, 406, 691 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Magain, P. 1984, A&A, 134, 189 [NASA ADS] [Google Scholar]
  40. Marino, A. F., Villanova, S., Milone, A. P., et al. 2011, ApJ, 730, L16 [NASA ADS] [CrossRef] [Google Scholar]
  41. Mashonkina, I. I., Shimanskii, V. V., & Sakhibullin, N. A. 2000, Astron. Rep., 44, 790 [NASA ADS] [CrossRef] [Google Scholar]
  42. Momany, Y., Bedin, L. R., Cassisi, S., et al. 2004, A&A, 420, 605 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Norris, J., Cottrell, P. L., Freeman, K. C., & Da Costa, G. S. 1981, ApJ, 244, 205 [NASA ADS] [CrossRef] [Google Scholar]
  44. Osborn, W. 1971, Obs., 91, 223 [Google Scholar]
  45. Pace, G., Recio-Blanco, A., Piotto, G., & Momany, Y. 2006, A&A, 452, 493 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Pasquini, L., Castillo, R., Dekker, H., et al. 2004, SPIE, 5492, 136 [Google Scholar]
  47. Peterson, R. C., Rood, R. T., & Crocker, D. A. 1995, ApJ, 453, 214 [NASA ADS] [CrossRef] [Google Scholar]
  48. Piotto, G. 2008, MSAIt, 79, 334 [Google Scholar]
  49. Piotto, G., Bedin, L. R., Anderson, J., et al. 2007, ApJ, 661, L53 [NASA ADS] [CrossRef] [Google Scholar]
  50. Recio-Blanco, A., Piotto, G., Aparicio, A., & Renzini, A. 2004, A&A, 417, 597 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Recio-Blanco, A., Aparicio, A., Piotto, G., De Angeli, F., & Djorgovski, S. G. 2006, A&A, 452, 875 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Sandage, A., & Wildey, R. 1967, ApJ, 150, 469 [NASA ADS] [CrossRef] [Google Scholar]
  53. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  54. Takeda, Y. 1997, PASJ, 49, 471 [NASA ADS] [Google Scholar]
  55. van den Bergh, S. 1967, AJ, 72, 70 [NASA ADS] [CrossRef] [Google Scholar]
  56. Ventura, P., D’Antona, F., Mazzitelli, I., & Gratton, R. 2001, ApJ, 550, L65 [NASA ADS] [CrossRef] [Google Scholar]
  57. Villanova, S., Piotto, G., & Gratton, R. G. 2009, A&A, 499, 755 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Villanova, S., Geisler, D., Piotto, G., & Gratton, R. 2011, ApJ, submitted [Google Scholar]

All Tables

Table 1

Basic data for program stars.

Table 2

Atmospheric parameters and abundances.

Table 3

Sensitivity and error analysis.

All Figures

thumbnail Fig. 1

Portions of the spectra of the RHB star #8342 (left column) and of the BHB star #7189 (right column). From top to bottom: the regions including the O I triplet at 7771-74 Å, the Na D lines at 5890-96 Å, and the Na doublet at 8183-94 Å. The spectra in the region of the O I triplet have been smoothed by a Gaussian with FWHM = 0.3 Å. The Na D interstellar lines are clearly visible in the median panels.

In the text
thumbnail Fig. 2

Colour magnitude diagram of the inner region of NGC 2808. Filled circles are confirmed cluster members, while open circles (13 and 14 red HB and red giants respectively) are non cluster members. Different colours are used for stars on the RGB (green), RHB (red), and BHB (blue). Stars not observed are shown as dots.

In the text
thumbnail Fig. 3

Left panel: Na-O anticorrelation for HB stars of NGC 2808. Filled circles are BHB stars; open circles are RHB stars. Right panel: the same, but for red giants (from Carretta et al. 2006).

In the text
thumbnail Fig. 4

Left: run of the [Na/O] abundance ratio with B − V colour along the HB of NGC 2808. Filled circles and open circles are for blue and red HB stars, respectively. Right: the same, but with Teff rather than colour.

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.