Search

Cepheid masses continue to be important tests of evolutionary tracks for intermediate mass stars as well as important predictors of their future fate. For systems where the secondary is a B star, {\it Hubble Space Telescope} ultraviolet spectra have been obtained. From these spectra a temperature can be derived, and from this a mass of the companion M$_2$. Once {\it Gaia} DR4 is available, proper motions can be used to determine the inclination of the orbit. Combining mass of the companion, M$_2$, the mass function from the ground-based orbit of the Cepheid and the inclination produces the mass of the Cepheid, M$_1$. The Cepheid system FN Vel is used here to demonstrate this approach and what limits can be put on the Cepheid mass for inclination between 50 and 130$^o$., Comment: accepted by AJ

The 30 year orbit of the Cepheid Polaris has been followed with observations by the CHARA Array (Center for High Angular Resolution Astronomy) from 2016 through 2021. An additional measurement has been made with speckle interferometry at the Apache Point Observatory. Detection of the companion is complicated by its comparative faintness--an extreme flux ratio. Angular diameter measurements appear to show some variation with pulsation phase. Astrometric positions of the companion were measured with a custom grid-based model-fitting procedure and confirmed with the CANDID software. These positions were combined with the extensive radial velocities discussed by Torres (2023) to fit an orbit. Because of the imbalance of the sizes of the astrometry and radial velocity datasets, several methods of weighting are discussed. The resulting mass of the Cepheid is 5.13$\pm$ 0.28 $M_\odot$. Because of the comparatively large eccentricity of the orbit (0.63), the mass derived is sensitive to the value found for the eccentricity. The mass combined with the distance shows that the Cepheid is more luminous than predicted for this mass from evolutionary tracks. The identification of surface spots is discussed. This would give credence to the identification of photometric variation with a period of approximately 120 days as a rotation period. Polaris has some unusual properties (rapid period change, a phase jump, variable amplitude, unusual polarization). However, a pulsation scenario involving pulsation mode, orbital periastron passage (Torres 2023), and low pulsation amplitude can explain these characteristics within the framework of pulsation seen in Cepheids., Comment: ApJ in press

The Cepheid AW Per is a component in a multiple system with a long period orbit. The radial velocities of Griffin (2016) cover the 38 year orbit well. An extensive program of interferometry with the CHARA array is reported here, from which the long period orbit is determined. In addition, a {\it Hubble Space Telescope} high resolution spectrum in the ultraviolet demonstrates that the companion is itself a binary with nearly equal mass components. These data combined with a distance from {\it Gaia} provide a mass of the Cepheid (primary) of M$_1$ = 6.79 $\pm$ 0.85 $M_\odot$. The combined mass of the secondary is M$_S$ = 8.79 $\pm$ 0.50 $M_\odot$. The accuracy of the mass will be improved after the fourth Gaia data release expected in approximately two years., Comment: Accepted for ApJ

An inner companion has recently been discovered orbiting the prototype of classical Cepheids, delta Cep, whose orbital parameters are still not fully constrained. We collected new precise radial velocity measurements of delta Cep in 2019 using the HARPS-N spectrograph mounted at the Telescopio Nazionale Galileo. Using these radial velocity measurements, we aimed to improve the orbital parameters of the system. We considered a template available in the literature as a reference for the radial velocity curve of the pulsation of the star. We then calculated the residuals between our global dataset (composed of the new 2019 observations plus data from the literature) and the template as a function of the pulsation phase and the barycentric Julian date. This provides the orbital velocity of the Cepheid component. Using a Bayesian tool, we derived the orbital parameters of the system. Considering priors based on already published Gaia constraints, we find for the orbital period a maximum a posteriori probability of Porb=9.32+/-0.03 years (uncertainties correspond to the 95% highest density probability interval), and we obtain an eccentricity e=0.71+/-0.02, a semimajor axis a=0.029 +/-0.003 arcsecond, and a center-of-mass velocity V0=-17.28+/-0.08 km/s, among other parameters. In this short analysis we derive the orbital parameters of the delta Cep inner binary system and provide a cleaned radial velocity curve of the pulsation of the star, which will be used to study its Baade-Wesselink projection factor in a future publication., Comment: accepted for publication in A&A Letters

Binary Cepheids with giant companions are crucial for studying the physical properties of Cepheid variables, providing the best means to measure their masses. Systems composed of two Cepheids are even more important but to date, only one such system in the Large Magellanic Cloud (LMC) was known. Our current aim is to increase the number of these systems tenfold and provide their basic characteristics. The final goal is to obtain the physical properties of the component Cepheids, including their masses and radii, and to learn about their evolution in the multiple systems, also revealing their origin. We started a spectroscopic monitoring of nine unresolved pairs of Cepheids from the OGLE catalog, to check if they are gravitationally bound. Two of these so-called double Cepheids are located in the LMC, five in the Small Magellanic Cloud (SMC), and two in the Milky Way (MW). We report the spectroscopic detection of binarity of all 9 of these double Cepheids with orbital periods from 2 to 18 years. This increases the number of known binary double (BIND) Cepheids from 1 to 10 and triples the number of all confirmed double-lined binary (SB2) Cepheids. For five BIND Cepheids disentangled pulsational light curves of the components show anti-correlated phase shifts due to orbital motion. We show the first empirical evidence that typical period-luminosity relations (PLRs) are rather binary Cepheid PLRs that include the companion's light. The statistics of pulsation period ratios of BIND Cepheids do not agree with those expected for pairs of the same-age Cepheids. These ratios together with the mass ratios far from unity suggest merger-origin of at least one component for about half of the systems. The SMC and MW objects are the first found in SB2 systems composed of giants in their host galaxies. The Milky Way BIND Cepheids are also the closest such systems, being located at about 11 and 26 kpc., Comment: 13 pages, 7 figures, 4 tables, accepted for publication in Astronomy & Astrophysics

The book consists of a number of short articles that present achievements of the Araucaria members, collaborators, and friends, in various aspects of distance determinations and related topics. It celebrates the 20-year anniversary of the Araucaria Project, acknowledges the people who worked for its success, and popularises our methods and results among broader readership. This book is a part of a project that has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 695099., Comment: 114 pages, book published in 2021 on behalf of the Nicolaus Copernicus Astronomical Center of the Polish Academy of Sciences, to celebrate 20 years of the Arauria Project

