Multiplicity of Galactic Cepheids from long-baseline interferometry

Multiplicity of Galactic Cepheids from long-baseline interferometry

I. CHARA/MIRC detection of the companion of V1334 Cygni
A. Gallenne Universidad de Concepción, Departamento de Astronomía, Casilla 160-C, Concepción, Chile    J. D. Monnier Astronomy Department, University of Michigan, 1034 Dennison Bldg, Ann Arbor, MI 48109-1090, USA    A. Mérand European Southern Observatory, Alonso de Córdova 3107, Casilla 19001, Santiago 19, Chile    P. Kervella LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Université Paris Diderot, 5 Place Jules Janssen, F-92195 Meudon, France    S. Kraus Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, MS-78, Cambridge, MA 02138, USA    G. H. Schaefer The CHARA Array of Georgia State University, Mount Wilson CA 91023, USA    W. Gieren Universidad de Concepción, Departamento de Astronomía, Casilla 160-C, Concepción, Chile    G. Pietrzyński Universidad de Concepción, Departamento de Astronomía, Casilla 160-C, Concepción, Chile Warsaw University Observatory, Al. Ujazdowskie 4, 00-478, Warsaw, Poland    L. Szabados Konkoly Observatory, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences, H-1121 Budapest, Konkoly Thege Miklós út 15-17, Hungary    X. Che Astronomy Department, University of Michigan, 1034 Dennison Bldg, Ann Arbor, MI 48109-1090, USA    F. Baron Astronomy Department, University of Michigan, 1034 Dennison Bldg, Ann Arbor, MI 48109-1090, USA    E. Pedretti School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews, Fife, KY16 9SS, UK    H. McAlister The CHARA Array of Georgia State University, Mount Wilson CA 91023, USA    T. ten Brummelaar The CHARA Array of Georgia State University, Mount Wilson CA 91023, USA    J. Sturmann The CHARA Array of Georgia State University, Mount Wilson CA 91023, USA    L. Sturmann The CHARA Array of Georgia State University, Mount Wilson CA 91023, USA    N. Turner The CHARA Array of Georgia State University, Mount Wilson CA 91023, USA    C. Farrington The CHARA Array of Georgia State University, Mount Wilson CA 91023, USA    N. Vargas The CHARA Array of Georgia State University, Mount Wilson CA 91023, USA
Received January 11, 2013; accepted February 6, 2013
Key Words.:
techniques: interferometric – techniques: high angular resolution – stars: variables: Cepheids – star: binaries: close
offprints: A. Gallenne

Context:More than 60 % of Cepheids are in binary or multiple systems. Studying such systems could lead to a better understanding of the age and evolution of Cepheids. These are also useful tools to estimate the mass of Cepheids, and constrain theoretical models of their pulsation and evolution.

Aims:We aim at determining the masses of Cepheids in binary systems, as well as their geometric distances and the flux contribution of the companions. The combination of interferometry with spectroscopy will offer a unique and independent estimate of the Cepheid masses.

Methods:Using long-baseline interferometry at visible and infrared wavelengths, it is possible to spatially resolve binary systems containing a Cepheid down to milliarcsecond separations. Based on the resulting visual orbit and radial velocities, we can then derive the fundamental parameters of these systems, particularly the masses of the components and the geometric distance. We therefore performed interferometric observations of the first-overtone mode Cepheid V1334 Cyg with the CHARA/MIRC combiner.

Results:We report the first detection of a Cepheid companion using long-baseline interferometry. We detect the signature of a companion orbiting V1334 Cyg at two epochs. We measure a flux ratio between the companion and the Cepheid  %, giving an apparent magnitude  mag. The combination of interferometric and spectroscopic data have enabled the unique determination of the orbital elements.   days, ,  mas, , , , and . We derive a minimal distance  pc, a minimum mass for both stars of , with a spectral type earlier than B5.5V for the companion star. Our measured flux ratio suggests that radial velocity detection of the companion using spectroscopy is within reach, and would provide an orbital parallax and model-free masses.


1 Introduction

