Free Access
Issue
A&A
Volume 593, September 2016
Article Number A45
Number of page(s) 17
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201528005
Published online 12 September 2016

© ESO, 2016

1. Introduction

The Baade-Wesselink (BW) method for the distance determination of Cepheids is, in its first version, purely spectro-photometric (Lindemann 1918; Baade 1926; Wesselink 1946). Results by Fouqué et al. (2007) and Storm et al. (2011a,b) illustrate the central and current role of this method in the distance scale calibration, even if recent calibrations of the Hubble constant rely exclusively on the trigonometric parallaxes of a few Galactic Cepheids (Riess et al. 2011; Benedict et al. 2007). The first interferometric version of the BW method was attempted in the visible by Mourard et al. (1997), and soon afterward in the infrared by Kervella et al. (1999, 2001) and Lane et al. (2000). Since then, the infrared interferometric BW method has been applied to a significant number of Cepheids, twelve in total (see Table 1), while among these stars only four have been observed by visible interferometers, and the pulsation could actually be resolved for only one of them, Car (Davis et al. 2009).

The principle of the interferometric version of the BW method is simple. Interferometric measurements lead to angular diameter estimations over the whole pulsation period, while the stellar radius variations can be deduced from the integration of the pulsation velocity. The latter is linked to the observational velocity deduced from spectral line profiles by the projection factor p (Nardetto et al. 2004; Mérand et al. 2005; Nardetto et al. 2007, 2009).

Table 1

Cepheids for which an interferometric BW method has been applied.

There are several underlying assumptions to the BW method. First, the limb-darkening of the star is assumed to be constant during the pulsation cycle. This has no impact on the interferometric analysis, at least in the infrared (Kervella et al. 2004), while in the visible Davis et al. (2009) reported the need to take into account a limb-darkening variation from Sydney University Stellar Interferometer (SUSI) observations. On the theoretical side, Nardetto et al. (2006a) found via hydrodynamical simulation that considering a constant limb darkening in the visible leads to a systematic shift of about 0.02 in phase on the angular diameter curve, which basically means that there is no impact on the derived distance (because the amplitude of the angular diameter curve is unchanged). Second, the projection-factor, mostly dominated by the limb-darkening calculated in the visible domain, is also assumed to be constant, which seems reasonable, at least in theory (Nardetto et al. 2004). Third, when applying the BW method, visible spectroscopy (e.g. Nardetto et al. 2006b) is often combined with infrared interferometric data, or even – in the recent photometric version of the BW method – with various photometric bands (Breitfelder et al. 2015). In the distance determination, we implicitly assume that the angular and linear diameters correspond to the same physical layer in the star. In this context, testing these hypotheses using visible interferometric observations seems to be of prime importance. This can be done first by simply verifying the internal consistency of the BW distances derived from visible and infrared interferometry, and then comparing them with available precise parallaxes (Benedict et al. 2007; Majaess et al. 2012).

This paper is the first in a series which investigates Cepheids with the Visible spEctroGraph and polArimeter (VEGA) beam combiner (Mourard et al. 2009) operating at the focus of the Center for High Angular Resolution Astronomy (CHARA) array (ten Brummelaar et al. 2005) located at the Mount Wilson Observatory (California, USA). The study focuses on the prototype δ Cep star. The VEGA data are first reduced (Sect. 2) and then analysed in terms of uniform disk angular diameters (Sect. 3). In Sect. 4, we show that δ Cep is clearly surrounded by a resolved structure, while some evidence in our interferometric data point toward an additional physical effect. We explore two hypotheses: a circumstellar reverberation and a strong limb-darkening variation. We draw our conclusions in Sect. 5.

Table 2

List and properties of calibration stars selected with the SearchCal software provided by the Jean-Marie Mariotti Center (JMMC) (Bonneau et al. 2006, 2011).

2. VEGA/CHARA observations of δ Cep

The CHARA array consists of six telescopes of 1 m in diameter, spread in a Y-shaped configuration, which offers 15 different baselines from 34 m to 331 m. These baselines can achieve a spatial resolution up to 0.3 mas in the visible. The interferometric observations were secured with the VEGA/CHARA instrument. The journal of observations is presented in Tables A.1 to A.2. Given the large angular size of δ Cep, we used only two short baselines, i.e. S1S2 and E1E2, with projected baselines ranging from about 27 to 31 m for S1S2 and from 52 to 66 m for E1E2. The data were processed using the standard VEGA pipeline (Mourard et al. 2009, 2011; Ligi et al. 2013) considering different spectral bands from 3 nm (in high spectral resolution mode) to 20 nm (in medium spectral resolution mode), and with reference wavelengths ranging from 500 nm to 745 nm (the minimum and maximum wavelengths of these bands are given in Tables A.1 to A.2).

In order to calibrate the squared visibilities, we considered the references stars HD 214734 (C1) and HD 213558 (C2). The first was used to calibrate the data of β Cep presented in Nardetto et al. (2011). Several nights of VEGA observations in 2008 and 2009 were also devoted to controlling the quality of these two calibrators by comparing their transfer functions (defined as the ratio of the actual to expected squared visibility) with those of six other calibrators listed in Table 2 (hereafter C3 to C8). The consistency among the eight calibrators is shown in Fig. A.1, while the general quality of VEGA/CHARA data is illustrated in Fig. A.2.

