A&A 436, 253-262 (2005)
DOI: 10.1051/0004-6361:20042075

VLTI/VINCI diameter constraints on the evolutionary status of $\delta $ Eri, $\xi $ Hya, $\eta $ Boo

F. Thévenin1 - P. Kervella2 - B. Pichon 1 - P. Morel 1 - E. Di Folco3 - Y. Lebreton4


1 - Laboratoire Cassiopée, UMR 6202 CNRS, Observatoire de la Côte d'Azur, BP 4229, 06304 Nice Cedex 4, France
2 - Laboratoire d'Études Spatiales et d'Instrumentation Astrophysique (LESIA), UMR 8109 du CNRS, Observatoire de Paris, Section de Meudon, 5 place Jules Janssen, 92195 Meudon Cedex, France
3 - European Southern Observatory, Karl-Schwarzschild Strasse 2, 85748 Garching, Germany
4 - GEPI, UMR 8111 CNRS, Observatoire de Paris, Section de Meudon, 5 place Jules Janssen, 92195 Meudon Cedex, France

Received 27 September 2004 / Accepted 16 January 2005

Abstract
Using VLTI/VINCI angular diameter measurements, we constrain the evolutionary status of three asteroseismic targets: the stars $\delta $ Eri, $\xi $ Hya, $\eta $ Boo. Our predictions of the mean large frequency spacing of these stars are in agreement with published observational estimations. Looking without success for a companion of $\delta $ Eri, we doubt its classification as an RS CVn star.

Key words: stars: evolution - stars: fundamental parameters - techniques: interferometric

1 Introduction

After two years of operation, the commissioning instrument VINCI of the VLTI has provided valuable stellar diameter measurements. Among the impact of these diameters are the studies of main sequence stars, where diameters combined with asteroseismic frequencies can be used to constrain evolutionary status and mass. Several papers have been published (Ségransan et al. 2003; Kervella et al. 2003a,b; 2004a; Di Folco et al. 2004) with important results on stellar fundamental parameters prior to the use of the dedicated VLTI light combiner AMBER (Petrov et al. 2003). The aim of the present paper is to complete previous studies using VINCI to measure the diameter of three subgiant and giant stars that are among selected asteroseismic targets for ground-based observations and space missions: $\delta $ Eri, $\xi $ Hya, $\eta $ Boo. We perform a preliminary study of their evolutionary status by constraining their mass, their helium content and their age. One of the purposes of this paper is to show that in the future, the use of stellar diameters will be a significant constraint for evolutionary models for given input physics. We first describe the characteristics of each of the three stars (Sect. 2) and then we present diameter measurements (Sect. 3) for each star. We construct evolutionary models satisfying spectro-photometric observable constraints and we compare asteroseismic large frequencies with measured ones. We present these models (Sect. 4) and we draw some conclusions on the classification and fundamental parameters of the three stars.

   
2 Global characteristics of the stars

The first part of Table 1 presents the observational data of the three stars. The second part of this table summarizes some input parameters and output data of the models.

Table 1: Observable characteristics of the stars and best model reproducing them. The subscripts " $_{\rm ini}$'' and " $_{\rm surf}$'' respectively refer to initial values and current surface quantities. Note that the presented errors of VLTI/VINCI angular diameters are the statistical ones followed by the systematical ones. Note also that $D/D_\odot $ is equal to $R/R_\odot $.

2.1 $\delta $ Eri