We aim to measure very precise and accurate model-independent masses and distances of detached binary stars. Precise masses at the $< 1$% level are necessary to test and calibrate stellar interior and evolution models, while precise and independent orbital parallaxes are essential to check for the next Gaia data releases. We combined RV measurements with interferometric observations to determine orbital and physical parameters of ten double-lined spectroscopic systems. We report new relative astrometry from VLTI/GRAVITY and, for some systems, new VLT/UVES spectra to determine the radial velocities of each component. We measured the distance of ten binary systems and the mass of their components with a precision as high as 0.03% (average level 0.2%). They are combined with other stellar parameters (effective temperatures, radii, flux ratios, etc.) to fit stellar isochrones and determine their evolution stage and age. We also compared our orbital parallaxes with Gaia and showed that half of the stars are beyond $1\sigma$ with our orbital parallaxes; although, their RUWE is below the frequently used cutoff of 1.4 for reliable Gaia astrometry. By fitting the telluric features in the GRAVITY spectra, we also estimated the accuracy of the wavelength calibration to be $\sim 0.02$% in high and medium spectral resolution modes. We demonstrate that combining spectroscopic and interferometric observations of binary stars provides extremely precise and accurate dynamical masses and orbital parallaxes. As they are detached binaries, they can be used as benchmark stars to calibrate stellar evolution models and test the Gaia parallaxes., Comment: Accepted for publication in A&A

The Baade-Wesselink (BW) method of distance determination of Cepheids is used to calibrate the distance scale. Various versions of this method are mainly based on interferometry and/or the surface-brightness color relation (SBCR). We quantify the impact of the SBCR, its slope, and its zeropoint on the projection factor. This quantity is used to convert the pulsation velocity into the radial velocity in the BW method. We also study the impact of extinction and of a potential circumstellar environment on the projection factor. We analyzed HARPS-N spectra of eta Aql to derive its radial velocity curve using different methods. We then applied the inverse BW method using various SBCRs in the literature in order to derive the BW projection factor. We find that the choice of the SBCR is critical: a scatter of about 8% is found in the projection factor for different SBCRs in the literature. The uncertainty on the coefficients of the SBCR affects the statistical precision of the projection factor only little (1-2\%). Confirming previous studies, we find that the method with which the radial velocity curve is derived is also critical, with a potential difference on the projection factor of 9%. An increase of 0.1 in E(B-V) translates into a decrease in the projection factor of 3%. A 0.1 magnitude effect of a circumstellar envelope (CSE) in the visible domain is rather small on the projection factor, about 1.5%. However, we find that a 0.1 mag infrared excess in the K band due to a CSE can increase the projection factor by about 6%. The impact of the surface-brightness color relation on the BW projection factor is found to be critical. Efforts should be devoted in the future to improve the SBCR of Cepheids empirically, but also theoretically, taking their CSE into account as well., Comment: Accepted for publication in Astronomy & Astrophysics

Planet(s) in binaries are unique architectures for testing predictions of planetary formation and evolution theories in very hostile environments. We used the IRDIS dual-band imager of SPHERE at VLT, and the speckle interferometric camera HRCAM of SOAR, to acquire high-angular resolution images of HD 196885 AB between 2015 and 2020. Radial velocity observations have been extended over almost 40 yr extending the radial velocity measurements HD 196885 A and resolving both the binary companion and the inner giant planet HD 196885 Ab. Finally, we took advantage of the exquisite astrometric precision of the dual-field mode of VLTI/GRAVITY (down to 30 {\mu}as) to monitor the relative position of HD 196885 A and B to search for the 3.6 yr astrometric wobble of the circumprimary planet Ab imprinted on the binary separation. Our observations enable to accurately constrain the orbital properties of the binary HD 196885 AB, seen on an inclined and retrograde orbit (iAB = 120.43 deg) with a semi-major axis of 19.78 au, and an eccentricity of 0.417. The GRAVITY measurements confirm for the first time the nature of the inner planet HD 196885 Ab by rejecting all families of pole-on solutions in the stellar or brown dwarf masses. The most favored island of solutions is associated with a Jupiter-like planet (MAb = 3.39 MJup), with moderate eccentricity (eAaAb = 0.44), and inclination close to 143.04 deg. This results points toward a significant mutual inclination (Phi = 24.36 deg) between the orbital planes (relative to the star) of the binary companion B and the planet Ab. Our dynamical simulations indicate that the system is dynamically stable over time. Eccentricity and mutual inclination variations could be expected for moderate von Zipele Kozai Lidov cycles that may affect the inner planet., Comment: 12 pages, 6 figures, accepted in A&A

Kernel phase imaging (KPI) enables the direct detection of substellar companions and circumstellar dust close to and below the classical (Rayleigh) diffraction limit. We present a kernel phase analysis of JWST NIRISS full pupil images taken during the instrument commissioning and compare the performance to closely related NIRISS aperture masking interferometry (AMI) observations. For this purpose, we develop and make publicly available the custom "Kpi3Pipeline" enabling the extraction of kernel phase observables from JWST images. The extracted observables are saved into a new and versatile kernel phase FITS file (KPFITS) data exchange format. Furthermore, we present our new and publicly available "fouriever" toolkit which can be used to search for companions and derive detection limits from KPI, AMI, and long-baseline interferometry observations while accounting for correlated uncertainties in the model fitting process. Among the four KPI targets that were observed during NIRISS instrument commissioning, we discover a low-contrast (~1:5) close-in (~1 $\lambda/D$) companion candidate around CPD-66~562 and a new high-contrast (~1:170) detection separated by ~1.5 $\lambda/D$ from 2MASS~J062802.01-663738.0. The 5-$\sigma$ companion detection limits around the other two targets reach ~6.5 mag at ~200 mas and ~7 mag at ~400 mas. Comparing these limits to those obtained from the NIRISS AMI commissioning observations, we find that KPI and AMI perform similar in the same amount of observing time. Due to its 5.6 times higher throughput if compared to AMI, KPI is beneficial for observing faint targets and superior to AMI at separations >325 mas. At very small separations (<100 mas) and between ~250-325 mas, AMI slightly outperforms KPI which suffers from increased photon noise from the core and the first Airy ring of the point-spread function., Comment: 34 pages, 17 figures, accepted for publication in PASP

Achernar, the closest and brightest classical Be star, presents rotational flattening, gravity darkening, occasional emission lines due to a gaseous disk, and an extended polar wind. It is also a member of a close binary system with an early A-type dwarf companion. We aim to determine the orbital parameters of the Achernar system and to estimate the physical properties of the components. We monitored the relative position of Achernar B using a broad range of high angular resolution instruments of the VLT/VLTI (VISIR, NACO, SPHERE, AMBER, PIONIER, GRAVITY, and MATISSE) over a period of 13 years (2006-2019). These astrometric observations are complemented with a series of more than 700 optical spectra for the period from 2003 to 2016. We determine that Achernar B orbits the Be star on a seven-year period, eccentric orbit (e = 0.7255 +/- 0.0014) which brings the two stars within 2 au at periastron. The mass of the Be star is found to be mA = 6.0 +/- 0.6 Msun for a secondary mass of mB = 2.0 +/- 0.1 Msun. We find a good agreement of the parameters of Achernar A with the evolutionary model of a critically rotating star of 6.4 Msun at an age of 63 million years. We also identify a resolved comoving low-mass star, which leads us to propose that Achernar is a member of the Tucana-Horologium moving group. Achernar A is presently in a short-lived phase of its evolution following the turn-off, during which its geometrical flattening ratio is the most extreme. Considering the orbital parameters, no significant interaction occurred between the two components, demonstrating that Be stars may form through a direct, single-star evolution path without mass transfer. Since component A will enter the instability strip in a few hundred thousand years, Achernar appears to be a promising progenitor of the Cepheid binary systems., Comment: 27 pages, Astronomy & Astrophysics, in press. The full set of continuum normalized high resolution spectra of Achernar is available at https://doi.org/10.5281/zenodo.6977303

