A&A 449, 127-134 (2006)
DOI: 10.1051/0004-6361:20053770
H. W. Zhang1 - G. Zhao2
1 - Department of Astronomy, School of Physics, Peking University,
Beijing 100871, PR China
2 -
National Astronomical Observatories, Chinese Academy of
Sciences, Beijing 100012, PR China
Received 22 July 2005 / Accepted 26 October 2005
Abstract
Context. The formation scenario of the Galactic thick disk is an unresolved problem. Chemical abundances in long-lived dwarf stars of the thin and thick disks provide information of the Galactic disk formation.
Aims. We present photospheric abundances of the O, Na, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Ni, and Ba elements for 32 mildly metal-poor stars with
.
According to their kinematics, age, and [
/Fe], sample stars are identified to thin disk, thick disk, and halo population memberships. Element abundances for sample stars are discussed as a function of metallicity.
Methods. High resolution and high signal-to-noise ratio spectra were obtained with the Coudé Echelle Spectrograph mounted on the 2.16 m telescope at the National Astronomical Observatories (Xinglong, China). Effective temperatures were estimated from colour indices, and surface gravities from Hipparcos parallaxes. Stellar abundances were determined from a differential LTE analysis. The kinematics parameters were calculated from the parallax, proper motion, and radial velocity. Stellar ages were determined from theoretical stellar evolution tracks.
Results. The average age of the thick disk stars is older than the thin disk stars. Our element abundance results extend and confirm previous works. The oxygen and other -elements (Mg, Si, Ca, and Ti) abundances of thin and thick disk stars show distinct trends at
.
The [Al/Fe] behaviour is exactly as an
-element, although the separation for [Na/Fe] of thin and thick disk stars is not clear. The elements V, Cr, and Ni follow Fe very closely, and there is no offset between thin and thick disk stars, but the Sc and Mn abundance trends of the thin and thick disk stars are different, and [Ba/Fe] of thin disk and thick disk stars shows different behaviour.
Key words: stars: abundances - stars: kinematics - Galaxy: disk
The existence of a population of stars with kinematics, ages, and chemical abundances that lie between the characteristic values for the halo and the disk populations is a long-standing problem in studies of Galactic structure and evolution.
Observational works from the 1980's showed that Galaxy as having two disk-like components: the thin and the thick disks (Gilmore & Reid 1983). The formation scenario of the Galactic thick disk is still unresolved (Majewski 1993). There are essentially two major formation scenarios for the Galactic thick disk: (i) the pre-thin disk (top-down) models, where the formation of the thick disk precedes the formation of the thin disk; and (ii) the post-thin disk (bottom-up) models, where the thick disk is the result of some action on or by the thin disk. Much more work on stellar ages, kinematics, and abundances has to be carried out before we can be certain about the basic scenario for the formation and evolution of our Galaxy.
The study of chemical abundance in long-lived dwarf stars of the thin and thick disks provides information about the Galactic disk formation and has become an active research area in recent years. A strong indication of the Galactic thin and thick disks as discrete populations with respect to kinematics and age comes from Edvardsson et al. (1993). The following works show systematic differences between the chemical composition of the thin and thick disk stars: Fuhrmann (1998, 2004), Gratton et al. (2000), Prochaska et al. (2000), Mashonkina & Gehren (2001), Reddy et al. (2003), Gehren et al. (2004, 2005), Brewer & Carney (2005). Recently, Bensby et al. (2003, 2004, 2005) selected two samples with a high kinematical probability of belonging to either the thin and thick disk and found distinct abundance trends at sub-solar metallicities. For a comprehensive review, the reader is referred to Nissen (2003).
Chen et al. (2000) studied the chemical composition of 90 F and G dwarf stars, and do not find any clear [/Fe] separation
between thin and thick stars. As pointed out by Prochaska et al.
(2000), this may, however, be due to the fact that they selected
dwarf stars in the temperature range
K. Hence, the old thick disk stars with
K were not included.
In order to extend the work of Chen et al. (2000), we study a sample
of candidates of thick disk stars. In this paper the results of 32 mildly metal-poor stars with
are presented. In
Sect. 2 we describe the observation and the data reduction, and in
Sect. 3 the abundance determination. In Sect. 4 we determine the
kinematic parameter, age, and population membership of sample stars.
In Sect. 5 we discuss the abundance trends with metallicity,
conclusions followed by in the last section.
The observation was performed on 10 nights during 1999 using the
Coudé Echelle Spectrograph attached to the 2.16 m telescope at
the National Astronomical Observatories (Xinglong, China). The red
arm of the spectrograph with a 31.6 grooves/mm grating was used in
combination with a prism as cross-disperser, thus providing a good
separation for the different echelle orders. With a 0.5 mm slit
(1.06
), the resolving power was 37 000 in the
middle-focus camera system. The detector was a Tek CCD (
pixels with
m2 each in size).
The signal-to-noise ratio of spectra at 6400 Å was about 150 pixel-1. The wavelength coverage was from 5600 Å to 8300 Å with some gaps. A detailed description of the technical aspects
of the spectrograph can be found in Zhao & Li (2001). During the
observation run, some spectra of bright, rapidly rotating,
spectral-type B stars were observed. These spectra were used to
remove the telluric O2 features.
A reduction of two-dimensional echelle spectral data was performed
using the ESO-MIDAS software package. The data reduction includes locating
the echelle order on the multi-order two-dimensional spectrum,
subtracting the background, and extracting the orders by summation
along the slit. The pixel-to-pixel variation was corrected by
dividing with a flatfield taken at the same night. The wavelength
calibration was based on a thorium-argon lamp. The radial velocities
were measured from about 20 intermediately strong and unblended
lines with an accuracy of 1.5 km s-1. The spectra were then
normalized by a continuum function determined by fitting a spline
curve to a set of pre-selected continuum windows (typically 20-30 per order) taken from the solar atlas.
The equivalent widths (EWs) were measured by direct integration and Gaussian fitting, depending on which method gave the best fit of the line profile. The EWs of all lines are given in Table 1. The accuracy of the EWs was estimated by comparing them to the independent measurements by Chen et al. (2000) for the 9 stars we have in common, which is shown in Fig. 1. The systematic difference between the two sets of measurements is small, and a linear least squares fitting gives:
![]() |
(1) |
The standard deviations around the relation are 5.0 mÅ (for 676 lines in common with Chen et al. 2000). We estimate an rms error of about 3 mÅ in our EWs.
![]() |
Figure 1: A comparison of equivalent widths measured in this work with those in Chen et al. (2000) for 9 stars in common. |
Open with DEXTER |
Table 2: Stellar parameters of sample stars.
The effective temperature was determined from the b-y and V-Kcolour indices using the IRFM calibrations of Alonso et al. (1996).
According to the discussion in Nissen et al. (2004), the error in
the effective temperature is expected to be around 70 K. With the
effective temperature and absolute magnitude, the stellar mass was
then determined from the star's position in the
diagram by interpolating the evolutionary tracks of
VandenBerg et al. (2000). The surface gravity was calculated from
Hipparcos parallax (ESA 1997) using the method described in Chen et al. (2000). We estimated that the error of log g is about 0.10 dex.
The microturbulence (
)
was obtained by requiring a zero
slope of [Fe/H] vs. EW.
The whole procedure of deriving
,
log g, and
was iterated to consistency when spectroscopic [Fe/H] was obtained.
The stellar parameters of sample stars are given in Table 2.
The abundance analysis was based on a net of flux constant, homogeneous, LTE model atmospheres interpolated from the extensive grids in Kurucz (1993). With the code ABONTEST8, which was kindly provided by Dr. Pierre Magain, abundances were derived by requiring that the calculated equivalent width from the model should match the observed one.
The details for line selection, atomic line data, calculation of the uncertainties, and the contribution from systematic errors are the same as those discussed in Chen et al. (2000). Solar abundances, calculated from the daylight spectrum, were used to derive stellar abundances relative to solar values. As expected, errors are around 0.05 dex for relative abundance, [X/Fe], and of the order of 0.10 dex for [Fe/H]. Finally, the abundances of the Fe, O, Mg, Si, Ca, Ti, Na, Al, Sc, V, Cr, Mn, Ni, Ba elements of 32 sample stars are given in Table 3.
We compared our results with recent studies that have some stars in common with ours. We found 9 stars in our sample are common to Chen et al. (2000), 5 stars to Bensby et al. (2003, 2005), and 4 stars to Reddy et al. (2003). The results of these comparisons are displayed in Table 4. We note that there is close agreement over the adopted atmospheric parameters and derived abundances.
Table 4: Average differences in atmospheric parameters and abundances.
Table 5: Kinematic properties, age, and population membership of sample stars.
The space velocities (U, V, W) and orbital parameters (
,
,
,
e) are calculated from the Hipparcos
parallax, proper motion, and radial velocity that we measured.
Calculation of the space velocity with respect to the Sun is based on the method presented by Johnson & Soderblom (1987). The correction of space velocity to the local standard of rest (LSR) is ( -10.0, +5.2, +7.2) km s-1 in (U, V, W), as derived by Dehnen & Binney (1998). The rotational velocity of the LSR with respect to the Galaxy is set to 225.2 km s-1 and the galactocentric distance is set to 8.5 kpc. The numerical integration of the Galactic orbits is done by adopting the Galactic mass model of Allen & Santillan (1991). A detailed discussion can be found in Zhang & Zhao (2005).
Using stellar effective temperature together with the absolute
magnitude based on the Hipparcos parallax approximate stellar age
can be interpolated according to [Fe/H] and [/Fe] from
theoretical stellar evolution tracks.
In this work we adopt the Yonsei-Yale isochrones (Yi et al. 2003),
which are calculated with new OPAL opacities and Kurucz model
atmospheres for a set of metallicities
Z = 0.00001, 0.0001, 0.0004,
0.001, 0.004, 0.007, 0.01, 0.02, 0.04, 0.06, 0.08, and [/Fe] = 0.0, 0.3, 0.6. The full set of stellar models and a FORTRAN
package that works for metallicity and
-enhancement
interpolation are available from the authors. We used different
-enhancement for sample stars according to our
determinations.
The identification of the population membership of a individual star is difficult. There are essentially two ways: the pure kinematical approach (e.g. Bensby et al. 2003, 2005) or by looking at a combination of kinematics, abundances, and stellar ages (e.g. Fuhrmann 1998, 2004).
In this work, we first calculated the relative probabilities for
the thick-disk-to-thin-disk (TD/D) and thick-disk-to-halo (TD/H)
membership using a method described by Bensby et al. (2003, 2005).
The population definition was mainly based on the probabilities,
but there are two exceptions: HD 108076 and HD 215257 have
intermediate membership probabilities. HD 108076 was re-defined to
thick disk star according to its age and -element
abundance, and HD 215257 to thin disk star according to its youth
and
-element abundance.
Famaey et al. (2005) identified several kinematical subgroups in a local sample of giants, representative of the solar neighbourhood. Among them, the Hercules stream is crucial because its motion can be confused with that of the thick disk. We have computed each star's probability as belonging to the Hercules stream and find this probability to be low, i.e. all stars are disk or halo stars.
Finally, 14 stars were defined as thin disk stars, 14 as thick disk
stars and 2 as halo stars. Of the remaining stars, HD 200580 recently
was identified as a visual binary with an angular distance of the
two stellar components of only 0.13 arcsec (Mason et al. 2001), and
HD 106516 is a single-lined spectroscopic binary and a blue straggler
candidate (Fuhrmann 2004). Table 4 gives the kinematic properties,
age, -elements abundance, and population membership for the
sample stars. In Fig. 2 the kinematic status of our sample stars was
shown in the Toomre diagram, while all our sample stars are
consistent with the Fuhrmann (1998, 2004) definition. Figure 3 shows
that thick disk stars have smaller
,
the minimum distance
in stellar orbit, than thin disk stars, so it is another feature to
use in identifying thin and thick disk stars.
Stellar age is a key feature for identifying stellar populations. For
all our thick disk stars, stellar age is larger than 10.0 Gyr and
the average is
Gyr. For comparison, the average age of
thin disk stars is
Gyr. However, the average ages of the
thin and thick disk stars given here are the results of our mildly
metal-poor samples, but the metallicity distributions of the thin
and thick disks peak at different metallicities, so our sample
probably includes many old thin disk stars.
![]() |
Figure 2:
Toomre diagram for the sample stars. Thin disk stars are
denoted by open circles, thick disk stars are filled circles, halo
stars are asterisk. The two circles delineate constant total space
velocities with respect to the LSR of
![]() |
Open with DEXTER |
![]() |
Figure 3:
![]() ![]() |
Open with DEXTER |
The element-to-iron ratio of each element is plotted as a function of [Fe/H] in Fig. 4. The discussions are given in the following subsections.
Both the quality of the spectra and the abundance analysis method of
Chen et al. (2000) are similar to our work and the results of two
works agree closely (see Table 4). In Fig. 5 we plot
-element abundances [
/Fe], which is defined as the
average abundance of Mg, Si, Ca, and Ti with respect to the Fe of
the two works. The population definition of the sample stars of Chen
et al. (2000) was made by using same method as in this paper.
As seen in Figs. 4 and 5, the -elements abundances of thin
and thick disk stars show distinct abundance trends at
.
The value of [
/Fe] in the thick disk
stars remains constant at a level of about +0.3 dex, whereas
[
/Fe] in the thin disk stars shows an increasing trend with
decreasing [Fe/H] and reaches thick disk values at
.
However, we noted a small overlap in [Fe/H]
between our thin and thick disk sample stars, so the abundance
trends must be checked in thick disk stars with higher [Fe/H] and
thin disk stars with lower [Fe/H].
For our two halo stars, [/Fe] of HD 194598 is +0.22
(
), lower than thick disk stars, and HD 219617 is +0.36
(
), similar to thick disk stars and to the
characteristic halo values.
The result is consistent with previous works. Fuhrmann (1998, 2004)
found the thin and thick disks to be chemically distinct for the Mg element. Bensby et al. (2003, 2005) showed that [/Fe] in the
thick disk stars is constant at
,
Above
,
[
/Fe] declines and the two disks merge together.
Soubiran & Girard (2005) and Brewer & Carney (2005) also found
similar behaviour.
This behaviour is interpreted as being due to the time delay between the type II SNe and long-lived type Ia SNe in the enrichment of the interstellar medium. Hence, star formation in the thick disk went on long enough that type Ia SNe started to enrich the gas out of which subsequent generations of thick disk stars formed.
It should be kept in mind that the population definition in most
works has only been based on stellar kinematics; whether the thick
disk stars with
are true thick disk stars or
contamination of thin disk stars needs further study. As pointed out
by Gehren et al. (2004, 2005), the combination of stellar
kinematics, age, and abundance makes the individual discrimination
between stars of the different populations possible.
Oxygen is also an -element that only forms in the interiors
of massive stars through hydrostatic burning of mainly He, C, and
Ne. Because oxygen is, next to hydrogen and helium, the most
abundant element in the Universe, its abundance has extra importance
for models of Galactic evolution. We discuss oxygen abundance in a
separate section.
The determination of oxygen abundance is often troublesome. The two sets of features most often used for oxygen abundance determinations are two forbidden [O I] lines at 6300 and 6363 Å and the permitted high-excitation O I triplet in the 7771-7775 Å region. The formation of the [O I] lines is described well by LTE calculations, but the lines are weak in the spectra of metal-poor solar-type dwarfs. On the other hand, O I triplet lines are strong in the metal-poor stars; however, the lines are known to be very sensitive to effective temperature and NLTE effects (see e.g. Kiselman 1991). In recent years there has been a lot of work on O I NLTE analysis (e.g. Gratton et al. 1999; Asplund 2004; Takeda & Honda 2005).
Our oxygen abundances are based on the O I triplet lines, which give systematically higher abundances than do forbidden lines. We use Eqs. (1) and (2) of Bensby et al. (2004) to scale the oxygen abundances derived from the triplet to those derived from the forbidden [O I] line.
As seen in Fig. 4, the oxygen abundances of thin and thick disk
stars also show distinct abundance trends, which agrees well with
what we found for other -elements (Mg, Si, Ca, and Ti).
![]() |
Figure 4: Abundance patterns. The symbols are the same as in Fig. 2. |
Open with DEXTER |
![]() |
Figure 5:
The ![]() |
Open with DEXTER |
Our result shows that [O/Fe] for thick disk stars with
is nearly constant, although the dispersions are
large. Bensby et al. (2004) present oxygen abundances for a large
sample of stars. They find distinctly different oxygen trends in the
thin and thick disks; at
the [O/Fe] trend is nearly
flat for thick disk stars, but the results of an extended sample
(Bensby et al. 2005) show an increasing trend with decreasing
metallicity. Jonsell et al. (2005) and Soubiran & Girard (2005)
also find that [O/Fe] possibly increases beyond
.
The
behaviour is different from other
-elements (Mg, Si, Ca, and
Ti). Bensby et al. (2004, 2005) also find different trends for
oxygen compared to other
-elements at
.
The
different behaviour between oxygen and other
-elements may
reflect the nucleosynthesis process of these elements, Bensby et al.
(2004) claim that oxygen is only produced in type II SNe. However,
this conclusion remains tentative, because the NLTE and other
effects make oxygen abundance determination with larger
uncertainties as compared with other
-elements.
Both sodium and aluminium can be subject to NLTE effects. Our sodium
abundances were determined from weak 6154/6160 lines. The NLTE
calculation by Shi et al. (2004) shows 6154/6160 lines giving small
NLTE effects (-0.1 dex) for the disk dwarfs. The [Al/Fe] is
slightly more affected by deviations from LTE, and the NLTE
correction is
+0.2 dex for our sample stars (Gehren et al.
2004).
Both Na and Al are mainly produced during the C and N burning in
massive stars, and they would therefore be ejected by SNe II, from a
phenomenological point of view, so that Al and Na could be
classified as mild -elements. As shown in Fig. 4, [Al/Fe]
behaves exactly like an
-element.
Our LTE sodium abundance trend is consistent with the NLTE result of Shi et al. (2004), although our [Na/Fe] systematically higher by 0.1 dex. If NLTE corrections are adopted, the differences will cancel.
Separation for [Na/Fe] of thin and thick disk stars is not quite as
clear as seen for [/Fe] and [Al/Fe]. The "merged''
appearance is also found by Bensby et al. (2003), but Bensby et al.
(2005) find that Na of the thin disk stars seems to be more abundant
than the thick disk stars. Brewer & Carney (2005) find that the Na
abundance remains near solar over the entire metallicity range for
both the thin and thick disk stars, although the average thick disk
abundance is slightly higher (0.04 dex) than that of the thin disk
stars.
The NLTE study of Gehren et al. (2004, 2005) shows gradual decline of [Na/Mg] with [Mg/Fe] for thin and thick disk stars, but no such trend for [Al/Mg] with [Mg/Fe]. The comparison of [Na/Mg] and [Al/Mg] will help us to understand the nucleosynthesis process of these two elements, but it is still not clear, while the yields of both Na and Al in SNe II are under debate.
As seen from Fig. 4, among the iron-group elements, V, Cr, and Ni follow Fe very closely, and there is no offset between thin and thick disk stars.
Hyperfine structure (HFS) correction of Sc and Mn lines are
significant, and our analysis adopts the hyperfine structure data of
Prochaska & McWilliam (2000). Previous studies of Sc show different
results. Zhao & Magain (1990) suggested that Sc was enhanced by
+0.25 dex in metal-poor stars. However, the result of Gratton
& Sneden (1991) showed a solar ratio for [Sc/Fe]. Nissen et al.
(2000) supports the high values of [Sc/Fe]. As seen from Fig. 4,
[Sc/Fe] of thin disk stars shows a decreasing trend with increasing
[Fe/H], thick disk stars with
have nearly constant
value,
,
but no significant enhancement at
.
Our [Sc/Fe] trend above
shows that Sc seems to follow the even-Z
elements like Mg, Si, Ca, and Ti; however, this similarity may not hold for stars with
,
which was also found by Prochaska & McWilliam (2000).
As seen from Fig. 4, there is a steplike change in [Mn/Fe] at
.
Thick disk stars with
have
,
whereas thin disk stars with
have
.
Nissen et al.
(2000) published a detailed study of the trend of [Mn/Fe] in disk
and metal-rich halo stars, but applied outdated data for the
hyperfine structure of the Mn lines. Using modern hyperfine
structure data, Prochaska & McWilliam (2000) found significant
corrections to the [Mn/Fe] values of Nissen et al. (2000). Their
revised data were plotted in Fig. 1.7 of Nissen (2003), and the
[Mn/Fe] trend is very similar to our work. The trend of [Mn/Fe]
almost mirrors that of [
/Fe] with respect to the
line, which suggests that type Ia SNe is a main source for the
production of Mn.
Barium is thought to be formed during the main s-process that
primarily occurs in low-mass AGB stars during He-shell burning. As
shown in Fig. 4, the [Ba/Fe] of thin and thick disk stars also shows
different behaviour. The slight declining trend of [Ba/Fe] with
increasing metallicity for thick disk stars claimed by Mashonkina et al. (2003) is not clear in our LTE results, although [Ba/Fe] of two
thick disk stars with
is larger than others with
higher [Fe/H].
Bensby et al. (2005) found the [Ba/Fe] trends are different for
the thin and thick disks. For the thick disk stars, the [Ba/Fe] is
flat, lying on a solar ratio, though the thin disk [Ba/Fe] trend
shows a prominent rise from the lowest [Fe/H] until reaching solar
metallicities, after which it starts to decline. Reddy et al.
(2003) found that [Ba/Fe] may be underabundant in thick disk stars
relative to in thin disk. Brewer & Carney (2005) found that
[Ba/Fe] is slightly enhanced in the thin disk stars compared to
the thick disk stars for
.
A series of very interesting papers on Ba abundance have been published by Mashonkina & Gehren (2000, 2001) and Mashonkina et al. (2003). Their results are obtained from NLTE differential analysis. Based on [Eu/Ba] ratio they estimate that the thick disk population formed on a timescale between 1.1 to 1.6 Gyr from the beginning of the protogalactic collapse.
Chemical abundances of 14 elements for 32 mildly metal-poor stars
were determined. According to kinematics, age, and [/Fe],
stars were identified to the thin disk, thick disk, and halo
population memberships. Element abundances for the different samples
were discussed as a function of metallicities. The main results are
summarized as follows:
The oxygen and other -element (Mg, Si, Ca, and Ti)
abundances of thin and thick disk stars show distinct trends at
.
For the odd-z elements, [Al/Fe] behaviour is exactly like an
-element, but separation for [Na/Fe] of thin and thick disk
stars is not clear.
In the iron-group elements, V, Cr, and Ni follow Fe very closely, and there is no offset between thin and thick disk stars. The Sc and Mn abundance trends of the thin and thick disk stars are different.
The [Ba/Fe] of thin disk and thick disk stars shows different behaviour.
Our abundance, kinematics, and age results of a sample of mildly metal-poor stars confirm and extend previous studies of Galactic thin and thick disk stars.
Although there has been much effort in this field in the last
decade, the abundance differences (e.g. -elements) for the
thin and thick disks stars are quite compelling, and one should note
that most of the studies have selected stars with extreme kinematics
to make it possible to classify stars as belonging to either the
thin or thick disk populations. Thus, the conclusions may be
affected by a kinematical bias. It is necessary to study large
samples in a relatively homogeneous way.
It would be very interesting to make in situ studies of the abundances and kinematics for stars in the various places of the thin and thick disks. The Large sky Area Multi-Object fiber Spectroscopic Telescope (LAMOST project) of China, which combines a large aperture with a wide field of view and unprecedented multi-objects spectrographs, provides a good chance to study a huge number of thin and thick disk stars (Zhao et al. 2005). The large scale spectroscopic survey is expected to lead to major improvements in understanding the chemical and dynamical aspects of Galactic thin and thick disk populations and should lead to a model of Galactic formation and evolution.
Acknowledgements
This research is supported by the National Natural Science Foundation of China under grant number 10433010 and NKBRSF No. G1999075406. It made use of the SIMBAD database, operated at the CDS, Strasbourg, France.