Search

Among the different types of evolved stars, some post-AGBs experience fast shocks that disrupt their circumstellar envelopes activating a complex chemistry, resembling that of young stellar objects, which is poorly characterized. The CHIPS project aims to understand the complex chemistry of these objects through observations in the (sub)mmand chemical modeling.

Binary interactions dominate the evolution of massive stars, but their role is less clear for low- and intermediate-mass stars. The evolution of a spherical wind from an asymptotic giant branch (AGB) star into a nonspherical planetary nebula (PN) could be due to binary interactions. We observed a sample of AGB stars with the Atacama Large Millimeter/submillimeter Array (ALMA) and found that their winds exhibit distinct nonspherical geometries with morphological similarities to planetary nebulae (PNe). We infer that the same physics shapes both AGB winds and PNe; additionally, the morphology and AGB mass-loss rate are correlated. These characteristics can be explained by binary interaction. We propose an evolutionary scenario for AGB morphologies that is consistent with observed phenomena in AGB stars and PNe., 19 pages main journal, 97 pages Supplementary Information

S-type asymptotic giant branch (AGB) stars are thought to be intermediates in the evolution of oxygen- to carbon-rich AGB stars. The chemical compositions of their circumstellar envelopes are also intermediate, but have not been studied in as much detail as their carbon- and oxygen-rich counterparts. We aim to determine the abundances of AlCl and AlF from rotational lines, which have been observed for the first time towards an S-type AGB star, W Aql. In combination with models based on PACS observations, we aim to update our chemical kinetics network based on these results. We analyse ALMA observations towards W Aql of AlCl in the ground and first two vibrationally excited states and AlF in the ground vibrational state. Using radiative transfer models, we determine the abundances and spatial abundance distributions of Al$^{35}$Cl, Al$^{37}$Cl, and AlF. We also model HCl and HF emission and compare these models to PACS spectra to constrain the abundances of these species. AlCl is found in clumps very close to the star, with emission confined within 0.1$^{\prime\prime}$ of the star. AlF emission is more extended, with faint emission extending 0.2$^{\prime\prime}$ to 0.6$^{\prime\prime}$ from the continuum peak. We find peak abundances, relative to H$_2$, of $1.7\times 10^{-7}$ for Al$^{35}$Cl, $7\times 10^{-8}$ for Al$^{37}$Cl and $1\times 10^{-7}$ for AlF. From the PACS spectra, we find abundances of $9.7\times 10^{-8}$ and $\leq 10^{-8}$, relative to H$_2$, for HCl and HF, respectively. The AlF abundance exceeds the solar F abundance, indicating that fluorine synthesised in the AGB star has already been dredged up to the surface of the star and ejected into the circumstellar envelope. From our analysis of chemical reactions in the wind, we conclude that AlF may participate in the dust formation process, but we cannot fully explain the rapid depletion of AlCl seen in the wind., Accepted for publication in Astronomy & Astrophysics

Observation of CO emission around asymptotic giant branch (AGB) stars is the primary method to determine gas mass-loss rates. While radiative transfer models have shown that molecular levels of CO can become mildly inverted, causing maser emission, CO maser emission has yet to be confirmed observationally. High-resolution observations of the CO emission around AGB stars now have the brightness temperature sensitivity to detect possible weak CO maser emission. We used high angular resolution observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the small-scale structure of CO $J=3-2$ emission around the oxygen-rich AGB star W Hya. We find CO maser emission amplifying the stellar continuum with an optical depth $\tau\approx-0.55$. The maser predominantly amplifies the limb of the star because CO $J=3-2$ absorption from the extended stellar atmosphere is strongest towards the centre of the star. The CO maser velocity corresponds to a previously observed variable component of high-frequency H$_2$O masers and with the OH maser that was identified as the amplified stellar image. This implies that the maser originates beyond the acceleration region and constrains the velocity profile since we find the population inversion primarily in the inner circumstellar envelope. We find that inversion can be explained by the radiation field at 4.6 $\mu$m and that the existence of CO maser emission is consistent with the estimated mass-loss rates for W Hya. However, the pumping mechanism requires a complex interplay between absorption and emission lines in the extended atmosphere. Excess from dust in the circumstellar envelope of W Hya is not sufficient to contribute significantly to the required radiation field at 4.6 $\mu$m. The interplay between molecular lines that cause the pumping can be constrained by future multi-level CO observations., Comment: Accepted for publication in A\&A, 14 pages, 14 figures

25 pags., 15 figs., 4 tabs., Aims. We aim to constrain the sizes of, and investigate deviations from spherical symmetry in, the CO circumstellar envelopes (CSEs) of 16 S-type stars, along with an additional 7 and 4 CSEs of C-type and M-type AGB stars, respectively. Methods. We map the emission from the CO J = 2-1 and 3-2 lines observed with the Atacama Compact Array (ACA) and its total power (TP) antennas, and fit with a Gaussian distribution in the uv- and image planes for ACA-only and TP observations, respectively. The major axis of the fitted Gaussian for the CO(2-1) line data gives a first estimate of the size of the CO-line-emitting CSE. We investigate possible signs of deviation from spherical symmetry by analysing the line profiles and the minor-to-major axis ratio obtained from visibility fitting, and by investigating the deconvolved images. Results. The sizes of the CO-line-emitting CSEs of low-mass-loss-rate (low-MLR) S-type stars fall between the sizes of the CSEs of C-stars, which are larger, and those of M-stars, which are smaller, as expected because of the differences in their respective CO abundances and the dependence of the photodissociation rate on this quantity. The sizes of the low-MLR S-type stars show no dependence on circumstellar density, as measured by the ratio of the MLR to terminal outflow velocity, irrespective of variability type. The density dependence steepens for S-stars with higher MLRs. While the CO(2-1) brightness distribution size of the low-density S-stars is in general smaller than the predicted photodissociation radius (assuming the standard interstellar radiation field), the measured size of a few of the high-density sources is of the same order as the expected photodissociation radius. Furthermore, our results show that the CO CSEs of most of the S-stars in our sample are consistent with a spherically symmetric and smooth outflow. For some of the sources, clear and prominent asymmetric features are observed which are indicative of intrinsic circumstellar anisotropy. Conclusions. As the majority of the S-type CSEs of the stars in our sample are consistent with a spherical geometry, the CO envelope sizes obtained in this paper will be used to constrain detailed radiative transfer modelling to directly determine more accurate MLR estimates for the stars in our sample. For several of our sources that present signs of deviation from spherical symmetry, further high-resolution observations would be necessary to investigate the nature of, and the physical processes behind, these asymmetrical structures. This will provide further insight into the mass-loss process and its related chemistry in S-type AGB stars., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2018.1.01434.S; ADS/JAO.ALMA#2017.1.00595.S and ADS/JAO.ALMA#2012.1.00524.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. MA acknowledges support from the Nordic ALMA Regional Centre (ARC) node based at Onsala Space Observatory. The Nordic ARC node is funded through Swedish Research Council grant No 2017-00648. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme under grant agreements No. 883867 [EXWINGS] and 730562 [RadioNet].