Surface brightness -- colour relations (SBCRs) are very useful tools for predicting the angular diameters of stars. They offer the possibility to calculate very precise spectrophotometric distances by the eclipsing binary method or the Baade-Wesselink method. Double-lined Detached Eclipsing Binary stars (SB2 DEBs) with precisely known trigonometric parallaxes allow for a calibration of SBCRs with unprecedented precision. In order to improve such calibrations, it is important to enlarge the calibration sample of suitable eclipsing binaries with very precisely determined physical parameters. We carefully chose a sample of ten SB2 DEBs in the solar neighbourhood which contain inactive main-sequence components. The components have spectral types from early A to early K. All systems have high-precision parallaxes from the Gaia mission. We analysed high precision ground- and space-based photometry simultaneously with the radial velocity curves derived from HARPS spectra. We used spectral disentangling to obtain the individual spectra of the components and used these to derive precise atmospheric parameters and chemical abundances. For almost all components, we derived precise surface temperatures and metallicities. We derived absolute dimensions for 20 stars with an average precision of 0.2% and 0.5% for masses and radii, respectively. Three systems show slow apsidal motion. One system, HD 32129, is most likely a triple system with a much fainter K6V companion. Also three systems contain metallic-line components and show strong enhancements of barium and ittrium. The components of all systems compare well to the SBCR derived before from the detached eclipsing binary stars. With a possible exception of HD 32129, they can be used to calibrate SBCRs with a precision better than 1% with available Gaia DR3 parallaxes., Comment: Accepted in Astronomy & Astrophysics, 18 pages

AU Mic is a young and nearby M-dwarf star harbouring a circumstellar debris disk and one recently discovered planet on an 8d orbit. Large-scale structures within the disk were also discovered and are moving outward at high velocity. We aim at studying this system with the highest spatial resolution in order to probe the innermost regions and to search for additional low-mass companion or set detection limits. The star was observed with two different techniques probing complementary spatial scales. We obtained new SAM observations with SPHERE, which we combined with data from NACO, PIONIER and GRAVITY. We did not detect additional companions within 0.02-7au from the star. We determined magnitude upper limits for companions of H~9.8mag within 0.02-0.5au, Ks~11.2mag within 0.4-2.4au and L'~10.7mag within 0.7-7au. Using theoretical isochrones, we converted into mass upper limits of ~17Mjup, ~12Mjup and ~9jup, respectively. The PIONIER observations allowed us to determine the angular diameter of AU Mic, 0.825+/-0.050mas, which converts to R = 0.862+/-0.052Rsun. We did not detect the newly discovered planets, but we derived upper limit masses for the innermost region of AU Mic. We do not have any detection with a significance beyond 3sigma, the most significant signal with PIONIER being 2.9sigma and with SPHERE being 1.6\sigma. We applied the pyMESS2 code to estimate the detection probability of companions by combining radial velocities, SPHERE imaging and our interferometric detection maps. We show that 99% of the companions down to ~0.5Mjup can be detected within 0.02au or 1Mjup down to 0.2au. The low-mass planets orbiting at <0.11au will not be directly detectable with the current AO and interferometric instruments due to its close orbit and very high contrast (~10e-10 in K). It will be also below the angular resolution and contrast limit of the next ELT IR imaging instruments., Comment: Accepted for publication in A&A. 15 pages. New version includes two new co-authors

We present time-series photometry of 21 nearby Type II Cepheids in the near-infrared J, H and Ks passbands. We use this photometry, together with the Third Gaia Early Data Release parallaxes, to determine for the first time period-luminosity relations (PLRs) for Type II Cepheids from field representatives of these old pulsating stars in the near-infrared regime. We found PLRs to be very narrow for BL Herculis stars, which makes them candidates for precision distance indicators. We then use archival photometry and the most accurate distance obtained from eclipsing binaries to recalibrate PLRs for Type II Cepheids in the Large Magellanic Cloud (LMC). Slopes of our PLRs in the Milky Way and in the LMC differ by slightly more than 2{\sigma} and are in a good agreement with previous studies of the LMC, Galactic Bulge and Galactic Globular Clusters Type II Cepheids samples. We use PLRs of Milky Way Type II Cepheids to measure the distance to the LMC and we obtain a distance modulus of 18.540$\pm$0.026(stat.)$\pm$0.034(syst.)mag in the WJKs Wesenheit index. We also investigate the metallicity effect within our Milky Way sample and we find rather significant value of about -0.2mag/dex in each band meaning that more metal-rich Type II Cepheids are intrinsically brighter than their more metal-poor counterparts, in agreement with the value obtained from Type II Cepheids in Galactic Globular Clusters. The main source of systematic error on our Milky Way PLRs calibration and the LMC distance is the current uncertainty of the Gaia parallax zero point., Comment: Accepted for ApJ

Classical Be stars are possible products of close binary evolution, in which the mass donor becomes a hot, stripped O or B-type subdwarf (sdO/sdB), and the mass gainer spins up and grows a disk to become a Be star. While several Be+sdO binaries have been identified, dynamical masses and other fundamental parameters are available only for a single Be+sdO system, limiting the confrontation with binary evolution models. In this work, we present direct interferometric detections of the sdO companions of three Be stars 28 Cyg, V2119 Cyg, and 60 Cyg, all of which were previously found in UV spectra. For two of the three Be+sdO systems, we present first orbits and preliminary dynamical masses of the components, revealing that one of them could be the first identified progenitor of a Be/X-ray binary with a neutron star companion. These results provide new sets of fundamental parameters that are crucially needed to establish the evolutionary status and origin of Be stars., Comment: accepted to ApJ