$\delta $ Eri (HD 23249, HR 1136, HIP 17378) has been thoroughly studied by photometry and spectroscopy and is classified as a K0 IV star (Keenan & Pitts 1980). It belongs to the group of the nearest stars with an accurate Hipparcos parallax of $ \rm 110.58 \pm 0.88 ~ mas$ (Perryman et al. 1997). The star has been classified as a weakly active and X-ray soft source (Huensch et al. 1999) after a lengthy search for its activity. Wilson & Bappu (1957) concluded that the possible detection of emission in the H&K lines is "exceedingly weak'' - so weak that it is questionable. It took more than 20 years to inconclusively detect its activity with Copernicus, revealing a weak emission in MgII (Weiler & Oegerle 1979). Fisher et al. (1983) tried to detect a periodic variation in the photometric data and concluded that, if it exists, the amplitude is below $ \pm $0.02 mag with a period of 10 days. They suggested that $\delta $ Eri could be classified as an RS CVn star. An RS CVn is defined as a F-G binary star having a period shorter than 14 days, with chromospheric activity and with a period of rotation synchronized with its orbital period (Linsky 1984), giving the star high rotational velocity inducing strong activity. This is in contrast with the low level of activity detected for $\delta $ Eri making doubtful its classification as an RS CVn star. $\delta $ Eri has a projected rotational velocity of $v \sin i = 1.0 ~ \rm km ~ s^{-1}$ (de  Meideros & Mayor 1999) and the hypothetical RS CVn classification forces us to conclude that the binary is seen pole-on therefore explaining the lack of photometric variation and also of any variation of the radial velocity (Santos et al. 2004). To reveal the presence of a close companion around $\delta $ Eri, we set several VLTI/VINCI observations at different baselines (see Sect. 3).

We estimate its bolometric luminosity $L_{\star}/L_{\odot} = 3.19 \pm 0.06 $ using the Alonso et al. (1999) empirical bolometric corrections (BC, $\rm BC = -0.24 \pm 0.01 $ for giants, this is the dominant source of uncertainty on the luminosity). We adopt the Santos et al. (2004) values for the effective temperature $ T_{\rm eff} = 5074 \pm 60 ~ \rm K $, logarithmic surface gravity $\log g = 3.77 \pm 0.16 $ and surface iron abundance $ \rm [Fe/H] = 0.13 \pm 0.03 $. These parameters are different from - but within the error bars of - the parameters proposed by Pijpers (2003) for this star, except the metallicity which is 0.24 dex higher. Bouchy & Carrier (2003) have measured a mean large frequency spacing of $ \rm 43.8~ \mu Hz $ that we will try to reproduce with our model. We recall that the large frequency spacing is defined as the difference between frequencies of modes with consecutive radial order n: $\Delta \nu_l (n) = \nu_{n,l}- \nu_{n-1,l} $. In the high frequency range, i.e. large radial orders, $\Delta \nu_l (n) $ is almost constant with a mean value strongly related to the square root of the mean density of the star. To obtain the mean large frequency separation, we average over l = 0 - 2.

2.2 $\xi $ Hya

$\xi $ Hya (HD 100407, HR 4450, HIP 56343) is a giant star (G7 III) which has been considered by Eggen (1977) as a spurious member of the Hyades group because it departs slightly from the regression line of giant stars in the colour diagrams (b-y, R-I) and (M1, R-I) of that stellar group.

Its Hipparcos parallax is $ \rm 25.23 \pm 0.83 ~ mas $. We estimate its bolometric luminosity $L_{\star}/L_{\odot} = 60.7 \pm 4.1 $ using BC ( $\rm BC = -0.26 \pm 0.01 $) from Alonso et al. (1999). We adopt the spectroscopic parameters derived by Mc William (1990): effective temperature $ T_{\rm eff} = 5010 \pm 100 ~ \rm K $, $ \log g = 2.93 \pm 0.30 $ and $ \rm [Fe/H] = -0.04 \pm 0.12 $. These parameters are different from - but within the error bars of - the parameters adopted by Frandsen et al. (2002) for this star. The star belongs to the HR diagram at the lowest part of the giant branch corresponding to an evolved star with a mass around $ 3~ M_{\odot} $. Using a set of CORALIE spectra, Frandsen et al. (2002) detected solar-like oscillations suggesting radial modes with the largest amplitudes almost equidistant around $ \rm 7.1~ \mu Hz $. That important detection opens the possibility of better constraining the model of this star for which the mass is not well-known.