Classical Cepheid stars have been considered as reliable tools to estimate distances in the universe for more than a century (see e.g. Leavitt_1912_03_1; Fernie 1969; Sandage & Tammann 2006; Barnes 2009; Bono et al. 2010). Their Period-Luminosity (P–L) relation makes them valuable to determine extragalactic distances and to calibrate secondary distance indicators. In addition to the determination of distances, Cepheids are also powerful astrophysical laboratories that provide fundamental clues for studying the pulsation and evolution of intermediate-mass stars (see e.g. Prada Moroni et al. 2012; Bono et al. 2006; Caputo et al. 2005).

The occurrence of Cepheids in binary (multiple) systems seems to be as high as 60 % for the brightest Cepheids (Szabados 2003), and is often neglected in the determination of the P–L relation. The position of a Cepheid in the P–L diagram could be biased by the presence of close bright companions. If the difference in magnitude between the Cepheid and its companion is small, the apparent magnitude of the Cepheid will be overestimated. This has a particular impact on the use of a surface brightness (SB) method to estimate the radius and luminosity (e.g. Fouqué & Gieren 1997; Gieren et al. 1998). The radial velocity measurements can also be altered because of orbital effects. This leads to a bias in the distance estimate, since radial velocities are necessary in the Baade-Wesselink method to evaluate the radius. The knowledge of potential photometric and astrometric biases caused by the companions is therefore important for the distance scale.

Primary (Cepheid) Secondaryfrom Evans (2000). .
from Klagyivik & Szabados (2009). from the 2MASS catalogue (Cutri et al. 2003). from the 2MASS catalogue (Cutri et al. 2003). Sp. Typefrom Samus et al. (2009). from Samus et al. (2009). from the surface brightness relation of Kervella et al. (2004a) from the P–L relation of Bono et al. (2002). Sp. Type
(days) (mas) (pc) (days) (days) (AU) (rad) (M)
5.87 4.46 4.66 F2Ib 3.333 0.534 683 1937.5 2 443 607 0.197 2.46 3.95 0.529
111, , : mean apparent , and magnitudes. Sp. Type: spectral type. : period of pulsation. : mean angular diameter. : distance. : orbital period. : time passage through periastron. : eccentricity of the orbit. : projected semi-major axis of the orbit of the Cepheid about the center of mass of the system. : argument of periastron. : spectroscopic mass function.
Table 1: Parameters of the Cepheid V1334 Cyg and its close companion.

Cepheids belonging to binary systems also offer the unique opportunity to make progress in resolving the Cepheid mass problem. For many years, stellar evolutionary models have predicted Cepheid masses larger than those derived from pulsation models (Neilson et al. 2011; Keller 2008; Bono et al. 2006, 1999). To investigate the origin of this discrepancy, the combination of spectroscopic and interferometric measurements will allow us to derive orbital elements and dynamical masses. Recent analyses for two eclipsing binary Cepheids in the LMC were carried out by Pietrzyński et al. (2010, 2011) which yielded masses with an accuracy  %. These dynamical masses are in agreement with the ones calculated from stellar pulsation models, suggesting that the pulsational theory provides the true current masses (see also Prada Moroni et al. 2012). However, this first conclusion needs additional mass estimates to better constrain the two models, particularly for wide binaries for which no significant physical interaction between the stars is expected.

A number of Galactic Cepheids are known to have companions closer than 30″ (e.g. Evans et al. 2008; Remage-Evans_2011_12_0), but most of them are located too close to the Cepheid (1–20 mas) to be observed with a 10–meter class telescopes. The already existing orbit measurements were estimated only from IUE spectrum or from the radial velocity variations. The only actual way to spatially resolve such systems is to use long-baseline interferometry. We therefore have started a long-term interferometric observing program that aims at studying a sample of seven northern and southern binary Cepheids. The first goal is to determine the angular separation and the apparent brightness ratio from the interferometric visibility and closure phase measurements. Our long-term objective, which needs a good sampling of the orbital period to get a reliable fit (several years), is to determine the full set of orbital elements, absolute masses and geometric distances.

The binary (or multiple) Cepheid systems were selected according to two main criteria. Firstly, the angular separations (mainly estimated from spectroscopy so far) have to be resolvable by the existing long-baseline interferometers. Secondly, the contrast between the Cepheid and its companion should not be too large, to ensure the flux contribution of the companion is detectable in the data. We therefore selected systems with an angular separation  mas, and a dynamic range .

