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.

Stars with initial masses between $\sim0.8$ and 8~$M_\odot$ present copious mass loss during the asymptotic giant branch (AGB) at the end of their lives. The accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the extended AGB stellar atmospheres. The details of this process are not yet well understood, however. Using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL, we observed light polarised by grains around W\,Hya, SW\,Vir, and R\,Crt, which have mass-loss rates between 10$^{-7}$ and 10$^{-6}~M_\odot~{\rm yr^{-1}}$. We find the distribution of dust to be asymmetric around 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 directly inferred from the ZIMPOL data. The dust grains polarise light more efficiently at 0.65~$\mu$m for R\,Crt and SW\,Vir and at 0.82~$\mu$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~\mu$m, while those around W\,Hya were larger, with sizes $\gtrsim 0.1~\mu$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. The radial profile of the polarised light observed around W\,Hya reveal a steep dust density profile. We find the wind-acceleration region of W\,Hya to extend to at least $\sim 7~R_\star$, in agreement with theoretical predictions of acceleration up to $\sim 12~R_\star$.

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.

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.

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.

The outflows of oxygen-rich asymptotic giant branch (AGB) stars are thought to be driven by radiation pressure by photon scattering on grains with sizes of tenths of microns. The details of the formation of dust in the extended atmospheres of these stars and the mass-loss process is still not well understood. We obtained quasi-simultaneous observations of the AGB star Mira using ALMA and ZIMPOL to probe the distribution of gas and large dust grains, respectively. The polarized light images show dust around Mira~A, the companion (Mira~B) and in a trail that connects the two sources. ALMA reveals that dust around Mira~A is contained in a high-gas-density region with a significant fraction of the polarized light arising from its edge. We constrained the gas density, temperature, and velocity within a few stellar radii from the star by modelling the CO $v=1, J=3-2$ line. We find a mass $(\sim 3.8 \pm 1.3) \times 10^{-4}~M_\odot$ to be contained between the stellar millimetre photosphere, $R^{\rm 338~GHz}_\star$, and $4~R^{\rm 338~GHz}_\star$. Our best-fit models with lower masses also reproduce the $^{13}$CO $v=0, J=3-2$ line emission from this region. We find TiO$_2$ and AlO abundances corresponding to 4.5% and $< 0.1$% of the total titanium and aluminium expected for a solar-composition gas. The low abundance of AlO allows for efficient Al depletion into dust already very close to the star, as expected from thermal dust emission observations and theoretical calculations of Mira variables. We constrain the presence of aluminium oxide grains based on the scattered light observations and our gas-phase model. We find that aluminium oxide grains can account for a significant fraction of the total aluminium atoms in this region only if the grains have sizes $\lesssim 0.02~\mu$m. This is an order of magnitude smaller than the maximum sizes predicted by dust-formation and wind-driving models.

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.

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.

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.

