A&A 478, 487-496 (2008)
DOI: 10.1051/0004-6361:20078523

A spectroscopic study of southern (candidate) $\gamma $ Doradus stars

II. Detailed abundance analysis and fundamental parameters[*]

H. Bruntt1,2 - P. De Cat3,4 - C. Aerts4,5


1 - Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen Ø, Denmark
2 - School of Physics A28, University of Sydney, 2006 NSW, Australia
3 - Royal Observatory of Belgium, Ringlaan 3, 1180 Brussel, Belgium
4 - Katholieke Universiteit Leuven, Celestijnenlaan 200B, 3001 Leuven, Belgium
5 - Department of Astrophysics, Radboud University Nijmegen, 6500 GL Nijmegen, The Netherlands

Received 21 August 2007 / Accepted 23 October 2007

Abstract
Context. The $\gamma $ Doradus stars are a recent class of variable main sequence F-type stars located on the red edge of the Cepheid instability strip. They pulsate in gravity modes, and this makes them particularly interesting for detailed asteroseismic analysis, which can provide fundamental knowledge of properties near the convective cores of intermediate-mass main sequence stars.
Aims. To improve current understanding of $\gamma $ Dor stars through theoretical modelling, additional constraints are needed. Our aim is to estimate the fundamental atmospheric parameters and determine the chemical composition of these stars. Detailed analyses of single stars have previously suggested links to Am and $\lambda$ Boo stars, so we wish to explore this interesting connection between chemical peculiarity and pulsation.
Methods. We analysed a sample of $\gamma $ Dor stars for the first time, including nine bona fide and three candidate members of the class. We determined the fundamental atmospheric parameters and compared the abundance pattern with other A-type stars. We used the semi-automatic software package VWA for the analysis. This code relies on the calculation of synthetic spectra and thus takes line-blending into account. This is important because of the fast rotation in some of the sample stars, and we made a thorough analysis of how VWA performs when increasing $v \sin i$. We obtained good results in agreement with previously derived fundamental parameters and abundances in a few selected reference stars with properties similar to the $\gamma $ Dor stars.
Results. We find that the abundance pattern in the $\gamma $ Dor stars is not distinct from the constant A- and F-type stars we analysed.

Key words: stars: fundamental parameters - stars: abundances - stars: chemically peculiar - stars: early-type - stars: oscillations - stars: variables: general

1 Introduction

The members of the $\gamma $ Dor class of variable stars are found near the main sequence on the cool edge of the Cepheid instability strip with spectral types A7-F5. They thus share properties with $\delta $ Scuti star variables, but the $\gamma $ Dor periods are an order of magnitude longer, indicative of g mode pulsation. The $\gamma $ Dor phenomenon was first identified by Balona et al. (1994), and Krisciunas & Handler (1995) presented the first list of six candidates. Henry et al. (2007) presents a list of 66 $\gamma $ Dor stars, and the group continues to grow as new members are discovered both among field stars (De Cat et al. 2006; Mathias et al. 2004; Henry & Fekel 2005) and in open clusters (Arentoft et al. 2007). Only a few stars show pulsations characteristic of $\gamma $ Dor and $\delta $ Scuti stars simultaneously (King et al. 2007; Rowe et al. 2006; Henry & Fekel 2005).

The $\gamma $ Dor stars have given new hope for a deeper understanding of main sequence stars with masses around 2 $M_\odot$ through asteroseismic analyses. Several $\delta $ Scuti stars have been studied extensively through both photometry and spectroscopy, and dozens of individual modes are now known in a few field stars (Breger et al. 2005) and also in members of open clusters (Bruntt et al. 2007). Whilst observational work has been very successful, comparison with theoretical models has so far not been able to provide a fully adequate description of all the observations (see Zima et al. 2006, for recent developments). While the driving in $\delta $ Scuti stars is well understood in terms of the opacity or $\kappa$ mechanism, the link to predicting observed mode amplitudes is weak. Furthermore, theoretical models of $\delta $ Scuti stars show that even moderate rotation leads to significant shifts in the mode frequencies (Suárez et al. 2006a,b), which complicates the confrontation of observations and models.

The theoretical framework for interpreting the observed pulsation in $\gamma $ Dor stars is well under way. Pulsations are thought to be driven by a flux blocking mechanism near the base of their convective envelopes (Dupret et al. 2004,2006). Moya et al. (2005) investigated a method of constraining the models of $\gamma $ Dor using frequency ratios. This method has been attempted on individual $\gamma $ Dor stars (Rodríguez et al. 2006; Moya et al. 2005). While the method is indeed very useful for providing constraints, no unique models that fit all the observations were found. An improvement would be to better constrain the fundamental atmospheric parameters including metallicity. The star studied by Moya et al. (2005), HD 12901, is included in our sample.

In the current study we carry out a detailed abundance analysis of a sample of $\gamma $ Dor stars described by (De Cat et al. 2006, hereafter Paper I). Thus, the current work is the second part of our detailed spectroscopic analysis of a sample of southern candidate $\gamma $ Dor stars. In Paper I we made a detailed analysis of the spectra to study binarity and the pulsation properties. We identified 10 new bona fide $\gamma $ Dor stars of which 40% are binary stars.

Detailed abundance analyses of $\gamma $ Dor stars have only been done for a few individual stars. Bruntt et al. (2002) analysed the $\gamma $ Dor star HD 49434 and found a metallicity slightly below solar, but the analysis was hampered by the high $v \sin i$ = 85 km s-1. Sadakane (2006) analysed HD 218396 and found solar abundance of C and O (but not S) and abundances of iron peak elements of -0.5 dex, thus suggesting a $\lambda$ Boo nature for this star. Henry & Fekel (2005) found evidence that the $\gamma $ Dor star HD 8801 is an Am star based on the strength of the Ca K line. Our aim is to shed light on these intriguing links that have been suggested between the $\gamma $ Dor variables and the chemically peculiar $\lambda$ Boo and Am-type stars (Sadakane 2006; Gray & Kaye 1999).