Cycle-to-cycle variations have never been detected for δ Cep, although they have been for a few other Cepheids (Anderson 2014; Anderson et al. 2016). The data are thus recomposed into a unique cycle and the data corresponding to nights that are close in pulsation phase are merged (see phases 0.180, 0.526, 0.627 in Tables A.1 to A.2). Delta Cephei was finally observed on 20 separate nights corresponding to 17 different pulsation phases. As the observations over a given night can be spread over several hours, we first derive the phase for each observation during the night and then calculate the average and the corresponding standard deviation. The standard deviations of the pulsation phases are rather low and range from 0.001 to 0.017, except for φ = 0.999 where we get an error of 0.025 because the star was observed at the very beginning and at the very end of the night. The pulsation phases can be found in Table 3 together with their respective errors, the number of visibility measurements (N), and the baseline used. For each calibrated visibility, the statistical and systematic calibration errors are given separately. The systematic calibration errors, owing to the uncertainty on the estimation of the diameter of the reference star, were found to be negligible compared to the statistical values. Therefore, we only considered the statistical uncertainties in the model fitting. In a few cases, the statistical uncertainty was clearly underestimated. Thus, we fixed the uncertainty on the calibrated squared visibility to 0.05 (Mourard et al. 2012) for the following nights: 2012 September 23, 2013 October 26, 2013 November 26, 2014 July 02, 2014 July 05, and 2014 July 09.

Table 3

Pulsation phases, number of visibility measurements (N), resulting uniform disk angular diameters, and reduced χ2.

thumbnail Fig. 1

Uniform disk diameters θUD derived from the VEGA/CHARA E1E2 data plotted as a function of the pulsation phase (blue crosses). The S1S2 data (corresponding to four different pulsation phases, see Table 3) were merged and the corresponding mean θUD is represented by a black horizontal line together with the 1σ uncertainty (blue dashed zone). The uniform disk angular diameter curve from FLUOR/CHARA in the K band (Mérand et al. 2005) is also shown for comparison. The red crosses and the light blue squares show the K band and R band UD angular diameter curves, respectively, as predicted from the SPIPS code (Merand et al. 2015a).

3. Uniform disk (UD) angular diameters

In Figs. A.3, the calibrated visibilities corresponding to E1E2 measurements are plotted as a function of the spatial frequency together with the corresponding (u,v) coverage and for each pulsation phase. We fit these calibrated visibilities by a uniform disk using a JMMC1 tool, LITpro (Tallon-Bosc et al. 2008). The results are given in Table 3 and plotted in Fig. 1 with blue crosses. The data corresponding to S1S2 measurements are shown in Fig. 2a. Owing to their low spatial frequencies, these S1S2 measurements are almost not sensitive to the pulsation of the star. This can be seen in Fig. 2b. Moreover, we have only a few S1S2 measurements per pulsation phase (four at most). Consequently, as a point of comparison with the E1E2 data, we merge the S1S2 data and fit them with a uniform disk angular diameter. The result is given in Table 3 and plotted in Fig. 1 by a horizontal dashed zone. For the data of 2014 October 24, which include both E1E2 and S1S2 baselines, we made two fits, one including all the S1S2 data and another with the E1E2 data. We overplot the very precise K-band uniform disk angular diameter curve obtained by Mérand et al. (2005) (magenta dots). To analyse our VEGA/CHARA angular diameter measurements, we consider the Spectroscopy-Photometry-Interferometry for Pulsating Stars algorithm (SPIPS; Merand et al. 2015a). The SPIPS code combines all the available observables of δ Cep: radial velocimetry (Bersier et al. 1994; Storm et al. 2004), interferometry (FLUOR/CHARA data), and photometry in the V (Berdnikov & Turner 2002; Engle et al. 2014; Kiss 1998; Moffett & Barnes 1984), J, H, and K bands (Barnes et al. 1997) in order to estimate the projection factor, the variation of the effective temperature, and the H and K band excesses. The data used for the fit for δ Cep have been published in digital form (Merand et al. 2015b). In the list of these outputs, we are particularly interested in the uniform disk angular diameter curves, calculated at 2200 nm and 800 nm, and derived directly from the framework of Merand et al. (2015a). The data are given in Table 4 and are shown in Fig. 1 by red crosses and light blue squares, respectively. In the infrared, the K-band UD angular diameter curve of SPIPS is consistent with the FLUOR/CHARA measurements. In the visible domain, the reduced χ2 between the VEGA (E1E2) and SPIPS UD angular diameters is 9, while it increases slightly to 11 when we replace the SPIPS angular diameter variation by a constant corresponding to the average of the VEGA UD angular diameters, i.e. 1.376 ± 0.030 mas. We thus detect a pulsation in the visible band, but make three remarks:

Table 4

Angular diameters as derived from the application of the SPIPS algorithm to δ Cep, using photometry, interferometry, and the cross-correlated radial velocity curves as inputs (see Merand et al. 2015a).

thumbnail Fig. 2