We report the detection of two lines at millimetre wavelengths towards the immediate surroundings of the post-giant and most likely post-common-envelope star HD101584 using high-angular-resolution ALMA observations. The circumstellar environment of this object is rich in different molecular species, but we find no viable identifications in terms of molecular lines. We aim to determine whether or not these lines can be attributed to the Rydberg transitions -- X30alpha and X26alpha -- of neutral atoms of elements heavier than carbon. A simple model in strict local thermodynamic equilibrium for a warm-gas environment of the moderate-temperature star (T_eff about 8500 K) was constructed to corroborate our findings. A geometrically thin, disc-like geometry seen face-on was chosen and a distance of 1 kpc. The observed flux densities of the lines and the continuum at 232 and 354 GHz can be reproduced using 10^(-3) M_sun of gas at a temperature of about 2800 K and a hydrogen density of about 10^(12) cm(-3), assuming solar abundances for the elements. The gas lies within a distance of about 5 au from the star (assuming a distance of 1 kpc). The ionisation fraction is low, about 3x10^(-5). The origin of such a region is not clear, but it may be related to a common-envelope-evolution phase. With these conditions, the line emissions are dominated by Rydberg transitions within the stable isotopes of Mg. A turbulent velocity field in the range 5.5 - 7.5 km s^(-1) is required to fit the Gaussian line shapes. An upper limit to the average magnetic field in the line-emitting region of 1G is set using the Zeeman effect in these lines. We speculate that Rydberg transitions of heavy elements may be an interesting probe for the close-in environments of other moderate-temperature objects like AGB stars, red supergiants, yellow hypergiants, and binaries of various types., Comment: 19 pages, 6 figures

One of the most poorly understood stellar evolutionary paths is that of binary systems undergoing common-envelope evolution, when the envelope of a giant star engulfs the orbit of a companion. The interaction that ensues leads to a great variety of astrophysical systems and associated phenomena, but happens over a very short timescale. Unfortunately, direct empirical studies of this momentous and complex phase are difficult at present because few objects experiencing, or having just experienced, common-envelope evolution are known. Here we present Atacama Large Millimeter/submillimeter Array observations of minor CO isotopologues towards a sample of sources known as water fountains, which reveal that almost all of them recently lost a substantial fraction of their initial mass over a timescale of less than a few tens to a few hundreds of years. The only known mechanism able to explain such rapid mass ejection, corresponding to a large fraction of the stellar mass, is the common-envelope evolution. A stellar population analysis shows that the number of water-fountain sources in the Milky Way is comparable to the expected number of common-envelope events that involve low-mass evolved stars. Thus, the known sample of water-fountain sources accounts for a large fraction of the systems undergoing a common-envelope phase in our Galaxy. As one of the distinguishing characteristics of water-fountain sources is their fast bipolar outflow, we conclude that outflows and jets play an important role right before, during or immediately after the common-envelope phase. © 2021, The Author(s), under exclusive licence to Springer Nature Limited., This study made use of the following ALMA data: ADS/JAO.ALMA 2018.1.00250.S and 2016.1.01032.S. ALMA is a partnership of the European Southern Observatory (ESO) (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan) and KASI (Republic of Korea), and in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The authors acknowledge support from the Nordic ALMA Regional Centre (ARC) node based at the Onsala Space Observatory. The Nordic ARC node and Swedish observations on APEX are funded through the Swedish Research Council grant no. 2017-00648. This publication is based on data acquired with the Atacama Pathfinder Experiment (APEX) under programme ID 0104.F-9310(A). APEX is a collaboration between the Max-Planck-Institut für Radioastronomie, the ESO and the Onsala Space Observatory. T.K. is supported by Swedish Research Council starting grant no. 2019-03777. W.H.T.V. and T.K. are supported by Swedish Research Council grant no. 2014-05713. C.S.C. is supported by the Spanish MICINN through grant no. PID2019-105203GB-C22. J.F.G. is supported by MCIU-AEI through the ‘Center of Excellence Severo Ochoa’ award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709), from grants AYA2017-84390-C2-1-R and PID2020-114461GB-I00 of AEI (10.13039/501100011033), co-funded by FEDER, and from the Amanogawa Galaxy Astronomy Research Center (AGARC). H.I. is supported by the MEXT KAKENHI programme (16H02167) and i-LINK+2019 programme at IAA/CSIC. J.F.G. and H.I. are supported by the Invitation Program for Foreign Researchers of the Japan Society for Promotion of Science (grant no. S14128).