2.3 $\eta $ Boo

$\eta $ Boo (HD 121370, HR 5235, HIP 67927) is a subgiant (G0 IV) spectroscopic binary (SB1) studied recently by Di Mauro et al. (2003, 2004) and Guenther (2004). Its Hipparcos parallax is $ \rm 88.17 \pm 0.75 ~ mas $. Having large overabundances of Si, Na, S, Ni and Fe, it has been considered as super-metal-rich by Feltzing & Gonzales (2001). We adopt here a luminosity $ L_{\star}/L_{\odot} = 8.95 \pm 0.20 $ using BC ( $\rm BC = -0.06 \pm 0.01 $, this is the dominant source of uncertainty on the luminosity) from Vandenberg & Clem (2003) for this subgiant, with an effective temperature $ T_{\rm eff} = 6050 \pm 150 ~ \rm K $ representing the average of five effective temperature determinations in the [Fe/H] catalogue of Cayrel de Strobel et al. (2001) and the spectroscopic $ \log g = 3.66 \pm 0.20 $ and $ \rm [Fe/H] = 0.24 \pm 0.07 $ from Feltzing & Gonzales (2001). These parameters are different from - but within the error bars of - the parameters adopted by Di Mauro et al. (2003, 2004) for this star. Asteroseismic observations of $\delta $ Eri have been reported by Carrier et al. (2005) with $ \rm\Delta \nu_0 = 39.9 \pm 0.1 ~ \mu Hz $ and by Kjeldsen et al. (2003) with $ \rm\Delta \nu_0 = 40.47 \pm 0.05 ~ \mu Hz $.

   
3 Diameter interferometric measurements

3.1 VINCI and the VLTI

The European Southern Observatory's Very Large Telescope Interferometer (Glindemann et al. 2000) has been operated on top of the Cerro Paranal, in Northern Chile since March 2001. For the observations reported in this work, the light coming from two telescopes (two 0.35 m test siderostats or VLT/UT1-UT3) was combined coherently in VINCI, the VLT Interferometer Commissioning Instrument (Kervella et al. 2000). We used a regular K band filter ( $ \rm\lambda = 2.0 {-} 2.4 ~ \mu m$) for these observations.

3.2 Data reduction

We used an improved version of the standard VINCI data reduction pipeline (Kervella et al. 2004b), whose general principle is based on the original FLUOR algorithm (Coudé du Foresto et al. 1997). The two calibrated output interferograms are subtracted to remove residual photometric fluctuations. Instead of the classical Fourier analysis, we implemented a time-frequency analysis (Ségransan et al. 1999) based on a continuous wavelet transform.

The atmospheric piston effect between the two telescopes corrupts the amplitude and the shape of the fringe peak in the wavelet power spectrum. As described in Kervella et al. (2004b), the properties of the fringe peaks in the time and frequency domains are monitored automatically, in order to reject from the processing the interferograms that are strongly affected by the atmospheric piston. This selection reduces the statistical dispersion of the squared coherence factors ($\mu^2$) measurement, and avoids biases from corrupted interferograms. The final $\mu^2$ values are derived by integrating the average wavelet power spectral density (PSD) of the interferograms at the position and frequency of the fringes. The residual photon and detector noise backgrounds are removed using a linear least squares fit of the PSD at high and low frequency. The statistical error bars on $\mu^2$ are computed from the series of $\mu^2$ values obtained on each target star (typically a few hundred interferograms) using the bootstrapping technique.

   
3.3 Measured visibilities and angular diameters

The visibility values obtained on $\delta $ Eri, $\xi $ Hya and $\eta $ Boo are listed in Tables 2 to 5, and plotted in Figs. 1 to 3.