As primary anchors of the distance scale, Cepheid stars play a crucial role in our understanding of the distance scale of the Universe because of their period-luminosity relation. Determining precise and consistent parameters (radius, temperature, color excess, and projection factor) of Cepheid pulsating stars is therefore very important. With the high-precision parallaxes delivered by the early third Gaia data release, we aim to derive various parameters of Cepheid stars in order to calibrate the period-luminosity and period-radius relations and to investigate the relation of period to p-factor. We applied an implementation of the parallax-of-pulsation method through the algorithm called Spectro-Photo-Interferometry of Pulsating Stars, which combines all types of available data for a variable star in a global modeling of its pulsation. We present the SPIPS modeling of a sample of 63 Galactic Cepheids. Adopting Gaia EDR3 parallaxes as an input associated with the best available dataset, we derive consistent values of parameters for these stars such as the radius, multiband apparent magnitudes, effective temperatures, color excesses, period changes, Fourier parameters, and the projection factor. We then derive new calibrations of the period-luminosity and period-radius relations. After investigating the dependences of the p-factor on the parameters of the stars, we find a high dispersion of its values and no evidence of its correlation with the period or with any other parameters. Statistically, the p-factor has an average value of p=1.26$\pm$0.07, but with an unsatisfactory agreement. In absence of any clear correlation between the p-factor and other quantities, the best agreement is obtained under the assumption that the p-factor can take any value in a band with a width of 0.15. This result highlights the need for a further examination of the physics behind the p-factor., Comment: 29 pages, 23 figures, 5 tables. Accepted for publication in Astronomy & Astrophysiscs journal (30/08/2021)

Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have developed MATISSE, the Multi AperTure mid-Infrared SpectroScopic Experiment, to access high resolution imaging in a wide spectral domain and explore topics such: stellar activity and mass loss; planet formation and evolution in the gas and dust disks around young stars; accretion processes around super massive black holes in AGN. Methods:The instrument is a spectro-interferometric imager covering three atmospheric bands (L,M,N) from 2.8 to 13.0 mu, combining four optical beams from the VLTI's telscopes. Its concept, related observing procedure, data reduction and calibration approach are the product of 30 years of instrumental research. The instrument utilizes a multi-axial beam combination that delivers spectrally dispersed fringes. The signal provides the following quantities at several spectral resolutions: photometric flux, coherent fluxes, visibilities, closure phases, wavelength differential visibilities and phases, and aperture-synthesis imaging. Results:We provide an overview of the physical principle of the instrument and its functionalities, the characteristics of the delivered signal, a description of the observing modes and of their performance limits. An ensemble of data and reconstructed images are illustrating the first acquired key observations. Conclusion:The instrument has been in operation at Cerro Paranal, ESO, Chile since 2018, and has been open for science use by the international community since April 2019. The first scientific results are being published now., Comment: 24 pages, 26 figures, submitted to Astronomy & Astrophysics

HD 98800 is a young ($\sim10$ Myr old) and nearby ($\sim45$ pc) quadruple system, composed of two spectroscopic binaries orbiting around each other (AaAb and BaBb), with a gas-rich disk in polar configuration around BaBb. While the orbital parameters of BaBb and AB are relatively well constrained, this is not the case for AaAb. A full characterisation of this quadruple system can provide insights on the formation of such a complex system. The goal of this work is to determine the orbit of the AaAb subsystem and refine the orbital solution of BaBb using multi-epoch interferometric observations with the VLTI/PIONIER and radial velocities. The PIONIER observations provide relative astrometric positions and flux ratios for both AaAa and BaBb subsystems. Combining the astrometric points with radial velocity measurements, we determine the orbital parameters of both subsystems. We refined the orbital solution of BaBb and derived, for the first time, the full orbital solution of AaAb. We confirmed the polar configuration of the circumbinary disk around BaBb. From our solutions, we also inferred the dynamical masses of AaAb ($M_{Aa} = 0.93 \pm 0.09$ and $M_{Ab} = 0.29 \pm 0.02$ M$_{\odot}$). We also revisited the parameters of the AB outer orbit. Using the N-body simulation, we show that the system should be dynamically stable over thousands of orbital periods and that it made preliminary predictions for the transit of the disk in front of AaAb which is estimated to start around 2026. We discuss the lack of a disk around AaAb, which can be explained by the larger X-ray luminosity of AaAb, promoting faster photo-evaporation of the disk. High-resolution infrared spectroscopic observations would provide radial velocities of Aa and Ab (blended lines in contemporary observations), which would allow us to calculate the dynamical masses of Aa and Ab independently of the parallax of BaBb., Comment: Accepted for publication in Astronomy \& Astrophysiscs journal (01/09/2021)

X-ray bursts have recently been discovered in the Cepheids $\delta$ Cep and $\beta$ Dor modulated by the pulsation cycle. We have obtained an observation of the Cepheid $\eta$ Aql with the XMM-Newton satellite at the phase of maximum radius, the phase at which there is a burst of X-rays in $\delta$ Cep. No X-rays were seen from the Cepheid $\eta$ Aql at this phase, and the implications for Cepheid upper atmospheres are discussed. We have also used the combination of X-ray sources and Gaia and 2MASS data to search for a possible grouping around the young intermediate mass Cepheid. No indication of such a group was found., Comment: accepted by AJ

We aim to investigate the IR excess of 45 MW Cepheids combining different observables to constrain the presence of CSE. We used the SPIPS algorithm, a robust implementation of the parallax-of-pulsation method that combines photometry, angular diameter, stellar effective temperature, and RV measurements in a global modelling of the pulsation. We obtained new photometric measurements at mid-IR with the VLT/VISIR complemented with literature data. We then compared the mean magnitude from 0.5$\mu$m to 70$\mu$m with stellar atmosphere models to infer the IR excess, which we attribute to the presence of a CSE. We report that at least 29% of our sample have a detected IR excess. We estimated a mean excess of 0.08mag in K and 0.13mag in N. Other Cepheids possibly have IR excess, but they were rejected due to their low detection level compared to a single-star model. We do not see any correlation between the IR excess and the pulsation period as previously suspected for MW Cepheids, but a rather constant trend. We also do not find any correlation between the CO absorption and the presence of a CSE, but rather with the stellar effective temperature, which confirms that the CO features previously reported are mostly photospheric. No bias caused by the presence of the CSE is detected on the average distance estimates from a SPIPS analysis with a fitted colour excess. We also do not find correlation between the presence of IR excess and the evolution stage. We report a fraction of 29% of Cepheids with an IR excess likely produced by the CSE. Longer period Cepheids do not exhibit greater excess than short periods as previously suspected from observations and theoretical dusty-wind models. Other mechanisms such as free-free emission, among others, may be at the origin of their formation. We also show that not fitting the colour excess leads to a bias on the distance estimates in our Galaxy., Comment: Accepted for publication in A&A