28 pags., 10 figs., 7 tabs., Context. This is the first publication from the DEATHSTAR project. The overall goal of the project is to reduce the uncertainties of the observational estimates of mass-loss rates from evolved stars on the Asymptotic Giant Branch (AGB). Aim. The aim in this first publication is to constrain the sizes of the 12CO emitting region from the circumstellar envelopes around 42 mostly southern AGB stars, of which 21 are M-Type and 21 are C-Type, using the Atacama Compact Array (ACA) at the Atacama Large Millimeter/submillimeter Array. The symmetry of the outflows is also investigated. Methods. Line emission from 12CO J = 2→1 and 3→2 from all of the sources were mapped using the ACA. In this initial analysis, the emission distribution was fit to a Gaussian distribution in the uv-plane. A detailed radiative transfer analysis will be presented in a future publication. The major and minor axis of the best-fit Gaussian at the line center velocity of the 12CO J = 2→1 emission gives a first indication of the size of the emitting region. Furthermore, the fitting results, such as the Gaussian major and minor axis, center position, and the goodness of fit across both lines, constrain the symmetry of the emission distribution. For a subsample of sources, the measured emission distribution is compared to predictions from previous best-fit radiative transfer modeling results. Results. We find that the CO envelope sizes are, in general, larger for C-Type than for M-Type AGB stars, which is as expected if the CO/H2 ratio is larger in C-Type stars. Furthermore, the measurements show a relation between the measured (Gaussian) 12CO J = 2→1 size and circumstellar density that, while in broad agreement with photodissociation calculations, reveals large scatter and some systematic differences between the different stellar types. For lower mass-loss-rate irregular and semi-regular variables of both M-and C-Type AGB stars, the 12CO J = 2→1 size appears to be independent of the ratio of the mass-loss rate to outflow velocity, which is a measure of circumstellar density. For the higher mass-loss-rate Mira stars, the 12CO J = 2→1 size clearly increases with circumstellar density, with larger sizes for the higher CO-Abundance C-Type stars. The M-Type stars appear to be consistently smaller than predicted from photodissociation theory. The majority of the sources have CO envelope sizes that are consistent with a spherically symmetric, smooth outflow, at least on larger scales. For about a third of the sources, indications of strong asymmetries are detected. This is consistent with what was found in previous interferometric investigations of northern sources. Smaller scale asymmetries are found in a larger fraction of sources. Conclusions. These results for CO envelope radii and shapes can be used to constrain detailed radiative transfer modeling of the same stars so as to determine mass-loss rates that are independent of photodissociation models. For a large fraction of the sources, observations at higher spatial resolution will be necessary to deduce the nature and origin of the complex circumstellar dynamics revealed by our ACA observations., This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.2.00025.S and #2017.1.00595.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www. cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/ consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

The mass-loss mechanism of asymptotic giant branch stars has long been thought to rely on two processes: stellar pulsations and dust formation. The details of the mass-loss mechanism have remained elusive, however, because of the overall complexity of the dust formation process in the very dynamical pulsation-enhanced atmosphere. Recently, our understanding of AGB stars and the associated mass loss has evolved significantly, thanks both to new instruments which allow sensitive and high-angular-resolution observations and the development of models for the convective AGB envelopes and the dust formation process. ALMA and SPHERE/ZIMPOL on the VLT have been very important instruments in driving this advance in the last few years by providing high-angular resolution images in the sub-mm and visible wavelengths, respectively. I will present observations obtained using these instruments at the same epoch (2.5 weeks apart) of the AGB star Mira that resolve even the stellar disk. The ALMA data reveals the distribution and dynamics of the gas around the star, while the polarised light imaged using SPHERE shows the distribution of the dust grains expected to drive the outflows. Moreover, the observations show a central source surrounded by asymmetric distributions of gas and dust, with complementary structures seen in the two components. We model the observed CO v = 1, J = 3−2 line to determine the density, temperature and velocity of gas close to the star. This model is then used to estimate the abundance of AlO. Our results show that only a very small fraction of aluminium (≲0.1%) is locked in AlO molecules. We also calculate models to fit the observed polarised light based on the gas densities we find. The low level of visible-light polarisation detected using ZIMPOL implies that, at the time of the observations, aluminium atoms are either not efficiently depleted into dust or the aluminium-oxide grains are relatively small (≲0.02μm).

Isotopologue line intensity ratios of circumstellar molecules have been widely used to trace the photospheric elemental isotopic ratios of evolved stars. However, depending on the molecular species and the physical conditions of the environment, the circumstellar isotopologue ratio may deviate considerably from the stellar atmospheric value. In this paper, we aim to examine how the CO and HCN abundance ratios vary radially due to chemical reactions in the outflows of AGB stars and the effect of excitation and optical depth on the resulting line intensity ratios. We find that the circumstellar 12CO/13CO can deviate from its atmospheric value by up to 25-94% and 6-60% for C- and O-type CSEs, respectively. We show that variations of the intensity of the ISRF and the gas kinetic temperature can significantly influence the CO isotopologue ratio in the outer CSEs. On the contrary, the H12CN/H13CN ratio is stable for all tested mass-loss rates. The RT modeling shows that the integrated line intensity ratio of CO of different rotational transitions varies significantly for stars with intermediate mass-loss rates due to combined chemical and excitation effects. In contrast, the excitation conditions for the both HCN isotopologues are the same. We demonstrate the importance of using the isotopologue abundance profiles from chemical models as inputs to RT models in the interpretation of isotopologue observations. Previous studies of CO isotopologue ratios are based on multi-transition data for individual sources and it is difficult to estimate the errors in the reported values due to assumptions that are not entirely correct according to this study. If anything, previous studies may have overestimated the circumstellar 12CO/13CO abundance ratio. The use of the HCN as a tracer of C isotope ratios is affected by fewer complicating problems, provided one accounts corrections for high optical depths., Comment: 14 pages, 11 figures

The initial conditions for the mass loss during the asymptotic giant branch (AGB) phase are set in their extended atmospheres, where, among others, convection and pulsation driven shocks determine the physical conditions. High resolution observations of AGB stars at (sub)millimetre wavelengths can now directly determine the morphology, activity, density, and temperature close to the stellar photosphere. We used the Atacama Large Millimeter/submillimeter Array (ALMA) high angular resolution observations to resolve the extended atmospheres of four of the nearest AGB stars: W Hya, Mira A, R Dor and R Leo. We interpreted the observations using a parameterised atmosphere model. We resolve all four AGB stars and determine the brightness temperature structure between $1$ and $2$ stellar radii. For W Hya and R Dor we confirm the existence of hotspots with brightness temperatures $>3000$ to $10000$~K. All four stars show deviations from spherical symmetry. We find variations on a timescale of days to weeks, and for R Leo we directly measure an outward motion of the millimetre wavelength surface with a velocity of at least $10.6\pm1.4$~km~s$^{-1}$. For all objects but W Hya we find that the temperature-radius and size-frequency relations require the existence of a (likely inhomogeneous) layer of enhanced opacity. The ALMA observations provide a unique probe of the structure of the extended AGB atmosphere. We find highly variable structures of hotspots and likely convective cells. In the future, these observations can be directly compared to multi-dimensional chromosphere and atmosphere models that determine the temperature, density, velocity, and ionisation structure between the stellar photosphere and the dust formation region. However, our results show that for the best interpretation, both very accurate flux calibration and near-simultaneous observations are essential., 18 pages, 13 figures, Accepted to A&A, final version after language editing

