Heading Gaia to measure atmospheric dynamics in AGB stars

Context. Asymptotic giant branch (AGB) stars are characterised by complex stellar surface dynamics that affect the measurements and amplify the uncertainties on stellar parameters. The uncertainties in observed absolute magnitudes have been found to originate mainly from uncertainties in the parallaxes. The resulting motion of the stellar photocentre could have adverse effects on the parallax determination with Gaia. Aims. We explore the impact of the convection-related surface structure in AGBs on the photocentric variability. We quantify these effects to characterise the observed parallax errors and estimate fundamental stellar parameters and dynamical properties. Methods. We use three-dimensional (3D) radiative hydrodynamics simulations of convection with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps in the Gaia G band [325–1030 nm]. From those maps, we calculate the intensity-weighted mean of all emitting points tiling the visible stellar surface (i.e. the photocentre) and evaluate its motion as a function of time. We extract the parallax error from Gaia data-release 2 (DR2) for a sample of semi-regular variables in the solar neighbourhood and compare it to the synthetic predictions of photocentre displacements. Results. AGB stars show a complex surface morphology characterised by the presence of few large-scale long-lived convective cells accompanied by short-lived and small-scale structures. As a consequence, the position of the photocentre displays temporal excursions between 0.077 and 0.198 AU (≈5 to ≈11% of the corresponding stellar radius), depending on the simulation considered. We show that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error of our sample of semi-regular variables. Finally, we present evidence for a correlation between the mean photocentre displacement and the stellar fundamental parameters: surface gravity and pulsation. We suggest that parallax variations could be exploited quantitatively using appropriate radiation-hydrodynamics (RHD) simulations corresponding to the observed star