The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. We first use photometric observations in various bands and Spitzer Space Telescope spectroscopy to constrain the IR excess of $\ell$ Car. Then, we analyze the VLTI/MATISSE measurements at a specific phase of observation, in order to determine the flux contribution, the size and shape of the environment of the star in the L band. We finally test the hypothesis of a shell of ionized gas in order to model the IR excess. We report the first detection in the L band of a centro-symmetric extended emission around l Car, of about 1.7$R_\star$ in FWHM, producing an excess of about 7.0\% in this band. In the N band, there is no clear evidence for dust emission from VLTI/MATISSE correlated flux and Spitzer data. On the other side, the modeled shell of ionized gas implies a more compact CSE ($1.13\pm0.02\,R_\star$) and fainter (IR excess of 1\% in the L band). We provide new evidences for a compact CSE of $\ell$ Car and we demonstrate the capabilities of VLTI/MATISSE for determining common properties of CSEs. While the compact CSE of $\ell$ Car is probably of gaseous nature, the tested model of a shell of ionized gas is not able to simultaneously reproduce the IR excess and the interferometric observations. Further Galactic Cepheids observations with VLTI/MATISSE are necessary for determining the properties of CSEs, which may also depend on both the pulsation period and the evolutionary state of the stars., Comment: 13 pages, 8 figures, accepted in Astronomy and Astrophysics

The surface brightness -- colour relation (SBCR) is a basic tool in establishing precise and accurate distances within the Local Group. Detached eclipsing binary stars with accurately determined radii and trigonometric parallaxes allow for a calibration of the SBCRs with unprecedented accuracy. We analysed four nearby eclipsing binary stars containing late F-type main sequence components: AL Ari, AL Dor, FM Leo and BN Scl. We determined very precise spectroscopic orbits and combined them with high precision ground- and space-based photometry. We derived the astrophysical parameters of their components with mean errors of 0.1% for mass and 0.4% for radius. We combined those four systems with another 24 nearby eclipsing binaries with accurately known radii from the literature for which $Gaia$ EDR3 parallaxes are available, in order to derive the SBCRs. The resulting SBCRs cover stellar spectral types from B9 V to G7 V. For calibrations we used Johnson optical $B$ and $V$, $Gaia$ $G_{\rm BP}$ and $G$ and 2MASS $JHK$ bands. The most precise relations are calibrated using the infrared $K$ band and allow to predict angular diameters of A-, F-, and G-type dwarf and subgiant stars with a precision of 1%., Comment: Accepted for publication in Astronomy & Astrophysics, 17 pages

The current tension between the direct and the early Universe measurements of the Hubble Constant, $H_0$, requires detailed scrutiny of all the data and methods used in the studies on both sides of the debate. The Cepheids in the type Ia supernova (SNIa) host galaxy NGC 5584 played a key role in the local measurement of $H_0$. The SH0ES project used the observations of this galaxy to derive a relation between Cepheids' periods and ratios of their amplitudes in different optical bands of the Hubble Space Telescope (HST), and used these relations to analyse the light curves of the Cepheids in around half of the current sample of local SNIa host galaxies. In this work, we present an independent detailed analysis of the Cepheids in NGC 5584. We employ different tools for our photometric analysis and a completely different method for our light curve analysis, and we do not find a systematic difference between our period and mean magnitude measurements compared to those reported by SH0ES. By adopting a period-luminosity relation calibrated by the Cepheids in the Milky Way, we measure a distance modulus $\mu=31.810\pm0.047$ (mag) which is in agreement with $\mu=31.786\pm0.046$ (mag) measured by SH0ES. In addition, the relations we find between periods and amplitude ratios of the Cepheids in NGC 5584 are significantly tighter than those of SH0ES and their potential impact on the direct $H_0$ measurement will be investigated in future studies., Comment: Accepted for publication in The Astrophysical Journal on 17-Feb-2021

Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. We use geometric models to fit the data. Our models include a smoothed ring, a flat disk with inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydro-dynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is $3.5 \pm 0.2$. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the $r<0.3$ au disk region is time-variable. We propose that this asymmetric structure, located in or near the inner rim of the dusty disk, orbits the star. For the physical origin of the asymmetry, we tested a hypothesis where a vortex is created by Rossby wave instability, and we find that a unique large scale vortex may be compatible with our data. The half-light radius of the L-band emitting region is $0.33\pm 0.01$ au, the inclination is ${52^\circ}^{+5^\circ}_{-7^\circ}$, and the position angle is $143^\circ \pm 3^\circ$. Our models predict that a non-negligible fraction of the L-band disk emission originates inside the dust sublimation radius for $\mu$m-sized grains. Refractory grains or large ($\gtrsim 10\ \mu$m-sized) grains could be the origin for this emission., Comment: accepted for publication in A&A

We present a new study of late-type eclipsing binary stars in the Small Magellanic Cloud (SMC) undertaken with the aim of improving the distance determination to this important galaxy. A sample of 10 new detached, double-lined eclipsing binaries indentified from the OGLE variable star catalogues and consisting of F- and G-type giant components has been analysed. The absolute physical parameters of the individual components have been measured with a typical accuracy of better than 3%. All but one of the systems consist of young and intermediate population stars with masses in the range of 1.4 to 3.8 M_Sun. This new sample has been combined with five SMC eclipsing binaries previously published by our team. Distances to the binary systems were calculated using a surface brightness - color calibration. The targets form an elongated structure, highly inclined to the plane of the sky. The distance difference between the nearest and most-distant system amounts to 10 kpc with the line of sight depth reaching 7 kpc. We find tentative evidence of the existence of a spherical stellar sub-structure (core) in the SMC coinciding with its stellar center, containing about 40% of the young and intermediate age stars in the galaxy. The radial extension of this sub-structure is ~1.5 kpc. We derive a distance to the SMC center of D_SMC=62.44 +/- 0.47 (stat.) +/- 0.81 (syst.) kpc corresponding to a distance modulus (m-M)_SMC=18.977 +/- 0.016 +/- 0.028 mag, representing an accuracy of better than 2%., Comment: 18 pages, accepted to the Astrophysical Journal

Cepheids in multiple systems provide information on the outcome of the formation of massive stars. They can also lead to exotic end-stage objects. This study concludes our survey of 70 galactic Cepheids using the {\it Hubble Space Telescope\} (\HST) Wide Field Camera~3 (WFC3) with images at two wavelengths to identify companions closer than $5\arcsec$. In the entire WFC3 survey we identify 16 probable companions for 13 Cepheids. The seven Cepheids having resolved candidate companions within $2"$ all have the surprising property of themselves being spectroscopic binaries (as compared with a 29\% incidence of spectroscopic binaries in the general Cepheid population). That is a strong suggestion that an inner binary is linked to the scenario of a third companion within a few hundred~AU\null. This characteristic is continued for more widely separated companions. Under a model where the outer companion is formed first, it is unlikely that it can anticipate a subsequent inner binary. Rather it is more likely that a triple system has undergone dynamical interaction, resulting in one star moving outward to its current location. {\it Chandra\} and {\it Gaia\} data as well as radial velocities and \HSTSTIS and {\it IUE\} spectra are used to derive properties of the components of the Cepheid systems. The colors of the companion candidates show a change in distribution at approximately 2000~AU separations, from a range including both hot and cool colors for closer companions, to only low-mass companions for wider separations., Comment: Accepted by ApJ

