Search

Classical Cepheids are a cornerstone class of pulsators, fundamental to testing stellar evolution and pulsation theories. Their secular period changes, characterized through $O-C$ (Observed minus Calculated) diagrams, offer valuable insights into their evolution. While evolutionary period changes are well understood from both observational and theoretical perspectives, shorter timescale period changes (on the order of ($\sim$ 10$^{2}$-10$^{4}$ days) - known as non-evolutionary period changes are yet to be systematically explored. In this work, we present a detailed and comprehensive search for non-evolutionary period changes using $O-C$ analysis of Magellanic Cloud (MC) Cepheids, based on 20+ years of OGLE photometry data. Our sample includes both the Large Magellanic Cloud (LMC) and the Small Magellanic Cloud (SMC) Cepheids, focusing on single radial mode Cepheids (both fundamental (FU) and first overtone (FO) modes). The results are grouped into two phenomena: (a) Cepheids in binary systems (b) Non-linear period changes., Comment: 7 pages, 4 figures. Submitted to Contributions of the Astronomical Observatory Skalnate Pleso as proceedings for "Binary and Multiple Stars in the Era of Big Sky Surveys"

The circumstellar envelopes (CSE) of Cepheids are still not well characterized despite their potential impact on distance determination via both the period-luminosity relation and the parallax-of-pulsation method. This paper aims to investigate Galactic Cepheids across the instability strip in the mid-infrared with MATISSE/VLTI in order to constrain the geometry and physical nature (gas and/or dust) of their CSEs. We secured observations of eight Galactic Cepheids from short up to long period of pulsation, with MATISSE/VLTI in $L$, $M$ and $N$-bands. The SED analysis in the mid-IR confirms the absence of dust spectral signature for all the star sample. For each star in $L$, $M$ and $N$-band we observe closure phases which are consistent with centro-symmetric geometry for the different targets. Finally, the visibilities in $L$, $M$ and $N$ bands are in agreement with the expected star angular diameter, although the observations are compatible with the presence of compact CSEs within the uncertainties. We provide 2$\,\sigma$ upper limits on the CSE flux contribution based on model residuals for several CSE radius, which yield to exclude models simultaneously large and bright ($R_\mathrm{CSE}\approx10\,R_\star$ and $f_\mathrm{CSE}\approx10\%$) for all the stars of the sample. Last, the visibilities in the $N$-band rule out CSE models with significant amount of different type of dust. The MATISSE observations of eight Cepheids with different pulsation period (from 7 up to 38$\,$day) and evolution stage, provide for the first time a comprehensive picture of Cepheids from mid-IR interferometry. We present additional evidences that circumstellar dust emission is negligible or absent around Cepheids for a wide range of stellar parameters in the instability strip. Further interferometric observations in the visible and the near-infrared will be necessary to disentangle the star and the CSE., Comment: 22 pages, 17 figures, accepted in Astronomy and Astrophysics

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

Y~Ophiuchi (Y~Oph) is a classical Cepheid reported to be as dim as a Cepheid of about half its pulsation period, and exhibits a low radial velocity and light-curves amplitude. Our objective is to conduct hydrodynamical pulsation modeling of Y~Oph to derive its distance and provide physical insight to its low amplitude and luminosity, constrained by an extensive set of observations. We first perform a linear analysis on a grid of models using hydrodynamical pulsation code \texttt{MESA-RSP} in order to find the combinations of mass, metallicity, effective temperature and luminosity resulting in linear excitation of pulsations with period of about 17$\,$days. Then, for the best combinations of stellar parameters, we perform non-linear computations to obtain the full-amplitude pulsations of these models. Last, we compare the results to a complete set of observations along the pulsation cycle. We adjust simultaneously the distance, the color excess and circumstellar envelope (CSE) model to fit the light curves and the angular diameter. We find that all pulsation models at high effective temperature are in remarkable agreement with the observations along the pulsation cycle. This result suggests that the low amplitude of Y~Oph can be explained by its location close to the blue edge of the instability strip. We also find that a pulsational mass of about 7-8$\,\mathrm{M}_\odot$ is consistent with a non-canonical evolutionary model with moderate overshooting, PL relation and \textit{Gaia} parallax. However, a much lower mass below 5$\,$M$_\odot$ is required to match Baade-Wesselink (BW) distance measurements from the literature. We show that the combination of the impact of the CSE on the photometry together with a projection factor of about 1.5 explains the discrepant distance and luminosity values obtained from BW methods., Comment: Accepted in A&A. 23 pages, 15 figures

The primary goal of this paper is to derive precise Fourier parameters of the radial velocity (RV) curves for fundamental and first-overtone Galactic Cepheids. For each star, we carefully selected RV measurements available in the literature which yield the highest precision of Fourier parameters. We performed a Fourier decomposition of the RV curves. We subtracted RV modulations caused by binary motion and have removed other residual trends. Finally, we have displayed and analyzed qualitatively the progressions found for Fourier parameters. We applied a standard identification of their pulsation mode based on their Fourier phase $\phi_{21}$. Our final sample includes 178 fundamental-mode and 33 first overtone pulsators, as well as 7 additional Cepheids whose pulsation mode is uncertain or undetermined according to our criteria. For the fundamental-mode Cepheids, we improved the precision of Fourier parameters in comparison of previous results from the literature. We are able to firmly identify V495 Cyg as a new first-overtone Cepheid. We confirm first-overtone nature of several other stars. We also show that $\alpha$ UMi should be definitely classified as a first-overtone pulsator. In 3 objects (VY Per, AQ Pup and QZ Nor) we found significant $\gamma$-velocity variations that we attribute to spectroscopic binarity. Finally, the analysis of the F-mode Fourier parameters up to 7th order reveals tight progression of Fourier phases. For $P<10\,$day we find a well defined upper limit for the Fourier amplitude $A_1$. The pulsation period coverage and the precision obtained, in particular for Fourier phase $\phi_{21}$, will be useful for studying the dynamics of Cepheid pulsations with the help of hydrodynamical models. Further radial velocity measurements from modern high-resolution spectroscopic instruments will be important to improve these results., Comment: 37 pages, 7 figures. Revised version. Originally submitted to A&A on 25/08/2023

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

Estimating metallicity of classical Cepheids is of prime importance for studying metallicity effect on stellar evolution, chemical evolution of galaxies, and ultimately its impact on period-luminosity relation used in the extragalactic distance scale. We aim at establishing new empirical relations for estimating the iron content of classical Cepheids for short and long-periods based on Fourier parameters from the $V$-band light curves. We calibrate new interrelations of Fourier parameters to convert $V$-band empirical relations into the $I$-band. Then we apply these relation in $V$ and $I$-bands to Cepheids from Milky Way (MW), Small and Large Magellanic Clouds (SMC and LMC) available in the literature. Last, we map the metallicity distribution in these galaxies for investigating potential application in galactic archeology. These empirical relations in $V$ and $I$ bands are able to derive the mean metallicity of a sample of MW, SMC and LMC Cepheids in agreement with literature values within 1$\sigma$. We also show that these relations are precise enough to reconstruct the radial metallicity gradients within the MW from OGLE data. The empirical relations in the $V$ and $I$ bands calibrated in this paper for short and long-period Cepheids provide a new useful tool to estimate the metallicity of Cepheids which are not accessible by spectroscopy. The calibration can be improved with further high-resolution spectroscopic observations of metal-poor Cepheids and homogeneous photometries in $V$ and $I$ bands., Comment: Accepted in A&A, 20 pages, 20 figures

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 {\mu}m), M (4.5 to 5 {\mu}m), and N (8 to 13 {\mu}m) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics (RHD) simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes Optim3D and Radmc3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of the stellar surface and circumstellar environment across a very large spectral domain for the first time., Comment: Accepted for publication on Astronomy and Astrophysics