We report the detection and investigate the properties of high-excitation lambda-doubling line emission of hydroxyl (OH) detected towards three asymptotic giant branch (AGB) stars (W Hya, R Dor, and IK Tau) using ALMA. The OH lines are produced very close to the central stars and seem optically thin and with no maser effect. We analyse the molecular excitation using population diagrams and find rotational temperatures of $\sim 2500$ K and column densities of $\sim 10^{19}$ cm$^{-2}$ for both W Hya and R Dor. For W Hya, we observe emission from vibrationally excited H2O arising from the same region as the OH emission. Moreover, CO $v = 1; J = 3 - 2$ emission also shows a brightness peak in the same region. Considering optically thin emission and the rotational temperature derived for OH, we find a CO column density $\sim 15$ times higher than that of OH, within an area of ($92 \times 84$) mas$^2$ centred on the OH emission peak. These results should be considered tentative because of the simple methods employed. The observed OH line frequencies differ significantly from the predicted transition frequencies in the literature, and provide the possibility of using OH lines observed in AGB stars to improve the accuracy of the Hamiltonian used for the OH molecule. We predict stronger OH lambda-doubling lines at millimetre wavelengths than those we detected. These lines will be a good probe of shocked gas in the extended atmosphere and are possibly even suitable as probes of the magnetic field in the atmospheres of close-by AGB stars through the Zeeman effect., Published in Astronomy and Astrophysics Letters

Context. The condensation of inorganic dust grains in the winds of evolved stars is poorly understood. As of today, it is not yet known which molecular clusters form the first dust grains in oxygen-rich (C/O < 1) asymptotic giant branch (AGB) winds. Aluminium oxides and iron-free silicates are often put forward as promising candidates for the first dust seeds. Aims. We aim to constrain the dust formation histories in the winds of oxygen-rich AGB stars. Methods. We obtained Atacama Large Millimeter/sub-millimeter array (ALMA) observations with a spatial resolution of 120 150 mas tracing the dust formation region of the low mass-loss rate AGB star, R Dor, and the high mass-loss rate AGB star, IK Tau. We detected emission line profiles of AlO, AlOH, and AlCl in the ALMA data and used these line profiles to derive a lower limit of atomic aluminium incorporated in molecules. This constrains the aluminium budget that can condense into grains. Results. Radiative transfer models constrain the fractional abundances of AlO, AlOH, and AlCl in IK Tau and R Dor. We show that the gas-phase aluminium chemistry is completely different in both stars with a remarkable difference in the AlO and AlOH abundance stratification. The amount of aluminium locked up in these three molecules is small, 1:110 w.r.t. H2, for both stars, i.e. only 2%of the total aluminium budget. An important result is that AlO and AlOH, which are the direct precursors of alumina (Al2O3) grains, are detected well beyond the onset of the dust condensation, which proves that the aluminium oxide condensation cycle is not fully effcient. The ALMA observations allow us to quantitatively assess the current generation of theoretical dynamical-chemical models for AGB winds. We discuss how the current proposed scenario of aluminium dust condensation for low mass-loss rate AGB stars within a few stellar radii from the star, in particular for R Dor and W Hya, poses a challenge if one wishes to explain both the dust spectral features in the spectral energy distribution (SED) in interferometric data and in the polarized light signal. In particular, the estimated grain temperature of Al2O3 is too high for the grains to retain their amorphous structure. We advocate that large gas-phase (Al2O3)n clusters (n > 34) can be the potential agents of the broad 11 m feature in the SED and in the interferometric data and we propose potential formation mechanisms for these large clusters.

We report the serendipitous detection with ALMA of the vibrationally-excited pure-rotational CO transition $v=1, J=3-2$ towards five asymptotic giant branch (AGB) stars, $o$ Cet, R Aqr, R Scl, W Aql, and $\pi^1$ Gru. The observed lines are formed in the poorly-understood region located between the stellar surface and the region where the wind starts, the so-called warm molecular layer. We successfully reproduce the observed lines profiles using a simple model. We constrain the extents, densities, and kinematics of the region where the lines are produced. R Aqr and R Scl show inverse P-Cygni line profiles which indicate infall of material onto the stars. The line profiles of $o$ Cet and R Scl show variability. The serendipitous detection towards these five sources shows that vibrationally-excited rotational lines can be observed towards a large number of nearby AGB stars using ALMA. This opens a new possibility for the study of the innermost regions of AGB circumstellar envelopes., Comment: 6 pages, 2 figures, 2 tables, 2016MNRAS.463L..74K

Context: The observed orbital characteristics of post-asymptotic giant branch (post-AGB) and post-red giant branch (post-RGB) binaries are not understood. We suspect that the missing ingredients to explain them probably lie in the continuous interaction of the central binary with its circumstellar environment. Aims: We aim at studying the circumbinary material in these complex systems by investigating the connection between the innermost and large-scale structures. Methods: We perform high-angular resolution observations in the near-infrared continuum of HD101584, which has a complex structure as seen at millimeter wavelengths with a disk-like morphology and a bipolar outflow due to an episode of strong binary interaction. To account for the complexity of the target we first perform an image reconstruction and use this result to fit a geometrical model to extract the morphological and thermal features of the environment. Results: The image reveals an unexpected double-ring structure. We interpret the inner ring to be produced by emission from dust located in the plane of the disk and the outer ring to be produced by emission from dust that is located 1.6[D/1kpc] au above the disk plane. The inner ring diameter (3.94[D/1kpc] au), and temperature (T=1540$\pm$10K) are compatible with the dust sublimation front of the disk. The origin of the out-of-plane ring (with a diameter of 7.39[D/1kpc] au and a temperature of 1014$\pm10$K) could be due to episodic ejection or a dust condensation front in the outflow. Conclusion: The observed outer ring is possibly linked with the blue-shifted side of the large scale outflow seen by ALMA and is tracing its launching location to the central star. Such observations give morphological constraints on the ejection mechanism. Additional observations are needed to constrain the origin of the out-of-plane structure., Comment: Accepted to A&A. 14 pages, 13 figures

