P. Kervella1 - N. Nardetto2 - D. Bersier3 - D. Mourard2 - V. Coudé du Foresto4
1 - European Southern Observatory, Alonso de Cordova 3107, Casilla 19001, Vitacura, Santiago 19, Chile
2 - Département Fresnel, UMR CNRS 6528, Observatoire de la Côte d'Azur,
BP 4229, 06304 Nice Cedex 4, France
3 - Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA
4 - LESIA, Observatoire de Paris-Meudon, 5 place Jules Janssen, 92195 Meudon Cedex, France
Received 19 June 2003 / Accepted 26 November 2003
Abstract
We report the angular diameter measurements of seven classical Cepheids,
X Sgr, Aql, W Sgr,
Gem,
Dor, Y Oph and
Car
that we have obtained with the VINCI instrument, installed at ESO's VLT Interferometer (VLTI).
We also present reprocessed archive data obtained with the FLUOR/IOTA instrument on
Gem,
in order to improve the phase coverage of our observations.
We obtain average limb darkened angular diameter values of
mas,
mas,
mas,
mas,
mas,
mas, and
mas.
For four of these stars,
Aql, W Sgr,
Dor, and
Car,
we detect the pulsational variation of their angular diameter. This enables
us to compute directly their distances, using a modified version of the
Baade-Wesselink method:
pc,
pc,
pc,
pc.
The stated error bars are statistical in nature.
Applying a hybrid method, that makes use of the Gieren et al. (1998)
Period-Radius relation to estimate the linear diameters, we obtain the following distances
(statistical and systematic error bars are mentioned):
pc,
pc,
pc,
pc,
pc,
pc,
pc.
Key words: techniques: interferometric - stars: variables: Cepheids - stars: oscillations
There are various ways to calibrate the P-L relation. The avenue chosen by the HST Key-Project was to assume a distance to the Large Magellanic Cloud (LMC), thereby adopting a zero point of the distance scale. Freedman et al. (2001) present an extensive discussion of all available LMC distances, and note, with other authors (see e.g. Benedict et al. 2002), that the LMC distance is currently the weak link in the extragalactic distance scale ladder. Another avenue is to determine the zero point of the P-L relation with Galactic Cepheids, using for instance parallax measurements, Cepheids in clusters, or through the Baade-Wesselink (BW) method. We propose in this paper and its sequels (Papers II and III) to improve the calibration of the Cepheid P-R, P-L and surface brightness-color relations through a combination of spectroscopic and interferometric observations of bright Galactic Cepheids.
In the particular case of the P-L relation, the slope ais well known from Magellanic Cloud Cepheids (e.g. Udalski et al. 1999),
though Lanoix et al. (1999) have suggested
that a Malmquist effect (population incompleteness) could bias this value.
On the other hand, the calibration of the zero-point b (the hypothetic absolute
magnitude of a 1-day period Cepheid) requires measurement of the distance to a
number of nearby Cepheids with high precision. For this purpose,
interferometry enables a new version of the Baade-Wesselink method
(BW, Baade 1926; Wesselink 1946)
for which we do not need to measure the star's temperature, as we have directly
access to its angular diameter (Davis 1979; Sasselov & Karovska 1994).
Using this method, we derive directly the distances to the four nearby Cepheids
Aql, W Sgr,
Dor and
Car.
For the remaining three objects of our sample, X Sgr,
Gem and Y Oph,
we apply a hybrid method to derive their distances, based on published values of
their linear diameters.
After a short description of the VINCI/VLTI instrument (Sect. 2),
we describe the sample Cepheids that we selected (Sect. 3).
In Sects. 4 and 5, we report our new observations
as well as reprocessed measurements of Gem
retrieved from the FLUOR/IOTA instrument archive.
Section 6 is dedicated to the computation of the corresponding
angular diameter values, taking into account
the limb darkening and the bandwidth smearing effects.
In Sects. 7 and 8, we investigate
the application of the BW method to our data, and we
derive the Cepheid distances.
We will discuss the consequences of these results for the calibration of the Period-Radius (P-R), Period-Luminosity (P-L) and Barnes-Evans relations of the Cepheids in forthcoming papers (Papers II and III).
The European Southern Observatory's Very Large Telescope Interferometer
(Glindemann et al. 2000) is in operation on Cerro Paranal, in Northern Chile since
March 2001. For the observations reported in this paper, the beams from two Test Siderostats
(0.35 m aperture) or two Unit Telescopes (8 m aperture) were recombined coherently in VINCI,
the VLT INterferometer Commissioning Instrument (Kervella et al. 2000, 2003a).
We used a regular K band filter (
m) that gives an effective observation
wavelength of
m for the effective temperature of typical Cepheids
(see Sect. 6.4 for details).
Three VLTI baselines were used for this program: E0-G1, B3-M0 and
UT1-UT3 respectively 66, 140 and 102.5 m in ground length.
Figure 1 shows their positions on the VLTI platform.
X Sgr | ![]() |
W Sgr | ![]() |
![]() |
Y Oph | ![]() |
|
HD 161592 | HD 187929 | HD 164975 | HD 37350 | HD 52973 | HD 162714 | HD 84810 | |
mVa | 4.581 | 3.942 | 4.700 | 3.731 | 3.928 | 6.164 | 3.771 |
mKb | 2.56 | 1.966 | 2.82 | 1.959 | 2.11 | 2.682 | 1.091 |
Sp. Type | F5-G2II | F6Ib-G4Ib | F4-G2Ib | F4-G4Ia-II | F7Ib-G3Ib | F8Ib-G3Ib | F6Ib-K0Ib |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Min
![]() |
5670 | 5400 | 5355 | 5025 | 5150 | ||
Mean
![]() |
6150 | 5870 | 5769 | 5490 | 5430 | 5300 | 5090 |
Max
![]() |
6820 | 6540 | 6324 | 6090 | 5750 | ||
Min ![]() |
1.86 | 1.25 | 1.72 | 1.60 | |||
Mean
![]() |
2.14 | 1.49 | 1.82 | 1.83 | 1.50 | 1.50 | 1.50 |
Max ![]() |
2.43 | 1.73 | 2.02 | 2.06 | |||
![]() |
0.04 | 0.05 | -0.01 | -0.01 | 0.04 | 0.05 | 0.30 |
T0 (JD-
![]() |
723.9488 | 519.2477 | 726.8098 | 214.2153 | 210.7407 | 715.4809 | 290.4158 |
P (days)g | 7.013059 | 7.176769 | 7.594904 | 9.842425 | 10.150967 | 17.126908 | 35.551341 |
Intensity profilesh | |||||||
a1 | +0.7594 | +0.8816 | +0.8002 | +0.7969 | +0.8713 | +0.8549 | +0.8500 |
a2 | -0.4530 | -0.7418 | -0.5135 | -0.4596 | -0.6536 | -0.5602 | -0.4991 |
a3 | +0.0347 | +0.3984 | +0.1583 | +0.1341 | +0.3283 | +0.2565 | +0.2113 |
a4 | +0.0751 | -0.0778 | +0.0109 | +0.0082 | -0.0610 | -0.0437 | -0.0340 |
|
Despite their brightness, Cepheids are located at large distances, and the H IPPARCOS
satellite (Perryman et al. 1997) could only obtain a limited number of Cepheid
distances with a relatively poor precision.
If we exclude the peculiar first overtone Cepheid UMi (Polaris), the closest Cepheid
is
Cep, located at approximately 250 pc (Mourard et al. 1997;
Nordgren et al. 2000).
As described by Davis (1979) and Sasselov & Karovska (1994),
it is possible to derive directly the distance to the Cepheids for which we can measure the amplitude
of the angular diameter variation. Even for the nearby Cepheids, this requires
an extremely high resolving power, as the largest Cepheid in the sky,
Car,
is only 0.003'' in angular diameter. Long baseline interferometry is therefore
the only technique that allows us to resolve these objects.
As a remark, the medium to long period Cepheids (
)
in the Large Magellanic
Cloud (LMC) (
kpc) are so small (
as)
that they would require a baseline of 20 km to be resolved in the K band (5 km in the visible).
However, such a measurement is highly desirable, as it would provide a
precise geometrical distance to the LMC, a critical step in the
extragalactic distance ladder.
![]() |
Figure 1: Layout of the three baselines used for the VINCI/VLTI Cepheids observations, UT1-UT3 (102.5 m), E0-G1 (66 m) and B3-M0 (140 m). |
Mourard (1996) has highlighted the capabilities of the VLTI for
the observation of nearby Cepheids, as it provides long baselines (up to 202 m) and thus
a high resolving power.
Though they are supergiant stars, the Cepheids are
very small objects in terms of angular size.
A consequence of this is that the limit on the number of interferometrically
resolvable Cepheids is not set by the size of the light collectors, but
by the baseline length.
From photometry only, several hundred Cepheids can produce interferometric
fringes using the VLTI Auxiliary Telescopes (1.8 m in diameter). However,
in order to measure accurately their size, one needs to resolve their disk
to a sufficient level.
Kervella (2001a) has compiled a list of
more than 30 Cepheids that can be measured from Paranal using
the VINCI and AMBER (Petrov et al. 2000) instruments.
Considering the usual constraints in terms of sky coverage, limiting magnitude
and accessible resolution, we have selected seven bright Cepheids observable
from Paranal Observatory (latitude
):
X Sgr,
Aql, W Sgr,
Dor,
Gem, Y Oph and
Car.
The periods of these stars cover a wide range, from 7 to 35.5 days. This
coverage is important to properly constrain the P-R and P-L relations.
To estimate the feasibility of the observations, the angular diameters of these
stars were deduced from the BW studies by Gieren et al. ().
For
Gem and
Aql, previously published direct interferometric
measurements by Nordgren et al. (2000), Kervella et al. (2001b)
and Lane et al. (2002) already demonstrated the feasibility of the observations.
The relevant parameters of the seven Cepheids of our sample, taken from the literature,
are listed in Table 1.
We used a modified version (Kervella et al. 2003c) of the standard VINCI data reduction pipeline, whose general principle is based on the original algorithm of the FLUOR instrument (Coudé du Foresto et al. 1997, 1998a). The VINCI/VLTI commissioning data we used for this study are publicly available through the ESO Archive, and result from two proposals of our group, that were accepted for ESO Periods 70 and 71.
The goal of the raw data processing is to extract the value of the modulated power contained in the interferometric fringes. This value is proportional to the squared visibility V2 of the source on the observation baseline, which is in turn directly linked to the Fourier transform of the light distribution of the source through the Zernike-Van Cittert theorem.
One of the key advantages of VINCI is to use single-mode fibers to filter out the perturbations induced by the turbulent atmosphere. The wavefront that is injected in the fibers is only the mode guided by the fiber (Gaussian in shape, see Ruilier 1999 or Coudé du Foresto 1998b for details). The atmospherically corrupted part of the wavefront is not injected into the fibers and is lost into the cladding. Due to the temporal fluctuations of the turbulence, the injected flux changes considerably during an observation. However, VINCI derives two photometric signals that can be used to subtract the intensity fluctuations from the interferometric fringes and normalize them continuously. The resulting calibrated interferograms are practically free of atmospheric corruption, except the piston mode (differential longitudinal delay of the wavefront between the two apertures) that tends to smear the fringes and affect their visibility. Its effect is largely diminished by using a sufficiently high scanning frequency, as was the case for the VINCI observations.
After the photometric calibration has been achieved, the two interferograms from the two interferometric outputs of the VINCI beam combiner are subtracted to remove the residual photometric fluctuations. As the two fringe patterns are in perfect phase opposition, this subtraction removes a large part of the correlated fluctuations and enhances the interferometric fringes. 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 onto a sine wave (Fourier transform), but onto a function, i.e. the wavelet, that is localised 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 wavelet 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 backgrounds are removed
by making a least squares fit of the PSD at high and low frequency.
The calibration of the Cepheids' visibilities was achieved using
well-known calibrator stars that have been selected in the
Cohen et al. (1999) catalogue, with the exception of Ind.
This dwarf star was measured separately (Ségransan et al. 2004)
and used to calibrate one of the
Aql measurements.
The angular diameters of 39 Eri A, HR 4050 and HR 4546
(which belong to the Cohen et al. 1999 catalogue)
were also measured separately, as these stars appeared to give
a slightly inconsistent value of the interferometric efficiency.
For 39 Eri A and HR 4546, the measured angular diameters
we find are
and
mas, respectively.
These measured values are only
lower than the Cohen et al. (1999)
catalogue values of
and
mas.
A possible reason for this difference could be the presence of faint,
main sequence companions in orbit around these two giant stars.
The additional contribution of these objects would bias the diameter found
by spectrophotometry towards larger values, an effect consistent with what we observe.
For HR 4050, we obtained
mas,
only +1
away from the catalogue value of
mas.
The characteristics of the selected calibrators are listed in Table 2.
The limb-darkened disk (LD) angular diameters of these stars were converted into uniform disk
values using linear coefficients taken from Claret et al. (1995). As demonstrated by
Bordé et al. (), 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 calibrators were observed soon before and after the Cepheids,
in order to verify that the interferometric efficiency (IE) has not changed
significantly during the Cepheid observation itself. In some cases, and due to the technical
nature of commissioning observations, part of the Cepheid observations could not be
bracketed, but only immediately preceded or followed by a calibrator.
However, the stability of the IE has proved to be generally very good, and we
do not expect any significant bias from these single-calibrator observations.
Some observations included several calibrators to allow a cross-check of
of their angular sizes. The calibrators were chosen as close
as possible in the sky to our target Cepheids, in order to be able
to observe them with similar airmass.
This selection has taken into account the constraints in terms of limiting
magnitude and sky coverage imposed by the VLTI siderostats and delay lines.
The IE was computed from the coherence factor measurements obtained on
the calibrators, taking into account the bandwidth smearing effect (see Sect. 6.4)
and a uniform disk angular diameter model. This calibration process yielded the final
squared visibilities listed in Tables 3 to 9.
Name | mV | mK | Sp. Type |
![]() |
![]() |
![]() |
![]() |
![]() |
|
![]() |
HD 12524 | 5.16 | 1.52 | K5III | 3780 | 1.9 |
![]() |
![]() |
![]() |
39 Eri A | HD 26846 | 4.90 | 2.25 | K3III | 4210 | 2.2 |
![]() |
![]() |
![]() |
![]() |
HD 27442 | 4.44 | 1.97 | K2IVa | 4460 | 2.3 |
![]() |
![]() |
![]() |
HR 2533 | HD 49968 | 5.69 | 2.10 | K5III | 3780 | 1.9 |
![]() |
![]() |
![]() |
HR 2549 | HD 50235 | 5.00 | 1.39 | K5III | 3780 | 1.9 |
![]() |
![]() |
![]() |
![]() |
HD 55865 | 3.77 | 1.52 | K0III | 4720 | 2.6 |
![]() |
![]() |
![]() |
6 Pup | HD 63697 | 5.18 | 2.62 | K3III | 4210 | 2.2 |
![]() |
![]() |
![]() |
HR 3046 | HD 63744 | 4.70 | 2.31 | K0III | 4720 | 2.6 |
![]() |
![]() |
![]() |
HR 4050 | HD 89388 | 3.38 | 0.60 | K3IIa | 4335 | 2.3 | 4.43 ![]() |
![]() |
![]() |
HR 4080 | HD 89998 | 4.83 | 2.40 | K1III | 4580 | 2.5 |
![]() |
![]() |
![]() |
HR 4526 | HD 102461 | 5.44 | 1.77 | K5III | 3780 | 1.9 |
![]() |
![]() |
![]() |
HR 4546 | HD 102964 | 4.47 | 1.56 | K3III | 4210 | 2.2 |
![]() |
![]() |
![]() |
HR 4831 | HD 110458 | 4.67 | 2.28 | K0III | 4720 | 2.6 |
![]() |
![]() |
![]() |
![]() |
HD 145897 | 5.25 | 1.60 | K3III | 4210 | 2.2 |
![]() |
![]() |
![]() |
70 Aql | HD 196321 | 4.90 | 1.21 | K5II | 3780 | 1.9 |
![]() |
![]() |
![]() |
7 Aqr | HD 199345 | 5.50 | 2.00 | K5III | 3780 | 1.9 |
![]() |
![]() |
![]() |
![]() |
HD 209100 | 4.69 | 2.18 | K4.5V | 4580 | 4.5 |
![]() |
![]() |
![]() |
![]() |
HD 209688 | 4.48 | 1.68 | K3III | 4210 | 2.2 |
![]() |
![]() |
![]() |
HR 8685 | HD 216149 | 5.41 | 1.60 | M0III | 3660 | 1.4 |
![]() |
![]() |
![]() |
|
The VINCI processing pipeline produces a number of outputs to the user for the data quality control, including in particular the average wavelet power spectral density (WPSD) of the processed interferograms. This is an essential tool to verify that no bias is present in the calibrated and normalized fringe power peak. Figure 2 shows the average WPSD of a series of 302 interferograms obtained on X Sgr. No bias is present, and the residual background is very low. The power integration being done between 2000 and 8000 cm-1, the complete modulated power of the fringes is taken into account without bias.
X Sgr was observed 8 times on the B3-M0 baseline (140 m ground length), using
exclusively the two 0.35 m Test Siderostats (TS).
The projected baseline length varied between
118.4 and 139.7 m, and the observed squared visibilities were confined between
V2 = 56.9 and 71.1%. Thanks to its declination of
,
X Sgr culminates
almost at zenith over Paranal (-24
), and all the observations were
obtained at very low airmasses. It is located on the sky near two other Cepheids of our
sample, Y Oph and W Sgr, and these three stars share the same calibrator,
Sco.
The average signal to noise ratio (SNR) was typically 2 to 5 on the photometric outputs of
VINCI, and 4 to 6 on the interferometric channels, for a constant fringe frequency of 242 Hz.
A total of 4977 interferograms were processed by the pipeline.
The same remarks apply to W Sgr and Y Oph, as they have almost the same
magnitude and similar angular diameters. The number of processed interferograms
for these two stars was 4231 and 2182, respectively, during 9 and 4 observing sessions.
Aql was observed once on the E0-G1 baseline (66 m) and 10 times on the
B3-M0 baseline (140 m ground length). The total number of processed interferograms
is 5584. The SNRs were typically 4 and 7 on the photometric and interferometric
outputs at a fringe frequency of 242 to 272 Hz.
Due to its northern declination (
)
and to the limits of the TS, it was not possible to observe
Aql for
more than two hours per night, therefore limiting the number of interferograms and
the precision of the measurements.
Dor is a difficult target for observation with the TS, as it is partially hidden
behind the TS periscopes that are used to direct the light into the VLTI tunnels.
This causes a partial vignetting of the beams and therefore a loss in SNR.
The data from the TS are thus of intermediate quality, considering the
brightness of this star.
It is located at a declination of
,
relatively close to
Car, and
therefore these two stars share some calibrators. In addition to the 5 observations
with the TS, four measurements were obtained during three commissioning runs on the
UT1-UT3 baseline. A total of 8129 interferograms were processed, of which 5187 were
acquired with the 8 m Unit Telescopes (96 min spread over four nights were spent on
Dor using UT1 and UT3).
At a declination of ,
Gem is not accessible to the TS due to a
mechanical limitation. This is the reason why this star was observed only on two
occasions with UT1 and UT3, for a total of 3857 interferograms,
obtained during 41 min on the target. The average on-source SNRs were typically
50 for the interferometric channels and 30 for the photometric signals, at a
fringe frequency of 694 Hz.
The data obtained using the FLUOR/IOTA instrument are described in Kervella et al. (2001b). They were reprocessed using the latest release of the FLUOR software that includes a better treatment of the photon shot noise than the 2001 version. As the baseline of IOTA is limited to 38 m, the visibility of the fringes is very high, and the precision on the angular diameter is reduced compared to the 102.5 m baseline UT1-UT3.
As for Dor, the observation of
Car (
)
is made
particularly difficult by the vignetting of the TS beams.
Thanks to its brightness (
)
the SNRs are 15-20 on the interferometric
channels, and 10-15 on the photometric signals, using the TS and a fringe frequency of 242 Hz.
One observation was obtained on the E0-G1 baseline (66 m ground length),
and 19 measurements on the B3-M0 baseline.
Car is the most observed star in our
sample, with a total of 22 226 processed interferograms. Its average diameter of approximately
3 mas makes it an ideal target for observations with baselines of 100 to 200 m. On the
B3-M0 baseline, we achieved projected baselines of 89.7 to 135.0 m, corresponding
to V2 values of 8 to 42%. This range is ideal to constrain the visibility model
and derive precise values of the angular diameter.
Figure 3 shows the squared visibility points obtained
at two phases on Car. The change in angular diameter is clearly visible.
Thanks to the variation of the projected baseline on sky, we have sampled a
segment of the visibility curve.
![]() |
Figure 4:
Average intensity profiles computed from the four-parameter approximations
of Claret (2000) for X Sgr (thin line) and ![]() |
The object of this section is to derive the angular diameters of the Cepheids as a function of their pulsational phase. We discuss the different types of models that can be used to compute the angular diameter from the squared visibility measurements.
This very simple, rather unphysical model is commonly used for interferometric studies
as it is independent of any stellar atmosphere model. The relationship between the
visibility V and the uniform disk angular diameter (UD) is:
![]() |
(1) |
While the true stellar light distributions depart significantly from the UD model,
the UD angular diameters
given in
Tables 3 to 9
have the advantage that they can easily be converted to LD values using any stellar
atmosphere model. This is achieved by computing a conversion factor
from the chosen intensity profile
(see e.g. Davis et al. 2000 for details).
The visibility curve shape before the first minimum 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 chose were computed by Claret (2000), based on
model atmospheres by Kurucz (1992).
They consist of four-parameter approximations to the function
,
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.
The analytical expression of these approximations is given by:
![]() |
(2) |
The limb darkening is directly measurable by interferometry around the first minimum
of the visibility function, as demonstrated by several authors on giant stars
(Quirrenbach et al. 1996; Wittkowski et al. 2001).
Unfortunately, even for Car observed in the K band, this requires a baseline of more than 180 m
that was not available for the measurements reported here. It is intended in the near future
to measure directly the LD of a few nearby Cepheids, using the shorter wavelength bands of
AMBER (Petrov et al. 2000) and the longest baselines of the VLTI (up to 202 m).
As shown by Marengo et al. (2002),
the atmosphere of the Cepheids departs from that of a non-variable
giant with identical
and
,
due in particular to the
presence of energetic shock waves at certain phases of the pulsation.
However, this effect is enhanced at visible wavelengths compared to the
infrared, and appears to be negligible in the case of the VINCI observations.
Marengo et al. (2003) have derived in the H band a relative
variation of the limb darkening coefficient
of only 0.2%. This is below the precision of our measurements and is
neglected in the rest of this paper. Furthermore, the VINCI/VLTI measurement
wavelength being longer (2.18
m) than the H band, the LD correction
is even smaller, as is its expected variation.
From the results of Marengo et al. (2003) it appears clearly that the interferometers operating at infrared wavelengths are ideally suited for Cepheid measurements that aim at calibrating the P-R and P-L relations. On the other hand, as pointed out by these authors, the visible wavelength interferometers should be favored to study the dynamical evolution of the atmosphere (including the limb darkening) during the pulsation. The geometrical determination of the pulsation parallax is almost independant of the adopted atmosphere model in the K band, while this is not the case at shorter wavelengths.
The VINCI instrument bandpass corresponds to the K band filter, transparent between
and
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 bandwidth smearing
(Kervella et al. 2003b).
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.
In this paper, we assume that the limb darkening law does not change over the K band. This is
reasonable for a hot and compact stellar atmosphere, but is also coherent with the range of visibilities
measured on the Cepheids of our sample.
If necessary, this computation can easily be extended to a wavenumber dependant
intensity profile.
Following Davis et al. (2000), using a Hankel integral, we can derive the
visibility law
from the intensity profile:
![]() |
(3) |
![]() |
(4) |
![]() |
(5) |
![]() |
(7) |
The
model is adjusted numerically to the
observed (B, V2) data using a classical
minimization
process to derive
.
A single angular diameter is
derived per observation session, the fit being done directly on
the set of V2 values obtained during the session. The
systematic and statistical errors are considered separately in the
fitting procedure, to estimate the contribution
of the uncertainty of the calibrator diameter on the final error bar.
Each observation session was generally executed in less than 3 h, a short time compared to the pulsation periods of the Cepheids of our sample. Therefore, we do not expect any phase induced smearing from this averaging.
The reference epochs T0 and periods P for each Cepheid are given in
Table 1. N is the number of batches (500 interferograms)
recorded during the corresponding observing session.
For each angular diameter, the statistical and systematic calibration errors are given separately,
except for the FLUOR/IOTA measurements of Gem, for which the systematic calibration
error is negligible compared to the statistical uncertainty.
For each star we used radial velocity data found in the literature.
Specifically, we collected data from Bersier (2002) for Aql,
Car, and
Dor; from Bersier et al. (1994) for
Gem;
from Babel et al. (1989) for W Sgr. All these data have been
obtained with the CORAVEL radial velocity spectrograph
(Baranne et al. 1979).
We also obtained data from Evans & Lyons (1986) for Y Oph
and from Wilson et al. (1989) for X Sgr.
In theory, the linear diameter variation could be determined by direct integration of
pulsational velocities
(within the assumption that the
photosphere is comoving
with the atmosphere of the Cepheid during its pulsation).
However these velocities are deduced from the measured radial
velocities by the use of a projection factor p. The Cepheid's radii determined
from the BW method depend directly from a good knowledge of p.
Sabbey et al. (1995) and
Krockenberger et al. (1997) have studied
in detail the way to determine the p-factor.
We used a constant projection factor p = 1.36 in order to transform the radial
velocities into pulsation velocities. Burki et al. (1982) have
shown that this value is appropriate for the radial velocity measurements that
we used.
Star |
![]() |
d0 (pc) | ![]() |
X Sgr |
![]() |
![]() |
0.38 |
![]() |
![]() |
![]() |
3.98 |
W Sgr |
![]() |
![]() |
0.90 |
![]() |
![]() |
![]() |
1.31 |
![]() |
![]() |
![]() |
0.51 |
Y Oph |
![]() |
![]() |
0.16 |
(![]() |
![]() |
![]() |
23.2 |
Star |
![]() |
d1 (pc) | ![]() |
X Sgr |
![]() |
![]() |
1.36 |
![]() |
![]() |
![]() |
0.40 |
W Sgr |
![]() |
![]() |
0.42 |
![]() |
![]() |
![]() |
0.23 |
![]() |
![]() |
![]() |
0.88 |
Y Oph |
![]() |
![]() |
0.03 |
![]() |
![]() |
![]() |
0.71 |
Star |
![]() |
d2 (pc) | ![]() |
![]() |
![]() |
![]() |
0.43 |
W Sgr |
![]() |
![]() |
0.48 |
![]() |
![]() |
![]() |
0.25 |
![]() |
![]() |
![]() |
0.49 |
A large number of BW studies have been published, using both visible and
infrared wavelength observations. For Gem and
Aql, the pulsation
has been resolved using the Palomar Testbed Interferometer (Lane et al. 2000, 2002), thererefore giving a direct estimate of the diameter and distance of these
stars. Table 10 gives a list of the existing diameter
estimates for the Cepheids of our sample from the application of the classical
BW method ("B AADE-W ESSELINK'' section of the table).
From the many different P-R relations available, we chose the Gieren et al. (1998) version, as it is based on infrared colors for the determination of the temperature of the stars. Compared to visible colors, the infrared colors give a much less dispersed P-R relation. Indeed, this relation has a very good intrinsic precision of the order of 5 to 10% for the period range of our sample. Moreover, it is identical to the law determined by Laney & Stobie (1995). The compatibility with the individual BW diameter estimates is also satisfactory. The linear diameters deduced from this P-R law are mentioned in the "E MPIRICAL P-R" section of Table 10. We assume these linear diameter values in the following.
In the case of Car, the fit of a constant diameter results in a very high
value. This means that the average diameters
and
should not be used for further analysis.
The pulsation curve of this star is not sampled uniformly by our interferometric
observations, with more values around the maximum diameter.
This causes the larger diameter values to have more weight in the average diameter
computation, and this produces a significant positive bias. This remark does not
apply to the orders 1 and 2 fitting methods.
As a remark, no significant phase shift is detected at a level of
(14 min of time) between the predicted radius curve of
Car
and the observed angular diameter curve. The values of P and T0 used for
the fit are given in Table 3.
Figures 5 to 11 show the best models for each star,
together with the VINCI/VLTI angular diameter measurements for the seven
Cepheids of our sample. Figure 12 gives an enlarged view of
the maximum diameter of Car.
![]() |
Figure 6:
Order 2 model fit for ![]() |
![]() |
Figure 9:
Order 0 model fit for ![]() |
![]() |
Figure 12:
Detail of Fig. 11 showing the angular diameter
curve of ![]() |
From the NPOI (Armstrong et al. 2001; Nordgren et al. 2000),
PTI (Lane et al. 2002) and VINCI/VLTI measurements, we know the
average UD angular diameters of Aql and
Gem
at several effective wavelengths with high precision.
Table 14 gives the angular diameter values and the corresponding wavelengths.
Claret's (2000) linear limb darkening parameters u were used to
compute the expected conversion factors
.
To read the u table, we have considered the closest parameters to the average
values for
Aql and
Gem in Table 1,
and we computed
using the formula from Hanbury Brown et al. (1974):
![]() |
(9) |
We note that the value of
for
Aql that we derive for the
NPOI observation,
mas, is not identical to
the LD angular diameter originally given by Armstrong et al. (2001),
mas. There is a 1
difference, that may
be due to the different source of limb darkening coefficient that these authors
used for their modeling (Van Hamme 1993).
The resulting
values for the three observations are compatible
at the 2
level, but there is a slight trend that points towards an
underestimation of the limb darkening effect at shorter wavelengths,
or alternatively its overestimation at longer wavelengths.
Considering that the limb darkening is already small in the infrared,
the first hypothesis seems more plausible.
Marengo et al. (2002, 2003) have
shown that the Cepheids limb darkening can be significantly
different from stable giant stars, particularly at visible wavelengths.
This could explain the observed difference between the
0.73
m and K band diameters of
Aql and
Gem,
the latter being probably closer to the true LD diameters,
thanks to the lower limb darkening in the infrared.
In the case of Aql, another explanation
could be that the measurement at visible wavelengths
is biased by the blue companion of
Aql. However, it is 4.6 mag fainter than the Cepheid in the V band (Böhm-Vitense &
Proffitt 1985, see also Sect. 9.2),
and therefore should not contribute significantly to the visibility of the fringes.
Ref. | ![]() ![]() |
![]() |
![]() |
![]() |
![]() |
||||
(1) | 0.73 |
![]() |
1.048 |
![]() |
(2) | 1.65 |
![]() |
1.024 |
![]() |
(3) | 2.18 |
![]() |
1.021 |
![]() |
![]() |
||||
(1) | 0.73 |
![]() |
1.051 |
![]() |
(2) | 1.65 |
![]() |
1.027 |
![]() |
(3) | 2.18 |
![]() |
1.023 |
![]() |
As demonstrated by several authors (see Szabados 2003 for a complete database), binarity and multiplicity are common in the Cepheid class. Evans (1992) has observed that 29% of the Cepheids of her sample have detectable companions.
Our sample of Cepheids contains four confirmed binary Cepheids, out
of a total of seven stars.
As it is biased towards bright and nearby Cepheids,
this large fraction is an indication that many Cepheids currently believed to be
single could have undetected companions.
X Sgr (Szabados 1989b),
Aql (Böhm-Vitense & Proffitt 1985),
and W Sgr (Böhm-Vitense & Proffitt 1985; Babel et al. 1989)
are confirmed members of binary or multiple systems.
Gem is a visual binary star (Proust et al. 1981), but the
separated companion does not contribute to our observations.
Y Oph was once suspected to be a binary (Pel 1978),
but Evans (1992) has not confirmed the companion, and has
set an upper limit of A0 on its spectral type.
The physical parameters of the companions of Aql and W Sgr
have been derived by Böhm-Vitense & Proffitt (1985)
and Evans (1991), based on ultraviolet spectra.
The latter has derived spectral types of B9.8V and A0V, respectively.
The orbital parameters of the binary W Sgr were computed by
Babel et al. (1989) and Albrow & Cottrell (1996).
Based on IUE spectra, Evans (1992) has
set an upper limit of A0 on the spectral type of the companion of X Sgr.
The difference in V magnitude between these three Cepheids and their
companions is
.
The
is even larger due to the blue color of these stars,
.
Therefore, the effect on our visibility measurements is negligible, with
a potential bias of
.
For example, this translates into a
maximum error of
as on the average angular diameter of
Aql,
(a relative error of
), that is significantly smaller than our
error bars (
).
In the K band, the effect of the companions of the other
Cepheids is also negligible at the precision level of our measurements.
However, the presence of companions will have
to be considered for future measurements with angular diameter precisions
of a few
as. In this respect, long-period Cepheids, such
as
Car, are more reliable, as their intrinsic brightness is larger than
the short-period pulsators, and therefore they dominate their potential
companions even more strongly.
Fernie et al. (1995b) have found that the amplitude of the
light curve of Y Oph has been decreasing for a few decades. A similar behavior
has been observed only on Polaris (e.g. Evans et al. 2002).
The uncertainty on our
measurements has not allowed us to
detect unambiguously the pulsation of this star, but it is clearly an
important target for future observations using the Auxiliary Telescopes
(1.8 m) of the VLTI in order to estimate its parameters with high precision.
Interestingly, Gieren et al. () have studied the impact of binary Cepheids on their determination of the period-luminosity relation using 100 Cepheids, and they conclude that it is negligible. This is due to the very large intrinsic luminosity of the Cepheids that overshine by several orders of magnitude most of the other types of stars.
We have reported in this paper our long-baseline interferometric
observations of seven classical Cepheids using the VINCI/VLTI instrument.
For four stars ( Aql, W Sgr,
Dor and
Car),
we were able to apply a modified version of the BW method, resulting in an
independent estimate of their distance.
For all stars, we also derived their distances from lower order
fitting methods, that use an a priori estimate of their linear diameter from
the P-R relation of Gieren et al. (1998).
We would like to emphasize that the order 0/1 and order 2
error bars are different in nature, and they should be
treated differently in any further use of these results.
While the order 2 error bars can be treated as statistical
(i.e. reduced by averaging), the order 0/1 methods errors are
dominated by the systematic uncertainty introduced by the a priori
estimation of the linear radius.
The respective contributions of the statistical and systematic uncertainties
are given separately in Tables 11 and 12.
These values assume a constant value of the p-factor of 1.36, and can
be scaled linearly for other values.
We will use these distances in Paper II, together with previously published measurements,
to calibrate the zero points of the Period-Radius and Period-Luminosity relations.
In Paper III, we will calibrate the surface brightness-color relation, with a particular emphasis on the
evolution of Car in this diagram over its pulsation. These three empirical
relations are of critical importance for the extragalactic distance scale.
The direct measurement of the limb darkening of nearby Cepheids by interferometry is the next step of the interferometric study of these stars. It will allow a refined modeling of the atmosphere of these stars. This observation will be achieved soon using in particular the long baselines of the VLTI equipped with the AMBER instrument, and the CHARA array for the northern Cepheids. Another improvement of the interferometric BW methow will come from radial velocity measurements in the near infrared (see e.g. Butler & Bell 1997). They will avoid any differential limb darkening between the interferometric and radial velocity measurements, and therefore make the resulting distances more immune to limb darkening uncertainties.
Acknowledgements
DB acknowledges support from NSF grant AST-9979812. PK acknowledges support from the European Southern Observatory through a postdoctoral fellowship. Based on observations collected at the European Southern Observatory, Cerro Paranal, Chile, in the framework of ESO shared-risk programme 071.D-0425 and unreferenced commissioning programme in P70. The VINCI/VLTI public commissioning data reported in this paper have been retrieved from the ESO/ST-ECF Archive (Garching, Germany). This work has made use of the wavelet data processing technique, developed by D. Ségransan (Observatoire de Genève), and embedded in the VINCI pipeline. This research has made use of the SIMBAD database at CDS, Strasbourg (France). We are grateful to the ESO VLTI team, without whose efforts no observation would have been possible.
JD | Stations | Baseline | Phase |
![]() |
![]() |
N |
![]() |
Calibrators |
(m) | ![]() ![]() |
![]() ![]() |
||||||
2 452 741.903 | B3-M0 | 138.366 | 0.560 |
![]() |
![]() |
2 | 0.66 | ![]() |
2 452 742.885 | B3-M0 | 137.432 | 0.700 |
![]() |
![]() |
3 | 0.52 | ![]() |
2 452 743.897 | B3-M0 | 137.903 | 0.844 |
![]() |
![]() |
3 | 0.08 | ![]() |
2 452 744.868 | B3-M0 | 139.657 | 0.983 |
![]() |
![]() |
2 | 0.09 | ![]() |
2 452 747.848 | B3-M0 | 139.530 | 0.408 |
![]() |
![]() |
1 | - | ![]() |
2 452 749.832 | B3-M0 | 139.084 | 0.691 |
![]() |
![]() |
2 | 0.35 | ![]() |
2 452 766.811 | B3-M0 | 138.853 | 0.112 |
![]() |
![]() |
4 | 0.09 | ![]() |
2 452 768.877 | B3-M0 | 128.228 | 0.406 |
![]() |
![]() |
6 | 0.62 | ![]() |
JD | Stations | Baseline | Phase |
![]() |
![]() |
N |
![]() |
Calibrators |
(m) | ![]() ![]() |
![]() ![]() |
||||||
2 452 524.564 | E0-G1 | 60.664 | 0.741 |
![]() |
![]() |
3 | 0.08 | 70 Aql |
2 452 557.546 | B3-M0 | 137.625 | 0.336 |
![]() |
![]() |
1 | - | ![]() |
2 452 559.535 | B3-M0 | 138.353 | 0.614 |
![]() |
![]() |
1 | - | 7 Aqr, ![]() |
2 452 564.532 | B3-M0 | 136.839 | 0.310 |
![]() |
![]() |
3 | 0.42 | 7 Aqr, ![]() |
2 452 565.516 | B3-M0 | 138.495 | 0.447 |
![]() |
![]() |
3 | 0.13 | 7 Aqr |
2 452 566.519 | B3-M0 | 137.845 | 0.587 |
![]() |
![]() |
5 | 0.23 | 7 Aqr |
2 452 567.523 | B3-M0 | 137.011 | 0.727 |
![]() |
![]() |
2 | 0.62 | 7 Aqr |
2 452 573.511 | B3-M0 | 136.303 | 0.561 |
![]() |
![]() |
1 | - | ![]() |
2 452 769.937 | B3-M0 | 139.632 | 0.931 |
![]() |
![]() |
3 | 0.06 | ![]() |
2 452 770.922 | B3-M0 | 139.400 | 0.068 |
![]() |
![]() |
3 | 0.15 | ![]() |
2 452 772.899 | B3-M0 | 138.188 | 0.343 |
![]() |
![]() |
3 | 0.16 | 7 Aqr |
JD | Stations | Baseline | Phase |
![]() |
![]() |
N |
![]() |
Calibrators |
(m) | ![]() ![]() |
![]() ![]() |
||||||
2 452 743.837 | B3-M0 | 137.574 | 0.571 |
![]() |
![]() |
1 | - | ![]() |
2 452 744.915 | B3-M0 | 137.166 | 0.713 |
![]() |
![]() |
2 | 0.04 | ![]() |
2 452 749.868 | B3-M0 | 139.632 | 0.365 |
![]() |
![]() |
1 | - | ![]() |
2 452 751.866 | B3-M0 | 139.538 | 0.628 |
![]() |
![]() |
1 | - | ![]() |
2 452 763.888 | B3-M0 | 131.830 | 0.211 |
![]() |
![]() |
4 | 0.73 | ![]() |
2 452 764.856 | B3-M0 | 135.926 | 0.339 |
![]() |
![]() |
4 | 0.76 | ![]() |
2 452 765.880 | B3-M0 | 132.679 | 0.473 |
![]() |
![]() |
4 | 1.43 | ![]() |
2 452 767.867 | B3-M0 | 132.637 | 0.735 |
![]() |
![]() |
3 | 0.01 | ![]() |
2 452 769.914 | B3-M0 | 120.648 | 0.005 |
![]() |
![]() |
2 | 0.33 | ![]() |
JD | Stations | Baseline | Phase |
![]() |
![]() |
N |
![]() |
Calibrators |
(m) | ![]() ![]() |
![]() ![]() |
||||||
2 452 215.795 | U1-U3 | 89.058 | 0.161 |
![]() |
![]() |
3 | 0.03 | ![]() ![]() |
2 452 216.785 | U1-U3 | 89.651 | 0.261 |
![]() |
![]() |
7 | 0.10 | ![]() |
2 452 247.761 | U1-U3 | 83.409 | 0.408 |
![]() |
![]() |
5 | 0.40 | ![]() |
2 452 308.645 | U1-U3 | 75.902 | 0.594 |
![]() |
![]() |
5 | 1.01 | HD 63697 |
2 452 567.827 | B3-M0 | 134.203 | 0.927 |
![]() |
![]() |
1 | - | HR 2549 |
2 452 744.564 | B3-M0 | 89.028 | 0.884 |
![]() |
![]() |
2 | 0.09 | HR 3046, 4831 |
2 452 749.514 | B3-M0 | 98.176 | 0.387 |
![]() |
![]() |
3 | 0.11 | HR 3046 |
2 452 750.511 | B3-M0 | 98.189 | 0.488 |
![]() |
![]() |
2 | 0.24 | HR 3046 |
2 452 751.519 | B3-M0 | 95.579 | 0.591 |
![]() |
![]() |
3 | 0.03 | HR 3046 |
JD | Stations | B, SF | Phase |
![]() |
![]() |
N |
![]() |
Calibrators |
![]() ![]() |
![]() ![]() |
|||||||
2 452 214.879 | U1-U3 | 82.423 | 0.408 |
![]() |
![]() |
8 | 0.25 | 39 Eri |
2 452 216.836 | U1-U3 | 72.837 | 0.600 |
![]() |
![]() |
6 | 0.28 | 39 Eri, ![]() |
2 451 527.972 | IOTA-38m | 84.870 | 0.739 |
![]() |
![]() |
1 | - | HD 49968 |
2 451 601.828 | IOTA-38m | 83.917 | 0.014 |
![]() |
![]() |
3 | 0.02 | HD 49968 |
2 451 259.779 | IOTA-38m | 83.760 | 0.318 |
![]() |
![]() |
1 | - | HD 49968 |
2 451 262.740 | IOTA-38m | 84.015 | 0.610 |
![]() |
![]() |
2 | 0.13 | HD 49968 |
2 451 595.863 | IOTA-38m | 83.790 | 0.427 |
![]() |
![]() |
2 | 1.72 | HD 49968 |
2 451 602.764 | IOTA-38m | 85.010 | 0.107 |
![]() |
![]() |
2 | 0.02 | HD 49968 |
JD | Stations | Baseline | Phase |
![]() |
![]() |
N |
![]() |
Calibrators |
(m) | ![]() ![]() |
![]() ![]() |
||||||
2 452 742.906 | B3-M0 | 139.569 | 0.601 |
![]() |
![]() |
2 | 0.10 | ![]() |
2 452 750.884 | B3-M0 | 139.057 | 0.067 |
![]() |
![]() |
2 | 0.41 | ![]() |
2 452 772.831 | B3-M0 | 139.657 | 0.349 |
![]() |
![]() |
3 | 0.22 | ![]() |
2 452 782.186 | B3-M0 | 129.518 | 0.168 |
![]() |
![]() |
4 | 0.30 | ![]() |
JD | Stations | Baseline | Phase |
![]() |
![]() |
N |
![]() |
Calibrators |
(m) | ![]() ![]() |
![]() ![]() |
HR | |||||
2 452 453.498 | E0-G1 | 61.069 | 0.587 |
![]() |
![]() |
4 | 0.01 | 4050 |
2 452 739.564 | B3-M0 | 130.468 | 0.634 |
![]() |
![]() |
2 | 0.03 | 4526 |
2 452 740.569 | B3-M0 | 128.821 | 0.662 |
![]() |
![]() |
7 | 0.77 | 4526 |
2 452 741.717 | B3-M0 | 96.477 | 0.694 |
![]() |
![]() |
5 | 0.28 | 4526 |
2 452 742.712 | B3-M0 | 99.848 | 0.722 |
![]() |
![]() |
5 | 0.09 | 4526 |
2 452 743.698 | B3-M0 | 99.755 | 0.750 |
![]() |
![]() |
2 | 0.08 | 4831 |
2 452 744.634 | B3-M0 | 114.981 | 0.776 |
![]() |
![]() |
6 | 0.73 | 4831 |
2 452 745.629 | B3-M0 | 115.791 | 0.804 |
![]() |
![]() |
2 | 0.01 | 3046, 4546, 4831 |
2 452 746.620 | B3-M0 | 116.828 | 0.832 |
![]() |
![]() |
5 | 0.65 | 3046, 4546 |
2 452 747.599 | B3-M0 | 120.812 | 0.860 |
![]() |
![]() |
3 | 0.70 | 4546, 4831 |
2 452 749.576 | B3-M0 | 124.046 | 0.915 |
![]() |
![]() |
4 | 1.18 | 4546 |
2 452 751.579 | B3-M0 | 122.555 | 0.971 |
![]() |
![]() |
4 | 1.16 | 3046, 4831 |
2 452 755.617 | B3-M0 | 112.185 | 0.085 |
![]() |
![]() |
1 | - | 4831 |
2 452 763.555 | B3-M0 | 120.632 | 0.308 |
![]() |
![]() |
6 | 1.02 | 4546 |
2 452 765.555 | B3-M0 | 119.629 | 0.365 |
![]() |
![]() |
6 | 1.19 | 4546 |
2 452 766.550 | B3-M0 | 120.005 | 0.393 |
![]() |
![]() |
7 | 0.99 | 4546 |
2 452 768.566 | B3-M0 | 115.135 | 0.450 |
![]() |
![]() |
7 | 0.46 | 4546 |
2 452 769.575 | B3-M0 | 113.082 | 0.478 |
![]() |
![]() |
3 | 0.03 | 3046, 4831 |
2 452 770.535 | B3-M0 | 121.152 | 0.505 |
![]() |
![]() |
2 | 0.20 | 3046, 4831 |
2 452 771.528 | B3-M0 | 122.014 | 0.533 |
![]() |
![]() |
3 | 0.88 | 4831 |
X Sgr | ![]() |
W Sgr | ![]() |
![]() |
Y Oph | ![]() |
|
INTERFEROMETRY | |||||||
Kervella et al. (2001b)
![]() |
63+35-19 | ||||||
Lane et al. (2002) |
![]() |
![]() |
|||||
Nordgren et al. (2000)
![]() |
69+28-15 | 60+25-14 | |||||
BAADE-WESSELINK | |||||||
Bersier et al. (1997) |
![]() |
![]() |
|||||
Fouqué et al. (2003) |
![]() |
![]() |
|||||
Krockenberger et al. (1997) |
![]() |
69.1+5.5-4.8 | |||||
Laney & Stobie (1995) |
![]() |
![]() |
![]() |
||||
Moffett & Barnes (1987)a |
![]() |
![]() |
![]() |
![]() |
|||
Moffett & Barnes (1987)b |
![]() |
![]() |
![]() |
![]() |
|||
Sabbey et al. (1995)c |
![]() |
![]() |
![]() |
||||
Sabbey et al. (1995)d |
![]() |
![]() |
![]() |
||||
Sachkov et al. (1998) | ![]() |
||||||
Taylor et al. (1997) |
![]() |
||||||
Taylor & Booth (1998) |
![]() |
||||||
Turner & Burke (2002) |
![]() |
![]() |
|||||
Sasselov & Lester (1990) | ![]() |
![]() |
|||||
MEAN B-W (overall ![]() |
52.5 (11.4) | 59.9 (5.7) | 57.0 (3.4) | 65.8 (7.2) | 65.3 (9.8) | 92.2 (-) | 180 (-) |
EMPIRICAL P-R | |||||||
Gieren et al. (1998) |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |