W Aql is an asymptotic giant branch (AGB) star with an atmospheric elemental abundance ratio C/O$\approx$0.98 and reported circumstellar molecular abundances intermediate between those of M-type (C/O$<$1) and C-type (C/O$>$1) AGB stars. This intermediate status is considered typical for S-type stars, although our understanding of the chemical content of their circumstellar envelopes (CSEs) is currently rather limited. We performed observations in the frequency range 159-268 GHz with the APEX telescope and make abundance estimates through comparison to available spectra towards some well-studied AGB stars and based on rotational diagram analysis in the case of SiC2. We conclude that W Aql's CSE appears considerably closer to that of a C-type AGB star than to that of an M-type AGB star. In particular, we detect emission from C2H, SiC2, SiN, and HC3N, molecules previously only detected towards the CSEs of C-type stars. This conclusion, based on the chemistry of the gaseous component of the CSE, is further supported by reports in the literature on the presence of atmospheric molecular bands and spectral features of dust species typical for C-type AGB stars. Although our observations mainly trace species in the outer regions of the CSE, our conclusion matches closely that based on recent chemical equilibrium models for the inner wind of S-type stars: the atmospheric and circumstellar chemistry of S-type stars likely resembles that of C-type AGB stars much more closely than that of M-type AGB stars. Further observational investigation of the gaseous circumstellar chemistry of S-type stars is required to characterise its dependence on the atmospheric C/O. Non-equilibrium chemical models of the CSEs of AGB stars need to address the particular class of S-type stars and the chemical variety that is induced by the range in atmospheric C/O.

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 Origins Space Telescope (Origins) traces our cosmic history, from the formation of the first galaxies and the rise of metals to the development of habitable worlds and present-day life. Origins does this through exquisite sensitivity to infrared radiation from ions, atoms, molecules, dust, water vapor and ice, and observations of extra-solar planetary atmospheres, protoplanetary disks, and large-area extragalactic fields. Origins operates in the wavelength range 2.8 to 588 microns and is 1000 times more sensitive than its predecessors due to its large, cold (4.5 K) telescope and advanced instruments. Origins was one of four large missions studied by the community with support from NASA and industry in preparation for the 2020 Decadal Survey in Astrophysics. This is the final study report.

The AGB star $\pi^{1}$ Gruis has a known companion (at a separation of ~400 AU) which cannot explain the strong deviations from the spherical symmetry of the CSE. Recently, hydrodynamic simulations of mass transfer in closer binary systems have successfully reproduced the spiral-shaped CSEs found around a handful of sources. There is growing evidence for an even closer, undetected companion complicating the case of $\pi^{1}$ Gruis further. The improved spatial resolution allows for the investigation of the complex circumstellar morphology and the search for imprints on the CSE of the third component. We have observed the 12CO J=3-2 line emission from $\pi^{1}$ Gruis using both the compact and extended array of Atacama Large Millimeter/submillimeter Array (ALMA). The interferometric data has furthermore been combined with data from the ALMA total power (TP) array. The imaged brightness distribution has been used to constrain a non-local, non-LTE 3D radiative transfer model of the CSE. The high-angular resolution ALMA data have revealed the first example of a source on the AGB where both a faster bipolar outflow and a spiral pattern along the orbital plane can be seen in the gas envelope. The spiral can be traced in the low- to intermediate velocity, from 13 to 25 km s$^{-1}$, equatorial torus. The largest spiral-arm separation is $\approx$5".5 and consistent with a companion with an orbital period of $\approx$330 yrs and a separation of less than 70 AU. The kinematics of the bipolar outflow is consistent with it being created during a mass-loss eruption where the mass-loss rate from the system increased by at least a factor of 5 during 10-15 yrs. The spiral pattern is the result of an undetected companion. The bipolar outflow is the result of a rather recent mass-loss eruption event.

Ultraviolet (UV) photodissociation of CO controls the abundances and distribution of CO and its photodissociation products. This significantly influences the gas-phase chemistry in the circumstellar envelope (CSE) around evolved stars. A better understanding of CO photodissociation in outflows also provides a more precise estimate of mass-loss rates. We aim to update the CO photodissociation rate in an expanding spherical envelope assuming that the interstellar radiation field (ISRF) photons penetrate through the envelope. This will allow us to precisely estimate the CO abundance distributions in CSEs around evolved stars. We used the most recent CO spectroscopic data to precisely calculate the depth dependency of the photodissociation rate of each CO dissociating line. We calculated the CO self- and mutual-shielding functions in an expanding envelope. We investigated the dependence of the CO profile on the five fundamental parameters mass-loss rate, the expansion velocity, the CO initial abundance, the CO excitation temperature, and the strength of the ISRF. Our derived CO envelope size is smaller than the commonly used radius derived by Mamon et al. 1988. The difference between results varies from 1% to 39% and depends on the H2 and CO densities of the envelope. We list two fitting parameters for a large grid of models to estimate the CO abundance distribution. We demonstrate that the CO envelope size can differ between outflows with the same effective content of CO, but different CO abundance, mass-loss rate, and the expansion velocity as a consequence of differing amounts of shielding by H2 and CO. We find that the abundance of CO close to the star and the outflow density both can have a significant effect on the size of the molecular envelope. We also demonstrate that modest variations in the ISRF can cause measurable differences in the envelope extent.

Models of the chemical evolution of the interstellar medium, galaxies, and the Universe rely on our understanding of the amounts and chemical composition of the material returned by stars and supernovae. Stellar yields are obtained from stellar-evolution models, which currently lack predictive prescriptions of stellar mass loss, although it significantly affects stellar lifetimes, nucleosynthesis, and chemical ejecta. Galaxy properties are derived from observations of the integrated light of bright member stars. Stars in the late stages of their evolution are among the infrared-brightest objects in galaxies. An unrealistic treatment of the mass-loss process introduces significant uncertainties in galaxy properties derived from their integrated light. We describe current efforts and future needs and opportunities to characterize AGB outflows: driving mechanisms, outflow rates, underlying fundamental physical and chemical processes such as dust grain formation, and dependency of these on metallicity.

Evolved stars dominate galactic spectra, enrich the galactic medium, expand to change their planetary systems, eject winds of a complex nature, produce spectacular nebulae and illuminate them, and transfer material between binary companions. While doing this, they fill the HR diagram with diagnostic loops that write the story of late stellar evolution. Evolved stars sometimes release unfathomable amounts of energy in neutrinos, light, kinetic flow, and gravitational waves. During these late-life times, stars evolve complexly, with expansion, convection, mixing, pulsation, mass loss. Some processes have virtually no spatial symmetries, and are poorly addressed with low-resolution measurements and analysis. Even a "simple" question as how to model mass loss resists solution. However, new methods offer increasingly diagnostic tools. Astrometry reveals populations and groupings. Pulsations/oscillations support study of stellar interiors. Optical/radio interferometry enable 2-3d imagery of atmospheres and shells. Bright stars with rich molecular spectra and velocity fields are a ripe opportunity for imaging with high spatial and spectral resolution, giving insight into the physics and modeling of later stellar evolution.

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.

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). 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. We observed the carbon AGB stars U Ant, DR Ser, and V644 Sco at 870 micron 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 tested 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. We derive dust masses in the shells of a few 10e-5 Msun, assuming spherical, solid grains. The best-fit grain radii are comparatively large, and indicate the presence of grains between 0.1 micron-2 micron. 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.

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.

