ABSTRACT
The processes by which red supergiants lose mass are not fully understood thus-far and their mass-loss rates lack theoretical constraints. The ambient surroundings of the nearby M0.5 Iab ...star Antares offer an ideal environment to obtain detailed empirical information on the outflow properties at its onset, and hence indirectly, on the mode(s) of mass-loss. We present and analyse optical VLT/SPHERE/ZIMPOL polarimetric imaging with angular resolution down to 23 milliarcsec, sufficient to spatially resolve both the stellar disc and its direct surroundings. We detect a conspicuous feature in polarized intensity that we identify as a clump containing dust, which we characterize through 3D radiative transfer modelling. The clump is positioned behind the plane of the sky, therefore has been released from the backside of the star, and its inner edge is only 0.3 stellar radii above the surface. The current dust mass in the clump is $1.3^{+0.2}_{-1.0} \times 10^{-8} \, \mathrm{M}_{\odot }$, though its proximity to the star implies that dust nucleation is probably still ongoing. The ejection of clumps of gas and dust makes a non-negligible contribution to the total mass lost from the star that could possibly be linked to localized surface activity such as convective motions or non-radial pulsations.
Rotational transitions in vibrationally excited AlO and TiO-two possible precursors of dust-were observed in the 300 GHz range (1 mm wavelength) toward the oxygen-rich asymptotic giant branch (AGB) ...stars R Dor and IK Tau with ALMA, and vibrationally excited AlO was observed toward the red supergiant VY CMa with the SMA. The J = 11 → 10 transition of TiO in the v = 1 and 2 levels and the N = 9 → 8 transition in the v = 2 level of AlO were identified toward R Dor, the J = 11 → 10 line of TiO was identified in the v = 1 level toward IK Tau, and two transitions in the v = 1 and 2 levels of AlO were identified toward VY CMa. The newly derived high vibrational temperature of TiO and AlO in R Dor of 1800 200 K and prior measurements of the angular extent confirm that the majority of the emission is from a region within 2R of the central star. A full radiative transfer analysis of AlO in R Dor yielded a fractional abundance of ∼3% of the solar abundance of Al. From a similar analysis of TiO a fractional abundance of ∼78% of the solar abundance of Ti was found. The observations provide indirect evidence that TiO is present in a rotating disk close to the star. Further observations in the ground and excited vibrational levels are needed to determine whether AlO, TiO, and TiO2 are seeds of the Al2O3 dust in R Dor, and perhaps in the gravitationally bound dust shells in other AGB stars with low mass-loss rates.
We study the hydrodynamical behavior occurring in the turbulent interaction zone of a fast-moving red supergiant star, where the circumstellar and interstellar material collide. In this ...wind-interstellar-medium collision, the familiar bow shock, contact discontinuity, and wind termination shock morphology form, with localized instability development. Our model includes a detailed treatment of dust grains in the stellar wind and takes into account the drag forces between dust and gas. The dust is treated as pressureless gas components binned per grain size, for which we use 10 representative grain size bins. Our simulations allow us to deduce how dust grains of varying sizes become distributed throughout the circumstellar medium. We show that smaller dust grains (radius <0.045 Delta *mm) tend to be strongly bound to the gas and therefore follow the gas density distribution closely, with intricate fine structure due to essentially hydrodynamical instabilities at the wind-related contact discontinuity. Larger grains which are more resistant to drag forces are shown to have their own unique dust distribution, with progressive deviations from the gas morphology. Specifically, small dust grains stay entirely within the zone bound by shocked wind material. The large grains are capable of leaving the shocked wind layer and can penetrate into the shocked or even unshocked interstellar medium. Depending on how the number of dust grains varies with grain size, this should leave a clear imprint in infrared observations of bow shocks of red supergiants and other evolved stars.
Context.
The asymptotic giant branch (AGB) marks the final evolutionary stage of stars with initial masses between ~0.8 and 8
M
⊙
. During this phase, stars undergo copious mass loss, which ...contributes significantly to the enrichment of the interstellar medium. The well-accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the density-enhanced and extended AGB stellar atmospheres. The details of the mass-loss process are not yet well understood, however. For oxygen-rich AGB stars, which are the focus of this study, the dust grains that drive the wind are expected to scatter visible light very efficiently because their sizes are relative large.
Aims.
We study the distribution of dust in the inner wind of oxygen-rich AGB stars to advance our understanding of the wind-driving process.
Methods.
We observed light scattered off dust grains that form around three oxygen-rich AGB stars (W Hya, SW Vir, and R Crt) with mass-loss rates between 10
−7
and 10
−6
M
⊙
yr
−1
using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL with three filters centred at 0.65, 0.75 and 0.82
μ
m. We compared the observed morphologies and the spectral dependence of the scattered light between the three sources and determined the radial profile, per image octant, of the dust density distribution around the closest target, W Hya.
Results.
We find the distribution of dust to be asymmetric for the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be inferred from the ZIMPOL data, but we speculate that it might be the consequence of a circumstellar disc or of the action of strong magnetic fields. The dust grains polarise light more efficiently at 0.65
μ
m for R Crt and SW Vir and at 0.82
μ
m for W Hya. This indicates that at the time of the observations, the grains around SW Vir and R Crt had sizes <0.1
μ
m, while those around W Hya were larger, with sizes ≳0.1
μ
m. The asymmetric distribution of dust around R Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. We fit radiative transfer models to the radial profile of the polarised light observed around W Hya and find a steep dust density profile, with steepness varying considerably with direction. We find the wind-acceleration region of W Hya to extend to at least ~7
R
⋆
. This is in agreement with theoretical predictions of wind acceleration up to ~12
R
⋆
, and highlights that ZIMPOL observations probe the crucial region around AGB stars where dust forms and is accelerated.
Context. High-angular-resolution observations of asymptotic giant branch (AGB) stars often reveal non-spherical morphologies for the gas and dust envelopes. Aims. We aim to make a pilot study to ...quantify the impact of different geometries (spherically symmetric, spiral-shaped, and disc-shaped) of the dust component of AGB envelopes on spectral energy distributions (SEDs), mass estimates, and subsequent mass-loss rate (MLR) estimates. We also estimate the error made on the MLR if the SED is fitted by an inappropriate geometrical model. Methods. We use the three-dimensional Monte-Carlo-based radiative-transfer code RADMC-3D to simulate emission from dusty envelopes with different geometries (but fixed spatial extension). We compare these predictions with each other, and with the SED of the AGB star EP Aqr that we use as a benchmark since its envelope is disc-like and known to harbour spiral arms, as seen in CO. Results. The SEDs involving the most massive envelopes are those for which the different geometries have the largest impact, primarily on the silicate features at 10 and 18 μm. These different shapes originate from large differences in optical depths. Massive spirals and discs appear akin to black bodies. Optically thick edge-on spirals and discs (with dust masses of 10−4 and 10−5 M⊙) exhibit black-body SEDs that appear cooler than those from face-on structures and spheres of the same mass, while optically thick face-on distributions appear as warmer emission. We find that our more realistic models, combined spherical and spiral distributions, are 0.1 to 0.5 times less massive than spheres with similar SEDs. More extreme, less realistic scenarios give that spirals and discs are 0.01 to 0.05 times less massive than corresponding spheres. This means that adopting the wrong geometry for an AGB circumstellar envelope may result in a MLR that is incorrect by as much as one to two orders of magnitude when derived from SED fitting.
Aims. Cool, evolved stars have copious, enriched winds. Observations have so far not fully constrained models for the shaping and acceleration of these winds. We need to understand the dynamics ...better, from the pulsating stellar surface to ~10 stellar radii, where radiation pressure on dust is fully effective. Asymmetric nebulae around some red supergiants imply the action of additional forces. Methods. We retrieved ALMA Science Verification data providing images of sub-mm line and continuum emission from VY CMa. This enables us to locate water masers with milli-arcsec accuracy and to resolve the dusty continuum. Results. The 658, 321, and 325 GHz masers lie in irregular, thick shells at increasing distances from the centre of expansion. For the first time this is confirmed as the stellar position, coinciding with a compact peak offset to the NW of the brightest continuum emission. The maser shells overlap but avoid each other on scales of up to 10 au. Their distribution is broadly consistent with excitation models but the conditions and kinematics are complicated by wind collisions, clumping, and asymmetries.
Context. Equatorial density enhancements (EDEs) are a very common astronomical phenomenon. Studies of the circumstellar environments (CSE) of young stellar objects and of evolved stars have shown ...that these objects often possess these features. These are believed to originate from different mechanisms, ranging from binary interactions to the gravitational collapse of interstellar material. Quantifying the effect of the presence of this type of EDE on the observables is essential for a correct interpretation of high-resolution data. Aims. We seek to investigate the manifestation in the observables of a circumstellar EDE, to assess which properties can be constrained, and to provide an intuitive bedrock on which to compare and interpret upcoming high-resolution data (e.g. ALMA data) using 3D models. Methods. We develop a simplified analytical parametrised description of a 3D EDE, with possible substructure such as warps, gaps, and spiral instabilities. In addition, different velocity fields (Keplerian, radial, super-Keplerian, sub-Keplerian and rigid rotation) are considered. The effect of a bipolar outflow is also investigated. The geometrical models are fed into the 3D radiative transfer code LIME, that produces 3D intensity maps throughout velocity space. We investigate the spectral signature of the J= 3-2 up to J= 7-6 rotational transitions of CO in the models, as well as the spatial aspect of this emission by means of channel maps, wide-slit position-velocity (PV) diagrams, stereograms, and spectral lines. Additionally, we discuss methods of constraining the geometry of the EDE, the inclination, the mass-contrast between the EDE and the bipolar outflow, and the global velocity field. Finally, we simulated ALMA observations to explore the effects of interferometric noise and artefacts on the emission signatures. Results. The effects of the different velocity fields are most evident in the PV diagrams. These diagrams also enable us to constrain the EDE height and inclination. A level of degeneracy may occur in the shapes of individual PV diagrams for different global velocity fields. The orthogonal PV diagrams may completely eliminate this ambiguity. Information on the EDE substructure is evident in the channel maps, but cannot be recovered from the PV diagrams, nor from the spectral lines. However, stereograms enable the detection of warping. For most inclinations the spectral lines are relatively broad, making it difficult to distinguish from an eventual superposed bipolar outflow component. Only under low inclination angles can one distinguish between these structures. Simulations of synthetic ALMA observations show how emission is affected when the largest angular scale of an antenna configuration is exceeded. For a rotating EDE, the emission around zero velocity will first fade because of destructive interference.
Context. Red supergiant (RSG) stars exhibit significant mass loss by means of a slow, dense wind. They are often considered to be the more massive counterparts of Asymptotic Giant Branch (AGB) stars. ...While AGB mass loss is related to their strong pulsations, the RSG are often only weakly variable. This raises the question of whether their wind-driving mechanism and the dust composition of the wind are the same. Aims. We study the conditions at the base of the wind by determining the dust composition of a sample of RSG. The dust composition is assumed to be sensitive to the density, temperature, and acceleration at the base of the wind. We compare the derived dust composition with the composition measured in AGB star winds. Methods. We compile a sample of 27 RSG infrared spectra (ISO-SWS) and supplement these with photometric measurements to derive the full spectral energy distribution (SED). These data are modelled using a dust radiative-transfer code, taking into account the optical properties of the relevant candidate materials to search for correlations between mass-loss rate, density at the inner edge of the dust shell, and stellar parameters. Results. We find strong correlations between the dust composition, mass-loss rate, and the stellar luminosity, roughly in agreement with the theoretical dust condensation sequence. We identify the need for a continuous (near-)IR dust opacity and tentatively propose amorphous carbon, and we note significant differences with AGB star winds in terms of the presence of PAHs, absence of “the” 13 μm band, and a lack of strong water bands. Conclusions. Dust condensation in RSG is found to experience a freeze-out process that is similar to that in AGB stars. Together with the positive effect of the stellar luminosity on the mass-loss rate, this suggests that radiation pressure on dust grains is an important ingredient in the driving mechanism. Still, differences with AGB stars are manifold and thus the winds of RSG should be studied individually in further detail.
Context. The evolution of intermediate and low-mass stars on the asymptotic giant branch is dominated by their strong dust-driven winds. More massive stars evolve into red supergiants with a similar ...envelope structure and strong wind. These stellar winds are a prime source for the chemical enrichment of the interstellar medium. Aims: We aim to (1) set up simple and general analytical expressions to estimate mass-loss rates of evolved stars, and (2) from those calculate estimates for the mass-loss rates of the asymptotic giant branch, red supergiant, and yellow hypergiant stars in our galactic sample. Methods: The rotationally excited lines of carbon monoxide (CO) are a classic and very robust diagnostic in the study of circumstellar envelopes. When sampling different layers of the circumstellar envelope, observations of these molecular lines lead to detailed profiles of kinetic temperature, expansion velocity, and density. A state-of-the-art, nonlocal thermal equilibrium, and co-moving frame radiative transfer code that predicts CO line intensities in the circumstellar envelopes of late-type stars is used in deriving relations between stellar and molecular-line parameters, on the one hand, and mass-loss rate, on the other. These expressions are applied to our extensive CO data set to estimate the mass-loss rates of 47 sample stars. Results: We present analytical expressions for estimating the mass-loss rates of evolved stellar objects for 8 rotational transitions of the CO molecule and thencompare our results to those of previous studies. Our expressions account for line saturation and resolving of the envelope, thereby allowing accurate determination of very high mass-loss rates. We argue that, for estimates based on a single rotational line, the CO(2-1) transition provides the most reliable mass-loss rate. The mass-loss rates calculated for the asympotic giant branch stars range from 4 \times 10^-8 M_ȯ yr^-1 up to 8 \times 10^-5 M_ȯ yr^-1. For red supergiants they reach values between 2 \times 10^-7 M_ȯ yr^-1 and 3 \times 10^-4 M_ȯ yr^-1. The estimates for the set of CO transitions allow time variability to be identified in the mass-loss rate. Possible mass-loss-rate variability is traced for 7 of the sample stars. We find a clear relation between the pulsation periods of the asympotic giant branch stars and their derived mass-loss rates, with a levelling off at ~3 \times 10^-5 M_ȯ yr^-1 for periods exceeding 850 days. Conclusions: Appendices are only available in electronic form at http://www.aanda.org
The relation between ices in the envelopes and disks surrounding young stellar objects (YSOs) and those in the quiescent interstellar medium (ISM) is investigated. For a sample of 31 stars behind ...isolated dense cores, ground-based and Spitzer spectra and photometry in the 1-25 Delta *mm wavelength range are combined. The baseline for the broad and overlapping ice features is modeled, using calculated spectra of giants, H2O ice and silicates. The adopted extinction curve is derived empirically. Its high resolution allows for the separation of continuum and feature extinction. The extinction between 13 and 25 Delta *mm is ~50% relative to that at 2.2 Delta *mm. The strengths of the 6.0 and 6.85 Delta *mm absorption bands are in line with those of YSOs. Thus, their carriers, which, besides H2O and CH3OH, may include NH+ 4, HCOOH, H2CO, and NH3, are readily formed in the dense core phase, before stars form. The 3.53 Delta *mm C-H stretching mode of solid CH3OH was discovered. The CH3OH/H2O abundance ratios of 5%-12% are larger than upper limits in the Taurus molecular cloud. The initial ice composition, before star formation occurs, therefore depends on the environment. Signs of thermal and energetic processing that were found toward some YSOs are absent in the ices toward background stars. Finally, the peak optical depth of the 9.7 Delta *mm band of silicates relative to the continuum extinction at 2.2 Delta *mm is significantly shallower than in the diffuse ISM. This extends the results of Chiar et al. to a larger sample and higher extinctions.