The symbiotic star R Aqr is part of a small sample of binary AGB stars observed with ALMA. The sample stars are: R Aqr, Mira, W Aql, and pi1 Gru. The sample covers a range in binary separation and wind properties, where R Aqr is the source with the smallest separation. The R Aqr binary pair consists of an M-type AGB star and a white dwarf at a separation of 45 mas, equivalent to about 10 AU at 218 pc. The aim of the ALMA study is to investigate the dependence of the wind shaping on the binary separation and to provide constraints for hydrodynamical binary interaction models. R Aqr is particularly interesting as the source with the smallest separation and a complex circumstellar environment that is strongly affected by the interaction between the two stars and by the high-energy radiation resulting from this interaction and from the hot white dwarf companion. The CO(J=3-2) line emission has been observed at 0.5" spatial resolution. The CO envelope around the binary pair is marginally resolved, showing what appears to be a rather complex distribution. The outer radius of the CO emitting region is estimated from the data and found to be about a factor of 10 larger than previously thought. This implies an average mass-loss rate during the past 100 yr of 2x10-7 Msun/yr, a factor of 45 less than previous estimates. The channel maps are presented and the molecular gas distribution is discussed and set into the context of what was previously known about the system. Additional molecular line emission detected within the bandwidth covered by the ALMA observations is also presented. Because of the limited extent of the emission, firm conclusions about the dynamical evolution of the system will have to wait for higher spatial resolution observations. However, the data presented here support the assumption that the mass-loss rate from the Mira star strongly varies and is focused on the orbital plane., 6 pages, 4 figures, accepted in A&A

(abridged) Our aim is to determine the radial abundance profile of SiO and HCN throughout the stellar outflow of R Dor, an oxygen-rich AGB star with a low mass-loss rate. We have analysed molecular transitions of CO, SiO, and HCN measured with the APEX telescope and all three instruments on the Herschel Space Observatory, together with literature data. Photometric data and the infrared spectrum measured by ISO-SWS were used to constrain the dust component of the outflow. Using both continuum and line radiative transfer methods, a physical envelope model of both gas and dust was established. We have performed an analysis of the SiO and HCN molecular transitions in order to calculate their abundances. We have obtained an envelope model that describes the dust and the gas in the outflow, and determined the abundance of SiO and HCN throughout the region of the outflow probed by our molecular data. For SiO, we find that the initial abundance lies between $5.5 \times 10^{-5}$ and $6.0 \times 10^{-5}$ w.r.t. H$_2$. The abundance profile is constant up to $60\ \pm 10\ R_*$, after which it declines following a Gaussian profile with an $e$-folding radius of $3.5 \pm 0.5 \times 10^{13}$ cm. For HCN, we find an initial abundance of $5.0 \times 10^{-7}$ w.r.t. H$_2$. The Gaussian profile that describes the decline starts at the stellar surface and has an $e$-folding radius $r_e$ of $1.85 \pm 0.05 \times 10^{15}$ cm. We cannot to unambiguously identify the mechanism by which SiO is destroyed at $60\ \pm 10\ R_*$. The initial abundances found are larger than previously determined (except for one previous study on SiO), which might be due to the inclusion of higher-$J$ transitions. The difference in abundance for SiO and HCN compared to high mass-loss rate Mira star IK Tau might be due to different pulsation characteristics of the central star and/or a difference in dust condensation physics., Comment: 16 pages, 13 figures, 7 tables, accepted for publication in Astronomy & Astrophysics

High resolution observations of the extended atmospheres of asymptotic giant branch (AGB) stars can now directly confront the theories that describe stellar mass loss. Using Atacama Large Millimeter/submillimeter Array (ALMA) high angular resolution ($30\times42$~mas) observations we have, for the first time, resolved stellar rotation of an AGB star, R~Dor. We measure an angular rotation velocity of $\omega_R\sin{i}=(3.5\pm0.3)\times10^{-9}$~rad~s$^{-1}$ which indicates a rotational velocity of $|\upsilon_{\rm rot}\sin{i}|=1.0\pm0.1$~km~s$^{-1}$ at the stellar surface ($R_*=31.2$~mas at $214$~GHz). The rotation axis projected on the plane of the sky has a position angle $\Phi=7\pm6^\circ$. We find that the rotation of R Dor is two orders of magnitude faster than expected for a solitary AGB star that will have lost most of its angular momentum. Its rotational velocity is consistent with angular momentum transfer from a close companion. As a companion has not been directly detected we thus suggest R~Dor has a low-mass, close-in, companion. The rotational velocity approaches the critical velocity, set by the local sound speed in the extended envelope, and is thus expected to affect the mass loss characteristics of R~Dor., Comment: 12 pages, 15 figures, Accepted for publication in Astronomy & Astrophysics

Our current understanding of the chemistry and mass-loss processes in solar-like stars at the end of their evolution depends critically on the description of convection, pulsations and shocks in the extended stellar atmosphere. Three-dimensional hydrodynamical stellar atmosphere models provide observational predictions, but so far the resolution to constrain the complex temperature and velocity structures seen in the models has been lacking. Here we present submillimeter continuum and line observations that resolve the atmosphere of the asymptotic giant branch star W Hya. We show that hot gas with chromospheric characteristics exists around the star. Its filling factor is shown to be small. The existence of such gas requires shocks with a cooling time larger than commonly assumed. A shocked hot layer will be an important ingredient in the models of stellar convection, pulsation and chemistry that underlie our current understanding of the late stages of stellar evolution., 30 pages, 9 figures, including Supplementary information. Author manuscript version before editorial/copyediting by Nature Astronomy. Journal version available via http://rdcu.be/xUWV

We present ALMA band 7 data of the extreme OH/IR star, OH 26.5+0.6. In addition to lines of CO and its isotopologues, the circumstellar envelope also exhibits a number of emission lines due to metal-containing molecules, e.g., NaCl and KCl. A lack of C18O is expected, but a non-detection of C17O is puzzling given the strengths of H217O in Herschel spectra of the star. However, a line associated with Si17O is detected. We also report a tentative detection of a gas-phase emission line of MgS. The ALMA spectrum of this object reveals intriguing features which may be used to investigate chemical processes and dust formation during a high mass-loss phase.