Eleven nearby (<300 pc), short-period (50-130 days) asymptotic giant branch (AGB) stars were observed in the CO J = (2-1) line. Detections were made towards objects that have evidence for dust production (Ks-[22] >~ 0.55 mag; AK Hya, V744 Cen, RU Crt, alpha Her). Stars below this limit were not detected (BQ Gem, eps Oct, NU Pav, II Hya, CL Hyi, ET Vir, SX Pav). Ks-[22] colour is found to trace mass-loss rate to well within an order of magnitude. This confirms existing results, indicating a factor of 100 increase in AGB-star mass-loss rates at a pulsation period of ~60 days, similar to the known "superwind" trigger at ~300 days. Between ~60 and ~300 days, an approximately constant mass-loss rate and wind velocity of ~3.7 x 10^-7 solar masses per year and ~8 km/s is found. While this has not been corrected for observational biases, this rapid increase in mass-loss rate suggests a need to recalibrate the treatment of AGB mass loss in stellar evolution models. The comparative lack of correlation between mass-loss rate and luminosity (for L <~ 6300 solar luminosities) suggests that the mass-loss rates of low-luminosity AGB-star winds are set predominantly by pulsations, not radiation pressure on dust, which sets only the outflow velocity. We predict that mass-loss rates from low-luminosity AGB stars, which exhibit optically thin winds, should be largely independent of metallicity, but may be strongly dependent on stellar mass.

The Origins Space Telescope, one of four large Mission Concept studies sponsored by NASA for review in the 2020 US Astrophysics Decadal Survey, will open unprecedented discovery space in the infrared, unveiling our cosmic origins. We briefly describe in this article the key science themes and architecture for OST. With a sensitivity gain of up to a factor of 1,000 over any previous or planned mission, OST will open unprecedented discovery space, allow us to peer through an infrared window teeming with possibility. OST will fundamentally change our understanding of our cosmic origins - from the growth of galaxies and black holes, to uncovering the trail of water, to life signs in nearby Earth-size planets, and discoveries never imagined. Built to be highly adaptable, while addressing key science across many areas of astrophysics, OST will usher in a new era of infrared astronomy.

We surveyed 20 AGB stars of different chemical types using the APEX telescope, and combined this with an IRAM 30 m and APEX survey of CS and SiS emission towards over 30 S-type stars. For those stars with detections, we performed radiative transfer modelling to determine abundances and abundance distributions. We detect CS towards all the surveyed carbon stars, some S-type stars, and the highest mass-loss rate oxygen-rich stars ($> 5\times 10^{-6}$ Msol yr$^{-1}$). SiS is detected towards the highest mass-loss rate sources of all chemical types ($> 8\times 10^{-7}$ Msol yr$^{-1}$). We find CS peak fractional abundances ranging from ~ $ 4\times 10^{-7}$ to ~ $2\times 10^{-5}$ for the carbon stars, from ~ $ 3\times 10^{-8}$ to ~ $1\times 10^{-7}$ for the oxygen-rich stars and from ~ $ 1\times 10^{-7}$ to ~ $8\times 10^{-6}$ for the S-type stars. We find SiS peak fractional abundances ranging from ~ $ 9\times 10^{-6}$ to ~ $ 2\times 10^{-5}$ for the carbon stars, from ~ $ 5\times 10^{-7}$ to ~ $ 2\times 10^{-6}$ for the oxygen-rich stars, and from ~ $ 2\times 10^{-7}$ to ~ $ 2\times 10^{-6}$ for the S-type stars. We derived Si$^{32}$S/Si$^{34}$S = 11.4 for AI Vol, the only star for which we had a reliable isotopologue detection. Overall, we find that wind density plays an important role in determining the chemical composition of AGB CSEs. It is seen that for oxygen-rich AGB stars both CS and SiS are detected only in the highest density circumstellar envelopes and their abundances are generally lower than for carbon-rich AGB stars by around an order of magnitude. For carbon-rich and S-type stars SiS was also only detected in the highest density circumstellar envelopes, while CS was detected consistently in all surveyed carbon stars and sporadically among the S-type stars.

Context. With the outstanding spatial resolution and sensitivity of the Atacama Large Millimeter/sub-millimeter Array (ALMA), molecular gas other than the abundant CO can be observed and resolved in circumstellar envelopes (CSEs) around evolved stars, such as the binary S-type Asymptotic Giant Branch (AGB) star W Aquilae. Aims. We aim to constrain the chemical composition of the CSE and determine the radial abundance distribution, the photospheric peak abundance, and isotopic ratios of a selection of chemically important molecular species in the innermost CSE of W Aql. The derived parameters are put into the context of the chemical evolution of AGB stars and are compared with theoretical models. Methods. We employ one-dimensional radiative transfer modeling - with the accelerated lambda iteration (ALI) radiative transfer code - of the radial abundance distribution of a total of five molecular species (CS, SiS, 30SiS, 29SiO and H13CN) and determine the best fitting model parameters based on high-resolution ALMA observations as well as archival single-dish observations. The additional advantage of the spatially resolved ALMA observations is that we can directly constrain the radial profile of the observed line transitions from the observations. Results. We derive abundances and e-folding radii for CS, SiS, 30SiS, 29SiO and H13CN and compare them to previous studies, which are based only on unresolved single-dish spectra. Our results are in line with previous results and are more accurate due to resolution of the emission regions.

We present the detection of neutral atomic carbon CI(3 P1 - 3 P0 ) line emission towards omi Cet. This is the first time that CI is detected in the envelope around an oxygen-rich M-type asymptotic giant branch (AGB) star. We also confirm the previously tentative CI detection around V Hya, a carbon-rich AGB star. As one of the main photodissociation products of parent species in the circumstellar envelope (CSE) around evolved stars, CI can be used to trace sources of ultraviolet (UV) radiation in CSEs. The observed flux density towards omi Cet can be reproduced by a shell with a peak atomic fractional abundance of $2.4 \times 10^{-5}$ predicted based on a simple chemical model where CO is dissociated by the interstellar radiation field. However, the CI emission is shifted by $\sim$ 4 km/s from the stellar velocity. Based on this velocity shift, we suggest that the detected CI emission towards omi Cet potentially arises from a compact region near its hot binary companion. The velocity shift could, therefore, be the result of the orbital velocity of the binary companion around omi Cet. In this case, the CI column density is estimated to be $1.1 \times 10^{19}$ cm$^{-2}$. This would imply that strong UV radiation from the companion and/or accretion of matter between two stars is most likely the origin of the CI enhancement. However, this hypothesis can be confirmed by high-angular resolution observations.

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 insights into the circumstellar chemistry of asymptotic giant branch (AGB) stars are largely based on studies of carbon-rich stars and stars with high mass-loss rates. In order to expand the current molecular inventory of evolved stars we present a spectral scan of the nearby, oxygen-rich star R Dor, a star with a low mass-loss rate ($\sim2\times10^{-7}M_{\odot}$/yr). We carried out a spectral scan in the frequency ranges 159.0-321.5GHz and 338.5-368.5 GHz (wavelength range 0.8-1.9mm) using the SEPIA/Band-5 and SHeFI instruments on the APEX telescope and we compare it to previous surveys, including one of the oxygen-rich AGB star IK Tau, which has a high mass-loss rate ($\sim5\times10^{-6}M_{\odot}$/yr). The spectrum of R Dor is dominated by emission lines of SO$_2$ and the different isotopologues of SiO. We also detect CO, H$_2$O, HCN, CN, PO, PN, SO, and tentatively TiO$_2$, AlO, and NaCl. Sixteen out of approximately 320 spectral features remain unidentified. Among these is a strong but previously unknown maser at 354.2 GHz, which we suggest could pertain to H$_2$SiO, silanone. With the exception of one, none of these unidentified lines are found in a similarly sensitive survey of IK Tau performed with the IRAM 30m telescope. We present radiative transfer models for five isotopologues of SiO ($^{28}$SiO, $^{29}$SiO, $^{30}$SiO, Si$^{17}$O, Si$^{18}$O), providing constraints on their fractional abundance and radial extent. We derive isotopic ratios for C, O, Si, and S and estimate that R Dor likely had an initial mass in the range 1.3-1.6$M_{\odot}$, in agreement with earlier findings based on models of H$_2$O line emission. From the presence of spectral features recurring in many of the measured thermal and maser emission lines we tentatively identify up to five kinematical components in the outflow of R Dor, indicating deviations from a smooth, spherical wind.