The calibration of the visibilities obtained on $\delta $ Eri and $\eta $ Boo was done using well-known calibrator stars that were selected from the Cohen et al. (1999) catalogue. The uniform disk (UD) angular diameter of these stars was converted into a limb darkened value and then to a K band uniform disk angular diameter using the recent non-linear law coefficients taken from Claret et al. (2000). As demonstrated by Bordé et al. (2002), the star diameters in this list have been measured very homogeneously to a relative precision of approximately 1%.

Table 2: $\delta $ Eri squared visibilities.

Table 3: $\delta $ Eri squared visibilities (continued from Table 2).

Table 4: $\xi $ Hya squared visibilities.

Table 5: $\eta $ Boo squared visibilities.


  \begin{figure}
\par\includegraphics[width=8.3cm,clip]{2075_f01.ps}
\end{figure} Figure 1: Squared visibility measurements obtained on $ \rm \delta ~ Eri$. The solid line is a limb darkened disk model with $ \rm \theta_{LD} = 2.394 \pm 0.014 \pm 0.025 ~ mas$ (statistical and systematic errors).
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{2075_f02.ps}
\end{figure} Figure 2: Squared visibility measurements obtained on $ \rm \xi ~Hya$. The solid line is a limb darkened disk model with $ \rm \theta_{LD} = 2.386 \pm 0.009 \pm 0.019 ~ mas$ (statistical and systematic errors).
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{2075_f03.ps}
\end{figure} Figure 3: Squared visibility measurements obtained on $\rm \eta ~ Boo$. The solid line is a limb darkened disk model with $ \rm \theta_{LD} = 2.200 \pm 0.027 \pm 0.016 ~ mas$ (statistical and systematic errors).
Open with DEXTER

The VINCI instrument has no spectral dispersion and its bandpass corresponds to the K band filter (2-2.4 $\mu$m). It is thus important to compute the precise effective wavelength of the instrument in order to determine the angular resolution at which we are observing the targets. The effective wavelength differs from the filter mean wavelength because of the detector quantum efficiency curve, the fiber beam combiner transmission and the object spectrum. It is only weakly variable as a function of the spectral type.

To derive the effective wavelength of our observations, we computed a model taking into account the star spectrum and the VLTI transmission. The instrumental transmission of VINCI and the VLTI was first modeled taking into account all known effects and then calibrated based on several bright reference star observations with the UTs (see Kervella et al. 2003b, for details).

Taking the weighted average wavelength of this model spectrum gives an effective wavelength of $\rm\lambda_{eff} = 2.178 \pm 0.003 ~ \mu m $ for $\delta $ Eri, $\xi $ Hya and $\eta $ Boo. The visibility fits were computed taking into account the limb darkening of the stellar disk of each stars. We used power law intensity profiles derived from the limb darkening models of Claret (2000) in the K band.

The resulting limb darkened diameters for the three program stars are given in Table 1. The statistical error bars were computed from the statistical dispersion of the series of $\mu^2$ values obtained on each star (typically a few hundred), using the bootstrapping technique. The systematic error bars come from the uncertainties on the angular diameters of the calibrators that were used for the observation. They impact the precision of the interferometric transfer function measurement, and thus affect the final visibility value. Naturally, these calibration error bars do not get smaller when the number of observations increases, as the statistical errors do. The detailed methods and hypothesis used to compute these error bars are given in Kervella et al. (2004b).

3.4 Search for a companion to $\delta $ Eri

