Through an innovative combination of multiple observing techniques and mod- eling, we are assembling a comprehensive understanding of the pulsation and close environment of Cepheids. We developed the SPIPS modeling tool that combines all observables (radial velocimetry, photometry, angular diameters from interferometry) to derive the relevant physical parameters of the star (effective temperature, infrared ex- cess, reddening,...) and the ratio of the distance and the projection factor d/p. We present the application of SPIPS to the long-period Cepheid RS Pup, for which we derive p = 1.25 +/- 0.06. The addition of this massive Cepheid consolidates the existing sample of p-factor measurements towards long-period pulsators. This allows us to conclude that p is constant or mildly variable around p = 1.29 +/- 0.04 (+/-3%) as a function of the pulsation period. The forthcoming Gaia DR2 will provide a considerable improvement in quantity and accuracy of the trigonometric parallaxes of Cepheids. From this sample, the SPIPS modeling tool will enable a robust calibration of the Cepheid distance scale.

The projection factor (p-factor) is an essential component of the classical Baade-Wesselink (BW) technique, that is commonly used to determine the distances to pulsating stars. It is a multiplicative parameter used to convert radial velocities into pulsational velocities. As the BW distances are linearly proportional to the p-factor, its accurate calibration for Cepheids is of critical importance for the reliability of their distance scale. We focus on the observational determination of the p-factor of the long-period Cepheid RS Pup (P = 41.5 days). This star is particularly important as this is one of the brightest Cepheids in the Galaxy and an analog of the Cepheids used to determine extragalactic distances. An accurate distance of 1910 +/- 80 pc (+/- 4.2%) has recently been determined for RS Pup using the light echoes propagating in its circumstellar nebula. We combine this distance with new VLTI/PIONIER interferometric angular diameters, photometry and radial velocities to derive the p-factor of RS Pup using the code Spectro-Photo-Interferometry of Pulsating Stars (SPIPS). We obtain p = 1.250 +/- 0.064 (+/-5.1%), defined for cross-correlation radial velocities. Together with measurements from the literature, the p-factor of RS Pup confirms the good agreement of a constant p = 1.293 +/- 0.039 (+/-3.0%) model with the observations. We conclude that the p-factor of Cepheids is constant or mildly variable over a broad range of periods (3.7 to 41.5 days).

The B-W method is used to determine the distance of Cepheids and consists in combining the angular size variations of the star, as derived from infrared surface-brightness relations or interferometry, with its linear size variation, as deduced from visible spectroscopy using the projection factor. While many Cepheids have been intensively observed by infrared beam combiners, only a few have been observed in the visible. This paper is part of a project to observe Cepheids in the visible with interferometry as a counterpart to infrared observations already in hand. Observations of delta Cep itself were secured with the VEGA/CHARA instrument over the full pulsation cycle of the star. These visible interferometric data are consistent in first approximation with a quasi-hydrostatic model of pulsation surrounded by a static circumstellar environment (CSE) with a size of theta_cse=8.9 +/- 3.0 mas and a relative flux contribution of f_cse=0.07+/-0.01. A model of visible nebula (a background source filling the field of view of the interferometer) with the same relative flux contribution is also consistent with our data at small spatial frequencies. However, in both cases, we find discrepancies in the squared visibilities at high spatial frequencies (maximum 2sigma) with two different regimes over the pulsation cycle of the star, phi=0.0-0.8 and phi=0.8-1.0. We provide several hypotheses to explain these discrepancies, but more observations and theoretical investigations are necessary before a firm conclusion can be drawn. For the first time we have been able to detect in the visible domain a resolved structure around delta~Cep. We have also shown that a simple model cannot explain the observations, and more work will be necessary in the future, both on observations and modelling.