2 Observations and selection of targets

We have obtained high-resolution spectra with the échelle spectrograph CORALIE attached to the 1.2-m Euler telescope (La Silla, Chile) for a sample of 37 known and candidate $\gamma $ Dor stars. For the details of the observations and the data reduction, we refer to Paper I. CORALIE covers the 3880-6810 Å region in 68 orders with a spectral resolution of 50 000. The typical S/N in the spectra is 100-150. For the abundance analysis, we selected the spectrum with the highest S/N. The wavelength calibrated spectra were rebinned to a step size of $\simeq$0.02 Å. Each order was normalised by fitting low-order polynomials to continuum windows identified in a synthetic spectrum. The orders were then merged to a single spectrum while making sure the overlapping orders agreed.

Stars with projected rotational velocities $v \sin i>70$ km s-1 have not been analysed due to two reasons: only very few unblended lines are available and we found that incorrect normalization of the spectra would introduce large systematic errors (Erspamer & North 2003). We also did not analyse the double-lined spectroscopic binaries from Paper I.

Table 1: Fundamental atmospheric parameters of the target stars as determined from photometric indices and parallaxes.

  
3 Aims and methods

The comparison of observed pulsation frequencies with theoretical models of $\gamma $ Dor stars will provide important insight into main sequence stars with convective cores. It is important to constrain the model space by putting constraints on $T_{\rm eff}$, $\log g$, and metallicity from observations (Moya et al. 2005). Except for one star, our targets are single stars and their fundamental parameters must be estimated by indirect methods like the calibration of photometric indices and spectroscopic analysis. In the present analyses we rely on the former as input for a detailed abundance analysis, and we seek to improve the estimates of the fundamental parameters.

  
3.1 Realistic uncertainty estimates

The fundamental physical parameters of stars are mass, radius, and luminosity. These relate to the effective temperature, $T_{\rm eff}$, and surface gravity, $\log g$. These two parameters, along with the metallicity, describe the properties of the applied model atmosphere in the abundance analysis. We adjust the parameters of the model to obtain a consistent result, specifically by requiring that the abundances measured for Fe lines of neutral and ionized species - and lines formed at different depths in the atmosphere - yield the same result. This is done by adjusting not only  $T_{\rm eff}$, $\log g$, but also the microturbulence ( $\xi _{\rm t}$), which is a crude parametrization of small-scale motions in the gas. We point out that in these analyses it is actually the temperature structure we adjust, but the final result we quote will be the $T_{\rm eff}$ and $\log g$ of the best-fitting model. It is thus not a direct measurement of the true values; therefore, it is important to assess how the fitted model parameters relate to the parameters of truly fundamental stars, i.e. stars where the properties are determined by model-independent means (Smalley 2005). The primary reference star is the Sun, but binary stars where accurate masses and radii have been determined are important for extending the validity of the method to higher temperatures. In our analysis we will test our method on the Sun and the visual binary HD 110379 (cf. Sect. 6).

Realistic uncertainties on $T_{\rm eff}$ and $\log g$ are especially important for our target stars, since some of them will likely be the targets for detailed asteroseismic studies. Our classical analysis of spectral lines yields intrinsic uncertainties on $T_{\rm eff}$ and $\log g$ in the range 50-100 K and 0.08-0.12 dex for slowly rotating stars. We estimate that at least 100 K and 0.1 dex must be added due to the limitations of the model atmosphere alone. This will lead to uncertainties on the derived abundances of the order of 0.08 dex, which must be added to the measured intrinsic scatter.

  
3.2 Model atmospheres

We used model atmospheres interpolated in the fine grid published by Heiter et al. (2002). These models are based on the original ATLAS9 code by Kurucz (1993) but use a more advanced convection description (Kupka 1996) based on Canuto & Mazzitelli (1992). Our sample consists of A- and F-type stars that have shallow convection zones.

One of the physical assumptions in the models is local thermodynamical equilibrium (LTE), but deviations from LTE start to become important for the hotter stars. We have not included the NTLE corrections in the present analysis but will estimate the importance of the effect here. According to Rentzsch-Holm (1996), the correction for neutral iron is [Fe I/H] $_{\rm NLTE}=$ [Fe I/H] $_{\rm LTE}+0.1$ dex for stars with solar metallicity and $T_{\rm eff}$ =7300 K. When this correction is applied, Fe II (unaffected by NLTE) must be increased by adding +0.2 to $\log g$. When extrapolating from Fig. 5 in Rentzsch-Holm (1996), the NTLE effect becomes negligible for stars cooler than about 6000 K.

The initial model atmosphere used for the abundance analysis of each star has $T_{\rm eff}$ from the V-K colour and $\log g$ from the HIPPARCOS estimates, except for HD 110379 where we used $\log g$ from the binary orbit (Smalley et al. 2002). For the metallicity we used the estimate from the Strömgren m1 index and solar metallicity for the two stars that did not have this index (HD 65526 and HD 218225). We note that the photometric amplitudes in variable targets are so tiny that they will not affect the applied indices. We used an initial microturbulence of 1.5 km s-1 for all stars. The adopted values for the fundamental parameters used for the initial models are printed in bold face in Table 1.

  
4 Abundance analysis with VWA

The software package VWA (Bruntt et al. 2004,2002) was used to measure abundances in the spectra and to constrain $T_{\rm eff}$ and $\log g$ for the slowly rotating stars. We have expanded VWA so it now has a graphical user interface (GUI), which allows the user to investigate the spectra in detail, pick lines manually, inspect the quality of fitted lines, etc.

