The photospheric radius is one of the fundamental parameters governing the radiative equilibrium of a star. We report new observations of the nearest solar-type stars alpha Centauri A (G2V) and B ...(K1V) with the VLTI/PIONIER optical interferometer. The combination of four configurations of the VLTI enable us to measure simultaneously the limb darkened angular diameter theta sub(LD) and the limb darkening parameters of the two solar-type stars in the near-infrared H band (lambda= 1.65 mu m). We obtain photospheric angular diameters of theta sub(LD)(A) = 8.502 + or - 0.038mas (0.43%) and theta sub(LD)(B) = 5.999 + or - 0.025mas (0.42%), through the adjustment of a power law limb darkening model. We find H band power law exponents of alpha (A) = 0.1404 + or - 0.0050 (3.6%) and alpha (B) = 0.1545 + or - 0.0044 (2.8%), which closely bracket the observed solar value ( alpha sub(middot in circle)= 0.15027). Combined with the parallax pi= 747.17 + or - 0.61mas previously determined, we derive linear radii of R sub(A)= 1.2234 + or - 0.0053 R sub(middot in circle)(0.43%) and R sub(B)= 0.8632 + or - 0.0037 R sub(middot in circle)(0.43%). The power law exponents that we derive for the two stars indicate a significantly weaker limb darkening than predicted by both 1D and 3D stellar atmosphere models. As this discrepancy is also observed on the near-infrared limb darkening profile of the Sun, an improvement of the calibration of stellar atmosphere models is clearly needed. The reported PIONIER visibility measurements of alpha Cen A and B provide a robust basis to validate the future evolutions of these models.
Context. Clouds are ubiquitous in exoplanet atmospheres and they represent a challenge for the model interpretation of their spectra. When generating a large number of model spectra, complex cloud ...models often prove too costly numerically, whereas more efficient models may be overly simplified. Aims. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. Methods. We used our radiative transfer code petitRADTRANS for generating the spectra, which we coupled to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. Results. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE, and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that C/O = 0.60 −0.08 +0.07 . Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: either cloudy or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, without constraining C/O. Conclusions. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO 2 or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.
Context. Theory surrounding the origin of the dust-laden winds from evolved stars remains mired in controversy. Characterizing the formation loci and the dust distribution within approximately the ...first stellar radius above the surface is crucial for understanding the physics that underlie the mass-loss phenomenon. Aims. By exploiting interferometric polarimetry, we derive the fundamental parameters that govern the dust structure at the wind base of a red supergiant. Methods. We present near-infrared aperture-masking observations of Betelgeuse in polarimetric mode obtained with the NACO/SAMPol instrument. We used both parametric models and radiative transfer simulations to predict polarimetric differential visibility data and compared them to SPHERE/ZIMPOL measurements. Results. Using a thin dust shell model, we report the discovery of a dust halo that is located at only 0.5 R⋆ above the photosphere (i.e. an inner radius of the dust halo of 1.5 R⋆). By fitting the data under the assumption of Mie scattering, we estimate the grain size and density for various dust species. By extrapolating to the visible wavelengths using radiative transfer simulations, we compare our model with SPHERE/ZIMPOL data and find that models based on dust mixtures that are dominated by forsterite are most favored. Such a close dusty atmosphere has profound implications for the dust formation mechanisms around red supergiants.
Abstract
Emission lines from Rydberg transitions are detected for the first time from a region close to the surface of Betelgeuse. The H30
α
line is observed at 231.905 GHz, with an FWHM ∼42 km s
−1
...and extended wings. A second line at 232.025 GHz (FWHM ∼21 km s
−1
), is modeled as a combination of Rydberg transitions of abundant low first ionization potential metals. Both H30
α
and the Rydberg combined line X30
α
are fitted by Voigt profiles, and collisional broadening with electrons may be partly responsible for the Lorentzian contribution, indicating electron densities of a few 10
8
cm
−3
. X30
α
is located in a relatively smooth ring at a projected radius of 0.9× the optical photospheric radius
R
⋆
, whereas H30
α
is more clumpy, reaching a peak at ∼1.4
R
⋆
. We use a semiempirical thermodynamic atmospheric model of Betelgeuse to compute the 232 GHz (1.29 mm) continuum and line profiles making simple assumptions. Photoionized abundant metals dominate the electron density, and the predicted surface of continuum optical depth unity at 232 GHz occurs at ∼1.3
R
⋆
, in good agreement with observations. Assuming a Saha–Boltzmann distribution for the level populations of Mg, Si, and Fe, the model predicts that the X30
α
emission arises in a region of radially increasing temperature and turbulence. Inclusion of ionized C and non-LTE effects could modify the integrated fluxes and location of emission. These simulations confirm the identity of the Rydberg transition lines observed toward Betelgeuse and reveal that such diagnostics can improve future atmospheric models.
Rapid rotation is a fundamental characteristic of classical Be stars and a crucial property allowing for the formation of their circumstellar disks. Past evolution in a mass and angular momentum ...transferring binary system offers a plausible solution to how Be stars attained their fast rotation. Although the subdwarf remnants of mass donors in such systems should exist in abundance, only a few have been confirmed due to tight observational constraints. An indirect method of detecting otherwise hidden companions is offered by their effect on the outer parts of Be star disks, which are expected to be disrupted or truncated. In the context of the infrared and radio continuum excess radiation originating in the disk, the disk truncation can be revealed by a turndown in the spectral energy distribution due to reduced radio flux levels. In this work, we search for signs of spectral turndown in a sample of 57 classical Be stars with radio data, which include new data for 23 stars and the longest-wavelength detections so far (λ 10 cm) for two stars. We confidently detect the turndown for all 26 stars with sufficient data coverage (20 of which are not known to have close binary companions). For the remaining 31 stars, the data are inconclusive as to whether the turndown is present or not. The analysis suggests that many if not all Be stars have close companions influencing their outer disks. If confirmed to be subdwarf companions, the mass transfer spin-up scenario might explain the existence of the vast majority of classical Be stars.
Context. The precise determinations of stellar mass at ≲1% provide important constraints on stellar evolution models. Accurate parallax measurements can also serve as independent benchmarks for the ...next Gaia data release. Aims. We measured the masses and distance of binary systems with a precision level better than 1% using a fully geometrical and empirical method. Methods. We obtained the first interferometric observations for the eclipsing systems AI Phe, AL Dor, KW Hya, NN Del, ψ Cen and V4090 Sgr with the VLTI/PIONIER combiner, which we combined with radial velocity measurements to derive their three-dimensional orbit, masses, and distance. Results. We determined very precise stellar masses for all systems, ranging in precision from 0.04% to 3.3%. We combined these measurements with the stellar effective temperature and linear radius to fit stellar isochrones models and determined the age of the systems. We also derived the distance to the systems with a precision level of 0.4%. Conclusions. The comparison of theoretical models with stellar parameters shows that stellar models are still deficient in simultaneously fitting the stellar parameters (Teff, R and M) with this level of precision on individual masses. This stresses the importance of precisely measuring the stellar parameters to better calibrate stellar evolution models. The precision of our model-independent orbital parallaxes varies from 24 μas as to 70 μas and the parallaxes provide a unique opportunity to verify whether the future Gaia measurements have systematic errors.
The availability of a number of new interferometric measurements of Main Sequence and subgiant stars makes it possible to calibrate the surface brightness relations of these stars using exclusively ...direct angular diameter measurements. These empirical laws make it possible to predict the limb darkened angular diameters theta sub(LD) of dwarfs and subgiants using their dereddened Johnson magnitudes, or their effective temperature. The smallest intrinsic dispersions of sigma less than or equal to 1% in theta sub(LD) are obtained for the relations based on the K and L magnitudes, for instance log theta sub(LD) = 0.0502 (B - L) + 0.5133 - 0.2L or log theta sub(LD) = 0.0755 (V - K) + 0.5170 - 0.2K. Our calibrations are valid between the spectral types A0 and M2 for dwarf stars (with a possible extension to later types when using the effective temperature), and between A0 and K0 for subgiants. Such relations are particularly useful for estimating the angular sizes of calibrators for long-baseline interferometry from readily available broadband photometry.
Context. Betelgeuse is an M-type supergiant that presents a circularly polarized (Stokes V) signal in its line profiles, interpreted in terms of a surface magnetic field. Aims. The weak circular ...polarization signal has been monitored over 7.5 years in order to follow its evolution on different timescales, and eventually to determine its physical origin. Linear polarization measurements have also been obtained regularly in the last few years. Methods. We used both the ESPaDOnS and Narval spectropolarimeters to obtain high signal-to-noise ratio spectra, which were processed by means of the least-squares deconvolution method. In order to ensure the reality of the very weak circular polarization, special care has been taken to limit instrumental effects. In addition, several tests were performed on the Stokes V signal to establish its stellar and Zeeman origin. Results. We confirm the magnetic nature of the circular polarization, pointing to a surface magnetic field of the order of 1 G. The Stokes V profiles present variations over different timescales, the most prominent one being close to the long secondary period (LSP; around 2000 d for Betelgeuse) often invoked in red evolved stars. This long period is also dominant for all the other Stokes parameters. The circular polarization is tentatively modeled by means of magnetic field concentrations mimicking spots, showing in particular that the velocity associated with each “spot” also follows the long timescale, and that this signal is nearly always slightly redshifted. Conclusions. From the coupled variations of both linear and circular polarization signatures in amplitude, velocity and timescale, we favour giant convection cells as the main engine at the origin of polarization signatures and variations in all the Stokes parameters. This strengthens support for the hypothesis that large convective cells are at the origin of the LSP.
Context.The universality of the Cepheid period-luminosity (PL) relations has been under discussion since metallicity effects were assumed to play a role in the value of the intercept and, more ...recently, of the slope of these relations. Aims.The goal of the present study is to calibrate the Galactic PL relations in various photometric bands (from B to K) and to compare the results to the well-established PL relations in the LMC. Methods.We use a set of 59 calibrating stars, the distances of which are measured using five different distance indicators: Hubble Space Telescope and revised Hipparcos parallaxes, infrared surface brightness and interferometric Baade-Wesselink parallaxes, and classical Zero-Age-Main-Sequence-fitting parallaxes for Cepheids belonging to open clusters or OB stars associations. A detailed discussion of absorption corrections and projection factor to be used is given. Results.We find no significant difference in the slopes of the PL relations between LMC and our Galaxy. Conclusions.We conclude that the Cepheid PL relations have universal slopes in all photometric bands, not depending on the galaxy under study (at least for LMC and Milky Way). The possible zero-point variation with metal content is not discussed in the present work, but an upper limit of 18.50 for the LMC distance modulus can be deduced from our data.
We provide an update for the astrometric solution for the Type II supernova progenitor Betelgeuse using the revised Hipparcos Intermediate Astrometric Data (HIAD) of van Leeuwen, combined with ...existing VLA and new e-MERLIN and ALMA positions. The 2007 Hipparcos refined abscissa measurements required the addition of so-called Cosmic Noise of 2.4 mas to find an acceptable 5-parameter stochastic solution. We find that a measure of radio Cosmic Noise should also be included for the radio positions because surface inhomogeneities exist at a level significant enough to introduce additional intensity centroid uncertainty. Combining the 2007 HIAD with the proper motions based solely on the radio positions leads to a parallax of mas ( pc), smaller than the Hipparcos 2007 value of 6.56 0.83 mas ( pc). Furthermore, combining the VLA and new e-MERLIN and ALMA radio positions with the 2007 HIAD, and including radio Cosmic Noise of 2.4 mas, leads to a nominal parallax solution of 4.51 0.80 mas ( pc), which, while only different from the 2008 solution of Harper et al., is 2.6 different from the solution of van Leeuwen. An accurate and precise parallax for Betelgeuse is always going to be difficult to obtain because it is small compared to the stellar angular diameter ( mas). We outline an observing strategy utilizing future mm and sub-mm high-spatial resolution interferometry that must be used if substantial improvements in the precision and accuracy of the parallax and distance are to be achieved.