Context: We describe the new SEPIA (Swedish-ESO PI Instrument for APEX) receiver, which was designed and built by the Group for Advanced Receiver Development (GARD), at Onsala Space Observatory (OSO) in collaboration with ESO. It was installed and commissioned at the APEX telescope during 2015 with an ALMA Band 5 receiver channel and updated with a new frequency channel (ALMA Band 9) in February 2016. Aims: This manuscript aims to provide, for observers who use the SEPIA receiver, a reference in terms of the hardware description, optics and performance as well as the commissioning results. Methods: Out of three available receiver cartridge positions in SEPIA, the two current frequency channels, corresponding to ALMA Band 5, the RF band 158--211 GHz, and Band 9, the RF band 600--722 GHz, provide state-of-the-art dual polarization receivers. The Band 5 frequency channel uses 2SB SIS mixers with an average SSB noise temperature around 45K with IF (intermediate frequency) band 4--8 GHz for each sideband providing total 4x4 GHz IF band. The Band 9 frequency channel uses DSB SIS mixers with a noise temperature of 75--125K with IF band 4--12 GHz for each polarization. Results: Both current SEPIA receiver channels are available to all APEX observers.

For the carbon AGB star R Sculptoris, the uncertain distance significantly affects the interpretation of observations regarding the evolution of the stellar mass loss during and after the most recent thermal pulse. We aim to provide a new, independent measurement of the distance to R Sculptoris, reducing the absolute uncertainty of the distance estimate to this source. R Scl is a semi-regular pulsating star, surrounded by a thin shell of dust and gas created during a thermal pulse approximately 2000 years ago. The stellar light is scattered by the dust particles in the shell at a radius of 19 arcsec. The variation in the stellar light affects the amount of dust-scattered light with the same period and amplitude ratio, but with a phase lag that depends on the absolute size of the shell. We measured this phase lag by observing the star R Scl and the dust-scattered stellar light from the shell at five epochs between June - December 2017. By observing in polarised light, we imaged the shell in the plane of the sky, removing any uncertainty due to geometrical effects. The phase lag gives the absolute size of the shell, and together with the angular size of the shell directly gives the absolute distance to R Sculptoris. We measured a phase lag between the stellar variations and the variation in the shell of 40.0 +/- 4.0 days. The angular size of the shell is measured to be 19.1 arcsec +/- 0.7 arcsec. Combined, this gives an absolute distance to R Sculptoris of 361 +/- 44 pc. We independently determined the absolute distance to R Scl with an uncertainty of 12%. The estimated distance is consistent with previous estimates, but is one of the most accurate distances to the source to date. In the future, using the variations in polarised, dust-scattered stellar light, may offer an independent possibility to measure reliable distances to AGB stars.

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.

The S-type asymptotic giant branch (AGB) star $\pi^{1}$ Gruis has a known companion at a separation of $\approx$400 AU. The envelope structure, including an equatorial torus and a fast bipolar outflow, is rarely seen in the AGB phase and is particularly unexpected in such a wide binary system. Therefore a second, closer companion has been suggested, but the evidence is not conclusive. The new ALMA $^{12}$CO and $^{13}$CO $J$=3-2 data, together with previously published $^{12}$CO $J$=2-1 data from the Submillimeter Array (SMA), and the $^{12}$CO $J$=5-4 and $J$=9-8 lines observed with Herschel/Heterodyne Instrument for the Far-Infrared (HIFI), is modeled with the 3D non-LTE radiative transfer code SHAPEMOL. The data analysis clearly confirms the torus-bipolar structure. The 3D model of the CSE that satisfactorily reproduces the data consists of three kinematic components: a radially expanding torus with velocity slowly increasing from 8 to 13 km s$^{-1}$ along the equator plane; a radially expanding component at the center with a constant velocity of 14 km s$^{-1}$; and a fast, bipolar outflow with velocity proportionally increasing from 14 km s$^{-1}$ at the base up to 100 km s$^{-1}$ at the tip, following a linear radial dependence. The results are used to estimate an average mass-loss rate during the creation of the torus of 7.7$\times$10$^{-7}$ M$_{\odot}$ yr$^{-1}$. The total mass and linear momentum of the fast outflow are estimated at 7.3$\times$10$^{-4}$ M$_{\odot}$ and 9.6$\times$10$^{37}$ g cm s$^{-1}$, respectively. The momentum of the outflow is in excess (by a factor of about 20) of what could be generated by radiation pressure alone, in agreement with recent findings for more evolved sources. The best-fit model also suggests a $^{12}$CO/$^{13}$CO abundance ratio of 50. Possible shaping scenarios for the gas envelope are discussed

The CO(J=3-2) emission from the CSE of the binary S-type AGB star W Aql has been observed at subarcsecond resolution using ALMA. The aim of this paper is to investigate the wind properties of the AGB star and to analyse how the known companion has shaped the CSE. The average mass-loss rate during the creation of the detected CSE is estimated through modelling, using the ALMA brightness distribution and previously published single-dish measurements as observational constraints. The ALMA observations are presented and compared to the results from a 3D smoothed particle hydrodynamics (SPH) binary interaction model with the same properties as the W Aql system and with two different orbital eccentricities. Three-dimensional radiative transfer modelling is performed and the response of the interferometer is modelled and discussed. The estimated average mass-loss rate of W~Aql agrees with previous results. The size of the emitting region is consistent with photodissociation models. The CO(J=3-2) emission is dominated by a smooth component overlayed with two weak arc patterns with different separations. The larger pattern is predicted by the binary interaction model with separations of 10" and therefore likely due to the known companion. It is consistent with a binary orbit with low eccentricity. The smaller separation pattern is asymmetric and coincides with the dust distribution, but the separation timescale (200 yrs) is not consistent with any known process of the system. The separation of the known companions of the system is large enough to not have a very strong effect on the circumstellar morphology. The density contrast across the envelope of a binary with an even larger separation will not be easily detectable, even with ALMA, unless the orbit is strongly asymmetric or the AGB star has a much larger mass-loss rate.