In this paper, we present the first results for the Cepheid V1334 Cyg, observed with the MIRC combiner in the band at the CHARA Array. The paper is organized as follows. We first introduce the observed Cepheid, with some knowledge about the primary and its companion. In Sect. 3, the instrument configurations and the data reduction are detailed. The data analysis, including the interferometric models used and the fitting steps, is discussed in Sect. 4. We then combined our interferometric results with spectroscopic data in Sect. 5 to derive the parameters of the V1334 Cyg system.

2 V1334 Cyg

V1334~Cyg (HD 203156, HR 8157) was the first Cepheid observed for our program. We present in Table 1 the known parameters of this system from the literature. The information about the Cepheid companion were derived from spectroscopic observations (Evans_1995_05_0; Evans 2000).

This short-period Cepheid is an interesting system because it has been studied for many years but the binary or triple nature of the system is still debated. It has been suspected to be a member of a visual binary system for decades (see e.g. Millis 1969; Abt & Levy 1970), with a separation between 0.1–0.2″, but it has not been spatially resolved so far. Early radial velocity measurements also showed strong evidence of a spectroscopic binary, but the orbital period could not be fully constrained. Abt & Levy (1970) found a period of  years for the visual component, but they also noticed a variation in the centre of mass velocity () of a shorter time-scale. The same variation was also observed by Szabados (1991), and is likely linked to a second closer companion. From the International Ultraviolet Exporer (IUE) low-resolution spectra, Evans_1995_05_0 detected the hottest star in the system, and derived the spectral type to be a B7.0V star. In a subsequent work, Evans (2000) concluded that the hottest star is the visual companion. From 30 years of radial velocity measurements, the same author solved for the orbital parameters of the companion and found a period of 1937 days (see Table 1), that is significantly shorter than the value derived by Abt & Levy (1970), strongly suggesting the presence of a third component. They also derived a projected semi-major axis for the orbit of the Cepheid around the center of mass  AU (see Table 1). Kiss & Vinkó (2000) also detected the change in , and hypothesized a yellow-bright close companion.

Many authors attempted to resolve the wide component, but we still do not have a firm conclusion because of its intermittent detection. For instance, Evans et al. (2006) set an upper limit of  mas using the Hubble Space Telescope in the far-ultraviolet (although depending on the brightness of the stars). Scardia et al. (2008) mentioned the detection of the companion from speckle interferometry in the band, and measured a separation of 160 mas, while most of previous speckle observations failed to resolve the system (from 1976 to 2005, see e.g. McAlister 1978; Hartkopf & McAlister 1984; Evans et al. 2006, see the last reference for a more detailled discussion), setting a separation  mas. As argued by Evans et al. (2006), these intermittent detections could be linked to a peculiar orbit, making the visual companion undetectable at some orbital phase, or it could even not exist.

3 Observations and data reduction

Figure 1: plane coverage of our observations of V1334 Cyg.

Our observations were performed using the Michigan InfraRed Combiner (MIRC), installed at the CHARA array (ten Brummelaar et al. 2005) located on Mount Wilson, California. The array consists of six 1 m aperture telescopes with an Y-shaped configuration (two telescopes on each branch), oriented to the east (E1, E2), west (W1,W2) and south (S1, S2), and so offering a good coverage of the plane. The baselines range from 34 m to 331 m, providing an angular resolution down to 0.5 mas in .

MIRC (Monnier et al. 2004, 2010) combines the light coming from all six telescopes in the or bands, with three spectral resolutions ( and 400). The recombination of six telescopes gives simultaneously 15 fringe visibilities and 20 closure phase measurements, that are our primary observables.

Our observations were carried out on UT 2012 July 27 and October 1 using an band filter (m) and either four, five or six telescopes. We used the low spectral resolution mode, where a prism splits the light on the detector into 8 narrow spectral channels. The plane coverage for these nights is presented in Fig. 1. We followed a standard observing procedure, i.e we monitored the interferometric transfer function by observing a calibrator before and after our Cepheids. The calibrators were selected using the SearchCal222Available at software (Bonneau et al. 2006, 2011) provided by the JMMC. The journal of the observations is presented in Table 2, and the corresponding calibrators are listed in Table 3.