$\delta $ Eri is classified as an RS CVn variable (Kholopov et al. 1998), and has shown a small amplitude photometric variability ( mV = 3.51 to 3.56). Fisher et al. (1983) have also reported photometric variations with an amplitude $\Delta m_V = 0.02$ over a period of 10 days. This small amplitude and the apparent absence of periodic radial velocity modulation lead these authors to propose that $\delta $ Eri is a close binary star seen nearly pole on ($i \le 5$ deg). Following this idea, we can suggest three hypotheses to explain the observed photometric variations:
1.
The main star is ellipsoidal. This would result in a modulation of its projected surface along the line of sight during its rotation. This deformation would be caused by the close gravitational interaction of the main star with the unseen companion.
2.
The companion creates a hot spot on the hemisphere of the main star that is facing it. It is changing in apparent surface when the system rotates, probably synchronously.
3.
The pole of the main component has a dark spot that changes in apparent surface during the rotation of the star.
The period of the photometric variations, if attributed to the presence of an orbiting companion, allows us to deduce the distance between the two components through Kepler's third law. At the distance of $\delta $ Eri, this corresponds to an angular separation of approximately 9 mas, easily resolvable using the moderately long baselines of the VLTI. Using the B3-D1 stations of the VLTI, we have taken advantage of the fact that the azimuth of the projected baseline is almost constant for observations of $\delta $ Eri to monitor the evolution of its visibility over a period of 13 nights. The projected length is also very well suited to the expected separation. Our interferometric data (Fig. 4) does not show any systematic deviation from the uniform disk model fit obtained using the longer baselines, at a level of $0.2 \pm 0.3$%, consistent with zero. From these measurements, we conclude that no companion is detected at a level of about $\pm 2$% of the luminosity of the primary star. This result is consistent with the fact that $\delta $ Eri does not deviate significantly from the surface-brightness relations determined by Kervella et al. (2004c).


  \begin{figure}
\par\includegraphics[width=7.8cm,clip]{2075_f04.ps}
\end{figure} Figure 4: Observed deviation of the squared visibilities of $ \rm \delta ~ Eri$ (B3-D1 baseline only) with respect to the visibility model of a $ \rm \theta_{UD} = 2.394 ~ mas $ uniform disk model. The dashed line represents the average deviation over all observations (0.21%).
Open with DEXTER

   
4 Models and results

In order to obtain a rapid estimate of the improvements brought by the new interferometric constraints on the radius on the determination of the mass and age of the three stars, we have calculated evolutionary stellar models that we compare to observations. In these models we have adopted a given set of standard input physics and the observational parameters described in Sect. 2 and Table 1. We do not intend to examine in detail the effects on the uncertainties of the details of the models (envelope, convection, overshooting or other extra mixing) on the results presented here.

The parameters used to construct our CESAM (Morel 1997) evolutionary models are summarized in Table 1. The convection is described by Canuto & Mazitelli's theory (1991, 1992) and the atmospheres are restored on the basis of Kurucz's atlas models (1992). The other input physics are identical to those adopted for the star Procyon (see Kervella et al. 2004a). The adopted metallicity Z/X, which is an input parameter for the evolutionary computations, is given by the iron abundance measured in the atmosphere with the help of the following approximation: $ \log\left(\frac ZX\right) \simeq [{\rm Fe/H]} + \log\left(\frac ZX\right)_\odot $. We use the solar mixture of Grevesse & Noels (1993): $\left(\frac Z X\right)_\odot = 0.0245 $.

The evolutionary tracks are initialized at the Pre-Main Sequence stage. Note that the age is counted from the ZAMS. In CESAM, the ZAMS is defined as the stage of the end of the Pre-Main Sequence where the gravitational energy release is equal to the nuclear one. We have computed models with and without microscopic diffusion of chemical species.

To fit observational data (effective temperature $ T_{\rm eff}$ , luminosity L and surface metallicity $ [Z/X]_{\rm surf}$) with corresponding results of various computations, we adjust the main stellar modeling parameters: mass, age and metallicity. In figures representing the zoom of HR diagram (Figs. 68, 10 and 12), the (rectangular) error boxes are derived from the values and accuracies of the stellar parameters quoted in Table 1. The present (new) values of radii, presented in this paper, select sub-areas in these error boxes and hence the new measures of diameters are used to discriminate between our models (see Table 1). Our best model is the one that satisfies first the luminosity and radius constraint and second the effective temperature constraint. On the zooms of the HR diagrams (see Figs. 6, 8, 10 and 12), the measured radius and its confidence interval appear as diagonal lines. We notice that the addition of the radius measurement reduces significantly the uncertainty domain, and in some cases tightens the allowed range for ages by a factor of three (see below). We have computed models that include overshooting of the convective core (radius $ R_{\rm co}$) over the distance $ O_{\rm v} = A_{\rm ov} \min(H_{\rm p}, R_{\rm co})$ where $ R_{\rm co}$ is the core radius, following the prescriptions of Schaller et al. (1992).

 