Aims. Our goal is to characterize the intermediate age, detached shell carbon star U Antliae morphologically and physically in order to study the mass-loss evolution after a possible thermal pulse. Methods. High spatial resolution ALMA observations of unprecedented quality in thermal CO lines allow us to derive first critical spatial and temporal scales and constrain modeling efforts to estimate mass-loss rates for both the present day as well as the ejection period of the detached shell. Results. The detached shell is remarkably thin, overall spherically symmetric, and shows a barely resolved filamentary substructure possibly caused by instabilities in the interaction zone of winds with different outflow velocities. The expansion age of the detached shell is of the order of 2700 years and its overall width indicates a high expansion-velocity and high mass-loss period of only a few hundred years at an average mass-loss rate of $\approx$ 10$^{-5}$ $M_\odot$ yr$^{-1}$. The post-high-mass-loss-rate-epoch evolution of U Ant shows a significant decline to a substantially lower gas expansion velocity and a mass-loss rate amounting to 4$\times$10$^{-8}$ $M_\odot$ yr$^{-1}$, at present being consistent with evolutionary changes as predicted for the period between thermal pulses.

Through a survey of (sub-)millimetre emission lines of various sulphur-bearing molecules, we aim to determine which molecules are the primary carriers of sulphur in different types of AGB stars. In this paper, the first in a series, we investigate the occurrence of H$_2$S in AGB circumstellar envelopes and determine its abundance, where possible. We have surveyed 20 AGB stars with a range of mass-loss rates and of different chemical types using the APEX telescope to search for rotational transition lines of five key sulphur-bearing molecules: CS, SiS, SO, SO$_2$ and H$_2$S. Here we present our results for H$_2$S, including detections, non-detections and detailed radiative transfer modelling of the detected lines. We compare results based on different descriptions of the molecular excitation of H$_2$S and different abundance distributions, including those derived from chemical modelling results. We detected H$_2$S towards five AGB stars, all of which have high mass-loss rates of $\dot{M}\geq 5\times 10^{-6}M_\odot$ yr$^{-1}$ and are oxygen-rich. H$_2$S was not detected towards the carbon or S-type stars that fall in a similar mass-loss range. For the stars in our sample with detections, we find peak o-H$_2$S abundances relative to H$_2$ between $4\times 10^{-7}$ and $2.5\times 10^{-5}$. Overall, we conclude that H$_2$S can play a significant role in oxygen-rich AGB stars with higher mass-loss rates, but is unlikely to play a key role in stars of other chemical types or the lower mass-loss rate oxygen-rich stars. For two sources, V1300 Aql and GX Mon, H$_2$S is most likely the dominant sulphur-bearing molecule in the circumstellar envelope.