a) S1S2 observed calibrated squared visibilities secured at four different pulsation phases plotted as a function of the spatial frequency together with the best fit of uniform disk model (solid blue line). In the upper right corner we show the corresponding (u,v) coverage. b) Same data fitted with a two-component model (open symbols) including a uniform pulsating disk and a CSE (see text). As a comparison, the solid lines are the corresponding one-component models, i.e. a uniform pulsating disk (in the R band) without the CSE.

  • 1.

    The VEGA S1S2 measurements provide a mean θUD (with a reduced χ2 of 0.9), which is significantly larger than the E1E2 angular diameter curve by about 10σ. This suggests the presence of a large angular structure around the pulsating Cepheid.

  • 2.

    The SPIPS predictions are consistent with the VEGA θUD angular diameters from 0.0 to 0.8 in phase (with a reduced χ2 of 4). This basically means that δ Cep is pulsating in a quasi-hydrostatic way over this range of pulsation phase, even if the two measurements (at φ = 0.741 and φ = 0.783) are slightly below the SPIPS curve (around 2σ).

  • 3.

    Five measurements between 0.8 and 1.0 are clearly from 2 to 6σ above the SPIPS curve (with a reduced χ2 of 8).

4. Detection of a circumstellar environment or a nebulae in the visible

In Fig. 2b, the S1S2 data are fitted with a model (open symbols) composed of the UD derived from the SPIPS UD angular diameter curve at the corresponding phases of VEGA observations and a second static disk describing the circumstellar environment (CSE). The orientation of the S1S2 measurements in the (u,v) coverage are between 148 and 171 degrees, preventing us from investigating a non-centrosymmetric CSE. The best-fit parameters are θCSE = 8.9 ± 3.0 mas (the size of the CSE) and fcsef=0.07±0.01\hbox{$ \frac{f_\mathrm{cse}}{f_\mathrm{\star}}=0.07\pm0.01$} (the relative flux contribution of the CSE, where fcse and f are the fluxes of the CSE and the star in the R band, respectively) with a reduced χ2 of 0.6. The two parameters, however, are significantly correlated (0.8) owing to the low number of S1S2 observations. Interestingly, we can also fit the S1S2 data (with the same reduced χ2 of 0.6) using a model composed of a pulsating disk, as in the previous case, together with a background contribution of 0.07 ± 0.01 in flux filling the field of view of the interferometer. This model physically corresponds to the visible counterpart of the infrared nebulae discovered by Marengo et al. (2010) and the HI nebula found with the Very Large Array (VLA) (Matthews et al. 2012). Unfortunately, our S1S2 data are too scarce to distinguish between these two hypotheses.

thumbnail Fig. 3

Panel a), the VEGA squared visibilities plotted as a function of x=πBp[m]θUD[mas]λ[nm]\hbox{$x = \frac{\pi B_{\rm p} \mathrm{[m]} \theta_\mathrm{UD} \mathrm{[mas]} }{\lambda \mathrm{[nm]}}$} in order to allow a comparison between data of different phases (i.e. different θUD) and effective wavelengths. Five models are overplotted: (1) the UD pulsating disk without CSE from the SPIPS algorithm (blue solid line); (2) the SPIPS UD pulsating disk + the CSE from Eq. (1) (red dotted line); (3) the SPIPS UD pulsating disk + the CSE from Eq. (1) but two times brighter (dashed line); (4) the SPIPS UD pulsating disk + the CSE or Eq. (1) but two times fainter (dash-dotted line); (5) the SPIPS UD pulsating disk + a background contribution filling the field of view of the interferometer (magenta dotted line). The residual between the observations and models 1 to 4 are shown in panels b), c), d), and e), respectively. In panels d) and e) we show only the data corresponding to phases intervals φ = 0.0−0.8 and φ = 0.8−1.0, respectively.

We have fitted the S1S2 data with a two-component model without considering the E1E2 data. Hence, we now have to verify whether this two-component model is consistent with our E1E2 data, in particular at minimum radius. In Sect. 3, we find that the E1E2 measurements are not consistent with a quasi-hydrostatic pulsating disk at the minimum radius of the star (between phase 0.8 and 1.0).

In order to address this issue, we plot in panel a) of Fig. 3 all the VEGA visibilities measured as a function of x=πBp[m]θUD[mas]λ[nm]\hbox{$x=\frac{\pi B_{\rm p} \mathrm{[m]} \theta_\mathrm{UD} \mathrm{[mas]}}{\lambda \mathrm{[nm]}}$}, where θUD is fixed and interpolated from the SPIPS UD angular diameter curve at the corresponding pulsation phase of VEGA observations. The data are rescaled in such a way that they can be compared, despite different θUD (or pulsation phases) and different wavelengths of observation. This is possible since the reference is the SPIPS θUD semi-theoretical curve represented by the solid blue line expressed by V2(x)=|2J1(x)x|\hbox{$V^2(x)= | \frac{2 J_1(x)}{x} |$}, where J1 is the Bessel function of the first order. This concept of pseudo-baseline was first introduced by Mérand et al. (2006). Using this approach, we confirm the three statements found in Sect. 3:

  • 1.

    The S1S2 measurements, whatever the pulsation phase, are significantly lower than in the rescaled uniform disk model (see the S1S2 measurements, i.e. data with x from approximately 0.5 to 1.4, Fig. 3b). This deviation is removed as soon as we consider a resolved structure around δ Cep (see below).

  • 2.

    The data with a pulsating phase from φ = 0.0 to φ = 0.8 are consistent with the uniform disk model for E1E2 measurements (blue crosses with x from 1.5 to 2.8 in Fig. 3b).

  • 3.

    The data with a pulsating phase from φ = 0.8 to φ = 1.0 are not consistent with the uniform disk model for E1E2 measurements (green circles around x = 1.7 in Fig. 3b).