Abundance analysis with VWA relies on the calculation of synthetic spectra. We use the SYNTH code by Valenti & Piskunov (1996), which works with ATLAS9 models and atomic parameters and line-broadening coefficients from the VALD database (Kupka et al. 1999). Compared to classical abundance analyses based on equivalent widths, our analysis has two important advantages:

These effects gradually become stronger when going to shorter wavelengths.

In our experience when $v \sin i$ becomes larger than about 50 km s-1, we cannot simultaneously constrain microturbulence, $T_{\rm eff}$, or $\log g$. This is because increased line blending and improper normalization of the continuum will introduce relatively large systematic errors (Erspamer & North 2003; see also Sect. 6.4).

Each line is fitted by iteratively changing the abundance to match the equivalent width (EW) of the observed and calculated spectrum. The EW is computed in a wavelength interval equal to the full-width half-maximum (FWHM) of the line. In some cases, e.g. if the line is partially blended in one wing of the line, the range for fitting the EW must be changed manually in the GUI. On a modern computer (3.2 Ghz Pentium IV), it takes about one hour to fit 250 lines for a star with low $v \sin i$. The fitted lines are inspected in the GUI, problems with the continuum level or asymmetries in the line are readily identified, and these lines are discarded. This is done automatically by calculating the $\chi^2$ of the fit in the core and the wings of the lines. This is followed by a manual inspection of the fitted lines.

An example of 12 lines fitted with VWA is shown in Fig. 1 for the star $\gamma $ Dor candidate HD 110379. The star has a moderately high projected rotational velocity of $v \sin i$ $\simeq$ 25 km s-1. It is seen that a few of the lines are affected by blends from strong neighbouring lines. As an example of the line lists we have used, we list the atomic parameters of the spectral lines we used for HD 110379 in Table B.1 in Appendix B.


  \begin{figure}
\par\includegraphics[width=17.6cm]{8523f01.eps}
\end{figure} Figure 1: Twelve Fe II lines in the candidate $\gamma $ Dor star HD 110379 fitted by VWA (continuous line). The wavelengths of the fitted lines are given in the bottom right corner of each panel.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=9.4cm]{8523f02a.eps} %
\includegraphics[width=9.4cm]{8523f02b.eps} %
\vskip 0.25cm %\end{figure} Figure 2: The two top panels show the abundances of Fe in the reference spectrum of the Sun and the target star HD 22001. The open symbols are Fe I lines and solid symbols are Fe II lines. In the bottom panel the abundances in HD 22001 are measured relative to the same lines in the Sun. As a consequence the rms scatter decreases by 40%. The solid lines are linear fits and the dashed lines are 95% confidence limits.
Open with DEXTER

4.1 Correcting log g values: relative abundances

In addition to the sample of stars in Table 1, we analysed a reference spectrum of the Sun (Hinkle et al. 2000), which has high resolution and high signal-to-noise (S/N $\simeq$ 1000). Using the results for the solar spectrum allows us to make a more precise differential abundance analysis. The abundances of Fe I and Fe II lines measured in the Sun and HD 22001 are compared in the two top panels in Fig. 2. The abundances are plotted against equivalent width and excitation potential in the left and right panels, respectively. There are 446 Fe I lines in the Sun but only 108 lines are available for HD 22001, mainly because its spectrum has lower S/N, the star has higher $v \sin i$ (i.e. fewer unblended lines), and is about 1200 K hotter than the Sun. The rms scatter of the Fe I abundance is about 0.15 dex for both stars.

It is seen that for the solar spectrum (top panels in Fig. 2), the abundance of Fe from neutral and ionized species do not agree and there is a significant positive correlation with excitation potential. The former could mean that $\log g$ is too high, while the latter indicates that the temperature of the model is too low. Since $\log g$ and $T_{\rm eff}$ are well-known for the Sun, we can make a first-order correction of the atmosphere models by measuring abundances in the target stars relative to the Sun. When doing this line-by-line any erroneous oscillator strengths, $\log gf$, are also corrected. This procedure has been used previously in detailed abundance studies of solar-like stars (Gonzalez 1998) and also stars of earlier type (Gillon & Magain 2006). To give an idea of the magnitude of the $\log gf$ corrections, we quote the rms of the corrections for a few elements: C/Sc/Ni: 0.11 dex, O/Ca/Fe lines: 0.18 dex, S/Ti/Cr: 0.23 dex, and for Si: 0.46 dex.

The result of the differential abundance analysis for HD 22001 is shown in the bottom panel in Fig. 2. It is seen that the rms scatter in the Fe I and Fe II lines is lower by about 40%. While the differential analysis improves the internal precision of the measured abundances significantly, one should note that our targets are 300-1500 K hotter than the Sun, and therefore systematic errors could be the dominant source of uncertainty on the abundances and the fundamental parameters. The amount of convection will be quite different in the sample stars compared to the Sun, and the temperature structure in these model atmospheres may not describe the observed stars correctly (Heiter et al. 2002). The fact that the rms scatter decreases significantly gives us some confidence in the differential analysis, but systematic effects on $T_{\rm eff}$ or $\log g$ could be introduced. To explore this caveat, we analysed some secondary and tertiary reference stars that have spectral types similar to our targets.

4.2 New results for the abundance in the Sun

Based on time-dependent 3D hydrodynamical models, updated atomic line parameters, and NLTE corrections, Grevesse et al. (2007) recently revised the abundances in the Sun. The overall metallicity, Z, has decreased significantly from previous estimates (Grevesse & Sauval 1998), i.e. from 0.017 to 0.0122 (-30%), mainly due to the new C, N, and O abundances. This result has vast implications in many fields of astrophysics. This includes detailed asteroseismology of $\gamma $ Dor stars, which the current work will provide important input to in terms of fundamental parameters and abundances. Our analysis was initiated before the new results, so they rely on the previous solar abundances from Grevesse & Sauval (1998). We are convinced that the analysis is still valid since it is carried out differentially with respect to an observed spectrum of the Sun. We have repeated the analysis applying the new solar abundance for one of our reference stars, HD 110379, as an explicit check. The mean abundance is 0.02 dex higher using either Fe  I and Fe  II lines, which is certainly within the uncertainty on the metallicity. The derived values for the microturbulence, $T_{\rm eff}$, and $\log g$ are unchanged.

  
5 Parameters from photometry and parallaxes

  
5.1 Strömgren indices and 2MASS (V - K)