We have performed a study of the characteristics of the circumstellar environment of the binary object HD101584, that provides information on a likely evolutionary scenario. We have obtained and analysed ALMA observations, complemented with observations using APEX, of a large number of molecular lines. An analysis of the spectral energy distribution has also been performed. Emissions from 12 molecular species (not counting isotopologues) have been observed, and most of them mapped with angular resolutions in the range 0.1" to 0.6". Four circumstellar components are identified: i) a central compact source of size 0.15", ii) an expanding equatorial density enhancement (a flattened density distribution in the plane of the orbit) of size 3", iii) a bipolar high-velocity outflow (150 km/s), and iv) an hourglass structure. The outflow is directed almost along the line of sight. There is evidence of a second bipolar outflow. The mass of the circumstellar gas is 0.5[D/1 kpc]^2 Msun, about half of it lies in the equatorial density enhancement. The dust mass is 0.01[D/1 kpc]^2 Msun, and a substantial fraction of this is in the form of large-sized, up to 1 mm, grains. The estimated kinetic age of the outflow is 770[D/1 kpc] yr. The kinetic energy and the scalar momentum of the accelerated gas are estimated to be 7x10^(45)[D/1 kpc]^2 erg and 10^(39)[D/1 kpc]^2 g cm/s, respectively. We provide good evidence that the binary system HD101584 is in a post-common-envelope-evolution phase, that ended before a stellar merger. Isotope ratios combined with stellar mass estimates suggest that the primary star's evolution was terminated already on the first red giant branch (RGB). Most of the energy required to drive the outflowing gas was probably released when material fell towards the companion., Comment: Accepted for publication in A&A

Aims. Polarisation observations of circumstellar dust and molecular (thermal and maser) lines provide unique information about dust properties and magnetic fields in circumstellar envelopes of evolved stars. Methods. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 5 science verification observations of the red supergiant VY CMa to study the polarisation of SiO thermal/maser lines and dust continuum at ∼1.7 mm wavelength. We analyse both linear and circular polarisation and derive the magnetic field strength and structure, assuming the polarisation of the lines originates from the Zeeman effect, and that of the dust originates from aligned dust grains. We also discuss other effects that could give rise to the observed polarisation. Results. We detect, for the first time, significant polarisation (∼3%) of the circumstellar dust emission at millimeter wavelengths. The polarisation is uniform with an electric vector position angle of ∼8°. Varying levels of linear polarisation are detected for the J = 4 - 328SiO v = 0, 1, 2, and 29SiO v = 0, 1 lines, with the strongest polarisation fraction of ∼30% found for the 29SiO v = 1 maser. The linear polarisation vectors rotate with velocity, consistent with earlier observations. We also find significant (up to ∼1%) circular polarisation in several lines, consistent with previous measurements. We conclude that the detection is robust against calibration and regular instrumental errors, although we cannot yet fully rule out non-standard instrumental effects. Conclusions. Emission from magnetically aligned grains is the most likely origin of the observed continuum polarisation. This implies that the dust is embedded in a magnetic field >13 mG. The maser line polarisation traces the magnetic field structure. The magnetic field in the gas and dust is consistent with an approximately toroidal field configuration, but only higher angular resolution observations will be able to reveal more detailed field structure. If the circular polarisation is due to Zeeman splitting, it indicates a magnetic field strength of ∼1-3 Gauss, consistent with previous maser observations.

We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We present observations in filters V, cntH$\alpha$, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntH$\alpha$, we are able to see variability on the stellar photosphere. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntH$\alpha$, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntH$\alpha$ show that the morphology of R Dor changes completely in a timespan of 48 days to a more axisymmetric and compact configuration. The polarised intensity is asymmetric in all epochs and varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, $\rho(r)=\rho_\circ r^{-n}$. On average, we find exponents of $- 4.5 \pm 0.5$ that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of $\sim 2\times10^{-4}$ at 5.0 $R_\star$. Given the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of $\sim 3.3\times10^{-4}$. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light., Comment: 10 pages, 8 figures, 4 tables

[Context] Context. The recent detection of warm H2O vapor emission from the outflows of carbon-rich asymptotic giant branch (AGB) stars challenges the current understanding of circumstellar chemistry. Two mechanisms have been invoked to explain warm H2O vapor formation. In the first, periodic shocks passing through the medium immediately above the stellar surface lead to H2O formation. In the second, penetration of ultraviolet interstellar radiation through a clumpy circumstellar medium leads to the formation of H2O molecules in the intermediate wind., [Aims] We aim to determine the properties of H2O emission for a sample of 18 carbon-rich AGB stars and subsequently constrain which of the above mechanisms provides the most likely warm H2O formation pathway., [Methods] Using far-infrared spectra taken with the PACS instrument onboard the Herschel telescope, we combined two methods to identify H2O emission trends and interpreted these in terms of theoretically expected patterns in the H2O abundance. Through the use of line-strength ratios, we analyzed the correlation between the strength of H2O emission and the mass-loss rate of the objects, as well as the radial dependence of the H2O abundance in the circumstellar outflow per individual source. We computed a model grid to account for radiative-transfer effects in the line strengths., [Results] We detect warm H2O emission close to or inside the wind acceleration zone of all sample stars, irrespective of their stellar or circumstellar properties. The predicted H2O abundances in carbon-rich environments are in the range of 10-6 up to 10-4 for Miras and semiregular-a objects, and cluster around 10-6 for semiregular-b objects. These predictions are up to three orders of magnitude greater than what is predicted by state-of-the-art chemical models. We find a negative correlation between the H2O/CO line-strength ratio and gas mass-loss rate for Ṁg> 5 × 10-7 M⊙ yr-1, regardless of the upper-level energy of the relevant transitions. This implies that the H2O formation mechanism becomes less efficient with increasing wind density. The negative correlation breaks down for the sources of lowest mass-loss rate, the semiregular-b objects., [Conclusions] Observational constraints suggest that pulsationally induced shocks play an important role in warm H2O formation in carbon-rich AGB stars, although photodissociation by interstellar UV photons may still contribute. Both mechanisms fail in predicting the high H2O abundances we infer in Miras and semiregular-a sources, while our results for the semiregular-b objects are inconclusive., R.L. acknowledges financial support from the Fund for Scientific Research – Flanders (FWO) under grant number ZKB5757-04-W01, from the Department of Physics and Astronomy of the KU Leuven, and from the Belgian Federal Science Policy Office via the PRODEX Program of ESA under grant number C90371. L.D. acknowledges financial support from the FWO. P.R., N.C., J.D., J.B., M.G., and B.V. acknowledge support from the Belgian Federal Science Policy Office via the PRODEX Programme of ESA. F.K. is supported by the FWF project P23586 and the ffg ASAP project HIL. E.G-A is a Research Associate at the Harvard-Smithsonian CfA, and thanks the Spanish Ministerio de Economía y Competitividad for support under projects AYA2010-21697-C05-0 and FIS2012-39162-C06-01 and partial support from NHSC/JPL RSA 1455432. This development has been supported by the funding agencies BMVIT (Austria), ESA-PRODEX (Belgium), CEA/CNES (France), DLR (Germany), ASI (Italy), and CICT/MCT (Spain). For the computations we used the infrastructure of the VSC – Flemish Supercomputer Center, funded by the Hercules Foundation and the Flemish Government – department EWI.