UT Star Configuration
2012 Jul. 27 09:52 HD~200577 S1-S2-E1-E2-W2
2012 Jul. 27 10:23 V1334 Cyg S1-S2-E1-E2-W2
2012 Jul. 27 11:08 HD~214200 S1-S2-E1-E2-W2
2012 Jul. 27 11:44 V1334 Cyg S1-S2-E1-E2
2012 Oct. 01 02:50 HD~185395 S1-S2-E1-E2-W1-W2
2012 Oct. 01 03:40 V1334 Cyg S1-S2-E1-E2-W2
2012 Oct. 01 05:01 HD~199956 S1-S2-E1-E2-W1-W2
2012 Oct. 01 05:44 V1334 Cyg S1-S2-E1-E2-W1-W2
2012 Oct. 01 06:28 HD~218470 S1-S2-E1-E2-W1-W2
2012 Oct. 01 07:11 V1334 Cyg S1-S2-E1-E2
2012 Oct. 01 07:48 HD~207978 S1-S2-E1-E2
Table 2: Journal of the observations.
Calibrator Sp. Type
(HD) (mas) ()
200577 6.1 4.1 G8III 3.2
214200Redetermined from internal MIRC calibration with HD 200577. 6.1 4.2 K0 15.5
185395Redetermined from internal MIRC calibration with HD 200577. 4.5 3.7 F4V 18.2
199956 6.6 4.5 K0 7.1
218470 5.7 4.7 F5V 15.5
207978 5.5 4.4 F6IV 21.6
333: magnitudes in and bands. : uniform disk angular diameter in band. : angular distance to the Cepheid.
Table 3: Calibrators used for our observations.

The data were reduced with the standard MIRC pipeline (Monnier et al. 2007). The main procedure is to compute squared visibilities and triple products for each baseline and spectral channel, and to correct for photon and readout noises. A recent upgrade includes a simultaneous measurement of the photometric channels with the fringes, enhancing the accuracy of MIRC down to 3 % in visibilities (Che et al. 2010, 2012). For the July observations, we used a coherent integration of 75 ms to improve the signal-to-noise ratio of the closure phase, at some expense of the visibility calibration, but for October observations, we used a more standard 17 ms integration time. We applied the same calibration error model as described in Monnier et al. (2012). We then did an incoherent average of 15 min for the final data. We note that this average prevent us from detecting a periodic signal from a potential third component, and the detection of any incoherent light would then be limited by the uncertainty in the visibilities.

4 Model fitting

To model the squared visibilities, triple amplitude and closure phase signals, we used the LITpro444LITpro software available at model fitting software (Tallon-Bosc et al. 2008), based on the Levenberg-Marquardt algorithm. It provides a set of elementary models that can be combined all together. The software also contains a tool allowing the search for the global minimum to solve for the problem of multiple minima.

The two epoch observations (July and October) were reduced separately in order to detect the changing position of the companion.

4.1 The models

As a first step, the primary component (the Cepheid) was modeled with a uniform disk (UD) angular diameter. The complex visibility model is:

with the first-order Bessel function, , the spatial frequencies, the UD angular diameter, and the wavelength.

The choice of a UD diameter instead of a limb-darkened (LD) disk for the fitting procedure is justified because the angular diameters of the Cepheids are small compared to the angular resolution of the interferometer, and the limb darkening effects are therefore undetectable. The conversion from UD to LD angular diameter was done by using a linear-law parametrization , with the LD coefficient from Claret & Bloemen (2011). The conversion is then given by the approximation (Hanbury Brown et al. 1974):

It is worth mentioning that the uncertainty of the limb-darkened coefficient has a small impact on the angular diameter conversion. A variation of 20 % of gives a LD diameter difference of less than 0.5 %.

As a second step, we fitted to the data a model with one companion, assumed to be unresolved by the array. The corresponding model for this binary system is:

where is the flux ratio between the companion and the Cepheid, and is the complex visibility model of an unresolved source:

with the relative position of the companion w.r.t the Cepheid.

The closure phase is then estimated from the modulus and argument of the bispectrum, , for each closed baseline triangle and spectral channel:

We assumed that the variation of the angular diameter between different acquisitions in the same night is negligible compared to our level of accuracy.

4.2 Fitting steps

We performed a least-squares model fit simultaneously with the squared visibility, triple amplitude and closure phase measurements. Our search strategy was the following. We first proceeded to a grid search in the space to determine the approximate position and brightness ratio of the companion. A first grid search between  mas with a 0.2 mas spacing was performed with various flux ratio (from with 0.005 steps). Then a second grid of  mas with a 0.01 mas spacing around the most likely position and brightness ratio was used to refine the position. Finally all parameters were fitted using the refined values. The first guess for the angular diameter was taken from Table 1.