Context: FS Canis Majoris (FS CMa, HD 45677) is an unclassified B[e] star surrounded by an inclined dust disk. The evolutionary stage of FS CMa is still debated. Perpendicular to the circumstellar disk, a bipolar outflow was detected. Infrared aperture-synthesis imaging provides us with a unique opportunity to study the disk structure. Aims: Our aim is to study the intensity distribution of the disk of FS CMa in the mid-infrared L and N bands. Methods: We performed aperture-synthesis imaging of FS CMa with the MATISSE instrument (Multi AperTure mid-Infrared SpectroScopic Experiment) in the low spectral resolution mode to obtain images in the L and N bands. We computed radiative transfer models that reproduce the L- and N-band intensity distributions of the resolved disks. Results: We present L- and N-band aperture-synthesis images of FS CMa reconstructed in the wavelength bands of 3.4-3.8 and 8.6-9.0 micrometer. In the L-band image, the inner rim region of an inclined circumstellar disk and the central object can be seen with a spatial resolution of 2.7 milliarcsec (mas). An inner disk cavity with an angular diameter of 6x12mas is resolved. The L-band disk consists of a bright northwestern (NW) disk region and a much fainter southeastern (SE) region. The images suggest that we are looking at the bright inner wall of the NW disk rim, which is on the far side of the disk. In the N band, only the bright NW disk region is seen. In addition to deriving the inclination and the inner disk radius, fitting the reconstructed brightness distributions via radiative transfer modeling allows one to constrain the innermost disk structure, in particular the shape of the inner disk rim., Comment: accepted for publication in A&A

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

