E. Carretta1 - R. G. Gratton1 - A. Bragaglia2 - P. Bonifacio3 - L. Pasquini4
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 -
INAF - Osservatorio Astronomico di Trieste, via Tiepolo 11, 34131 Trieste, Italy
4 -
ESO - European Southern Observatory
Received 22 September 2003 / Accepted 5 November 2003
Abstract
We used the UVES spectrograph on Kueyen (VLT UT2) to perform
abundance analysis of stars at the base of the giant branch (RGB) and at the
main sequence turn-off in the globular cluster 47 Tuc (NGC 104). High
dispersion spectra (
)
for 3 dwarfs and 9 subgiants were
analyzed. We found an overall [Fe/H] value of -0.67
0.01
0.04,
rms = 0.05 dex. The mean value obtained from stars at the base of the RGB is
virtually the same as that obtained from the dwarfs. Na and O abundances show a
star-to-star scatter and are anti-correlated, not unlike
previous findings
in more metal-poor clusters. The extension of the anticorrelation is not
as extreme as found in other clusters, but it is clearly present in unevolved
stars. Al abundances do not show a significant spread; this could be indicative
that p-captures in high temperature regions are less likely to occur in
metal-rich clusters. The
[
/Fe] ratios suggest
a slight excess of Ti with respect to field
stars of similar metallicities.
We found an enhanced odd-even
effect for Fe-group elements. In particular,
the [Mn/Fe] ratio is almost 0.2 dex
deficient with respect to the mean value of the dissipative component of field
stars in the solar neighborhood. While the theoretical and observational
framework for Mn is still poorly understood, the chemistry of 47 Tuc, coupled
with its kinematics and its somewhat younger age, might give hints favouring
an origin in a formely independent, possibly larger, subsystem of our Galaxy.
Key words: stars: abundances - stars: atmospheres - stars: evolution - stars: population II - Galaxy: globular clusters: general - Galaxy: globular clusters: individual: 47 Tuc
The galactic globular clusters (GCs) are important for a large set of astrophysical problems, ranging from stellar evolution to cosmology. Their chemical compositon is a fundamental probe of the early chemical evolution in the galactic environment. However, before using globular clusters as tracers of the early chemical conditions in the Galaxy, one has to disentangle the original compositions from any effect which could have resulted in a pollution of their photospheres. Possible effects include (i) possible self-enrichment of the protoclouds from which the GCs originate (Cayrel 1986; Brown et al. 1991, 1995); (ii) changes in abundances derived from the internal evolution of the stars (both canonical or non-canonical mixing, see Kraft 1994; Gratton et al. 2000b); and (iii) modifications of the composition due to a possible second generation of stars born from the low velocity wind of intermediate mass globular cluster stars (Parmentier et al. 1999; Parmentier & Gilmore 2001).
A growing body of evidence shows that the usual assumption of globular cluster stars being exactly coeval and chemically homogeneous (the basis of their classification as a pure simple stellar population) is unreliable. Useful records of the first phases of formation and early evolution of GCs are hidden in the photospheric abundances of stars presently observed. Once correctly interpreted, they may provide an understanding of the origin and evolution of GCs which may not be achieved by other means.
The most noticeable of these anomalies are those involving the light element C, N, O (for exhaustive reviews see e.g. Smith 1987; Smith & Mateo 1990; Kraft 1994 and references therein), as detected from molecular abundances of CN, CH, NH and OH, often accompanied by correlated anomalies in atomic species such as O, Na, Al and Mg (see the survey by the Lick-Texas group).
The observed abundance pattern may be explained by H-burning at high temperatures, where reactions and products of the CNO, NeNa and MgAl cycles occur. Denisenkov & Denisenkova (1990), Langer et al. (1993) have shown that production of 23Na and 27Al from proton captures on 22Ne and 25Mg is possible even in low mass stars, in regions close to the ON-burning shell.
To explain the complex observational picture, there are presently two main lines of thought. The first (see Cottrell & Da Costa 1981; Brown & Wallerstein 1992) considers primordial inhomogeneities, i.e. already existing in the gas from which the cluster stars formed, or released from the very early generation of massive stars born and evolved in GCs. This scenario is based on the belief that Na and Al cannot be synthetized in the low mass stars presently observed in GCs, so that the overabundances of these elements, and the correlated CN-enhancements, must be due to the yields of more massive stars that polluted the next generation of cluster stars.
The second scenario explains the CN and anticorrelated anomalies as due to internal deep mixing, that brings up to the surface fresh 14N from CN processes, depleting at the same time the abundance of 12C of the surface layers. The constant decrease of the isotopic ratio 12C/13C up to values near the equilibrium ratio of 3.5 (see the review by Suntzeff 1993) supports this mechanism. All stars experience a mixing phase (the standard first dredge-up) during the first ascent on the giant branch; however, to explain the fact that only some of the stars are CN-strong, one has to assume the existence of a mechanism of extra mixing.
Following the model of Sweigart & Mengel (1979), this additional mechanism
could be represented by
meridional circulation currents driven by core rotation. To explain
the observations, these currents should be able
to reach the region of the ON process, bringing up to the surface processed
N-rich and O-poor material that contributes to the CN enhancements even in the
case of large C depletions.
However, in the model of Sweigart & Mengel (1979) a gradient of mean molecular
weight does inhibit the formation of circulation currents (Tassoul & Tassoul
1984), and so the additional mixing, below a critical level (log
2.2) on the RGB. This is the point of maximum penetration of the
convective envelope, succesively crossed by the H-burning shell propagating
outward (the so-called RGB-bump).
Observations show the opposite: anomalous CN abundances exist at considerably lower luminosities than predicted. Furthermore, a radial dependence in the distributions of CN-strong and CN-weak stars (the first being more centrally concentrated) seems to be present in 47 Tuc and possibly in other clusters (Norris 1987), and this is not easily explained by the mixing scenario.
Recently we definitively settled this issue by finding a clear Na-O anticorrelation among unevolved (turn-off) stars in the globular cluster NGC 6752 (Gratton et al. 2001, hereafter Paper I). These stars do not meet the constraints of high central temperature and extended envelope required to self-produce and bring to the surface the synthetized proton-capture elements. The same seems to hold also in the more metal-poor cluster NGC 6397 (Carretta et al. 2003a, in preparation).
The Na-O anticorrelation is apparently limited to the dense environments of GCs: there are almost no stars showing CN-strong or Na-rich composition in the field, where the chances for pollution are much less than in the denser GCs. Gratton et al. (2000b) convincingly showed that the abundances of elements heavier than CN (e.g. Na and O) are unaltered in field stars of whatever evolutionary stages, from pre-main sequence up to the horizontal branch (HB) phase.
There are several reasons that make 47 Tuc a preferred target for studying the complex issue of chemical anomalies in globular cluster stars:
The favourable location of the ESO VLT complex in the southern hemisphere, coupled with the high efficiency of the UVES spectrograph, now allow us to obtain high-resolution spectra at a S/N adequate for a precise, detailed abundance analysis of several elements in 47 Tuc dwarfs.
![]() |
Figure 1: Positions of program stars (large symbols) on the Strömgren colour-magnitude diagram of 47 Tuc. Photometry is from Grundahl et al. (1999). |
| Open with DEXTER | |
Observations were carried out in three different runs (September 2000; August
and October 2001) with the spectrograph UVES at the second 8m Unit Telescope of
the VLT (Paranal ESO
Observatory); each run consisted of 6 nights assigned to the ESO Large
Program 165.L-0263.
Due to poor weather conditions in the September 2000 run, we were not able to
complete the acquisition of the spectra of 47 Tuc. Bad weather/seeing also
plagued the August 2001 and partially the October 2001 runs. In total, 5 candidate turn-off stars and 9 candidate subgiant stars were observed.
However, one of the turn-off stars turned out to be non-member, both from radial
velocity (see Lucatello & Gratton 2003) and elemental abundances, and
for another we were able to gather only a low-S/N (
15) spectrum, useless
apart from the membership determination.
Hence we have useful observational material (high-resolution, high S/N spectra)
only for 9 stars at the base of the RGB and 3 stars at the main
sequence turn-off in 47 Tuc.
Following Paper I, we tried to observe stars likely having both strong and weak CN bands, selecting objects on the basis of Strömgren photometry, using the index c1 (increasing with the strength of CN bands, Grundahl et al. 1999).
Selected stars are indicated in Fig. 1, and their relevant data are listed in Table 1.
Table 1: Data for observed stars in 47 Tuc (identifications from Grundahl et al. 1999).
As for NGC 6397 and NGC 6752 (see Paper I), data were acquired using the
dichroic beamsplitter #2 at UVES. The CD2, centered at 420 nm (with
spectral coverage
356-484 nm), was used in the blue arm of the
spectrograph; and the CD4 centered at 750 nm (covering
555-946 nm) was adopted in the red arm. Slit length was always 8 arcsec;
the slit width was mostly set at 1 arcsec (corresponding to a resolution of 43 000 at the order centers). However, this value was slightly modified downward
or upward in a few cases, according to the seeing conditions.
Typical exposure times were 2 h for each subgiant and about 4 h for each turn-off star, split in multiple exposures. Values of the S/N per pixel (about 5 per resolution element) measured at about 6700 Å are listed in Table 1.
Following the procedure adopted in Paper I, we compared effective temperatures
from observed colours (both Johnson B-V and Strömgren b-y were used) with
spectroscopic temperatures derived from comparison of Balmer line profiles (namely H
)
with those obtained from spectral synthesis.
An example of the fit is shown in Fig. 2 for star 1081.
![]() |
Figure 2:
Derivation of temperature from the H |
| Open with DEXTER | |
Since (i) these last temperatures are reddening-free and (ii) they are derived
following a uniform procedure for field stars (see Paper I), this
approach allows us to compare directly our derived abundances for 47 Tuc with
those obtained for field stars. This is a key requisite for accurate distances
from the Main Sequence Fitting Method (see Gratton et al. 2003a)
.
However, some differences do exist with respect to the analysis performed for stars in NGC 6397 and NGC 6752.
Estimates of the microturbulent velocity
for each star were derived, as
usual, by eliminating the trend of abundances with expected line strengths;
for this, we used abundances of Fe I, for which a large number of lines were
measured in each star. Unlike Paper I, we found that, for stars at the base of
the giant branch, the star-to-star scatter in Fe abundances was reduced if
individual values of microturbulent velocity were used.
Finally, the overall model metallicity [A/H] was chosen as that of the model
atmospheres extracted from the grid of ATLAS models
with the overshooting option switched off computed by Castelli
, that best reproduce the measured equivalent widths (EW).
Final adopted atmospheric parameters are listed in Table 2.
Table 2: Adopted atmospheric parameters and derived Iron abundances for observed stars in 47 Tuc.
| |
Figure 3: Comparison between EWs measured on the spectra of two dwarfs in 47 Tuc ( upper panel; stars 1012 and 1081) and two subgiants ( lower panel, stars 433 and 435). The dashed lines indicate equality, and the rms scatters about the average differences amount to 6.8 mÅ in the case of subgiants and 5.5 mÅ for TO-stars, in this example. |
| Open with DEXTER | |
Automatic equivalent width measurements of spectral lines were performed exploiting a recent version of the spectrum analysis package developed in Padova and partially described in Bragaglia et al. (2001) and Carretta et al. (2002). The number of useful lines depends on the metallicity of star and on the S/N. The metal abundance being the same for dwarfs and early subgiants in 47 Tuc, we were able to measure similar numbers of Iron lines for the two groups, i.e. from 50 to 60 Fe I lines for the SGB and from 30 to 40 Fe I lines for dwarfs, and about 10 Fe II lines in each group, securing a good database with statistical significance for the following analysis.
From the classical formula by Cayrel (1989) widely used to give an estimate of the error in measured EW as a function of the full width half-maximum, we expect an error of 3.8 mÅ, due to the characteristics of our spectra and their typical S/N ratios (see Table 1). This value can be compared with observations.
An empirical estimate of internal errors in the equivalent widths can be obtained by comparing stars of similar physical status and with spectra of similar quality. We performed this exercise by cross-comparing the sets of EWs measured for subgiants (discarding from the comparison only the 2 stars with lower S/N) and turn-off dwarfs separately. From 15 comparisons of subgiants, we derived an average rms scatter of 6.4 mÅ for the average difference EW(star1)-EW(star2), based typically on about 140 lines in common between each pair of stars. In the same way, the 3 inter-comparisons of the three dwarfs resulted in an average rms scatter of 7.9 mÅ. If we suppose that errors can be attributed in equal proportions to the stars, we obtain with typical errors in EWs of 4.5 mÅ for subgiants and 5.6 mÅ for the turn-off stars. Figure 3 shows two of these comparisons, one for a pair of subgiant stars and one for two dwarfs.
Comparing the observed typical errors in EWs with the expectations of the Cayrel' formula, we see that another error source is missing, giving residual uncertainties of 2.5 mÅ (quadratic sum) for SGB stars and 4.1 mÅ for TO stars. This can be likely attributed to uncertainties in the location of the continuum (neglected in the Cayrel formula); with our procedure (relationship between FWHM and central depth), and considering for simplicity a triangular shape for the lines, we estimate that errors in the (automatic) continuum tracing at a level of 1% (for subgiants) or 2% (for dwarfs) can explain the residual discrepancy.
Table 3 shows the sensitivity of the derived abundances to variations in the adopted atmospheric parameters for Fe and other elements; this is obtained by re-iterating the analysis while varying each time only one of the parameters. This exercise was repeated for a subgiant (star 435) and a dwarf (star 1081).
Note that in some cases the entities of variations listed in
Table 3 are likely to be overestimates of the actual errors in
the atmospheric parameters. For instance, uncertainties in temperatures were
estimated in Paper I to be
90 K for the dwarfs and
60 K
for the subgiants. The same values are applicable here.
Errors in the surface gravities can be estimated by taking into account
uncertainties in distance moduli and bolometric corrections (affecting
luminosities:
0.06 dex), temperatures (
0.03 dex), and masses
(
0.04 dex). If we sum in quadrature all these contributions, we
find that the adopted gravities have internal errors not larger than 0.08 dex.
Typical errors in the overall metal abundance are less than 0.10 dex, and have a small effect in the final error budget.
To estimate the proper error bars in the microturbulent velocity values, we
started from the original atmospheric parameters adopted for the subgiant 435 and
the dwarf 1081; then we used the same set of Fe lines to repeat the analysis
changing the
value until the 1
value from the slope of the
abundance/line strength relation was reached. A simple comparison allows us to
give an estimate of 1
internal errors associated to
:
they are
about 0.7 km s-1 for subgiants and 0.2 km s-1 for TO stars.
Column 7 of Table 3 allows to estimate the effect of errors in the EW; this was obtained by weighting the average error from a single line (obtained separately from Fe for the 9 subgiants and the 3 dwarfs) with the square root of the mean number of lines (listed in Col. 6 of the Table) measured for each element in subgiants and dwarf respectively. Finally, the last column reports the total error bars as the quadratic sum of all the contributions.
For 47 Tuc we derived from our entire sample of 12 stars an average Fe abundance of [Fe/H] =-0.67
0.01
0.04 dex, where the first is the
internal error (standard deviation of the mean) and the second represents the
systematic errors, as the weighted average of errors found for TO stars and
subgiants. As stated in Paper I, systematic errors in our case are
mainly due to uncertainties in the temperature scale, which is related to the
scale defined by field stars analyzed in the very same way of our cluster stars.
Table 3: Sensitivities of abundance ratios to errors in the atmospheric parameters and in the equivalent widths.
Note that the Fe solar abundances from the reference analysis made
using the solar model from the Kurucz (1995) grid, with no overshooting, are
slightly different from those used in Paper I, due to the
different treatment of collisional damping (see Gratton et al. 2003b).
The present adopted
values are
(Fe) = 7.54 for neutral iron and 7.49 for singly ionized iron.
We found no difference in the average Fe abundance between subgiant and dwarf
stars in 47 Tuc: [Fe/H] =-0.67
0.02 (9 stars) and [Fe/H] =-0.67
0.01 (3 stars), respectively. As far as Fe is concerned, 47 Tuc seems to be a very
homogeneous cluster: the rms scatter from star to star is not larger than 0.05 dex, i.e. no more than 12% for this cluster. This value can be entirely
explained by internal errors, suggesting a very small real star-to-star scatter
in Fe abundances. This scatter well agrees with the upper limit
of
0.04 dex estimated by Hesser et al. (1987) from the photometric
intrinsic width of the main sequence in 47 Tuc.
The new average metallicity agrees well (within the quoted uncertainties) with
the value derived in Carretta & Gratton (1997) from 5 red giants:
[Fe/H] =-0.70
0.03. As a further comparison, the revised cluster abundance
based on individual stellar abundances of Fe II, as given recently by Kraft &
Ivans (2003) is [Fe/H] =-0.63 dex.
At face value, the stars in this cluster share the
same amount of iron in their photospheres, irrespective of their evolutionary
status. A similar result has been obtained for other clusters as well (NGC 6397:
Gratton et al. 2001; M 5: Ramirez & Cohen 2003; M 71: Ramirez & Cohen 2002).
We can check the uncertainties in Fe abundances using the
ionization and excitation equilibria as probes, since temperature and gravities were not
derived from our spectra, but independently from colours and theoretical
isochrones (see also the discussion in Carretta et al. 2002).
In Fig. 4 we plot the differences between
abundances
obtained from neutral and singly ionized Fe lines and the slopes of the
abundances from neutral lines
(Fe I) with respect to the excitation
potential
as a function of effective temperature.
![]() |
Figure 4:
Upper panel: differences between the abundances of Fe from singly
ionized and neutral lines as a function of effective temperature for stars in
47 Tuc. The solid line is the average difference for the entire sample; the
lower dashed line is the average difference for the 3 dwarfs and the upper
dashed line represents the average difference for the subgiant stars. The difference for
the Sun, obtained from our reference analysis, is also plotted. Lower panel:
the slopes of the relationship resulting between the Fe I abundances and the
excitation potential as a function of adopted
|
| Open with DEXTER | |
From the upper panel of this figure we see that on average the lines of singly
ionized Fe give slightly larger abundances than neutral lines (mean difference
+0.08 dex, rms scatter
,
12 stars). The existence of a trend
with temperature is questionable if we take into consideration the error bars.
The largest differences are found among the subgiants, whereas the average
difference for the 3 dwarfs is opposite in sign (
dex), and could as well be a reflection of the same
difference in our solar reference values (-0.05 dex, see above and
Fig. 4, upper panel).
Do we have to take into account possible departures from the LTE assumption?
The dominating effect being overionization, this could explain the difference
observed among subgiants; however, the good agreement found when comparing
abundances for stars in different evolutionary phases does not support this
hypothesis. On the other hand, errors in atmospheric parameters may affect this
difference.
We started from quoted errors in temperature (
60 K and
90 K for
subgiants and dwarfs respectively), and used the approximate slope of the
isochrone in the subgiant/dwarf region
/
(about
0.05 dex/100 K) to estimate the corresponding errors in gravity. We then
read from Table 3 the changes in surface gravity and
overall metal abundance [A/H] induced by typical errors in temperature, computing the
resulting changes in Fe I and Fe II abundance.
The average observed value of
of +0.08 dex (and
more than half of the average difference for subgiants) is reproduced for
temperatures that are lower by 60 (and 90 K) than those presently adopted.
Most of the difference therefore could simply be due to a temperature scale too
low, although well within the error bars.
Further supporting evidence comes from the observed rms scatter (0.10 dex)
observed in the difference
.
From the last
Column of Table 3, the predicted scatters are 0.06 dex and
0.08 dex for subgiant and dwarfs, respectively. Again, a significant fraction
of the observed scatter is explained only by the combination of errors, both in
atmospheric parameters and measured EWs.
Finally, the lower panel of Fig. 4 confirms that our
adopted temperature scale cannot be grossly wrong, since there is no noticeable
trend of derived abundances as a function of excitation potential. The average
value of the slope
we derive (-0.019 dex/eV)
compares well with that derived from the solar analysis
(-0.023,
dex/eV) by Carretta et al. (2002), using a very
similar line list.
In the present work we will study only the abundances of the light elements O, Na, Al, the group of the
-elements and the Fe-groups elements.
Abundances of neutron capture elements in 47 Tuc (as well as in NGC 6752 and
NGC 6397) will be presented in James et al. (2003, in preparation).
The carbon and nitrogen abundances, as well as the 12C/13C isotopic
ratios, derived from spectral synthesis of CH and CN molecular bands for stars
in 47 Tuc will be presented and discussed in a forthcoming paper (Carretta et al. 2003b, in preparation).
Abundances of lithium in stars of 47 Tuc will be discussed in detail in a forthcoming paper. Here we point out that abundances of Li and Na seem to be anticorrelated among 47 Tuc dwarfs. The dwarf star with the largest equivalent width of the Li 6707 Å doublet (39.1 mÅ) is also the one with the lowest measured Na abundance ([Na/Fe] =-0.34 dex, star 1081, see below). The dwarf star with the lowest measured equivalent width of the Li doublet (14.1 mÅ) is also the dwarf with the highest Na abundance (star 975, [Na/Fe] = +0.22). In summary, our preliminary data seem to indicate that in 47 Tuc, turn-off stars with the largest Na abundance also show the largest Li depletion in their photospheres.
Table 4: Abundances of O, Na and Al in stars of 47 Tuc.
Oxygen abundances in these warm stars were derived exclusively from the
permitted near-IR triplet at 7771-75 Å. In cooler subgiants, only some
of the lines were measured by our automatic procedure, while in warmer dwarfs
all the lines of the triplet were successfully measured. Line measurements were
further checked by eye on the spectra. Only upper limits were derived for three
SGB stars.
Final abundances and upper limits are given in
Table 4, corrected for non-LTE effects as described in Gratton et al. (1999), from statistical equilibrium calculations based on
empirically calibrated collisional H I cross sections. As stated in Carretta
et al. (2000b), these corrections are on average rather small and their
adoption brings the O abundances derived from the permitted triplet in good agreement with those from the forbidden [OI] line at 6300 Å, as far as stars with
K are considered, which is clearly our case (see also the
discussion in Gratton et al. 2003b).
We carefully looked for the forbidden [OI] line at 6300 Å in the spectra of all observed subgiants. Unfortunately, due to the cluster low radial velocity, in most of stars the [OI] line falls under the telluric emission line. However, in star 201075 the radial velocity shifts the emission slightly away, so that the stellar [OI] line could be measured. Using direct integration, we obtained an EW of 10 mÅ. To evaluate the O abundance, we must take into proper account the contribution of the Ni I line at 6300.33 Å (see Allende-Prieto et al. 2001; Johanson et al. 2003). With the adopted atmospheric parameters and the derived Ni abundance (see below), the Ni line should contribute to 1 mÅ in this star. The resulting EW of 9 mÅ for the [OI] line gives an oxygen abundance of [O/Fe] =+0.34 dex, that compares well with that of +0.41 dex from the IR permitted triplet.
Table 5:
Abundance of
-elements in stars of 47 Tuc.
47 Tuc is about 0.8 dex more metal-rich than NGC 6752, studied in Paper I.
Hence, we were able to measure a number of additional Na I lines,
beside those of the strong infrared doublet at 8183-8194 Å. Lines
of the doublet at 6154-6160 Å were observed in all subgiants and several
measurements of the weak line at 4751 Å were also made.
The doublet at 5682-88 Å is not covered by our spectra.
These lines are
listed as clean in the solar analysis by Holweger (1971).
Corrections for departures from LTE were taken from Gratton et al. (1999),
using collisional cross sections calibrated as in Carretta et al. (2000b).
These corrections are not large (on average,
[Na/Fe]
- [Na/Fe]
=
,
dex for 12 stars), and they improve the agreement between
abundances derived from different indicators. For example, the average
difference in Na abundances obtained from the line at 8183 Å and
those derived as the mean from the doublet lines at 6154-60 Å is only
-0.07 dex (
,
5 stars). Results corrected for the non-LTE effects
are given in Table 4.
Table 6: Abundances of Fe-group elements in stars of 47 Tuc: Sc to Cr.
Table 7: Abundances of Fe-group elements in stars of 47 Tuc: Mn to Zn.
The strong AlI resonance lines at 3944/61 Å were not measured; they are heavily saturated and likely affected by large departures from LTE (Baumüller & Gehren 1997). Moreover, these lines are located on the wings of the strong Ca II H and K lines, in a region where the S/N ratios of our spectra is rather low. Instead, we updated the line list of Paper I (that was limited to the doublet at 6696-98 Å) including also the doublets at 7835-36 Å and 8772-73 Å, with oscillator strengths from Lambert & Warner (1968). Our derived Al abundances rest on this set of lines. The non-LTE effects are smaller for high excitation Al lines as the ones used here (see Baümuller & Gehren 1997). Since the amount of correction for the subordinate doublets at 6696/6698 Å, 7835/7836 Å and 8772/8773 Å increases toward decreasing metal abundance, we do not include any correction for the Al abundances derived for the metal-rich cluster 47 Tuc, even if a constant offset of about 0.15 dex cannot be excluded (see Table 2 in Baumuller & Gehren 1997). Al abundances are summarized in Table 4.
A number of lines for several elements produced by
-capture were
measured. Results for Mg, Si, Ca and Ti are given in Table 5,
together with the number of lines measured and the rms values obtained for
each species. Abundances from Ti lines in two different stages of ionization
allow us to check for additional evidence of departures from the LTE assumption. From Table 5 there is a
small overabundance of singly ionized Ti with respect to the abundances from
neutral lines (on average [Ti/Fe]
- [Ti/Fe]
= 0.12 dex,
with
dex, 12 stars), with a hint of differences decreasing
with increasing temperature among the subgiant stars. However, due to the
rather limited sample size and the small range in
covered, we
do not draw any further conclusions.
Lines of several elements belonging to the Fe-group (Sc, V, Cr, Mn, Ni and Zn)
were measured. Computations of corrections due to the hyperfine structure (HFS)
was performed in the relevant cases (Sc, V and Mn); references are given in
Gratton et al. (2003b).
Another test of the adopted atmospheric parameters and of the LTE assumption
is possible here since we
detected lines of both neutral and singly ionized chromium. The average
difference we found is [Cr/Fe]
- [Cr/Fe]
(
) dex (11 stars), with no trend with temperature.
Even with large scatter (mainly due to the few lines of ionized chromium
measured), this finding supports our adopted scale for the atmospheric
parameters.
A summary of abundances obtained for the Fe-group elements is given in Tables 6 and 7.
Mean abundances for individual elements are summarized in Table 8, where we used only the subgiants, to have a more homogeneous sample of stars all in the same evolutionary phase. In Col. 3, the average from available stars (Col. 2) is computed for each element ratio, and the observed star-to-star scatter is given in Col. 4. The total uncertainty, as derived from the values of Table 3, is listed in Col. 5.
Note that
can be regarded as the predicted total error
expected as a consequence of uncertainties in the adopted atmospheric
parameters, combined with errors in EW measurements. The differences give an
idea of what elements show a scatter that likely exceeds the
predicted observational errors in stars of 47 Tuc.
As one can see from this table, apart from cases in which the scatter is
large due to the small number of lines measured (like Sc II and Cr II),
two elements clearly stand out: O and Al
.
The scatter observed in Na abundances would be
much larger if the three turn-off stars were considered.
![]() |
Figure 5: [Na/Fe] ratio as a function of [O/Fe]. Panel a): filled circles are turn-off stars and open circles are subgiant in 47 Tuc from the present study. The arrows represent the stars with only upper limits in derived O abundances. Panel b): overimposed (open squares) to stars of this study are 3 red giant stars from Carretta (1994). Panel c): 10 red giants (open triangles) of the study of Sneden et al. (1994) in M 71 (NGC 6838), from the re-analysis by Carretta (1994), are added to the stars of the previous sources. |
| Open with DEXTER | |
In panel (a) of Fig. 5 we plot the abundance ratios of oxygen and sodium from the present study. The three dwarfs and the subgiants show a well defined trend, with [O/Fe] abundances anticorrelated with [Na/Fe] ones. Upper limits in O abundances derived for three other subgiants also contribute to clearly define the anticorrelation. This is the first time that the Na-O anticorrelation is detected among scarcely evolved stars in this cluster.
In panel (b) of this figure we added three red giants analyzed by Carretta (1994). O abundances for these stars were derived from spectral synthesis of the forbidden [OI] line at 6300 Å in high resolution spectra acquired with the CASPEC spectrograph at the ESO 3.6 m telescope. Na abundances include the corrections for departures from LTE, using the same prescription as in Gratton et al. (1999) and in the present study.
Table 8: Mean abundance ratios for early subgiant stars in 47 Tuc and in field stars of similar metallicities.
At face value, and despite the small sample size, there seems to be no significant difference between the pattern shown by giants and less evolved stars. This impression is confirmed by panel (c) of Fig. 5, where we plot as open triangles results from a reanalysis (Carretta 1994) of 10 RGB stars observed by the Lick-Texas group (Sneden et al. 1994) in the globular cluster M 71 (NGC 6838), a "twin" of 47 Tuc, as far as the overall metal abundance is concerned. Again, the locus involved in the O-Na anticorrelation seems to be the same for stars in different evolutionary phases. We prefer not to use the recent results by Ramirez & Cohen (2002) to plot only Na and O values obtained through homogeneous analyses. Their Fig. 11 is in very good agreement with the present results and lead us to conclude that, whatever the phenomenon producing the Na-O anticorrelation, it is active also in metal-rich clusters, down to scarcely evolved stars.
The extension of the Na-O anticorrelation is not as large as that of more metal-poor cluster, but is undoubtly present.
![]() |
Figure 6: Run of the [Mg/Fe] ratio as a function of [Al/Fe], for stars in 47 Tuc. Symbols are as in the previous figure. |
| Open with DEXTER | |
On the other hand, elements that are produced by p-captures at higher temperatures, such as Al and Mg, do not show a clear anticorrelation, as evident in Fig. 6.
The spread in [Mg/Fe] is not astonishing, among the
subgiants, while it is somewhat larger in [Al/Fe]. However, this spread
is only marginally significant (see Table 8),
slightly exceeding the one expected from uncertainties due to
measurement errors in the EWs and to the effect of adopted atmospheric
parameters
.
On the other hand, Al is not correlated with Na.
We note that this is not at odds with the analysis of Al in stars of M 71 (Ramirez & Cohen 2002), where the existing spread could be explained purely as due to a combination of uncertainties in analysis (see their Table 7). This could suggest that products of p-captures in the MgAl-cycle are less likely to occur in metal-rich clusters like 47 Tuc and in M 71. Theoretical models seem to support this finding: in fact, the efficiency of the Hot Bottom Burning (HBB) in AGB stars (D'Antona 2003, private communication) is expected to decreases with increasing metal abundance and the temperature at the base of the convective envelope is likely not high enough to allow the burning of Mg in Al in the metal-rich clusters. The limited size of the observed samples precludes any further conclusion.
In summary, the pattern of light elements observed in unevolved turn-off stars and in stars at the base of the red giant branch in 47 Tuc seems to be not very different to that observed in other more metal-poor globular clusters. The range in Na and O abundances spanned by turn-off stars in 47 Tuc is almost comparable to that observed by Gratton et al. (2001) in NGC 6752; admittedly, the number of stars sampled in 47 Tuc is small.
The overall distribution of light elements shows that processes of proton-capture are at work. In the atmospheres of the studied stars we are seeing exactly the products of these reactions. In this case, the line of thought is the same as presented in Gratton et al. (2001): turn-off stars do not reach the temperature regime where the ON and NeNa cycles required to produce the Na-O anticorrelation are active, and moreover these stars have too small convective envelopes to have efficiently mixed the ashes of these nuclear processing up to the surface. The same conclusion holds also for subgiants.
We are seeing products of nuclear burning and dredge-up in other stars, that are now not observable, but that returned their elements to the intracluster medium or directly to the surface of presently observed stars.
Aside from changes in the pattern of light elements, we might expect that
cluster stars could significantly differ from field stars of the same
metallicity, if the protocluster was able to retain the yields of
sustained and independent chemical evolution.
Having obtained an extensive set of abundances of
-elements in star of
47 Tuc, we may test this hypothesis by looking at the pattern of elements
possibly enriched by the ejecta of core-collapse supernovae (Cayrel 1986; Brown
et al. 1991, 1995; Parmentier et al. 1999; Parmentier & Gilmore 2001).
This suggestion is inspired by the fact
that in our own Galaxy the metallicity distribution of GCs is apparently quite
different from that of field stars (e.g., Carney 1993).
Observationally, we might then expect that GCs have a larger [
/Fe] ratio than field stars of similar metal abundance. This signature is likely subtle, and may have escaped detection in the study by Carney (1996); it should be more evident among the metal-rich globular clusters,
because observations of field stars indicate that at metallicity [Fe/H]
there is already a significant contribution by type Ia SNe, which has lowered the [
/Fe] ratio from
0.4 to
0.3 (see Fig. 7).
![]() |
Figure 7:
Run of the overabundances of individual |
| Open with DEXTER | |
The elements to be considered here are those not involved in the O-Na anticorrelation (Mg, Si, Ti, and Ca), for which very accurate abundances have been obtained. Results for 47 Tuc are compared in Fig. 7 with the overabundances measured in field stars of similar metallicity in the dissipative component of our Galaxy (Gratton et al. 2003b). Values are reported in Col. 6 of Table 8, while in Col. 8 the differences (47 Tuc minus field) between the abundances are listed. The analysis is completely consistent with the present one.
The result for Mg is less significant, since, on average, Mg might be somewhat
reduced due to its participation in the Mg-Al anticorrelation phenomenon,
although there is no evidence that this is the case for 47 Tuc.
Offsets for individual elements are +0.04 dex for Si, -0.01 dex
for Ca and 0.16 dex for Ti (average of Ti I and Ti II).
On average, 47 Tuc seems to be overabundant in
-elements by about 0.06
0.04 dex with respect to halo/thick disk stars of similar metallicities.
This result is not significant.
On the other hand, this figure seems to indicate that some excess of Ti may indeed be present (a similar indication being obtained for M4 by Ivans et al. 1999), pointing out possibly that GCs produced a significant fraction of their metals. More extended observations are needed to firmly establish this point.
A much more puzzling result concerns the abundance of odd Fe-group elements. More than thirty years ago, Helfer et al. (1959) noted that among metal-poor stars an enhanced odd-even effect for Fe-group elements existed.
![]() |
Figure 8: Differences for several elements between the clusters analyzed in Paper I and in the present work, and field stars of similar metallicities (from Gratton et al. 2003b). Differences were arbitrarily shifted for display purposes. Solid lines connect Fe-group elements having consecutive atomic number. Notice the enhanced odd-even effect, in particular in 47 Tuc: Sc, V and Mn are clearly deficient with respect to Ti, Cr and Fe. |
| Open with DEXTER | |
A graphical representation of the odd-even effect in globular clusters is given in Fig. 8, where the differences of abundances with respect to field stars of similar metallicities are displayed for the three clusters analyzed so far in the Large Program 165.L-0263. In GCs, the odd-number Fe-group elements show abundances that are more deficient with respect to the even-number Fe-group elements than in field stars. This odd-even effect seems to increase with increasing overall metallicity, reaching as much as about 0.2 dex in 47 Tuc. This exceeds the observational errors.
Here we concentrate on manganese, as the run of this element has been studied for long time in several stars; moreover, it is the one showing the most marked underabundance. On the other hand, the major shortcoming is that the nucleosynthetic site of production of Mn is not yet well known, although Mn is believed to be formed mostly in type Ia SNe (e.g. Nakamura et al. 1999).
Gratton (1989) pointed out that the run of the ratio [Mn/Fe] in metal-poor
stars is mirroring in some way the run of
-elements, being almost flat
at very low metallicities, then rising with increasing [Fe/H], starting at
about [Fe/H] =-1.
The high, constant [
/Fe] ratio in metal-poor stars is commonly
interpreted as due to the delay with which elements are injected in the
interstellar medium (ISM) as a consequence of different typical lifetimes of the
producers. Massive stars ending their lifes as SN II enrich in very short times
the ISM with
-enriched yields; when the longer-lived
precursors of SN Ia
undergo thermonuclear explosion the major source of Fe is activated, thus decreasing the ratio [
/Fe]. In this framework,
Gratton (1989) suggested that the
observed run of [Mn/Fe] vs. [Fe/H]
could indicate that the fraction of manganese produced by SN II is lower
than the fraction that SN Ia are allowed to inject in the ISM.
However, the production of elements with odd atomic numbers in the Fe-group depends on the available neutron excess. Arnett (1971) proposed that the yields of Mn produced in SN II are then a function of the metallicity, through the neutron excess. In this view, the rise at higher metallicity would be due to the occurrence of more metal-rich supernovae.
This approach was recently revived by McWilliam et al. (2003; hereinafter McW03). They used new observations of giants in the galactic bulge and in the dwarf spheroidal Sagittarius to constrain the run of [Mn/Fe] at low metal abundance and at [Fe/H] values typical of the solar neighborhood. Using these two environments to pinpoint the underabundances of Mn, they claimed that manganese yields from both SN Ia and SN II are metallicity dependent.
However, since their results are based on the still few reliable high resolution spectra obtainable for these rather distant systems, we gathered Mn abundances for several kinds of environments, to better clarify the observational framework. Results are displayed in Fig. 9.
![]() |
Figure 9:
Run of the [Mn/Fe] ratio as a function of [Fe/H] for
stars in different environments. Light crosses ([Fe/H]
|
| Open with DEXTER | |
To have a reference, we included field stars from both the accretion
and dissipative components as found by Gratton et al. (2003b), covering the
range of metallicity up to [Fe/H]
.
These data were complemented at
higher metallicity with thin disk stars from the study of Reddy et al. (2003),
after an offset of +0.12 dex was added to their [Mn/Fe] ratios to take into
account differences in the adopted solar reference abundances. Also, we plotted
data from Feltzing & Gustaffson (1998), after subtracting 0.04 dex to mantain
the offset of 0.16 dex with Reddy et al. data, as adopted in McW03.
Data points for red giants in Sgr were taken from Bonifacio et al. (2000) and were read from figures in the McW03 paper, like those for the galactic bulge.
Several features in this diagram are worth of discussion.
1) Field stars show the same pattern noted by Gratton (1989): decreasing
[Mn/Fe] for decreasing metallicity, from super-solar values down to about
[Fe/H]
,
where a plateau, although with large scatter, can be
envisioned around [Mn/Fe]
.
Both the quoted
interpretations of SN Ia producing a higher
fraction of Mn with respect to SN II and/or more metal-rich SN Ia yielding more Mn than metal-poor SN Ia can explain this result.
2) Most of the red giants surveyed to date in dwarf spheroidal galaxies
(dSphs, open circles in Fig. 9, from Shetrone et al. 2001,
2003) span a range in metallicity restricted between -2 < [Fe/H] <-1, and
the overwhelming majority is scattered around a mean value of [Mn/Fe]
,
with no trend increasing the observed [Mn/Fe] ratios over 1 full dex in [Fe/H]. Stars in present day dSphs seem to behave differently from those in the cannibalized dwarf Sgr (see
below). Our choice of the field reference sample does not affect our
conclusions; it would hold also using the sample assembled by McW03.
3) The locus occupied by stars in dSphs seems to be coincident with that populated by the accretion component selected by Gratton et al. (2003b), in the [Mn/Fe]-[Fe/H] diagram. This field population is composed of non-rotating or even counter-rotating halo stars, with distinctive chemistry with respect to the dissipative component of the halo/thick disk.
Hence, as far as Mn is concerned, it is tantalizing to conclude that these stars, likely relics of objects accreted or captured in the past by our Galaxy, share the same kind of chemical signature as present-day dSphs, i.e. stellar systems with star formation occurring at a uniformly low rate, even if at different times (Tolstoy et al. 2003).
On the other hand, it is well known that the pattern of
-elements differs between stars in dSphs and in the galactic field. The
low (even subsolar) ratios observed in giant stars of present-day dwarfs
suggest that the action of very massive stars was not strong enough to affect
the chemical evolution in these objects (Tolstoy et al. 2003).
The stars of the galactic accretion component (Gratton et al. 2003b), however,
show higher [
/Fe] ratios; it is possible that in the
primordial satellites that were likely accreted in the Galaxy, leaving behind
this accreted component, stars with initial masses >
were
able to contribute a higher content of
-elements.
4) The disrupting dwarf spheroidal Sgr seems to represent a different case. The five stars analyzed by McW03 and the two by Bonifacio et al. (2000) span a much more extended range in metallicity and reach a higher overall metallicity, maybe reflecting a higher mass of the host galaxy. Sgr was probably able to retain a higher amount of gas, reaching up to solar values of [Fe/H]. The most metal poor component of Sgr seems to be intermingled with field stars, and as well with the most metal-poor component of the bulge population.
Plotting in the same panel both bulge and Sgr stars analyzed by McW03, it seems that in the high metallicity regime these populations seem to differ from the galactic field population. However considering also the two stars analyzed by Bonifacio et al. (2000) this distinction seems less clear; Mn abundances for the larger sample presented by Bonifacio et al. (2003) will allow a clearer view of the situation. The conclusion of McW03, that type Ia SNe produce metallicity-dependent yields, as far as Mn is concerned, is considerably weakened if the stars of Bonifacio et al. (2000) are also considered.
5) Concerning the bulge population, the stars analyzed in McW03
seem to bracket the field population. Using the new data, however, and by
plotting side-by-side bulge and Sgr stars, it is quite
difficult to distinguish the different populations, at the metal poor end.
In our opinion, in the high
metallicity regime ([Fe/H]
), we can only conclude that data encompass
a large range of [Mn/Fe] ratios. This impression is
strengthened by the average
abundance of the bulge cluster NGC 6528 (Carretta et al. 2001), which is based
on four cluster stars. We therefore caution that
it may be dangerous to argue about
implications from Mn in metal-rich populations, since the scatter at solar
metallicity is as large as 0.6 dex in [Mn/Fe].
6) The run of globular clusters in the [Mn/Fe] vs. [Fe/H] plane seems to be more
or less flat, with some scatter. NGC 6528, as already noted, stands out, since
its [Mn/Fe] ratio is about 0.4 dex below the mean value of stars in the solar
neighborhood. This cluster is also rather peculiar for its kinematics: its high
radial velocity rules out the possibility that NGC 6528 is a disk cluster, but
we cannot exclude that we are just seeing an inner halo cluster presently
passing through the bulge (see discussion in Carretta et al. 2001).
47 Tuc, on the other hand, seems to lie under the mean locus defined by the
field stars at [Fe/H]
.
7) Finally, a recent paper by Ivans et al. (2003) shows that the
situation is
much more complex. They studied in detail three well known
-poor
stars (five-pointed stars in Fig. 9; we dropped in this
figure star BD+24
1676, whose Mn abundance is based on only 1 line).
While these stars have
low amounts of
-elements also shown by stars
in dSphs,
their Mn content stands out about 1 dex above the mean level of stars
of similar (very low) metallicity. Moreover, it is likely that these peculiar
stars formed in sub-systems that were distinct from the global galactic
population and subsequently accreted, as hinted by their orbital characteristics
(see Ivans et al. 2003). The interplay between the
abundances (yields) of
-elements and those of Mn is still poorly
understood. The position of these stars
argues that the yield of Mn and Fe are quite different in SN II and SN Ia, in
agreement with theoretical predictions (e.g. Tsujimoto et al. 1995).
8) The situation for Mn is somehow mirrored by that of Ni, an even element:
while in the Galactic field stars Ni tracks iron very closely ([Ni/Fe]
)
down to extremely low metallicities ([Fe/H]
,
Cayrel et al.
2004), the GCs Ruprecht 106, Pal 12 (Brown et al. 1997) and Terzan 7 (Sbordone
et al. 2004) as well as the Sgr field stars (Bonifacio et al. 2002) show
[Ni/Fe] < 0. Pal 12 and Terzan 7 are associated with Sgr. In all these
objects [
/Fe] is solar or sub-solar. 47 Tuc behaves for Ni like the field stars and contrasts quite strikingly with
Terzan 7, which has a similar metallicity. Any theory which seeks to explain
anomalous Mn abundances should at the same time be able to explain the Ni abundances.
In this framework, it seems rather difficult to use our derived abundances for 47 Tuc to put strong constraints on the production of Mn. The commonly adopted hypothesis are: (i) the contribution from SN Ia is larger than that from SN II; in this case, our results imply that type Ia SNe contributed less to the building up of elements in 47 Tuc with respect to field stars of similar metallicity; and (ii) Mn yields are metallicity-dependent (via the neutron excess) and the SNe polluting the gas that formed 47 Tuc were more metal-poor (by about 1 dex) than those enriching the field stars.
In both cases, however, we might interpret these data as hints of a peculiar chemical history for 47 Tuc. Two tentative scenarios can be envisioned:
(a) the proto-cluster cloud that formed 47 Tuc was able to build up a significant fraction of its own metals;
(b) 47 Tuc might come from a larger object having a chemical evolution quite distinct from that of our Galaxy. In the scenario of an accretion of a rather large satellite in the past of our Galaxy (see e.g. Gratton et al. 2000a; Freeman & Bland-Hawthorne 2002), "inflating" the primordial disk to produce the present thick disk, 47 Tuc could have formed out of the material of the infalling satellite during the accretion phase. To support this, we note that the age of this cluster is about 2 Gyr less than the mean age of typical old galactic clusters (Rosenberg et al. 1999; Gratton et al. 2003a). Moreover, among the metal-rich ([Fe/H] >-0.8) and more massive clusters, 47 Tuc is the object that lies furthest from the galactic plane (-3 kpc, from the February 2003 update of Harris 1996). Besides, there are very few stars of similar metallicity in the dissipative component whose maximum height above the galactic plane outrank the distance at which 47 Tuc is presently observed below this plane (Gratton et al. 2003b).
We caution that the observational framework for the Mn abundances is still quite complex and poorly understood. Is it possible that 47 Tuc, one of the most massive globular clusters in the Galaxy, has been formed in a formerly independent subsystem? While at present it is not possible to give a clearcut answer with the available data, we suggest that some hints do exist from both its chemical composition and orbital features.
Using the high resolution UVES spectra we firmly established the
metallicity of 47 Tuc to be [Fe/H] =-0.67
0.05, in substantial agreement
with previous results in the literature. As for the more metal-poor
clusters NGC 6397 and NGC 6752, we find that TO and SGB stars share
the same iron abundance and there is no star-to-star scatter in excess
of that expected from observational errors.
As in NGC 6752, we have been able to demonstrate the existence of a Na-O anticorrelation which must be the signature of a previous generation of stars, since TO and SGB stars do not have high enough internal temperatures to produce such patterns, no matter how deep the mixing.
Some scatter has been found also for Li, Mg and Al; no clear correlations have been found for the last two elements.
For most of the other elements, 47 Tuc seems to behave as field stars of similar metallicity, with two possible exceptions: Ti seems to be marginally higher than in field stars, Mn seems to be marginally lower. The errors in our analysis and the scatter among the field stars are such that none of these anomalies is highly significant; however, if real, they could point towards the origin of 47 Tuc as part of an independent, possibly larger, system later captured by the Galaxy.
Acknowledgements
This research has made use of the SIMBAD data base, operated at CDS, Strasbourg, France, and was funded by COFIN2001028897 by Ministero Università e Ricerca Scientifica, Italy. We wish to thank F. Grundahl for providing photometry of stars in 47 Tuc. We thank the ESO staff at Paranal (Chile) for their help during observing runs and Elena Pancino for reading the manuscripts and making useful comments.