We have derived effective temperature ( $T_{\rm eff}$), surface gravity ($\log g$), and metallicity ([Fe/H]) using both Strömgren indices and the V-K colour from the 2MASS point source catalog (Cutri et al. 2003). Furthermore, we used parallaxes from HIPPARCOS to determine $\log g$. The fundamental atmospheric parameters of the sample of $\gamma $ Dor stars we analysed are given in Table 1.

We used the TEMPLOGG software (Rogers 1995) to derive the fundamental atmospheric parameters from the Strömgren indices. The on-line version of TEMPLOGG provides uncertainty estimates that are too optimistic, since they are solely based on the uncertainties of the photometric indices. In Table 1 we quote uncertainties on $T_{\rm eff}$, $\log g$, and [Fe/H] of 250 K, 0.2 dex, and 0.1 dex (Kupka & Bruntt 2001; Rogers 1995).

We used Strömgren colour indices from the compilation of Hauck & Mermilliod (1998), but they were not available for two stars: HD 65526 and HD 218225. The H$_\beta$ index was not available for HD 12901 and HD 14940 so interstellar reddening, E(b-y), could not be determined. However, in all cases, E(b-y)<0.01 with the exception of HD 125081, which has $E(b-y)=0.045\pm0.005$. This star also has a high m1 index, indicating the star is quite metal rich. From TEMPLOGG, we get [Fe/H]  $=+0.54\pm0.10$, but this is based on an extrapolation from the calibration by Olsen (1988). The Strömgren indices for HD 110379 listed in SIMBAD are incorrect, as also noted by Scardia et al. (2007). Instead we used the average of the indices listed in Crawford et al. (1966), Cameron (1966), and Olsen (1983), which are all in good agreement. The Strömgren indices in Hauck & Mermilliod (1998) for HD 7455 are the mean of Stetson (1991) and Perry (1991), which are not in agreement: the difference in the b-y index is 0.086. Using the indices from Stetson (1991) yields $T_{\rm eff}$ =5800 K and $\log g$ =4.9, while Perry (1991) gives $T_{\rm eff}$ =6460 K and $\log g$ =4.0.

The V-K colours from the 2MASS catalogue were used to estimate $T_{\rm eff}$ with the calibrations from Masana et al. (2006). We adopted the interstellar reddening from TEMPLOGG and assumed E(b-y)=0 for the five stars where it was not available. We used E(V-K)=3.8 E(b-y) using Cardelli et al. (1989). The V-K calibration only has a weak dependence on [Fe/H] and $\log g$, so we assumed [Fe/H]  $=-0.2\pm0.2$ and $\log g$  $=4.0\pm0.3$ for all stars. This is a valid assumption since the maximum change in $T_{\rm eff}$ is 40 K when changing either [Fe/H] or $\log g$ by 2$\sigma$. The four brightest stars, HD 22001, HD 27290, HD 33262, and HD 110379, have V<5, and their 2MASS K band magnitudes are based on saturated images. For this reason the errors are large, i.e. $\sigma_K=0.23$ mag, instead of $\simeq$0.02 mag for the other stars. For these stars the uncertainty on $T_{\rm eff}$ from the V-K calibration is around 500 K, while for the other stars it is around 80-100 K. We find that $T_{\rm eff}$ from the Strömgren indices and V-K agree within the uncertainties except for HD 110379. This star is part of the visual binary system $\gamma $ Vir, and the companion, which has equal brightness, is within 4 arcsec (Scardia et al. 2007).

  
5.2 log g from HIPPARCOS parallaxes

The Hertzsprung-Russell diagram for the stars is shown in Fig. 3. It is based on the adopted $T_{\rm eff}$ and luminosities calculated from the HIPPARCOS parallaxes and a solar bolometric magnitude of $M_{\rm bol, \odot}=4.75$. In Fig. 3 we also show the $\gamma $ Dor instability strip from Dupret et al. (2005) based on models with a mixing-length parameter $\alpha=2.0$ along with evolution tracks from Lejeune & Schaerer (2001) for metallicities Z=0.008 and Z=0.02, which bracket the range for our targets. Each track is marked by the mass in solar units. From these tracks we can estimate the masses of the stars to be in the range 1.6-2.0 $M_\odot$. Since the target stars lie in the region of the "hook'' of the evolution tracks, the uncertainty on the masses is $\simeq$0.2 $M_\odot$. We have assumed a common mass of $1.8\pm0.2$ $M_\odot$ except for the well-studied binary star HD 110379, which has a known mass $M/M_\odot=1.4$ (Scardia et al. 2007).

  \begin{figure}
\par\includegraphics[width=8cm,clip]{8523f03.eps}
\end{figure} Figure 3: Hertzsprung-Russell diagram for the sample of $\gamma $ Dor candidate stars analysed. Circles, boxes, and triangle symbols are used for stars with low, moderate, and high $v \sin i$. Three evolution tracks from Lejeune & Schaerer (2001) are shown for metallicity Z=0.008 (solid lines) and one track with $M/M_\odot =2.0$ for Z=0.02 (dashed line). The dotted lines mark the $\gamma $ Dor instability strip predicted from theoretical models by Dupret et al. (2005).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=90,width=14cm,clip]{8523f04.eps}
\end{figure} Figure 4: Comparison of the spectra for the Sun and HD 110379 for which $v \sin i$ is 2 and 24 km s-1. Atomic lines selected for abundance analysis of the Sun are marked with solid vertical lines, while the neighbouring spectral lines are marked with dotted lines.
Open with DEXTER