In Fig. 3a, we overplot a red dotted line to the squared visibility curve corresponding to a uniform disk surrounded by the CSE using the formula V2(x)=(f|2J1(x)x|+fcse|2J1(srx)srx|)2,\begin{equation} V^2(x)=\left(f_{\star} \left| \frac{2 J_1(x)}{x} \right|+f_{\mathrm{cse}} \left| \frac{2 J_1({s_{\rm r}} x)}{{s_{\rm r}} x} \right| \right)^2 , \end{equation}(1)where fcse = 0.065 and f = 0.935 (corresponding to fcsef=0.07±0.01\hbox{$ \frac{f_\mathrm{cse}}{f_\mathrm{\star}}=0.07\pm0.01$} and fcse + f = 1), sr is the mean size ratio between the CSE and the pulsating disk. Considering the mean UD angular diameter of the SPIPS curve (θmean = 1.39 mas), we find sr = 6.32. The result is unchanged (i.e. within the width of the line) if we use the 6.75 and 6.07 values corresponding to the minimum and maximum UD angular diameters.

We find that the S1S2 measurements are properly fitted by this two-component model, as expected (Fig. 3c). However, around x = 1.7 the E1E2 measurements corresponding to φ = 0.0−0.8 and φ = 0.8−1.0 are respectively above and below the curve (by about 2σ). All the E1E2 measurements are partly above the curve (about 1σ) at larger frequencies (x> 2.3).

If we consider the model composed of a uniform disk surrounded by a background, we obtain V2(x)=(f|2J1(x)x|)2fcse2,\begin{equation} V^2(x)=\left(f_{\star} \left| \frac{2 J_1(x)}{x} \right|\right)^2-f_{\mathrm{cse}}^2\,, \end{equation}(2)and the corresponding curve is plotted by a magenta dotted line in Fig. 3a. We arrive at the same conclusions as the model of the pulsating uniform disk surrounded by a CSE. In particular, we find the presence of two regimes in phase, even considering a CSE or a background. We explore two hypotheses.

First, we cannot fit all the VEGA measurements at the same time if we consider a CSE two times fainter or brighter (see black dotted and dot-dashed lines in Fig. 3a). Qualitatively, however, it seems that considering a CSE two times fainter between phase 0.0 and 0.8 is a good compromise to fit the S1S2 and E1E2 measurements together, even if not totally satisfactory (see Fig. 3d). To consider a CSE two times brighter between phase 0.8 and 1.0 (when the size of the star is at a minimum, and is hot and bright) helps to fit the data except at large spatial frequencies (see Fig. 3e). If this hypothesis is correct, then the δ Cep would light up its environment (CSE or background) differently at maximum and minimum radius.

The second possibility is to consider that the CSE (or the background) has a constant brightness and that the discrepancy found for E1E2 measurements comes mainly from an additional structure angularly smaller (φ = 0.0−0.8) or larger (φ = 0.8−1.0) in size than the SPIPS θUD stellar pulsating disk, for instance a strong limb-darkening variation. The limb-darkening coefficient, defined as k=θUDθLD\hbox{$ k=\frac{\theta_{\mathrm{UD}}}{\theta_{\mathrm{LD}}}$}, is used to convert the θUD diameters into limb-darkened θLD angular diameters. The k coefficient variation in the visible band (at 600 nm) and over the cycle of pulsation of the star is found to be from the hydrodynamical model of 0.015 (Nardetto et al. 2006a), which corresponds to about 1%. Similarly, Marengo et al. (2002, 2003) have found a variation for the k-parameter of about 0.01 and 0.02 for ζ Gem and at 570 nm using quasi-hydrostatic and hydrodynamical models, respectively. Also for comparison, Davis et al. (2009) used a phase-dependent k-parameter to analyse their SUSI data of Car, and found a variation of 0.012. This means typically a 1% difference in θUD, and thus a 1% difference in x, which corresponds to 0.02 in absolute value around x = 1.7 in Fig. 3. Consequently, the limb-darkening variation of δ Cep is one order of magnitude lower than the discrepancies found for the E1E2 measurements. In the case of δ Cep, the theoretical average value of k is of 0.954 (Nardetto et al. 2006a). So far, the only published direct measurement of the limb-darkening coefficient of a Cepheid (in the visible) has been performed by Pilecki et al. (2013) for a Cepheid (P = 3.80 d) in an eclipsing binary system in the LMC. They find the limb-darkening effect to be much stronger than expected for a star of that temperature, giving a possible range for k from 0.91 to 0.93, assuming however no temperature dependence. Recently, Gieren et al. (2015) studied another Cepheid in the LMC (P = 2.99 d) and found a similar rather strong limb-darkening effect. If we consider the lower k value found by Pilecki et al. (2013), all the E1E2 measurements in Fig. 3 should translate towards lower values of x − 4% or about 0.06. The measurements in blue crosses in the figure would be better fitted, while green open circles would not (in particular near x ≃ 1.7). The limb-darkening effect is probably not the key issue to explain the discrepancy found for the E1E2 measurements.

5. Discussion