2012-07-27 2012-10-01
Single star model
1.63 2.08
Binary model
0.34 1.24
555 and are the uniform and limb-darkened disk angular diameter, respectively. , and correspond to the flux ratio and position of the companion. is the reduced of the corresponding fitted model.
Table 4: Summary of the parameters estimated from the model fit.

4.3 Results

Figure 2: Closure phase signal for the first epoch. The color-coded squares are the data, while the black dots represent the binary model.
Figure 3: Closure phase signal for the second epoch. The color-coded squares are the data, while the black dots represent the binary model.
Figure 4: Squared visibility measurements for the two epochs. The data are represented in blue, while the red dots are the fitted binary model.

The closure phase signal and squared visibilities, presented in Figs. 2, 3 and 4, clearly show a departure from a single star with a symmetric brightness distribution. These variations reveal the presence of at least one companion.

We chose as first guess  mas, and performed the fitting procedure with a grid search as explained in Sect. 4.2. The probability maps for our two epochs of observations are shown in Fig. 5, and the fitted parameters are reported in Table 4. The model is also represented graphically with black dots in Figs. 2, 3 and red dots in Fig. 4. We notice a good agreement between the model and the data. The fitted model give the most probable location of this companion at an angular separation  mas and a position angle for the first epoch, and  mas and for the second. The measured flux ratio is also particularly consistent between the two epochs.

We estimated the limb-darkened angular diameter to be  mas for the first epoch, using the limb-darkening coefficient (Claret & Bloemen 2011), chosen from the stellar parameters:  K, , [Fe/H] = 0.0 and  m s (Luck et al. 2008, at pulsation phase ). For the second epoch, the photospheric parameters did not change significantly (because it is low-amplitude pulsation Cepheid), leading to a same to find  mas (at pulsation phase ). These measurements are in agreement with the averaged angular diameters predicted from a surface brightness relation ( mas, Kervella et al. 2004b; Moskalik & Gorynya 2005), taking into account the small amplitude variation. As no IR photometric curves are available, it was not possible to compare our value with an angular diameter derived from a SB technique for the actual phase of our measurements.

Figure 5: Probability maps for the two observed epochs, showing the maximum obtained at each point of the search region with a 3 % flux ratio.

5 Combination of spectroscopy and interferometry

5.1 Preliminary orbit

Figure 6: Left: fitted (solid line) and measured (blue dots) orbital velocity. Middle: fitted (solid line) and measured (blue dots) pulsation velocity. Right: orbit of V1334 Cyg Ab. The data points are the MIRC results from Table 4.

Because of its intermittent detections, the existence of the wide component is rather uncertain, and might not exist. We claim that the companion detected with MIRC and IUE (Evans_1995_05_0) can be identified as the spectroscopic one. Therefore, our measured separations should be compatible with the spectroscopic orbit, and it should be possible to derive all the orbital parameters.

We therefore combined our astrometric measurements with the single-line radial velocity data gathered by Evans (2000), in order to derive the complete orbital elements for V1334 Cyg Ab. We stress that this analysis is preliminary as we have only two astrometric measurements so far. We computed the orbital solutions through the formalism developed by Wright & Howard (2009), slightly modified to include the pulsation velocity of the form:

where and are the fitted Fourier coefficients, and denotes the pulsation phase.

We solved for all the orbital elements () and pulsation parameters (). The reference epoch of maximum light for the Cepheid, , was held fixed to the value given by Samus et al. (2009). The initial values for the orbit and the pulsation were chosen from Evans (2000) and Samus et al. (2009). The final elements derived from our combined fit are listed in Table 5. The fitted parameters are in good agreement with those from Evans (2000). The quoted uncertainties for the elements derived from spectroscopy () were estimated using the bootstrapping technique (with replacement and 500 bootstrap samples). For the remaining elements derived from interferometry ( and ), we refitted the orbits 500 times, each time adding Gaussian noise to each astrometric point according to their uncertainties. The standard deviation from these trials is then used as the uncertainty. Fig. 6 (left and middle) shows the orbital and pulsation velocities disentangled from the radial velocity measurements. The solid black lines denote our fitted curves. The final best fit orbit of V1334 Cyg Ab is also plotted in Fig. 6 (right) with our MIRC measurements marked by the blue dot symbols.

