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.

Red supergiants are cool massive stars and are the largest and the most luminous stars in the universe. They are characterized by irregular or semi-regular photometric variations, the physics of which is not clearly understood. The paper aims at deriving the velocity field in the red supergiant star $\mu$ Cep and relating it to the photometric variability with the help of the tomographic method. The tomographic method allows to recover the line-of-sight velocity distribution over the stellar disk and within different optical-depth slices. The method is applied to a series of high-resolution spectra of $\mu$ Cep, and these results are compared to those obtained from 3D radiative-hydrodynamics CO5BOLD simulations of red supergiants. Fluctuations in the velocity field are compared with photometric and spectroscopic variations, the latter being derived from the TiO band strength and serving (at least partly) a proxy of the variations in effective temperature. The tomographic method reveals a phase shift between the velocity and spectroscopic/photometric variations. This phase shift results in a hysteresis loop in the temperature - velocity plane, with a timescale of a few hundred days, similar to the photometric one. The similarity between the hysteresis loop timescale measured in $\mu$ Cep and the timescale of acoustic waves disturbing the convective pattern suggests that such waves play an important role in triggering the hysteresis loops.

A tomographic method, aiming at probing velocity fields at depth in stellar atmospheres, is applied to the red supergiant star {\mu} Cep and to snapshots of 3D radiative-hydrodynamics simulation in order to constrain atmospheric motions and relate them to photometric variability.

Cool giant and supergiant star atmospheres are characterized by complex velocity fields originating from convection and pulsation processes which are not fully understood yet. The velocity fields impact the formation of spectral lines, which thus contain information on the dynamics of stellar atmospheres. The tomographic method allows to recover the distribution of the component of the velocity field projected on the line of sight at different optical depths in the stellar atmosphere. The computation of the contribution function to the line depression aims at correctly identifying the depth of formation of spectral lines in order to construct numerical masks probing spectral lines forming at different optical depths. The tomographic method is applied to 1D model atmospheres and to a realistic 3D radiative hydrodynamics simulation performed with CO5BOLD in order to compare their spectral line formation depths and velocity fields. In 1D model atmospheres, each spectral line forms in a restricted range of optical depths. On the other hand, in 3D simulations, the line formation depths are spread in the atmosphere mainly because of temperature and density inhomogeneities. Comparison of CCF profiles obtained from 3D synthetic spectra with velocities from the 3D simulation shows that the tomographic method correctly recovers the distribution of the velocity component projected on the line of sight in the atmosphere.