On the binarity of the classical Cepheid X Sgr from interferometric observations
Key Words.:Stars: variables: Cepheids – Binaries: close – Techniques: interferometric – Instrumentation: interferometers – Infrared: stars
Context:Optical-infrared interferometry can provide direct geometrical measurements of the radii of Cepheids and/or reveal unknown binary companions of these stars. Such information is of great importance for a proper calibration of Period-Luminosity relations and for determining binary fraction among Cepheids.
Aims:We observed the Cepheid X Sgr with VLTI/AMBER in order to confirm or disprove the presence of the hypothesized binary companion and to directly measure the mean stellar radius, possibly detecting its variation along the pulsation cycle.
Methods:From AMBER observations in MR mode we performed a binary model fitting on the closure phase and a limb-darkened model fitting on the visibility.
Results:Our analysis indicates the presence of a point-like companion at a separation of 10.7 mas and 5.6 mag fainter than the primary, whose flux and position are sharply constrained by the data. The radius pulsation is not detected, whereas the average limb-darkened diameter results to be 1.480.08 mas, corresponding to 533 R. at a distance of 333.3 pc.
Classical Cepheids are fundamental astrophysical objects, being the most popular primary distance indicators in the nearby Universe. These pulsating stars obey well defined Period-Luminosity (PL) relations, so that it is possible to derive their distance moduli by measuring periods and apparent magnitudes.
Cepheids can be used to trace the evolutionary properties of Helium burning intermediate-mass stars, as they are crossing the instability strip along the so-called “blue loops” in the HR diagram. However, the comparison between pulsation and evolutionary masses disclosed the problem known as “Cepheid mass discrepancy” (Cox 1980), due to the pulsation masses resulting systematically smaller than the evolutionary ones (see e.g. Bono et al. (2001); Beaulieu et al. (2001); Caputo et al. (2005); Keller & Wood (2006)).
Recently two double eclipse binary Cepheids have been identified in the Large Magellanic Cloud (Pietrzyński et al. (2011), Pietrzyński et al. (2011)), providing the first measurements of the dynamical mass of classical Cepheids with 1-3% accuracy, whose values support the pulsation masses, or the Mass-Luminosity relations based on evolutionary models that account for mild convective core overshooting (Cassisi & Salaris (2011); Prada Moroni et al. (2012)).
For non-eclipsing binaries, direct orbital determination of the companion by means of optical interferometry would provide dynamical mass determinations able to disentangle the mass discrepancy problem.
The binary fraction of the Cepheids has an impact on the Cepheid distance scale in the optical bands, since the companions are typically main sequence (MS) stars. The binarity analysis has been applied to a limited sample (Szabados 2003), but current estimates indicate that the binary fraction among MS B-type stars is of the order of 60–70% (Brott et al. (2011); Chini et al. (2012)).
In order to improve the current observational scenario we have observed the Cepheid X Sgr with the VLTI. This particular target was selected because: i) its trigonometric parallax has been recently measured by Benedict et al. (2007); ii) accurate radial velocity curves are available (Mathias et al. (2006); Storm et al. (2011)), and iii) it was suggested by Szabados (1990), on the basis of temporal changes in the -velocity, that X Sgr might be a binary, although this finding has always been controversial (Mathias et al. (2006); Groenewegen (2008)).
2 Observations and data reduction
2.1 Data acquisition
We observed X Sgr (HD161592, R.A. 17:47:33.62, Dec. -27:49:50.83, J2000) using VLTI/AMBER (Petrov et al. 2007) with the Auxiliary Telescopes. We used medium spectral resolution in the -band (MR, R1500, 2.1m) and the fringe tracker FINITO (Gai et al. 2004) with a frame exposure time of 1 s. Observations were carried out with a maximum projected baseline length of m, corresponding to mas spatial resolution, with the atmospheric seeing varied in the range.
|Date||Target sequence||Baselines||of frames for||of frames for|
Data from JMMC Stellar Diameters Catalogue - JSDC (Lafrasse et al. 2010)
2.2 Data reduction
We used the amdlib v3
Each AMBER file contains many individual frames, which must be quality-selected to reject tracking losses in both FINITO and telescope guiding. Such events are associated to large values in the optical path difference (OPD) and cause a poor signal to noise ratio (SNR) of the fringes. We adopted the amdlib “PISTON+SNR” selection criteria, limiting the OPD variation below 8m and keeping only the frames with the highest fringe SNR. The adopted percentage of frames, in the range for visibility and for closure phase , has been optimized for each night separately (see Table 1).
Hence, we computed, for each file, the average values of and , weighted by the error estimate provided by amdlib; the same was done for the calibrators.
Finally, a night transfer function (TF) was computed using the calibrators listed in Table 1.
3 Results and Discussion
3.1 Visibilities and closure phases
We see that the source is resolved, because is less than unity, and that significant deviations from a pure Uniform Disk (UD) are evident. This could indicate a possible detection of the Cepheid pulsation, or might suggest a more complex morphology, like that of a binary. In any case, in the plot we only see small deviations from zero, which is an indication that the companion, if present, should be faint and close to our detection limit.
Both plots, in particular the , show that the error bars are clearly larger than the dispersion of the data points, indicating a correlated noise caused by the low-frequency variations (on time-scales of minutes) of the TF during the night.
3.2 Model fitting
Szabados (1990) suggested, on the basis of changes in the -velocity, that X Sgr might be a binary system with a period of 507 days. The complex nature of X Sgr was further supported by Sasselov & Lester (1990) who spectrally found that line asymmetries and line-doubling were well beyond the typical Cepheids behaviour; in particular, the two absorption components showed different Doppler shifts, but they did not display secular velocity changes compared with the -velocity. This complex dynamical structure was also supported by other high-resolution (R120,000) optical spectra by Mathias et al. (2006). They hypothesized a binary, or even triple, nature of X Sgr. The data available in the literature were analyzed by Feast et al. (2008) and by Groenewegen (2008) who also found evidence of binarity, these latter suggesting a slightly longer orbital period (570 days).
In order to check the companion hypothesis, we performed a binary model fit to the closure phase . We decided to use only this observable, for its low sensitivity on the unknown radius variation, for which we adopted a crude UD estimate of mas, provided by the values at the largest baselines, namely the data from nights 2011-04-24, 2011-05-03 and 2012-03-31.
We started with Run A alone assuming no temporal variation in the binary configuration within the two months, which is a plausible assumption, considering the mentioned orbital period estimates. Therefore we fit the Run A data all together in order to increase accuracy. Run B was not included here, since it was acquired 10 months later.
We adopted a grid search on a wide parameter space in and , where are the R.A. and Dec. offsets of the companion (these are searched within a radius of 100 mas around the primary, see Absil et al. (2011)), is the flux ratio of the companion to the primary, and is the companion UD radius.
We obtain the best fit result for a flux ratio of (corresponding to a magnitude difference), while tends to zero, thus confirming the presence of a faint and unresolved companion, whose position is constrained by a single and well defined minimum at 7.15 S, 8.00 E mas from the primary (see Fig 3). The parameters of the companion are summarized in Table 2.
The fact that the best fit reduced is significantly smaller than one is due to the large error bars associated to the TF, already described in Sec. 3.1.
To double check our finding, we built a series of synthetic binary models with different flux ratios and the same error bars of the real data, from which we derived a sensitivity limit of , i.e. 3 times lower than the measured flux ratio.
To refine this result we removed our assumption of fixed binary configuration, in order to check if we could detect possible displacements among different nights.
Therefore, we repeated the fitting separately for each night of Run A, to find out if the location of the minimum moves as a function of time (nights 2011-05-03/2011-05-07 and 2011-06-11/2011-06-25 were fitted as single points).
The corresponding maps (Fig. 3, panel A) show several local minima within the 99 confidence interval, all of them having the same probability to be the true position of the companion. This shows us that a single baseline configuration is not sufficient to unambiguously locate such a faint object.
However, each map shows a particular local minimum in the same position found using the entire data set, and even more importantly the little displacements of these minima are in the correct temporal sequence, as disclosed by Fig. 3, panel B.
This empirical evidence further confirms our detection of the postulated companion of X Sgr, and although current data do not allow us to constrain the orbit, we have a hint on the direction of its motion.
The binary fit for the single night of Run B also show the same kind of multiple minima map (not displayed), but with no minima close to the quoted position, as expected after a time interval of ten months.
|8.000.02 E||7.150.04 S||0.5||(61)|
The pulsation phase was estimated using the formula provided by Szabados (1989) for constraining the pulsation period and the epoch of maximum of the different data sets.
For the 2012-03-31 night a pure LD model is adopted, given the undetermined binary parameters.
3.3 Angular diameter
To determine the mean Cepheid diameter, and possibly its variation, we fit the separately for each night by fixing the binary parameters (Table 2) and varying only the primary radius . Following Kervella et al. (2004), we adopted the Limb Darkening (LD) intensity profile by Claret (2000).
The radius only depends on the mean
value of each visibility spectrum (), whose error bar we computed as
the quadratic sum of the weighted average error plus the correlated component due to the TF variation
For the diameter determination we only took into account the nights 2011-04-24, 2011-05-03 and 2012-03-31, which have sufficiently long baselines and good quality. The night of 2011-05-03 plays a key role in this context, since it shows a pure random error with a negligible correlated component.
In the table we distinguish the random errors and the systematics coming from the uncertainties in the calibrator diameters (Table 1). Note that the latter are larger than the former only for our best data set, namely the night 2011-05-03, showing the importance of using many different calibrators when a very accurate diameter measurement is required.
|This paper||53 3|
|This paper||55 4|
|Kervella et al. (2004)||53 3|
|Storm et al. (2011)||50 3|
Solution obtained by neglecting the presence of the companion.
The differences in the three diameter estimates, within from each other, do not allow us to positively detect the radius pulsation, and are probably due to the mentioned TF temporal variation. The same limitation was also found by Kervella et al. (2004) using interferometric data collected with VINCI at VLTI.
Thus we only computed a mean physical radius of 533 R by adopting the trigonometric parallax of 3.00.18 mas from the HST (Benedict et al. 2007). In Table 3 we compared it with recent estimates from both optical interferometry (Kervella et al. 2004) and Infra Red Surface Brightness (IRSB, Groenewegen (personal communication) and Storm et al. (2011)), having used the same distance for all of these measurements. Finally, we also report the radius obtained by neglecting the binary companion, which has a 3% impact on the mean radius.
We found evidence of the binary nature of the Cepheid X Sgr by analysing new interferometric data collected with AMBER at VLTI. The main conclusions of our work can be summarized as follows:
i) Our analysis of the closure phase indicates that both the flux and the position of the companion can be constrained by our data sets. We find that the angular separation is 10.7 mas and that the companion is 5.6 mag fainter than the primary in the K-band.
ii) The data allow us to directly measure an average Limb-Darkening diameter of 1.480.08 mas (i.e. a physical radius of 533 R) which agrees, at 1 level, with similar interferometric and IRSB mean radii for X Sgr available in the literature.
iii) We find that MR AMBER closure phases obtained from 3 baseline configurations have the sensitivity needed to detect companions around nearby classical Cepheids. This evidence becomes even more compelling if we take into account that current data were collected in medium quality seeing conditions.
iv) The average visibilities of the 2001-05-03 night show that a non systematic error as small as 1% can be obtained with AMBER in MR mode in stable TF conditions. This would be sufficient to clearly trace the differential diameter pulsation in future measurements, provided that the same calibrator star is used for all the observations.
v) We find that, in the case of a binary Cepheid, the effect introduced by the companion on the primary radius measurement (3%) is comparable with the expected radius variation (5 10%).
Acknowledgements.We are grateful to the Nice Observatory (France), and in particular to F. Millour, N. Nardetto and P. Cruzalebes for helpful discussions, and to M.A.T. Gronewegen and J. Storm for many useful insights concerning Cepheids and for providing us with their recent radius measurements of X Sgr.
- thanks: Based on observations collected at the European Southern Observatory, Paranal, Chile
- offprints: Gianluca Li Causi,
- Available at http://www.jmmc.fr/amberdrs
- where the latter was estimated, for each visibility spectrum, as the squared difference between the r.m.s. of the error bars and the r.m.s. of the residuals with respect to a smooth polynomial fitting.
- Absil, O., Le Bouquin, J.-B., Berger, J.-P., et al. 2011, A&A, 535, A68
- Beaulieu, J. P., Buchler, J. R., & Kolláth, Z. 2001, A&A, 373, 164
- Benedict, G. F., McArthur, B. E., Feast, M. W., et al. 2007, AJ, 133, 1810
- Bono, G., Gieren, W. P., Marconi, M., Fouqué, P., & Caputo, F. 2001, ApJ, 563, 319
- Brott, I., Evans, C. J., Hunter, I., et al. 2011, A&A, 530, A116
- Caputo, F., Bono, G., Fiorentino, G., Marconi, M., & Musella, I. 2005, ApJ, 629, 1021
- Cassisi, S. & Salaris, M. 2011, ApJ, 728, L43
- Chelli, A., Utrera, O. H., & Duvert, G. 2009, A&A, 502, 705
- Chini, R., Hoffmeister, V. H., Nasseri, A., Stahl, O., & Zinnecker, H. 2012, MNRAS, 424, 1925
- Claret, A. 2000, A&A, 363, 1081
- Cox, A. N. 1980, ARA&A, 18, 15
- Feast, M. W., Laney, C. D., Kinman, T. D., van Leeuwen, F., & Whitelock, P. A. 2008, MNRAS, 386, 2115
- Gai, M., Menardi, S., Cesare, S., et al. 2004, in New Frontiers in Stellar Interferometry, Proceedings of SPIE Volume 5491. Edited by Wesley A. Traub. Bellingham, WA: The International Society for Optical Engineering, 2004., p.528, ed. W. A. Traub, 528–+
- Groenewegen, M. A. T. 2008, A&A, 488, 25
- Keller, S. C. & Wood, P. R. 2006, ApJ, 642, 834
- Kervella, P., Nardetto, N., Bersier, D., Mourard, D., & Coudé du Foresto, V. 2004, A&A, 416, 941
- Lafrasse, S., Mella, G., Bonneau, D., et al. 2010, VizieR Online Data Catalog, 2300, 0
- Mathias, P., Gillet, D., Fokin, A. B., et al. 2006, A&A, 457, 575
- Petrov, R. G., Malbet, F., Weigelt, G., et al. 2007, A&A, 464, 1
- Pietrzyński, G., Thompson, I. B., Graczyk, D., et al. 2011, ApJ, 742, L20
- Prada Moroni, P. G., Gennaro, M., Bono, G., et al. 2012, ApJ, 749, 108
- Sasselov, D. D. & Lester, J. B. 1990, ApJ, 362, 333
- Storm, J., Gieren, W., Fouqué, P., et al. 2011, A&A, 534, A95
- Szabados, L. 1989, Commmunications of the Konkoly Observatory Hungary, 94, 1
- Szabados, L. 1990, MNRAS, 242, 285
- Szabados, L. 2003, Information Bulletin on Variable Stars, 5394, 1
- Tatulli, E., Millour, F., Chelli, A., et al. 2007, A&A, 464, 29