4.1 $\delta $ Eri

First, we adopt an initial helium content similar to the Sun, $ Y_{\rm ini} = 0.28 $ and $ [Z/X]_{\rm ini} = 0.148 $, both stars having similar ages and abundances (this will be confirmed hereafter).

Then, with mass and metallicity as free parameters, we have computed a grid of evolutionary tracks in order to reproduce observational data. Our best model without diffusion and without overshooting gives $ M = 1.215 ~ M_\odot $ and an age (from the ZAMS) of $ \rm 6196 ~ Myr$. Our best model with diffusion and an overshooting value of $ A_{\rm ov} = 0.15 $ in agreement with the results of Ribas et al. (2000) gives $ M = 1.215 ~ M_\odot $, an age (from the ZAMS) of $ \rm 6194 ~ Myr $ and a diameter of $ D = 2.328 ~ D_\odot $. See Figs. 5 and 6.

The mean large frequency splitting found for our best model is $ \rm 45.27 ~ \mu Hz$. This result is in agreeement within two per cent with the value of $ \rm 43.8~ \mu Hz $ of the mean large frequency splitting reported by Carrier et al. (2003).

4.2 $\xi $ Hya

We have computed a grid of evolutionary tracks (with and without diffusion) in order to reproduce observational data. Hence, we derived the following parameters: $ M = 2.65 ~ M_\odot $, $ Y_{\rm ini} = 0.275 $ and $ [Z/X]_{\rm ini} \equiv 0.0 $. Our best model with diffusion and an overshooting value of $ A_{\rm ov} = 0.20 $ in agreement with the results of Ribas et al. (2000) gives us an age (from the ZAMS) of $ \rm 509.5 ~ Myr $ and a diameter of $ D = 10.3 ~ D_\odot $. To improve the modeling, a better precision of the diameter is required as it is the case for the two other stars discussed in this paper, for which the accuracy is better by an order of magnitude. See Figs. 7 and 8.


  \begin{figure}
\par\includegraphics[angle=-90, width=8.5cm,clip]{2075_f05.ps}
\end{figure} Figure 5: Evolutionary tracks in the H-R diagram for $\delta $ Eri from label "A'' (0 Myr) to label "G'' (6000 Myr), shown by upper case letters and squares with time steps of 1000 Myr; from label "h'' (6100 Myr) to label "j'' (6300 Myr), shown by lower case letters and triangles with time steps of 100 Myr.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90, width=8.5cm,clip]{2075_f06.ps}
\end{figure} Figure 6: Zoom of the evolutionary tracks in the H-R diagram for $\delta $ Eri from label "a'' (6140 Myr) to label "j'' (6230 Myr), shown by lower case letters and triangles with time steps of 10 Myr (except label "G'' at 6200 Myr shown by an upper case letter and a square). Our best model is close to label "f'' at 6194 Myr (see Table 1).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90, width=8.5cm,clip]{2075_f07.ps}
\end{figure} Figure 7: Evolutionary tracks in the H-R diagram for $\xi $ Hya from label "A'' (0 Myr) to "F'' (500 Myr), shown by upper case letters and squares with time steps of 100 Myr; from label "g'' (502 Myr) to "p'' (511 Myr), show by lower case letters and triangles with time steps of 1 Myr.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90, width=8.5cm,clip]{2075_f08.ps}
\end{figure} Figure 8: Zoom of the evolutionary tracks in the H-R diagram for $\xi $ Hya. from label "a'' (509.2 Myr) to "g'' (509.8 Myr), shown by lower case letters and triangles with time steps of 0.1 Myr. Our best model is close to label "d'' at 509.5 Myr (see Table 1).
Open with DEXTER