Using both dynamical and chemical modelling, we derive an accurate abundance profile for the molecule SiO in the stellar wind of R Dor, an O-rich AGB star. SiO plays a key role in the dust formation process in O-rich AGB stars. This method will be applied to additional molecules, with the aim to achieve a detailed overview of the molecular abundance pattern in the wind of R Dor.

Context. Asymptotic giant branch (AGB) stars experience strong mass loss driven by dust particles formed in the upper atmospheres. The dust is released into the interstellar medium, and replenishes galaxies with synthesised material from the star. The dust grains further act as seeds for continued dust growth in the diffuse medium of galaxies. As such, understanding the properties of dust produced during the asymptotic giant branch phase of stellar evolution is important for understanding the evolution of stars and galaxies. Recent observations of the carbon AGB star R Scl have shown that observations at far-infrared and submillimetre wavelengths can effectively constrain the grain sizes in the shell, while the total mass depends on the structure of the grains (solid vs. hollow or fluffy). Aims. We aim to constrain the properties of the dust observed in the submillimetre in the detached shells around the three carbon AGB stars U Ant, DR Ser, and V644 Sco, and to investigate the constraints on the dust masses and grain sizes provided by far-infrared and submm observations. Methods. We observed the carbon AGB stars U Ant, DR Ser, and V644 Sco at 870 μm using LABOCA on APEX. Combined with observations from the optical to far-infrared, we produced dust radiative transfer models of the spectral energy distributions (SEDs) with contributions from the stars, present-day mass-loss and detached shells. We assume spherical, solid dust grains, and test the effect of different total dust masses and grain sizes on the SED, and attempted to consistently reproduce the SEDs from the optical to the submm. Results. We derive dust masses in the shells of a few 10−5 M ⊙. The best-fit grain radii are comparatively large, and indicate the presence of grains between 0.1 μm and 2 μm. The LABOCA observations suffer from contamination from 12CO (3 − 2), and hence gives fluxes that are higher than the predicted dust emission at submm wavelengths. We investigate the effect on the best-fitting models by assuming different degrees of contamination and show that far-infrared and submillimetre observations are important to constrain the dust mass and grain sizes in the shells. Conclusions. Spatially resolved observations of the detached shells in the far-infrared and submillimetre effectively constrain the temperatures in the shells, and hence the grain sizes. The dust mass is also constrained by the observations, but additional observations are needed to constrain the structure of the grains.

Context. On the asymptotic giant branch (AGB), Sun-like stars lose a large portion of their mass in an intensive wind and enrich the surrounding interstellar medium with nuclear processed stellar material in the form of molecular gas and dust. For a number of carbon-rich AGB stars, thin detached shells of gas and dust have been observed. These shells are formed during brief periods of increased mass loss and expansion velocity during a thermal pulse, and open up the possibility to study the mass-loss history of thermally pulsing AGB stars. Aims. We study the properties of dust grains in the detached shell around the carbon AGB star R Scl and aim to quantify the influence of the dust grain properties on the shape of the spectral energy distribution (SED) and the derived dust shell mass. Methods. We modelled the SED of the circumstellar dust emission and compared the models to observations, including new observations of Herschel/PACS and SPIRE (infrared) and APEX/LABOCA (sub-millimeter). We derived present-day mass-loss rates and detached shell masses for a variation of dust grain properties (opacities, chemical composition, grain size, and grain geometry) to quantify the influence of changing dust properties to the derived shell mass. Results. The best-fitting mass-loss parameters are a present-day dust mass-loss rate of 2 × 10−10 M⊙ yr−1 and a detached shell dust mass of (2.9 ± 0.3) × 10−5 M⊙. Compared to similar studies, the uncertainty on the dust mass is reduced by a factor of 4. We find that the size of the grains dominates the shape of the SED, while the estimated dust shell mass is most strongly affected by the geometry of the dust grains. Additionally, we find a significant sub-millimeter excess that cannot be reproduced by any of the models, but is most likely not of thermal origin.

Low- and intermediate-mass stars go through a period of intense mass-loss at the end of their lives in a phase known as the asymptotic giant branch (AGB). During the AGB a significant fraction of their initial mass is expelled in a stellar wind. This process controls the final stages of their evolution and contributes to the chemical evolution of galaxies. However, the wind-driving mechanism of AGB stars is not yet well understood, especially so for oxygen-rich sources. Characterizing both the present-day mass-loss and wind structure and the evolution of the mass-loss rate of such stars is paramount to advancing our understanding of this processes. We modelled the dust envelope of W Hya using an advanced radiative transfer code. The dust model was analysed in the light of a previously calculated gas-phase wind model and compared to measurements available in the literature, such as infrared spectra, infrared images, and optical scattered light fractions. We find that the dust spectrum of W Hya can partly be explained by a gravitationally bound dust shell that probably is responsible for most of the amorphous Al$_2$O$_3$ emission. The composition of the large ($\sim$\,0.3\,$\mu$m) grains needed to explain the scattered light cannot be constrained, but probably is dominated by silicates. Silicate emission in the thermal infrared was found to originate from beyond 40 AU from the star and we find that they need to have substantial near-infrared opacities to be visible at such large distances. The increase in near-infrared opacity of the dust at these distances roughly coincides with a sudden increase in expansion velocity as deduced from the gas-phase CO lines. Finally, the recent mass loss of W Hya is confirmed to be highly variable and we identify a strong peak in the mass-loss rate that occurred about 3500 years ago and lasted for a few hundred years., Comment: 15 pages, 13 figures