We used this mass estimate, the adopted $T_{\rm eff}$, and the HIPPARCOS parallaxes to determine $\log g$ values. In particular we used[*] $\log g = 4 [T_{\rm eff}] + [M] + 2 \log \pi + 0.4 (V + BC_V + 0.26) + 4.44$, where $[T_{\rm eff}] = \log~ (T_{\rm eff}/T_{{\rm eff}~\odot})$ and $[M] = \log~ (M/M_\odot)$. We used bolometric corrections (BCV) from the tables by Bessell et al. (1998). If we assume the mass is known to 10% and $T_{\rm eff}$ to 4% for all stars in the sample, the uncertainty on $\log g$ will depend on the uncertainty of the parallax: 13 out of 18 stars have uncertainties below 7%, while five stars, HD 7455, HD 26298, HD 125081, HD 126516, and HD 218225, have uncertainties around 15%. The uncertainty on $\log g$ is 0.13 dex and 0.20 dex for these two groups of stars. The $\log g$ values from the Strömgren c1 index and the HIPPARCOS parallaxes agree, but the estimated uncertainty on the latter is significantly lower: typical uncertainties on $\log g$ are 0.2 and 0.1 dex, respectively.

  
6 Fundamental parameters from spectroscopy

The two important aims of the current work are to determine the fundamental parameters of the $\gamma $ Dor stars and to compare their abundance pattern with other stars of similar spectral type. We first analyse a few stars with well-known parameters. We then analyse synthetic spectra in order to estimate uncertainties and make sure our method can be used to reliably constrain $T_{\rm eff}$, $\log g$, and metallicity. Based on these results we will decide on the approach for the detailed analysis of the target stars.

VWA can automatically adjust the microturbulence, $T_{\rm eff}$, and $\log g$ either simultaneously or any parameter can be fixed. This part of the analysis is based only on Fe lines, which are the most numerous in the stellar spectra. The iterative process of adjusting the parameters is to:

Unfortunately, microturbulence and $T_{\rm eff}$ are not independent. Therefore, only the most sensitive lines are used to adjust the microturbulence (typically ${\it EW} < 80$ mÅ), while we also use stronger lines ( ${\it EW}< 150$ mÅ) to constrain $T_{\rm eff}$. For the slowly rotating stars, VWA needs to run 5-8 iterations before it converges. In each iteration, up to 120 lines are used and the CPU time is typically 1-3 h for each star.

An important limitation to detailed spectroscopic analyses arises when the lines are broadened due to rotation. In Fig. 4 we show part of the spectrum for the CORALIE spectra of the Sun and the selected of lines for the abundance analysis. The spectrum of HD 110379, which has $v \sin i$ =25 km s-1, is also shown for comparison. It is seen that line blending is worse, which illustrates that correct placement of the continuum can be difficult as $v \sin i$ increases. We will assess the importance of rotation by analysing synthetic spectra with increasing $v \sin i$ below.

  
6.1 Primary reference: solar spectrum from CORALIE

We used VWA to analyse a spectrum of the Sun measured with CORALIE. The S/N was 180, which is slightly higher than the spectra for the target stars. We ran the software with four models with parameters offset by $\pm$300 K in $T_{\rm eff}$ and $\pm$0.4 dex in $\log g$, respectively. In Table 2 we compare the results for the derived fundamental parameters. The results are very close to the canonical values of $T_{\rm eff}$ =5777 K, $\log g$ =4.44, and [Fe/H] =0.00. The quoted uncertainties on $T_{\rm eff}$ and $\log g$ were estimated as is described in Sect. 7 and the quoted uncertainty on [Fe/H] is the rms value of the abundance determined from the Fe  I lines. Uncertainties on $T_{\rm eff}$ are rounded off to 10 K. The uncertainties in Table 2 do not include any contribution from the uncertainty on the atmosphere models (cf. Sect. 3.1).

  
6.2 Secondary reference: the visual binary star HD 110379

HD 110379 is the A component in the visual binary system $\gamma $ Vir, which has two identical components (Popper 1980). From the orbital mass, the HIPPARCOS parallax, and measured spectrophotometric fluxes, constraints can be placed on $\log g$ and $T_{\rm eff}$. Thus, Smalley et al. (2002) include HD 110379 in their sample of fundamental stars and derived $T_{\rm eff}$  $=7143\pm450$ K and $\log g$  $=4.21\pm0.02$.

We used HD 110379 to test the robustness of results from VWA. The S/N in the spectrum is 120, and the star has $v \sin i$ =24 km s-1 (part of the spectrum is shown in Fig. 4). For the initial model we used the $T_{\rm eff}$ and $\log g$ estimate from the Strömgren indices and the HIPPARCOS parallax, i.e. $T_{\rm eff}$ =6860 K and $\log g$ =4.39 (cf. Table 1). We also perturbed the initial guess for the fundamental parameters and converged at the parameters listed in Table 2. One of the results (p$\log g$) is in excellent agreement with $\log g$ from the binary orbit, while the other result (p $T_{\rm eff}$) has lower values of both $T_{\rm eff}$ and $\log g$. This could be an indication that the determined values of $\log g$ and $T_{\rm eff}$ are not independent, although the results for the tertiary reference stars, except perhaps for HD 37594, do not indicate that this is a general problem.

Table 2: Fundamental parameters obtained with VWA for the reference stars after perturbing the initial value of either $T_{\rm eff}$ or $\log g$.

  
6.3 Tertiary references: HD 32115, HD 37594, HD 49933