Solar-like oscillations of that star were discovered by Frandsen et al. (2002) with a mean spacing of $ \rm 7.1~ \mu Hz $; see also Teixeira et al. (2003). From our model, we computed a value of $ \rm 7.2 ~ \mu Hz $ similar to the theoretical value presented by Frandsen et al. or Teixeira et al.

4.3 $\eta $ Boo

Concerning the values of $ T_{\rm eff}$ and its corresponding uncertainty, we have chosen conservative values based upon various determinations: Feltzing & Gonzales (2001) gives $ T_{\rm eff} = 6000 \pm 100 ~ \rm K $ whereas Cayrel de Strobel (2001) gives a range between 5943 and 6219 K . We notice that DiMauro et al. adopt $ T_{\rm eff} = 6028 \pm 45 ~ \rm K $ but in our study, we take advantage of the constraint given by the new diameter value which reduces the uncertainty as shown in Figs. 10 or 12.

In a first attempt to characterize this star, DiMauro et al. (2003) limit the range of mass between $ 1.64 ~ M_\odot $ and $ 1.75 ~ M_\odot $. Guenther (2004) adopted in his conclusion a mass of $ 1.706 ~ M_\odot $ with an initial chemical composition: $ X_{\rm ini} = 0.71 $ , $ Y_{\rm ini} = 0.25 $ and $ Z_{\rm ini} = 0.04 $. In the present study, we have computed a grid of models and it appears that the best fitting parameters are $ M = 1.70 ~ M_\odot $ with an initial chemical composition $ X_{\rm ini} = 0.70 $, $ Y_{\rm ini} = 0.26 $ and $ \rm Z_{ini} = 0.04 $. A first set of models have been computed with the simplest available reliable physics (and therefore without diffusion, as probably done by the previously cited authors). A second set of models have also been computed with improved physics. Thus, we include convective overshooting (with $ A_{\rm ov} = 0.15 $, see previous discussion), diffusion and radiative diffusivity (see Morel & Thévenin 2002) which controls diffusion of chemical elements in intermediate mass stars. The two sets of results give similar results except for the ages: the age of the best model with diffusion (2738.5 Myr) is higher than the age of the best model without diffusion (2355.0 Myr).

