Classical Cepheids (CCs) are at the heart of the empirical extragalactic distance ladder. Milky Way CCs are the only stars of this class accessible to trigonometric parallax measurements. Until recently, the most accurate trigonometric parallaxes of Milky Way CCs were the HST/FGS measurements collected by Benedict et al. (2002, 2007) and HST/WFC3 measurements by Riess et al. (2018). Unfortunately, the second Gaia data release (GDR2) has not yet delivered reliable parallaxes for Galactic CCs, failing to replace the HST as the foundation of the Galactic calibrations of the Leavitt law. We aim at calibrating independently the Leavitt law of Milky Way CCs based on the GDR2 catalog of trigonometric parallaxes. As a proxy for the parallaxes of a sample of 23 Galactic CCs, we adopt the GDR2 parallaxes of their spatially resolved companions. As the latter are unsaturated, photometrically stable stars, this novel approach allows us to bypass the GDR2 bias on the parallax of the CCs that is induced by saturation and variability. We present new Galactic calibrations of the Leavitt law in the J, H, K, V, Wesenheit WH and Wesenheit WVK bands based on the GDR2 parallaxes of the CC companions. We show that the adopted value of the zero point of the GDR2 parallaxes, within a reasonable range, has a limited impact on our Leavitt law calibration.

Classical Cepheids provide the foundation for the empirical extragalactic distance ladder. Milky Way Cepheids are the only stars of this class accessible to trigonometric parallax measurements. However, the parallaxes of Cepheids from the second Gaia data release (GDR2) are affected by systematics due to the absence of chromaticity correction and occasionally by saturation. As a proxy for the parallaxes of 36 Galactic Cepheids, we adopt either the GDR2 parallaxes of their spatially resolved companions or the GDR2 parallax of their host open cluster. This novel approach allows us to bypass the systematics on the GDR2 Cepheids parallaxes that is induced by saturation and variability. We adopt a GDR2 parallax offset of 0.046 mas with an uncertainty of 0.015 mas that covers most of the recent estimates. We present new Galactic calibrations of the Leavitt law in the V, J, H, K_S and Wesenheit W_H bands. We compare our results with previous measurements from the Hubble Space Telescope (HST) and compute a revised value for the Hubble constant anchored to Milky Way Cepheids. From an initial Hubble constant of 76.18 +/- 2.37 km/s/Mpc based on parallax measurements without Gaia, we derive a revised value by adopting companion and average cluster parallaxes in place of direct Cepheid parallaxes and we find H0 = 73.07 +/- 1.75 (stat. + syst.) +/- 1.88 (ZP) km/s/Mpc when all Cepheids are considered, and H0 = 73.51 +/- 1.76 (stat. + syst.) +/- 1.91 (ZP) km/s/Mpc for fundamental mode pulsators only.

Despite observational evidences, InfraRed (IR) excess of classical Cepheids are seldom studied and poorly understood, but probably induces systematics on the Period-Luminosity (PL) relation used in the calibration of the extragalactic distance scale. This study aims to understand the physical origin of the IR excess found in the spectral energy distribution (SED) of 5 Cepheids : RS Pup (P=41.46d), zeta Gem (P=10.15d), eta Aql (P=7.18d), V Cen (P=5.49d) and SU Cyg (P=3.85d). A time series of atmospheric models along the pulsation cycle are fitted to a compilation of data, including optical and near-IR photometry, Spitzer spectra (secured at a specific phase), interferometric angular diameters, effective temperature and radial velocity measurements. Herschel images in two bands are also analyzed qualitatively. In this fitting process, based on the SPIPS algorithm, a residual is found in the SED, whatever the pulsation phase, and for wavelengths larger than about $1.2\mu$m, which corresponds to the so-determined infrared excess of Cepheids. This IR excess is then corrected from interstellar medium absorption in order to infer or not the presence of dust shells, and is finally used in order to fit a model of a shell of ionized gas. For all Cepheids, we find a continuum IR excess increasing up to about -0.1 magnitudes at 30$\mu$m, which cannot be explained by a hot or cold dust model of CircumStellar Environment (CSE). We show, for the first time, that the IR excess of Cepheids can be explained by free-free emission from a thin shell of ionized gas, with a thickness of about 15% of the star radius, a mass of $10^{-9}-10^{-7}$ Msol and a temperature ranging from 3500 to 4500K. This result has to be tested with interferometers operating in visible, in the mid-IR or in the radio domain. The impact of such CSEs of ionized gas on the PL relation of Cepheids needs also more investigations.