It has been shown recently that the infrared emission of Cepheids, which is constant over the pulsation cycle, might be due to a pulsating shell of ionized gas of about 15\% of the stellar radius, which could be attributed to the chromospheric activity of Cepheids. The aim of this paper is to investigate the dynamical structure of the chromosphere of Cepheids along the pulsation cycle and quantify its size. We present H$\alpha$ and Calcium Near InfraRed triplet (Ca IR) profile variations using high-resolution spectroscopy with the UVES spectrograph of a sample of 24 Cepheids with a good period coverage from $\approx$ 3 to 60 days. After a qualitative analysis of the spectral lines profiles, we quantify the Van Hoof effect (velocity gradient between the H$\alpha$ and Ca IR) as a function of the period of the Cepheids. Then, we use the Schwarzschild mechanism (a line doubling due to a shock wave) to quantify the size of the chromosphere. We find a significant Van Hoof effect for Cepheids with period larger than $P=10$ days, in particular H$\alpha$ lines are delayed with a velocity gradient up to $\Delta v \approx$30 km/s compared to Ca IR. We find that the size of the chromosphere of long-period Cepheids is of at least $\approx$ 50\% of the stellar radius, which is consistent at first order with the size of the shell made of ionized gas previously found from the analysis of infrared excess. Last, for most of the long-period Cepheids in the sample, we report a motionless absorption feature in the H$\alpha$ line that we attribute to a circumstellar envelope that surrounds the chromosphere. Analyzing the Ca~IR lines of Cepheids is of importance to potentially unbias the period-luminosity relation from their infrared excess, particularly in the context of forthcoming observations from the Radial Velocity Spectrometer (RVS) on board \textit{Gaia}, that could be sensitive to their chromosphere., Comment: 25 pages, 30 figures

The surface brightness - color relationship (SBCR) is a poweful tool for determining the angular diameter of stars from photometry. It was for instance used to derive the distance of eclipsing binaries in the Large Magellanic Cloud (LMC), which led to its distance determination with an accuracy of 1%. We calibrate the SBCR for red giant stars in the 2.1 < V-K < 2.5 color range using homogeneous VEGA/CHARA interferometric data secured in the visible domain, and compare it to the relation based on infrared interferometric observations, which were used to derive the distance to the LMC. Observations of eight G-K giants were obtained with the VEGA/CHARA instrument. The derived limb-darkened angular diameters were combined with a homogeneous set of infrared magnitudes in order to constrain the SBCR. The average precision we obtain on the limb-darkened angular diameters of the eight stars in our sample is 2.4%. For the four stars in common observed by both VEGA/CHARA and PIONIER/VLTI, we find a 1 sigma agreement for the angular diameters. The SBCR we obtain in the visible has a dispersion of 0.04 magnitude and is consistent with the one derived in the infrared (0.018 magnitude). The consistency of the infrared and visible angular diameters and SBCR reinforces the result of 1% precision and accuracy recently achieved on the distance of the LMC using the eclipsing-binary technique. It also indicates that it is possible to combine interferometric observations at different wavelengths when the SBCR is calibrated., Comment: 9 pages, in press in A&A revue

Classical Cepheids provide the foundation for the empirical extragalactic distance ladder. Milky Way Cepheids are the only stars in this class accessible to trigonometric parallax measurements. However, the parallaxes of Cepheids from the second Gaia data release (GDR2) are affected by systematics because of 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 zero-point (ZP) 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 calibrations based on non-Gaia measurements 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 H_0 = 72.8 +/- 1.9 (statistical + systematics) +/- 1.9 (ZP) km/s/Mpc when all Cepheids are considered and H0 = 73.0 +/- 1.9 (statistical + systematics) +/- 1.9 (ZP) km/s/Mpc for fundamental mode pulsators only., Comment: 15 pages, 20 figures. Accepted for publication in A&A

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., Comment: This paper was updated and replaced by arXiv:2006.08763

Context. The precise determinations of stellar mass at $\sim$1% provide important constraints on stellar evolution models. Accurate parallax measurements can also serve as independent benchmarks for the next Gaia data release. Aims. We aim at measuring the masses and distance of binary systems with a precision level better than 1% using a fully geometrical and empirical method. Methods. We obtained the first interferometric observations for the eclipsing systems AI Phe, AL Dor, KW Hya, NN Del, $\psi$ Cen and V4090 Sgr with the VLTI/PIONIER combiner, which we combined with radial velocity measurements to derive their three-dimensional orbit, masses, and distance. Results. We determined very precise stellar masses for all systems, ranging from 0.04% to 3.3 % precision level. We combined these measurements with stellar effective temperature and linear radius to fit stellar isochrones models and determined the age of the systems. We also derived the distance to the systems with a precision level as high as 0.4%. Conclusions. The comparison of theoretical models with stellar parameters shows that stellar models are still deficient in simultaneously fitting the stellar parameters (Teff , R and M) with such level of precision on individual masses. This stresses the importance of precisely measuring the stellar parameters to better calibrate stellar evolution models. The precision of our model-independent orbital parallaxes varies from 24$\mu$as as to 70$\mu$as and they provide a unique opportunity to check on the future Gaia measurements for possible systematic errors., Comment: Accepted for publication in A&A

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., Comment: Accepted in A&A, 21 pages, 17 figures

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., Comment: 24 pages, 18 figures, 4 tables. Formally accepted for publication in Astronomy & Astrophysics on August 5th, 2019

Context. The multiplicity of classical Cepheids (CCs) and RR Lyrae stars (RRLs) is still imperfectly known, particularly for RRLs. Aims. In order to complement the close-in short orbital period systems presented in Paper I, our aim is to detect the wide, spatially resolved companions of the targets of our reference samples of Galactic CCs and RRLs. Methods. Angularly resolved common proper motion pairs were detected using a simple progressive selection algorithm to separate the most probable candidate companions from the unrelated field stars. Results. We found 27 resolved, high probability gravitationally bound systems with CCs out of 456 examined stars, and one unbound star embedded in the circumstellar dusty nebula of the long-period Cepheid RS Pup. We found seven spatially resolved, probably bound systems with RRL primaries out of 789 investigated stars, and 22 additional candidate pairs. We report in particular new companions of three bright RRLs: OV And (companion of F4V spectral type), RR Leo (M0V), and SS Oct (K2V). In addition, we discovered resolved companions of 14 stars that were likely misclassified as RRLs. Conclusions. The detection of resolved non-variable companions around CCs and RRLs facilitates the validation of their Gaia DR2 parallaxes. The possibility to conduct a detailed analysis of the resolved coeval companions of CCs and old population RRLs will also be valuable to progress on our understanding of their evolutionary path., Comment: 42 pages, published in Astronomy & Astrophysics