(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.

(abbreviated) We aim to constrain the dust formation histories in the winds of oxygen-rich AGB stars. We have obtained ALMA observations with a spatial resolution of 120x150 mas tracing the dust formation region of a low mass-loss rate and a high mass-loss rate AGB star, R Dor and IK Tau. Emission line profiles of AlO, AlOH and AlCl are detected and are used to derive a lower limit of atomic aluminium incorporated in molecules. 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 3 molecules is small, <=1.1e-7, for both stars, i.e. only <=2% of the total aluminium budget. This leaves ample of room for aluminium to be incorporated in grains. A fundamental result is that AlO and AlOH, being the direct precursors of alumina grains, are detected well beyond the onset of the dust condensation proving that the aluminium oxide condensation cycle is not fully efficient. 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 at a distance of ~1.5 Rstar, in particular for the stars 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 polarized light signal. In particular, the estimated grain temperature of Al2O3 is too high for the grains to retain their amorphous structure. We propose that large gas-phase (Al2O3)n-clusters (n>34) can be the potential agents of the broad 11 micron feature in the SED and in the interferometric data and we explain how these large clusters can be formed.

We use Atacama Large Millimeter/submillimeter Array Band 5 science verification observations of the red supergiant VY CMa to study the polarization of SiO thermal/masers lines and dust continuum at ~1.7 mm wavelength. We analyse both linear and circular polarization and derive the magnetic field strength and structure, assuming the polarization 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 polarization. We detect, for the first time, significant polarization (~3%) of the circumstellar dust emission at millimeter wavelengths. The polarization is uniform with an electric vector position angle of $\sim8^\circ$. Varying levels of linear polarization are detected for the J=4-3 28SiO v=0, 1, 2, and 29SiO v=0, 1 lines, with the strongest polarization fraction of ~30% found for the 29SiO v=1 maser. The linear polarization vectors rotate with velocity, consistent with earlier observations. We also find significant (up to ~1%) circular polarization 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. Emission from magnetically aligned grains is the most likely origin of the observed continuum polarization. This implies that the dust is embedded in a magnetic field >13 mG. The maser line polarization 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 polarization is due to Zeeman splitting, it indicates a magnetic field strength of ~1-3 Gauss, consistent with previous maser observations.

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.

Using the PACS and SPIRE spectrometers on-board the Herschel Space Observatory, we obtained spectra of two red supergiants (RSGs) in the Large Magellanic Cloud (LMC). Multiple rotational CO emission lines (J=6-5 to 15-14) and 15 H2O lines were detected from IRAS 05280-6910, and one CO line was detected from WOH G64. This is the first time CO rotational lines have been detected from evolved stars in the LMC. Their CO line intensities are as strong as those of the Galactic RSG, VY CMa. Modelling the CO lines and the spectral energy distribution results in an estimated mass-loss rate for IRAS 05280-6910 of 3x10^-4 Msun per yr. The model assumes a gas-to-dust ratio and a CO-to-H2 abundance ratio is estimated from the Galactic values scaled by the LMC metallicity ([Fe/H]~-0.3), i.e., that the CO-to-dust ratio is constant for Galactic and LMC metallicities within the uncertainties of the model. The key factor determining the CO line intensities and the mass-loss rate found to be the stellar luminosity.

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.

Context. The recent detection of warm H$_2$O 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 H$_2$O vapor formation. In the first, periodic shocks passing through the medium immediately above the stellar surface lead to H$_2$O formation. In the second, penetration of ultraviolet interstellar radiation through a clumpy circumstellar medium leads to the formation of H$_2$O molecules in the intermediate wind. Aims. We aim to determine the properties of H$_2$O emission for a sample of 18 carbon-rich AGB stars and subsequently constrain which of the above mechanisms provides the most likely warm H$_2$O formation pathway. Methods, Results, and Conclusions. See paper.

Context: AGB stars lose a large percentage of their mass in a dust-driven wind. This creates a circumstellar envelope, which can be studied through thermal dust emission and molecular emission lines. In the case of high mass-loss rates, this study is complicated by the high optical depths and the intricate coupling between gas and dust radiative transfer characteristics. An important aspect of the physics of gas-dust interactions is the strong influence of dust on the excitation of several molecules, including H2O. Aims: The dust and gas content of the envelope surrounding the high mass-loss rate OH/IR star OH 127.8+0.0, as traced by Herschel observations, is studied, with a focus on the H2O content and the dust-to-gas ratio. We report detecting a large number of H2O vapor emission lines up to J = 9 in the Herschel data, for which we present the measured line strengths. Methods: The treatments of both gas and dust species are combined using two numerical radiative transfer codes. The method is illustrated for both low and high mass-loss-rate sources. Specifically, we discuss different ways of assessing the dust-to-gas ratio: 1) from the dust thermal emission spectrum and the CO molecular gas line strengths; 2) from the momentum transfer from dust to gas and the measured gas terminal velocity; and 3) from the determination of the required amount of dust to reproduce H2O lines for a given H2O vapor abundance. These three diagnostics probe different zones of the outflow, for the first time allowing an investigation of a possible radial dependence of the dust-to-gas ratio. Results: ... Continued in article.

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.

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, thus 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 a dust emissivity spectral index of $\beta =-0.1$ at its peak, implying that it is optically thick at these frequencies with a cool core of $T_{d}\lesssim 100$ K. 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 dust emissivity spectral index of $\beta =0.7$. We find that at least $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 anisotropic morphology of the dust indicates a continuous, directed mass loss over a few decades, suggesting that this mass loss cannot be driven by large convection cells alone.

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. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the $^{12}$CO molecule. Here we broaden this study by modelling an extensive set of H$_{2}$O and $^{28}$SiO lines observed by the three instruments on board Herschel using a state-of-the-art molecular excitation and radiative transfer code. The oxygen isotopic ratios and the $^{28}$SiO 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. The modelling of H$_{2}$O and $^{28}$SiO confirms the properties of the envelope model of W Hya derived from $^{12}$CO lines. We find an H$_2$O ortho-to-para ratio of 2.5\,$^{+2.5}_{-1.0}$, consistent with what is expected for an AGB wind. The O$^{16}$/O$^{17}$ ratio indicates that W Hya has an initial mass of about 1.5 M$_\odot$. Although the ortho- and para-H$_{2}$O lines observed by HIFI appear to trace gas of slightly different physical properties, a turbulence velocity of $0.7\pm0.1$ km s$^{-1}$ fits the HIFI lines of both spin isomers and those of $^{28}$SiO well. The ortho- and para-H$_2^{16}$O and $^{28}$SiO abundances relative to H$_{2}$ are $(6^{+3}_{-2}) \times 10^{-4}$, $(3^{+2}_{-1}) \times 10^{-4}$, and $(3.3\pm 0.8)\times 10^{-5}$, respectively. Assuming a solar silicon-to-carbon ratio, the $^{28}$SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles.

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.