Spectroscopy only This work
(Evans 2000)
From Samus et al. (2009), and held fixed when fitting (HJD)
666: orbital period. : time passage through periastron. : eccentricity. : radial velocity semi-amplitude of the primary. : systemic velocity. : argument of periastron. : position angle of the ascending node. : semi-major axis. : orbital inclination. : apparent magnitude in . : pulsation period. : reference epoch of maximum light. : Fourier parameters.
Table 5: Fitted orbital and pulsation parameters of V1334 Cyg Ab.

5.2 Apparent magnitude, spectral type and mass of the companion

Combining the -band magnitude given by the 2MASS catalog (Cutri et al. 2003) with our averaged measured flux ratio  %, we derived for the companion a magnitude , and  mag. As no -band light curve is available to estimate the Cepheid magnitude at our pulsation phase, an additional uncertainty of 3 % was quadratically added to take into account the phase mismatch. The choice of these 3 % is based on the amplitude variation of the light curve in (Klagyivik & Szabados 2009, that is surely to be lower in ).

The absolute magnitude, , can be estimated knowing the distance to the system. However, there is no accurate determination of the distance for this Cepheid. The Hipparcos data give a distance  pc (van Leeuwen 2007). The use of a -band P–L relation for first overtone (FO) pulsators (Bono et al. 2002, non-canonical model) gives  pc, while converting the overtone period to the fundamental one with the period ratio data from Alcock et al. (1995) and using a -band P–L relation for fundamental (F) mode pulsator (Storm et al. 2011) gives  pc. As there is no optimum value, we plotted in Fig. 7 (top panel) the spectral type vs. the distance. The previous cited distance ranges set a spectral type for the companion between a B8.0V and B4.0V star. The extinction was assumed negligible at our observing wavelength ( mag, estimated using the total-to-selective absorption ratios , from Fouqué et al. 2003, and the average color excess from Evans_1995_05_0 and Kovtyukh et al. 2008)

Figure 7: Absolute magnitude, spectral type of the companion (top panel), and mass of each component (bottom panel) as function of the distance. The shaded gray areas denote the distance range of Hipparcos, F mode P–L relation (left), and FO mode P–L relation (right). The vertical dotted line represents the minimum distance where .

We can also estimate the mass ratio, , for a range of distances by combining the Kepler’s third law with the spectroscopic mass function,


which yields:

with in arcseconds, in pc, in yr, and in . Combining Eq. 1 and 2 gives the mass of each component:


We plotted Eq. 3 and 4 vs. the distance in Fig. 7 (bottom panel). The shaded gray areas denote the distance range of Hipparcos, F mode P–L relation, and FO mode P–L relation. As we expect a lower mass companion (Evans_1995_05_0), we can set a lower limit to the distance, implies  pc, plotted as a dotted vertical line in Fig. 7. This limit is consistent with the distance derived from the FO P–L relation ( pc, Bono et al. 2002), and would mean a mass ratio equal to 1 for this system. With this limit, Eq. 3 yields a minimum mass for the companion (and so the Cepheid) , and a spectral type earlier than a B5.5V star.

6 Conclusion

We presented new multiple telescope interferometric observations of the classical Cepheid V1334 Cyg in the band. For the first time, we were able to spatially resolve the companion at two epochs. We derived the limb-darkened angular diameter for the Cepheid at their corresponding pulsation phase, the relative positions, and the flux ratio of the companion. We combined our accurate astrometric measurements with existing spectroscopic data and derived preliminary orbital solutions for the V1334 Cyg system. We also determined a minimal distance to the system to be  pc, which yields to a minimum mass for the component . We also found that its spectral type is earlier than a B5.5V star.

Our work, using multi-telescope recombination, provided unique and useful informations, both on the Cepheid and its companion. These innovative results show the capabilities of long-baseline interferometry to study Cepheids in binary or multiple systems. This is particularly important as most of the companions are located too close to the star to be spatially resolved by a single telescope. This technique provides new observables that can be efficiently combined with spectroscopic results to provide innovative constraints on the system properties. We plan to continue our observing program in the next years to increase the sample of binary Cepheids with well determined orbital elements.