We observed δ Cep intensively on twenty nights. The initial purpose of these observations was to apply the BW method, but instead we detected for the very first time a static resolved structure around δ Cep contributing to about 7% of the total flux in the visible. This is not totally surprising as envelopes around Cepheids have already been discovered by long-baseline interferometry in the K band with VLTI and CHARA (Kervella et al. 2006; Mérand et al. 2006). In addition, four Cepheids have also been observed in the N band with VISIR and MIDI (Kervella et al. 2009; Gallenne et al. 2013) and one with NACO (Gallenne et al. 2011, 2012). Some evidence has also been found using high-resolution spectroscopy (Nardetto et al. 2008). From these observations, the typical size of the envelope of Cepheids seems to be around 3 stellar radii and the flux contribution from 2% to 10% of the continuum in the K band for medium- and long-period Cepheids, respectively, while it is around 10% or more in the N band. However, Mérand et al. (2005) do not mention any contribution from a CSE in the case of δ Cep in their infrared FLUOR data, while Mérand et al. (2006) found an improved agreement with a larger set of data when considering a model with a CSE. They used a ring model 3.54 mas in size, 0.5 mas in width, and with a relative brightness of 1.5%, which is confirmed in Merand et al. (2015a). The processes at work in infrared and in the visible regarding the CSE are different. We expect thermal emission in the infrared and scattering in the visible. It is also worth mentioning that the mid-infrared flux in excess is not necessarily the evidence for mass loss (see e.g. Schmidt 2015) and that the resolved structure we observe around δ Cep, instead of a CSE, could be simply a visible nebulae (also contributing to 7% of the total flux), as reported in the infrared by Marengo et al. (2010) and in the radio domain by Matthews et al. (2012).

Our second result is the presence of an additional second-order discrepancy between the observations and the models (pulsating disk + CSE or pulsating disk + background) at high spatial frequencies, which is not clearly understood. One possibility could be that the star is lighting up its environment differently at minimum (φ ≃ 0.8−1.0) and maximum (φ ≃ 0.0−0.8) radius. This reverberation effect would then be more important in the visible than in the infrared since the contribution in flux of the CSE (or the background) in the visible band is about 7% compared to 1.5% in the infrared. Interestingly, the phase interval from 0.8 to 1.0 is usually disregarded in the infrared surface brightness version of the BW method owing to general poor agreement between spectroscopic and photometric angular diameters (Storm et al. 2011a,b). The other possibility is to consider a static environment with a constant brightness in time, but with an additional dynamical effect. However, from the hydrodynamical simulations and observations of Cepheids in eclipsing binary systems, it seems that the limb-darkening effect in the visible band does not vary strongly enough to reproduce the second-order discrepancy found in our interferometric observations. Also possible is a strong compression or shockwave occurring near the photosphere at the minimum radius, which is indeed the layer probed by photometry and interferometry. However, such a shock is not seen in the hydrodynamical code (Nardetto et al. 2004) or in the spectroscopic data, conversely to X Sgr, an atypical Cepheid in which a shockwave has indeed been reported (Mathias et al. 2006). Why this dynamical effect is seen in the visible band and not in the infrared is another pending question in this hypothesis. We also exclude non-radial pulsation, since very accurate space photometry does not support the detection of such modes in Galactic classical Cepheids (e.g. Poretti et al. 2015), and the previous spectroscopic campaigns would have shown up quite easily. Furthermore, Anderson et al. (2015) found the signature of a companion in recent spectroscopic data of δ Cep. The expected flux contribution is supposed to be lower than 1% in the visible, clearly not detectable by VEGA/CHARA. Finally, reverberation seems to be the more plausible hypothesis. More interferometric data on Cepheids in the visible might shed light on it. In particular, covering the spatial frequencies properly (x from 0.2 to 3) at a unique pulsating phase (maximum and minimum radii) is probably required to go further in the analysis.


Acknowledgments