In the era of precision cosmology, it is essential to empirically determine the Hubble constant with an accuracy of one per cent or better. At present, the uncertainty on this constant is dominated by the uncertainty in the calibration of the Cepheid period - luminosity relationship (also known as Leavitt Law). The Large Magellanic Cloud has traditionally served as the best galaxy with which to calibrate Cepheid period-luminosity relations, and as a result has become the best anchor point for the cosmic distance scale. Eclipsing binary systems composed of late-type stars offer the most precise and accurate way to measure the distance to the Large Magellanic Cloud. Currently the limit of the precision attainable with this technique is about two per cent, and is set by the precision of the existing calibrations of the surface brightness - colour relation. Here we report the calibration of the surface brightness-colour relation with a precision of 0.8 per cent. We use this calibration to determine the geometrical distance to the Large Magellanic Cloud that is precise to 1 per cent based on 20 eclipsing binary systems. The final distane is 49.59 +/- 0.09 (statistical) +/- 0.54 (systematic) kiloparsecs., Comment: Published in Nature 567, 200 (2019)

Classical Cepheids (CCs) and RR Lyrae stars (RRLs) are important classes of variable stars used as standard candles to estimate galactic and extragalactic distances. Their multiplicity is imperfectly known, particularly for RRLs. Astoundingly, to date only one RRL has convincingly been demonstrated to be a binary, TU UMa, out of tens of thousands of known RRLs. Our aim is to detect the binary and multiple stars present in a sample of Milky Way CCs and RRLs. In the present article, we combine the Hipparcos and Gaia DR2 positions to determine the mean proper motion of the targets, and we search for proper motion anomalies (PMa) caused by close-in orbiting companions. We identify 57 CC binaries from PMa out of 254 tested stars and 75 additional candidates, confirming the high binary fraction of these massive stars. For 28 binary CCs, we determine the companion mass by combining their spectroscopic orbital parameters and astrometric PMa. We detect 13 RRLs showing a significant PMa out of 198 tested stars, and 61 additional candidates. We determine that the binary fraction of CCs is likely above 80%, while that of RRLs is at least 7%. The newly detected systems will be useful to improve our understanding of their evolutionary states. The discovery of a significant number of RRLs in binary systems also resolves the long-standing mystery of their extremely low apparent binary fraction., Comment: 18 pages, 6 figures, Astronomy & Astrophysics, in press

Using a sample of 81 galactic, detached eclipsing binary stars we investigated the global zero-point shift of their parallaxes with the Gaia Data Release 2 (DR2) parallaxes. The stars in the sample lay in a distance range of 0.04-2 kpc from the Sun. The photometric parallaxes {\pi}_Phot of the eclipsing binaries were determined by applying a number of empirical surface brightness - color (SBC) relations calibrated on optical-infrared colors. For each SBC relation we calculated the individual differences d{\pi}_i = ({\pi}_Gaia - {\pi}_Phot)_i and then we calculated unweighted and weighted means. As the sample covers the whole sky we interpret the weighted means as the global shifts of the Gaia DR2 parallaxes with respect to our eclipsing binary sample. Depending on the choice of the SBC relation the shifts vary from -0.094 mas to -0.025 mas. The weighted mean of the zero-point shift from all colors and calibrations used is d{\pi} = -0.054 +/- 0.024 mas. However, the SBC relations based on (B-K) and (V-K) colors, which are the least reddening dependent and have the lowest intrinsic dispersions, give a zero-point shift of d{\pi} = -0.031 +/- 0.011 mas in full agreement with results obtained by Lindegren et al. and Arenou et al. Our result confirms the global shift of Gaia DR2 parallaxes of d{\pi} = -0.029 mas reported by the Gaia team, but we do not confirm the larger zero-point shift reported by a number of follow-up papers., Comment: Like original version but with updated grant numbers

We aim at detecting and characterizing the main-sequence companions of a sample of known and suspected Galactic binary Cepheids. We used the multi-telescope interferometric combiners MIRC and PIONIER to detect and measure the astrometric positions of the high-contrast companions orbiting 16 bright Galactic Cepheids. We made use of the CANDID algorithm to search for the companions and set detection limits. We also present new high-precision RVs which were used to fit radial pulsation and orbital velocities. We report the detection of the companions orbiting the Cepheids U Aql, BP Cir, and S Mus, possible detections for FF Aql, Y Car, BG Cru, X Sgr, V350 Sgr, and V636 Sco, while no component is detected around U Car, YZ Car, T Mon, R Mus, S Nor, W Sgr and AH Vel. For U Aql and S Mus, we performed a preliminary orbital fit combining astrometry with new high-precision RVs, providing the full set of orbital and pulsation parameters. Assuming the distance, we estimated preliminary masses of M(U Aql) = 4.97+/-0.62Msun and M(S Mus) = 4.63+/-0.99Msun. For YZ Car, W Sgr, V350 Sgr, and V636 Sco, we revised the spectroscopic orbits using new high-precision RVs, while we updated the pulsation parameters for BP Cir, BG Cru, S Nor and AH Vel. Our interferometric observations also provide measurements of the angular diameters, that can be used in a Baade-Wesselink type analysis. We have now several astrometric detections of Cepheid companions. When RVs of the companions will be available, such systems will provide accurate and independent masses and distances. Orbital parallaxes with an accuracy <5% will be particularly useful for a better calibration of the P-L relation. The final Gaia parallaxes will be also particularly helpful for single-line spectroscopic systems, where mass and distance are degenerate. Mass measurements are necessary for a better understanding of the age and evolution of Cepheids, Comment: 29 pages, 11 figures, 23 tables. Accepted for publication in A&A

Cepheid stars play a considerable role as extragalactic distances indicators, thanks to the simple empirical relation between their pulsation period and their luminosity. They overlap with that of secondary distance indicators, such as Type Ia supernovae, whose distance scale is tied to Cepheid luminosities. However, the Period-Luminosity (P-L) relation still lacks a calibration to better than 5%. Using an original combination of interferometric astrometry with optical and ultraviolet spectroscopy, we measured the geometrical distance d = 720.35+/-7.84 pc of the 3.33 d period Cepheid V1334 Cyg with an unprecedented accuracy of +/-1 %, providing the most accurate distance for a Cepheid. Placing this star in the P-L diagram provides an independent test of existing period-luminosity relations. We show that the secondary star has a significant impact on the integrated magnitude, particularly at visible wavelengths. Binarity in future high precision calibrations of the P-L relations is not negligible, at least in the short-period regime. Subtracting the companion flux leaves V1334 Cyg in marginal agreement with existing photometric-based P-L relations, indicating either an overall calibration bias or a significant intrinsic dispersion at a few percent level. Our work also enabled us to determine the dynamical masses of both components, M1 = 4.288 +/- 0.133 Msun (Cepheid) and M2 = 4.040 +/- 0.048 Msun (companion), providing the most accurate masses for a Galactic binary Cepheid system., Comment: Accepted for publication in ApJ