Accurate radial velocities ($v_{\rm rad}$) of Cepheids are mandatory within the context of distance measurements via the Baade-Wesselink technique. The most common $v_{\rm rad}$ derivation method consists in cross-correlating the observed spectrum with a binary template and measuring a velocity on the resulting profile. Yet for Cepheids, the spectral lines selected within the template as well as the way of fitting the cross-correlation function (CCF) have a significant impact on the measured $v_{\rm rad}$. We detail the steps to compute consistent Cepheid CCFs and $v_{\rm rad}$, and we characterise the impact of Cepheid spectral properties and $v_{\rm rad}$ computation method on the resulting line profiles. We collected more than 3900 high-resolution spectra from seven different spectrographs of 64 classical Cepheids. These spectra were standardised through a single process on pre-defined wavelength ranges. We built six correlation templates selecting un-blended lines of different depths from a synthetic Cepheid spectrum, on three different wavelength ranges from 390 to 800 nm. Each spectrum was cross-correlated with these templates to build the corresponding CCFs. We derived a set of line profile observables as well as three different $v_{\rm rad}$ measurements from each CCF. This study confirms that both the template wavelength range, its mean line depth and width, and the $v_{\rm rad}$ computation method significantly impact the $v_{\rm rad}$. Deriving more robust Cepheid $v_{\rm rad}$ time series require to minimise the asymmetry of the line profile and its impact on the $v_{\rm rad}$. Centroid $v_{\rm rad}$, that exhibit slightly smaller amplitudes but significantly smaller scatter than Gaussian or biGaussian $v_{\rm rad}$, should thus be favoured. Stronger lines are also less asymmetric and lead to more robust $v_{\rm rad}$ than weaker lines.

The dynamical structure of the atmosphere of Cepheids has been well studied in the optical. Several authors have found very interesting spectral features in the J band, but little data have been secured beyond 1.6 um. However, such observations can probe different radial velocities and line asymmetry regimes, and are able to provide crucial insights into stellar physics. Our goal was to investigate the infrared line-forming region in the K-band domain, and its impact on the projection factor and the k-term of Cepheids. We secured CRIRES observations for the long-period Cepheid l Car, with a focus on the unblended spectral line NaI2208.969 nm. We measured the corresponding radial velocities (by using the first moment method) and the line asymmetries (by using the bi-Gaussian method). These quantities are compared to the HARPS visible spectra we previously obtained on l Car. The optical and infrared radial velocity curves show the same amplitude (only about 3% of difference), with a slight radial velocity shift of about 0.5 +/- 0.3 km s^-1 between the two curves. Around the minimum radius (phase ~ 0.9) the visible radial velocity curve is found in advance compared to the infrared one (phase lag), which is consistent with an infrared line forming higher in the atmosphere (compared to the visible line) and with a compression wave moving from the bottom to the top of the atmosphere during maximum outward velocity. The asymmetry of the K-band line is also found to be significantly different from that of the optical line.

We present the modeling tool we developed to incorporate multi-technique observations of Cepheids in a single pulsation model: the Spectro-Photo-Interferometry of Pulsating Stars (SPIPS). The combination of angular diameters from optical interferometry, radial velocities and photometry with the coming Gaia DR2 parallaxes of nearby Galactic Cepheids will soon enable us to calibrate the projection factor of the classical Parallax-of-Pulsation method. This will extend its applicability to Cepheids too distant for accurate Gaia parallax measurements, and allow us to precisely calibrate the Leavitt law's zero point. As an example application, we present the SPIPS model of the long-period Cepheid RS Pup that provides a measurement of its projection factor, using the independent distance estimated from its light echoes.

We focus on empirically measure the p-factor of a homogeneous sample of 29 LMC and 10 SMC Cepheids for which an accurate average LMC/SMC distance were estimated from eclipsing binary systems. We used the SPIPS algorithm, which is an implementation of the BW method. As opposed to other conventional use, SPIPS combines all observables, i.e. radial velocities, multi-band photometry and interferometry into a consistent physical modeling to estimate the parameters of the stars. The large number and their redundancy insure its robustness and improves the statistical precision. We successfully estimated the p-factor of several MC Cepheids. Combined with our previous Galactic results, we find the following P-p relation: -0.08(log P-1.18)+1.24. We find no evidence of a metallicity dependent p-factor. We also derive a new calibration of the P-R relation, logR=0.684(log P-0.517)+1.489, with an intrinsic dispersion of 0.020. We detect an IR excess for all stars at 3.6 and 4.5um, which might be the signature of circumstellar dust. We measure a mean offset of $\Delta m_{3.6}=0.057$mag and $\Delta m_{4.5}=0.065$mag. We provide a new P-p relation based on a multi-wavelengths fit, and can be used for the distance scale calibration from the BW method. The dispersion is due to the MCs width we took into account because individual Cepheids distances are unknown. The new P-R relation has a small intrinsic dispersion, i.e. 4.5% in radius. Such precision will allow us to accurately apply the BW method to nearby galaxies. Finally, the IR excesses we detect raise again the issue on using mid-IR wavelengths to derive P-L relation and calibrate the $H_0$. These IR excesses might be the signature of circumstellar dust, and are never taken into account when applying the BW method at those wavelengths. Our measured offsets may give an average bias of 2.8% on the distances derived through mid-IR P-L relations.

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).