As shown, for example, in Fig. 10, without the constraint given by the diameter, the age would range from 2295 Myr (between label "b'' and label "c'') to 2410 Myr (close to label "n''), with a derived uncertainty of 115 Myr. For a given set of input physics, the constraint on diameter reduces the uncertainty on the age by about a factor of three: the age would be ranging from 2323 Myr (close to label "e'') to 2370 Myr (close to label "j''), corresponding to a (reduced) uncertainty of 47 Myr (Figs. 9-12). Note that our model for $\eta $ Boo with diffusion (Figs. 11 and 12) has the star in a very short-lived phase of evolution (which is, of course, possible but with a small, but non zero, probability).

  \begin{figure}
\par\includegraphics[angle=-90, width=8.5cm,clip]{2075_f09.ps}
\end{figure} Figure 9: Evolutionary tracks in the H-R diagram for $\eta $ Boo (model without diffusion) from label "A'' (0 Myr) to "E'' (2000 Myr), shown by upper case letters and squares with time steps of 500 Myr; from label "f'' (2200 Myr) to "o'' (2650 Myr), shown by lower case letters and triangles with time steps of 50 Myr (except label "L'' at 2500 Myr shown by an upper case letter and a square).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90, width=8.5cm,clip]{2075_f10.ps}
\end{figure} Figure 10: Zoom of the evolutionary tracks in the H-R diagram for $\eta $ Boo (model without diffusion) from label "a'' (2280 Myr) to "n'' (2410 Myr), shown by lower case letters and triangles with time steps of 10 Myr (except labels "C'' at 2300 Myr and label "M'' a 2400 Myr shown by an upper case letters and squares). Our best model is close to label "h'' at 2350 Myr (see Table 1).
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90, width=8.3cm,clip]{2075_f11.ps}
\end{figure} Figure 11: Evolutionary tracks in the H-R diagram for $\eta $ Boo (model with diffusion) from label "A'' (0 Myr) to label "F'' (2500 Myr), shown by upper case letters and squares with a time step of 500 Myr; label "g'' at 2700 Myr shown by a triangle; from label "h'' (2750 Myr) to label "k'' (2900 Myr), shown by lower case letters and triangles with a time step of 50 Myr.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=-90, width=8.5cm,clip]{2075_f12.ps}
\end{figure} Figure 12: Zoom of the evolutionary tracks in the H-R diagram for $\eta $ Boo (model with diffusion) from label "a'' (2738 Myr) to label "g'' (2744 Myr), shown by lower case letters and triangles with time steps of 1 Myr (except label "C'' at 2740 Myr shown by a square); from label "H'' (2745 Myr) to label "K'' (2760 Myr), shown by upper case letters and squares with time steps of 5 Myr. Our best model is between label "a'' (at 2738 Myr) and label "b'' (at 2739 Myr) (see Table 1).
Open with DEXTER

5 Concluding remarks

We have measured with the instrument VLTI/VINCI the angular diameters of three subgiant and giant stars and used them as an additive constraint to the spectro-photometric and asteroseismic ones to perform a study of the stellar evolutionary status.

Owing the position of the three stars in the HR diagram, the determination of the modeling parameters, in particulary the age, is very sensitive to the input physics, due to the rapidity of the stellar evolution compared to the size of the error boxes.

With our input physics and observational constraints, $\delta $ Eri is a star at the end of the subgiant phase ( $ M = 1.215 ~ M_\odot $) with an age of 6.2 Gyr. We attempt without success to detect a close companion forcing us to conclude that the classification of $\delta $ Eri as an RS CVn star is doubtful.

$\xi $ Hya has been constrained with success with a model adopting a mass of $ 2.65 ~ M_\odot $ and an age of 510 Myr.

$\eta $ Boo is a subgiant slightly more evolved than Procyon with a similar age of 2.7 Gyr. With a mass of at $ M = 1.7 ~ M_\odot $ (similar to the mass adopted by Di Mauro et al. 2003), we were able to reproduce the VLTI/VINCI radius. We notice that because of the short evolutionary time scales of a model crossing rather large error boxes, the results of the models - in particular the age - are very sensitive to the input physics (for instance, the core mixing. Some progress in the asteroseismic observations is now required to better constrain the evolutionary state of giant stars for which the frequency spacings (Bouchy & Carrier 2003; Bedding & Kjeldsen 2003) are still relatively imprecise. The improvement of the angular diameter estimations in the future will further tighten the uncertainty domain in the HR diagram, especially as detailed modeling of the atmosphere will be required. This improvement will naturally require a higher precision on the parallax value to derive the linear diameters.

Acknowledgements
The VINCI public commissioning data reported in this paper has been retrieved from the ESO/ST-ECF Archive. The VINCI pipeline includes the wavelets processing technique, developed by D. Ségransan (Obs. de Genève). No VLTI observation would have been possible without the efforts of the ESO VLTI team, to whom we are grateful. This work has been performed using the computing facilities provided by the program Simulations Interactives et Visualisation en Astronomie et Mécanique (SIVAM) at the computer center of the Observatoire de la Côte d'Azur. This research has made use of the Simbad database operated at CDS, Strasbourg, France. We thank the referee, T. R. Bedding, for his suggested improvements of this paper.

References

 

Copyright ESO 2005