The authors acknowledge the support of the French Agence Nationale de la Recherche (ANR), under grant ANR-15-CE31-0012- 01 (project UnlockCepheids). We acknowledge financial support from “Programme National de Physique Stellaire” (PNPS) of CNRS/INSU, France. The CHARA Array is funded by the National Science Foundation through NSF grants AST-0606958 and AST-0908253 and by Georgia State University through the College of Arts and Sciences, as well as the W. M. Keck Foundation. W.G. gratefully acknowledges financial support for this work from the BASAL Centro de Astrofisica y Tecnologias Afines (CATA) PFB-06/2007, and from the Millenium Institute of Astrophysics (MAS) of the Iniciativa Cientifica Milenio del Ministerio de Economia, Fomento y Turismo de Chile, project IC120009. We acknowledge financial support for this work from ECOS-CONICYT grant C13U01. Support from the Polish National Science Center grant MAESTRO 2012/06/A/ST9/00269 is also acknowledged. EP and MR acknowledge financial support from PRIN INAF-2014. N.N., P.K., A.G., and W.G. acknowledge support from the French-Chilean exchange program ECOS- Sud/CONICYT (C13U01). This project was partially supported by the Polish Ministry of Science grant Ideas Plus. This research has made use of the SIMBAD and VIZIER databases at CDS, Strasbourg, France (http://cdsweb.u- strasbg.fr/), the Jean-Marie Mariotti Center Aspro service (http://www.jmmc.fr/aspro), and the electronic bibliography maintained by the NASA/ADS system. This research has made use of the Jean-Marie Mariotti Center SearchCal service (http://www.jmmc.fr/searchcal) co-developed by FIZEAU and LAOG/IPAG, and CDS Astronomical Databases SIMBAD and VIZIER (http://cdsweb.u-strasbg.fr/). This research has made use of the Jean-Marie Mariotti Center LITpro service co-developed by CRAL, IPAG, and LAGRANGE (LITpro software available at http://www.jmmc.fr/litpro). This research has made use of the Jean-Marie Mariotti Center OiDB service (http://oidb.jmmc.fr).

References

  1. Anderson, R. I. 2014, A&A, 566, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Anderson, R. I., Sahlmann, J., Holl, B., et al. 2015, ApJ, 804, 144 [NASA ADS] [CrossRef] [Google Scholar]
  3. Anderson, R. I., Mérand, A., Kervella, P., et al. 2016, MNRAS, 455, 4231 [NASA ADS] [CrossRef] [Google Scholar]
  4. Armstrong, J. T., Nordgren, T. E., Germain, M. E., et al. 2001, AJ, 121, 476 [NASA ADS] [CrossRef] [Google Scholar]
  5. Baade, W. 1926, Astron. Nachr., 228, 359 [NASA ADS] [CrossRef] [Google Scholar]
  6. Barnes, III, T. G., Fernley, J. A., Frueh, M. L., et al. 1997, PASP, 109, 645 [NASA ADS] [CrossRef] [Google Scholar]
  7. Benedict, G. F., McArthur, B. E., Fredrick, L. W., et al. 2002, AJ, 124, 1695 [NASA ADS] [CrossRef] [Google Scholar]
  8. Benedict, G. F., McArthur, B. E., Feast, M. W., et al. 2007, AJ, 133, 1810 [NASA ADS] [CrossRef] [Google Scholar]
  9. Berdnikov, L. N., & Turner, D. G. 2002, VizieR Online Data Catalog J/ApJS/137/209 [Google Scholar]
  10. Bersier, D., Burki, G., Mayor, M., & Duquennoy, A. 1994, A&AS, 108, 25 [NASA ADS] [Google Scholar]
  11. Bonneau, D., Clausse, J.-M., Delfosse, X., et al. 2006, A&A, 456, 789 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bonneau, D., Delfosse, X., Mourard, D., et al. 2011, A&A, 535, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Breitfelder, J., Kervella, P., Mérand, A., et al. 2015, A&A, 576, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Breitfelder, J., Mérand, A., Kervella, P., et al. 2016, A&A, 587, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Davis, J., Jacob, A. P., Robertson, J. G., et al. 2009, MNRAS, 394, 1620 [NASA ADS] [CrossRef] [Google Scholar]
  16. Engle, S. G., Guinan, E. F., Harper, G. M., Neilson, H. R., & Remage Evans, N. 2014, ApJ, 794, 80 [NASA ADS] [CrossRef] [Google Scholar]
  17. Fouqué, P., Arriagada, P., Storm, J., et al. 2007, A&A, 476, 73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Gallenne, A., Mérand, A., Kervella, P., & Girard, J. H. V. 2011, A&A, 527, A51 [Google Scholar]
  19. Gallenne, A., Kervella, P., & Mérand, A. 2012, A&A, 538, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Gallenne, A., Mérand, A., Kervella, P., et al. 2013, A&A, 558, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Gieren, W., Pilecki, B., Pietrzyński, G., et al. 2015, ApJ, 815, 28 [NASA ADS] [CrossRef] [Google Scholar]
  22. Kervella, P., Coudé du Foresto, V., Traub, W. A., & Lacasse, M. G. 1999, in Working on the Fringe: Optical and IR Interferometry from Ground and Space, eds. S. Unwin, & R. Stachnik, ASP Conf. Ser., 194, 22 [Google Scholar]
  23. Kervella, P., Coudé du Foresto, V., Perrin, G., et al. 2001, A&A, 367, 876 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Kervella, P., Nardetto, N., Bersier, D., Mourard, D., & Coudé du Foresto, V. 2004, A&A, 416, 941 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Kervella, P., Mérand, A., Perrin, G., & Coudé du Foresto, V. 2006, A&A, 448, 623 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Kervella, P., Mérand, A., & Gallenne, A. 2009, A&A, 498, 425 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Kiss, L. L. 1998, in A Half Century of Stellar Pulsation Interpretation, eds. P. A. Bradley, & J. A. Guzik, ASP Conf. Ser., 135, 173 [Google Scholar]
  28. Kukarkin, B. V., Kholopov, P. N., Pskovsky, Y. P., et al. 1971, in General Catalogue of Variable Stars, 3rd edn. [Google Scholar]
  29. Lafrasse, S., Mella, G., Bonneau, D., et al. 2010, VizieR Online Data Catalog: II/300 [Google Scholar]
  30. Lane, B. F., Kuchner, M. J., Boden, A. F., Creech-Eakman, M., & Kulkarni, S. R. 2000, Nature, 407, 485 [NASA ADS] [CrossRef] [Google Scholar]
  31. Lane, B. F., Creech-Eakman, M. J., & Nordgren, T. E. 2002, ApJ, 573, 330 [NASA ADS] [CrossRef] [Google Scholar]
  32. Ligi, R., Mourard, D., Nardetto, N., & Clausse, J.-M. 2013, J. Astron. Instrum., 2, 40003 [Google Scholar]
  33. Lindemann, F. A. 1918, MNRAS, 78, 639 [NASA ADS] [CrossRef] [Google Scholar]
  34. Majaess, D., Turner, D., & Gieren, W. 2012, ApJ, 747, 145 [NASA ADS] [CrossRef] [Google Scholar]
  35. Marengo, M., Sasselov, D. D., Karovska, M., Papaliolios, C., & Armstrong, J. T. 2002, ApJ, 567, 1131 [NASA ADS] [CrossRef] [Google Scholar]
  36. Marengo, M., Karovska, M., Sasselov, D. D., et al. 2003, ApJ, 589, 968 [NASA ADS] [CrossRef] [Google Scholar]
  37. Marengo, M., Evans, N. R., Barmby, P., et al. 2010, ApJ, 725, 2392 [NASA ADS] [CrossRef] [Google Scholar]
  38. Mathias, P., Gillet, D., Fokin, A. B., et al. 2006, A&A, 457, 575 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Matthews, L. D., Marengo, M., Evans, N. R., & Bono, G. 2012, ApJ, 744, 53 [NASA ADS] [CrossRef] [Google Scholar]
  40. Mérand, A., Kervella, P., Coudé du Foresto, V., et al. 2005, A&A, 438, L9 [Google Scholar]
  41. Mérand, A., Kervella, P., Coudé du Foresto, V., et al. 2006, A&A, 453, 155 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Mérand, A., Aufdenberg, J. P., Kervella, P., et al. 2007, ApJ, 664, 1093 [NASA ADS] [CrossRef] [Google Scholar]
  43. Merand, A., Kervella, P., Breitfelder, J., et al. 2015a, A&A, 584, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Merand, A., Kervella, P., Breitfelder, J., et al. 2015b, VizieR Online Data Catalog: J/A+A/584/A80 [Google Scholar]
  45. Moffett, T. J., & Barnes, III, T. G. 1984, ApJS, 55, 389 [NASA ADS] [CrossRef] [Google Scholar]
  46. Mourard, D., Bonneau, D., Koechlin, L., et al. 1997, A&A, 317, 789 [Google Scholar]
  47. Mourard, D., Clausse, J. M., Marcotto, A., et al. 2009, A&A, 508, 1073 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Mourard, D., Bério, P., Perraut, K., et al. 2011, A&A, 531, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Mourard, D., Challouf, M., Ligi, R., et al. 2012, in SPIE Conf. Ser., 8445, 0 [Google Scholar]
  50. Nardetto, N., Fokin, A., Mourard, D., et al. 2004, A&A, 428, 131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Nardetto, N., Fokin, A., Mourard, D., & Mathias, P. 2006a, A&A, 454, 327 [Google Scholar]
  52. Nardetto, N., Mourard, D., Kervella, P., et al. 2006b, A&A, 453, 309 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Nardetto, N., Mourard, D., Mathias, P., Fokin, A., & Gillet, D. 2007, A&A, 471, 661 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Nardetto, N., Groh, J. H., Kraus, S., Millour, F., & Gillet, D. 2008, A&A, 489, 1263 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Nardetto, N., Gieren, W., Kervella, P., et al. 2009, A&A, 502, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Nardetto, N., Mourard, D., Tallon-Bosc, I., et al. 2011, A&A, 525, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Nordgren, T. E., Armstrong, J. T., Germain, M. E., et al. 2000, ApJ, 543, 972 [NASA ADS] [CrossRef] [Google Scholar]
  58. Pilecki, B., Graczyk, D., Pietrzyński, G., et al. 2013, MNRAS, 436, 953 [NASA ADS] [CrossRef] [Google Scholar]
  59. Poretti, E., Le Borgne, J. F., Rainer, M., et al. 2015, MNRAS, 454, 849 [NASA ADS] [CrossRef] [Google Scholar]
  60. Riess, A. G., Macri, L., Casertano, S., et al. 2011, ApJ, 730, 119 [NASA ADS] [CrossRef] [Google Scholar]
  61. Schmidt, E. G. 2015, ApJ, 813, 29 [NASA ADS] [CrossRef] [Google Scholar]
  62. Storm, J., Carney, B. W., Gieren, W. P., et al. 2004, A&A, 415, 531 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  63. Storm, J., Gieren, W., Fouqué, P., et al. 2011a, A&A, 534, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Storm, J., Gieren, W., Fouqué, P., et al. 2011b, A&A, 534, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Tallon-Bosc, I., Tallon, M., Thiébaut, E., et al. 2008, in SPIE Conf. Ser., 7013 [Google Scholar]
  66. ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 453 [NASA ADS] [CrossRef] [Google Scholar]
  67. Wesselink, A. J. 1946, Bull. Astron. Inst. Netherlands, 10, 91 [Google Scholar]

Appendix A: VEGA observations: Log, transfer functions, and visibility curves

Table A.1

Observing log with the E1E2 baseline.

Table A.2

Observing log with the S1S2 baseline.

thumbnail Fig. A.1

Two calibrators, C1 and C2, used to calibrate the VEGA/CHARA data of δ Cep are compared in terms of transfer functions with six other calibrators, C3 to C8 (listed in Table 2), on six different nights during 2008 and 2009. The consistency between the calibrators, which is shown by a linear fit (black line), is independent of the two-telescope baseline considered, S1S2 or W1W2.

thumbnail Fig. A.2

Transfer function (calculated either from the C1 or C2 calibrator) is shown together with the raw visibilities of the target (T) for four different nights in our sample, using different spectral configurations and two baselines, E1E2 and S1S2. For nights 2008 July 30 and 2014 October 24, the squared visibilities (or transfer functions) associated with calibrators C1 and C2 cannot be consistently compared as they correspond to different spectral configurations, but they are shown to illustrate the quality of the data. The 2014 October 24 observing night is the only one for which we have successive data from E1E2 (from around 3.5h to 6h in UTC) and S1S2 (from around 6.5h to 8h in UTC). For nights 2013 July 25 and 2013 December 17, the instrument configuration is the same over the entire observing time, and the stability of the transfer function is evident (black line).

thumbnail Fig. A.3

Observed calibrated squared visibilities (red dots) are plotted as a function of the spatial frequency for each pulsation phase (indicated in the lower left corner of each panel) together with the best fit of uniform disk model (solid blue line). In the upper right corner (inset panel) we show the corresponding (u,v) coverage.

thumbnail Fig. A.3

continued.

All Tables

Table 1

Cepheids for which an interferometric BW method has been applied.

Table 2

List and properties of calibration stars selected with the SearchCal software provided by the Jean-Marie Mariotti Center (JMMC) (Bonneau et al. 2006, 2011).

Table 3

Pulsation phases, number of visibility measurements (N), resulting uniform disk angular diameters, and reduced χ2.

Table 4

Angular diameters as derived from the application of the SPIPS algorithm to δ Cep, using photometry, interferometry, and the cross-correlated radial velocity curves as inputs (see Merand et al. 2015a).

Table A.1

Observing log with the E1E2 baseline.

Table A.2

Observing log with the S1S2 baseline.

All Figures

thumbnail Fig. 1

Uniform disk diameters θUD derived from the VEGA/CHARA E1E2 data plotted as a function of the pulsation phase (blue crosses). The S1S2 data (corresponding to four different pulsation phases, see Table 3) were merged and the corresponding mean θUD is represented by a black horizontal line together with the 1σ uncertainty (blue dashed zone). The uniform disk angular diameter curve from FLUOR/CHARA in the K band (Mérand et al. 2005) is also shown for comparison. The red crosses and the light blue squares show the K band and R band UD angular diameter curves, respectively, as predicted from the SPIPS code (Merand et al. 2015a).

In the text
thumbnail Fig. 2

a) S1S2 observed calibrated squared visibilities secured at four different pulsation phases plotted as a function of the spatial frequency together with the best fit of uniform disk model (solid blue line). In the upper right corner we show the corresponding (u,v) coverage. b) Same data fitted with a two-component model (open symbols) including a uniform pulsating disk and a CSE (see text). As a comparison, the solid lines are the corresponding one-component models, i.e. a uniform pulsating disk (in the R band) without the CSE.