The distance to pulsating stars is classically estimated using the parallax-of-pulsation (PoP) method, which combines spectroscopic radial velocity measurements and angular diameter estimates to derive the distance of the star. An important application of this method is the determination of Cepheid distances, in view of the calibration of their distance scale. However, the conversion of radial to pulsational velocities in the PoP method relies on a poorly calibrated parameter, the projection factor (p-factor). We aim to measure empirically the value of the p-factors of a homogeneous sample of nine Galactic Cepheids for which trigonometric parallaxes were measured with the Hubble Space Telescope Fine Guidance Sensor. We use the SPIPS algorithm, a robust implementation of the PoP method that combines photometry, interferometry, and radial velocity measurements in a global modeling of the pulsation. We obtained new interferometric angular diameters using the PIONIER instrument at the Very Large Telescope Interferometer, completed by data from the literature. Using the known distance as an input, we derive the value of the p-factor and study its dependence with the pulsation period. We find the following p-factors: 1.20 $\pm$ 0.12 for RT Aur, 1.48 $\pm$ 0.18 for T Vul, 1.14 $\pm$ 0.10 for FF Aql, 1.31 $\pm$ 0.19 for Y Sgr, 1.39 $\pm$ 0.09 for X Sgr, 1.35 $\pm$ 0.13 for W Sgr, 1.36 $\pm$ 0.08 for $\beta$ Dor, 1.41 $\pm$ 0.10 for $\zeta$ Gem, and 1.23 $\pm$ 0.12 for $\ell$ Car. These values are consistently close to p = 1.324 $\pm$ 0.024. We observe some dispersion around this average value, but the observed distribution is statistically consistent with a constant value of the p-factor as a function of the pulsation period. The error budget of our determination of the p-factor values is presently dominated by the uncertainty on the parallax, a limitation that will soon be waived by Gaia.

Baade-Wesselink-type (BW) techniques enable geometric distance measurements of Cepheid variable stars in the Galaxy and the Magellanic clouds. The leading uncertainties involved concern projection factors required to translate observed radial velocities (RVs) to pulsational velocities and recently discovered modulated variability. We carried out an unprecedented observational campaign involving long-baseline interferometry (VLTI/PIONIER) and spectroscopy (Euler/Coralie) to search for modulated variability in the long-period (P $\sim$ 35.5 d) Cepheid Carinae. We determine highly precise angular diameters from squared visibilities and investigate possible differences between two consecutive maximal diameters, $\Delta_{\rm{max}} \Theta$. We characterize the modulated variability along the line-of-sight using 360 high-precision RVs. Here we report tentative evidence for modulated angular variability and confirm cycle-to-cycle differences of $\ell$ Carinae's RV variability. Two successive maxima yield $\Delta_{\rm{max}} \Theta$ = 13.1 $\pm$ 0.7 (stat.) {\mu}as for uniform disk models and 22.5 $\pm$ 1.4 (stat.) {\mu}as (4% of the total angular variation) for limb-darkened models. By comparing new RVs with 2014 RVs we show modulation to vary in strength. Barring confirmation, our results suggest the optical continuum (traced by interferometry) to be differently affected by modulation than gas motions (traced by spectroscopy). This implies a previously unknown time-dependence of projection factors, which can vary by 5% between consecutive cycles of expansion and contraction. Additional interferometric data are required to confirm modulated angular diameter variations. By understanding the origin of modulated variability and monitoring its long-term behavior, we aim to improve the accuracy of BW distances and further the understanding of stellar pulsations.

The parallax of pulsation, and its implementations such as the Baade-Wesselink method and the infrared surface bright- ness technique, is an elegant method to determine distances of pulsating stars in a quasi-geometrical way. However, these classical implementations in general only use a subset of the available observational data. Freedman & Madore (2010) suggested a more physical approach in the implementation of the parallax of pulsation in order to treat all available data. We present a global and model-based parallax-of-pulsation method that enables including any type of observational data in a consistent model fit, the SpectroPhoto-Interferometric modeling of Pulsating Stars (SPIPS). We implemented a simple model consisting of a pulsating sphere with a varying effective temperature and a combina- tion of atmospheric model grids to globally fit radial velocities, spectroscopic data, and interferometric angular diameters. We also parametrized (and adjusted) the reddening and the contribution of the circumstellar envelopes in the near-infrared photometric and interferometric measurements. We show the successful application of the method to two stars: delta Cep and eta Aql. The agreement of all data fitted by a single model confirms the validity of the method. Derived parameters are compatible with publish values, but with a higher level of confidence. The SPIPS algorithm combines all the available observables (radial velocimetry, interferometry, and photometry) to estimate the physical parameters of the star (ratio distance/ p-factor, Teff, presence of infrared excess, color excess, etc). The statistical precision is improved (compared to other methods) thanks to the large number of data taken into account, the accuracy is improved by using consistent physical modeling and the reliability of the derived parameters is strengthened thanks to the redundancy in the data.