(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 investigate the use of 183 GHz H2O masers for characterization of the physical conditions and mass loss process in the circumstellar envelopes of evolved stars. We used APEX SEPIA Band 5 to observe the 183 GHz H2O line towards 2 Red Supergiant and 3 Asymptotic Giant Branch stars. Simultaneously, we observed lines in 28SiO v0, 1, 2 and 3, and for 29SiO v0 and 1. We detected the 183 GHz H2O line towards all the stars with peak flux densities greater than 100 Jy, including a new detection from VY CMa. Towards all 5 targets, the water line had indications of being due to maser emission and had higher peak flux densities than for the SiO lines. The SiO lines appear to originate from both thermal and maser processes. Comparison with simulations and models indicate that 183 GHz maser emission is likely to extend to greater radii in the circumstellar envelopes than SiO maser emission and to similar or greater radii than water masers at 22, 321 and 325 GHz. We speculate that a prominent blue-shifted feature in the W Hya 183 GHz spectrum is amplifying the stellar continuum, and is located at a similar distance from the star as mainline OH maser emission. From a comparison of the individual polarizations, we find that the SiO maser linear polarization fraction of several features exceeds the maximum fraction allowed under standard maser assumptions and requires strong anisotropic pumping of the maser transition and strongly saturated maser emission. The low polarization fraction of the H2O maser however, fits with the expectation for a non-saturated maser. 183 GHz H2O masers can provide strong probes of the mass loss process of evolved stars. Higher angular resolution observations of this line using ALMA Band 5 will enable detailed investigation of the emission location in circumstellar envelopes and can also provide information on magnetic field strength and structure.

The detection and study of molecular gas in born-again stars would be of great importance to understand their composition and chemical evolution. In addition, the molecular emission would be an invaluable tool to explore the physical conditions, kinematics and formation of asymmetric structures in the circumstellar envelopes of these evolved stars. However, until now, all attempts to detect molecular emission from the cool material around born-again stars have failed. We carried out observations using the APEX and IRAM 30m telescopes to search for molecular emission toward four well studied born-again stars, V4334 Sgr, V605 Aql, A30 and A78, that are thought to represent an evolutionary sequence. We detected for the first time emission from HCN and H$^{13}$CN molecules toward V4334 Sgr, and CO emission in V605 Aql. No molecular emission was detected above the noise level toward A30 and A78. A first estimate of the H$^{12}$CN/H$^{13}$CN abundance ratio in the circumstellar environment of V4334 Sgr is $\approx$3, which is similar to the value of the $^{12}$C/$^{13}$C ratio measured from other observations. We derived a rotational temperature of $T_{\rm rot}$=13$\pm1$ K, and a total column density of $N_{{\rm HCN}}$=1.6$\pm0.1\times$10$^{16}$ cm$^{-2}$ for V4334 Sgr. This result sets a lower limit on the amount of hydrogen that was ejected into the wind during the born-again event of this source. For V605 Aql, we obtained a lower limit for the integrated line intensities $I_{^{12}\rm C}$/$I_{^{13}\rm C}$>4.

Abridged. The 12CO/13CO ratio in the circumstellar envelope (CSE) of asymptotic giant branch (AGB) stars has been extensively used as the tracer of the photospheric 12C/13C ratio. However, spatially-resolved ALMA observations of R Scl, a carbon rich AGB star, have shown that the 12CO/13CO ratio is not consistent over the entire CSE. Hence, it can not necessarily be used as a tracer of the 12C/13C ratio. The most likely hypothesis to explain the observed discrepancy between the 12CO/13CO and 12C/13C ratios is CO isotopologue selective photodissociation by UV radiation. Unlike the CO isotopologue ratio, the HCN isotopologue ratio is not affected by UV radiation. Therefore, HCN isotopologue ratios can be used as the tracer of the atomic C ratio in UV irradiated regions. We have performed a detailed non-LTE excitation analysis of circumstellar H12CN and H13CN line emission around R Scl, observed with ALMA and APEX, using a radiative transfer code, ALI. The spatial extent of the molecular distribution for both isotopologues is constrained based on the spatially resolved H13CN(4-3) ALMA observations. We find fractional abundances of H12CN/H2 = (5.0 +\- 2.0) x 10^{-5} and H13CN/H2 = (1.9 +\- 0.4) x 10^{-6} in the inner wind (r < (2.0 +\- 0.25) x 10^{15} cm) of R Scl. The derived circumstellar isotopologue ratio of H12CN/H13CN = 26.3 +\- 11.9 is consistent with the photospheric ratio of 12C/13C ~ 19 \pm 6. We show that the circumstellar H12CN/H13CN ratio traces the photospheric 12C/13C ratio. These results support the previously proposed explanation that CO isotopologue selective-shielding is the main factor responsible for the observed discrepancy between 12C/13C and 12CO/13CO ratios in the inner CSE of R Scl. This indicates that UV radiation impacts on the CO isotopologue ratio.

Aluminium monoxide, AlO, is likely efficiently depleted from the gas around oxygen-rich evolved stars to form alumina clusters and dust seeds. Its presence in the extended atmospheres of evolved stars has been derived from optical spectroscopy. More recently, AlO gas was also detected at long wavelengths around the supergiant VY CMa and the oxygen-rich asymptotic giant branch (AGB) star o Cet (Mira A). In search of AlO, we mined data obtained with APEX, the IRAM 30m telescope, Herschel/HIFI, SMA, and ALMA, which were primarily aimed at studying other molecular species. We report here on observations of AlO towards a sample of eight oxygen-rich AGB stars in different rotational transitions, up to seven for some stars. We present definite detections of one rotational transition of AlO for o Cet and R Aqr, and tentative detections of one transition for R Dor and o Cet, and two for IK Tau and W Hya. The presented spectra of WX Psc, R Cas, and TX Cam show no signature of AlO. For o Cet, R Aqr, and IK Tau, we find that the AlO(N=9-8) emission likely traces the inner parts of the wind, out to only a few tens of AU, where the gas has not yet reached its terminal velocity. The conclusive detections of AlO emission in the case of o Cet and R Aqr confirm the presence of AlO gas in outflows of AGB stars. The tentative detections further support this. Since most of the observations presented in this study were obtained with stronger emission from other species than AlO in mind, observations with higher sensitivity in combination with high angular resolution will improve our understanding of the presence and behaviour of AlO. From the current data sets we cannot firmly conclude whether there is a direct correlation between the wind properties and the detection rate of AlO emission. We hope that this study can serve as a stimulus to perform sample studies in search of AlO in oxygen-rich outflows.

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.

Herschel FIR imaging and spectroscopy were taken at several epochs to probe the central point source and the extended environment of V838 Mon. PACS and SPIRE maps were used to obtain photometry of the near and far dust around V838 Mon. Fitting reveals 0.5-0.6 solar masses of ~19K dust in the environs (~2.7pc) surrounding the star. The surface-integrated infrared flux (signifying the thermal light echo) and derived dust properties do not vary significantly between the epochs. We also fit the SED of the point source. As the peak of the SED lies outside the Herschel spectral range, it is only by incorporating data from other observatories and epochs that we can perform useful fitting; with this we explicitly assume no evolution of the point source between the epochs. We find warm dust with a temperature of ~300K distributed over a radius of 150-200AU. PACS and SPIRE spectra were also used to detect emission lines from the extended environment around the star. We fit the far-infrared lines of CO arising from the point source, from an extended environment around V838 Mon. Assuming a model of a spherical shell for this gas, we find that the CO appears to arise from two temperature zones: a cold zone (Tkin ~18K) that could be associated with the ISM or possibly with a cold layer in the outermost part of the shell, and a warm (Tkin ~400K) zone that is associated with the extended environment of V838 Mon within a region of radius of ~210AU. The SiO lines arise from a warm/hot zone. We did not fit the lines of H2O as they are far more dependent on the model assumed.

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.

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

We aim to constrain the temperature and velocity structures, and H2O abundances in the winds of a sample of M-type AGB stars. We further aim to determine the effect of H2O line cooling on the energy balance in the inner circumstellar envelope. We use two radiative-transfer codes to model molecular emission lines of CO and H2O towards four M-type AGB stars. We focus on spectrally resolved observations of CO and H2O from HIFI. The observations are complemented by ground-based CO observations, and spectrally unresolved CO and H2O observations with PAC. The observed line profiles constrain the velocity structure throughout the circumstellar envelopes (CSEs), while the CO intensities constrain the temperature structure in the CSEs. The H2O observations constrain the o-H2O and p-H2O abundances relative to H2. Finally, the radiative-transfer modelling allows to solve the energy balance in the CSE, in principle including also H2O line cooling. The fits to the line profiles only set moderate constraints on the velocity profile, indicating shallower acceleration profiles in the winds of M-type AGB stars than predicted by dynamical models, while the CO observations effectively constrain the temperature structure. Including H2O line cooling in the energy balance was only possible for the low-mass-loss-rate objects in the sample, and required an ad hoc adjustment of the dust velocity profile in order to counteract extreme cooling in the inner CSE. H2O line cooling was therefore excluded from the models. The constraints set on the temperature profile by the CO lines nevertheless allowed us to derive H2O abundances. The derived H2O abundances confirm previous estimates and are consistent with chemical models. However, the uncertainties in the derived abundances are relatively large, in particular for p-H2O, and consequently the derived o/p-H2O ratios are not well constrained.

The sulphur compounds SO and SO$_2$ have not been widely studied in the circumstellar envelopes of asymptotic giant branch (AGB) stars. By presenting and modelling a large number of SO and SO$_2$ lines in the low mass-loss rate M-type AGB star R Dor, and modelling the available lines of those molecules in a further four M-type AGB stars, we aim to determine their circumstellar abundances and distributions. We use a detailed radiative transfer analysis based on the accelerated lambda iteration method to model circumstellar SO and SO$_2$ line emission and molecular data files for both SO and SO$_2$ that are more extensive than those previously available. Using 17 SO lines and 98 SO2 lines to constrain our models for R Dor, we find an SO abundance of 6.7x10$^{-6}$ and an SO$_2$ abundance of 5x10$^{-6}$ with both species having high abundances close to the star. We also modelled $^{34}$SO and found an abundance of 3.1x10$^{-7}$, giving an $^{32}$SO/$^{34}$SO ratio of 21.6. We derive similar results for the circumstellar SO and SO$_2$ abundances and their distributions for the low mass-loss rate object W Hya. For these stars, the circumstellar SO and SO$_2$ abundances are much higher than predicted by chemical models and these two species may account for all available sulphur. For the higher mass-loss rate stars, we find shell-like SO distributions with peak abundances that decrease and peak abundance radii that increase with increasing mass-loss rate. The positions of the peak SO abundance agree very well with the photodissociation radii of H$_2$O. We find evidence that SO is most likely through the photodissociation of H$_2$O and the subsequent reaction between S and OH. The S-bearing parent molecule appears not to be H$_2$S. The SO$_2$ models suggest an origin close to the star for this species, also disagreeing with current chemical models.

Thermal pulses are fundamental to the chemical evolution of AGB stars and their circumstellar envelopes. A further consequence of thermal pulses is the formation of detached shells of gas and dust around the star. We aim to determine the physical properties of the detached gas shell around R Sculptoris, in particular the shell mass and temperature, and to constrain the evolution of the mass-loss rate during and after a thermal pulse. We analyse CO(1-0), CO(2-1), and CO(3-2) emission, observed by. The spatial resolution of the ALMA data allows us to separate the detached shell emission from the extended emission inside the shell. We perform radiative transfer modelling of both components to determine the shell properties and the post-pulse mass-loss properties. The ALMA data show a gas shell with a radius of 19.5" expanding at 14.3km/s. The different scales probed by the ALMA Cycle 0 array show that the shell must be entirely filled with gas, contrary to the idea of a detached shell. The comparison to single-dish spectra and radiative transfer modelling confirms this. We derive a shell mass of 4.5e-3 Msun with a temperature of 50K. Typical timescales for thermal pulses imply a pulse mass-loss rate of 2.3e-5 Msun/yr. For the post-pulse mass-loss rate, we find evidence for a gradual decline of the mass-loss rate, with an average value of 1.6e-5 Msun/yr. The total amount of mass lost since the last thermal pulse is 0.03 Msun, a factor four higher compared to classical models, with a sharp decline in mass-loss rate immediately after the pulse. We find that the mass-loss rate after a thermal pulse has to decline more slowly than generally expected from models of thermal pulses. This may cause the star to lose significantly more mass during a thermal pulse cycle, which affects the chemical evolution of the star and the interstellar medium.

Titanium dioxide, TiO$_2$, is a refractory species that could play a crucial role in the dust-condensation sequence around oxygen-rich evolved stars. To date, gas phase TiO$_2$ has been detected only in the complex environment of the red supergiant VY CMa. We aim to constrain the distribution and excitation of TiO$_2$ around VY CMa in order to clarify its role in dust formation. We analyse spectra and channel maps for TiO$_2$ extracted from ALMA science verification data. We detect 15 transitions of TiO$_2$, and spatially resolve the emission for the first time. The maps demonstrate a highly clumpy, anisotropic outflow in which the TiO$_2$ emission likely traces gas exposed to the stellar radiation field. A roughly east-west oriented, accelerating bipolar-like structure is found, of which the blue component runs into and breaks up around a solid continuum component. A distinct tail to the south-west is seen for some transitions, consistent with features seen in the optical and near-infrared. We find that a significant fraction of TiO$_2$ remains in the gas phase outside the dust-formation zone and suggest that this species might play only a minor role in the dust-condensation process around extreme oxygen-rich evolved stars like VY CMa.

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.

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.

Cool, evolved stars have copious, enriched winds. The structure of these winds and the way they are accelerated is not well known. We need to improve our understanding by studying the dynamics from the pulsating stellar surface to about 10 stellar radii, where radiation pressure on dust is fully effective. Some red supergiants have highly asymmetric nebulae, implicating additional forces. We retrieved ALMA Science Verification data providing images of sub-mm line and continuum emission from VY CMa. This enables us to locate water masers with milli-arcsec precision and resolve the dusty continuum. The 658-, 321- and 325-GHz masers lie in irregular, thick shells at increasing distances from the centre of expansion. For the first time this is confirmed as the stellar position, coinciding with a compact peak offset to the NW of the brightest continuum emission. The maser shells (and dust formation zone) overlap but avoid each other on tens-au scales. Their distribution is broadly consistent with excitation models but the conditions and kinematics appear to be complicated by wind collisions, clumping and asymmetries.

We report on the discovery of strong intensity variations in the high rotational lines of abundant molecular species towards the archetypical circumstellar envelope of IRC+10216. The observations have been carried out with the HIFI instrument on board \textit{Herschel}\thanks{\textit{Herschel} is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA} and with the IRAM\thanks{This work was based on observations carried out with the IRAM 30-meter telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain)} 30-m telescope. They cover several observing periods spreading over 3 years. The line intensity variations for molecules produced in the external layers of the envelope most probably result from time variations in the infrared pumping rates. We analyze the main implications this discovery has on the interpretation of molecular line emission in the envelopes of Mira-type stars. Radiative transfer calculations have to take into account both the time variability of infrared pumping and the possible variation of the dust and gas temperatures with stellar phase in order to reproduce the observation of molecular lines at different epochs. The effect of gas temperature variations with stellar phase could be particularly important for lines produced in the innermost regions of the envelope. Each layer of the circumstellar envelope sees the stellar light radiation with a different lag time (phase). Our results show that this effect must be included in the models. The sub-mm and FIR lines of AGB stars cannot anymore be considered as safe intensity calibrators.

