Convection-related Surface Structures Research Articles

Probing red supergiant dynamics through photo-center displacements measured by Gaia

Context. Red supergiant (RSGs) are cool massive stars in a late phase of their evolution when the stellar envelope becomes fully convective. They are the brightest stars in the universe at infrared light and can be detected in galaxies far beyond the Local Group, allowing for accurate determination of chemical composition of galaxies. The study of their physical properties is extremely important for various phenomena including the final fate of massive stars as type II supernovae and gravitational wave progenitors. Aims. We explore the well-studied nearby young stellar cluster χ Per, which contains a relatively large population of RSG stars. Using Gaia EDR3 data, we find the distance of the cluster (d = 2.260 ± 0.020 kpc) from blue main sequence stars and compare with RSG parallax measurements analysing the parallax uncertainties of both groups. We then investigate the variability of the convection-related surface structure as a source for parallax measurement uncertainty. Methods. We use state-of-the-art three-dimensional radiative hydrodynamics simulations of convection with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps in the Gaia G photometric system. We calculate the variabiltiy, as a function of time, of the intensity-weighted mean (or the photo-center) from the synthetic maps. We then select the RSG stars in the cluster and compare their uncertainty on parallaxes to the predictions of photocentre displacements. Results. The synthetic maps of RSG show extremely irregular and temporal variable surfaces due to convection-related dynamics. Consequentially, the position of the photo-center varies during Gaia measurements between 0.033 and 0.130 AU (≈1 to ≈5% of the corresponding simulation stellar radius). We argue that the variability of the convection-related surface structures accounts for a substantial part of the Gaia EDR3 parallax error of the RSG sample of χ Per. Conclusions. We suggest that the variation of the uncertainty on Gaia parallax could be exploited quantitatively using appropriate RHD simulations to extract, in a unique way, important information about the stellar dynamics and parameters of RSG stars.

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

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–4 μm), M (4.5–5 μm), and N (8–13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics 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 VX Sgr’s stellar surface and circumstellar environment across a very large spectral domain for the first time.

Optical interferometry and Gaia measurement uncertainties reveal the physics of asymptotic giant branch stars

Context. Asymptotic giant branch (AGB) stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics, which amplifies the uncertainties on stellar parameters and distances. 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 (RHD) 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 using the average 3D intensity profile and, finally, 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 μm, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths (3.299 ± 0.005 mas and 3.053 ± 0.006 mas at 1.52 and 1.70 μm, respectively) and recovered the radius (R = 307 ± 41 and R = 284 ± 38 R⊙) 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.

Measuring stellar granulation during planet transits

Context. Stellar activity and convection-related surface structures might cause bias in planet detection and characterization that use these transits. Surface convection simulations help to quantify the granulation signal. Aims. We used realistic three-dimensional (3D) radiative hydrodynamical (RHD) simulations from the Stagger grid and synthetic images computed with the radiative transfer code Optim3D to model the transits of three prototype planets: a hot Jupiter, a hot Neptune, and a terrestrial planet. Methods. We computed intensity maps from RHD simulations of the Sun and a K-dwarf star at different wavelength bands from optical to far-infrared that cover the range of several ground- and space-based telescopes which observe exoplanet transits. We modeled the transit using synthetic stellar-disk images obtained with a spherical-tile imaging method and emulated the temporal variation of the granulation intensity generating random images covering a granulation time-series of 13.3 h. We measured the contribution of the stellar granulation on the light curves during the planet transit. Results. We identified two types of granulation noise that act simultaneously during the planet transit: (i) the intrinsic change in the granulation pattern with timescale (e.g., 10 min for solar-type stars assumed in this work) is smaller than the usual planet transit (~hours as in our prototype cases); and (ii) the fact that the transiting planet occults isolated regions of the photosphere that differ in local surface brightness as a result of convective-related surface structures. First, we showed that our modeling approach returns granulation timescale fluctuations that are comparable with what has been observed for the Sun. Then, our statistical approach shows that the granulation pattern of solar and K-dwarf-type stars have a non-negligible effect of the light curve depth during the transit, and, consequentially on the determination of the planet transit parameters such as the planet radius (up to 0.90% and ~0.47% for terrestrial and gaseous planets, respectively). We also showed that larger (or smaller) orbital inclination angles with respect to values corresponding to transit at the stellar center display a shallower transit depth and longer ingress and egress times, but also granulation fluctuations that are correlated to the center-to-limb variation: they increase (or decrease) the value of the inclination, which amplifies the fluctuations. The granulation noise appears to be correlated among the different wavelength ranges either in the visible or in the infrared regions. Conclusions. The prospects for planet detection and characterization with transiting methods are excellent with access to large amounts of data for stars. The granulation has to be considered as an intrinsic uncertainty (as a result of stellar variability) on the precise measurements of exoplanet transits of planets. The full characterization of the granulation is essential for determining the degree of uncertainty on the planet parameters. In this context, the use of 3D RHD simulations is important to measure the convection-related fluctuations. This can be achieved by performing precise and continuous observations of stellar photometry and radial velocity, as we explained with RHD simulations, before, after, and during the transit periods.

New view on exoplanet transits

Context. An important benchmark for current observational techniques and theoretical modeling of exoplanet atmospheres is the transit of Venus (ToV). Stellar activity and, in particular, convection-related surface structures, potentially cause fluctuations that can affect the transit light curves. Surface convection simulations can help interpreting the ToV as well as other transits outside our solar system.

Planet transit and stellar granulation detection with interferometry

Context. Stellar activity and, in particular, convection-related surface structures, potentially cause bias in planet detection and characterisation. In the latter, interferometry can help disentangle the signal of the transiting planet.

Stellar granulation and interferometry

Stars are not smooth. Their photosphere is covered by a granulation pattern associated with the heat transport by convection. The convection-related surface structures have different size, depth, and temporal variations with respect to the stellar type. The related activity (in addition to other phenomena such as magnetic spots, rotation, dust, etc.) potentially causes bias in stellar parameters determination, radial velocity, chemical abundances determinations, and exoplanet transit detections. The role of long-baseline interferometric observations in this astrophysical context is crucial to characterize the stellar surface dynamics and correct the potential biases. In this Chapter, we present how the granulation pattern is expected for different kind of stellar types ranging from main sequence to extremely evolved stars of different masses and how interferometric techniques help to study their photospheric dynamics.

Radiative hydrodynamics simulations of red supergiant stars

Context. The red supergiant (RSG) Betelgeuse is an irregular variable star. Convection may play an important role in understanding this variability. Interferometric observations can be interpreted using sophisticated simulations of stellar convection. Aims. We compare the visibility curves and closure phases obtained from our 3D simulation of RSG convection with CO5BOLD to various interferometric observations of Betelgeuse from the optical to the H band in order to characterize and measure the convection pattern on this star. Methods. We use 3D radiative-hydrodynamics (RHD) simulation to compute intensity maps in different filters and we thus derive interferometric observables using the post-processing radiative transfer code OPTIM3D. The synthetic visibility curves and closure phases are compared to observations. Results. We provide a robust detection of the granulation pattern on the surface of Betelgeuse in the optical and in the H band based on excellent fits to the observed visibility points and closure phases. Moreover, we determine that the Betelgeuse surface in the H band is covered by small to medium scale (5-15 mas) convection-related surface structures and a large (30 mas) convective cell. In this spectral region, H2O molecules are the main absorbers and contribute to the small structures and to the position of the first null of the visibility curve (i.e. the apparent stellar radius).

Atmospheric Dynamics of Red Supergiant Stars and Interferometry

We developed a 3D pure LTE radiative transfer code to derive observables expected for RSGs, with emphasis on small scale structures, from radiative-hydrodynamic (RHD) simulations of red supergiant stars (RSGs) carried out with CO5BOLD (Freytag et al. 2002). We show that the convection-related surface structures are observable with today's interferometers. Moreover, the RHD simulations are a great improvement over parametric models for the interpretation of interferometric observations.