It is yet well understood how mass-loss rates from evolved stars depend on metallicities. With a half of the solar metallicity and the distance of only 50 kpc, the evolved stars of the Large Magellanic Cloud (LMC) are an ideal target for studying mass loss at low metallicity. We have obtained spectra of red-supergiants in the LMC, using the Hershel Space Observatory, detecting CO thermal lines fro J=6–5 up to 15–14 lines. Modelling CO lines with non-LTE Radiative transfer code suggests that CO lines intensities can be well explained with high gas-to-dust ratio, with no obvious reduction in mass-loss rate at the LMC. We conclude that the luminosities of the stars are primary factors on mass-loss rates, rather than the metallicity.

Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. We modeled the CO emission from the wind of the low mass-loss rate oxygen-rich AGB star W Hya using data obtained by the HIFI, PACS, and SPIRE instruments onboard the Herschel Space Observatory and ground-based telescopes. $^{12}$CO and $^{13}$CO lines are used to constrain the intrinsic $^{12}$C/$^{13}$C ratio from resolved HIFI lines. The acceleration of the outflow up to about 5.5 km/s is quite slow and can be represented by a beta-type velocity law with index 5. Beyond this point, acceleration up the terminal velocity of 7 km/s is faster. Using the J=10-9, 9-8, and 6-5 transitions, we find an intrinsic $^{12}$C/$^{13}$C ratio of $18\pm10$ for W Hya, where the error bar is mostly due to uncertainties in the $^{12}$CO abundance and the stellar flux around 4.6 $\mu$m. To match the low-excitation CO lines, these molecules need to be photo-dissociated at about 500 stellar radii. The radial dust emission intensity profile measured by PACS images at 70 $\mu$m shows substantially stronger emission than our model predicts beyond 20 arcsec. The initial slow acceleration of the wind implies inefficient wind driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for W Hya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M$_\odot$ or more. However, the uncertainty in the main-sequence mass derived based on this isotopologic ratio is large., Comment: 13 pages, 13 figures

Context. The evolution of low- and intermediate-mass stars on the asymptotic giant branch (AGB) is mainly controlled by the rate at which these stars lose mass in a stellar wind. Understanding the driving mechanism and strength of the stellar winds of AGB stars and the processes enriching their surfaces with products of nucleosynthesis are paramount to constraining AGB evolution and predicting the chemical evolution of galaxies.Aims. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the 12CO molecule. Here we broaden this study by including an extensive set of H2O and 28SiO lines. It is the first time such a comprehensive study is performed for this source. The oxygen isotopic ratios and the 28SiO abundance profile can be connected to the initial stellar mass and to crucial aspects of dust formation at the base of the stellar wind, respectively.Methods. We model the molecular emission observed by the three instruments on board Herschel Space Observatory using a state-of-the-art molecular excitation and radiative transfer code. We also account for the dust radiation field in our calculations.Results. We find an H2O ortho-to-para ratio of 2.5 +2.5-1.0, consistent with what is expected for an AGB wind. The O16/O17 ratio indicates that W Hya has an initial mass of about 1.5 M⊙. Although the ortho- and para-H2O lines observed by HIFI appear to trace gas of slightly different physical properties, we find that a turbulence velocity of 0.7 ± 0.1 km s-1 fits the HIFI lines of both spin isomers and those of 28SiO well.Conclusions. The modelling of H2O and 28SiO confirms the properties of the envelope model of W Hya, as derived from 12CO lines, and allows us to constrain the turbulence velocity. The ortho- and para-H216O and 28SiO abundances relative to H2 are (6+ 3-2)×10-4, (3+ 2-1)×10-4, and (3.3 ± 0.8) × 10-5, respectively, in agreement with expectations for oxygen-rich AGB outflows. Assuming a solar silicon-to-carbon ratio, the 28SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles.

The processes leading to dust formation and the subsequent role it plays in driving mass-loss in cool evolved stars is an area of intense study. Here, we present high resolution ALMA Science Verification data of the continuum emission around the highly evolved oxygen-rich red supergiant VY CMa. These data enable us to study the dust in its inner circumstellar environment at a spatial resolution of 129 mas at 321 GHz and 59 mas at 658 GHz, allowing us to trace dust on spatial scales down to 11 R$_{\star}$ (71 AU). Two prominent dust components are detected and resolved. The brightest dust component, C, is located 334 mas (61 R$_{\star}$) south-east of the star and has a dust mass of at least $2.5\times 10^{-4} $M$_{\odot}$. It has an emissivity spectral index of $\beta =-0.1$ at its peak, implying that it is either optically thick at these frequencies with a cool core of $T_{d}\lesssim 100$ K, and/or contains very large dust grains. Interestingly, not a single molecule in the ALMA data has emission close to the peak of this massive dust clump. The other main dust component, VY, is located at the position of the star and contains a total dust mass of $4.0 \times 10^{-5} $M$_{\odot}$. It also contains a weaker dust feature extending over $60 $R$_{\star}$ to the north with the total component having a typical emissivity spectral index of $\beta =0.7$. We find that $>17%$ of the dust mass around VY CMa is located in clumps ejected within a more quiescent roughly spherical stellar wind, with a quiescent dust mass loss rate of $5 \times 10^{-6}$ M$_{\odot} $yr$^{-1}$. The observations suggest a continuous preferentially directed mass-loss from the star over many decades and do not support current models of convective driven only mass loss in red supergiant stars. We thus suggest other forces, i.e., MHD disturbances, are also needed to explain the observations. ispartof: Astronomy & Astrophysics vol:573 pages:5- status: published