We analysed CORALIE spectra of three slowly rotating A- and F-type stars for which detailed analyses have been published: HD 32115, HD 37594 (Bikmaev et al. 2002), and HD 49933 Gillon & Magain (2006). The fundamental parameters for these stars are constrained by photometric indices and spectroscopic analysis. Therefore, $\log g$ is known to about 0.15 dex from the parallax, which is about an order of magnitude worse than for the primary and secondary reference. The S/N in the three spectra from CORALIE is 200, 220, and 140 and the stars have $v \sin i$ =9, 17, and 10 km s-1.

As for the primary and secondary reference stars, we offset the initial parameters to test the convergence of VWA. The results are shown in Table 2. It is encouraging that for each star, the results agree within the error bars. The two slowly rotating A-type stars HD 32115 and HD 37594 were analysed by Bikmaev et al. (2002), who adopted a fixed value for $T_{\rm eff}$ based on Strömgren indices and the H$\alpha$ line and $\log g$ from the HIPPARCOS parallax. The F-type star HD 49933 was analysed by Gillon & Magain (2006), who used an approach similar to VWA to fit $T_{\rm eff}$ and $\log g$ as part of the analysis.

Our results are in acceptable agreement with previous studies. The metallicities agree within 0.1 dex, while the largest difference in $T_{\rm eff}$ and $\log g$ are 200 K and 0.2. The differences are largest for HD 32115 and HD 37594, but we recall that $T_{\rm eff}$ and $\log g$ were not adjusted as part of the abundance analysis by Bikmaev et al. (2002). On the other hand, the agreement is good for HD 49933, in which case the VWA analysis is quite similar to the approach of Gillon & Magain (2006). We recall that we did not include NTLE effects, although for the two hottest stars, the effect on $\log g$ would be about +0.2 dex. However, the studies we are comparing with here also did not include any correction.

  
6.4 Abundance analysis of synthetic spectra

We tested VWA's ability to determine $T_{\rm eff}$, $\log g$, and metallicity by using synthetic spectra with the SYNTH code. This is the "ideal'' case for abundance analysis since all $\log gf$ values are known and the spectrum is correctly normalized by design. Also, the input fundamental parameters are known. To mimic the quality of the observed data, we added random noise corresponding to S/N =100 in the continuum. We calculated spectra with $T_{\rm eff}$ =6750 and 7250 K and $\log g$ =4.3. For the cooler model, we used a range in $v \sin i$ of 10-60 km s-1, and for the hotter one, we used $v \sin i=10$ km s-1. The spectra were calculated in the range 4500-5600 Å where most of the lines are present. For slow rotation ( $v \sin i <20$), we used about 100 and 20 lines of Fe I and Fe II in the analysis. For the fast rotation (40 and 60 km s-1), only half as many lines were used.

We offset the initial models in $T_{\rm eff}$ ($\pm$500 K) or $\log g$ (+0.4 dex) and let VWA determine the best parameters. For slow and moderate rotation, $v \sin i=10$-40 km s-1, we found that the models converged satisfactorily: the largest difference in $T_{\rm eff}$, $\log g$, and [Fe/H] were 30 K, 0.05 dex, and 0.03 dex. For stars with high $v \sin i =60$ km s-1, we found $\Delta T_{\rm eff}=60$ K, $\Delta \log g=0.05$, $\Delta {\rm [Fe/H]}=0.08$.

We also calculated the uncertainties on $T_{\rm eff}$ and $\log g$ from the analysis of the synthetic spectra, and we list the results in Table 3. Uncertainties on the fundamental parameters from the models are not included (cf. Sect. 3.1). The last column gives the rms value of the Fe I and Fe II abundance. In comparison, the uncertainties for the observed secondary and tertiary reference stars listed in Table 2 are roughly twice as large. The reasons are likely a combination of imperfect continuum normalization, the remaining errors in the oscillator strengths, and differences in the temperature structure in the models and the real stars. The reference stars are all slowly rotating stars, but we may, as a first approximation, scale the uncertainties for the ideal case in Table 3 by a factor two. Thus, for stars with $v \sin i>40$ km s-1, the uncertainties become larger than the estimates from photometric indices or the HIPPARCOS parallax. We have therefore chosen not to use $T_{\rm eff}$ and $\log g$ as free parameters for stars with $v \sin i>40$ km s-1. We use this result when defining our strategy for the analysis of the target stars.

Table 3: Uncertainties on the fundamental parameters from the VWA analysis of synthetic spectra with two values of $T_{\rm eff}$ and with increasing $v \sin i$ for one value of $T_{\rm eff}$.

Table 4: The three groups of stars, depending on their $v \sin i$.

Table 5: Fundamental parameters determined with VWA for the 18 target stars and the four reference stars.

  
6.5 Not all stars are equal: low, moderate, and high v sin i

Based on the analysis of the reference stars and the synthetic spectra, we have put the stars in three groups depending on their $v \sin i$ value in Table 4: low (4-25 km s-1), moderate (35-40 km s-1), and high $v \sin i$ (50-70 km s-1). Furthermore, the stars are sorted according to the variability type from Paper I, i.e. either constant, candidate $\gamma $ Dor, bona fide $\gamma $ Dor, or bona fide $\delta $ Scuti stars.

The procedure for analysis with VWA depends on which group the star belongs to:

  
7 Fundamental parameters

In Table 5 we list the final fundamental parameters of the 18 targets stars and the four reference stars. We list the average values from Table 2 for the reference stars. Note that for the moderate and fast rotators some of the parameters were held fixed and these are marked by F in Table 5, e.g. $\log g$ from the HIPPARCOS parallax.

