A&A 436, 253-262 (2005)
DOI: 10.1051/0004-6361:20042075
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
Eri,
Hya,
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
Eri, we doubt its
classification as an RS CVn star.
Key words: stars: evolution - stars: fundamental parameters - techniques: interferometric
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 "
'' and "
'' 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
is equal to
.
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
(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
0.02 mag
with a period of 10 days.
They suggested that
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
Eri making doubtful its classification as an RS CVn star.
Eri has a projected rotational velocity of
(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
Eri,
we set several VLTI/VINCI observations at different baselines (see Sect. 3).
We estimate its bolometric luminosity
using the Alonso et al. (1999)
empirical bolometric corrections (BC,
for giants,
this is the dominant source of uncertainty on the luminosity).
We adopt the Santos et al. (2004) values for the effective temperature
,
logarithmic surface gravity
and surface
iron abundance
.
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
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:
.
In the high frequency range, i.e. large radial orders,
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.
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
.
We estimate its bolometric luminosity
using BC (
)
from Alonso et al. (1999).
We adopt the spectroscopic parameters derived by Mc William (1990):
effective temperature
,
and
.
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
.
Using a set of CORALIE spectra, Frandsen et al. (2002)
detected solar-like oscillations suggesting radial modes with the largest amplitudes
almost equidistant around
.
That important detection opens the possibility of better
constraining the model of this star for which the mass is not well-known.
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
.
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
using BC
(
,
this is the dominant source of
uncertainty on the luminosity) from Vandenberg & Clem (2003) for this subgiant,
with an effective temperature
representing the average of five effective temperature determinations in the [Fe/H]
catalogue of Cayrel de Strobel et al. (2001) and the spectroscopic
and
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
Eri have been reported
by Carrier et al. (2005) with
and by Kjeldsen et al. (2003) with
.
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 (
)
for these observations.
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 ()
measurement, and avoids biases from corrupted interferograms.
The final
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
are computed from the series of
values obtained on each target star (typically a few hundred interferograms)
using the bootstrapping technique.
The calibration of the visibilities obtained on Eri and
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:
Eri squared visibilities.
Table 3:
Eri squared visibilities
(continued from Table 2).
Table 4:
Hya squared visibilities.
Table 5:
Boo squared visibilities.
![]() |
Figure 1:
Squared visibility measurements obtained on
![]() ![]() |
Open with DEXTER |
![]() |
Figure 2:
Squared visibility measurements obtained on
![]() ![]() |
Open with DEXTER |
![]() |
Figure 3:
Squared visibility measurements obtained on
![]() ![]() |
Open with DEXTER |
The VINCI instrument has no spectral dispersion and its bandpass
corresponds to the K band filter (2-2.4 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
for
Eri,
Hya and
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
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).
![]() |
Figure 4:
Observed deviation of the squared visibilities of
![]() ![]() |
Open with DEXTER |
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:
.
We use the solar mixture of Grevesse & Noels (1993):
.
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
,
luminosity L and surface metallicity
)
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. 6, 8, 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
)
over the distance
where
is the core radius,
following the prescriptions of Schaller et al. (1992).
First, we adopt an initial helium content similar to the Sun,
and
,
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
and an age (from the ZAMS) of
.
Our best model with diffusion and an overshooting value of
in agreement with the results of Ribas et al. (2000) gives
,
an age (from the ZAMS) of
and a diameter of
.
See Figs. 5 and 6.
The mean large frequency splitting found for our best model is
.
This result is in agreeement within two per cent with the value of
of the mean large frequency splitting reported by Carrier et al. (2003).
We have computed a grid of evolutionary tracks (with and without diffusion)
in order to reproduce observational data.
Hence, we derived the following parameters:
,
and
.
Our best model with diffusion and an overshooting value of
in agreement with the results of Ribas et al. (2000) gives us
an age (from the ZAMS) of
and a diameter of
.
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.
![]() |
Figure 5:
Evolutionary tracks in the H-R diagram for ![]() |
Open with DEXTER |
![]() |
Figure 6:
Zoom of the evolutionary tracks in the H-R diagram for ![]() |
Open with DEXTER |
![]() |
Figure 7:
Evolutionary tracks in the H-R diagram for ![]() |
Open with DEXTER |
![]() |
Figure 8:
Zoom of the evolutionary tracks in the H-R diagram for ![]() |
Open with DEXTER |
Solar-like oscillations of that star were discovered by Frandsen et al. (2002)
with a mean spacing of
;
see also Teixeira et al. (2003).
From our model, we computed a value of
similar to the theoretical value presented by
Frandsen et al. or Teixeira et al.
Concerning the values of
and its corresponding uncertainty, we
have chosen conservative values based upon various
determinations:
Feltzing & Gonzales (2001) gives
whereas
Cayrel de Strobel (2001) gives a range between 5943 and 6219 K .
We notice that DiMauro et al. adopt
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
and
.
Guenther (2004) adopted in his conclusion a mass of
with an initial chemical composition:
,
and
.
In the present study, we have computed a grid of models and it appears that the
best fitting parameters are
with an initial chemical composition
,
and
.
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
,
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
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).
![]() |
Figure 9:
Evolutionary tracks in the H-R diagram for ![]() |
Open with DEXTER |
![]() |
Figure 10:
Zoom of the evolutionary tracks in the H-R diagram for ![]() |
Open with DEXTER |
![]() |
Figure 11:
Evolutionary tracks in the H-R diagram for ![]() |
Open with DEXTER |
![]() |
Figure 12:
Zoom of the evolutionary tracks in the H-R diagram for ![]() |
Open with DEXTER |
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,
Eri is a star at the end of the subgiant phase
(
)
with an age of 6.2 Gyr.
We attempt without success to detect a close companion forcing us to conclude
that the classification of
Eri as an RS CVn star is doubtful.
Hya has been constrained with success with a model adopting a mass of
and an age of 510 Myr.
Boo is a subgiant slightly more evolved than Procyon with a similar
age of 2.7 Gyr.
With a mass of at
(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.