skip to main content

Title: 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 O PTIM 3D 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.  more » « less
Award ID(s):
1715788 1636624 1506540
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; « less
Date Published:
Journal Name:
Astronomy & Astrophysics
Page Range / eLocation ID:
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. Stars evolving along the asymptotic giant branch (AGB) can become carbon rich in the final part of their evolution. The detailed description of their spectra has led to the definition of several spectral types: N, SC, J, and R. To date, differences among them have been partially established only on the basis of their chemical properties. Aims. An accurate determination of the luminosity function (LF) and kinematics together with their chemical properties is extremely important for testing the reliability of theoretical models and establishing on a solid basis the stellar population membership of the different carbon star types. Methods. Using Gaia Data Release 2 ( Gaia DR2) astrometry, we determine the LF and kinematic properties of a sample of 210 carbon stars with different spectral types in the solar neighbourhood with measured parallaxes better than 20%. Their spatial distribution and velocity components are also derived. Furthermore, the use of the infrared Wesenheit function allows us to identify the different spectral types in a Gaia -2MASS diagram. Results. We find that the combined LF of N- and SC-type stars are consistent with a Gaussian distribution peaking at M bol  ∼ −5.2 mag. The resulting LF, however, shows two tails at lower and higher luminosities more extended than those previously found, indicating that AGB carbon stars with solar metallicity may reach M bol  ∼ −6.0 mag. This contrasts with the narrower LF derived in Galactic carbon Miras from previous studies. We find that J-type stars are about half a magnitude fainter on average than N- and SC-type stars, while R-hot stars are half a magnitude brighter than previously found, although fainter in any case by several magnitudes than other carbon types. Part of these differences are due to systematically lower parallaxes measured by Gaia DR2 with respect to H IPPARCOS values, in particular for sources with parallax ϖ < 1 mas. The Galactic spatial distribution and velocity components of the N-, SC-, and J-type stars are very similar, while about 30% of the R-hot stars in the sample are located at distances greater than ∼500 pc from the Galactic plane, and show a significant drift with respect to the local standard of rest. Conclusions. The LF derived for N- and SC-type in the solar neighbourhood fully agrees with the expected luminosity of stars of 1.5−3 M ⊙ on the AGB. On a theoretical basis, the existence of an extended low-luminosity tail would require a contribution of extrinsic low-mass carbon stars, while the high-luminosity tail would imply that stars with mass values up to ∼5 M ⊙ may become carbon stars on the AGB. J-type stars differ significantly not only in their chemical composition with respect to the N- and SC-types, but also in their LF, which reinforces the idea that these carbon stars belong to a different type whose origin is still unknown. The derived luminosities of R-hot stars means that it is unlikely that these stars are in the red-clump, as previously claimed. On the other hand, the derived spatial distribution and kinematic properties, together with their metallicity values, indicate that most of the N-, SC-, and J-type stars belong to the thin disc population, while a significant fraction of R-hot stars show characteristics compatible with the thick disc. 
    more » « less

    Westerlund 1 (Wd1) is potentially the largest star cluster in the Galaxy. That designation critically depends upon the distance to the cluster, yet the cluster is highly obscured, making luminosity-based distance estimates difficult. Using Gaia Data Release 2 (DR2) parallaxes and Bayesian inference, we infer a parallax of $0.35^{+0.07}_{-0.06}$ mas corresponding to a distance of $2.6^{+0.6}_{-0.4}$ kpc. To leverage the combined statistics of all stars in the direction of Wd1, we derive the Bayesian model for a cluster of stars hidden among Galactic field stars; this model includes the parallax zero-point. Previous estimates for the distance to Wd1 ranged from 1.0 to 5.5 kpc, although values around 5 kpc have usually been adopted. The Gaia DR2 parallaxes reduce the uncertainty from a factor of 3 to 18 per cent and rules out the most often quoted value of 5 kpc with 99 per cent confidence. This new distance allows for more accurate mass and age determinations for the stars in Wd1. For example, the previously inferred initial mass at the main-sequence turn-off was around 40 M⊙; the new Gaia DR2 distance shifts this down to about 22 M⊙. This has important implications for our understanding of the late stages of stellar evolution, including the initial mass of the magnetar and the LBV in Wd1. Similarly, the new distance suggests that the total cluster mass is about four times lower than previously calculated.

    more » « less
  3. Abstract Increasing main-sequence stellar luminosity with stellar mass leads to the eventual dominance of radiation pressure in stellar-envelope hydrostatic balance. As the luminosity approaches the Eddington limit, additional instabilities (beyond conventional convection) can occur. These instabilities readily manifest in the outer envelopes of OB stars, where the opacity increase associated with iron yields density and gas-pressure inversions in 1D models. Additionally, recent photometric surveys (e.g., TESS) have detected excess broadband low-frequency variability in power spectra of OB star lightcurves, called stochastic low-frequency variability (SLFV). This motivates our novel 3D Athena++ radiation hydrodynamical (RHD) simulations of two 35 M ⊙ star envelopes (the outer ≈15% of the stellar radial extent), one on the zero-age main sequence and the other in the middle of the main sequence. Both models exhibit turbulent motion far above and below the conventional iron-opacity peak convection zone (FeCZ), obliterating any “quiet” part of the near-surface region and leading to velocities at the photosphere of 10–100 km s −1 , directly agreeing with spectroscopic data. Surface turbulence also produces SLFV in model lightcurves with amplitudes and power-law slopes that are strikingly similar to those of observed stars. The characteristic frequencies associated with SLFV in our models are comparable to the thermal time in the FeCZ (≈3–7 day −1 ). These ab initio simulations are directly validated by observations and, though more models are needed, we remain optimistic that 3D RHD models of main-sequence O-star envelopes exhibit SLFV originating from the FeCZ. 
    more » « less
  4. Context . The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R ⊕ and 3.0 M ⊕ . It is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets. Aims . To prepare for future studies, we aim to thoroughly characterise the planetary system with new accurate and precise data collected with state-of-the-art photometers from space and spectrometers and interferometers from the ground. Methods . We collected light curves of seven new transits observed with the CHEOPS space mission and new radial velocities obtained with MAROON-X at the 8.1 m Gemini North telescope and CARMENES at the 3.5 m Calar Alto telescope, together with previously published spectroscopic and photometric data from the two spectrographs and TESS. We also performed near-infrared interferometric observations with the CHARA Array and new photometric monitoring with a suite of smaller telescopes (AstroLAB, LCOGT, OSN, TJO). This extraordinary and rich data set was the input for our comprehensive analysis. Results . From interferometry, we measure a limb-darkened disc angular size of the star Gl 486 at θ LDD = 0.390 ± 0.018 mas. Together with a corrected Gaia EDR3 parallax, we obtain a stellar radius R * = 0.339 ± 0.015 R ⊕ . We also measure a stellar rotation period at P rot = 49.9 ± 5.5 days, an upper limit to its XUV (5-920 A) flux informed by new Hubble /STIS data, and, for the first time, a variety of element abundances (Fe, Mg, Si, V, Sr, Zr, Rb) and C/O ratio. Moreover, we imposed restrictive constraints on the presence of additional components, either stellar or sub-stellar, in the system. With the input stellar parameters and the radial-velocity and transit data, we determine the radius and mass of the planet Gl 486 b at R p = 1.343 −0.062 +0.063 R ⊕ and M p = 3.00 −0.12 +0.13 M ⊕ , with relative uncertainties of the planet radius and mass of 4.7% and 4.2%, respectively. From the planet parameters and the stellar element abundances, we infer the most probable models of planet internal structure and composition, which are consistent with a relatively small metallic core with respect to the Earth, a deep silicate mantle, and a thin volatile upper layer. With all these ingredients, we outline prospects for Gl 486 b atmospheric studies, especially with forthcoming James Webb Space Telescope ( Webb ) observations. 
    more » « less
  5. Abstract Observations indicate that turbulent motions are present on most massive star surfaces. Starting from the observed phenomena of spectral lines with widths that are much larger than their thermal broadening (e.g., micro- and macroturbulence), and considering the detection of stochastic low-frequency variability (SLFV) in the Transiting Exoplanet Survey Satellite photometry, these stars clearly have large-scale turbulent motions on their surfaces. The cause of this turbulence is debated, with near-surface convection zones, core internal gravity waves, and wind variability being proposed. Our 3D gray radiation hydrodynamic (RHD) models previously characterized the convective dynamics of the surfaces, driven by near-surface convection zones, and provided reasonable matches to the observed SLFV of the most luminous massive stars. We now explore the complex emitting surfaces of these 3D RHD models, which strongly violate the 1D assumption of a plane-parallel atmosphere. By post-processing the gray RHD models with the Monte Carlo radiation transport code Sedona , we synthesize stellar spectra and extract information from the broadening of individual photospheric lines. The use of Sedona enables the calculation of the viewing angle and temporal dependence of spectral absorption line profiles. By combining uncorrelated temporal snapshots together, we compare the turbulent broadening from the 3D RHD models to the thermal broadening of the extended emitting region, showing that our synthesized spectral lines closely resemble the observed macroturbulent broadening from similarly luminous stars. More generally, the new techniques that we have developed will allow for systematic studies of the origins of turbulent velocity broadening from any future 3D simulations. 
    more » « less