A&A 408, 681-688 (2003)
DOI: 10.1051/0004-6361:20030994
P. Kervella1 - F. Thévenin2 - P. Morel2 - P. Bordé3 - E. Di Folco4
1 - European Southern Observatory, Alonso de Cordova 3107, Casilla 19001, Vitacura, Santiago 19, Chile
2 - Département Cassini, UMR CNRS 6529, Observatoire de la Côte
d'Azur, BP 4229, 06304 Nice Cedex 4, France
3 - LESIA, Observatoire de Paris-Meudon, 5 place Jules Janssen, 92195 Meudon Cedex, France
4 - European Southern Observatory, Karl-Schwarzschild-str. 2, 85748 Garching, Germany
Received 10 April 2003 / Accepted 23 June 2003
Abstract
The interferometric observations of dwarf stars in the solar neighbourhood,
combined with HIPPARCOS parallaxes provide very precise values
of their linear diameters. In this paper, we report the direct measurement of the
angular diameter of the bright star Sirius A with the VINCI/VLTI instrument.
We obtain a uniform disk angular diameter of
mas in the K band
and a limb darkened value of
mas.
In combination with the HIPPARCOS parallax of
mas,
this translates into a linear diameter of
.
Using the VINCI/VLTI interferometric diameter and the published properties of
Sirius A, we derive internal structure models corresponding to
ages between 200 and
Myr. This range is defined mainly
by the hypothesis on the mass of the star, the overshoot and the metallicity.
Key words: stars: individual: Sirius - stars: fundamental parameters - stars: evolution - techniques: interferometric
The brightest star in the sky, Sirius A (HR 2491, HD 48915), is an A1V type dwarf. It is also a member of the fifth nearest binary system. The existence of a massive companion, Sirius B, was predicted by Bessel (1844) based on the apparent motion of Sirius A. It was later discovered visually by A. G. Clark in 1862, and is still today among the most massive white dwarfs (WD) known. Chandrasekhar (1935) published the first model to explain its characteristics. Sirius has been reported several times to be a red star in the ancient times, whereas it is now almost white. The reason behind this color change is still not understood, though several scenarios have been proposed in the last decades (Schlosser & Bergmann 1985; Bonnet-Bideau & Gry 1991; Whittet 1999). During the past century, many efforts were directed towards compiling the properties of both stars like their masses or chemical abundances. Sirius A exhibits enhanced metal abundances in its photosphere and has therefore been classified as Am by Abt (1979). As the orbital period is relatively short (50 years), two revolutions of the Sirius A-B binary have already been measured, allowing Gatewood & Gatewood (1978) to derive precisely the masses of the two components (Sect. 2.1). Moreover, Sirius A has been the first main sequence star to have its diameter measured interferometrically by Hanbury Brown et al. (1967, 1974a).
The chief object in our study is to model the internal structure of Sirius A in order to reproduce its macroscopic characteristics. In this process, we use the angular diameter of this star from our new VINCI/VLTI measurements and the HIPPARCOS parallax to derive its linear diameter. From our models, we compute the age and the initial helium content of Sirius A. We summarize in Sect. 2 the fundamental parameters of this star, while Sects. 3 to 5 are dedicated to the VINCI data processing and analysis. This section includes in particular the fit of the limb darkened disk model to the measured visibilities. Finally, we discuss in Sect. 6 the macroscopic constraints that we apply to our CESAM numerical models (Morel 1997) of the internal structure of Sirius A, and we predict the asteroseismic large frequency spacing of this star.
Table 1:
Relevant parameters of CMa and its calibrators from the literature.
The age of the Sirius system has been poorly debated in the literature.
Based on Wood's (1995) sequences of cooling age of WD
(carbon-oxygen model with thick H layer),
Sirius B is found to be a young WD of about 160 (Holberg et al. 1998) to
Myr (Holberg, private communication).
The revisited mass of Sirius B (1.034
)
corresponds
to a progenitor mass of
,
i.e.
the approximate mass of a B5 main sequence star,
from the empirical initial to final mass relation
of Weidemann (2002).
Such a massive star evolves to a WD in
Myr.
Consequently, assuming a simultaneous origin of both A & B components,
the Sirius system is expected to be between 200 and 250 Myr old.
Sirius shows deep metallic spectral lines, and we propose that the corresponding enhancement
in metal abundance is a consequence of the strong radiative forces at work in the interior of this star.
This hypothesis is supported by the very low rotational velocity (Royer et al. 2002)
km s-1, that favors the action of these forces.
In their comparison with the star Vega, Qiu et al. (2001) concluded that
Sirius A appears to be enhanced in abundances by about 1.0 dex on
average while C, O, Sc and Ca are underabundant.
In consequence, for such a slow rotator and intermediate mass,
it is very hazardous to estimate the mean metal content in the interior of the
star from its surface metal abundance.
Therefore, we propose to focus only on the abundances of CNO,
that are believed to be little modified by radiative forces (Richer et al. 2000),
to estimate the value of
at present day.
We adopt an average abundance at present day of
in our models, that includes the
microscopic diffusion of the elements.
is very important for the internal stellar opacity of the
computed model and is a key parameter with respect to the age
of the star (see Sect. 6).
It results, when using these elements as tracers, that the average metal
abundance of Sirius A is in fact about half that of the Sun.
While appearing to be low, the adopted
abundance is comparable with Vega's observed value. Vega is a fast
rotator and therefore should not be affected by radiative accelerations
masking its true metal abundance.
This corresponds to an initial value at zero age of
.
In order to validate the previous assumption that the metallicity of outerlayers
does not represent the metallicity of the whole star, we also present, in Sect. 6,
a model of Sirius A with the observed metallicity
,
corresponding to the observed mean surface abundances.
The photometric properties of Sirius are summarized in Table 1.
The magnitude and the parallax give an intrinsic luminosity of
.
Following Lemke (1989),
we adopt an effective temperature of
K and
a logarithm of the spectroscopic surface gravity of
.
The European Southern Observatory's Very Large Telescope Interferometer
(Glindemann et al. 2000) is operated on top of the Cerro Paranal, in Northern Chile since
March 2001. For the observations reported in this paper, the light coming from two test siderostats
(0.35 m aperture) was recombined coherently in VINCI, the VLT INterferometer Commissioning Instrument
(Kervella et al. 2000). We used a regular K band filter (
m)
but VINCI can also be operated in the H band (
m) using an
integrated optics beam combiner (Berger et al. 2001; Kervella et al. 2003a).
Two VLTI baselines were used for this program: D1-B3 and E0-G1, respectively 24 and 66 m
in ground length.
Table 2: Sirius squared visibilities.
We used a modified version of the standard VINCI data reduction pipeline (Kervella et al. 2003c), 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 the continuous wavelet transform (Farge 1992). In this approach, the projection of the signal is not done onto a sine wave (Fourier transform), but on a function, i.e. the wavelet, that is localized in both time and frequency. We used as a basis the Morlet wavelet, a Gaussian envelope multiplied by a sine wave. With the proper choice of the number of oscillations inside the Gaussian envelope, this wavelet closely matches a VINCI interferogram. It is therefore very efficient at localizing the signal in both time and frequency.
The differential piston corrupts the amplitude and the shape of the
fringe peak in the wavelets power spectrum.
A selection based on the shape of fringe peak in the
time-frequency domain is used to remove "pistonned'' and false detection interferograms.
Squared coherence factors
are then derived by integrating the
wavelet power spectral density (PSD)
of the interferograms at the position and frequency of the fringes.
The residual photon and detector noise are removed
by making a least squares fit of the PSD at high and low frequency.
The resulting measurement stability is satisfactory: on a good series of 500 interferograms of Sirius
(10 min) the standard deviation is typically 2.0%, and the final bootstrapped
statistical error on the average
is only 0.16%.
The
values are converted into calibrated visibilities V2 through
the observation of calibrator stars (Sect. 4).
The calibration of Sirius visibilities was achieved using well-known calibrator stars that have been selected in the Cohen et al. (1999) catalogue. The characteristics of the selected calibrators are listed in Table 1. The limb-darkened disk (LD) angular diameter of these stars was converted into a uniform disk value using linear coefficients taken from Claret et al. (1995). As demonstrated by Bordé et al. (2002), the star diameters in the Cohen et al. (1999) list have been measured very homogeneously to a relative precision of approximately 1% and agree well with other angular diameter estimation methods. The interferometric efficiency is computed from the coherence factors obtained on these calibrators, taking into account the bandwidth smearing effect (Sect. 5.2) and a uniform disk angular diameter model. This process yields the calibrated squared visibilities that are used for the model fit.
In Table 2, two error bars are given for each V2 value:
From the visibility curve shape before the first minimum, it is almost impossible
to distinguish between a uniform disk (UD) and limb darkened (LD) model. Therefore, it
is necessary to use a model of the stellar disk limb darkening to deduce the photospheric angular
size of the star, from the observed visibility values.
The intensity profiles that we used were computed by Claret (2000), based on
model atmospheres by Kurucz (1992).
We chose the four parameters
law of this author, where
is the cosine of the azimuth of a surface element of the star. They are
accurate approximations of the numerical results from the ATLAS modeling code.
To compare the different types of approximated limb darkening laws,
we chose the following approximate parameters for Sirius:
K (for Claret et al. 1995 and Claret 2000)
or
K (for Claret 1998),
,
,
km s-1.
These were chosen as close as possible to the Cohen et al. (1992) values.
As noted by Claret (2000), the sensitivity of the limb
darkening of a
K star to the metallicity is negligible,
so the true metal content of Sirius is not critical.
Different
square root and four parameters laws for Sirius, from Claret et al. (1995)
and Claret (1998, 2000) are plotted in Fig. 3.
Though no formal error bars are provided for these limb darkening models, the small differences between
the four curves demonstrates the good internal consistency of the different types of laws.
The limb darkening is directly measurable by interferometry beyond the first minimum of the visibility function, as demonstrated by several authors on giant stars (Quirrenbach et al. 1996; Wittkowski et al. 2001). Unfortunately, for stars of the angular size of Sirius observed in the K band, this requires baseline of 80 to 120 m that were not available for the measurements reported here. It is intended in the near future to measure directly the LD of a number of nearby stars, using the VINCI and AMBER (Petrov et al. 2000) instruments on the long baselines of the VLTI. With lengths of up to 202 m, they will allow the exploration of the secondary and higher order lobes of the visibility function.
The VINCI instrument bandpass corresponds to the K band filter (2-2.4 m).
An important effect of this relatively large spectral bandwidth
is that several spatial frequencies are simultaneously
observed by the interferometer. This effect is known as the bandwidth smearing.
It is usually negligible for
%, but this is not the case when visibilities
come closer to the first minimum of the visibility function.
In the case of Sirius A observed with the E0-G1 baseline, V2 is already low (13%) and we have to consider it.
To account for the bandwidth smearing, the model visibility is computed for regularly spaced
wavenumber spectral bins over the K band, and then integrated to obtain the model visibility.
We assume in this paper a limb darkening model following the four parameters
law of Claret (2000):
![]() |
(1) |
![]() |
(2) |
![]() |
(3) |
![]() |
(4) |
![]() |
(5) |
![]() |
(6) |
The
model is adjusted numerically to the observed (B, V2) data using a classical
minimization process to derive
.
The reduced
of the global fit is 0.30, characteristic of a satisfactory correspondence between
the model and the measured V2 values. Figures 1 and 2
show the position of the VINCI points on the squared visibility curve of the chosen model.
![]() |
Figure 1: Overview of the VINCI/VLTI visibility measurements obtained on Sirius. The solid line is the best fit visibility model, that takes into account the limb darkening and the bandwidth smearing. |
Open with DEXTER |
![]() |
Figure 2:
Visibility measurements obtained on Sirius on the E0-G1 baseline. The thin line is the best fit model,
and the dashed lines mark the limits of the 1![]() |
Open with DEXTER |
![]() |
Figure 3: Published LD laws for Sirius in the K band, from Claret et al. (1995), Claret (1998) and Claret (2000). The square root law is already a good approximation to the Kurucz atmosphere models. The solid line shows results with the most recent four parameters law introduced by Claret (2000). This last version was chosen for the Sirius LD fit of this paper. |
Open with DEXTER |
The
minimization gives a limb darkened angular diameter
mas for Sirius A, while a simpler uniform
disk model yields
mas.
The conversion formula between the LD angular diameter and the
linear diameter D (in solar unit) is:
Our LD angular diameter value can be compared to the value published by
Hanbury Brown et al. (1974a).
They found, using the Narrabri intensity interferometer,
mas.
The difference with our value of
mas
is
,
making the two results statistically compatible within their error bars.
Davis & Tango (1986) have obtained a
value of
mas using the amplitude
interferometer, and revised the Hanbury Brown et al. (1974a)
value to
mas. Using a limb darkening
coefficient
u = 0.5900 (B band from Claret et al. 2000) and
the conversion factor formula from Hanbury Brown et al. (1974b):
![]() |
(8) |
Using spectro-photometric observations, Cohen et al. (1992) have derived the angular diameter of Sirius and found 6.04 mas (with Kurucz model atmospheres), in remarkable agreement with our direct measurement.
Due to the relatively slow rotational velocity of Sirius (
km s-1,
from Royer et al. 2002), we do not expect any detectable flattening of its disk.
In addition, the flux of Sirius B is totally negligible compared to A,
in particular in the K band.
We therefore do not foresee an asymmetry of the visibility function in
azimuth due to these two contributors.
Hanbury Brown et al. (1974a) have observed a small asymmetry of the Sirius visibility function, proposing tentatively that this may come from a dust disk around this star. In the error bars, our measurements do not show this asymmetry, as the D1-C3 and E0-G1 points are in good agreement with our single disk visibility model. Nevertheless, it should be noted that we have probed only two azimuth values, namely 80 and 160 degrees (D1-C3 and E0-G1 baselines, respectively), and that our error bars on the diameter measurements on the shorter baseline are relatively large. Therefore, we cannot exclude formally the presence of a disk accounting for approximately 1% of the stellar flux.
As a remark, Benest & Duvent (1995) have proposed that
a late M5 dwarf could be orbiting Sirius A, but our measurements do not have the necessary
sensitivity to record the visibility modulation that would be produced by a star
of this type, with .
We computed a number of evolutionary models of Sirius A using the CESAM
code (Morel 1997) including the pre-main sequence evolution.
The ordinary assumptions of stellar modeling are made, i.e. spherical symmetry,
no rotation, no magnetic field, no mass loss. The relevant nuclear reaction rates
are taken from the NACRE compilation (Angulo et al. 1999). The equation of
state adopted is EFF (Eggleton et al. 1973), and
the OPAL opacities are from Iglesias & Rogers (1996) with the
Grevesse & Noels (1993) mixture.
The microscopic diffusion is described using the formalism of
Burgers (1969) with the resistance coefficients of Paquette et al.
(1986). We take into account the radiative diffusivity as
recommended by Morel & Thévenin (2002), that limits the efficiency
of the microscopic diffusion in outerlayers of stars with intermediate masses.
We have neglected the radiative accelerations, as well as
the changes of abundance ratios between metals within Z.
The atmosphere is restored using Hopf's law (Mihalas 1978).
The adopted radius of the star is the bolometric radius,
where
.
In the convection zones the temperature gradient is
computed according to the
convection
theory with a mixing length parameter of
(Canuto & Mazzitelli 1991, 1992).
Following the prescriptions of Schaller et al. (1992)
we have computed models that include overshooting of the
convective core (radius
)
over the distance
.
We adopted a value A=0.15in agreement with the results by Ribas et al. (2000) for
a 2
star.
In addition, we have also computed stellar models without overshooting (A=0.00)
of the convective core and with A=0.20 to estimate the effect of this
parameter on the age of the star.
As a comparison, we have also built a model with enhanced
,
corresponding to the atmospheric abundances published
by Qiu et al. (2001). It is intended to demonstrate the strong effect of
this overestimated value of the internal abundance on the evolution stage of Sirius.
Each model is described by about 760 mass shells.
We considered that our model was representative of Sirius A when
we reached the interferometric radius and the observed luminosity and effective
temperature of this star (within their respective error bars).
The main characteristics of the computed models are presented in Table 3.
Table 3:
Characteristics of our CMa A models
remaining within the uncertainty box of the
observed properties of this star.
The subscripts i and s refer respectively to
the initial content at zero age and the surface content at the end of
the computed evolution. The adopted luminosity is
with
K.
The observed abundance of oxygen is
dex (see text).
The precision of the age is
12 Myr corresponding to the crossing
time of the shaded parallelogram.
The HR evolution diagram is
presented in Fig. 4.
We adopt the standard definition of the ZAMS, i.e. where nuclear reactions
begin to dominate gravitation as the primary source of energy by at least 50%.
For the mass of Sirius A this corresponds to an elapsed time of 4.2 Myr.
Note that 100% of energy comes from nuclear reactions only after 10 Myr.
The evolution of models starts with the quasi-static contraction
of a cloud with a central temperature of 0.5 MK.
According to these definitions, our model b of Sirius, which fits the VINCI/VLTI radius
and satisfies the adopted luminosity and
,
is aged
of
Myr. The error bar is fixed by the uncertainty on the
VINCI/VLTI radius, that is itself largely dominated by the HIPPARCOS
parallax error bar.
Because the age of the WD can reach 210 Myr (Holberg, private communication)
we also computed models with a lower mass 2.07
which fit both the
observed radius and the
.
We had to change sligthly
and
the value of
within its error bar, see model d,
resulting in an age of 243 Myr for Sirius.
We can produce a greater age of 250 Myr (model e) when increasing the overshooting
parameter A to the value A=0.20.
At the end of the pre-main sequence, Sirius A has generated
a convective core due to the equilibration of the CNO nuclear cycle.
At the expected age of Sirius A for model dthe hydrogen mass fraction at center is
with a temperature of
K. The base of the external convective zone is situated at 0.9926
and
.
Removing the overshooting from our model reduces the age of the binary by 10%
(model a)
while the model c with enhanced abundances leads to a too young age (25 Myr) for
Sirius A, in significant disagreement with the age of Sirius B.
At the present status of our observational constraints on Sirius A, it is impossible
to decide if overshoot exists or not in the central convective core.
Note that our model gives a surface abundance of
dex
which agrees with the abundance of the oxygen found by
Kamp et al. (2002).
Finally, if Sirius A is able to excite radial modes of oscillation,
we predict the mean large frequency spacing
to be ranging
between 81 and 82
Hz
This is the primary observable of asteroseismology, and corresponds to
the difference between the frequencies of oscillation modes with consecutive
radial order n (see Table 3). This large spacing is sensitive
primarily to the stellar density. It could therefore give a direct, high accuracy
estimation of the mass of Sirius A when combined with the
observed radius reported in this paper.
![]() |
Figure 4:
Evolutionary tracks in the HR diagram of models of Sirius A. The dashed rectangle
delimits the uncertainty domain in luminosity and effective temperature, while the shaded area
represents the uncertainty on the interferometric radius. The continuous line corresponds to
model b with overshoot and the dashed-dot line to model a without overshoot,
both with a mass of 2.12 ![]() ![]() |
Open with DEXTER |
The accuracy on this age is greatly strengthened by the VINCI/VLTI radius, thus
encouraging further studies to improve our knowledge of the diameter of nearby stars.
Based on our model, we predict that the asteroseismic
large frequency spacing of Sirius A should be 81 to 82 Hz,
if it exhibits radial oscillations.
A more complete modeling of the external layers of this star for asterosismic
frequencies prediction, including radiative accelerations, is currently
in preparation.
We would like to emphasize here the need for very high accuracy parallax values,
beyond the HIPPARCOS precision. We have shown that the uncertainty on our
determination of the linear radius of Sirius A is largely dominated by the error
on the distance, despite its proximity (2.7 pc).
The next generation of interferometric instruments, such as the AMBER
beam combiner (Petrov et al. 2000), will provide high precision
angular diameters, together with their wavelength dependence.
In order to constrain the stellar structure models, that output linear values,
it will be necessary to match the precision on the distance to the angular diameter
uncertainty. It is expected that future space based astrometric missions
will be able to give a significant improvement over the HIPPARCOS
parallaxes, and are therefore crucial in this respect.
This work on Sirius A and the study of Cen A & B by Kervella
et al. (2003b) have demonstrated that for the nearest stars, thanks to
the very small error bar on the linear radius, we have advantage to replace the
classical error bars in the HR diagram, e.g. luminosity and effective
temperature, by the couple luminosity and radius.
The direct measurement of the limb darkening of Sirius A is the next step of the interferometric study of this star. It will allow a refined modeling of its atmosphere and will complete the calibration of the radius that we presented in this paper. This observation will be achieved soon using the long baselines of the VLTI.
Acknowledgements
These interferometric measurements have been obtained using the VINCI instrument installed at the VLTI. The VLTI is operated by the European Southern Observatory at Cerro Paranal, Chile. This work has made use of the wavelets data processing technique, developed by D. Ségransan (Observatoire de Genève) and integrated in the VINCI pipeline. Observations with the VLTI are only made possible through the efforts of the VLTI team, for which we are grateful. The VINCI public commissioning data reported in this paper have been retrieved from the ESO/ST-ECF Archive (Garching, Germany). This research has made use of the SIMBAD database at CDS, Strasbourg (France).