Carbon is one of the most abundant components in the Universe. While silicates have been the main focus of solid phase studies in protoplanetary discs (PPDs), little is known about the solid carbon content especially in the planet-forming regions ($\sim $0.1 to 10 au). Fortunately, several refractory carbonaceous species present C-H bonds (such as hydrogenated nano-diamond and amorphous carbon as well as polycyclic aromatic hydrocarbons (PAHs)), which generate infrared (IR) features that can be used to trace the solid carbon reservoirs. The new mid-IR instrument MATISSE, installed at the Very Large Telescope Interferometer (VLTI), can spatially resolve the inner regions ($\sim$ 1 to 10 au) of PPDs and locate, down to the au-scale, the emission coming from carbon grains. Our aim is to provide a consistent view on the radial structure, down to the au-scale, as well as basic physical properties and the nature of the material responsible for the IR continuum emission in the inner disk region around HD 179218. We implemented a temperature-gradient model to interpret the disk IR continuum emission, based on a multiwavelength dataset comprising a broadband spectral energy distribution (SED) and VLTI H-, L-, and N-bands interferometric data obtained in low spectral resolution. Then, we added a ring-like component, representing the carbonaceous L-band features-emitting region, to assess its detectability in future higher spectral resolution observations employing mid-IR interferometry.

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

Surface brightness-colour relations (SBCRs) are used to derive the stellar angular diameters from photometric observations. They have various astrophysical applications, such as the distance determination of eclipsing binaries or the determination of exoplanet parameters. However, strong discrepancies between the SBCRs still exist in the literature, in particular for early and late-type stars. We aim to calibrate new SBCRs as a function of the spectral type and the luminosity class of the stars. Our goal is also to apply homogeneous criteria to the selection of the reference stars and in view of compiling an exhaustive and up-to-date list of interferometric late-type targets. We implemented criteria to select measurements in the JMMC Measured Diameters Catalog (JMDC). We then applied additional criteria on the photometric measurements used to build the SBCRs, together with stellar characteristics diagnostics. We built SBCRs for F5/K7-II/III, F5/K7-IV/V, M-II/III and M-V stars, with respective RMS of $\sigma_{F_{V}} = 0.0022$ mag, $\sigma_{F_{V}} = 0.0044$ mag, $\sigma_{F_{V}} = 0.0046$ mag, and $\sigma_{F_{V}} = 0.0038$ mag. This results in a precision on the angular diameter of 1.0\%, 2.0\%, 2.1\%, and 1.7\%, respectively. These relations cover a large $V-K$ colour range of magnitude, from 1 to 7.5. Our work demonstrates that SBCRs are significantly dependent on the spectral type and the luminosity class of the star. Through a new set of interferometric measurements, we demonstrate the critical importance of the selection criteria proposed for the calibration of SBCR. Finally, using the Gaia photometry for our samples, we obtained (G-K) SBCRs with a precision on the angular diameter between 1.1\% and 2.4\%., Comment: 19 pages

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

Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations., Comment: Accepted for publication on Astronomy and Astrophysics

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