A&A 440, 901-908 (2005)
DOI: 10.1051/0004-6361:20053019
R. G. Gratton1 - A. Bragaglia2 - E. Carretta2 - F. De Angeli3 - S. Lucatello1 - G. Piotto3 - A. Recio Blanco3,4
1 - INAF - Osservatorio Astronomico di Padova, Vicolo
dell'Osservatorio 5, 35122 Padova, Italy
2 -
INAF - Osservatorio Astronomico di Bologna, Via Ranzani 1,
40127 Bologna, Italy
3 -
Dipartimento di Astronomia, Universitá di Padova, Vicolo
dell'Osservatorio 5, 35122 Padova, Italy
4 -
Observatoire Astronomique de la Côte d'Azur,
Boulevard de l'Observatoire, BP 4229, 06304 Nice Cedex 4, France
Received 8 March 2005 / Accepted 4 May 2005
Abstract
Accurate reddenings for Globular Clusters could be obtained by comparing the
colour-temperature relation obtained using temperatures from a reddening-free
indicator (
)
with that given by standard colour-temperature
calibrations. The main difficulty in such derivations is the large errors in
temperatures for individual stars due to uncertainties on the removal of
instrumental signature for each individual star. The large multiplexing
opportunity offered by FLAMES at VLT2 allowed us to obtain spectra centred on
at a resolution of R = 6000 and 5<S/N<50 for 120 stars near the
turn-off of NGC 6752 with GIRAFFE from a single 1300 seconds exposure. This set
of spectra was used to derive effective temperatures from fittings of
profiles with typical errors of about
200 K and reddening
estimates with individual errors of 0.05 mag. Averaging all individual
reddenings, a high precision reddening estimate was obtained for the cluster:
E(B-V)=0.046
0.005. The same exposure provided UVES spectra of seven
stars near the red giant branch bump at a resolution of 40 000 and
20<S/N<40. These spectra, combined with temperatures from colours (corrected
for our high precision reddening value), provided Fe abundances with internal
errors of 0.026 dex and with average metallicity [Fe/H] = -1.48
0.01
0.06 dex (random + systematic). Abundances were obtained for several other
elements, allowing e.g. an accurate estimate of the ratio between the
-elements and Fe ([
/Fe] = +0.27
0.01). The O-Na anticorrelation is evident from our UVES data, in agreement with past results.
This analysis shows the power of FLAMES for analysis of globular clusters: the
accurate reddenings and metal abundances obtained by a procedure like that
described here, combined with distance determinations from cluster dynamics or
main sequence fitting, and high quality colour-magnitude diagrams, could allow
derivation of ages with errors below 1 Gyr for individual globular clusters.
Key words: stars: abundances - stars: evolution - stars: Population II - Galaxy: globular clusters: general - Galaxy: formation
Accurate estimates of the absolute ages of (the oldest) Globular Clusters (GCs), coupled with determinations of the Hubble constants H0 from the spectrum of fluctuations of the microwave background determined from either the WMAP experiment (Spergel et al. 2003), or from the HST Key Project (Freedman et al. 2001) may provide a stringent lower limit to the age of the Universe and constrains the exponent w of the equation of state for dark energy (Jimenez et al. 2003), independently of type Ia SN observations (see e.g. Gratton et al. 2003a).
Furthermore, within the framework of a standard CDM model, where the
age of the Universe is accurately fixed at 13.7
0.2 Gyr by the WMAP results (Spergel et al. 2003), the age of GCs can be used to constrain the
epoch of formation of the Galaxy, linking the local Universe to the distant
one (see Carretta et al. 2000; Gratton et al. 2003a).
Relative ages are fundamental to describing the early history of our Galaxy. In this framework it should be noted that galactic GCs divide into two main groups: halo and thick disk (or bulge) GCs (Zinn 1985). The differential ages method suggests that these two groups might have ages that differ by about 2 Gyr (Rosenberg et al. 1999), a result that seems supported by absolute ages (Gratton et al. 2003a). This has implications for both cosmology, where only the oldest GCs are of interest, and galactic evolution. Observations of GC systems in other galaxies suggest a link between GC formation and strong dynamical interactions (Peebles & Dicke 1968; Schweizer & Seitzer 1993). The oldest group of GCs might then be related to the very early phases of the galactic collapse, while the second one may instead trace a later accretion event (see Freeman & Bland-Hawthorn 2002), possibly related to the end of the thick disk phase indicated by chemistry (Gratton et al. 1996; Gratton et al. 2000; Fuhrmann 1998).
An important goal is then to derive absolute ages with internal errors of 1 Gyr for a sample of GCs. Ages for GCs with such small errors may be derived
using only the luminosity of the turn-off (TO), which in turn requires
accurate distances, with errors <
.
In perspective, most accurate and
robust distances (error <2%) for a few GCs will be obtained using
geometrical methods (Piotto et al. 2004). At present, distances with errors of
3-5% can be obtained for a larger sample of GCs using the Main Sequence
Fitting Method (MSFM), exploiting local subdwarfs as standard candles (see
Gratton et al. 1997; Pont et al. 1998; Carretta et al. 2000; Gratton et al. 2003a). Main sources of errors in MSFM are possible systematic differences in
reddenings and metallicities between field and GC stars, while accurate
initial He abundances in GCs have been determined by Cassisi et al. (2003; see
also Salaris et al. 2004). Both of them can be reduced to within the required
accuracy if a reddening-free temperature indicator is used for both field and
GC stars of similar evolutionary phases. Analysis of the results of the ESO LP 165.L-0263 (Gratton et al. 2003a) relative to three GCs (NGC 6397, NGC 6752 and 47 Tuc), spanning almost the total metallicity range of galactic GCs, showed that this approach may provide reddenings accurate to
mag, metallicities accurate to
0.04 dex, distances accurate to 4%, and ages with errors of about
1 Gyr. Also, geometrical distances,
which determines the true distance modulus toward a cluster, will take
advantage of accurate reddening and metallicity determinations, since apparent
distance moduli are required to derive ages.
Note that here we are only interested in relative reddening and metallicity
determinations: the adopted scale may be tied to that of field stars
exploiting the clusters observed within the LP 165.L-0263 (47 Tuc, NGC 6397,
and NGC 6752: Gratton et al. 2003a). The same temperature indicator may be
adopted (H profile).
In this paper we describe a pilot program on NGC 6752, which exploits the
multiplexing capabilities of FLAMES, the VLT multifibre facility (Pasquini et al. 2002). The large number of spectra that could be obtained using GIRAFFE
allowed proper reduction of the major source of errors in temperatures
derived from H:
flat fielding. On the other hand, low resolution and S/N were not too critical in such observations, allowing use of faint turn-off
stars. The simultaneous acquisition of spectra of a few red giants with UVES
additionally allowed accurate determination of the chemical composition.
For this purpose we preferred to use relatively warm stars, for which the
analysis is expected to be robust. This procedure's success suggests the
usefulness of an extensive program on other globular clusters. Ages within
1 Gyr are now fully within reach for a substantial sample in both Zinn's groups.
Table 1: Summary of observations.
Data used in this paper are based on a single 1300 s exposure obtained
on June 24th, 2004, with FLAMES at Kueyen (=VLT2) in service mode. The
observations were obtained at a rather large airmass (about 1.8) and in very
good seeing conditions (FWHM = 0
62 at zenith and 5000 Å) (see also
Table 1).
120 fibres feeding the GIRAFFE spectrograph were centred on stars slightly
brighter than the turn-off of NGC 6752 in the magnitude range
16.7<V<17.2.
Stars were carefully selected from high quality photometric observations obtained with the Wide-Field Imager (WFI) at the 2.2 m ESO-MPI telescope (total field of view of 34
33 arcmin2). For a detailed representation of the data reduction and calibration of this data set we refer the reader to Momany et al. (2004).
The astrometric calibration of NGC 6752 reference images employed over 7000 stars from the GSC2.2 catalogue (Loomis et al. 2004) using the IRAF MSCRED
package. The internal accuracy of the astrometry has been estimated to be
about 0.15 arcsec, well within the FLAMES requirements (0.2 arcsec). To
further confirm fulfilment of the requirements, we matched our
astrometrically calibrated NGC 6752 catalogue with that of UCAC2 (Zacharias et al. 2004) and estimated the positional residuals for 550 stars in common. The
residuals show a Gaussian distribution with an rms of
0.05 arcsec in both coordinates.
Only uncrowded stars were considered, that is, stars not showing any companion
brighter than
mag within 2.5 arcsec, or brighter than
mag within 10 arcsec. The targets were selected to lie
close to the cluster mean loci in the colour-magnitude diagram. A posteriori,
radial velocities confirmed membership of all but two of the observed stars.
Eight fibres were additionally used to monitor sky background; they were
pointed toward carefully selected empty sky regions.
The GIRAFFE spectrograph was used with the LR06 grating; the spectra cover the
wavelength range 6400-7100 Å at a resolution of about
.
Pixel-to-pixel S/N of the spectra (measured from the scatter of individual
spectral points in the wavelength range 6660-6670 Å) ranged from 5 to 50,
with typical values around 20. The S/N values were generally lower for stars
in the outer regions of the cluster, most likely because these stars were not
well centred on the fibre heads. This can be attributed to the effects of
differential refraction at the rather large airmass of observation.
The average of the eight sky spectra were used to subtract telluric emissions
(in particular, emission in H)
from GIRAFFE spectra. Appropriate
scaling factors were evaluated to take into account the transmission of
individual fibres.
Seven fibres feeding the UVES spectrograph were centred on stars close to the RGB bump ( 13.2<V<14.2), while one was dedicated to the sky. The spectra cover the wavelength range 4700-6900 Å and have 20<S/N<40.
The two sets of spectra were reduced using the dedicated FLAMES pipelines (BLDRS Python software 0.5.3 version for the GIRAFFE spectra; uves/2.1.1 version for UVES spectra). We found that this UVES pipeline does not accurately subtract the background between orders in the green-yellow part of the spectra. Only a few lines measured on these portions of the spectra were considered in the present analysis.
The following procedure was used to derive accurate temperatures from the
H profile. First, instrumental fluxes within 10 narrow bands of
5 Å width in the region including
were measured on the GIRAFFE spectra by integrating the observed spectra, after shifting them in wavelength
for the geocentric radial velocity of each star. Cosmic ray hits were removed
before evaluating the fluxes. The list of the bands used is given in Table 2.
The fluxes measured in each band were then normalized to a pseudocontinuum given by a straight line connecting the average of the first two bands with the average of the two last bands. The normalized fluxes for all the stars observed with GIRAFFE are given in Table 3 (available only in electronic form at the CDS). The second column of this table gives also the S/N ratio for each spectrum, computed from the spectral region 6660-6670 Å, where there is no significant feature.
Since the
profiles are expected to be fairly symmetric and since
bands are defined symmetrically with respect to the line center, errors in
these normalized fluxes can be obtained by comparing fluxes measured on the
blue and red side of the
.
Eliminating a few outliers, the
comparisons are as follows:
![]() |
(1) |
![]() |
(2) |
Table 2:
Definition of the bands used to derive
's.
Individual heliocentric radial velocities measured by the FLAMES pipeline are
presented in the third column of Table 3. In a few cases (7 stars out of 120), these measures of the radial velocities by the automatic
routine in the pipeline were obviously wrong, perhaps due to the strong
telluric signal present in this wavelength range. Radial velocities for these
stars were measured by fitting H,
and zeroing the radial velocity on
the telluric bands consistently with the other stars. All but one (29049) of
the stars appear to be members of the cluster on the basis of radial velocity;
the velocity for star 39462 is discrepant, too, but this spectrum has very low S/N, so we suspend judgement about it. However we did exclude these two stars
from our estimates of reddening of NGC 6752. The mean radial velocity is
-32.0
0.6 km s-1, with an rms scatter of individual values of
6.0 km s-1 (118 stars), with no obvious correlation of the spread with
neither S/N nor distance from the cluster centre.
The average value of the radial velocity agrees very well with those estimated
by Webbink (1988: -32.2
3.2 km s-1) and Dubath et al. (1997:
-32.0
1.6 km s-1), while it is slightly larger than that given by
Rutledge et al. (1997: -27.4
2.7 km s-1). For comparison, the seven stars observed with UVES provided a slightly lower average velocity (-23.8
2.1 km s-1, see Table 6).
Errors in the GIRAFFE radial velocities should be roughly 5-6 km s-1, as given by the FLAMES pipeline. For comparison, the radial velocity
error is expected to be roughly
km s-1 when
using the formula by Landman et al. (1982) and assuming that the radial
velocity signal is given by H
alone. For the typical S/N of the
program spectra, errors are then expected in the range 0.7-6 km s-1.
Given these large uncertainties on the errors attached to these radial
velocities, they are of little use in estimating the internal velocity
dispersion in NGC 6752.
![]() |
Figure 1:
Top panel: average H![]() ![]() |
Open with DEXTER |
Effective temperatures were derived by comparing the normalized average fluxes
in bands 5 and 6 with those expected from Kurucz (1992) model atmospheres
(with the overshooting option switched off) of different temperatures.
Gravities and metal abundances adopted for these model atmospheres are those
appropriate for the program stars. The
absorption and broadening
were modelled using the same assumptions as Castelli et al. (1997). The
theoretical profiles were further broadened by convolution with a Gaussian
profile mimicking the instrumental profile.
To put these effective temperature derivations in a clearer perspective, we
show in Fig. 1 the average spectrum for the stars with Giraffe
spectra with S/N>25 in the region around H.
Superimposed are
synthetic spectra computed for gravities and metal abundances appropriate for
the program stars and various values of
(=5600, 5800, 6000,
6200, 6400, 6600 K). The limits of the bands used are also shown as vertical
marks.
Values for the
's are given in Table 4 (available
only in electronic form at the CDS).
Reddening estimates for individual stars can be derived by comparing the observed colours with those predicted from the reddening-free temperatures and an appropriate colour-temperature relation.
The B-V colours were given by the WFI photometry (Momany et al. 2004). They
were corrected blueward by 0.020 mag to put them on the same scale of Thompson
et al. (1999) used in Gratton et al. (2003a). Based on the scatter around the
mean relation, we expect that errors in colours for individual stars are of
0.025 mag at this magnitude. We found clear correlations between the
's derived from
,
the B-V colours, and the V magnitudes within our sample (see Fig. 2).
Reddenings toward NGC 6752 can finally be evaluated by comparing the observed colour-temperature relation with that expected from models. Each individual star provided a reddening estimate, as listed in the last column of Table 4. No clear trend in these reddening estimates with e.g. V magnitude or location on the field could be discerned. The individual values were then averaged together to provide a best estimate. Values obtained in this way are listed in Table 5. For comparison, the value obtained by using a similar procedure from spectra of 20 stars taken at higher resolution (Gratton et al. 2003a) is also given, as well as the values from the compilation by Harris (1996) and from the reddening maps of Schlegel et al. (1998).
![]() |
Figure 2: Top panel: effective Temperature - B-V colour diagram for the stars observed with GIRAFFE; Left panel: effective Temperature - V magnitude diagram for the same stars; right panel: B-V colour - V magnitude diagram for the same stars. |
Open with DEXTER |
There are several aspects in this procedure that may introduce systematic
errors in our reddening derivations: systematic uncertainties in flat
fielding, inappropriate modelling of H,
errors in the model
atmospheres, photometric errors etc. However, we recall here that what
matters for the age derivations are not the absolute values of the reddenings,
but instead the relative values between cluster and field stars. The good
agreement with the determination of Gratton et al. (2003a), which is on a
scale consistent with that adopted for the field subdwarfs, supports the
present technique.
Table 5: Various reddening estimates.
Table 6: Photometry and spectrum data for stars observed with UVES.
Table 6 gives the main parameters for the stars observed with UVES. The colours have been corrected blueward by 0.008 mag to put them onto the same scale of Thompson et al. (1999). Note that this correction is slightly different from that found appropriate for the TO-stars, which suggests the presence of a colour term in one of the two photometry.
The EWs were measured on the spectra using the ROSA code (Gratton 1988; see
Table 7) with Gaussian fittings to the measured profiles, which
exploit a linear relation between EWs and FWHM of the lines, derived from a
subset of lines characterized by cleaner profiles. Since the observed stars
span a very limited parameter range, errors in these EWs can be computed by
comparing values derived from individual stars with the average value for the
whole sample. Typical errors obtained using this procedure are listed in the
last column of Table 6. They are roughly reproduced by the
formula
mÅ. Considering the resolution and
sampling of the spectra, the errors in the EWs are in agreement with
expectations based on photon noise statistics (Cayrel 1988). Finally, we
notice that due to the problems in background subtraction in the green-yellow
part of the spectra, only lines with wavelength >5900 Å were considered,
save for Na, Mg, and Si, for which also lines in the 5600-5750 Å region were considered.
Effective temperatures were derived from dereddened B-V colour using the
calibration by Alonso et al. (1999): we interpolated the values at
[Fe/H] = -1.5 from the tables for [Fe/H] = -1 and -2. Rather than directly
using the individual stellar colours, we preferred to use the colours of
the mean loci at the same magnitudes of the program stars. Individual stellar
colours would have produced larger individual errors (0.014 mag,
corresponding to 34 K, rather than <
K with the procedure adopted
here: see below).
Internal uncertainties in these temperatures can be obtained by considering
the errors in the V magnitudes (<0.05 mag) and the slope of the magnitude
temperature relationship, which is 193 K/mag in the range of interest for the
program stars. We get internal uncertainties of <10 K, corresponding to
about 0.011 dex in [Fe/H]. Systematic errors are larger. The uncertainty in
the reddening (0.005 mag) multiplied for the slope of the
colour-temperature relation (about -2400 K/mag) yields a systematic error of
12 K, which is about 0.013 dex in [Fe/H]. Much larger is the uncertainty
in the adopted temperature scale, most likely in the range 50-100 K,
producing possible errors in the Fe abundances in the range 0.055-0.11 dex.
Surface gravities were obtained from the location of the stars in the
colour-magnitude diagram. This procedure requires assumptions about the
distance modulus (from Gratton et al. 2003a), the bolometric corrections (from
Alonso et al. 1999), and about the masses (we assumed a mass of 0.9 ,
close to the value given by isochrones fittings). Uncertainties in these
gravities are small; we estimate a total error of about 0.15 dex, dominated by
systematic effects in the temperature scale.
Microturbulent velocities
were determined by eliminating trends in the
relation between expected line strength and abundances (see Magain 1984). To
estimate errors in these values we notice that we found that the error in the
EWs contributes 57% of the variance of the errors in the abundances for
individual lines. Given the typical uncertainties in the slope of expected
line strength vs. abundances, this implies an expected random error in the
microturbulent velocities of
0.09 km s-1. This value coincides with
the star-to-star scatter in microturbulent velocities.
Finally, model metal abundances were set in agreement with derived Fe abundance. The adopted model atmosphere parameters are listed in Table 8.
Table 8: Atmospheric parameters for stars observed with UVES.
Table 9: Iron abundances for stars observed with UVES.
Individual [Fe/H] values are listed in Table 9, as well as
averages over the whole sample. Reference solar abundances are as in Gratton
et al. (2003b). The average Fe abundance from all stars is [Fe/H] = -1.48
0.02 (error of the mean) with an rms scatter of 0.038 dex from 7 stars.
If we consider only the five stars with S/N>25, we have: [Fe/H] = -1.48
0.01 with an rms scatter of 0.028 dex.
Table 10: Uncertainties in Fe abundances for stars observed with UVES.
There is a small offset of 0.07 dex between abundances given by neutral and
singly ionized Fe I lines. This might be attributed to the use of very few
lines for Fe II, but might also indicate some errors (0.15 dex) in the
surface gravities, as well as a systematic error of about 50 K in the
effective temperatures. For comparison, Gratton et al. (2001) found
[Fe/H] = -1.44 for neutral iron, and -1.55 for singly ionized iron.
Table 10 lists the impact of various uncertainties on the
derived Fe abundances. Variations in parameters of the model atmospheres
(effective temperatures
,
surface gravities
,
model metal
abundances [A/H], microturbulent velocities
)
were obtained by changing
each of the parameters at a time. The second column gives the variation of the
parameter used to estimate the changes in the abundances from neutral (Col. 3) and singly ionized (Col. 4) Fe lines. Columns 5 and 6 give the random
(i.e. appropriate to each star) and systematic (scale errors for all stars)
uncertainties in the various parameters, while Cols. 7 and 8 give the
corresponding errors in the Fe abundances. The last row gives total errors,
which have been obtained by combining errors due to the various parameters.
The observed star-to-star scatter in Fe abundances is very small, in particular if only higher quality (S/N>25) spectra are considered. In spite of this, one may wonder if there is some evidence for real star-to-star scatter in the Fe abundances, or at least put some upper limit on this scatter, even though the sample of stars observed with UVES is not extensive. In Sect. 5.3 we have seen that the expected star-to-star scatter in Fe abundances due to the adopted temperatures is only 0.008 dex. More relevant is the error due to the microturbulent velocities.
To evaluate this source of error, we first note that considering only spectra
with S/N>25, the error in abundances from individual lines from each
spectrum is 0.068 dex, while the line-to-line rms scatter of the average
abundances from the 5 spectra is 0.097 dex. This indicates that only part of
the variance in the internal abundance errors is due to random errors in the
EWs, variable from star-to-star. The remaining contribution can be attributed
to systematic errors proper to each line (oscillator strengths, blends and
systematic effects on positioning of the continuum level). Given these facts,
the typical internal error in the abundances due to EWs can be estimated to be
dex. Also, in the same way we may distribute the
measured errors in the microturbulent velocities (determined from the
uncertainty in the slope of the expected line strength vs.
abundance fit) between random (i.e. star-to-star variable) and
systematic (i.e. constant throughout the analysis of all stars) errors. Only
the random (star-to-star variable) error should be considered when discussing
the star-to-star abundance variations. By combining the various sources of
random errors quadratically, we get a prediction of 0.026 dex for the
star-to-star spread in the Fe abundances. This compares very well with the
measured star-to-star scatter of 0.028 dex. The conclusion is that there is
very scarce evidence for an intrinsic star-to-star scatter in the abundances;
a one-sided 1-
upper limit is 0.017 dex.
Table 11: Element-to-element abundance ratios for individual stars, averaged over the whole sample. For comparison, we also give the values found by Gratton et al. (2003a) and James et al. (2004).
The O-Na anticorrelation in NGC 6752 is well known from previous observations of both TO and subgiant stars (Gratton et al. 2001; Carretta et al. 2005), as well as of red giants (Carretta 1994; Norris & Da Costa 1995; Yong et al. 2003). It is fully confirmed by the present data for RGB bump stars (see Fig. 3). Figure 3 also collects all data available up to now, which clearly shows how an extensive O-Na anticorrelation can be seen along all evolutionary phases.
![]() |
Figure 3: [Na/Fe] ratio as a function of [O/Fe] for stars in NGC 6752. Red filled circles are our RGB bump stars in the present study. Green-filled and open circles are subgiant and turn-off stars, respectively, from Carretta et al. (2005). Literature data are as follow: green diamonds with crosses inside are bright red giants from the extensive study by Yong et al. (2003), and open green triangles are red giant stars from Norris & Da Costa (1995; 6 stars) and Carretta (1994; 4 stars). |
Open with DEXTER |
Table 11 lists the average abundances obtained for various
elements. Also in this case, solar abundances were as in Gratton et al.
(2003a). In the same table we also compare the abundances obtained in this
paper with the analysis of Gratton et al. (2003a) and James et al. (2004). The
comparison is very good for the best determined elements: we found clear
overabundances of the -elements ([
/Fe] = +0.27
0.01), a
small deficiency of Ni, and slight overabundance of Ba. We notice that NGC 6752
closely traces the composition of the dissipative component considered by
Gratton et al. (2003b), in agreement with its kinematics (Dinescu et al.
1999). The overabundance of
-elements looks quite similar to those of
other globular clusters (see e.g. Gratton et al. 2004).
We have shown that a single 1300 s exposure with FLAMES at VLT2 can
provide accurate estimates of the reddening toward NGC 6752, as well as of the
chemical composition of the cluster (errors of 0.005 mag and of 0.02 dex,
respectively). Similar analyses may provide results of comparable accuracy for
other globular clusters, too, on a uniform scale. While results of similar accuracy have already been obtained for a
few clusters (including NGC 6752), use of FLAMES allows such high accuracy with
much less observing time (about a factor of 20). An extensive program over a
large number of clusters may reduce errors in age determination for those
clusters for which accurate distances could be obtained from either the main
sequence fitting method, or even better from dynamical methods.
Acknowledgements
This research was funded by PRIN 2003029437 "Continuità e discontinuità nella formazione della nostra Galassia'' by the Italian MIUR.