Spectra, taken with the heterodyne instrument, HIFI, aboard the Herschel Space Observatory, of O-rich asymptotic giant branch (AGB) stars which form part of the guaranteed time key program HIFISTARS are presented. The aim of this program is to study the dynamical structure, mass-loss driving mechanism, and chemistry of the outflows from AGB stars as a function of chemical composition and initial mass. We used the HIFI instrument to observe nine AGB stars, mainly in the H2O and high rotational CO lines We investigate the correlation between line luminosity, line ratio and mass-loss rate, line width and excitation energy. A total of nine different molecules, along with some of their isotopologues have been identified, covering a wide range of excitation temperature. Maser emission is detected in both the ortho- and para-H2O molecules. The line luminosities of ground state lines of ortho- and para-H2O, the high-J CO and NH3 lines show a clear correlation with mass-loss rate. The line ratios of H2O and NH3 relative to CO J=6-5 correlate with the mass-loss rate while ratios of higher CO lines to the 6-5 is independent of it. In most cases, the expansion velocity derived from the observed line width of highly excited transitions formed relatively close to the stellar photosphere is lower than that of lower excitation transitions, formed farther out, pointing to an accelerated outflow. In some objects, the vibrationally excited H2O and SiO which probe the acceleration zone suggests the wind reaches its terminal velocity already in the innermost part of the envelope, i.e., the acceleration is rapid. Interestingly, for R Dor we find indications of a deceleration of the outflow in the region where the material has already escaped from the star.

We present the first results of a high-spectral-resolution survey of the carbon-rich evolved star IRC+10216 that was carried out with the HIFI spectrometer onboard Herschel. This survey covers all HIFI bands, with a spectral range from 488 to 1901GHz. In this letter we focus on the band-1b spectrum, in a spectral range 554.5-636.5GHz, where we identified 130 spectral features with intensities above 0.03 K and a signal-to-noise ratio >5. Detected lines arise from HCN, SiO, SiS, CS, CO, metal-bearing species and, surprisingly, silicon dicarbide (SiC2). We identified 55 SiC2 transitions involving energy levels between 300 and 900 K. By analysing these rotational lines, we conclude that SiC2 is produced in the inner dust formation zone, with an abundance of ~2x10^-7 relative to molecular hydrogen. These SiC2 lines have been observed for the first time in space and have been used to derive an SiC2 rotational temperature of ~204 K and a source-averaged column density of ~6.4x10^15 cm^-2. Furthermore, the high quality of the HIFI data set was used to improve the spectroscopic rotational constants of SiC2.