Uncertainties on the fundamental parameters were estimated by evaluating the sensitivity to the changes in microturbulence, $T_{\rm eff}$, and $\log g$. In Fig. 5 we show examples for HD 110379. The top panel is for the final parameters and the following panels are for increased $\xi _{\rm t}$, $\log g$, and $T_{\rm eff}$, respectively. Following Gillon & Magain (2006), the uncertainty on $T_{\rm eff}$ is found by multiplying the change in $T_{\rm eff}$ by the ratio of the uncertainty of the slope and the change in slope, s, i.e. $\sigma(T_{\rm eff}) = \Delta T_{\rm eff} \cdot {{\sigma(s)}/{\Delta(s)}}$. This uncertainty is added quadratically to the estimated uncertainty from the model atmospheres as was discussed in Sect. 3.1.

The metallicity, [M/H], is computed as the average of the five metals Ca, Sc, Ti, Cr and Fe for both neutral and ionized lines, with the requirement that at least five lines were used for any element.

In Fig. 6 we show the differences between the parameters from VWA and the initial parameters (cf. Table 1). It is seen that some of the moderate and fast rotators have $\Delta T_{\rm eff}=0$ K because they are not very sensitive to changes in $T_{\rm eff}$. For the slowly rotating stars, we find that in most cases $T_{\rm eff}$ and $\log g$ found by VWA is close to the initial model. A few exceptions are found that illustrate the importance of using more than one method to estimate the fundamental parameters.

The largest deviation is for HD 7455 where $T_{\rm eff}$ was 600 K lower and $\log g$ 0.7 dex higher than the initial model. Our result resolves the dispute over the discrepant Strömgren indices from the two different sources mentioned in Sect. 5.1. We find a large discrepancy for HD 125081, where we get a $\log g$ and metallicity that is 0.4 dex and 0.7 lower, respectively. This is the most evolved star and is also the only star with a significant interstellar reddening. If there was no interstellar reddening, [Fe/H] would be lower but not as low as we find from the abundance analysis. We find a high $T_{\rm eff}$ and high metallicity for HD 167858, but the uncertainty on $T_{\rm eff}$ is quite large.


  \begin{figure}
\par\includegraphics[width=8.8cm,clip]{8523f05.eps}
\end{figure} Figure 5: Abundance of Fe in HD 110379 for four different atmosphere models. The top panel is for the adopted parameters, while the next three panels are for increased $\xi _{\rm t}$, $T_{\rm eff}$, and $\log g$ (as indicated in the right panels). Fe I and Fe II lines are plotted with open and solid points, respectively. The solid line is a linear fit to the Fe I lines and the dashed lines indicate the 95% confidence limit of the fit.
Open with DEXTER


  \begin{figure}
\par\hskip -0.3cm
\includegraphics[width=9.6cm]{8523f06.eps}
\end{figure} Figure 6: Comparison of the fundamental parameters found from VWA and photometric indices. Circles, squares, and triangles are used for the slow, moderate, and fast rotators. Star symbols are used for the three tertiary reference stars. The HD numbers of the slow rotators are indicated in the middle panel.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=4cm,angle=90]{8523f07a.eps} \includegraphics[width=4cm,angle=90]{8523f07b.eps}
\end{figure} Figure 7: Abundance pattern in the bona fide and three candidate $\gamma $ Dor stars ( top panel) and reference stars and constant stars ( bottom panel). Only the results for slow and moderate rotators are shown. The abundances are plotted in the same order from left to right as the HD numbers in Tables A.1 and A.2. Note that the abundances are measured relative to the abundance of Fe I measured in each star.
Open with DEXTER

8 Abundances

The abundances for the bona fide $\gamma $ Dor and candidate $\gamma $ Dor stars are shown in the top panel in Fig. 7, and results for the constant and reference stars are shown in the bottom panel. Results are only shown for the slow and moderate rotators. For each element, each point corresponds to the HD numbers in the same order from left to right in Tables A.1 and A.2. Note that the abundance of each element has been offset by the abundance of Fe I, which is our primary metallicity indicator. When this offset is applied we see the abundance pattern quite clearly, especially for the reference stars where nearly all abundances lie within $0.0\pm0.2$ dex. In Appendix A we list the individual abundances for all the stars we have analysed.

8.1 Abundance pattern in the constant stars

The bottom panel in Fig. 7 shows the results for the reference stars and the constant stars in the sample. We see systematic offsets of about -0.5 to +0.2 for C, Mn, Cu, Zn, and Ba, which to some extent may be explained by the assumption of LTE. However, for C the LTE correction for stars around 7000 K is negative and of the order of -0.1 dex (Rentzsch-Holm 1996). We included only the line transitions available in VALD, while for certain elements, like Ba (McWilliam 1998), hyperfine structure is important.

For the Fe-peak elements the scatter is quite low, while the scatter from star to star is higher for the lighter elements, C, O, Na, Mg, and S. For the light elements, C to S, typically 1-5 lines are available for each element (cf. Table A.2) and so the uncertainties are quite large because of systematic errors due to erroneous continuum placement, blends, etc.

8.2 Abundance pattern in the $\gamma $ Doradus stars

The top panel in Fig. 7 shows the abundance pattern for the bona fide $\gamma $ Dor and candidate $\gamma $ Dor stars. We see systematic offsets for C, Mn, Cu, Zr, and Ba, which is also seen for some elements in the reference and constant stars. For a given element we see that the scatter is larger than for the reference stars. However, four of the eight stars have moderate $v \sin i$, so fewer lines are available.

One of our goals of the present study is to find evidence of a link between chemical peculiarity and the $\gamma $ Dor stars. In particular we searched for evidence of the following patterns:

We do not find clear evidence of these patterns in any of the stars. It was noted by Handler (1999) that all $\gamma $ Dor candidates have metallicities close to the solar value, based on the Strömgren m1 index. This is supported by our detailed spectroscopic analyses, since [M/H] lies in the range -0.27 to +0.22for 11 $\gamma $ Dor candidates (see Table 5). The exception is the fast rotator HD 218225, but its value also has the largest uncertainty, [M/H]  $=0.57\pm0.20$.