The extragalactic distance scale builds on the Cepheid period-luminosity (PL) relation. In this paper, we want to carry out a strictly differential comparison of the absolute PL relations obeyed by classical Cepheids in the Milky Way (MW), LMC and SMC galaxies. Taking advantage of the substantial metallicity difference among the Cepheid populations in these three galaxies, we want to establish a possible systematic trend of the PL relation absolute zero point as a function of metallicity, and determine the size of such an effect in optical and near-infrared photometric bands. We are using the IRSB Baade-Wesselink type method as calibrated by Storm et al. to determine individual distances to the Cepheids in our samples in MW, LMC and SMC. For our analysis, we use a greatly enhanced sample of Cepheids in the SMC (31 stars) as compared to the small sample (5 stars) available in our previous work. We use the distances to determine absolute Cepheid PL relations in optical and near-infrared bands in each of the three galaxies.} {Our distance analysis of 31 SMC Cepheids with periods from 4-69 days yields tight PL relations in all studied bands, with slopes consistent with the corresponding LMC and MW relations. Adopting the very accurately determined LMC slopes for the optical and near-infrared bands, we determine the zero point offsets between the corresponding absolute PL relations in the 3 galaxies. We find that in all bands the metal-poor SMC Cepheids are intrinsically fainter than their more metal-rich counterparts in the LMC and MW. In the $K$ band the metallicity effect is $-0.23\pm0.06$~mag/dex while in the $V,(V-I)$ Wesenheit index it is slightly stronger, $-0.34\pm0.06$~mag/dex. We find some evidence that the PL relation zero point-metallicity relation might be nonlinear, becoming steeper for lower metallicities., Comment: 20 pages, accepted for publication by A&A

While planet formation is thought to occur early in the history of a protoplanetary disk, the presence of planets embedded in disks, or of other processes driving disk evolution, might be traced from their imprints on the disk structure. We observed the T Tauri star HD 143006, located in the 5-11 Myr-old Upper Sco region, in polarized scattered light with VLT/SPHERE at near-infrared wavelengths, reaching an angular resolution of ~0.037" (~6 au). We obtained two datasets, one with a 145 mas diameter coronagraph, and the other without, enabling us to probe the disk structure down to an angular separation of ~0.06" (~10 au). In our observations, the disk of HD 143006 is clearly resolved up to ~0.5" and shows a clear large-scale asymmetry with the eastern side brighter than the western side. We detect a number of additional features, including two gaps and a ring. The ring shows an overbrightness at a position angle (PA) of ~140 deg, extending over a range in position angle of ~60 deg, and two narrow dark regions. The two narrow dark lanes and the overall large-scale asymmetry are indicative of shadowing effects, likely due to a misaligned inner disk. We demonstrate the remarkable resemblance between the scattered light image of HD 143006 and a model prediction of a warped disk due to an inclined binary companion. The warped disk model, based on the hydrodynamic simulations combined with 3D radiative transfer calculations, reproduces all major morphological features. However, it does not account for the observed overbrightness at PA~140 deg. Shadows have been detected in several protoplanetary disks, suggesting that misalignment in disks is not uncommon. However, the origin of the misalignment is not clear. As-yet-undetected stellar or massive planetary companions could be responsible for them, and naturally account for the presence of depleted inner cavities., Comment: Accepted for publication in A&A

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., Comment: A&A, in press

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., Comment: 4 pages, 3 figures, Proceedings of the IAU Symposium 330 "Astrometry and Astrophysics in the Gaia sky" (Nice, France, 24-28 April 2017), pp. 305-308

Hi-5 is a high-contrast (or high dynamic range) infrared imager project for the VLTI. Its main goal is to characterize young extra-solar planetary systems and exozodiacal dust around southern main-sequence stars. In this paper, we present an update of the project and key technology pathways to improve the contrast achieved by the VLTI. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5~$\mu$m) and advanced fringe tracking strategies. We also address the strong exoplanet science case (young exoplanets, planet formation, and exozodiacal disks) offered by this wavelength regime as well as other possible science cases such as stellar physics (fundamental parameters and multiplicity) and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed., Comment: 15 pages, 3 figures, SPIE proceedings, based on arXiv:1801.04148

As part of a program to determine dynamical masses of Cepheids, we have imaged the nearest and brightest Cepheid, Polaris, with the Hubble Space Telescope Wide Field Planetary Camera 2 and Wide Field Camera 3. Observations were obtained at three epochs between 2007 and 2014. In these images, as in HST frames obtained in 2005 and 2006, which we discussed in a 2008 paper, we resolve the close companion Polaris Ab from the Cepheid Polaris Aa. Because of the small separation and large magnitude difference between Polaris Aa and Ab, we used PSF deconvolution techniques to carry out astrometry of the binary. Based on these new measurements, we have updated the elements for the 29.59 yr orbit. Adopting the distance to the system from the recent Gaia Data Release 2, we find a dynamical mass for the Cepheid of 3.45 +/- 0.75 Msun, although this is preliminary, and will be improved by CHARA measurements covering periastron. As is the case for the recently determined dynamical mass for the Cepheid V1334 Cyg, the mass of Polaris is significantly lower than the "evolutionary mass" predicted by fitting to evolutionary tracks in the HR diagram. We discuss several questions and implications raised by these measurements, including the pulsation mode, which instability-strip crossing the stars are in, and possible complications such as rotation, mass loss, and binary mergers. The distant third star in the system, Polaris B, appears to be older than the Cepheid, based on isochrone fitting. This may indicate that the Cepheid Polaris is relatively old and is the result of a binary merger, rather than being a young single star., Comment: Accepted by Astrophysical Journal

Observations of 48 red-clump stars were obtained in the H band with the PIONIER instrument installed at the Very Large Telescope Interferometer. Limb-darkened angular diameters were measured by fitting radial intensity profile I(r) to square visibility measurements. Half the angular diameters determined have formal errors better than 1.2%, while the overall accuracy is better than 2.7%. Average stellar atmospheric parameters (effective temperatures, metallicities and surface gravities) were determined from new spectroscopic observations and literature data and combined with precise Gaia parallaxes to derive a set of fundamental stellar properties. These intrinsic parameters were then fitted to existing isochrone models to infer masses and ages of the stars. The added value from interferometry imposes a better and independent constraint on the R-Teff plane. Our derived values are consistent with previous works, although there is a strong scatter in age between various models. This shows that atmospheric parameters, mainly metallicities and surface gravities, still suffer from a non-accurate determination, limiting constraints on input physics and parameters of stellar evolution models., Comment: Accepted for publication in A&A, 12 pages