We have mapped the CO(3-2) line emission around the Mira AB system at 0.5 resolution using the Atacama Large Millimeter/submillimeter Array (ALMA). The CO map shows amazing complexity. The circumstellar gas has been shaped by different dynamical actors during the evolution of the system and several morphological components can be identified. The companion is marginally resolved in continuum emission and is currently at 0.487$\pm$0.006 separation. In the main line component, centered on the stellar velocity, spiral arcs around Mira A are found. The spiral appears to be relatively flat and oriented in the orbital plane. An accretion wake behind the companion is clearly visible and the projected arc separation is of order 5''. In the blue wing of the line emission, offset from the main line, several large ($\sim$5-10''), opposing arcs are found. We tentatively suggest that this structure is created by the wind of Mira B blowing a bubble in the expanding envelope of Mira A.

To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. These observations are part of the Herschel Guaranteed Time Key Program HIFISTARS. We detected the J=6-5, J=10-9, and J=16-15 lines of 12CO and 13CO at about 100, 300, and 750K above the ground state (and the 13CO J=9-8 line). These lines are crucial for improving the modelling of the internal layers of the envelope around VY CMa. We also detected 27 lines of H2O and its isotopomers, and 96 lines of species such as NH3, SiO, SO, SO2 HCN, OH and others, some of them originating from vibrationally excited levels. Three lines were not unambiguously assigned. Our observations confirm that VY CMa's envelope must consist of two or more detached components. The molecular excitation in the outer layers is significantly lower than in the inner ones, resulting in strong self-absorbed profiles in molecular lines that are optically thick in this outer envelope, for instance, low-lying lines of H2O. Except for the most abundant species, CO and H2O, most of the molecular emission detected at these sub-mm/FIR wavelengths arise from the central parts of the envelope. The spectrum of VY CMa is very prominent in vibrationally excited lines, which are caused by the strong IR pumping present in the central regions. Compared with envelopes of other massive evolved stars, VY CMa's emission is particularly strong in these vibrationally excited lines, as well as in the emission from less abundant species such as H13CN, SO, and NH3.