Long-baseline interferometry is an important technique to spatially resolve binary or multiple systems in close orbits. By combining several telescopes together and spectrally dispersing the light, it is possible to detect faint components around bright stars. Aims. We provide a rigorous and detailed method to search for high-contrast companions around stars, determine the detection level, and estimate the dynamic range from interferometric observations. We developed the code CANDID (Companion Analysis and Non-Detection in Interferometric Data), a set of Python tools that allows us to search systematically for point-source, high-contrast companions and estimate the detection limit. The search pro- cedure is made on a N x N grid of fit, whose minimum needed resolution is estimated a posteriori. It includes a tool to estimate the detection level of the companion in the number of sigmas. The code CANDID also incorporates a robust method to set a 3{\sigma} detection limit on the flux ratio, which is based on an analytical injection of a fake companion at each point in the grid. We used CANDID to search for the companions around the binary Cepheids V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul. First, we showed that our previous discoveries of the components orbiting V1334 Cyg and AX Cir were detected at > 13 sigmas. The companion around AW Per is detected at more than 15 sigmas with a flux ratio of f = 1.22 +/- 0.30 %. We made a possible detection of the companion orbiting RT Aur with f = 0.22 +/- 0.11 %. It was detected at 3.8{\sigma} using the closure phases only, and so more observations are needed to confirm the detection. We also set the detection limit for possible undetected companions. We found that there is no companion with a spectral type earlier than B7V, A5V, F0V, B9V, A0V, and B9V orbiting V1334 Cyg, AX Cir, RT Aur, AW Per, SU Cas, and T Vul, respectively.

The distances of pulsating stars, in particular Cepheids, are commonly measured using the parallax of pulsation technique. The different versions of this technique combine measurements of the linear diameter variation (from spectroscopy) and the angular diameter variation (from photometry or interferometry) amplitudes, to retrieve the distance in a quasi-geometrical way. However, the linear diameter amplitude is directly proportional to the projection factor (hereafter p-factor), which is used to convert spectroscopic radial velocities (i.e., disk integrated) into pulsating (i.e., photospheric) velocities. The value of the p-factor and its possible dependence on the pulsation period are still widely debated. Our goal is to measure an observational value of the p-factor of the type-II Cepheid kappa Pavonis, whose parallax was measured with an accuracy of 5% using HST/FGS. We used this parallax as a starting point to derive the p-factor of kappa Pav, using the SPIPS technique, which is a robust version of the parallax-of-pulsation method that employs radial velocity, interferometric and photometric data. We applied this technique to a combination of new VLTI/PIONIER optical interferometric angular diameters, new CORALIE and HARPS radial velocities, as well as multi-colour photometry and radial velocities from the literature. We obtain a value of p = 1.26 +/- 0.07 for the p-factor of kappa Pav. This result agrees with several of the recently derived Period-p-factor relationships from the literature, as well as previous observational determinations for Cepheids. Individual estimates of the p-factor are fundamental to calibrating the parallax of pulsation distances of Cepheids. Together with previous observational estimates, the projection factor we obtain points to a weak dependence of the p-factor on period.

As one of the most luminous Cepheids in the Milky Way, the 41.5-day RS Puppis is an analog of the long-period Cepheids used to measure extragalactic distances. An accurate distance to this star would therefore help anchor the zero-point of the bright end of the period-luminosity relation. But, at a distance of about 2 kpc, RS Pup is too far away for measuring a direct trigonometric parallax with a precision of a few percent with existing instrumentation. RS Pup is unique in being surrounded by a reflection nebula, whose brightness varies as pulses of light from the Cepheid propagate outwards. We present new polarimetric imaging of the nebula obtained with HST/ACS. The derived map of the degree of linear polarization pL allows us to reconstruct the three-dimensional structure of the dust distribution. To retrieve the scattering angle from the pL value, we consider two different polarization models, one based on a Milky Way dust mixture and one assuming Rayleigh scattering. Considering the derived dust distribution in the nebula, we adjust a model of the phase lag of the photometric variations over selected nebular features to retrieve the distance of RS Pup. We obtain a distance of 1910 +/- 80 pc (4.2%), corresponding to a parallax of 0.524 +/- 0.022 mas. The agreement between the two polarization models we considered is good, but the final uncertainty is dominated by systematics in the adopted model parameters. The distance we obtain is consistent with existing measurements from the literature, but light echoes provide a distance estimate that is not subject to the same systematic uncertainties as other estimators (e.g. the Baade-Wesselink technique). RS Pup therefore provides an important fiducial for the calibration of systematic uncertainties of the long-period Cepheid distance scale.