8.3 Comparison with previous studies

Abundance studies have been done previously for three of the stars in our sample, HD 48501, HD 110379, and HD 167858. Boesgaard & Tripicco (1987) analysed HD 48501 and found [Fe/H] =+0.01 and [Ca/H] =+0.20, while we get [Fe/H] =-0.08 and [Ca/H] +0.23, which is in very good agreement. Our results also roughly agree with Boesgaard & Tripicco (1986), who found high metallicity in HD 167858 at [Fe/H] =+0.15 and [Ca/H] =+0.17, while we find [Fe/H] =+0.27 and [Ca/H] +0.32. These two studies were based on relatively few lines in a limited optical range. Erspamer & North (2003) analysed several elements in HD 110379, and we have good agreement with differences below 0.1 dex, although two elements, Mg and Sc, differ by 0.2 dex.

Detailed asteroseismic modelling was attempted for the bona fide $\gamma $ Dor star HD 12901 by Moya et al. (2005). Their analysis was hampered by the uncertain metallicity of about [Fe/H] =-0.4found from the Strömgren m1 index. We find a metallicity of ${\rm [M/H]}=-0.13\pm0.16$, which we recommend using in future asteroseismic analyses of HD 12901. The star is a fast rotator with $v \sin i \simeq 64$ km s-1, so we cannot constrain $T_{\rm eff}$ and $\log g$ based on our analysis with VWA.

Two stars in our sample are included in the catalogue of Ap and Am stars compiled by Renson et al. (1991): HD 125081 is listed as a chemically peculiar star with abnormal abundances of Sr, Cr, and Eu. HD 167858 is noted as having a "doubtful nature'', but the source of this claim is not given. Paunzen & Maitzen (1998) did not find any strong chemical peculiarity in these two stars based on their measurements of the $\Delta p$ photometric index. Our present analyses of the stars support this result.

9 Conclusions

We have presented a detailed abundance analysis of a group of bona fide and candidate $\gamma $ Dor stars. In addition we analysed a number of constant stars with similar stellar parameters. There seems to be larger scatter in the abundances for the $\gamma $ Dor stars, but we find no strong evidence that the overall abundance pattern is different from other A- and F-type stars. Furthermore, the metallicity is quite close to the solar value in all cases. We have constrained the fundamental parameters of 18 single field stars from Paper I, of which about half are potential $\gamma $ Dor stars.

We also analysed a few reference stars in order to thoroughly test the performance of the VWA software package. The software gives reliable results for the value of $T_{\rm eff}$ and $\log g$ for our primary and secondary reference stars, i.e. the Sun and the astrometric binary HD 110379; the latter has a well-determined $\log g$ but poorly determined $T_{\rm eff}$. Our results also agree well with previous analyses of three tertiary reference A- and F-type stars, although these single field stars do not have well-determined values of $T_{\rm eff}$ and $\log g$. Our analysis of synthetic spectra with increasing rotational velocity shows that, for stars with $v \sin i > 50$ km s-1, our method cannot be used to constrain the microturbulence, $T_{\rm eff}$, and $\log g$ simultaneously. For the slowly rotating stars with $v \sin i < 25$ km s-1, we can constrain $T_{\rm eff}$, $\log g$, and [Fe/H] to about 120 K, 0.13 dex, and 0.09 dex including estimated uncertainties of the applied model atmospheres. These results are certainly an improvement over photometric uncertainties, which are typically at least twice as large. We expect that our results will be useful in future asteroseismic studies of $\gamma $ Dor stars.

Acknowledgements
The project was supported by the Australian and Danish Research Councils and by the Research Council of Leuven University under grant GOA-2003/04. This research has made use of the SIMBAD database, operated at the CDS, Strasbourg, France. We used atomic data extracted from the VALD data base made available through the Institute of Astronomy in Vienna, Austria. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. We are grateful to Friedrich Kupka for his very useful suggestions and to Fabien Carrier for providing the spectrum of the Sun from CORALIE. We thank the referee, Patrick François, for his very useful comments.

References

 

  
Online Material

  
Appendix A: Detailed abundance results

In Tables A.1-A.3, we list the abundances of individual elements of the target stars. The abundances are measured differentially line-by-line with respect to an observed spectrum of the Sun published by Hinkle et al. (2000). The tables list the uncertainty on the mean value (in parenthesis) and the number of lines used in the analysis. For example, the abundance of Carbon in HD 14940 is measured to be $\log N_{\rm C}/N_{\rm total}-(\log N_{\rm C}/N_{\rm total})_\odot =-0.38\pm0.03$from four lines. The quoted uncertainty is an internal value and does not include contributions from the uncertainties on the adopted fundamental parameters or shortcomings of the applied model atmosphere, which contribute by about 0.08 dex on the uncertainty of the abundances (cf. Sect. 3.1).

Table A.1: Abundances for the bona fide and candidate $\gamma $ Dor stars. The results are shown in the top panel in Fig. 7.

Table A.2: Abundances for the constant stars, the bona fide $\delta $ Scuti star HD 125081, and the tertiary reference stars (HD 32115, HD 37594, and HD 49933). The results are shown in the bottom panel in Fig. 7.

Table A.3: Abundances in the target stars with high $v \sin i$.

  
Appendix B: Line list for HD 110379

In Table B.1 we list the lines used in the abundance analysis of HD 110379 with atomic parameters extracted from the VALD database (Kupka et al. 1999). It represents a typical example of lines used in the analysis for the slowly rotating stars in our sample.

Table B.1: The atomic number, element name, wavelength, and $\log gf$ from the VALD database for the lines used in the analysis of HD 110379.



Copyright ESO 2008