A&A 475, 1003-1009 (2007)
DOI: 10.1051/0004-6361:20078233
S. Hekker1 - J. Meléndez2
1 - Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands
2 -
Research School of Astronomy & Astrophysics, Mount Stromlo Observatory, Cotter Road, Weston Creek, ACT 2611, Australia
Received 6 July 2007 / Accepted 6 September 2007
Abstract
Context. A radial velocity survey of about 380 G and K giant stars is ongoing at Lick observatory. For each star we have a high signal to noise ratio template spectrum, which we use to determine spectroscopic stellar parameters.
Aims. The aim of this paper is to present spectroscopic stellar parameters, i.e. effective temperature, surface gravity, metallicity and rotational velocity for our sample of G and K giant stars.
Methods. Effective temperatures, surface gravities and metallicities are determined from the equivalent width of iron lines, by imposing excitation and ionisation equilibrium through stellar atmosphere models. Rotational velocities are determined from the full width at half maximum (FWHM) of moderate spectral lines. A calibration between the FWHM and total broadening (rotational velocity and macro turbulence) is obtained from stars in common between our sample and the sample from Gray (1989, ApJ, 347, 1021). Macro turbulence is determined from the macro turbulence vs. spectral type relations from Gray (2005, The Observation and Analysis of Stellar Photospheres).
Results. The metallicity we derive is essentially equal to the literature values, while the effective temperature and surface gravity are slightly higher by 56 K and 0.15 dex, respectively. A method comparison is performed with 72 giants in common with Luck & Heiter (2007, AJ, 133, 2464), which shows that both methods give similar results. Our rotational velocities are comparable with those obtained by Gray (1989, ApJ, 347, 1021), but somewhat higher than those of de Medeiros & Mayor (1999, A&AS, 139, 433), which is consistent with the different diagnostics used to determine them.
Conclusions. We are able to determine spectroscopic stellar parameters for about 380 G and K giant stars uniformly (112 stars are being analysed spectroscopically for the first time). For stars available in the literature, we find reasonable agreement between literature values and values determined in the present work. In addition, we show that the metallicity enhancement of companion hosting stars might also be valid for giant stars, with the planet hosting giants being 0.13
0.03 dex (i.e. 35
10%) more metal rich than our total sample of stars.
Key words: stars: abundances - stars: fundamental parameters - stars: rotation - methods: observational - techniques: spectroscopic
For the determination of spectroscopic stellar parameters, one needs high resolution spectra with high signal to noise ratio. These spectra are available from radial velocity surveys and are often used to determine stellar parameters. For instance, properties of cool stars from the Keck, Lick and AAT planet search are described by Valenti & Fischer (2005). Atmospheric parameters for stars observed by the N2K consortium (Fischer et al. 2005) are described by Robinson et al. (2007a). Santos et al. (2005,2004) present stellar parameters and metallicities from the planet search using ESO facilities and the ELODIE spectrograph at the 1.93 m telescope at the Observatoire de Haute Provence. Also, basic stellar parameters for 72 evolved stars, previously studied for radial velocity variations, are presented by da Silva et al. (2006). Some of these results are not only interesting in terms of the stellar parameters, but also reveal which stars are most likely to harbour sub-stellar companions. As first shown by Gonzalez (1997), and confirmed with larger samples by Fischer & Valenti (2005) and Santos et al. (2005), metal rich stars are more likely to harbour companions than metal poor ones.
Spectroscopic stellar parameters are most commonly determined by fitting the observed spectrum directly, see for instance Valenti & Fischer (2005), or by imposing excitation and ionisation equilibrium for metal lines, using an LTE analysis and a grid of model atmospheres, see for instance Santos et al. (2005,2004), da Silva et al. (2006), Takeda et al. (2002) and Luck & Heiter (2007).
Rotational velocity and macro turbulence can only be determined directly with the Fourier transform technique, see for instance Gray (1989). Benz & Mayor (1981) have shown that accurate rotational velocities can also be deduced for dwarfs from a cross correlation function, by performing a calibration with the direct measurements of Gray (1989). de Medeiros & Mayor (1999) extended this technique for giant stars. Fekel (1997) used the full width at half maximum (FWHM) of weak to moderate spectral lines to determine rotational velocities, also by performing a calibration with the results of Gray (1989).
In 1999, a radial velocity survey of about 180 K giant stars was started at UCO/Lick Observatory, USA. This ongoing survey has recently been expanded to about 380 G and K giants. From the initial sample of 180 stars, companions have been announced for Draconis (Frink et al. 2002) and Pollux (Reffert et al. 2006). Stars with radial velocity variations of less than 20 ms-1 were presented as stable stars by Hekker et al. (2006), and an investigation into the mechanism(s) causing the radial velocity variations is presented by Hekker et al. (2007). Some binaries discovered with this survey, as well as an extensive overview of the sample, will be presented in forthcoming papers.
In this paper, we determine stellar parameters, i.e. effective temperature ( ), surface gravity () and metallicity ([Fe/H]), as well as rotational velocity ( ) for all stars in the sample. In Sect. 2, we describe the observations. In Sects. 3 and 4, we present the methods used, and results for the stellar parameters and rotational velocity, respectively. In Sect. 5 a summary of our results is presented.
For the radial velocity survey, giants were selected from the Hipparcos catalog (Perryman & ESA 1997), based on the criteria described by Frink et al. (2001). The selected stars are all brighter than 6 mag, are presumably single and have photometric variations <0.06 mag. These criteria are the same for the initial sample (K1 and later giants) as well as for the extension (G and K0, K1 giants). The survey started in 1999 at Lick observatory using the Coude Auxiliary Telescope (CAT) in conjunction with the Hamilton echelle spectrograph (R = 60 000). The radial velocity measurements are performed with an iodine cell in the light path as described by Marcy & Butler (1992) and Valenti et al. (1995). Radial velocities are determined from the comparison of a stellar spectrum obtained with an iodine cell in the light path, and the convolution of a template iodine spectrum and a template stellar spectrum obtained without an iodine cell in the light path (Butler et al. 1996). For each target star we have a high signal to noise ratio template spectrum. These templates are used for the determination of the stellar parameters described in this paper. Thorium-argon images taken at the beginning and end of each night are used for wavelength calibration.
Table 1: Iron lines considered in our analysis.
Spectroscopic stellar parameters ( , and [Fe/H]) are determined by measuring the equivalent width (EW) of iron lines. The iron lines used in this work are listed in Table 1. The lines were carefully selected to avoid blends by atomic and CN lines. CN blends were visually inspected by comparing a synthetic spectrum computed with laboratory CN lines (Meléndez & Barbuy 1999) with the high resolution visible atlas of the cool giant Arcturus (Hinkle et al. 2000). The values are based on laboratory works, in some cases with small adjustments using the Arcturus atlas. For Fe I, they are from the Oxford group (e.g. Blackwell et al. 1995), Hannover group (e.g. O'Brian et al. 1991; May et al. 1974; Bard & Kock 1994; and Milford et al. 1994). For Fe II, the values are from the laboratory normalisation performed by Meléndez et al. (2006).
Figure 1: EWs obtained with the ARES software as a function of EWs computed using IRAF. The dots indicate a star with = 4050 K ("cool'') and the asterisks a star with = 4900 K ("hot''). The solid line is a1 to 1 relation. | |
Open with DEXTER |
It is very time consuming to determine EWs for about 380 stars by hand, using for instance the "splot'' routine from IRAF, we therefore used the publicly available Automatic Routine for line Equivalent widths in stellar Spectra (ARES) (Sousa et al. 2007). To check for possible differences between the EWs determined with ARES and those obtained with IRAF, we plot the EWs obtained with ARES vs. those obtained using IRAF (see Fig. 1). This comparison is done for a "hot'' ( = 4900 K) and a "cool'' ( = 4050 K) star in the sample. The mean differences between the EWs measured with ARES and IRAF are: = 1.6 mÅ and 2.2 mÅ, with standard deviations of 6.6 mÅ and 2.6 mÅ, for the "cool'' and "hot'' star respectively. The nearly 1 to 1 relation between the EWs obtained with both methods shows that it is reasonable to use EWs obtained with ARES. For some stars one or more lines appeared to be too strong (stronger than 200 mÅ) for a reliable parameter estimate. These lines are discarded.
Figure 2: Difference between our spectroscopic values and literature values as a function of literature values for the effective temperature ( left), the logarithm of the surface gravity in cgs units ( centre) and metallicity ( right). | |
Open with DEXTER |
From the EWs, stellar parameters are determined by imposing excitation and ionisation equilibrium through stellar atmosphere models. The micro turbulence (vt) was obtained by requiring no dependence of Fe I against equivalent width. We performed a spectroscopic LTE analysis using the 2002 version of MOOG (Sneden 1973) and Kurucz model atmospheres, which include overshooting (Castelli et al. 1997). The resulting stellar parameters for each star are listed in Table 4 (only available in the online version). The reference solar iron abundance used in this study is A(Fe) = 7.49 and was obtained using the same grid of Kurucz models.
Based on the scatter of the Fe lines, we obtained observational errors (standard errors) of 0.024 dex for Fe I and 0.037 dex for Fe II. Furthermore, we compute internal errors for a change of +80 K, +0.2 dex and +0.2 km s-1 in , and micro turbulence, respectively. The errors due to change in stellar parameters are shown for three stars with different temperatures in Table 2.
We compare our spectroscopic stellar parameters with values obtained from the literature. These literature values are from an updated version (Ramírez & Meléndez 2005) of the Cayrel de Strobel et al. (2001) [Fe/H] catalogue (254 stars in common), including the Luck & Heiter (2007) catalogue. If possible, suspicious literature values were corrected according to the normalisation suggested by Taylor (1999). In Fig. 2 the difference between our spectroscopic and literature values are plotted for the temperature, surface gravity and metallicity. We find the following trimean difference and pseudo-sigma for the stellar parameters:
Table 2: Internal errors due to changes in stellar parameters for three stars with different temperatures. The total error in the last column includes the observational errors (0.024 dex for Fe I and 0.037 dex for Fe II).
The difference in our spectroscopic metallicity and the literature value is essentially zero and we conclude that our metallicity scale is correct. Furthermore, our values are in good agreement with the literature, with a scatter of only 84 K, and a zero point difference of 56 K, our values being higher than the values in the literature. From the left panel of Fig. 2 one can see that the difference is largest for the coolest stars in the sample. This might be due to the fact that the models are less accurate for low temperatures. In addition, the number of spectral lines increases with decreasing temperature, the spectra might be too crowded at lower temperatures, and also the lines get stronger and more dependent on the micro turbulence. Our results below 4000 K should be interpreted with caution.
The spectroscopic gravities we derived also agree well with the literature, with a scatter of only 0.22 dex, and a zero point difference of 0.15 dex. We checked whether the enhanced values from our spectroscopic analyses are related to the higher temperatures we obtained, compared to literature values. We therefore performed a test for three stars with 4170 K, 4445 K and 4980 K, respectively. We increased with 100 K and determined , while keeping the micro turbulence fixed. For all three stars we obtained higher values for the increased temperatures. This reveals that the higher values for , compared to the literature values, are related to the higher effective temperatures.
Figure 3: Difference between our spectroscopic values and values from Luck & Heiter (2007) as a function of Luck & Heiter (2007) values for the effective temperature ( left), the logarithm of the surface gravity in cgs units ( centre) and metallicity ( right) for the 72 stars in common between the samples. | |
Open with DEXTER |
Table 3: Companion hosting giants with their mean metallicities and literature sources.
Recently, Luck & Heiter (2007) presented a homogeneous spectroscopic analysis of 298 giants in the local region, using between 300 and 400 Fe I lines for each star and MARCS stellar models (Gustafsson et al. 2003) on spectra with R=60 000. We used the 72 stars in common between Luck & Heiter (2007) and our sample, to see how well we can determine spectroscopic stellar parameters with just two dozen carefully selected iron lines instead of a few hundred iron lines. Note that different models are used for the two analyses, however spectra have the same resolution. A comparison for each parameter is shown in Fig. 3.
We find the following trimean difference and pseudo-sigma:
By now 20 sub-stellar companions are announced as orbiting giant stars. Some recent publications, e.g. Sadakane et al. (2005) and Pasquini et al. (2007), suggest that giant stars with companions are metal poor, which is quite different from the known metallicity enhancement in dwarf stars hosting companions. Schuler et al. (2005) and da Silva et al. (2006) argue that giant stars with companions may be metal poor, due to a stellar mass-companion relation instead of a metallicity-companion relation. Indeed, Fischer & Valenti (2005) also find a relation between stellar mass and companions, but conclude that this is likely spurious. Here, we look at the metallicities of the giants with announced companions and compare these with the metallicities of the giants in our sample. For consistency with our work, the companion hosting giants are selected to encompass the same range in stellar parameters as our sample of giants, i.e. stars with 750 K and 2.1 1.5.
The companion hosting giants and their metallicities adopted for the comparison are listed in Table 3. In order to perform a homogeneous comparison between the metallicities of the total and the companion hosting giant sample we determined the metallicity zero-points of the literature works we are using for the companion hosting giant stars. We found the following zero-points: -0.05 dex for Luck & Heiter (2007), -0.10 dex for Schuler et al. (2005), -0.02 dex for da Silva et al. (2006) and Doellinger et al. (2007), -0.07 dex for Santos et al. (2003,2005) and Sousa et al. (2006), -0.11 dex for Johnson et al. (2007), Valenti & Fischer (2005) and Robinson et al. (2007b), 0.00 dex for Mishenina et al. (2006), -0.05 dex for Sadakane et al. (2005), Sato et al. (2007) and Takeda et al. (2005), and -0.85 dex for Twarog et al. (1997). Note that Twarog et al. (1997) only provide a mean metallicity for the cluster NGC 2423 based on photometry. The normalisation factor is based on the Hyades. For Niedzielski et al. (2007) we do not have any stars in common, so we adopted a zero-point of +0.01 dex, which we obtained from the global comparison with the literature. In the event that more individual measurements of a star were available in the literature, the mean metallicity was adopted. The literature sources are mentioned in the last column of Table 3.
Figure 4: Distribution of metallicities of all stars in our sample with mean, median and trimean values of -0.12, -0.09 and -0.095 dex, respectively. The metallicity of 20 giants with an announced companion in the literature (see text for selection criteria) are plotted with the dashed histogram. The latter distribution is normalised to the total number of giants in our sample. These giants have mean, median and trimean metallicity values of -0.05, +0.025 and -0.015 dex, respectively. | |
Open with DEXTER |
In Fig. 4 the metallicity distribution of all stars in the present sample is shown together with the metallicity distribution of giant stars with announced companions. The companion hosting giant star distribution shows a gap at -0.1 dex. We do not know whether this gap is real or due to low number statistics. The mean, median and trimean metallicities are -0.12, -0.09 and -0.095 dex and -0.05, +0.025 and -0.015 dex for the total and companions hosting sample, respectively, while the peaks of the histograms are at -0.05 dex and +0.05 dex. Gaussians fitted through the two distribution have their centres at -0.06 dex and +0.09 dex (+0.085 dex in case one Gauss is fitted and +0.09 for 2 Gaussian fits) for the total and companion hosting sample, respectively. Therefore, the metallicity enhancement for companion hosting giants is 0.13 0.03 dex. This is similar to the metallicity enhancement found by Fischer & Valenti (2005). Their comparison between metallicities of all stars in the sample and companion hosting dwarfs reveals that companion hosting dwarfs are more metal rich by 0.13 dex. If they compare the metallicity enhancement as a function of stellar mass, they find also that, independent of mass, the metallicity distribution of dwarfs with companions is 0.12 dex higher than the average metallicity of all stars in the sample.
Although the determination of precise stellar masses, for both our sample of field giant stars and the companion hosting giant stars, is beyond the scope of the present paper, we stress that overall both samples have comparable masses. Using the masses given by Allende Prieto & Lambert (1999), we estimate a typical mass of 2.0 for our sample, with the bulk of the field giant sample in the range 1.4-2.9 (first to third quartiles), which is the typical range covered by companion hosting giants according to Sadakane et al. (2005) and Johnson et al. (2007).
The metallicity enhancement for companion hosting giants should be taken with caution. First, it is still based on low number statistics. Second, for nearly 200 stars in our sample, we do not have a long enough time span of radial velocity observations to detect companions, in case these are present. Third, there is still some discussion ongoing about some of the companion hosting giants with companions in (nearly) circular orbits. The observed radial velocity variations could in principle also be due to a mechanism intrinsic to the star.
The present conclusion that companion hosting giants are on average metal rich is rather different compared to other studies. Sadakane et al. (2005) and Pasquini et al. (2007) both agree that the giant stars with companions are typically not much more metal rich than [Fe/H] = 0.0. Sadakane et al. (2005) analysed only a few giants, therefore their results may be due to statistics based on small numbers. However, Pasquini et al. (2007) used a total of 14 giants, slightly less than our sample of 20 giants. The two samples have only 10 stars in common (HD 11977, HD 13189, HD 28305, HD 47536, HD 62509, HD 73108, HD 104985, HD 122430, HD 137759, HD 222404), which include the most metal poor giants in the sample, but does not contain some of the more metal rich ones. Since their histogram for the companion hosting giants does not include stars with [Fe/H] > +0.2 dex, we suspect that the four candidates which are not available in the literature are also not as metal rich as for example HD 27442 and HD 177830, both with [Fe/H] > +0.3 (Santos et al. 2003; Valenti & Fischer 2005). Note that Pasquini et al. (2007) did not correct for potential systematic errors in the literature [Fe/H] values for giant stars with companions.
Finally, caution should be taken when comparing dwarfs with giants, since systematic errors are expected due to differences in their atmospheric structures and parameters (Asplund 2005). Indeed, in a large study of several hundreds of giants and dwarfs, Luck & Heiter (2007) show that the abundances of Fe, Mn and Ba in giants may be affected by systematic uncertainties when compared with dwarf stars. A consistent study such as ours, comparing field giants to companion hosting giants, is therefore more reliable.
We computed rotational velocities for our sample of giant stars, using the method described by Fekel (1997). The FWHM for moderate spectral lines at 6432.68, 6452.33, 6454.99, 6455,60, 6456.38, 6469.15 and 6471.66 Å is determined and averaged. The dispersion of the FWHM resulting from individual lines is typically 0.025 mÅ. The instrumental broadening is determined from thorium-argon (ThAr) images taken at the beginning and end of each night. The FWHM of several ThAr lines, in the same spectral region as the stellar lines, are determined and averaged. The intrinsic stellar broadening is computed as .
The intrinsic stellar broadening is converted to rotational velocity
,
using the results from Gray (1989). For the 51 stars in common (excluding 2 outliers, for which we find higher velocities than Gray 1989), the intrinsic broadening is plotted as a function of the total broadening (
)
determined by Gray (1989), as shown in Fig. 5. A second order polynomial is fitted:
Figure 5: The intrinsic stellar FWHM of spectral lines as a function Gray's total broadening, for 51 stars in common between our sample and Gray (1989). The best fit (Eq. (1)) is shown as the solid line, while the best fit obtained by Fekel (1997) is shown as the dashed line. | |
Open with DEXTER |
Figure 6: Log g vs. for all stars in our sample. The solid line shows the best fit (Eq. (2)). | |
Open with DEXTER |
The macro turbulence is derived from the spectral type as shown in Fig. 17.10 from Gray (2005). Each luminosity class has its own relation. We estimate the luminosity class from a vs.
relation (Fig. 6). Most stars in the sample are luminosity class III stars, and therefore, the second order best fit relation, shown in Eq. (2), is used for class III stars. This relation has a robust sigma scatter of 0.25 dex.
(3) | |||
(4) | |||
(5) |
We checked our final values by comparing the values of the 51 stars in common between our sample and Gray (1989), see Fig. 7. The values are located around the 1 to 1 relation indicated by the solid line, which shows that the results of both samples are consistent. We also have 184 stars in common with de Medeiros & Mayor (1999) and compare our values with theirs in Fig. 8. Our values are on average higher than those obtained by de Medeiros & Mayor (1999). This is probably due to the different diagnostics used. de Medeiros & Mayor (1999) show that the relation between their values, and those obtained by Gray (1989) for class III and IV, has an offset of 1.15 and a correlation coefficient of 1.18. We plotted the 1 to 1 relation, solid line, as well as the relation between obtained by Gray (1989) and de Medeiros & Mayor (1999) in Fig. 8, dashed line. The data are located around this latter relation. This indicates that the difference between the results obtained here and from de Medeiros & Mayor (1999) are due to the different diagnostics used to determine . Also, Luck & Heiter (2007) find that the CORAVEL values may suffer from systematic differences with respect to values derived from other techniques. For all stars and are listed in Table 4 (only available at the CDS).
Figure 7: obtained here vs. obtained by Gray (1989). The solid line is a 1 to 1 relation. | |
Open with DEXTER |
Figure 8: obtained here vs. obtained by de Medeiros & Mayor (1999). The solid line is a 1 to 1 relation and the dashed line is the relation between obtained by Gray (1989) and de Medeiros & Mayor (1999). The arrows indicate upper limits. | |
Open with DEXTER |
We have determined spectroscopic stellar parameters for a sample of about 380 G and K giant stars. Among these, 112 stars are analysed spectroscopically for the first time. Our metallicities agree with values found in the literature and we conclude that our metallicity scale is not severely affected by systematic errors. Our temperatures are 50 K higher compared to those from the literature. The difference is largest for stars with lowest temperatures. This is probably due to the lower accuracy of atmosphere models in this temperature range, the increased number and strength of spectral lines and increasing dependence on micro turbulence in cooler stars. An increase in temperature causes an increase in surface gravity and our values are 0.15 dex higher compared to the literature values.
The comparison between the mean metallicity of our total sample of giant stars and giant stars with announced companions reveals that the companion hosting stars have a 0.13 0.03 dex higher metallicity than the mean metallicity of our total sample. This is in agreement with the enhanced metallicity of companion hosting dwarf stars, but is based on low number statistics.
Rotational velocities are determined using the method described by Fekel (1997). Stars in common between our sample and that observed by Gray (1989) are used to convert FWHM of moderate lines [Å] to total line broadening [km s-1]. We used a vs. correlation to determine the luminosity class of the stars. This luminosity class was subsequently used to calculate the macro turbulence, which has a different relation with temperature for different classes. Our data are in agreement with those obtained by Gray (1989), but are on average larger than the values obtained by de Medeiros & Mayor (1999). This is due to the different diagnostics used to determine .
Acknowledgements
We are very grateful to Andreas Quirrenbach, Sabine Reffert and Dave Mitchell, who started the radial velocity survey of giant stars, to allow us to use the template spectra for the present analysis. We also want to thank Conny Aerts and Ignas Snellen for carefully reading the manuscript and useful discussions. We would like to thank the referee for carefully reading the manuscript and constructive comments.