The stellar population in the central parsec of the Galaxy is dominated by an old (several Gyr) population, but young, massive stars dominate the luminosity function. We have studied the most luminous of these stars, GCIRS 7, in order to constrain the age of the recent star formation event in the Galactic Centre and to characterise it as an interferometric reference for observations of the Galactic Centre with the instrument GRAVITY, which will equip the Very Large Telescope Interferometer in the near future. We present the first H-band interferometric observations of GCIRS 7, obtained using the PIONIER visitor instrument on the VLTI using the four 8.2-m unit telescopes. In addition, we present unpublished K-band VLTI/AMBER data, build JHKL light-curves based on data spanning 4 decades, and measured the star's effective temperature using SINFONI spectroscopy. GCIRS 7 is marginally resolved at H-band (in 2013: uniform-disk diameter=1.076+/-0.093mas, R=960+/-92Rsun at 8.33+/-0.35kpc). We detect a significant circumstellar contribution at K-band. The star and its environment are variable in brightness and in size. The photospheric H-band variations are well modelled with two periods: P0~470+/-10 days (amplitude ~0.64mag) and long secondary period LSP~2700-2850 days (~1.1mag). As measured from CO equivalent width, <Teff>=3600+/-195K. The size, periods, luminosity (<Mbol>=-8.44+/-0.22) and effective temperature are consistent with an M1 supergiant with an initial mass of 22.5+/-2.5Msun and an age of 6.5-10Myr (depending on rotation). This age is in remarkable agreement with most estimates for the recent star formation event in the central parsec. Caution should be taken when using this star as an interferometric reference as it is variable in size, is surrounded by a variable circumstellar environment and large convection cells may form on its photosphere.

Aims: We aim at detecting and characterizing the main-sequence companion of the Cepheid AX Cir ($P_\mathrm{orb} \sim $ 18 yrs). The long-term objective is to estimate the mass of both components and the distance to the system. Methods: We used the PIONIER combiner at the VLT Interferometer to obtain the first interferometric measurements of the short-period Cepheid AX Cir and its orbiting component. Results: The companion is resolved by PIONIER at a projected separation $\rho = 29.2 \pm 0.2$ mas and projection angle $PA = 167.6 \pm 0.3^{\circ}$. We measured $H$-band flux ratios between the companion and the Cepheid of $0.90 \pm 0.10$ % and $0.75 \pm 0.17$ %, respectively at a pulsation phase for the Cepheid $\phi = 0.24$ and 0.48. The lower contrast at $\phi = 0.48$ is due to increased brightness of the Cepheid compared to the $\phi = 0.24$. This gives an average apparent magnitude $m\mathrm{_H (comp)} = 9.06 \pm 0.24$ mag. The limb-darkened angular diameter of the Cepheid at the two pulsation phases was measured to be $\theta_\mathrm{LD} = 0.839 \pm 0.023$ mas and $\theta_\mathrm{LD} = 0.742 \pm 0.020$ mas, respectively at $\phi = 0.24$ and 0.48. A lower limit on the total mass of the system was also derived based on our measured separation, we found $M_\mathrm{T} \geq 9.7 \pm 0.6 M_\odot$.

Aims. We study the close environment of nearby Cepheids using high spatial resolution observations in the mid-infrared with the VLTI/MIDI instrument, a two-beam interferometric recombiner. Methods. We obtained spectra and visibilities for the classical Cepheids X Sgr and T Mon. We fitted the MIDI measurements, supplemented by B, V, J, H, K literature photometry, with the numerical transfer code DUSTY to determine the dust shell parameters. We used a typical dust composition for circumstellar environments. Results. We detect an extended dusty environment in the spectra and visibilities for both stars, although T Mon might suffer from thermal background contamination. We attribute this to the presence of a circumstellar envelope (CSE) surrounding the Cepheids. This is optically thin for X Sgr (tau(0.55microns) = 0.008), while it appears to be thicker for T Mon (tau(0.55micron) = 0.15). They are located at about 15-20 stellar radii. Following our previous work, we derived a likely period-excess relation in the VISIR PAH1 filter, f(8.6micron)[%]= 0.81(+/-0.04)P[day]. We argue that the impact of CSEs on the mid-IR period-luminosity (P-L) relation cannot be negligible because they can bias the Cepheid brightness by up to about 30 %. For the K-band P-L relation, the CSE contribution seems to be lower (< 5 %), but the sample needs to be enlarged to firmly conclude that the impact of the CSEs is negligible in this band.