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.
Abstract
We obtained 16 VLT/X-shooter observations of GX 339−4 in quiescence during the period 2016 May–September and detected absorption lines from the donor star in its NIR spectrum. This allows us ...to measure the radial velocity curve and projected rotational velocity of the donor for the first time. We confirm the 1.76 day orbital period and we find that
K
2
= 219 ± 3 km s
−1
,
γ
= 26 ± 2 km s
−1
, and
v
sin
i
= 64 ± 8 km s
−1
. From these values we compute a mass function
f
(
M
) = 1.91 ± 0.08
M
⊙
, a factor ∼3 lower than previously reported, and a mass ratio
q
= 0.18 ± 0.05. We confirm the donor is a K-type star and estimate that it contributes
∼
4
%
–
50
%
of the light in the
J
- and
H
-bands. We constrain the binary inclination to 37°
<
i
<
78° and the black hole (BH) mass to
2.3
M
⊙
<
M
BH
<
9.5
M
⊙
. GX 339−4 may therefore be the first BH to fall in the “mass-gap” of 2–5
M
⊙
.
Context. Spectroscopy of exoplanet atmospheres at high-resolving powers is rapidly gaining popularity for measuring the presence of atomic and molecular species. While this technique is particularly ...robust against contaminant absorption in the Earth’s atmosphere, the non-stationary stellar spectrum, in the form of either Doppler shift or distortion of the line profile during planetary transits, creates a non-negligible source of noise that can alter or even prevent detection. Aims. Our aim was to use state-of-the art three-dimensional stellar simulations to directly remove the signature of the star from observations prior to cross correlation with templates for the planet’s atmosphere, which are commonly used to extract the faint exoplanet signal from noisy data. Methods. We computed synthetic spectra from 3D simulations of stellar convection resolved both spatially and temporally, and we coupled them with an analytical model reproducing the correct geometry of a transiting exoplanet. We applied the method to the early K-dwarf, HD 189733, and re-analyzed transmission and emission spectroscopy of its hosted exoplanet. In addition, we also analyzed emission spectroscopy of the non transiting exoplanet 51 Pegasi b, orbiting a solar-type star. Results. We find a significant improvement in planet detectability when removing the stellar spectrum with our method. In all cases, we show that the method is superior to a simple parametrisation of the stellar line profile or to the use of 1D stellar models. We show that this is due to the intrinsic treatment of convection in 3D simulations, which allows us to correctly reproduce asymmetric and blue-shifted spectral lines, and intrinsically model center-to-limb variation and Rossiter-McLaughlin effect potentially altering the interpretation of exoplanet transmission spectra. In the case of 51 Pegasi b, we succeed in confirming a previous tentative detection of the planet’s K-band spectrum due to the improved suppression of stellar residuals. Conclusions. Future high-resolution observations will benefit from the synergy with stellar spectroscopy and can be used to test the correct modeling of physical processes in stellar atmospheres. We highlight key improvements in modeling techniques and knowledge of opacity sources to extend this work to shorter wavelengths and later-type stars.
Context.
Despite being the best studied red supergiant star in our Galaxy, the physics behind the photometric variability and mass loss of Betelgeuse is poorly understood. Moreover, recently the star ...has experienced an unusual fading with its visual magnitude reaching a historical minimum. The nature of this event was investigated by several studies where mechanisms, such as episodic mass loss and the presence of dark spots in the photosphere, were invoked.
Aims.
We aim to relate the atmospheric dynamics of Betelgeuse to its photometric variability, with the main focus on the dimming event.
Methods.
We used the tomographic method which allowed us to probe different depths in the stellar atmosphere and to recover the corresponding disk-averaged velocity field. The method was applied to a series of high-resolution HERMES observations of Betelgeuse. Variations in the velocity field were then compared with photometric and spectroscopic variations.
Results.
The tomographic method reveals that the succession of two shocks along our line-of-sight (in February 2018 and January 2019), the second one amplifying the effect of the first one, combined with underlying convection and/or outward motion present at this phase of the 400 d pulsation cycle, produced a rapid expansion of a portion of the atmosphere of Betelgeuse and an outflow between October 2019 and February 2020. This resulted in a sudden increase in molecular opacity in the cooler upper atmosphere of Betelgeuse and, thus, in the observed unusual decrease of the star’s brightness.
Aims. We compute the emergent stellar spectra from the UV to far infrared for different viewing angles using realistic 3D model atmospheres for a large range in stellar parameters to predict the ...stellar limb darkening. Methods. We have computed full 3D LTE synthetic spectra based on 3D radiative hydrodynamic atmosphere models from the STAGGER-grid in the ranges: T sub(eff) from 4000 to 7000 K, log g from 1.5 to 5.0, and Fe/H, from -4.0 to +0.5. From the resulting intensities, we derived coefficients for the standard limb darkening laws considering a number of often-used photometric filters. Furthermore, we calculated theoretical transit light curves, in order to quantify the differences between predictions by the widely used 1D model atmosphere and our 3D models. Results. The 3D models are often found to predict steeper darkening towards the limb compared to the 1D models, mainly due to the temperature stratifications and temperature gradients being different in the 3D models compared to those predicted with 1D models based on the mixing length theory description of convective energy transport. The resulting differences in the transit light curves are rather small; however, these can be significant for high-precision observations of extrasolar transits, and are able to lower the residuals from the fits with 1D limb darkening profiles. Conclusions. We advocate the use of the new limb darkening coefficients provided for the standard four-parameter non-linear power law, which can fit the limb darkening more accurately than other choices.
Context. It has been proposed that convection in red supergiant (RSG) stars produces large-scale granules causing observable surface inhomogeneities. This convection is also extremely vigorous and is ...suspected to be one of the main causes of mass-loss in RSGs. It should thus be understood in detail. Evidence has accumulated of asymmetries in the photospheres of RSGs, but detailed studies of granulation are still lacking. Interferometric observations provide an innovative way of addressing this question, but they are still often interpreted using smooth symmetrical limb-darkened intensity distributions, or simple, spotted, ad hoc models. Aims. We explore the impact of the granulation on visibility curves and closure phases using the radiative transfer code OPTIM3D. We simultaneously assess how 3D simulations of convection in RSG with CO5BOLD can be tested by comparing with these observations. Methods. We use 3D radiative hydrodynamical (RHD) simulations of convection to compute intensity maps at various wavelengths and time, from which we derive interferometric visibility amplitudes and phases. We study their behaviour with time, position angle, and wavelength, and compare them to observations of the RSG α Ori. Results. We provide average limb-darkening coefficients for RSGs. We describe the prospects for the detection and characterization of granulation (i.e., contrast, size) on RSGs. We demonstrate that our RHD simulations provide an excellent fit to existing interferometric observations of α Ori, in contrast to limb darkened disks. This confirms the existence of large convective cells on the surface of Betelgeuse.
Context. Red supergiants are massive evolved stars that contribute extensively to the chemical enrichment of our Galaxy. It has been shown that convection in those stars produces large granules that ...cause surface inhomogeneities and shock waves in the photosphere. The understanding of their dynamics is crucial for unveiling the unknown mass-loss mechanism, their chemical composition, and their stellar parameters. Aims: We present a new generation of red supergiant simulations with a more sophisticated opacity treatment performed with 3D radiative-hydrodynamics code CO5BOLD. Methods: In the code the coupled equations of compressible hydrodynamics and non-local radiation transport are solved in the presence of a spherical potential. The stellar core is replaced by a special spherical inner boundary condition, where the gravitational potential is smoothed and the energy production by fusion is mimicked by a simply producing heat corresponding to the stellar luminosity. All outer boundaries are transmitting for matter and light. The post-processing radiative transfer code OPTIM3D is used to extract spectroscopic and interferometric observables. Results: We show that if one relaxes the assumption of frequency-independent opacities, this leads to a steeper mean thermal gradient in the optical thin region that strongly affects the atomic strengths and the spectral energy distribution. Moreover, the weaker temperature fluctuations reduce the incertitude on the radius determination with interferometry. We show that 1D models of red supergiants must include a turbulent velocity that is calibrated on 3D simulations to obtain the effective surface gravity that mimics the effect of turbulent pressure on the stellar atmosphere. We provide an empirical calibration of the ad hoc micro- and macroturbulence parameters for 1D models using the 3D simulations: we find that there is no clear distinction between the different macroturbulent profiles needed in 1D models to fit 3D synthetic lines.
We present the STAGGER-grid, a comprehensive grid of time-dependent, threedimensional (3D), hydrodynamic model atmospheres for late-type stars with realistic treatment of radiative transfer, covering ...a wide range in stellar parameters. This grid of 3D models is intended for various applications besides studies of stellar convection and atmospheres per se, including stellar parameter determination, stellar spectroscopy and abundance analysis, asteroseismology, calibration of stellar evolution models, interferometry, and extrasolar planet search. In this introductory paper, we describe the methods we applied for the computation of the grid and discuss the general properties of the 3D models as well as of their temporal and spatial averages. For the first time, we are able to quantify systematically over a broad range of stellar parameters the uncertainties of 1D models arising from the simplified treatment of physics, in particular convective energy transport. In agreement with previous findings, we find that the differences can be rather significant, especially for metal-poor stars.