Phosphorus-bearing compounds have only been studied in the circumstellar environments (CSEs) of the asymptotic giant branch (AGB) star IRC +10216 and the protoplanetary nebula CRL 2688, both C-rich objects, and the O-rich red supergiant VY CMa. The current chemical models cannot reproduce the high abundances of PO and PN derived from observations of VY CMa. No observations have been reported of phosphorus in the CSEs of O-rich AGB stars. We aim to set observational constraints on the phosphorous chemistry in the CSEs of O-rich AGB stars, by focussing on the Mira-type variable star IK Tau. Using the IRAM 30m telescope and the Submillimeter Array (SMA), we observed four rotational transitions of PN (J=2-1,3-2,6-5,7-6) and four of PO (J=5/2-3/2,7/2-5/2,13/2-11/2,15/2-13/2). The IRAM 30m observations were dedicated line observations, while the SMA data come from an unbiased spectral survey in the frequency range 279-355 GHz. We present the first detections of PN and PO in an O-rich AGB star and estimate abundances X(PN/H2) of about 3x10^-7 and X(PO/H2) in the range 0.5-6.0x10^-7. This is several orders of magnitude higher than what is found for the C-rich AGB star IRC +10216. The diameter (<=0.7") of the PN and PO emission distributions measured in the interferometric data corresponds to a maximum radial extent of about 40 stellar radii. The abundances and the spatial occurrence of the molecules are in very good agreement with the results reported for VY CMa. We did not detect PS or PH3 in the survey. We suggest that PN and PO are the main carriers of phosphorus in the gas phase, with abundances possibly up to several 10^-7. The current chemical models cannot account for this, underlining the strong need for updated chemical models that include phosphorous compounds.

We report the detection in emission of various highly excited rotational transitions of HNC (J = 6-5 through J =12-11) toward the carbon-star envelope IRC +10 216 using the HIFI instrument on-board the Herschel Space Observatory. Observations of the J = 1-0 and J = 3-2 lines of HNC with the IRAM 30-m telescope are also presented. The lines observed with HIFI have upper level energies corresponding to temperatures between 90 and 340 degrees Kelvin, and trace a warm and smaller circumstellar region than that seen in the interferometric maps of the J = 1-0 transition, whose emission extends up to a radius of 20". After a detailed chemical and radiative transfer modeling, we find that the presence of HNC in the circumstellar envelope of IRC +10 216 is consistent with formation from the precursor ion HCNH+, which in turn is produced through several proton transfer reactions which are triggered by the cosmic-ray ionization. We also find that the radiative pumping through 21 um photons to the first excited state of the bending mode v2 plays a crucial role to populate the high-J HNC levels involved in the transitions observed with HIFI. Emission in these high-J rotational transitions of HNC is expected to be strong in regions which are warm and dense and/or have an intense infrared flux at wavelengths around 21 um.

The C-rich AGB star IRC+10216 undergoes strong mass loss, and quasi-periodic density enhancements in the circumstellar matter have been reported. CO is ubiquitous in the CSE, while CCH emission comes from a spatially confined shell. With the IRAM 30m telescope and Herschel/HIFI, we recently detected unexpectedly strong emission from the CCH N=4-3, 6-5, 7-6, 8-7, and 9-8 transitions, challenging the available chemical and physical models. We aim to constrain the physical properties of IRC+10216's CSE, including the effect of episodic mass loss on the observed emission. In particular, we aim to determine the excitation region and conditions of CCH and to reconcile these with interferometric maps of the N=1-0 transition. Via radiative-transfer modelling, we provide a physical description of the CSE, constrained by the SED and a sample of 20 high-resolution and 29 low-resolution CO lines. We further present detailed radiative-transfer analysis of CCH. Assuming a distance of 150pc, the SED is modelled with a stellar luminosity of 11300Lsun and a dust-mass-loss rate of 4.0\times10^{-8}Msun/yr. Based on the analysis of 20 high resolution CO observations, an average gas-mass-loss rate for the last 1000yrs of 1.5\times10^{-5}Msun/yr is derived. This gives a gas-to-dust-mass ratio of 375, typical for an AGB star. The gas kinetic temperature throughout the CSE is described by 3 powerlaws: it goes as r^{-0.58} for r<9R*, as r^{-0.40} for 9<=r<=65R*, and as r^{-1.20} for r>65R*. This model successfully describes all 49 CO lines. We show the effect of wind-density enhancements on the CCH-abundance profile, and the good agreement of the model with the CCH N=1-0 transition and with the lines observed with the 30m telescope and HIFI. We report on the importance of radiative pumping to the vibrationally excited levels of CCH and the significant effect this has on the excitation of all levels of the CCH-molecule.

In 2001, the discovery of circumstellar water vapour around the ageing carbon star IRC+10216 was announced. This detection challenged the current understanding of chemistry in old stars, since water vapour was predicted to be absent in carbon-rich stars. Several explanations for the occurrence of water vapour were postulated, including the vaporization of icy bodies (comets or dwarf planets) in orbit around the star, grain surface reactions, and photochemistry in the outer circumstellar envelope. However, the only water line detected so far from one carbon-rich evolved star can not discriminate, by itself, between the different mechanisms proposed. Here we report on the detection by the Herschel satellite of dozens of water vapour lines in the far-infrared and sub-millimetre spectrum of IRC+10216, including some high-excitation lines with energies corresponding to ~1000 K. The emission of these high-excitation water lines can only be explained if water vapour is present in the warm inner region of the envelope. A plausible explanation for the formation of warm water vapour appears to be the penetration of ultraviolet (UV) photons deep into a clumpy circumstellar envelope. This mechanism triggers also the formation of other molecules such as ammonia, whose observed abundances are much higher than hitherto predicted.

In the framework of the HIFISTARS guaranteed time key programme, we measured more than 70 molecular emission lines with high signal-to-noise ratio towards VY CMa using the high-resolution HIFI spectrometer on board the Herschel satellite. The kinematic information obtained from the measured water lines supports the hypothesis of multiple outflow components. The observed high-intensity maser lines give no indication for strong polarisation.

Observations of high-excitation molecular emission lines can greatly increase our understanding of AGB winds, as they trace the innermost regions of the circumstellar envelope. The PACS spectrometer on-board the Herschel Space Telescope, provides for the first time the spectral resolution and sensitivity necessary to trace these lines. We report on the first modelling efforts of a PACS spectral scan for the OH/IR star V669 Cas. Central to our methodology is the consistent treatment of both dust and gas by using a line radiative transfer and a continuum radiative transfer code conjointly. Water emission lines are found to be extremely sensitive to the dust-to-gas ratio, emphasizing the need of consistent modelling for dust and gas.

CONTEXT: Water vapour maser emission from evolved oxygen-rich stars remains poorly understood. Additional observations, including polarisation studies and simultaneous observation of different maser transitions may ultimately lead to greater insight. AIMS: We have aimed to elucidate the nature and structure of the VY CMa water vapour masers in part by observationally testing a theoretical prediction of the relative strengths of the 620.701 GHz and the 22.235 GHz maser components of ortho water vapour. METHODS: In its high-resolution mode (HRS) the Herschel Heterodyne Instrument for the Infrared (HIFI) offers a frequency resolution of 0.125 MHz, corresponding to a line-of-sight velocity of 0.06 km/s, which we employed to obtain the strength and linear polarisation of maser spikes in the spectrum of VY CMa at 620.701 GHz. Simultaneous ground based observations of the 22.235 GHz maser with the Max-Planck-Institut f\"ur Radioastronomie 100-meter telescope at Effelsberg, provided a ratio of 620.701 GHz to 22.235 GHz emission. RESULTS:We report the first astronomical detection to date of water vapour maser emission at 620.701 GHz. In VY CMa both the 620.701 and the 22.235 GHz polarisation are weak. At 620.701 GHz the maser peaks are superposed on what appears to be a broad emission component, jointly ejected asymmetrically from the star. We observed the 620.701 GHz emission at two epochs 21 days apart, both to measure the potential direction of linearly polarised maser components and to obtain a measure of the longevity of these components. Although we do not detect significant polarisation levels in the core of the line, they rise up to approximately 6% in its wings.

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.

We aim to (1) set up simple and general analytical expressions to estimate mass-loss rates of evolved stars, and (2) from those calculate estimates for the mass-loss rates of asymptotic giant branch (AGB), red supergiant (RSG), and yellow hypergiant stars in our galactic sample. Rotationally excited lines of CO are a very robust diagnostic in the study of circumstellar envelopes (CSEs). When sampling different layers of the CSE, observations of these molecular lines lead to detailed profiles of kinetic temperature, expansion velocity, and density. A state-of-the-art, nonlocal thermal equilibrium, and co-moving frame radiative transfer code that predicts CO line intensities in the CSEs of late-type stars is used in deriving relations between stellar and molecular-line parameters, on the one hand, and mass-loss rate, on the other. We present analytical expressions for estimating the mass-loss rates of evolved stellar objects for 8 rotational transitions of the CO molecule, apply them to our extensive CO data set covering 47 stars, and compare our results to those of previous studies. Our expressions account for line saturation and resolving of the envelope, thereby allowing accurate determination of very high mass-loss rates. We argue that, for estimates based on a single rotational line, the CO(2-1) transition provides the most reliable mass-loss rate. The mass-loss rates calculated for the AGB stars range from 4x10^-8 Msun/yr up to 8x10^-5 Msun/yr. For RSGs they reach values between 2x10^-7 Msun/yr and 3x10^-4 Msun/yr. The estimates for the set of CO transitions allow time variability to be identified in the mass-loss rate. Possible mass-loss-rate variability is traced for 7 of the sample stars. We find a clear relation between the pulsation periods of the AGB stars and their derived mass-loss rates, with a levelling off at approx. 3x10^-5 Msun/yr for periods exceeding 850 days.

Circumstellar envelopes (CSEs) of a variety of evolved stars have been found to contain ammonia (NH3) in amounts that exceed predictions from conventional chemical models by many orders of magnitude. The observations reported here were performed in order to better constrain the NH3 abundance in the CSEs of four, quite diverse, oxygen-rich stars using the NH3 ortho J_K = 1_0 - 0_0 ground-state line. We used the Heterodyne Instrument for the Far Infrared aboard Herschel to observe the NH3 J_K = 1_0 - 0_0 transition near 572.5 GHz, simultaneously with the ortho-H2O J_Ka,Kc = 1_1,0 -1_0,1 transition, toward VY CMa, OH 26.5+0.6, IRC+10420, and IK Tau. We conducted non-LTE radiative transfer modeling with the goal to derive the NH3 abundance in these objects' CSEs. For the latter two stars, Very Large Array imaging of NH3 radio-wavelength inversion lines were used to provide further constraints, particularly on the spatial extent of the NH3-emitting regions. Results. We find remarkably strong NH3 emission in all of our objects with the NH3 line intensities rivaling those obtained for the ground state H2O line. The NH3 abundances relative to H2 are very high and range from 2 x 10-7 to 3 x 10-6 for the objects we have studied. Our observations confirm and even deepen the circumstellar NH3 enigma. While our radiative transfer modeling does not yield satisfactory fits to the observed line profiles, it leads to abundance estimates that confirm the very high values found in earlier studies. New ways to tackle this mystery will include further Herschel observations of more NH3 lines and imaging with the Expanded Very Large Array.

During their asymptotic giant branch, evolution low-mass stars lose a significant fraction of their mass through an intense wind, enriching the interstellar medium with products of nucleosynthesis. We observed the nearby oxygen-rich asymptotic giant branch star IK Tau using the high-resolution HIFI spectrometer onboard Herschel. We report on the first detection of H_2^{16}O and the rarer isotopologues H_2^{17}O and H_2^{18}O in both the ortho and para states. We deduce a total water content (relative to molecular hydrogen) of 6.6x10^{-5}, and an ortho-to-para ratio of 3:1. These results are consistent with the formation of H_2O in thermodynamical chemical equilibrium at photospheric temperatures, and does not require pulsationally induced non-equilibrium chemistry, vaporization of icy bodies or grain surface reactions. High-excitation lines of 12CO, 13CO, 28SiO, 29SiO, 30SiO, HCN, and SO have also been detected. From the observed line widths, the acceleration region in the inner wind zone can be characterized, and we show that the wind acceleration is slower than hitherto anticipated.

The interstellar medium is enriched primarily by matter ejected from evolved low and intermediate mass stars. The outflows from these stars create a circumstellar envelope in which a rich gas-phase and dust-nucleation chemistry takes place. We observed the nearest carbon-rich evolved star, IRC+10216, using the PACS (55-210 {\mu}m) and SPIRE (194-672 {\mu}m) spectrometers on board Herschel. We find several tens of lines from SiS and SiO, including lines from the v=1 vibrational level. For SiS these transitions range up to J=124-123, corresponding to energies around 6700K, while the highest detectable transition is J=90-89 for SiO, which corresponds to an energy around 8400K. Both species trace the dust formation zone of IRC+10216, and the broad energy ranges involved in their detected transitions permit us to derive the physical properties of the gas and the particular zone in which each species has been formed. This allows us to check the accuracy of chemical thermodynamical equilibrium models and the suggested depletion of SiS and SiO due to accretion onto dust grains.

We report on the detection of anhydrous hydrochloric acid (hydrogen chlorine, HCl) in the carbon-rich star IRC+10216 using the spectroscopic facilities onboard the Herschel satellite. Lines from J=1-0 up to J=7-6 have been detected. From the observed intensities, we conclude that HCl is produced in the innermost layers of the circumstellar envelope with an abundance relative to H2 of 5x10^-8 and extends until the molecules reach its photodissociation zone. Upper limits to the column densities of AlH, MgH, CaH, CuH, KH, NaH, FeH, and other diatomic hydrides have also been obtained.

Aims: Little information exists on the circumstellar molecular abundance stratifications of many molecules. The aim is to study the circumstellar chemical abundance pattern of 11 molecules and isotopologs ($^{12}$CO, $^{13}$CO, SiS, $^{28}$SiO, $^{29}$SiO, $^{30}$SiO, HCN, CN, CS, SO, SO$_2$) in the oxygen-rich evolved star IK~Tau. Methods: We have performed an in-depth analysis of a large number of molecular emission lines excited in the circumstellar envelope around IK~Tau. The analysis is done based on a non-local thermodynamic equilibrium (non-LTE) radiative transfer analysis, which calculates the temperature and velocity structure in a self-consistent way. The chemical abundance pattern is coupled to theoretical outer wind model predictions including photodestruction and cosmic ray ionization. Not only the integrated line intensities, but also the line shapes, are used as diagnostic tool to study the envelope structure. Results: The deduced wind acceleration is much slower than predicted from classical theories. SiO and SiS are depleted in the envelope, possibly due to the adsorption onto dust grains. For HCN and CS a clear difference with respect to inner wind non-equilibrium predictions is found, either indicating uncertainties in the inner wind theoretical modeling or the possibility that HCN and CS (or the radical CN) participate in the dust formation. The low signal-to-noise profiles of SO and CN prohibit an accurate abundance determination; the modeling of high-excitation SO$_2$ lines is cumbersome, possibly related to line misidentifications or problems with the collisional rates. The SiO isotopic ratios ($^{29}$SiO/$^{28}$SiO and $^{30}$SiO/$^{28}$SiO) point toward an enhancement in $^{28}$SiO compared to results of classical stellar evolution codes. Predictions for H$_2$O lines in the spectral range of the Herschel/HIFI mission are performed. [abbreviated]