The authors would like to thank the CHARA Array and Mount Wilson Observatory staff for their support. Research conducted at the CHARA Array is funded by the National Science Foundation through NSF grant AST-0908253, by Georgia State University, the W. M. Keck Foundation, the Packard Foundation, and the NASA Exoplanet Science Institute. JDM acknowledges funding from the NSF grants AST-0707927 and AST-0807577. WG and GP gratefully acknowledge financial support for this work from the BASAL Centro de Astrofísica y Tecnologías Afines (CATA) PFB-06/2007. Support from the Polish National Science Centre grant MAESTRO DEC-2012/06/A/ST9/00269 and the Polish Ministry of Science grant Ideas Plus (awarded to G. P.) is also acknowledge. This research received the support of PHASE, the high angular resolution partnership between ONERA, Observatoire de Paris, CNRS, and University Denis Diderot Paris 7. A. G. acknowledges support from FONDECYT grant 3130361. LSz has been supported by the ESTEC Contract No.4000106398/12/NL/KML. This work made use of the SIMBAD and VIZIER astrophysical database from CDS, Strasbourg, France and the bibliographic informations from the NASA Astrophysics Data System. This research has made use of the Jean-Marie Mariotti Center LITpro service, co-developed by CRAL, LAOG and FIZEAU, and SearchCal service, co-developed by FIZEAU and LAOG/IPAG, and of CDS Astronomical Databases SIMBAD and VIZIER.


  • Abt & Levy (1970) Abt, H. A. & Levy, S. G. 1970, PASP, 82, 334
  • Alcock et al. (1995) Alcock, C., Allsman, R. A., Axelrod, T. S., et al. 1995, AJ, 109, 1653
  • Barnes (2009) Barnes, T. G. 2009, in AIPC Series, ed. J. A. Guzik & P. A. Bradley, Vol. 1170, 3–12
  • Bonneau et al. (2006) Bonneau, D., Clausse, J.-M., Delfosse, X., et al. 2006, A&A, 456, 789
  • Bonneau et al. (2011) Bonneau, D., Delfosse, X., Mourard, D., et al. 2011, A&A, 535, A53
  • Bono et al. (2006) Bono, G., Caputo, F., & Castellani, V. 2006, Mem. Soc. Astron. Italiana, 77, 207
  • Bono et al. (2010) Bono, G., Caputo, F., Marconi, M., & Musella, I. 2010, ApJ, 715, 277
  • Bono et al. (2002) Bono, G., Groenewegen, M. A. T., Marconi, M., & Caputo, F. 2002, ApJ, 574, L33
  • Bono et al. (1999) Bono, G., Marconi, M., & Stellingwerf, R. F. 1999, ApJS, 122, 167
  • Caputo et al. (2005) Caputo, F., Bono, G., Fiorentino, G., Marconi, M., & Musella, I. 2005, ApJ, 629, 1021
  • Che et al. (2012) Che, X., Monnier, J. D., Kraus, S., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. F. Delplancke, J. K. Rajagopal, & F. Malbet, Vol. 8445
  • Che et al. (2010) Che, X., Monnier, J. D., & Webster, S. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. W. C. Danchi, F. Delplancke, & J. K. Rajagopal, Vol. 7734, 77342V–77342V–9
  • Claret & Bloemen (2011) Claret, A. & Bloemen, S. 2011, A&A, 529, A75
  • Cutri et al. (2003) Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog, 2246, 0
  • Evans (2000) Evans, N. R. 2000, AJ, 119, 3050
  • Evans et al. (2006) Evans, N. R., Franz, O., Massa, D., et al. 2006, PASP, 118, 1545
  • Evans et al. (2008) Evans, N. R., Schaefer, G. H., Bond, H. E., et al. 2008, AJ, 136, 1137
  • Fernie (1969) Fernie, J. D. 1969, PASP, 81, 707
  • Fouqué & Gieren (1997) Fouqué, P. & Gieren, W. P. 1997, A&A, 320, 799
  • Fouqué et al. (2003) Fouqué, P., Storm, J., & Gieren, W. 2003, in Lecture Notes in Physics, Vol. 635, Stellar Candles for the Extragalactic Distance Scale, ed. D. Alloin & W. Gieren (Berlin Springer), 21
  • Gieren et al. (1998) Gieren, W. P., Fouque, P., & Gomez, M. 1998, ApJ, 496, 17
  • Hanbury Brown et al. (1974) Hanbury Brown, R., Davis, J., Lake, R. J. W., & Thompson, R. J. 1974, MNRAS, 167, 475
  • Hartkopf & McAlister (1984) Hartkopf, W. I. & McAlister, H. A. 1984, PASP, 96, 105
  • Keller (2008) Keller, S. C. 2008, ApJ, 483
  • Kervella et al. (2004a) Kervella, P., Bersier, D., Mourard, D., et al. 2004a, A&A, 428, 587
  • Kervella et al. (2004b) Kervella, P., Thévenin, F., Di Folco, E., & Ségransan, D. 2004b, A&A, 426, 297
  • Kiss & Vinkó (2000) Kiss, L. L. & Vinkó, J. 2000, MNRAS, 314, 420
  • Klagyivik & Szabados (2009) Klagyivik, P. & Szabados, L. 2009, A&A, 504, 959
  • Kovtyukh et al. (2008) Kovtyukh, V. V., Soubiran, C., Luck, R. E., et al. 2008, MNRAS, 389, 1336
  • Luck et al. (2008) Luck, R. E., Andrievsky, S. M., Fokin, A., & Kovtyukh, V. V. 2008, AJ, 136, 98
  • McAlister (1978) McAlister, H. A. 1978, PASP, 90, 288
  • Millis (1969) Millis, R. L. 1969, Lowell Observatory Bulletin, 7, 113
  • Monnier et al. (2010) Monnier, J. D., Anderson, M., Baron, F., et al. 2010, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. W. C. Danchi, F. Delplancke, & J. K. Rajagopal, Vol. 7734, 77340G–77340G–12
  • Monnier et al. (2004) Monnier, J. D., Berger, J.-P., Millan-Gabet, R., & ten Brummelaar, T. A. 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. W. A. Traub, Vol. 5491, 1370
  • Monnier et al. (2012) Monnier, J. D., Che, X., Zhao, M., et al. 2012, ApJ, 761, L3
  • Monnier et al. (2007) Monnier, J. D., Zhao, M., Pedretti, E., et al. 2007, Science, 317, 342
  • Moskalik & Gorynya (2005) Moskalik, P. & Gorynya, N. A. 2005, Acta Astronomica, 55, 247
  • Neilson et al. (2011) Neilson, H. R., Cantiello, M., & Langer, N. 2011, A&A, 529, L9
  • Pietrzyński et al. (2010) Pietrzyński, G., Thompson, I. B., Gieren, W., et al. 2010, Nature, 468, 542
  • Pietrzyński et al. (2011) Pietrzyński, G., Thompson, I. B., Graczyk, D., et al. 2011, ApJ, 742, L20
  • Prada Moroni et al. (2012) Prada Moroni, P. G., Gennaro, M., Bono, G., et al. 2012, ApJ, 749, 108
  • Samus et al. (2009) Samus, N. N., Durlevich, O. V., et al. 2009, VizieR Online Data Catalog: B/gcvs. Originally published in: Institute of Astronomy of Russian Academy of Sciences and Sternberg State Astronomical Institute of the Moscow State University, 1, 2025
  • Sandage & Tammann (2006) Sandage, A. & Tammann, G. A. 2006, ARA&A, 44, 93
  • Scardia et al. (2008) Scardia, M., Prieur, J.-L., Pansecchi, L., et al. 2008, Astronomische Nachrichten, 329, 54
  • Storm et al. (2011) Storm, J., Gieren, W., Fouqué, P., et al. 2011, A&A, 534, A94
  • Szabados (1991) Szabados, L. 1991, Commun. of the Konkoly Observatory Hungary, 96, 123
  • Szabados (2003) Szabados, L. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 298, GAIA Spectroscopy: Science and Technology, ed. U. Munari, 237
  • Tallon-Bosc et al. (2008) Tallon-Bosc, I., Tallon, M., Thiébaut, E., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7013, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series
  • ten Brummelaar et al. (2005) ten Brummelaar, T. A., McAlister, H. A., Ridgway, S. T., et al. 2005, ApJ, 628, 453
  • van Leeuwen (2007) van Leeuwen, F. 2007, A&A, 474, 653
  • Wright & Howard (2009) Wright, J. T. & Howard, A. W. 2009, ApJS, 182, 205
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description