In the text
thumbnail Fig. 3

Panel a), the VEGA squared visibilities plotted as a function of x=πBp[m]θUD[mas]λ[nm]\hbox{$x = \frac{\pi B_{\rm p} \mathrm{[m]} \theta_\mathrm{UD} \mathrm{[mas]} }{\lambda \mathrm{[nm]}}$} in order to allow a comparison between data of different phases (i.e. different θUD) and effective wavelengths. Five models are overplotted: (1) the UD pulsating disk without CSE from the SPIPS algorithm (blue solid line); (2) the SPIPS UD pulsating disk + the CSE from Eq. (1) (red dotted line); (3) the SPIPS UD pulsating disk + the CSE from Eq. (1) but two times brighter (dashed line); (4) the SPIPS UD pulsating disk + the CSE or Eq. (1) but two times fainter (dash-dotted line); (5) the SPIPS UD pulsating disk + a background contribution filling the field of view of the interferometer (magenta dotted line). The residual between the observations and models 1 to 4 are shown in panels b), c), d), and e), respectively. In panels d) and e) we show only the data corresponding to phases intervals φ = 0.0−0.8 and φ = 0.8−1.0, respectively.

In the text
thumbnail Fig. A.1

Two calibrators, C1 and C2, used to calibrate the VEGA/CHARA data of δ Cep are compared in terms of transfer functions with six other calibrators, C3 to C8 (listed in Table 2), on six different nights during 2008 and 2009. The consistency between the calibrators, which is shown by a linear fit (black line), is independent of the two-telescope baseline considered, S1S2 or W1W2.

In the text
thumbnail Fig. A.2

Transfer function (calculated either from the C1 or C2 calibrator) is shown together with the raw visibilities of the target (T) for four different nights in our sample, using different spectral configurations and two baselines, E1E2 and S1S2. For nights 2008 July 30 and 2014 October 24, the squared visibilities (or transfer functions) associated with calibrators C1 and C2 cannot be consistently compared as they correspond to different spectral configurations, but they are shown to illustrate the quality of the data. The 2014 October 24 observing night is the only one for which we have successive data from E1E2 (from around 3.5h to 6h in UTC) and S1S2 (from around 6.5h to 8h in UTC). For nights 2013 July 25 and 2013 December 17, the instrument configuration is the same over the entire observing time, and the stability of the transfer function is evident (black line).

In the text
thumbnail Fig. A.3

Observed calibrated squared visibilities (red dots) are plotted as a function of the spatial frequency for each pulsation phase (indicated in the lower left corner of each panel) together with the best fit of uniform disk model (solid blue line). In the upper right corner (inset panel) we show the corresponding (u,v) coverage.

In the text
thumbnail Fig. A.3

continued.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.