A&A 417, 293-300 (2004)
DOI: 10.1051/0004-6361:20031759
S. Ambika 1 - M. Parthasarathy 1 - W. Aoki 2 - T. Fujii 2,3 - Y. Nakada 4,5 - Y. Ita 4 - H. Izumiura 6
1 - Indian Institute of Astrophysics, Koramangala, Bangalore, 560034, India
2 - National Astronomical Observatory, Mitaka, Tokyo 181-8588, Japan
3 - Faculty of Science, Kagoshima University, 1-21-35 Korimoto, Kagoshima
890-0065, Japan
4 - Institute of Astronomy, The University of Tokyo, Mitaka, Tokyo, 181-0015,
Japan
5 - Kiso Observatory, Institute of Astronomy, University of Tokyo, Mitake,
Kiso, Nagano 397-0101, Japan
6 - Okayama Astrophysical Observatory, National Astronomical Observatory,
Kamogata, Okayama 719-0232, Japan
Received 20 May 2003 / Accepted 18 November 2003
Abstract
We present a detailed model-atmosphere analysis of ZNG 4, a
UV-bright star
in the globular cluster M 13. From the analysis of a high resolution
(
)
spectrum of the object, we derive the
atmospheric parameters to be
K,
log
and
.
Except
for magnesium, chromium and strontium, all other even Z elements are
enhanced with titanium and calcium being
overabundant by a factor of 0.8 dex. Sodium is enhanced by a factor of 0.2
dex. The luminosity of ZNG 4 and its position in the color-magnitude
diagram of the cluster indicate that it is a Supra Horizontal Branch (SHB)
(post-HB) star.
The underabundance of He and overabundances of Ca, Ti, Sc and Ba in the
photosphere of ZNG 4 indicate that diffusion and radiative
levitation of elements may be in operation in M 13 post-HB stars even at
of 8500 K.
Detailed and more accurate abundance analysis of post-HB stars in several
globular clusters is needed to further understand their abundance anomalies.
Key words: stars: abundances - stars: evolution - stars: Population II - stars: horizontal-branch - Galaxy: globular cluster: individual: M 13 - stars: individual: ZNG 4
The term "UV-bright stars" was introduced by Zinn et al. (1972) for stars in globular clusters that lie above the horizontal branch (HB) and are bluer than red giants. The name resulted from the fact that, in the U band, these stars were brighter than all other cluster stars. Further investigations showed that this group of stars consist of blue horizontal branch (BHB) stars, supra horizontal branch stars (SHB), post asymptotic giant branch stars (post-AGB), post-early AGB (P-EAGB) stars and AGB-manque stars (de Boer 1985, 1987; Sweigart et al. 1974; Brocato et al. 1990; Dorman et al. 1993; Gonzalez & Wallerstein 1994).
To derive the chemical composition of UV-bright stars in globular
clusters
and to understand their evolutionary stages, we started a program
to obtain high resolution spectra of these objects in selected globular
clusters with the High Dispersion Spectrograph (HDS, Noguchi et al. 2002) of
the 8.2 m Subaru
Telescope. We selected a few UV-bright stars in the
globular cluster M 13 from the papers of Zinn et al. (1972) and Harris
et al. (1983) to derive their chemical composition. In this paper we
report the analysis of a high resolution spectrum of the UV-bright star
ZNG 4 (RA (
)
and Dec (
)
(2000)) (Zinn et al. 1972) in M 13 as the first target of our program.
M 13 (NGC 6205) is a nearby wellstudied globular cluster with a
distance modulus of
and metallicity
of
(Kraft & Ivans 2003).
The position of ZNG 4 in the color-magnitude diagram of M 13 (Paltrinieri
et al. 1998) is shown in Fig. 1. Many of the globular
clusters show a prominent gap in the blue tail of the HB, which is
presumed to be due to differential mass loss on the Red Giant Branch (RGB).
In M 13 it is observed at
K (Ferraro et al. 1997).
High resolution spectroscopic studies of M 13 BHB stars lying on either
side of the gap were carried out by Peterson et al. (1983, 1995) and
Behr et al. (1999, 2000a). They found anomalous photospheric
abundances in BHB stars. These photospheric anomalies are most likely
due to diffusion - the gravitational settling of helium and radiative
levitation of the metal atoms in the stable atmosphere of hot stars.
They found variations in the photospheric abundances and rotational
velocities of BHB stars as a function of their effective temperatures.
![]() |
Figure 1: Color magnitude diagram (CMD) of globular cluster M 13 obtained by Paltrinieri et al. (1998). The arrow indicates the position of ZNG 4 in the CMD. |
Open with DEXTER |
We have obtained a high resolution (
)
spectrum of ZNG 4 on 15th (UT: 14h 45m) April
2001 with the Subaru/HDS. The spectrum covering the
wavelength range
4142
-6814
was obtained in an exposure time of 20 min.
There was no moon light problem during the observations and the sky
background in the data was close to zero. We neglected the sky background
in our data reduction.
The data was bias-subtracted, trimmed, flat-fielded to remove pixel to pixel variations, converted to a one-dimensional spectrum, and normalized to the continuum using standard CCD data reduction package (NOAO IRAF). The spectrum has an average signal to noise ratio of 35. The reference spectrum of thorium-argon was used for the wavelength calibration.
The various orders in our echelle spectrum of ZNG 4 have well defined continuum and the normalization of the continuum was carried out using the IRAF echelle spectra reduction programs. The continuum level in the adjacent echelle orders to those containing the Balmer lines was useful in defining the continuum in the Balmer line regions and the profiles were normalized with a polynomial fit.
The radial velocity of ZNG 4 was derived from the wavelength shifts of many
absorption lines. The average heliocentric velocity is found to be
which
is in agreement with the value derived by Zinn (1974) (-253 km s-1).
It is also in agreement with the heliocentric velocities of M 13 BHB stars
derived by Behr et al. (1999) and Moehler et al. (2003).
For the initial estimate of effective temperature, we looked for the
published CCD photometry of the star.
Recent CCD photometry of M 13 was carried out by Rey et al. (2001).
However the ZNG 4 area of the cluster was not included in their observations
(Rey, private communication). We used the published CCD photometry
of ZNG 4 by Paltrinieri et al. (1998), who give, B=14.096 and V=13.964.
(B-V) = 0.132 and
E(B-V) = 0.02 (Kraft & Ivans 2003) will yield
(B-V)0 = 0.112 which corresponds to
K
(Flower 1996). However, the (B-V)0 and
calibration
given by Flower (1996) is for Population I stars.
For our analysis, excitation potential and oscillator strengths of the lines
were taken from the Vienna Atomic Line Database
(http://www.astro.univie.ac.at/~vald/). We employed the latest (2002) version of MOOG, an LTE stellar line analysis
program (Sneden 1973) and Kurucz (1993) grid of ATLAS models.
MOOG has been used successfully in the analysis of the spectra of warmer
stars with
K (Preston & Sneden 2000).
We have also analyzed the spectra using the Kurucz WIDTH program (Kurucz CDROM 13, 1993) for verification. We used the line list obtained using version 43 of the Synspec code of Hubeny and Lanz which is distributed as part of their TLUSTY model atmosphere program. (http://tlusty.gsfc.nasa.gov/Synspec43/synspec-line.html) and also the information from the Kurucz linelist (http://kurucz.harvard.edu/linelists.html).
The value of effective temperature was obtained by the method of excitation balance, forcing the slope of abundances from Fe I lines versus excitation potential to be zero. The surface gravity was then set by ionization- equilibrium, forcing abundances obtained from neutral (Fe I) and ionized (Fe II) species to be equal. The microturbulent velocity was estimated by demanding that there should be no dependence of the Fe I abundance upon equivalent widths of Fe I lines.
![]() |
Figure 2: Top figure is the plot of the abundances from Fe lines versus excitation potential of the lines. Figure at the bottom is the plot of the abundances from Fe lines versus their equivalent widths. Triangles represent the Fe I lines and stars denote the Fe II lines. |
Open with DEXTER |
The plots of the abundances versus excitation potentials and abundances
versus equivalent widths in the case of Fe I and Fe II lines are shown
in Fig. 2. Such plots were made by varying the
,
log gand
in steps of 250 K, 0.5 and 0.5 km
respectively
to estimate the uncertainties in these parameters.
From our analysis, we find that
K, log g = 2.5 and
fit the data best (Fig. 2). From the
abovementioned method of analysis we find the uncertainties in
to be 250 K, log g = 0.5 dex and
km
.
Uncertainties in derived abundances as a result of errors in the
determination of the parameters and errors in the measurements
of equivalent widths are found to be of the order of 0.2 dex.
Table 1: Chemical composition of ZNG 4 in M 13.
Using the derived atmospheric parameters and abundances, a synthetic spectrum was generated and plotted over the observed spectrum for verification. The observed and synthetic spectra were found to match well with the above mentioned atmospheric parameters and the final abundances are given in Tables 1 and 3. A region of the observed and synthetic spectrum is shown in Fig. 3. The abundances derived using MOOG (Tables 1 and 3) are in good agreement with the abundances derived using the WIDTH (Table 1).
![]() |
Figure 3:
Synthetic spectrum calculated with the atmospheric parameters
(
![]() ![]() ![]() ![]() |
Open with DEXTER |
We tried to estimate the
and log g from
the analysis of Balmer lines in the spectrum of ZNG 4 using
the Kurucz spectral atlas for Balmer lines (Kurucz CDROM 13, 1993).
We could not get a satisfactory fit between the observed and theoretical
Balmer line profiles with the atmospheric parameters
K,
log g= 2.5,
km
and
.
We also tried models that take into account the alpha element enhancements
and no convective overshooting
(ANOVER models: http://kurucz.harvard.edu/grids/gridm15ANOVER/bm15ak2nover.dat).
They did not make any difference to the abovementioned atmospheric parameters
(Fig. 4).
The best fit to
profile was obtained with the parameters
K, log g= 2.0,
km
and
(Fig. 4).
However, the excitation balance and ionization equilibrium for Mg and Fe lines
could not be achieved with the above
parameters (see the last column in Table 1) and the abundances of Mg I, Mg II and Fe I, Fe II were found to differ
significantly (Table 1).
Therefore, we chose the model atmosphere determined from the
analysis of metal lines (i.e.
K, log g= 2.5,
km
and
)
to represent the atmosphere of the star.
The problem of fitting Balmer line profiles of HB stars has been mentioned by Grundahl et al. (1999) (and references therein). For stars being more luminous than HB stars, mass loss and/or extended atmosphere may influence the Balmer line profiles (Vink & Cassisi 2002).
![]() |
Figure 4:
Observed
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The mean abundances of ZNG 4 relative to the Sun (Anders & Grevesse 1989) are given in Table 1, together with the number of lines used in the analysis and the standard deviation of abundances estimated from individual species.
Table 2: Comparison of the abundances of ZNG 4 with the abundances of M 13 BHB and RGB stars.
![]() |
Figure 5: The abundances of different elements in ZNG 4 (represented by asterix) compared with M 13/ J11 (BHB star; triangles), M 13/WF4-3485 (hot BHB star; diamonds) and that of M 13/ L262 (RGB star; squares) with their error bars (Table 2). In the cases where only upper limit of the abundances are given (3rd and 4th column of Table 2), those elements are not shown in the above figure. Average abundances of Y and Ba in M 13 RGB stars are from Armosky et al. (1994). |
Open with DEXTER |
Analysis of Mg lines gives
which is the same as the M 13
cluster metallicity. (The equivalent widths of the Mg II lines at 4481
(Table 3)
were obtained by multiple gaussian fit to the lines in the observed spectrum).
Silicon is overabundant compared to iron
(
). Calcium and titanium are found to be overabundant
(
and
). There is a 0.5 dex
difference in the abundances derived from the Ca I lines at 4226.73
and
4454.78
(Table 3). However, the abundance of Ca derived from the Ca II line
at 5019.97
is in agreement with that derived from the Ca I 4454.78
line.
The reason for the deviation in the abundance derived from the Ca I 4226.73
line is not clear. It may be due to the relatively low signal to noise ratio
of the data around this wavelength range.
There seems to be no interstellar contribution to the Ca II line.
Since the star has a radial velocity of -257 km s-1, the stellar
lines are well
separated from the lines of interstellar origin.
The abundance of Cr and Fe (
and
)
are found to be close to the metallicity of the cluster.
On the other hand, Sc is found to be overabundant (
).
Na lines show an overabundance of 0.2 dex.
We have detected one line of Sr II, two lines of Y II and two lines of Ba II.
Sr seems to be underabundant (
), while Y and Ba are
overabundant (
and
).
We have detected the He I line at 4471.47 ,
which yields an abundance
of
which implies an underabundance of 0.55 dex compared to the
solar value. This is in agreement with the underabundane of He found in hot
BHB stars (Moehler 1999; Moehler et al. 2003).
We have not detected C, N and O lines in our spectrum of ZNG 4.
Assuming an equivalent width of 5 m
as the detectable limit in our
spectrum of ZNG 4, we find the upper limit of
to be +0.32 dex
(based on the C I 5052.17
line), that of
to be +1.15 dex
(based on the N I 4214.80
line) and that of
to be
+0.01 dex (based on the O I 6156.78
line).
Globular cluster stars show anticorrelation
of sodium and oxygen abundances (Kraft et al. 1997). In ZNG 4 we find
enhancement of sodium, therefore we expect an underabundance of oxygen.
Also star to star abundance variations in the
light elements C, N, O, Na, Mg and Al occur among the bright giants of a
number of globular clusters (Ivans et al. 1999). The absence of C, N and O
lines in our spectrum of ZNG 4 may be due to the underabundance of these
elements. A much higher resolution and high signal to noise ratio spectrum of
ZNG 4 may reveal the lines of C, N and O lines if present.
The chemical composition of ZNG 4 shows significant deviations in element abundances from the expected metallicity of M 13. Ti, Ca, Sc and Ba are found to be relatively overabundant (Fig. 5) compared to RGB stars and cool BHB stars of M 13, whereas the abundances of Mg, Cr and Fe are in agreement with the cluster metallicity.
From a study of 22 M 13 G-K giants, Kraft et al. (1993, 1997) found the
abundances of Fe, Sc, V and Ni to be
,
,
[V/Fe] = 0.00 and
.
They found Ca and Ti
to be mildly overabundant ([Ca/Fe] = +0.24 and [Ti/Fe] = +0.29). Si was
overabundant by +0.34 dex. Study of M 13 giants by Armosky
et al. (1994) yields the average abundance of Fe, Y and Ba:
,
and
,
whereas in ZNG 4 we find significant
overabundance of Ca, Sc, Ti, Y and Ba compared to that found in
M 13 giants. Also, overabundance of Na [+0.2 dex] and absence of O lines
support the anticorrelation of Na and O abundances found in M 13 giants
(Kraft et al. 1997).
Behr et al. (1999) have studied the BHB stars in M 13 on either side of the
HB gap.
They find the photospheric compositions and stellar rotation rates to
vary strongly as a function of
of the stars.
Among the cooler stars in their sample, at
of 8500 K, the metal
abundances are in rough agreement with the canonical cluster metallicity
and the hotter stars with
greater than 10 000 K show a deficiency
of He and enhancement of Fe, Ti, Cr by a factor of 300. However,
Mg remains at the canonical cluster metallicity. In ZNG 4 also, Mg abundance is
in agreement with the M 13 metallicity. Abundances similar to that found in
BHB stars of M 13 were also found in BHB stars
of M 15 (Behr et al. 2000b) and BHB stars of NGC 6752 (Glaspy et al. 1989;
Moehler et al. 1999). The abundance anomalies in these
BHB stars are most likely due to diffusion - the gravitational settling
of helium (Greestein et al. 1967) and radiative levitation of metal
atoms (Michaud et al. 1983). Rotational velocities (
)
appear to
have a bimodal distribution in cooler BHB stars, whereas the hotter
BHB stars with
greater than 10 000 K are found to be slow
rotators.
However, the abundances and rotational velocity of ZNG 4 are contrary to what
is expected for stars on the red side of the HB gap at 11 000 K. It is a slow
rotator (
km
). It shows underabundance of He and
overabundances of Ti, Ca, Sc and Ba.
In Table 2 and Fig. 5, we have compared the
abundances of ZNG 4 with the abundances of M 13 BHB stars of
7681 K (J 11)
and
12 750 K (WF4-3485) (Behr 2000c) and with the
abundances of the RGB star L262 (Cavallo & Nagar 2000). In Table 2, the last
column shows the mean RGB abundances in M 13, where elements from Na to Fe
are taken from Kraft et al. (1997) and Y and Ba abundances are from
Armosky et al. (1994). It is evident from the abundances listed in Table 2
and from Fig. 5 that ZNG 4 shows overabundance of metals
compared to that of a M 13 RGB star and also when compared to that of a M 13 BHB
star of similar temperature.
These results indicate that in ZNG 4, diffusion and radiative levitation
of elements may be in operation. Slowly
rotating HB stars are also seen on the cooler side of HB gap, but abundance
anomalies start from 11 000 K (Moehler et al. 1999; Behr et al. 2000b).
This implies that ZNG 4 may have the properties of the stars on
the blue side of the HB gap although it has
of 8500 K.
In this regard, more accurate determination of abundances of these elements
in ZNG 4 and similar stars in M 13 is needed to confirm our results and
conclusions.
This may be explained in two ways. One is that, for some stars in M 13,
the onset of diffusion seems to start at lower
(
8500 K). The other argument would be that the star has evolved from the blue
side of the HB
gap and is moving towards the red with higher luminosity as indicated by the
post-HB evolutionary tracks of Gingold (1976) and Dorman et al. (1993).
The BHB stars hotter than 11 500 K typically show strong photospheric helium depletions due to gravitational settling (Moehler et al. 2000, 2003). The calculations of Michaud et al. (1983) indicate that helium depletion should be accompanied by photospheric enhancement of metals, since the same stable atmosphere that permits gravitational settling also permits the levitation of elements with large radiative cross sections. The depletion of helium and overabudance of some of the metals in the photosphere of ZNG 4 is in qualitative agreement with the calculations of Michaud et al. (1983).
Recently, Turcotte et al. (1998) and Richer et al. (2000) made diffusion
simulations to explain the abundance patterns of chemically peculiar A and
F stars. Their predicted abundance patterns are qualitatively similar to
that found in ZNG 4. However, none of the recent diffusion studies
treated the cases of BHB stars and post-HB stars. This phenomenon
may be related to the disappearance of surface convection and
hence to the formation of a stable stellar atmospheres. HB stars
and post-HB stars cooler than
K have deep convective
envelopes (Sweigart 2002). Hotter than this temperature the envelope
convection breaks into distinct shells associated with the
ionization of H and He. Note that the surface convection
disappears at 11 000 K (Sweigart 2002) and BHB stars hotter than
this show moderate to severe abundance anomalies.
ZNG 4 has a V magnitude of 13.964 (Paltrinieri et al. 1998). Considering
the distance modulus of M 13 to be 14.42 and E(B-V) towards M 13
to be 0.02 (Kraft & Ivans 2003), we estimated the absolute magnitude
(Mv) of the star to be -0.522. For stars with
around 8500 K, the bolometric correction (BC) is negligible (Flower 1996).
Considering BC=0, we get the bolometric magnitude (
)
to
be -0.522, which corresponds to a luminosity of 127
[
].
Using the equation connecting the mass, effective temperature and bolometric
magnitude, we find the surface gravity, log g = 2.6 (assuming the mass of ZNG 4
to be 0.5
), which agrees well with the value estimated from
the analysis of the spectrum of ZNG 4.
The post-AGB star Barnard 29 [
]
which is a member of M 13 is more luminous than ZNG 4.
The abundance pattern of ZNG 4 is very different from
that of the post-AGB star Barnard 29 (Conlon et al. 1994; Moehler et al.
1998). The M 13 BHB stars with
around 8500 K have a luminosity
of about 40
,
whereas ZNG 4 is more luminous by about a
factor of 3, which indicates that ZNG 4 can be classified as a supra
horizontal branch star (post-HB).
Stars that lie 1.5 magnitude above the HB stars are classified
as supra horizontal branch (SHB) stars in the photometric studies of M 13
(Zinn 1974) and NGC 6522 (Shara et al. 1998).
No detailed abundance analysis of SHB stars in globular clusters
is available to compare with the abundances of ZNG 4.
Since ZNG 4 is a post-HB star and it has evolved from a hot BHB star stage
and may had severe abundance anomalies similar to those found in
the hot BHB stars of M 13 (Table 2). The present
K of
ZNG 4 indicates that thin layers of subsurface convection if present may
have diluted the severe abundance anomalies due to diffusion and radiative
levitation that took place during its hot BHB stage of evolution. It is
important to derive the chemical composition of a significant sample of
post-HB stars hotter than 11 000 K and much cooler than 11 000 K
to further understand the role of diffusion, radiative
levitation, rotation and convection during the post-HB stage of evolution.
Acknowledgements
We would like to thank Dr. B.Dorman and Dr. B. Behr for helpful information on M 13 BHB stars, Dr. F. R. Ferraro for kindly providing the table of B and V magnitudes of M 13 cluster stars and the referee Dr. S. Moehler for helpful comments.
Table 3: Data for spectral lines measured in the spectrum of ZNG 4 in M 13.