We present WASP-43b climate simulations with deep wind jets (down to 700~bar) that are linked to retrograde (westward) flow at the equatorial day side for $p<0.1$~bar. Retrograde flow inhibits efficient eastward heat transport and naturally explains the small hotspot shift and large day-night-side gradient of WASP-43b ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days) observed with Spitzer. We find that deep wind jets are mainly associated with very fast rotations ($P_{\text{rot}}=P_{\text{orb}}\leq 1.5$~days) which correspond to the Rhines length smaller than $2$ planetary radii. We also diagnose wave activity that likely gives rise to deviations from superrotation. Further, we show that we can achieve full steady state in our climate simulations by imposing a deep forcing regime for $p>10$~bar: convergence time scale $\tau_{\text{conv}}=10^6-10^8$~s to a common adiabat, as well as linear drag at depth ($p\geq 200$~bar), which mimics to first order magnetic drag. Lower boundary stability and the deep forcing assumptions were also tested with climate simulations for HD~209458b ($P_{\text{orb}}=P_{\text{rot}}=3.5$~days). HD~209458b simulations always show shallow wind jets (never deeper than 100~bar) and unperturbed superrotation. If we impose a fast rotation ($P_{\text{orb}}=P_{\text{rot}}=0.8135$~days), also the HD~209458b-like simulation shows equatorial retrograde flow at the day side. We conclude that the placement of the lower boundary at $p=200$~bar is justified for slow rotators like HD~209458b, but we suggest that it has to be placed deeper for fast-rotating, dense hot Jupiters ($P_{\text{orb}}\leq 1.5$~days) like WASP-43b. Our study highlights that the deep atmosphere may have a strong influence on the observable atmospheric flow in some hot Jupiters.

{We investigate the atmospheric composition of the long period ($P_{\rm orb}=$~10 days), eccentric exo-Saturn WASP-117b. WASP-117b could be in atmospheric temperature and chemistry similar to WASP-107b. In mass and radius WASP-117b is similar to WASP-39b, which allows a comparative study of these planets.} % methods heading (mandatory) { We analyze a near-infrared transmission spectrum of WASP-117b taken with Hubble Space Telescope/WFC3 G141, which was reduced with two independent pipelines to ensure a robust detection of the $3 \sigma$-water spectrum. High resolution measurements were taken with VLT/ESPRESSO in the optical.} % results heading (mandatory) {Using a 1D atmosphere model with isothermal temperature, uniform cloud deck and equilibrium chemistry, the Bayesian evidence of a retrieval analysis of the transmission spectrum indicates a preference for a high atmospheric metallicity ${\rm [Fe/H]}=2.58^{+0.26}_{-0.37}$ and clear skies. The data are also consistent with a lower-metallicity composition ${\rm [Fe/H]}<1.75$ and a cloud deck between $10^{-2.2} - 10^{-5.1}$~bar, but with weaker Bayesian preference. We retrieve a low CH$_4$ abundance of $<10^{-4}$ volume fraction within $1 \sigma$. We cannot constrain the temperature between theoretically imposed limits of 700 and 1000~K. Further observations are needed to confirm quenching of CH$_4$ with $K_{zz}\geq 10^8$~cm$^2$/s. We report indications of Na and K in the VLT/ESPRESSO high resolution spectrum with substantial Bayesian evidence in combination with HST data.}

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

We present and analyse SO and SO$_2$, recently observed with high angular resolution and sensitivity in a spectral line survey with ALMA, for two oxygen-rich AGB stars: the low mass-loss rate R Dor and high mass-loss rate IK Tau. We analyse 8 lines of SO detected towards both stars, 78 lines of SO$_2$detected towards R Dor and 52 lines of SO$_2$ detected towards IK Tau. We detect several lines of $^{34}$SO, $^{33}$SO and $^{34}$SO$_2$ towards both stars, and tentatively S$^{18}$O towards R Dor, and hence derive isotopic ratios for these species. The spatially resolved observations show us that the two sulphur oxides are co-located towards R Dor and trace out the same wind structures in the circumstellar envelope (CSE). Much of the emission is well reproduced with a Gaussian abundance distribution spatially centred on the star. Emission from the higher energy levels of SO and SO$_2$ towards R Dor provide evidence in support of a rotating inner region of gas identified in earlier work. The new observations allow us to refine the abundance distribution of SO in IK Tau derived from prior observations with single antennas, and confirm the distribution is shell-like with the peak in the fractional abundance not centred on the star. The confirmation of different types of SO abundance distributions will help fine-tune chemical models and allows for an additional method to discriminate between low and high mass-loss rates for oxygen-rich AGB stars.

Asteroseismology offers the possibility of probing stellar interiors and testing evolutionary and seismic models. Precise photometry and spectroscopy obtained during multi-site campaigns on young open clusters allows discovering rich samples of pulsating stars and using them in a simultaneous seismic modelling called ensemble asteroseismology. The aim of this study is to obtain the age of the open cluster NGC 6910 by means of ensemble asteroseismology of the early-type pulsating members, to derive their stellar parameters, and to classify the excited modes. We used time-series analysis, performed photometric and spectroscopic mode identification, and calculated grids of evolutionary and seismic models to apply the procedure of ensemble asteroseismology for nine pulsating members of NGC 6910. With two iterations of the procedure of ensemble asteroseismology applied to nine pulsating stars we derived an age of 10.6$^{+0.9}_{-0.8}$ Myr for NGC 6910. Of the nine pulsating stars examined in the paper, eight are $\beta$ Cep stars, including three that are hybrid $\beta$ Cep and slowly pulsating B-type (SPB) pulsators, and one is an SPB star. Interestingly, the least massive $\beta$ Cep star, NGC 6910-38, has a mass of about 5.6 M$_\odot$. The present theory does not predict unstable $p$ modes in B-type stars with such a low mass. The $g$ modes with relatively high frequencies ($>3.5$ d$^{-1}$), observed in three members of the cluster, are also stable according to seismic modelling. Both findings pose a challenge for theoretical calculations and prompt a revision of the opacities. The procedure of ensemble asteroseismology was found to be successful for NGC 6910 and $\chi$ Per on the basis of pulsating B-type stars and can therefore be applied to other young open clusters that are rich in such stars.

The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run a grid of 70 3D hydrodynamics simulations of a progressively accelerating wind propagating in the Roche potential formed by a mass-loosing evolved star in orbit with a main sequence companion. We resolve the flow structure both in the immediate vicinity of the secondary, where bow shocks, outflows and wind-captured disks form, and up to 40 orbital separations, where spiral arms, arcs and equatorial density enhancements develop. When the companion is deeply engulfed in the wind, the lower terminal wind speeds and more progressive wind acceleration around oxygen-rich AGB stars make them more prone than carbon-rich AGB stars to display more disturbed outflows, a disk-like structure around the companion and a wind concentrated in the orbital plane. In these configurations, a large fraction of the wind is captured by the companion which leads to a significant shrinking of the orbit over the mass-loss timescale, if the donor star is at least a few times more massive than its companion. Provided the companion has a mass of at least a tenth of the mass of the donor star, it can compress the wind in the orbital plane up to large distances. Our grid of models covers a wide scope of configurations function of the dust chemical content, the terminal wind speed relative to the orbital speed, the extension of the dust condensation region around the cool evolved star and the mass ratio. It provides a frame of reference to interpret high-resolution maps of the outflows surrounding cool evolved stars.

Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster.

Rate coefficients have been measured for the reaction of CH radicals with formaldehyde, CH$_{2}$O, over the temperature range 31 - 133 K using a pulsed Laval nozzle apparatus combined with pulsed laser photolysis and laser induced fluorescence spectroscopy. The rate coefficients are very large and display a distinct decrease with decreasing temperature below 70 K, although classical collision rate theory fails to reproduce this temperature dependence. The measured rate coefficients have been parameterized and used as input for astrochemical models for both dark cloud and AGB stellar outflow scenarios. The models predict a distinct change (up to a factor of two) in the abundance of ketene, H$_{2}$CCO, which is the major expected molecular product of the CH + CH$_{2}$O reaction.

Chemical modelling of AGB outflows is typically focused on either non-thermodynamic equilibrium chemistry in the inner region or photon-driven chemistry in the outer region. We include, for the first time, a comprehensive dust-gas chemistry in our AGB outflow chemical kinetics model, including both dust-gas interactions and grain-surface chemistry. The dust is assumed to have formed in the inner region, and follows an interstellar-like dust-size distribution. Using radiative transfer modelling, we obtain dust temperature profiles for different dust types in an O-rich and a C-rich outflow. We calculate a grid of models, sampling different outflow densities, drift velocities between the dust and gas, and dust types. Dust-gas chemistry can significantly affect the gas-phase composition, depleting parent and daughter species and increasing the abundance of certain daughter species via grain-surface formation followed by desorption/sputtering. Its influence depends on four factors: outflow density, dust temperature, initial composition, and drift velocity. The largest effects are for higher density outflows with cold dust and O-rich parent species, as these species generally have a larger binding energy. At drift velocities larger than $\sim 10$ km s$^{-1}$, ice mantles undergo sputtering; however, they are not fully destroyed. Models with dust-gas chemistry can better reproduce the observed depletion of species in O-rich outflows. When including colder dust in the C-rich outflows and adjusting the binding energy of CS, the depletion in C-rich outflows is also better reproduced. To best interpret high-resolution molecular line observations from AGB outflows, dust-gas interactions are needed in chemical kinetics models.

Unravelling the composition and characteristics of gas and dust lost by asymptotic giant branch (AGB) stars is important as these stars play a vital role in the chemical life cycle of galaxies. The general hypothesis of their mass loss mechanism is a combination of stellar pulsations and radiative pressure on dust grains. However, current models simplify dust formation, which starts as a microscopic phase transition called nucleation. Various nucleation theories exist, yet all assume chemical equilibrium, growth restricted by monomers, and commonly use macroscopic properties for a microscopic process. Such simplifications for initial dust formation can have large repercussions on the type, amount, and formation time of dust. By abandoning equilibrium assumptions, discarding growth restrictions, and using quantum mechanical properties, we have constructed and investigated an improved nucleation theory in AGB wind conditions for four dust candidates, TiO$_2$, MgO, SiO and Al$_2$O$_3$. This paper reports the viability of these candidates as first dust precursors and reveals implications of simplified nucleation theories. Monomer restricted growth underpredicts large clusters at low temperatures and overpredicts formation times. Assuming the candidates are present, Al$_2$O$_3$ is the favoured precursor due to its rapid growth at the highest considered temperatures. However, when considering an initially atomic chemical mixture, only TiO$_2$-clusters form. Still, we believe Al$_2$O$_3$ to be the prime candidate due to substantial physical evidence in presolar grains, observations of dust around AGB stars at high temperatures, and its ability to form at high temperatures and expect the missing link to be insufficient quantitative data of Al-reactions.

Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss which spreads heavy elements into the interstellar medium. The structure and the dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO $J = 2-1$ line from the red supergiant star $\mu$~Cep with the NOEMA interferometer. In the line the synthesized beam was $0.92 \times 0.72$~arcsec ($590 \times 462$~au at 641~pc). The continuum map shows only the unresolved contribution of the free-free emission of the star chromosphere. The continuum-subtracted channel maps reveal a very inhomogeneous and clumpy circumstellar environment. In particular, we detected a bright CO clump, as bright as the central source in the line, at 1.80~arcsec south-west from the star, in the blue channel maps. After a deprojection of the radial velocity assuming two different constant wind velocities, the observations were modelled using the 3D radiative transfer code \textsc{lime} to derive the characteristics of the different structures. We determine that the gaseous clumps observed around $\mu$~Cep are responsible for a mass loss rate of $(4.9 \pm 1.0) \times 10^{-7}~{\rm M}_\odot\,{\rm yr}^{-1}$, in addition to a spatially unresolved wind component with an estimated mass-loss rate of $2.0 \times 10^{-6}~{\rm M}_\odot\,{\rm yr}^{-1}$. Therefore, the clumps have a significant role in $\mu$~Cep's mass loss ($\ge 25 \%$). We cannot exclude that the unresolved central outflow may be made of smaller unresolved clumps.

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.

A molecular line survey has been carried out toward the carbon-rich asymptotic giant branch star CW Leo employing the HIFI instrument on board of the Herschel satellite. Numerous features from 480 GHz to beyond 1100 GHz could be assigned unambiguously to the fairly floppy SiC$_2$ molecule. However, predictions from laboratory data exhibited large deviations from the observed frequencies even after some lower frequency data from this survey were incorporated into a fit. Therefore, we present a combined fit of all available laboratory data together with data from radio-astronomical observations.

In 1981, the idea of a superwind that ends the life of cool giant stars was proposed. Extreme OH/IR-stars develop superwinds with the highest mass-loss rates known so far, up to a few 10^(-4) Msun/yr, informing our understanding of the maximum mass-loss rate achieved during the Asymptotic Giant Branch (AGB) phase. A condundrum arises whereby the observationally determined duration of the superwind phase is too short for these stars to become white dwarfs. Here, we report on the detection of spiral structures around two cornerstone extreme OH/IR-stars, OH26.5+0.6 and OH30.1-0.7, identifying them as wide binary systems. Hydrodynamical simulations show that the companion's gravitational attraction creates an equatorial density enhancement mimicking a short extreme superwind phase, thereby solving the decades-old conundrum. This discovery restricts the maximum mass-loss rate of AGB stars around the single-scattering radiation-pressure limit of a few 10^(-5) Msun/yr. This brings about crucial implications for nucleosynthetic yields, planet survival, and the wind-driving mechanism.

We aim to determine the distributions of molecular SiS and CS in the circumstellar envelopes of oxygen-rich asymptotic giant branch stars and how these distributions differ between stars that lose mass at different rates. In this study we analyse ALMA observations of SiS and CS emission lines for three oxygen-rich galactic AGB stars: IK Tau, with a moderately high mass-loss rate of $5\times10^{-6}$M$_\odot$ yr$^{-1}$, and W Hya and R Dor with low mass loss rates of $\sim1\times10^{-7}$M$_\odot$ yr$^{-1}$. These molecules are usually more abundant in carbon stars but the high sensitivity of ALMA allows us to detect their faint emission in the low mass-loss rate AGB stars. The high spatial resolution of ALMA also allows us to precisely determine the spatial distribution of these molecules in the circumstellar envelopes. We run radiative transfer models to calculate the molecular abundances and abundance distributions for each star. We find a spread of peak SiS abundances with $\sim10^{-8}$ for R Dor, $\sim10^{-7}$ for W Hya, and $\sim3\times10^{-6}$ for IK Tau relative to H$_2$. We find lower peak CS abundances of $\sim7\times10^{-9}$ for R Dor, $\sim7\times10^{-8}$ for W Hya and $\sim4\times10^{-7}$ for IK Tau, with some stratifications in the abundance distributions. For IK Tau we also calculate abundances for the detected isotopologues: C$^{34}$S, $^{29}$SiS, $^{30}$SiS, Si$^{33}$S, Si$^{34}$S, $^{29}$Si$^{34}$S, and $^{30}$Si$^{34}$S. Overall the isotopic ratios we derive for IK Tau suggest a lower metallicity than solar.

The material lost through stellar winds of Asymptotic Giant Branch (AGB) stars is one of the main contributors to the chemical enrichment of galaxies. The general hypothesis of the mass loss mechanism of AGB winds is a combination of stellar pulsations and radiative pressure on dust grains, yet current models still suffer from limitations. Among others, they assume chemical equilibrium of the gas, which may not be justified due to rapid local dynamical changes in the wind. This is important as it is the chemical composition that regulates the thermal structure of the wind, the creation of dust grains in the wind, and ultimately the mass loss by the wind. Using a self-consistent hydrochemical model, we investigated how non-equilibrium chemistry affects the dynamics of the wind. This paper compares a hydrodynamical and a hydrochemical dust-free wind, with focus on the chemical heating and cooling processes. No sustainable wind arises in a purely hydrodynamical model with physically reasonable pulsations. Moreover, temperatures are too high for dust formation to happen, rendering radiative pressure on grains impossible. A hydrochemical wind is even harder to initiate due to efficient chemical cooling. However, temperatures are sufficiently low in dense regions for dust formation to take place. These regions occur close to the star, which is needed for radiation pressure on dust to sufficiently aid in creating a wind. Extending this model self-consistently with dust formation and evolution, and including radiation pressure, will help to understand the mass loss by AGB winds.

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

Context: At the end of their lives AGB stars are prolific producers of dust and gas. The details of this mass-loss process are still not understood very well. Herschel PACS and SPIRE spectra offer a unique way of investigating properties of AGB stars in general and the mass-loss process in particular. Methods: The HIPE software with the latest calibration is used to process the available PACS and SPIRE spectra of 40 evolved stars. The spectra are convolved with the response curves of the PACS and SPIRE bolometers and compared to the fluxes measured in imaging data of these sources. Custom software is used to identify lines in the spectra, and to determine the central wavelengths and line intensities. Standard molecular line databases are used to associate the observed lines. Because of the limited spectral resolution of the spectrometers several known lines are typically potential counterparts to any observed line. To help identifications the relative contributions in line intensity of the potential counterpart lines are listed for three characteristic temperatures based on LTE calculations and assuming optically thin emission. Result: The following data products are released: the reduced spectra, the lines that are measured in the spectra with wavelength, intensity, potential identifications, and the continuum spectra, i.e. the full spectra with all identified lines removed. As simple examples of how this data can be used in future studies we have fitted the continuum spectra with three power laws and find that the few OH/IR stars seem to have significantly steeper slopes than the other oxygen- and carbon-rich objects in the sample. As another example we constructed rotational diagrams for CO and fitted a two-component model to derive rotational temperatures.

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.

We report on the succesful search for CO (2-1) and (3-2) emission associated with OH/IR stars in the Galactic Bulge. We observed a sample of eight extremely red AGB stars with the APEX telescope and detected seven. The sources were selected at sufficient high galactic latitude to avoid interference by interstellar CO, which hampered previous studies of inner galaxy stars. To study the nature of our sample and the mass loss we constructed the SEDs from photometric data and Spitzer IRS spectroscopy. In a first step we apply radiative transfer modelling to fit the SEDs and obtain luminosities and dust mass loss rates (MLR). Through dynamical modelling we then retrieve the total MLR and the gas-to-dust ratios. We derived variability periods of our stars. The luminosities range between approximately 4000 and 5500 Lsun and periods are below 700 days. The total MLR ranges between 1E-5 and 1E-4 Msun/yr. Comparison with evolutionary models shows that the progenitor mass is approximately 1.5 Msun, similar to the Bulge Miras and are of intermediate age (3 Gyr). The gas-to-dust ratios are between 100 and 400 and are similar to what is found for OH/IR stars in the galactic Disk. One star, IRAS 17347-2319, has a very short period of approximately 300 days which may be decreasing further. It may belong to a class of Mira variables with a sudden change in period as observed in some Galactic objects. It would be the first example of an OH/IR star in this class and deserves further follow-up observations.

High-resolution interferometric observations of the circumstellar environments of AGB stars show a variety of morphologies. Guided by the unusual carbon monoxide line profile of the AGB star EP Aquarii, we have observed its circumstellar environment with ALMA band 6 in cycle 4. We describe the morphological complexity of the CO, SiO and SO2 molecular emission. The CO emission exhibits the characteristics of a bi-conical wind with a bright nearly face-on spiral feature around the systemic velocity. This is the first convincing detection of a spiral morphology in an O-rich wind. Based on the offsets of the centers of the two bi-conical wind hemispheres, we deduce the position angle of the inclination axis to be ~150 deg measured counterclockwise from north. Based on the velocity width of the spiral signature we estimate the inclination angle of the system to be between 4 deg and 18 deg. The central emission zone exhibits a morphology that resembles simulations modelling the spiral-inducing wind-Roche-lobe-overflow mechanism. Though the spiral may be a (companion-induced) density enhancement in the stellar outflow, the extremely narrow width of the spiral signature in velocity space suggests that it may be a hydrodynamical perturbation in a face-on differentially rotating disk. The SiO emission does not show the spiral, but exhibits a local emission void approximately 0.5'' west of the continuum brightness peak. We hypothesise that this may be a local environment caused by the presence of a (stellar) companion with a mass of at most 0.1 Msol, based on its non-detection in the continuum...

A spectral scan of the circumstellar environment of the asymptotic giant branch (AGB) star R~Doradus was taken with ALMA in cycle 2 at frequencies between 335 and 362 GHz and with a spatial resolution of $\sim$150 milliarcseconds. Many molecular lines show a spatial offset between the blue and red shifted emission in the innermost regions of the wind. The position-velocity diagrams of this feature, in combination with previous SPHERE data and theoretical work point towards the presence of a compact differentially rotating disk, orientated nearly edge-on. We model the $^{\rm 28}$SiO ($v=1,~J=8\to7$) emission with a disk model. We estimate the disk mass and angular momentum to be $3 \times 10^{-6}$ Solar masses and $5 \times 10^{40}\ {\rm m^2 kg/s}$. The latter presents an `angular momentum problem' that may be solved by assuming that the disk is the result of wind-companion interactions with a companion of at least 2.5 earth masses, located at 6 AU, the tentatively determined location of the disk's inner rim. An isolated clump of emission is detected to the south-east with a velocity that is high compared to the previously determined terminal velocity of the wind. Its position and mean velocity suggest that it may be associated with a companion planet, located at the disk's inner rim.

The circumstellar environment of the AGB star L$_{\rm 2}$ Puppis was observed with ALMA in cycle 3, with a resolution of $15 \times 18 \rm\ mas$. The molecular emission shows a differentially rotating disk, inclined to a nearly edge-on position. In the first paper in this series (paper I) the molecular emission was analysed to accurately deduce the motion of the gas in the equatorial regions of the disk. In this work we model the optically thick $^{\rm 12}$CO $J=$3$-$2 and the optically thin $^{\rm 13}$CO $J=$3$-$2 rotational transition to constrain the physical conditions in the disk. To realise this effort we make use of the 3D NLTE radiative transfer code {\tt LIME}. The temperature structure and velocity structure show a high degree of complexity, both radially and vertically. The radial H$_{\rm 2}$ density profile in the disk plane is characterised by a power law with a slope of $-3.1$. We find a $^{\rm 12}$CO over $^{\rm 13}$CO abundance ratio of 10 inside the disk. Finally, estimations of the angular momentum in the disk surpass the expected available angular momentum of the star, strongly supporting the indirect detection of a compact binary companion reported in paper I. We estimate the mass of the companion to be around 1 Jupiter mass.

Equatorial density enhancements (EDEs) are a very common astronomical phenomenon. Studies of the circumstellar environments (CSE) of young stellar objects and of evolved stars have shown that these objects often possess these features. These are believed to originate from different mechanisms, ranging from binary interactions to the gravitational collapse of interstellar material. Quantifying the effect of the presence of this type of EDE on the observables is essential for a correct interpretation of high-resolution data. We seek to investigate the manifestation in the observables of a circumstellar EDE, to assess which properties can be constrained, and to provide an intuitive bedrock on which to compare and interpret upcoming high-resolution data (e.g. \emph{ALMA} data) using 3D models. We develop a simplified analytical parametrised description of a 3D EDE, with a possible substructure such as warps, gaps, and spiral instabilities. In addition, different velocity fields (Keplerian, radial, super-Keplerian, sub-Keplerian and rigid rotation) are considered. The effect of a bipolar outflow is also investigated. The geometrical models are fed into the 3D radiative transfer code {\tt LIME}, that produces 3D intensity maps throughout velocity space. We investigate the spectral signature of the $J$=3$-$2 up to $J$=7$-$6 rotational transitions of CO in the models, as well as the spatial aspect of this emission by means of channel maps, wide-slit position-velocity (PV) diagrams, stereograms, and spectral lines. Additionally, we discuss methods of constraining the geometry of the EDE, the inclination, the mass-contrast between the EDE and the bipolar outflow, and the global velocity field. Finally, we simulated \emph{ALMA} observations...

(abridged) In the inner regions of AGB outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium (TE) chemical models. The presence of the majority of these species can be explained by shock-induced non-TE chemical models, where shocks caused by the pulsating star take the chemistry out of TE in the inner region. Moreover, a non-uniform density structure has been detected in several AGB outflows. A detailed parameter study on the quantitative effects of a non-homogeneous outflow has so far not been performed. We implement a porosity formalism for treating the increased leakage of light associated with radiation transport through a clumpy, porous medium. The effects from the altered UV radiation field penetration on the chemistry, accounting also for the increased reaction rates of two-body processes in the overdense clumps, are examined. We present a parameter study of the effect of clumping and porosity on the chemistry throughout the outflow. Both the higher density within the clumps and the increased UV radiation field penetration have an important impact on the chemistry, as they both alter the chemical pathways. The increased amount of UV radiation in the inner region leads to photodissociation of parent species, releasing the otherwise deficient elements. We find an increased abundance in the inner region of all species not expected to be present assuming TE chemistry, such as HCN in O-rich outflows, H$_2$O in C-rich outflows, and NH$_3$ in both. Outflows whose clumps have a large overdensity and that are very porous to the interstellar UV radiation field yield abundances comparable to those observed in O- and C-rich outflows for most of the unexpected species investigated. The inner wind abundances of H$_2$O in C-rich outflows and of NH$_3$ in O- and C-rich outflows are however underpredicted.

A spectral line and imaging survey of the low mass-loss rate AGB star R Dor (Mdot ~ 1e-7 Msun/yr) and the high mass-loss rate AGB star IK Tau (Mdot ~5e-6 Msun/yr) was made with ALMA between 335 and 362 GHz at a spatial resolution of ~150 mas, corresponding to the locus of the main dust formation region of both targets. Some 200 spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states. Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2. We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v_inf, previously derived for both objects (36 km/s versus v_inf ~17.7 km/s for IK Tau and 23 km/s versus v_inf ~5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond ~400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ~3.7e-7 Msun for IK Tau and ~2e-8 Msun for R Dor.

Using ALMA, we observed the stellar wind of two oxygen-rich Asymptotic Giant Branch (AGB) stars, IK Tau and R Dor, between 335 and 362 GHz. One aim was to detect metal oxides and metal hydroxides (AlO, AlOH, FeO, MgO, MgOH), some of which are thought to be direct precursors of dust nucleation and growth. We report on the potential first detection of FeO (v=0, Omega=4, J=11-10) in RDor (mass-loss rate, Mdot, ~1e-7 Msun/yr). The presence of FeO in IK Tau (Mdot~5e-6 Msun/yr) cannot be confirmed due to a blend with 29SiS, a molecule that is absent in R Dor. The detection of AlO in R Dor and of AlOH in IK Tau was reported earlier by Decin et al. (2017). All other metal oxides and hydroxides, as well as MgS, remain undetected. We derive a column density N(FeO) of 1.1+/-0.9e15 cm^{-2} in R Dor, or a fractional abundance [FeO/H]~1.5e-8 accounting for non-LTE effects. The derived fractional abundance [FeO/H] is a factor ~20 larger than conventional gas-phase chemical kinetic predictions. This discrepancy may be partly accounted for by the role of vibrationally excited OH in oxidizing Fe, or may be evidence for other currently unrecognised chemical pathways producing FeO. Assuming a constant fractional abundance w.r.t. H_2, the upper limits for the other metals are [MgO/H_2] <5.5e-10 (R Dor) and <7e-11 (IK Tau), [MgOH/H_2] <9e-9 (R Dor) and <1e-9 (IK Tau), [CaO/H_2] <2.5e-9 (R Dor) and <1e-10 (IK Tau), [CaOH/H_2] <6.5e-9 (R Dor) and <9e-10 (IK Tau), and [MgS/H_$] <4.5e-10 (R Dor) and <6e-11 (IK Tau). The retrieved upper limit abundances for these latter molecules are in accord with the chemical model predictions.

Stratosphere circulation is important to interpret abundances of photo-chemically produced compounds like ozone that we aim to observe to assess habitability of exoplanets. We thus investigate a tidally locked ExoEarth scenario for TRAPPIST-1b, TRAPPIST-1d, Proxima Centauri~b and GJ 667 C~f with a simplified 3D atmosphere model and for different stratospheric wind breaking assumptions. These planets are representatives for different circulation regimes for orbital periods: $P_{orb}=1-100$~days. The circulation of exoplanets with $P_{orb} \leq $ 25~days can be dominated by the standing tropical Rossby wave in the troposphere and also in the stratosphere: It leads to a strong equatorial eastward wind jet and to 'Anti-Brewer-Dobson'-circulation that confines air masses to the stratospheric equatorial region. Thus, the distribution of photo-chemically produced species and aerosols may be limited to an 'equatorial transport belt'. In contrast, planets with $P_{orb}>25$~days, like GJ~667~C~f, exhibit efficient thermally driven circulation in the stratosphere that allows for a day side-wide distribution of air masses. The influence of the standing tropical Rossby waves on tidally locked ExoEarths with $P_{orb} \leq 25$~days can, however, be circumvented with deep stratospheric wind breaking alone - allowing for equator-to-pole transport like on Earth. For planets with $3 \leq P_{orb} \leq 6$~days, the extratropical Rossby wave acts as an additional safe-guard against the tropical Rossby wave in case of shallow wind breaking. Therefore, TRAPPIST-1d is less prone to have an equatorial transport belt in the stratosphere than Proxima~Centauri~b. Even our Earth model shows an equatorial wind jet, if stratosphere wind breaking is inefficient.

We observed Betelgeuse using ALMA's extended configuration in band 7 (f~340 GHz, {\lambda}~0.88 mm), resulting in a very high angular resolution of 18 mas. Using a solid body rotation model of the 28SiO(v=2,J=8-7) line emission, we show that the supergiant is rotating with a projected equatorial velocity of v_eq sin i = 5.47 +/- 0.25 km/s at the equivalent continuum angular radius R_star = 29.50 +/- 0.14 mas. This corresponds to an angular rotation velocity of {\omega} sin i = (5.6 +/- 1.3) x 10^(-9) rad/s. The position angle of its north pole is PA = 48.0 +/- 3.5{\deg}. The rotation period of Betelgeuse is estimated to P/sin i = 36 +/- 8 years. The combination of our velocity measurement with previous observations in the ultraviolet shows that the chromosphere is co-rotating with the star up to a radius of ~10 au (45 mas or 1.5x the ALMA continuum radius). The coincidence of the position angle of the polar axis of Betelgeuse with that of the major ALMA continuum hot spot, a molecular plume, and a partial dust shell (from previous observations) suggests that focused mass loss is currently taking place in the polar region of the star. We propose that this hot spot corresponds to the location of a particularly strong "rogue" convection cell, which emits a focused molecular plume that subsequently condenses into dust at a few stellar radii. Rogue convection cells therefore appear to be an important factor shaping the anisotropic mass loss of red supergiants.

We study the sub-Keplerian rotation and dust content of the circumstellar material around the asymptotic giant branch (AGB) star L$_2$ Puppis. We find that the thermal pressure gradient alone cannot explain the observed rotation profile. We find that there is a family of possible dust populations for which radiation pressure can drive the observed sub-Keplerian rotation. This set of solutions is further constrained by the spectral energy distribution (SED) of the system, and we find that a dust-to-gas mass ratio of $\sim10^{-3}$ and a maximum grain size that decreases radially outwards can satisfy both the rotation curve and SED. These dust populations are dynamically tightly coupled to the gas azimuthally. However grains larger than $\sim0.5\mu$m are driven outward radially by radiation pressure at velocities $\sim5$kms$^{-1}$, which implies a dust replenishment rate of $\sim3\times10^{-9}$M$_\odot$yr$^{-1}$. This replenishment rate is consistent with observational estimates to within uncertainties. Coupling between the radial motion of the dust and gas is weak and hence the gas does not share in this rapid outward motion. Overall we conclude that radiation pressure is a capable and necessary mechanism to explain the observed rotation profile of L$_2$ Pup, and offers other additional constraints on the dust properties.

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

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.

The mechanisms responsible for heating the extended atmospheres of early-M spectral-type supergiants are poorly understood. So too is the subsequent role these mechanisms play in driving the large mass-loss rates of these stars. Here we present ALMA long (i.e., $\sim$16 km) baseline 338 GHz (0.89 mm) continuum observations of the free-free emission in the extended atmosphere of the M2 spectral-type supergiant Betelgeuse. The spatial resolution of 14 mas exquisitely resolves the atmosphere, revealing it to have a mean temperature of 2760 K at $\sim$1.3 R$_{\star}$, which is below both the photospheric effective temperature ($T_{\textrm{eff}} = 3690$ K) and the temperatures at $\sim$$2 $R$_{\star}$. This is unambiguous proof for the existence of an inversion of the mean temperature in the atmosphere of a red supergiant. The emission is clearly not spherically symmetric with two notable deviations from a uniform disk detected in both the images and visibilities. The most prominent asymmetry is located in the north-east quadrant of the disk and is spatially resolved showing it to be highly elongated with an axis-ratio of 2.4 and occupying $\sim$$5\%$ of the disk projected area. Its temperature is approximately 1000 K above the measured mean temperature at 1.3 R$_{\star}$. The other main asymmetry is located on the disk limb almost due east of the disk center and occupies $\sim$$3\%$ of the disk projected area. Both emission asymmetries are clear evidence for localized heating taking place in the atmosphere of Betelgeuse. We suggest that the detected localized heating is related to magnetic activity generated by large-scale photospheric convection.

We provide an update for the astrometric solution for the Type II supernova progenitor Betelgeuse using the revised Hipparcos Intermediate Astrometric Data (HIAD) of van Leeuwen, combined with existing VLA and new e-MERLIN and ALMA positions. The 2007 Hipparcos refined abscissa measurements required the addition of so-called Cosmic Noise of 2.4 mas to find an acceptable 5-parameter stochastic solution. We find that a measure of radio Cosmic Noise should also be included for the radio positions because surface inhomogeneities exist at a level significant enough to introduce additional intensity centroid uncertainty. Combining the 2007 HIAD with the proper motions based solely on the radio positions leads to a parallax of $\pi = 5.27\pm 0.78$ mas ($190^{+33}_{-25}$ pc), smaller than the Hipparcos 2007 value of $6.56\pm 0.83$ mas ($152^{+22}_{-17}$ pc; van Leeuwen 2007). Furthermore, combining the VLA and new e-MERLIN and ALMA radio positions with the 2007 HIAD, and including radio Cosmic Noise of 2.4 mas, leads to a nominal parallax solution of $4.51 \pm 0.80$ mas ($222^{+48}_{-34}$ pc), which while only $0.7\sigma$ different from the 2008 solution of Harper et al. it is $2.6\sigma$ different from the solution of van Leeuwen. An accurate and precise parallax for Betelgeuse is always going to be difficult to obtain because it is small compared to the stellar angular diameter ($\theta=44$ mas). We outline an observing strategy, utilizing future mm and sub-mm high-spatial resolution interferometry that must be used if substantial improvements in precision and accuracy of the parallax and distance are to be achieved.

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

AIM: In this study we examine rotational emission lines of two isotopologues of water: H$_2$$^{17}$O and H$_2$$^{18}$O. By determining the abundances of these molecules, we aim to use the derived isotopologue --- and hence oxygen isotope --- ratios to put constraints on the masses of a sample of M-type AGB stars that have not been classified as OH/IR stars. METHODS: We use detailed radiative transfer analysis based on the accelerated lambda iteration method to model the circumstellar molecular line emission of H$_2$$^{17}$O and H$_2$$^{18}$O for IK Tau, R Dor, W Hya, and R Cas. The emission lines used to constrain our models come from Herschel/HIFI and Herschel/PACS observations and are all optically thick, meaning that full radiative transfer analysis is the only viable method of estimating molecular abundance ratios. RESULTS: We find generally low values of the $^{17}$O/$^{18}$O ratio for our sample, ranging from 0.15 to 0.69. This correlates with relatively low initial masses, in the range $\sim1.0$ to 1.5 M$_\odot$ for each source, based on stellar evolutionary models. We also find ortho-to-para ratios close to 3, which are expected from warm formation predictions. CONCLUSIONS: The $^{17}$O/$^{18}$O ratios found for this sample are at the lower end of the range predicted by stellar evolutionary models, indicating that the sample chosen had relatively low initial masses.

The aim of this paper is to investigate the $^{17}$O/$^{18}$O ratio for a sample of AGB stars, containing M-, S- and C-type stars. These ratios are evaluated in relation to fundamental stellar evolution parameters: the stellar initial mass and pulsation period. Circumstellar $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O line observations were obtained for a sample of nine stars with various single-dish long-wavelength facilities. Line intensity ratios are shown to relate directly to the surface $^{17}$O/$^{18}$O abundance ratio. Stellar evolution models predict the $^{17}$O/$^{18}$O ratio to be a sensitive function of initial mass and to remain constant throughout the entire TP-AGB phase for stars initially less massive than 5\,$M_{\odot}$. This makes the measured ratio a probe of the initial stellar mass. Observed $^{17}$O/$^{18}$O ratios are found to be well in the range predicted by stellar evolution models that do not consider convective overshooting. From this, accurate initial mass estimates are calculated for seven sources. For the remaining two sources two mass solutions result, though with a larger probability that the low-mass solution is the correct one. Finally, hints at a possible separation between M/S- and C-type stars when comparing the $^{17}$O/$^{18}$O ratio to the stellar pulsation period are presented.

Six billion years from now, while evolving on the asymptotic giant branch (AGB), the Sun will metamorphose from a red giant into a beautiful planetary nebula. This spectacular evolution will impact the Solar System planets, but observational confirmations of the predictions of evolution models are still elusive as no planet orbiting an AGB star has yet been discovered. The nearby AGB red giant L2 Puppis (d = 64 pc) is surrounded by an almost edge-on circumstellar dust disk. We report new observations with ALMA at very high angular resolution (18 x 15 mas) in band 7 (f ~ 350 GHz) that allow us to resolve the velocity profile of the molecular disk. We establish that the gas velocity profile is Keplerian within the central cavity of the dust disk, allowing us to derive the mass of the central star L2 Pup A, mA = 0.659 +/- 0.011 +/- 0.041 Msun (+/- 6.6%). From evolutionary models, we determine that L2 Pup A had a near-solar main sequence mass, and is therefore a close analog of the future Sun in 5 to 6 Gyr. The continuum map reveals the presence of a secondary source (B) at a radius of 2 AU contributing fB/ fA = 1.3 +/- 0.1% of the flux of the AGB star. L2 Pup B is also detected in CO emission lines at a radial velocity of vB = 12.2 +/- 1.0 km/s. The close coincidence of the center of rotation of the gaseous disk with the position of the continuum emission from the AGB star allows us to constrain the mass of the companion to mB = 12 +/- 16 MJup. L2 Pup B is most likely a planet or low mass brown dwarf with an orbital period around 5 years. Its continuum brightness and molecular emission suggest that it may be surrounded by an extended molecular atmosphere or an accretion disk. L2 Pup therefore emerges as a promising vantage point on the distant future of our Solar System.

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.

Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models.

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.

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.

More than 3000 exoplanets have been detected so far, and more and more spectroscopic observations of exoplanets are performed. Future instruments are eagerly awaited as they will be able to provide spectroscopic data with a greater accuracy and sensitivity than what is currently available. An important aspect to consider is temporal stellar atmospheric disturbances that can influence the planetary composition, and hence spectra, and potentially can lead to incorrect assumptions about the steady-state atmospheric composition of the planet. We focus on perturbations that come from the host star in the form of flare events that significantly increase the photon flux impingement on the exoplanet atmosphere. In some cases, and particularly for M stars, this sudden increase may last for several hours. We aim at answering the question to what extent a stellar flare is able to modify the chemical composition of the planetary atmosphere and, therefore influence the resulting spectra. We use a 1D thermo-photochemical model to study the neutral atmospheric composition of two hypothetic planets located around the star AD Leo. This active star has already been observed during a flare. We use the spectroscopic data from this flare event to simulate the evolution of the chemical composition of the atmospheres of the two hypothetic planets. We compute synthetic spectra to evaluate the implications for observations. The increase of the incoming photon flux affects the chemical abundances of some important species down to altitudes associated with an atmospheric pressure of 1 bar, that can lead to variations in planetary spectra if performed during transit.

New high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J=3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8+-0.5)x10^{-8} for ortho-NH3 and (3.2^{+0.7}_{-0.6})x10^{-8} for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1 sigma confidence level).

We investigate how night side cooling and surface friction impact surface temperatures and large scale circulation for tidally locked Earth-like planets. For each scenario, we vary the orbital period between $P_{rot}=1-100$~days and capture changes in climate states. We find drastic changes in climate states for different surface friction scenarios. For very efficient surface friction ($t_{s,fric}=$ 0.1 days), the simulations for short rotation periods ($P_{rot} \leq$ 10 days) show predominantly standing extra tropical Rossby waves. These waves lead to climate states with two high latitude westerly jets and unperturbed meridional direct circulation. In most other scenarios, simulations with short rotation periods exhibit instead dominance by standing tropical Rossby waves. Such climate states have a single equatorial westerly jet, which disrupts direct circulation. Experiments with weak surface friction ($t_{s,fric}=~10 -100$ days) show decoupling between surface temperatures and circulation, which leads to strong cooling of the night side. The experiment with $t_{s,fric}= 100$ days assumes climate states with easterly flow (retrograde rotation) for medium and slow planetary rotations $P_{rot}= 12 - 100$~days. We show that an increase of night side cooling efficiency by one order of magnitude compared to the nominal model leads to a cooling of the night side surface temperatures by 80-100~K. The day side surface temperatures only drop by 25~K at the same time. The increase in thermal forcing suppresses the formation of extra tropical Rossby waves on small planets ($R_P=1 R_{Earth}$) in the short rotation period regime ($P_{rot} \leq$ 10 days).

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.

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.

We investigate 3D atmosphere dynamics for tidally locked terrestrial planets with an Earth-like atmosphere and irradiation for different rotation periods ($P_{rot}=1-100$ days) and planet sizes ($R_P=1-2 R_{Earth}$) with unprecedented fine detail. We could precisely identify three climate state transition regions that are associated with phase transitions in standing tropical and extra tropical Rossby waves. We confirm that the climate on fast rotating planets may assume multiple states ($P_{rot}\leq 12$ days for $R_P=2 R_{Earth}$). Our study is, however, the first to identify the type of planetary wave associated with different climate states: The first state is dominated by standing tropical Rossby waves with fast equatorial superrotation. The second state is dominated by standing extra tropical Rossby waves with high latitude westerly jets with slower wind speeds. For very fast rotations ($P_{rot}\leq 5$~days for $R_P=2 R_{Earth}$), we find another climate state transition, where the standing tropical and extra tropical Rossby wave can both fit on the planet. Thus, a third state with a mixture of the two planetary waves becomes possible that exhibits three jets. Different climate states may be observable, because the upper atmosphere's hot spot is eastward shifted with respect to the substellar point in the first state, westward shifted in the second state and the third state shows a longitudinal 'smearing' of the spot across the substellar point. We show, furthermore, that the largest fast rotating planet in our study exhibits atmosphere features known from hot Jupiters like fast equatorial superrotation and a temperature chevron in the upper atmosphere.

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.

Aim: The late stages of stellar evolution are mainly governed by the mass of the stars. Low- and intermediate-mass stars lose copious amounts of mass during the asymptotic giant branch (AGB) which obscure the central star making it difficult to study the stellar spectra and determine the stellar mass. In this study, we present observational data that can be used to determine lower limits to the stellar mass. Method: Spectra of nine heavily reddened AGB stars taken by the Herschel Space Observatory display numerous molecular emission lines. The strongest emission lines are due to H2O. We search for the presence of isotopologues of H2O in these objects. Result: We detected the 16O and 17O isotopologues of water in these stars, but lines due to H2^{18}O are absent. The lack of 18O is predicted by a scenario where the star has undergone hot-bottom burning which preferentially destroys 18O relative to 16O and 17O. From stellar evolution calculations, this process is thought to occur when the stellar mass is above 5 Msun for solar metallicity. Hence, observations of different isotopologues of H2O can be used to help determine the lower limit to the initial stellar mass. Conclusion: From our observations, we deduce that these extreme OH/IR stars are intermediate-mass stars with masses of >= 5 Msun. Their high mass-loss rates of ~ 1.0e-4 Msun/yr may affect the enrichment of the interstellar medium and the overall chemical evolution of our Galaxy.

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.

Recent high-resolution observations have shown stellar winds to harbour complexities which strongly deviate from spherical symmetry, generally assumed as standard wind model. One such morphology is the archimedean spiral, generally believed to be formed by binary interactions, which has been directly observed in multiple sources. We seek to investigate the manifestation in the observables of spiral structures embedded in the spherical outflows of cool stars. We aim to provide an intuitive bedrock with which upcoming ALMA data can be compared and interpreted. By means of an extended parameter study, we model rotational CO emission from the stellar outflow of asymptotic giant branch stars. To this end, we develop a simplified analytical parametrised description of a 3D spiral structure. This model is embedded into a spherical wind, and fed into the 3D radiative transfer code LIME, which produces 3D intensity maps throughout velocity space. Subsequently, we investigate the spectral signature of rotational transitions of CO of the models, as well as the spatial aspect of this emission by means of wide-slit PV diagrams. Additionally, the potential for misinterpretation of the 3D data in a 1D context is quantified. Finally, we simulate ALMA observations to explore the impact of interefrometric noise and artifacts on the emission signatures. The spectral signatures of the CO rotational transition v=0 J=3-2 are very efficient at concealing the dual nature of the outflow. Only a select few parameter combinations allow for the spectral lines to disclose the presence of the spiral structure. The inability to disentangle the spiral from the spherical signal can result in an incorrect interpretation in a 1D context. Consequently, erroneous mass loss rates would be calculated...

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 discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Atmospheres with a high C/O ratio are expected to contain an important quantity of hydrocarbons, including heavy molecules (with more than 2 carbon atoms). To study correctly these C-rich atmospheres, a chemical scheme adapted to this composition is necessary. We have implemented a chemical scheme that can describe the kinetics of species with up to 6 carbon atoms. This chemical scheme has been developed with specialists of combustion and validated through experiments on a wide range of T and P. This chemical network is available on the online database KIDA. We have created a grid of 12 models to explore different thermal profiles and C/O ratios. For each of them, we have compared the chemical composition determined with a C0-C2 chemical scheme (species with up to 2 carbon atoms) and with the C0-C6 scheme. We found no difference in the results obtained with the two schemes when photolyses are not included in the model, whatever the temperature of the atmosphere. In contrast, when there is photochemistry, differences can appear in the upper atmosphere. These differences are found for all the tested PT profiles in the case that the C/O ratio is above 1. When the C/O ratio of the atmosphere is solar, differences are only found at temperatures lower than 1000K. The differences linked to the use of different chemical schemes do not have important influence on the synthetic spectra. However, we have confirmed that C2H2 and HCN as possible tracers of warm C-rich atmospheres. The use of this new chemical scheme is mandatory to model atmospheres with a high C/O ratio and, in particular, if one is interested in studying in details the photochemistry. If one is just interested in the synthetic spectra, the use of a smaller scheme may be sufficient.

Understanding the formation of planetary nebulae remains elusive because in the preceding asymtotic giant branch (AGB) phase these stars are heavily enshrouded in an optically thick dusty envelope. To further understand the morphology of the circumstellar environments of AGB stars we observe the closest carbon-rich AGB star IRC+10216 in scattered light. When imaged in scattered light at optical wavelengths, IRC+10216 surprisingly shows a narrow equatorial density enhancement, in contrast to the large-scale spherical rings that have been imaged much further out. We use radiative transfer models to interpret this structure in terms of two models: firstly, an equatorial density enhancement, commonly observed in the more evolved post-AGB stars, and secondly, in terms of a dust rings model, where a local enhancement of mass-loss creates a spiral ring as the star rotates. We conclude that both models can be used to reproduce the dark lane in the scattered light images, which is caused by an equatorially density enhancement formed by dense dust rather than a bipolar outflow as previously thought. We are unable to place constraints on the formation of the equatorial density enhancement by a binary system.

We developed a 1D photo-thermochemical model to study the atmosphere of warm exoplanets. The chemical scheme used in this model is completely new in planetology and has been constructed in collaboration with specialists of combustion. It has been validated as a whole through experiments on a large range of temperature (300 - 2500 K) and pressure (1 mbar - 100 bar), allowing to study a wide variety of exoplanets. We have used this chemical model to study the atmosphere of two warm Neptunes, GJ3470b and GJ436b, and the influence of different parameters (vertical mixing, metallicity, temperature, . . . ) on their chemical composition. We present here the results obtained in these studies.

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.

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

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.

(abbreviated) We aim to study the inner wind of the well-known AGB star CW Leo. Different diagnostics probing different geometrical scales have pointed toward a non-homogeneous mass-loss process: dust clumps are observed at milli-arcsec scale, a bipolar structure is seen at arcsecond-scale and multi-concentric shells are detected beyond 1". We present the first ALMA Cycle 0 band 9 data around 650 GHz. The full-resolution data have a spatial resolution of 0".42x0".24, allowing us to study the morpho-kinematical structure within ~6". Results: We have detected 25 molecular lines. The emission of all but one line is spatially resolved. The dust and molecular lines are centered around the continuum peak position. The dust emission has an asymmetric distribution with a central peak flux density of ~2 Jy. The molecular emission lines trace different regions in the wind acceleration region and suggest that the wind velocity increases rapidly from about 5 R* almost reaching the terminal velocity at ~11 R*. The channel maps for the brighter lines show a complex structure; specifically for the 13CO J=6-5 line different arcs are detected within the first few arcseconds. The curved structure present in the PV map of the 13CO J=6-5 line can be explained by a spiral structure in the inner wind, probably induced by a binary companion. From modeling the ALMA data, we deduce that the potential orbital axis for the binary system lies at a position angle of ~10-20 deg to the North-East and that the spiral structure is seen almost edge-on. We infer an orbital period of 55 yr and a binary separation of 25 au (or ~8.2 R*). We tentatively estimate that the companion is an unevolved low-mass main-sequence star. The ALMA data hence provide us for the first time with the crucial kinematical link between the dust clumps seen at milli-arcsecond scale and the almost concentric arcs seen at arcsecond scale.

Radiative equilibrium temperatures are calculated for the troposphere of a tidally locked Super-Earth based on a simple greenhouse model, using Solar System data as a guideline. These temperatures provide in combination with a Newtonian relaxation scheme thermal forcing for a 3D atmosphere model using the dynamical core of the Massachusetts Institute of Technology global circulation model (MITgcm). Our model is of the same conceptional simplicity than the model of Held & Suarez (1994) and is thus computationally fast. Furthermore, because of the coherent, general derivation of radiative equilibrium temperatures, our model is easily adaptable for different planets and atmospheric scenarios. As a case study relevant for Super-Earths, we investigate a Gl581g-like planet with Earth-like atmosphere and irradiation and present results for two representative rotation periods of $P_{rot}$ = 10 days and $P_{rot}$ = 36.5 days. Our results provide proof of concept and highlight interesting dynamical features for the rotating regime 3 < $P_{rot}$ < 100 days, which was shown by Edson et al. (2011) to be an intermediate regime between equatorial superrotation and divergence. We confirm that the $P_{rot}$ = 10 days case is more dominated by equatorial superrotation dynamics than the $P_{rot}$ = 36.5 days case, which shows diminishing influence of standing Rossby- Kelvin waves and increasing influence of divergence at the top of the atmosphere. We argue that this dynamical regime change relates to the increase in Rossby deformation radius, in agreement with previous studies. However, we also pay attention to other features that are not or only in partial agreement with other studies, like, e.g., the number of circulation cells and their strength, the role and extent of thermal inversion layers, and the details of heat transport

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.

We aim to examine the role of the interstellar magnetic field in shaping the extended morphologies of slow dusty winds of Asymptotic Giant-branch (AGB) stars in an effort to pin-point the origin of so-called eye shaped CSE of three carbon-rich AGB stars. In addition, we seek to understand if this pre-planetary nebula (PN) shaping can be responsible for asymmetries observed in PNe. Hydrodynamical simulations are used to study the effect of typical interstellar magnetic fields on the free-expanding spherical stellar winds as they sweep up the local interstellar medium (ISM). The simulations show that typical Galactic interstellar magnetic fields of 5 to 10 muG, are sufficient to alter the spherical expanding shells of AGB stars to appear as the characteristic eye shape revealed by far-infrared observations. The typical sizes of the simulated eyes are in accordance with the observed physical sizes. However, the eye shapes are of transient nature. Depending on the stellar and interstellar conditions they develop after 20,000 to 200,000yrs and last for about 50,000 to 500,000 yrs, assuming that the star is at rest relative to the local interstellar medium. Once formed the eye shape will develop lateral outflows parallel to the magnetic field. The "explosion" of a PN in the center of the eye-shaped dust shell gives rise to an asymmetrical nebula with prominent inward pointing Rayleigh-Taylor instabilities. Interstellar magnetic fields can clearly affect the shaping of wind-ISM interaction shells. The occurrence of the eyes is most strongly influenced by stellar space motion and ISM density. Observability of this transient phase is favoured for lines-of-sight perpendicular to the interstellar magnetic field direction. The simulations indicate that shaping of the pre-PN envelope can strongly affect the shape and size of PNe.

Many massive stars travel through the interstellar medium at supersonic speeds. As a result they form bow shocks at the interface between the stellar wind. We use numerical hydrodynamics to reproduce such bow shocks numerically, creating models that can be compared to observations. In this paper we discuss the influence of two physical phenomena, interstellar magnetic fields and the presence of interstellar dust grains on the observable shape of the bow shocks of massive stars. We find that the interstellar magnetic field, though too weak to restrict the general shape of the bow shock, reduces the size of the instabilities that would otherwise be observed in the bow shock of a red supergiant. The interstellar dust grains, due to their inertia can penetrate deep into the bow shock structure of a main sequence O-supergiant, crossing over from the ISM into the stellar wind. Therefore, the dust distribution may not always reflect the morphology of the gas. This is an important consideration for infrared observations, which are dominated by dust emission. Our models clearly show, that the bow shocks of massive stars are useful diagnostic tools that can used to investigate the properties of both the stellar wind as well as the interstellar medium.

Phosphine (PH3) was tentatively identified a few years ago in the carbon star envelopes IRC+10216 and CRL2688 from observations of an emission line at 266.9 GHz attributable to the J=1-0 rotational transition. We report the detection of the J=2-1 rotational transition of PH3 in IRC+10216 using the HIFI instrument on board Herschel, which definitively confirms the identification of PH3. Radiative transfer calculations indicate that infrared pumping to excited vibrational states plays an important role in the excitation of PH3 in the envelope of IRC+10216, and that the observed lines are consistent with phosphine being formed anywhere between the star and 100 R* from the star, with an abundance of 1e-8 relative to H2. The detection of PH3 challenges chemical models, none of which offers a satisfactory formation scenario. Although PH3 locks just 2 % of the total available phosphorus in IRC+10216, it is together with HCP, one of the major gas phase carriers of phosphorus in the inner circumstellar layers, suggesting that it could be also an important phosphorus species in other astronomical environments. This is the first unambiguous detection of PH3 outside the solar system, and a further step towards a better understanding of the chemistry of phosphorus in space.

For the Mass-loss of Evolved StarS (MESS) programme, the unprecedented spatial resolution of the PACS photometer on board the Herschel space observatory was employed to map the dusty environments of asymptotic giant branch (AGB) and red supergiant (RSG) stars. Among the morphologically heterogeneous sample, a small fraction of targets is enclosed by spherically symmetric detached envelopes. Based on observations in the 70 {\mu}m and 160 {\mu}m wavelength bands, we investigated the surroundings of the two carbon semiregular variables S Sct and RT Cap, which both show evidence for a history of highly variable mass-loss. S Sct exhibits a bright, spherically symmetric detached shell, 138" in diameter and co-spatial with an already known CO structure. Moreover, weak emission is detected at the outskirts, where the morphology seems indicative of a mild shaping by interaction of the wind with the interstellar medium, which is also supported by the stellar space motion. Two shells are found around RT Cap that were not known so far in either dust emission or from molecular line observations. The inner shell with a diameter of 188" shows an almost immaculate spherical symmetry, while the outer ~5' structure is more irregularly shaped. MoD, a modification of the DUSTY radiative transfer code, was used to model the detached shells. Dust temperatures, shell dust masses, and mass-loss rates are derived for both targets.

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

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.

Many evolved stars travel through space at supersonic velocities, which leads to the formation of bow shocks ahead of the star where the stellar wind collides with the interstellar medium (ISM). Herschel observations of the bow shock of $\alpha$-Orionis show that the shock is almost free of instabilities, despite being, at least in theory, subject to both Kelvin-Helmholtz and Rayleigh-Taylor instabilities. A possible explanation for the lack of instabilities lies in the presence of an interstellar magnetic field. We wish to investigate whether the magnetic field of the interstellar medium (ISM) in the Orion arm can inhibit the growth of instabilities in the bow shock of $\alpha$-Orionis. We used the code MPI-AMRVAC to make magneto-hydrodynamic simulations of a circumstellar bow shock, using the wind parameters derived for $\alpha$-Orionis and interstellar magnetic field strengths of $B\,=\,1.4,\, 3.0$, and $5.0\, \mu$G, which fall within the boundaries of the observed magnetic field strength in the Orion arm of the Milky Way. Our results show that even a relatively weak magnetic field in the interstellar medium can suppress the growth of Rayleigh-Taylor and Kelvin-Helmholtz instabilities, which occur along the contact discontinuity between the shocked wind and the shocked ISM. The presence of even a weak magnetic field in the ISM effectively inhibits the growth of instabilities in the bow shock. This may explain the absence of such instabilities in the Herschel observations of $\alpha$-Orionis.

We present an analysis of the far-infrared and submillimetre molecular emission line spectrum of the luminous M-supergiant VY CMa, observed with the SPIRE and PACS spectrometers aboard the Herschel Space Observatory. Over 260 emission lines were detected in the 190-650-micron SPIRE FTS spectra, with one-third of the observed lines being attributable to H2O. Other detected species include CO, 13CO, H2^18O, SiO, HCN, SO, SO2, CS, H2S, and NH3. Our model fits to the observed 12CO and 13CO line intensities yield a 12C/13C ratio of 5.6+-1.8, consistent with measurements of this ratio for other M supergiants, but significantly lower than previously estimated for VY CMa from observations of lower-J lines. The spectral line energy distribution for twenty SiO rotational lines shows two temperature components: a hot component at 1000 K, which we attribute to the stellar atmosphere and inner wind, plus a cooler ~200 K component, which we attribute to an origin in the outer circumstellar envelope. We fit the line fluxes of 12CO, 13CO, H2O and SiO, using the SMMOL non-LTE line transfer code, with a mass-loss rate of 1.85x10^-4 Msun yr^-1 between 9 R* and 350 R*. To fit the observed line fluxes of 12CO, 13CO, H2O and SiO with SMMOL non-LTE line radiative transfer code, along with a mass-loss rate of 1.85x10^-4 Msun yr^-1. To fit the high rotational lines of CO and H2O, the model required a rather flat temperature distribution inside the dust condensation radius, attributed to the high H2O opacity. Beyond the dust condensation radius the gas temperature is fitted best by an r^-0.5 radial dependence, consistent with the coolant lines becoming optically thin. Our H2O emission line fits are consistent with an ortho:para ratio of 3 in the outflow.

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 spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

Convection, pulsation and magnetic fields have all been suggested as mechanisms for the transport of mass and energy from the optical photosphere of red supergiants, out to the region where the stellar wind is launched. We imaged the red supergiant Betelgeuse at 0.06-0.18 arcsec resolution, using e-MERLIN at 5.5--6.0 GHz, with a sensitivity of ~0.01 mJy/beam. Most of the radio emission comes from within an ellipse (0.235x0.218) arcsec^2 (~5x the optical radius), with a flux density of 1.62 mJy, giving an average brightness temperature ~1250 K. This radio photosphere contains two hotspots of 0.53 and 0.79 mJy/beam, separated by 90 milli-arcsec, with brightness temperatures 5400+/-600 K and 3800+/-500 K. Similar hotspots, at more than double the distance from the photosphere of those seen in any other regime, were detected by the less-sensitive `old' MERLIN in 1992, 1995 and 1996 and many exceed the photospheric temperature of 3600 K. Such brightness temperatures are high enough to emanate from pockets of chromospheric plasma. Other possibilities include local shock heating, the convective dredge-up of hot material or exceptionally cool, low density regions, transparent down to the hottest layer at ~40 milliarcsec radius. We also detect an arc 0.2--0.3 arcsec to the SW, brightness temperature ~150 K, in a similar direction to extensions seen on both smaller and larger scales in the infra-red and in CO at mm wavelengths. These preliminary results will be followed by further e-MERLIN, VLA and ALMA observations to help resolve the problem of mass elevation from 1 to 10 R* in red supergiants.

We report the results of a sensitive search for the minor isotopologues of water, H2-17O and H2-18O, toward the carbon-rich AGB star IRC+10216 (a.k.a. CW Leonis) using the HIFI instrument on the Herschel Space Observatory. This search was motivated by the fact that any detection of isotopic enhancement in the H2-17O and H2-18O abundances would have strongly implicated CO photodissociation as the source of the atomic oxygen needed to produce water in a carbon-rich circumstellar envelope. Our observations place an upper limit of 1/470 on the H2-17O/H2-16O abundance ratio. Given the isotopic 17O/16O ratio of 1/840 inferred previously for the photosphere of IRC+10216, this result places an upper limit of a factor 1.8 on the extent of any isotope-selective enhancement of H2-17O in the circumstellar material, and provides an important constraint on any model that invokes CO photodissociation as the source of O for H2O production. In the context of the clumpy photodissociation model proposed previously for the origin of water in IRC+10216, our limit implies that 12C-16O (not 13C-16O or SiO) must be the dominant source of 16O for H2O production, and that the effects of self-shielding can only have reduced the 12C-16O photodissociation rate by at most a factor ~ 2.

We report the first detection of pure rotational transitions of TiO and TiO_2 at (sub-)millimeter wavelengths towards the red supergiant VY CMa. A rotational temperature, T_rot, of about 250 K was derived for TiO_2. Although T_rot was not well constrained for TiO, it is likely somewhat higher than that of TiO_2. The detection of the Ti oxides confirms that they are formed in the circumstellar envelopes of cool oxygen-rich stars and may be the "seeds" of inorganic-dust formation, but alternative explanations for our observation of TiO and TiO_2 in the cooler regions of the envelope cannot be ruled out at this time. The observations suggest that a significant fraction of the oxides is not converted to dust, but instead remains in the gas phase throughout the outflow.

The Mass loss of Evolved StarS (MESS) sample offers a selection of 78 Asymptotic Giant Branch (AGB) stars and Red Supergiants (RSGs) observed with the PACS photometer on-board Herschel at 70 and 160 {\mu}m. For most of these objects, the dusty AGB wind differs from spherical symmetry and the wind shape can be subdivided into four classes. In the present paper we concentrate on the influence of a companion on the morphology of the stellar wind. Literature was searched to find binaries in the MESS sample and these are subsequently linked to their wind-morphology class to assert that the binaries are not distributed equally among the classes. In the second part of the paper we concentrate on the circumstellar environment of the two prominent objects R Aqr and W Aql. Each shows a characteristic signature of a companion interaction with the stellar wind. For the symbiotic star R Aqr, PACS revealed two perfectly opposing arms which in part reflect the previously observed ring-shaped nebula in the optical. However, from the far-IR there is evidence that the emitting region is elliptical rather than circular. The outline of the wind of W Aql seems to follow a large Archimedean spiral formed by the orbit of the companion but also shows strong indications of an interaction with the interstellar medium. The nature of the companion of W Aql was investigated and the magnitude of the orbital period supports the size of the spiral outline.

Cool objects glow in the infrared. The gas and solid-state species that escape the stellar gravitational attraction of evolved late-type stars in the form of a stellar wind are cool, with temperatures typically $\la$1500\,K, and can be ideally studied in the infrared. These stellar winds create huge extended circumstellar envelopes with extents approaching $10^{19}$\,cm. In these envelopes, a complex kinematical, thermodynamical and chemical interplay determines the global and local structural parameters. Unraveling the wind acceleration mechanisms and deriving the complicated structure of the envelopes is important to understand the late stages of evolution of ~97% of stars in galaxies as our own Milky Way. That way, we can also assess the significant chemical enrichment of the interstellar medium by the mass loss of these evolved stars. The Herschel Space Observatory is uniquely placed to study evolved stars thanks to the excellent capabilities of the three infrared and sub-millimeter instruments on board: PACS, SPIRE and HIFI. In this review, I give an overview of a few important results obtained during the first two years of Herschel observations in the field of evolved low and intermediate mass stars, and I will show how the Herschel observations can solve some historical questions on these late stages of stellar evolution, but also add some new ones.

Context. The interaction between stellar winds and the interstellar medium (ISM) can create complex bow shocks. The photometers on board the Herschel Space Observatory are ideally suited to studying the morphologies of these bow shocks. Aims. We aim to study the circumstellar environment and wind-ISM interaction of the nearest red supergiant, Betelgeuse. Methods. Herschel PACS images at 70, 100, and 160 micron and SPIRE images at 250, 350, and 500 micron were obtained by scanning the region around Betelgeuse. These data were complemented with ultraviolet GALEX data, near-infrared WISE data, and radio 21 cm GALFA-HI data. The observational properties of the bow shock structure were deduced from the data and compared with hydrodynamical simulations. Results. The infrared Herschel images of the environment around Betelgeuse are spectacular, showing the occurrence of multiple arcs at 6-7 arcmin from the central target and the presence of a linear bar at 9 arcmin. Remarkably, no large-scale instabilities are seen in the outer arcs and linear bar. The dust temperature in the outer arcs varies between 40 and 140 K, with the linear bar having the same colour temperature as the arcs. The inner envelope shows clear evidence of a non-homogeneous clumpy structure (beyond 15 arcsec), probably related to the giant convection cells of the outer atmosphere. The non-homogeneous distribution of the material even persists until the collision with the ISM. A strong variation in brightness of the inner clumps at a radius of 2 arcmin suggests a drastic change in mean gas and dust density some 32 000 yr ago. Using hydrodynamical simulations, we try to explain the observed morphology of the bow shock around Betelgeuse. Conclusions: [abbreviated]

Some planetary systems harbour debris disks containing planetesimals such as asteroids and comets. Collisions between such bodies produce small dust particles, the spectral features of which reveal their composition and, hence, that of their parent bodies. A measurement of the composition of olivine crystals has been done for the protoplanetary disk HD100546 and for olivine crystals in the warm inner parts of planetary systems. The latter compares well with the iron-rich olivine in asteroids (x<0.29). In the cold outskirts of the Beta Pictoris system, an analogue to the young Solar System, olivine crystals were detected but their composition remained undetermined, leaving unknown how the composition of the bulk of Solar System cometary olivine grains compares with that of extrasolar comets. Here we report the detection of the 69-micrometre-wavelength band of olivine crystals in the spectrum of Beta Pictoris. Because the disk is optically thin, we can associate the crystals with an extrasolar proto-Kuiper belt a distance of 15-45 astronomical units from the star (one astronomical unit is the Sun-Earth distance), determine their magnesium-rich composition (x=0.01+/-0.001) and show that they make up 3.6+/-1.0 per cent of the total dust mass. These values are strikingly similar to those for the dust emitted by the most primitive comets in the Solar System, even though Beta Pictoris is more massive and more luminous and has a different planetary system architecture.

We present new Herschel/PACS images at 70, 100, and 160 micron of the well-known, nearby, carbon-rich asymptotic giant branch star IRC+10216 revealing multiple dust shells in its circumstellar envelope. For the first time, dust shells (or arcs) are detected until 320 arcsec. The almost spherical shells are non-concentric and have an angular extent between 40 deg and 200 deg. The shells have a typical width of 5 arcsec - 8 arcsec, and the shell separation varies in the range of 10 arcsec - 35 arcsec, corresponding to 500-1700 yr. Local density variations within one arc are visible. The shell/intershell density contrast is typically 4, and the arcs contain some 50% more dust mass than the smooth envelope. The observed (nested) arcs record the mass-loss history over the past 16 000 yr, but Rayleigh-Taylor and Kelvin-Helmholtz instabilities in the turbulent astropause and astrosheath will erase any signature of the mass-loss history for at least the first 200 000 yr of mass loss. Accounting for the bowshock structure, the envelope mass around IRC+10216 contains >2Msun of gas and dust mass. It is argued that the origin of the shells is related to non-isotropic mass-loss events and clumpy dust formation.

Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) are believed to be the high-mass counterparts of stars in the AGB and early post-AGB phases. We study the mass-loss in the post main-sequence evolution of massive stars, through the properties of their envelopes in the intermediate and warm gas layers. These are the regions where the acceleration of the gas takes place and the most recent mass-loss episodes can be seen. We used the HIFI instrument on-board the Herschel Space Observatory to observe sub-mm and FIR transitions of CO, water, and their isotopologues in a sample of two RSGs (NML Cyg and Betelgeuse) and two YHGs (IRC+10420 and AFGL 2343) stars. We present an inventory of the detected lines and analyse the information revealed by their spectral profiles. On the basis of the results presented in an earlier study, we model the CO and 13CO emission in IRC+10420 and compare it to a set of lines ranging from the mm, to the FIR. Red supergiants have stronger high-excitation lines than the YHGs, indicating that they harbour dense and hot inner shells contributing to these transitions. Consequently, these high-J lines in RSGs originate from acceleration layers that have not yet reached the circumstellar terminal velocity and have narrower profiles than their flat-topped lower-J counterparts. The YHGs tend to lack this inner component, in line with the picture of detached, hollow envelopes derived from studies at longer wavelengths. NH3 is only detected in two sources (NML Cyg, IRC+10420), which are also observed to be the strongest water-line emitters of the studied sample. In contrast, OH is detected in all sources and does not seem to correlate with the water line intensities. We show that the IRC+10420 model derived solely from mm low-J CO transitions is capable of reproducing the high-J transitions when the temperature in the inner shell is simply lowered by about 30%.

The broad 30 \mu m feature in carbon stars is commonly attributed to MgS dust particles. However, reproducing the 30 \mu m feature with homogeneous MgS grains would require much more sulfur relative to the solar abundance. Direct gas-phase condensation of MgS occurs at a low efficiency. Precipitation of MgS on SiC precursor grains provides a more efficient formation mechanism, such that the assumption of homogeneous MgS grains may not be correct. Using a Monte Carlo-based radiative transfer code, we aim to model the 30 \mu m feature of the extreme carbon star LL Peg with MgS dust particles. We find that for LL Peg this modeling is insensitive to the unknown MgS optical properties at \lambda < 10 \mu m. When MgS is allowed to be in thermal contact with amorphous carbon and SiC, the amount of MgS required to reproduce the strength of 30 \mu m feature agrees with the solar abundance of sulfur, thereby resolving the reported MgS mass problem. We conclude that MgS is a valid candidate to be the carrier of the 30 \mu m feature when it is part of a composite grain population that has optical properties representative of an ensemble of particle shapes.

CW Leo has been observed six times between October 2009 and June 2012 with the SPIRE instrument on board the Herschel satellite. Variability has been detected in the flux emitted by the central star with a period of 639 \pm 4 days, in good agreement with determinations in the literature. Variability is also detected in the bow shock around CW Leo that had previously been detected in the ultraviolet and Herschel PACS/SPIRE data. Although difficult to prove directly, our working hypothesis is that this variability is directly related to that of the central star. In this case, fitting a sine curve with the period fixed to 639 days results in a time-lag in the variability between bow shock and the central star of 402 \pm 37 days. The orientation of the bow shock relative to the plane of the sky is unknown (but see below). For an inclination angle of zero degrees, the observed time-lag translates into a distance to CW Leo of 130 \pm 13 pc, and for non-zero inclination angles the distance is smaller. Fitting the shape of the bow shock with an analytical model (Wilkin 1996), the effect of the inclination angle on the distance may be estimated. Making the additional assumption that the relative peculiar velocity between the interstellar medium (ISM) and CW Leo is determined entirely by the star space velocity with respect to the local standard of rest (i.e. a stationary ISM), the inclination angle is found to be (-33.3 \pm 0.8) degrees based on the observed proper motion and radial velocity. Using the Wilkin model, our current best estimate of the distance to CW Leo is 123 \pm 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 \pm 150 Lsol (internal error).

We report on the discovery of an infrared emission band in the Spitzer spectrum of the S-type AGB star NP Aurigae that is caused by TiO molecules in the circumstellar environment. We modelled the observed emission to derive the temperature of the TiO molecules (\approx 600 K), an upper limit on the column density (\approx 10^17.25 cm^{-2}) and a lower limit on the spatial extent of the layer that contains these molecules. (\approx 4.6 stellar radii). This is the first time that this TiO emission band is observed. A search for similar emission features in the sample of S-type stars yielded two additional candidates. However, owing to the additional dust emission, the identification is less stringent. By comparing the stellar characteristics of NP Aur to those of the other stars in our sample, we find that all stars with TiO emission show large-amplitude pulsations, s-process enrichment, and a low C/O ratio. These characteristics might be necessary requirements for a star to show TiO in emission, but they are not sufficient.

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.

S-type AGB stars are thought to be in the transitional phase between M-type and C-type AGB stars. Because of their peculiar chemical composition, one may expect a strong influence of the stellar C/O ratio on the molecular chemistry and the mineralogy of the circumstellar dust. In this paper, we present a large sample of 87 intrinsic galactic S-type AGB stars, observed at infrared wavelengths with the Spitzer Space Telescope, and supplemented with ground-based optical data. On the one hand, we derive the stellar parameters from the optical spectroscopy and photometry, using a grid of model atmospheres. On the other, we decompose the infrared spectra to quantify the flux-contributions from the different dust species. Finally, we compare the independently determined stellar parameters and dust properties. For the stars without significant dust emission, we detect a strict relation between the presence of SiS absorption in the Spitzer spectra and the C/O ratio of the stellar atmosphere. These absorption bands can thus be used as an additional diagnostic for the C/O ratio. For stars with significant dust emission, we define three groups, based on the relative contribution of certain dust species to the infrared flux. We find a strong link between group-membership and C/O ratio. We show that these groups can be explained by assuming that the dust-condensation can be cut short before silicates are produced, while the remaining free atoms and molecules can then form the observed magnesium sulfides or the carriers of the unidentified 13 and 20 micron features. Finally, we present the detection of emission features attributed to molecules and dust characteristic to C-type stars, such as molecular SiS, hydrocarbons and magnesium sulfide grains. We show that we often detect magnesium sulfides together with molecular SiS and we propose that it is formed by a reaction of SiS molecules with Mg.

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.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

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.

Far-infrared Herschel/PACS images at 70 and 160 micron of a sample of 78 Galactic evolved stars are used to study the (dust) emission structures, originating from stellar wind-ISM interaction. In addition, two-fluid hydrodynamical simulations of the coupled gas and dust in wind-ISM interactions are used to compare with the observations. Four distinct classes of wind-ISM interaction (i.e. "fermata", "eyes", "irregular", and "rings") are identified and basic parameters affecting the morphology are discussed. We detect bow shocks for ~40% of the sample and detached rings for ~20%. De-projected stand-off distances (R_0) -- defined as the distance between the central star and the nearest point of the interaction region -- of the detected bow shocks ("fermata" and "eyes") are derived from the PACS images and compared to previous results, model predictions and the simulations. All observed bow shocks have stand-off distances smaller than 1 pc. Observed and theoretical stand-off distances are used together to independently derive the local ISM density. Both theoretical (analytical) models and hydrodynamical simulations give stand-off distances for adopted stellar properties that are in good agreement with the measured de-projected stand-off distance of wind-ISM bow shocks. The possible detection of a bow shock -- for the distance limited sample -- appears to be governed by its physical size as set roughly by the stand-off distance. In particular the star's peculiar space velocity and the density of the ISM appear decisive in detecting emission from bow shocks or detached rings. Tentatively, the "eyes" class objects are associated to (visual) binaries, while the "rings" generally appear not to occur for M-type stars, only for C or S-type objects that have experienced a thermal pulse.

We present numerical simulations of the hydrodynamical interactions that produce circumstellar shells. These simulations include several scenarios, such as wind-wind interaction and wind-ISM collisions. In our calculations we have taken into account the presence of dust in the stellar wind. Our results show that, while small dust grains tend to be strongly coupled to the gas, large dust grains are only weakly coupled. As a result, the distribution of the large dust grains is not representative of the gas distribution. Combining these results with observations may give us a new way of validating hydrodynamical models of the circumstellar medium.

We report the detection of emission in the J=1-0 rotational transition of hydrogen fluoride (HF), together with observations of the J=1-0 to J=3-2 rotational lines of H35Cl and H37Cl, towards the envelope of the carbon star IRC +10216. High-sensitivity, high-spectral resolution observations have been carried out with the HIFI instrument on board Herschel, allowing us to resolve the line profiles and providing insights into the spatial distribution of the emission. Our interpretation of the observations, with the use of radiative transfer calculations, indicates that both HF and HCl are formed in the inner regions of the envelope close to the AGB star. Thermochemical equilibrium calculations predict HF and HCl to be the major reservoirs of fluorine and chlorine in the atmospheres of AGB stars. The abundances relative to H2 derived for HF and HCl, 8e-9 and 1e-7 respectively, are substantially lower than those predicted by thermochemical equilibrium, indicating that F and Cl are likely affected by significant depletion onto dust grains, although some chlorine may be in the form of atomic Cl. The H35Cl/H37Cl abundance ratio is 3.3 +/- 0.3. The low abundance derived for HF in IRC +10216 makes it likely that the fluorine abundance is not enhanced over the solar value by nucleosynthesis in the AGB star, although this conclusion may not be robust because the HF abundance we derive is a lower limit to the elemental abundance of F. These observations suggest that both HF and HCl should be detectable through low J rotational transitions in other evolved stars.

Herschel's PACS instrument observed the environment of the binary system Mira Ceti in the 70 and 160 micron bands. These images reveal bright structures shaped as five broken arcs and fainter filaments in the ejected material of Mira's primary star. The overall shape of the IR emission around Mira deviates significantly from the expected alignment with Mira's exceptionally high space velocity. The observed broken arcs are neither connected to each other nor are they of a circular shape; they stretch over angular ranges of 80 to 100 degrees. By comparing Herschel and GALEX data, we found evidence for the disruption of the IR arcs by the fast outflow visible in both Halpha and the far UV. Radial intensity profiles are derived, which place the arcs at distances of 6-85" (550 - 8000 AU) from the binary. Mira's IR environment appears to be shaped by the complex interaction of Mira's wind with its companion, the bipolar jet, and the ISM.

The asymptotic giant branch (AGB) stars X Her and TX Psc have been imaged at 70 and 160 microns with the PACS instrument onboard the Herschel satellite, as part of the large MESS (Mass loss of Evolved StarS) Guaranteed Time Key Program. The images reveal an axisymmetric extended structure with its axis oriented along the space motion of the stars. This extended structure is very likely to be shaped by the interaction of the wind ejected by the AGB star with the surrounding interstellar medium (ISM). As predicted by numerical simulations, the detailed structure of the wind-ISM interface depends upon the relative velocity between star+wind and the ISM, which is large for these two stars (108 and 55 km/s for X Her and TX Psc, respectively). In both cases, there is a compact blob upstream whose origin is not fully elucidated, but that could be the signature of some instability in the wind-ISM shock. Deconvolved images of X Her and TX Psc reveal several discrete structures along the outermost filaments, which could be Kelvin-Helmholtz vortices. Finally, TX Psc is surrounded by an almost circular ring (the signature of the termination shock?) that contrasts with the outer, more structured filaments. A similar inner circular structure seems to be present in X Her as well, albeit less clearly.

We study the hydrodynamical behavior occurring in the turbulent interaction zone of a fast moving red supergiant star, where the circumstellar and interstellar material collide. In this wind-interstellar medium collision, the familiar bow shock, contact discontinuity, and wind termination shock morphology forms, with localized instability development. Our model includes a detailed treatment of dust grains in the stellar wind, and takes into account the drag forces between dust and gas. The dust is treated as pressureless gas components binned per grainsize, for which we use ten representative grainsize bins. Our simulations allow to deduce how dust grains of varying sizes become distributed throughout the circumstellar medium. We show that smaller dust grains (radius <0.045 micro-meters) tend to be strongly bound to the gas and therefore follow the gas density distribution closely, with intricate finestructure due to essentially hydrodynamical instabilities at the wind-related contact discontinuity. Larger grains which are more resistant to drag forces are shown to have their own unique dust distribution, with progressive deviations from the gas morphology. Specifically, small dust grains stay entirely within the zone bound by shocked wind material. The large grains are capable of leaving the shocked wind layer, and can penetrate into the shocked or even unshocked interstellar medium. Depending on how the number of dust grains varies with grainsize, this should leave a clear imprint in infrared observations of bowshocks of red supergiants and other evolved stars.

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.

This paper presents spectra in the 2 to 20 micron range of quiescent cloud material located in the IC 5146 cloud complex. The spectra were obtained with NASA's Infrared Telescope Facility (IRTF) SpeX instrument and the Spitzer Space Telescope's Infrared Spectrometer. We use these spectra to investigate dust and ice absorption features in pristine regions of the cloud that are unaltered by embedded stars. We find that the H2O-ice threshold extinction is 4.03+/-0.05 mag. Once foreground extinction is taken into account, however, the threshold drops to 3.2 mag, equivalent to that found for the Taurus dark cloud, generally assumed to be the touchstone quiescent cloud against which all other dense cloud and embedded young stellar object observations are compared. Substructure in the trough of the silicate band for two sources is attributed to CH3OH and NH3 in the ices, present at the ~2% and ~5% levels, respectively, relative to H2O-ice. The correlation of the silicate feature with the E(J-K) color excess is found to follow a much shallower slope relative to lines of sight that probe diffuse clouds, supporting the previous results by Chiar et al. (2007).

The relation between ices in the envelopes and disks surrounding YSOs and those in the quiescent interstellar medium is investigated. For a sample of 31 stars behind isolated dense cores, ground-based and Spitzer spectra and photometry in the 1-25 um wavelength range are combined. The baseline for the broad and overlapping ice features is modeled, using calculated spectra of giants, H2O ice and silicates. The adopted extinction curve is derived empirically. Its high resolution allows for the separation of continuum and feature extinction. The extinction between 13-25 um is ~50% relative to that at 2.2 um. The strengths of the 6.0 and 6.85 um absorption bands are in line with those of YSOs. Thus, their carriers, which, besides H2O and CH3OH, may include NH4+, HCOOH, H2CO and NH3, are readily formed in the dense core phase, before stars form. The 3.53 um C-H stretching mode of solid CH3OH was discovered. The CH3OH/H2O abundance ratios of 5-12% are larger than upper limits in the Taurus molecular cloud. The initial ice composition, before star formation occurs, therefore depends on the environment. Signs of thermal and energetic processing that were found toward some YSOs are absent in the ices toward background stars. Finally, the peak optical depth of the 9.7 um band of silicates relative to the continuum extinction at 2.2 um is significantly shallower than in the diffuse interstellar medium. This extends the results of Chiar et al. (2007) to a larger sample and higher extinctions.

We report the preliminary results of a survey for water vapor in a sample of eight C stars with large mid-IR continuum fluxes: V384 Per, CIT 6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, and IRC+40540. This survey, performed using the HIFI instrument on board the Herschel Space Observatory, entailed observations of the lowest transitions of both ortho- and para-water: the 556.936 GHz 1(10)-1(01) and 1113.343 GHz 1(11)-0(00) transitions, respectively. Water vapor was unequivocally detected in all eight of the target stars. Prior to this survey, IRC+10216 was the only carbon-rich AGB star from which thermal water emissions had been discovered, in that case with the use of the Submillimeter Wave Astronomy Satellite (SWAS). Our results indicate that IRC+10216 is not unusual, except insofar as its proximity to Earth leads to a large line flux that was detectable with SWAS. The water spectral line widths are typically similar to those of CO rotational lines, arguing against the vaporization of a Kuiper belt analog (Ford & Neufeld 2001) being the general explanation for water vapor in carbon-rich AGB stars. There is no apparent correlation between the ratio of the integrated water line fluxes to the 6.3 micron continuum flux - a ratio which measures the water outflow rate - and the total mass-loss rate for the stars in our sample.

MESS (Mass-loss of Evolved StarS) is a Guaranteed Time Key Program that uses the PACS and SPIRE instruments on board the Herschel Space Observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. In total, of order 150 objects are observed in imaging and about 50 objects in spectroscopy. This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated using results obtained during Herschel's science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type stars.

We report the results of observations of ten rotational transitions of water vapor toward the carbon-rich AGB (asymptotic giant branch) star IRC+10216 (CW Leonis), carried out with Herschel's HIFI instrument. Each transition was securely detected by means of observations using the dual beam switch mode of HIFI. The measured line ratios imply that water vapor is present in the inner outflow at small distances (few x 1.E+14 cm) from the star, confirming recent results reported by Decin et al. from observations with Herschel's PACS and SPIRE instruments. This finding definitively rules out the hypothesis that the observed water results from the vaporization of small icy objects in circular orbits. The origin of water within the dense C-rich envelope of IRC+10216 remains poorly understood. We derive upper limits on the H2-17O/H2-16O and H2-18O/H2-16O isotopic abundance ratios of ~ 5.E-3 (3 sigma), providing additional constraints on models for the origin of the water vapor in IRC+10216.

Studying the composition of dust in the interstellar medium (ISM) is crucial in understanding the cycle of dust in our galaxy. The mid-infrared spectral signature of amorphous silicates, the most abundant dust species in the ISM, is studied in different lines-of-sight through the Galactic plane, thus probing different conditions in the ISM. We have analysed 10 spectra from the Spitzer archive, of which 6 lines-of-sight probe diffuse interstellar medium material and 4 probe molecular cloud material. The 9.7 um silicate absorption features in 7 of these spectra were studied in terms of their shape and strength. In addition, the shape of the 18 um silicate absorption features in 4 of the diffuse sightline spectra were analysed. The 9.7 um silicate absorption bands in the diffuse sightlines show a strikingly similar band shape. This is also the case for all but one of the 18 um silicate absorption bands observed in diffuse lines-of-sight. The 9.7 um bands in the 4 molecular sightlines show small variations in shape. These modest variations in the band shape are inconsistent with the interpretation of the large variations in {\tau}_9.7/E(J-K) between diffuse and molecular sightlines in terms of silicate grain growth. Instead, we suggest that the large changes in {\tau}_9.7 / E(J-K) must be due to changes in E(J-K).

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.

In this article we present the detection of the 69 {\mu}m band of the crystalline olivine forsterite within the MESS key program of Herschel. We determine the temperature of the forsterite grains by fitting the 69 {\mu}m band.

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.

We report on the initial analysis of Herschel/HIFI carbon monoxide (CO) observations of the Martian atmosphere performed between 11 and 16 April 2010. We selected the (7-6) rotational transitions of the isotopes ^{13}CO at 771 GHz and C^{18}O at 768 GHz in order to retrieve the mean vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO mixing ratio was determined as 980 \pm 150 ppm, in agreement with the 900 ppm derived from Herschel/SPIRE observations in November 2009.

We report on an initial analysis of Herschel/HIFI observations of hydrogen chloride (HCl), hydrogen peroxide (H_2O_2), and molecular oxygen (O_2) in the Martian atmosphere performed on 13 and 16 April 2010 (L_s ~ 77{\deg}). We derived a constant volume mixing ratio of 1400 +/- 120 ppm for O_2 and determined upper limits of 200 ppt for HCl and 2 ppb for H_2O_2. Radiative transfer model calculations indicate that the vertical profile of O_2 may not be constant. Photochemical models determine the lowest values of H_2O_2 to be around L_s ~ 75{\deg} but overestimate the volume mixing ratio compared to our measurements.

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.

A set of new, sensitive, and spectrally resolved, sub-millimeter line observations are used to probe the warm circumstellar gas around the S-type AGB star chi Cyg. The observed lines involve high rotational quantum numbers, which, combined with previously obtained lower-frequency data, make it possible to study in detail the chemical and physical properties of, essentially, the entire circumstellar envelope of chi Cyg. The data were obtained using the HIFI instrument aboard Herschel, whose high spectral resolution provides valuable information about the line profiles. Detailed, non-LTE, radiative transfer modelling, including dust radiative transfer coupled with a dynamical model, has been performed to derive the temperature, density, and velocity structure of the circumstellar envelope. We report the first detection of circumstellar H2O rotational emission lines in an S-star. Using the high-J CO lines to derive the parameters for the circumstellar envelope, we modelled both the ortho- and para-H2O lines. Our modelling results are consistent with the velocity structure expected for a dust-driven wind. The derived total H2O abundance (relative to H2) is (1.1 +/- 0.2)E-5, much lower than that in O-rich stars. The derived ortho-to-para ratio of 2.1 +/- 0.6 is close to the high-temperature equilibrium limit, consistent with H2O being formed in the photosphere.

We performed Herschel/HIFI observations of several CO lines in the far-infrared/sub-mm in the protoplanetary nebula CRL618. The high spectral resolution provided by HIFI allows measurement of the line profiles. Since the dynamics and structure of the nebula is well known from mm-wave interferometric maps, it is possible to identify the contributions of the different nebular components (fast bipolar outflows, double shells, compact slow shell) to the line profiles. The observation of these relatively high-energy transitions allows an accurate study of the excitation conditions in these components, particularly in the warm ones, which cannot be properly studied from the low-energy lines. The 12CO J=16-15, 10-9, and 6-5 lines are easily detected in this source. 13CO J=10-9 and 6-5 are also detected. Wide profiles showing spectacular line wings have been found, particularly in 12CO 16-15. Other lines observed simultaneously with CO are also shown. Our analysis of the CO high-J transitions, when compared with the existing models, confirms the very low expansion velocity of the central, dense component, which probably indicates that the shells ejected during the last AGB phases were driven by radiation pressure under a regime of maximum transfer of momentum. No contribution of the diffuse halo found from mm-wave data is identified in our spectra, because of its low temperature. We find that the fast bipolar outflow is quite hot, much hotter than previously estimated; for instance, gas flowing at 100 km/s must have a temperature higher than ~ 200 K. Probably, this very fast outflow, with a kinematic age < 100 yr, has been accelerated by a shock and has not yet cooled down. The double empty shell found from mm-wave mapping must also be relatively hot, in agreement with the previous estimate.

We report the discovery of water vapour toward the carbon star V Cygni. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1(11) - 0(00) para-water transition at 1113.3430 GHz in the upper sideband of the Band 4b receiver. The observed spectral line profile is nearly parabolic, but with a slight asymmetry associated with blueshifted absorption, and the integrated antenna temperature is 1.69 \pm 0.17 K km/s. This detection of thermal water vapour emission, carried out as part of a small survey of water in carbon-rich stars, is only the second such detection toward a carbon-rich AGB star, the first having been obtained by the Submillimeter Wave Astronomy Satellite toward IRC+10216. For an assumed ortho-to-para ratio of 3 for water, the observed line intensity implies a water outflow rate ~ (3 - 6) E-5 Earth masses per year and a water abundance relative to H2 of ~ (2-5) E-6. This value is a factor of at least 1E+4 larger than the expected photospheric abundance in a carbon-rich environment, and - as in IRC+10216 - raises the intriguing possibility that the observed water is produced by the vapourisation of orbiting comets or dwarf planets. However, observations of the single line observed to date do not permit us to place strong constraints upon the spatial distribution or origin of the observed water, but future observations of additional transitions will allow us to determine the inner radius of the H2O-emitting zone, and the H2O ortho-to-para ratio, and thereby to place important constraints upon the origin of the observed water emission.

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 high spectral resolution and sensitivity of Herschel/HIFI allows for the detection of multiple rotational water lines and accurate determinations of water production rates in comets. In this letter we present HIFI observations of the fundamental 110-101 (557 GHz) ortho and 111-000 (1113 GHz) para rotational transitions of water in comet 81P/Wild 2 acquired in February 2010. We mapped the extent of the water line emission with five point scans. Line profiles are computed using excitation models which include excitation by collisions with electrons and neutrals and solar infrared radiation. We derive a mean water production rate of $1.0 \times 10^{28}$ molecules s$^{-1}$ at a heliocentric distance of 1.61 AU about 20 days before perihelion, in agreement with production rates measured from the ground using observations of the 18-cm OH lines. Furthermore, we constrain the electron density profile and gas kinetic temperature, and estimate the coma expansion velocity by fitting the water line shapes.

We report on the initial analysis of a Herschel/PACS full range spectrum of Neptune, covering the 51-220 micrometer range with a mean resolving power of ~ 3000, and complemented by a dedicated observation of CH4 at 120 micrometers. Numerous spectral features due to HD (R(0) and R(1)), H2O, CH4, and CO are present, but so far no new species have been found. Our results indicate that (i) Neptune's mean thermal profile is warmer by ~ 3 K than inferred from the Voyager radio-occultation; (ii) the D/H mixing ratio is (4.5+/-1) X 10**-5, confirming the enrichment of Neptune in deuterium over the protosolar value (~ 2.1 X 10**-5); (iii) the CH4 mixing ratio in the mid stratosphere is (1.5+/-0.2) X 10**-3, and CH4 appears to decrease in the lower stratosphere at a rate consistent with local saturation, in agreement with the scenario of CH4 stratospheric injection from Neptune's warm south polar region; (iv) the H2O stratospheric column is (2.1+/-0.5) X 10**14 cm-2 but its vertical distribution is still to be determined, so the H2O external flux remains uncertain by over an order of magnitude; and (v) the CO stratospheric abundance is about twice the tropospheric value, confirming the dual origin of CO suspected from ground-based millimeter/submillimeter observations.

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 have obtained the first continuous disk averaged spectrum of Mars from 450 to 1550 Ghz using the Herschel-SPIRE Fourier Transform Spectrometer. The spectrum was obtained at a constant resolution of 1.4 GHz across the whole band. The flux from the planet is such that the instrument was operated in "bright source" mode to prevent saturation of the detectors. This was the first successful use of this mode and in this work we describe the method used for observing Mars together with a detailed discussion of the data reduction techniques required to calibrate the spectrum. We discuss the calibration accuracy obtained and describe the first comparison with surface and atmospheric models. In addition to a direct photometric measurement of the planet the spectrum contains the characteristic transitions of 12CO from J 5-4 to J 13-12 as well as numerous H2O transitions. Together these allow the comparison to global atmospheric models allowing the mean mixing ratios of water and 12CO to be investigated. We find that it is possible to match the observed depth of the absorption features in the spectrum with a fixed water mixing ratio of 1 x 10-4 and a 12CO mixing ratio of 9 x 10-4

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.

High-resolution far-infrared and sub-millimetre spectroscopy of water lines is an important tool to understand the physical and chemical properties of cometary atmospheres. We present observations of several rotational ortho- and para-water transitions in comet C/2008 Q3 (Garradd) performed with HIFI on Herschel. These observations have provided the first detection of the 2_{12}-1_{01} (1669 GHz) ortho and 1_{11}-0_{00} (1113 GHz) para transitions of water in a cometary spectrum. In addition, the ground-state transition 1_{10}-1_{01} at 557 GHz is detected and mapped. By detecting several water lines quasi-simultaneously and mapping their emission we can constrain the excitation parameters in the coma. Synthetic line profiles are computed using excitation models which include excitation by collisions, solar infrared radiation, and radiation trapping. We obtain the gas kinetic temperature, constrain the electron density profile, and estimate the coma expansion velocity by analyzing the map and line shapes. We derive water production rates of 1.7-2.8 x 10^{28} s^{-1} over the range r_h = 1.83-1.85 AU.

We present far-infrared and submillimetre spectra of three carbon-rich evolved objects, AFGL 2688, AFGL 618 and NGC 7027. The spectra were obtained with the SPIRE Fourier transform spectrometer on board the Herschel Space Observatory, and cover wavelengths from 195-670 um, a region of the electromagnetic spectrum hitherto difficult to study in detail. The far infrared spectra of these objects are rich and complex, and we measure over 150 lines in each object. Lines due to 18 different species are detected. We determine physical conditions from observations of the rotational lines of several molecules, and present initial large velocity gradient models for AFGL 618. We detect water in AFGL 2688 for the first time, and confirm its presence in AFGL 618 in both ortho and para forms. In addition, we report the detection of the J=1-0 line of CH+ in NGC 7027.

With a luminosity > 10^5 Lsun and a mass-loss rate of about 2.10-4 Msun/yr, the red supergiant VY CMa truly is a spectacular object. Because of its extreme evolutionary state, it could explode as supernova any time. Studying its circumstellar material, into which the supernova blast will run, provides interesting constraints on supernova explosions and on the rich chemistry taking place in such complex circumstellar envelopes. We have obtained spectroscopy of VYCMa over the full wavelength range offered by the PACS and SPIRE instruments of Herschel, i.e. 55 to 672 micron. The observations show the spectral fingerprints of more than 900 spectral lines, of which more than half belong to water. In total, we have identified 13 different molecules and some of their isotopologues. A first analysis shows that water is abundantly present, with an ortho-to-para ratio as low as 1.3:1, and that chemical non-equilibrium processes determine the abundance fractions in the inner envelope.

Detached circumstellar dust shells are detected around three carbon variables using Herschel-PACS. Two of them are already known on the basis of their thermal CO emission and two are visible as extensions in IRAS imaging data. By model fits to the new data sets, physical sizes, expansion timescales, dust temperatures, and more are deduced. A comparison with existing molecular CO material shows a high degree of correlation for TT Cyg and U Ant but a few distinct differences with other observables are also found.

Herschel PACS and SPIRE images have been obtained over a 30'x30' area around the well-known carbon star CW Leo (IRC +10 216). An extended structure is found in an incomplete arc of ~22' diameter, which is cospatial with the termination shock due to interaction with the interstellar medium (ISM) as defined by Sahai & Chronopoulos from ultraviolet GALEX images. Fluxes are derived in the 70, 160, 250, 350, and 550 um bands in the region where the interaction with the ISM takes place, and this can be fitted with a modified black body with a temperature of 25+-3 K. Using the published proper motion and radial velocity for the star, we derive a heliocentric space motion of 25.1 km/s. Using the PACS and SPIRE data and the analytical formula of the bow shock structure, we infer a de-projected standoff distance of the bow shock of R0 = (8.0+-0.3)x10^17 cm. We also derive a relative velocity of the star with respect to the ISM of (106.6+-8.7)/sqrt(n_ISM) km/s, where n_ISM is the number density of the local ISM.

Comet C/2006 W3 (Christensen) was observed in November 2009 at 3.3 AU from the Sun with Herschel. The PACS instrument acquired images of the dust coma in 70- and 160-micrometers filters, and spectra covering several H2O rotational lines. Spectra in the range 450-1550 GHz were acquired with SPIRE. The comet emission continuum from 70 to 672 micrometers was measured, but no lines were detected. The spectral energy distribution indicates thermal emission from large particles and provides a measure of the size distribution index and dust production rate. The upper limit to the water production rate is compared to the production rates of other species (CO, CH3OH, HCN, H2S, OH) measured with the IRAM 30-m and Nancay telescopes. The coma is found to be strongly enriched in species more volatile than water, in comparison to comets observed closer to the Sun. The CO to H2O production rate ratio exceeds 220%. The dust to gas production rate ratio is on the order of 1.

The aim of this paper is to study the molecular composition in the circumstellar envelope around the oxygen-rich star IK Tau. We observed IK Tau in several (sub)millimeter bands using the APEX telescope during three observing periods. To determine the spatial distribution of the $\mathrm{^{12}CO(3-2)}$ emission, mapping observations were performed. To constrain the physical conditions in the circumstellar envelope, multiple rotational CO emission lines were modeled using a non local thermodynamic equilibrium radiative transfer code. The rotational temperatures and the abundances of the other molecules were obtained assuming local thermodynamic equilibrium. An oxygen-rich Asymptotic Giant Branch star has been surveyed in the submillimeter wavelength range. Thirty four transitions of twelve molecular species, including maser lines, were detected. The kinetic temperature of the envelope was determined and the molecular abundance fractions of the molecules were estimated. The deduced molecular abundances were compared with observations and modeling from the literature and agree within a factor of 10, except for SO$_2$, which is found to be almost a factor 100 stronger than predicted by chemical models. From this study, we found that IK Tau is a good laboratory to study the conditions in circumstellar envelopes around oxygen-rich stars with (sub)millimeter-wavelength molecular lines. We could also expect from this study that the molecules in the circumstellar envelope can be explained more faithful by non-LTE analysis with lower and higher transition lines than by simple LTE analysis with only lower transition lines. In particular, the observed CO line profiles could be well reproduced by a simple expanding envelope model with a power law structure.

As a result of the variability survey in Chi Persei and NGC6910, the number of Beta Cep stars that are members of these two open clusters is increased to twenty stars, nine in NGC6910 and eleven in Chi Persei. We compare pulsational properties, in particular the frequency spectra, of Beta Cep stars in both clusters and explain the differences in terms of the global parameters of the clusters. We also indicate that the more complicated pattern of the variability among B type stars in Chi Persei is very likely caused by higher rotational velocities of stars in this cluster. We conclude that the sample of pulsating stars in the two open clusters constitutes a very good starting point for the ensemble asteroseismology of Beta Cep-type stars and maybe also for other B-type pulsators.

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]

Polycyclic aromatic hydrocarbons (PAHs) produce characteristic infrared emission bands that have been observed in a wide range of astrophysical environments, where carbonaceous material is subjected to ultraviolet (UV) radiation. Although PAHs are expected to form in carbon-rich AGB stars, they have up to now only been observed in binary systems where a hot companion provides a hard radiation field. In this letter, we present low-resolution infrared spectra of four S-type AGB stars, selected from a sample of 90 S-type AGB stars observed with the infrared spectrograph aboard the Spitzer satellite. The spectra of these four stars show the typical infrared features of PAH molecules. We confirm the correlation between the temperature of the central star and the centroid wavelength of the 7.9 {\mu}m feature, present in a wide variety of stars spanning a temperature range from 3 000 to 12 000 K. Three of four sources presented in this paper extend this relation towards lower temperatures. We argue that the mixture of hydrocarbons we see in these S-stars has a rich aliphatic component. The fourth star, BZ CMa, deviates from this correlation. Based on the similarity with the evolved binary TU Tau, we predict that BZ CMa has a hot companion as well.

We aim to determine the fundamental parameters of a sample of B stars with apparent visual magnitudes below 8 in the field-of-view of the CoRoT space mission, from high-resolution spectroscopy. We developed an automatic procedure for the spectroscopic analysis of B-type stars with winds, based on an extensive grid of FASTWIND model atmospheres. We use the equivalent widths and/or the line profile shapes of continuum normalized hydrogen, helium and silicon line profiles to determine the fundamental properties of these stars in an automated way. After thorough tests, both on synthetic datasets and on very high-quality, high-resolution spectra of B stars for which we already had accurate values of their physical properties from alternative analyses, we applied our method to 66 B-type stars contained in the ground-based archive of the CoRoT space mission. We discuss the statistical properties of the sample and compare them with those predicted by evolutionary models of B stars. Our spectroscopic results provide a valuable starting point for any future seismic modelling of the stars, should they be observed by CoRoT.

HD172189 is a spectroscopic eclipsing binary system with a rapidly-rotating pulsating delta Scuti component. It is also a member of the open cluster IC4756. These combined characteristics make it an excellent laboratory for asteroseismic studies. To date, HD172189 has been analysed in detail photometrically but not spectroscopically. For this reason we have compiled a set of spectroscopic data to determine the absolute and atmospheric parameters of the components. We determined the radial velocities (RV) of both components using four different techniques. We disentangled the binary spectra using KOREL, and performed the first abundance analysis on both disentangled spectra. By combining the spectroscopic results and the photometric data, we obtained the component masses, 1.8 and 1.7 solar masses, and radii, 4.0 and 2.4 solar radii, for inclination i = 73.2 degrees, eccentricity e = 0.28, and orbital period 5.70198 days. Effective temperatures of 7600 K and 8100 K were also determined. The measured vsini are 78 and 74 km/s, respectively, giving rotational periods of 2.50 and 1.55 days for the components. The abundance analysis shows [Fe/H] = -0.28 for the primary (pulsating) star, consistent with observations of IC4756. We also present an assessment of the different analysis techniques used to obtain the RVs and the global parameters.

We present IRS spectra and revised MIPS photometry for the 18 members of the Eta Chamaeleontis cluster. Aged 8 Myr, the Eta Cha cluster is one of the few nearby regions within the 5-10 Myr age range, during which the disk fraction decreases dramatically and giant planet formation must come to an end. For the 15 low-mass members, we measure a disk fraction ~50%, high for their 8 Myr age, and 4 of the 8 disks lack near-IR excesses, consistent with the empirical definition of "transition'' disks. Most of the disks are comparable to geometrically flat disks. The comparison with regions of different ages suggests that at least some of the "transition" disks may represent the normal type of disk around low-mass stars. Therefore, their flattened structure and inner holes may be related to other factors (initial masses of the disk and the star, environment, binarity), rather than to pure time evolution. We analyze the silicate dust in the disk atmosphere, finding moderate crystalline fractions (~10-30%) and typical grain sizes ~1-3 micron, without any characteristic trend in the composition. These results are common to other regions of different ages, suggesting that the initial grain processing occurs very early in the disk lifetime (<1 Myr). Large grain sizes in the disk atmosphere cannot be used as a proxy for age, but are likely related to higher disk turbulence. The dust mineralogy varies between the 8-12micron and the 20-30 micron features, suggesting high temperature dust processing and little radial mixing. Finally, the analysis of IR and optical data on the B9 star Eta Cha reveals that it is probably surrounded by a young debris disk with a large inner hole, instead of being a classical Be star.

Context: Among late-type red giants, an interesting change occurs in the structure of the outer atmospheric layers as one moves to later spectral types in the Hertzsprung-Russell diagram: a chromosphere is always present, but the coronal emission diminishes and a cool massive wind steps in. Aims: Where most studies have focussed on short-wavelength observations, this article explores the influence of the chromosphere and the wind on long-wavelength photometric measurements. Methods: The observational spectral energy distributions are compared with the theoretical predictions of the MARCS atmosphere models for a sample of 9 K- and M-giants. The discrepancies found are explained using basic models for flux emission originating from a chromosphere or an ionized wind. Results: For 7 out of 9 sample stars, a clear flux excess is detected at (sub)millimeter and/or centimeter wavelengths. The precise start of the excess depends upon the star under consideration. The flux at wavelengths shorter than about 1 mm is most likely dominated by an optically thick chromosphere, where an optically thick ionized wind is the main flux contributor at longer wavelengths. Conclusions: Although the optical to mid-infrared spectrum of the studied K- and M-giants is well represented by a radiative equilibrium atmospheric model, the presence of a chromosphere and/or ionized stellar wind at higher altitudes dominates the spectrum in the (sub)millimeter and centimeter wavelength ranges. The presence of a flux excess also has implications on the role of these stars as fiducial spectrophotometric calibrators in the (sub)millimeter and centimeter wavelength range.

Gamma Doradus stars are excellent targets for asteroseismology since the gravity modes present in these stars probe the deep stellar interiors. Mode identification will improve the knowledge of these stars considerably. A selected group of Gamma Doradus stars and some candidates were observed with the Mercator telescope to find and/or confirm the periodicities in the light variations and to derive reliable amplitude ratios in different pass bands. A frequency analysis was performed on all new data obtained in the Geneva photometric system. In order to have more reliable and accurate frequencies, the new data were combined with similar data from the literature and with Hipparcos observations. A set of frequencies that minimized the the residuals in a harmonic fit was searched for while allowing means and amplitudes to vary from one observation set to another. Frequencies and amplitudes in the photometric passbands of the Geneva system are given for 21 Gamma Doradus stars. We report the discovery of HD 74504 as a newly found Gamma Doradus star. Our study provides the first extensive multicolour database for the understanding of gravity modes in F-type stars.

Context: Red super-giant (RSG) stars exhibit significant mass loss through a slow and dense wind. They are often considered to be the more massive counter parts of Asymptotic Giant Branch (AGB) stars. While the AGB mass loss is linked to their strong pulsations, the RSG are often only weakly variable. Aim: To study the conditions at the base of the wind, by determining the dust composition in a sample of RSG. The dust composition is thought to be sensitive to the density, temperature and acceleration at the base of the wind. Method: We compile a sample of 27 RSG infrared spectra (ISO-SWS) and supplement these with photometric measurements to obtain the full spectral energy distribution (SED). These data are modelled using a dust radiative transfer code. The results are scrutinised for correlations. Results: We find (1) strong correlations between dust composition, mass-loss rate and stellar luminosity, roughly in agreement with the theoretical dust condensation sequence, (2) the need for a continuous (near-)IR dust opacity and tentatively propose amorphous carbon, and (3) significant differences with AGB star winds: presence of PAHs, absence of 'the' 13 micron band, and a lack of strong water bands. Conclusions: Dust condensation in RSG is found to experience a similar freeze-out process as in AGB stars. Together with the positive effect of the stellar luminosity on the mass-loss rate, this suggests that radiation pressure on dust grains is an important ingredient in the driving mechanism. Still, differences with AGB stars are manifold and thus the winds of RSG deserve separate studies.

We present spectral energy distribution modelling of 6875 stars in omega Centauri, obtaining stellar luminosities and temperatures by fitting literature photometry to state-of-the-art MARCS stellar models. By comparison to four different sets of isochrones, we provide a new distance estimate to the cluster of 4850 +/- 200 (random) +/- 120 (systematic error) pc, a reddening of E(B-V) = 0.08 +/- 0.02 +/- 0.02 mag and a differential reddening of Delta[E(B-V)] < 0.02 mag for an age of 12 Gyr. Several new post-early-AGB candidates are also found. Infra-red excesses of stars were used to measure total mass-loss rates for individual stars down to ~7 x 10^-8 Msun/yr. We find a total dust mass-loss rate from the cluster of 1.3 (+0.8/-0.5) x 10^-9 Msun/yr, with the total gas mass-loss rate being > 1.2 (+0.6/-0.5) x 10^-6 Msun/yr. Half of the cluster's dust production and 30% of its gas production comes from the two most extreme stars - V6 and V42 - for which we present new Gemini/T-ReCS mid-infrared spectroscopy, possibly showing that V42 has carbon-rich dust. The cluster's dust temperatures are found to be typically >~550 K. Mass loss apparently does not vary significantly with metallicity within the cluster, but shows some correlation with barium enhancement, which appears to occur in cooler stars, and especially on the anomalous RGB. Limits to outflow velocities, dust-to-gas ratios for the dusty objects and the possibility of short-timescale mass-loss variability are also discussed in the context of mass loss from low-metallicity stars. The ubiquity of dust around stars near the RGB-tip suggests significant dusty mass loss on the RGB; we estimate that typically 0.20--0.25 Msun of mass loss occurs on the RGB. From observational limits on intra-cluster material, we suggest the dust is being cleared on a timescale of <~10^5 years.

We determine an absolute calibration for the MIPS 24 microns band and recommend adjustments to the published calibrations for 2MASS, IRAC, and IRAS photometry to put them on the same scale. We show that consistent results are obtained by basing the calibration on either an average A0V star spectral energy distribution (SED), or by using the absolutely calibrated SED of the sun in comparison with solar-type stellar photometry (the solar analog method). After the rejection of a small number of stars with anomalous SEDs (or bad measurements), upper limits of ~ 1.5% (rms) are placed on the intrinsic infrared SED variations in both A dwarf and solar-type stars. These types of stars are therefore suitable as general-purpose standard stars in the infrared. We provide absolutely calibrated SEDs for a standard zero magnitude A star and for the sun to allow extending this work to any other infrared photometric system. They allow the recommended calibration to be applied from 1 to 25 microns with an accuracy of ~2 %, and with even higher accuracy at specific wavelengths such as 2.2, 10.6, and 24 microns, near which there are direct measurements. However, we confirm earlier indications that Vega does not behave as a typical A0V star between the visible and the infrared, making it problematic as the defining star for photometric systems. The integration of measurements of the sun with those of solar-type stars also provides an accurate estimate of the solar SED from 1 through 30 microns, which we show agrees with theoretical models.

Asymptotic Giant Branch (AGB) stars are typified by strong dust-driven, molecular outflows. For long, it was believed that the molecular setup of the circumstellar envelope of AGB stars is primarily determined by the atmospheric C/O ratio. However, recent observations of molecules such as HCN, SiO, and SO reveal gas-phase abundances higher than predicted by thermodynamic equilibrium (TE) models. UV-photon initiated dissociation in the outer envelope or non-equilibrium formation by the effect of shocks in the inner envelope may be the origin of the anomolous abundances. We aim at detecting (i.) a group of `parent' molecules (CO, SiO, HCN, CS), predicted by the non-equilibrium study of Cherchneff (2006) to form with almost constant abundances independent of the C/O ratio and the stellar evolutionary stage on the Asymptotic Giant Branch (AGB), and (ii.) few molecules, such as SiS and SO,which are sensitive to the O- or C-rich nature of the star. Several low and high excitation rotational transitions of key molecules are observed at mm and sub-mm wavelengths with JCMT and APEX in four AGB stars: the oxygen-rich Mira WX Psc, the S star W Aql, and the two carbon stars V Cyg and II Lup. A critical density analysis is performed to determine the formation region of the high-excitation molecular lines. We detect the four `parent' molecules in all four objects, implying that, indeed, these chemical species form whatever the stage of evolution on the AGB. High-excitation lines of SiS are also detected in three stars with APEX, whereas SO is only detected in the oxygen-rich star WX Psc. This is the first multi-molecular observational proof that periodically shocked layers above the photosphere of AGB stars show some chemical homogeneity, whatever the photospheric C/O ratio and stage of evolution of the star.

One of the key ingredients in establishing the relation between input signal and output flux from a spectrometer is accurate determination of the spectrophotometric calibration. In the case of spectrometers onboard satellites, the accuracy of this part of the calibration pedigree is ultimately linked to the accuracy of the set of reference SEDs that the spectrophotometric calibration is built on. In this paper, we deal with the spectrophotometric calibration of infrared (IR) spectrometers onboard satellites in the 2 to 200 micron range. We aim at comparing the different reference SEDs used for the IR spectrophotometric calibration. The emphasis is on the reference SEDs of stellar standards with spectral type later than A0, with special focus on the theoretical model atmosphere spectra. Using the MARCS model atmosphere code, spectral reference SEDs were constructed for a set of IR stellar standards (A dwarfs, solar analogs, G9-M0 giants). A detailed error analysis was performed to estimate proper uncertainties on the predicted flux values. It is shown that the uncertainty on the predicted fluxes can be as high as 10%, but in case high-resolution observational optical or near-IR data are available, and IR excess can be excluded, the uncertainty on medium-resolution SEDs can be reduced to 1-2% in the near-IR, to ~3% in the mid-IR, and to ~5% in the far-IR. Moreover, it is argued that theoretical stellar atmosphere spectra are at the moment the best representations for the IR fluxes of cool stellar standards. When aiming at a determination of the spectrophotometric calibration of IR spectrometers better than 3%, effort should be put into constructing an appropriate set of stellar reference SEDs based on theoretical atmosphere spectra for some 15 standard stars with spectral types between A0V and M0III.

Low and intermediate mass stars lose a significant fraction of their mass through a dust-driven wind during the Asymptotic Giant Branch (AGB) phase. Recent studies show that winds from late-type stars are far from being smooth. Mass-loss variations occur on different time scales, from years to tens of thousands of years. The variations appear to be particularly prominent towards the end of the AGB evolution. The occurrence, amplitude and time scale of these variations are still not well understood. The goal of our study is to gain insight into the structure of the circumstellar envelope (CSE) of WX Psc and map the possible variability of the late-AGB mass-loss phenomenon. We have performed an in-depth analysis of the extreme infrared AGB star WX Psc by modeling (1) the CO J=1-0 through 7-6 rotational line profiles and the full spectral energy distribution (SED) ranging from 0.7 to 1300 micron. We hence are able to trace a geometrically extended region of the CSE. Both mass-loss diagnostics bear evidence of the occurrence of mass-loss modulations during the last ~2000 yr. In particular, WX Psc went through a high mass-loss phase (Mdot~5e-5 Msun/yr) some 800 yr ago. This phase lasted about 600 yr and was followed by a long period of low mass loss (Mdot~5e-8 Msun/yr). The present day mass-loss rate is estimated to be ~6e-6 Msun/yr. The AGB star WX Psc has undergone strong mass-loss rate variability on a time scale of several hundred years during the last few thousand years. These variations are traced in the strength and profile of the CO rotational lines and in the SED. We have consistently simulated the behaviour of both tracers using radiative transfer codes that allow for non-constant mass-loss rates.

It is generally acknowledged that the mass loss of Asymptotic Giant Branch (AGB) stars undergoes variations on different time scales. We constructed models for the dust envelopes for a sample of AGB stars to assess whether mass-loss variations influence the spectral energy distribution (SED). To constrain the variability, extra observations at millimetre wavelengths (1.2 mm) were acquired. From the analysis of the dust models, two indications for the presence of mass-loss variations can be found, being (1) a dust temperature at the inner boundary of the dust envelope that is far below the dust condensation temperature and (2) an altered density distribution with respect to rho(r) \propto r^{-2} resulting from a constant mass-loss rate. For 5 out of the 18 studied sources a two-component model of the envelope is required, consisting of an inner region with a constant mass-loss rate and an outer region with a less steep density distribution. For one source an outer region with a steeper density distribution was found. Moreover, in a search for time variability in our data set at 1.2 mm, we found that WX Psc shows a large relative time variation of 34% which might partially be caused by variable molecular line emission.

A method is developed for fitting theoretically predicted astronomical spectra to an observed spectrum. Using a hierarchical Bayesian principle, the method takes both systematic and statistical measurement errors into account, which has not been done before in the astronomical literature. The goal is to estimate fundamental stellar parameters and their associated uncertainties. The non-availability of a convenient deterministic relation between stellar parameters and the observed spectrum, combined with the computational complexities this entails, necessitate the curtailment of the continuous Bayesian model to a reduced model based on a grid of synthetic spectra. A criterion for model selection based on the so-called predictive squared error loss function is proposed, together with a measure for the goodness-of-fit between observed and synthetic spectra. The proposed method is applied to the infrared 2.38--2.60 \mic ISO-SWS data (Infrared Space Observatory - Short Wavelength Spectrometer) of the star $\alpha$ Bootis, yielding estimates for the stellar parameters: effective temperature \Teff = 4230 $\pm$ 83 K, gravity $\log$ g = 1.50 $\pm$ 0.15 dex, and metallicity [Fe/H] = $-0.30 \pm 0.21$ dex.

Context. Abundance analysis of post-AGB objects as probes of AGB nucleosynthesis. Aims. A detailed photospheric abundance study is performed on the carbon-rich post-AGB candidate MACHO47.2496.8 in the LMC. Methods. High-resolution, high signal-to-noise ESO VLT-UVES spectra of MACHO47.2496.8 are analysed by performing detailed spectrum synthesis modelling using state-of-the-art carbon-rich MARCS atmosphere models. Results. The spectrum of MACHO47.2496.8 is not only dominated by bands of carbon bearing molecules, but also by lines of atomic transitions of s-process elements. The metallicity of [Fe/H]=-1.4 is surprisingly low for a field LMC star. The C/O ratio, however difficult to quantify, is greater than 2, and the s-process enrichment is large: the light s-process elements are enhanced by 1.2 dex compared to iron ([ls/Fe]=+1.2), while for the heavy s-process elements an even stronger enrichment is measured: [hs/Fe]=+2.1. The lead abundance is comparable to the [hs/Fe]. With its low intrinsic metallicity and its luminosity at the low end of the carbon star luminosity function, the star represents likely the final stage of a low initial mass star. Conclusions. The LMC RV Tauri star MACHO47.2496.8 is highly carbon and s-process enriched, and is most probable a genuine post-C(N-type) AGB star. This is the first detailed abundance analysis of an extragalactic post-AGB star to date.

Context: Mass loss plays a dominant role in the evolution of low mass stars while they are on the Asymptotic Giant Branch (AGB). The gas and dust ejected during this phase are a major source in the mass budget of the interstellar medium. Recent studies have pointed towards the importance of variations in the mass-loss history of such objects. Aims: By modelling the full line profile of low excitation CO lines emitted in the circumstellar envelope, we can study the mass-loss history of AGB stars. Methods: We have developed a non-LTE radiative transfer code, which calculates the velocity structure and gas kinetic temperature of the envelope in a self-consistent way. The resulting structure of the envelope provides the input for the molecular line radiative calculations which are evaluated in the comoving frame. The code allows for the implementation of modulations in the mass-loss rate. This code has been benchmarked against other radiative transfer codes and is shown to perform well and efficiently. Results: We illustrate the effects of varying mass-loss rates in case of a superwind phase. The model is applied to the well-studied case of VY CMa. We show that both the observed integrated line strengths as the spectral structure present in the observed line profiles, unambiguously demonstrate that this source underwent a phase of high mass loss (~ 3.2E-4 Msun/yr) some 1000 yr ago. This phase took place for some 100 yr, and was preceded by a low mass-loss phase (~ 1E-6 Msun/yr) taking some 800 yr. The current mass-loss rate is estimated to be in the order of 8E-5 Msun/yr. Conclusions: In this paper, we demonstrate that both the relative strength of the CO rotational line profiles and the (non)-occurrence of spectral structure in the profile offer strong diagnostics to pinpoint the mass-loss history.

In this paper we study the extended atmosphere of the late-type supergiant Alpha Orionis. Infrared spectroscopy of red supergiants reveals strong molecular bands, some of which do not originate in the photosphere but in a cooler layer of molecular material above it. Lately, these layers have been spatially resolved by near and mid-IR interferometry. In this paper, we try to reconcile the IR interferometric and ISO-SWS spectroscopic results on Alpha Orionis with a thorough modelling of the photosphere, molecular layer(s) and dust shell. From the ISO and near-IR interferometric observations, we find that Alpha Orionis has only a very low density water layer close above the photosphere. However, mid-IR interferometric observations and a narrow-slit N-band spectrum suggest much larger extra-photospheric opacity close to the photosphere at those wavelengths, even when taking into account the detached dust shell. We argue that this cannot be due to the water layer, and that another source of mid-IR opacity must be present. We show that this opacity source is probably neither molecular nor chromospheric. Rather, we present amorphous alumina (Al2O3) as the best candidate and discuss this hypothesis in the framework of dust-condensation scenarios.

We present near-IR interferometric measurements of the K1.5 giant Arcturus (alpha Bootis), obtained at the IOTA interferometer with the FLUOR instrument, in four narrow filters with central wavelengths ranging from 2.03 to 2.39 micron. These observations were expected to allow us to quantify the wavelength dependence of the diameter of a typical K giant. They are compared to predictions from both plane-parallel and spherical model atmospheres. Unexpectedly, neither can explain the observed visibilities. We show how these data suggest the presence of a companion, in accordance with the Hipparcos data on this star, and discuss this solution with respect to Arcturus' single star status.

We describe state-of-the-art MARCS-code model atmospheres generated for a group of A dwarf, G dwarf, and late-G to mid-K giant standard stars, selected to photometrically calibrate the Spitzer-IRS, and compare the synthetic spectra to observations of HR 6688, HR 6705, and HR 7891. The general calibration processes and uncertainties are briefly described, and the differences between various templated composite spectra of the standards are addressed. In particular, a contrast between up-to-date model atmospheres and previously published composite and synthetic spectra is illustrated for wavelength ranges around 8um (where the SiO Delta(v) = 1 band occurs for the cooler standards) and lambda greater than 20um, where the use of the Engelke function will lead to increasingly large discrepancies due to the neglect of gravity in cool stars. At this point, radiometric requirements are being met, absolute flux calibration uncertainties (1-sigma) are ~20% in the SH and LH, and ~15% in the SL and LL data, and order-to-order flux uncertainties are ~10% or less. Iteration between the MARCS model atmosphere inputs and the data processing will improve the S/N ratios and calibration accuracies.

presented. The 2.38 -- 4.08 micron wavelength-range of band 1 of ISO-SWS (Short-Wavelength Spectrometers on board of the Infrared Space Observatory) in which many different molecules -- with their own dependence on each of the stellar parameters -- are absorbing, enables us to estimate the effective temperature, the gravity, the microturbulence, the metallicity, the CNO-abundances, the 12C/13C-ratio and the angular diameter from the ISO-SWS data. Using the Hipparcos' parallax, the radius, luminosity and gravity-inferred mass are derived. The stellar parameters obtained are in good agreement with other published values, though also some discrepancies with values deduced by other authors are noted. For a few stars (Delta Dra, Xi Dra, Alpha Tuc, H Sco and Alpha Cet) some parameters -- e.g. the CNO-abundances -- are derived for the first time. By examining the correspondence between different ISO-SWS observations of the same object and between the ISO-SWS data and the corresponding synthetic spectrum, it is shown that the relative accuracy of ISO-SWS in band 1 (2.38 -- 4.0 micron) is better than 2% for these high-flux sources. The high level of correspondence between observations and theoretical predictions, together with a confrontation of the estimated Teff(ISO) value with Teff-values derived from colours -- which demonstrates the consistency between V-K, BC(K), Teff and the angular diameter derived from optical or IR data -- proves that both the used MARCS models to derive the stellar quantities and the flux calibration of the ISO-SWS detectors have reached a high level of reliability.

Vega, Sirius, Beta Leo, Alpha Car and Alpha Cen A belong to a sample of twenty stellar sources used for the calibration of the detectors of the Short-Wavelength Spectrometer on board the Infrared Space Observatory (ISO-SWS). While general problems with the calibration and with the theoretical modelling of these stars are reported in Decin et al. (2002), each of these stars is discussed individually in this paper. As demonstrated in Decin et al. (2002), it is not possible to deduce the effective temperature, the gravity and the chemical composition from the ISO-SWS spectra of these stars. But since ISO-SWS is absolutely calibrated, the angular diameter of these stellar sources can be deduced from their ISO-SWS spectra, which consequently yields the stellar radius (R), the gravity-inferred mass (M) and the luminosity (L) for these stars. These deduced parameters are confronted with other published values and the goodness-of-fit between observed ISO-SWS data and the corresponding synthetic spectrum is discussed.

The fine calibration of the ISO-SWS detectors (Infrared Space Observatory - Short Wavelength Spectrometer) has proven to be a delicate problem. We therefore present a detailed spectroscopic study in the 2.38 -- 12 micron wavelength range of a sample of 16 A0 -- M2 stars used for the calibration of ISO-SWS. By investigating the discrepancies between the ISO-SWS data of these sources, the theoretical predictions of their spectra, the high-resolution FTS-KP (Kitt Peak) spectrum of Alpha Boo and the solar FTS-ATMOS (Atmospheric Trace Molecule Spectroscopy) spectrum, both calibration problems and problems in computing the theoretical models and the synthetic spectra are revealed. The underlying reasons for these problems are sought for and the impact on the further calibration of ISO-SWS and on the theoretical modelling is discussed extensively.

We combine spectra of VY CMa obtained with the short- and long-wavelength spectrometers, SWS and LWS, on the Infrared Space Observatory to provide a first detailed continuum spectrum of this highly luminous star. The circumstellar dust cloud through which the star is observed is partially self-absorbing, which makes for complex computational modeling. We review previous work and comment on the range of uncertainties about the physical traits and mineralogical composition of the modeled disk. We show that these uncertainties significantly affect the modeling of the outflow and the estimated mass loss. In particular, we demonstrate that a variety of quite diverse models can produce good fits to the observed spectrum. If the outflow is steady, and the radiative repulsion on the dust cloud dominates the star's gravitational attraction, we show that the total dust mass-loss rate is $\sim 4\times 10^{-6}M_{\odot}$ yr$^{-1}$, assuming that the star is at a distance of 1.5 kpc. Several indications, however, suggest that the outflow from the star may be spasmodic. We discuss this and other problems facing the construction of a physically coherent model of the dust cloud and a realistic mass-loss analysis.

For the astronomical community analyzing ISO-SWS data, a first point to assess when judging and qualifying the observational data concerns the flux calibration accuracy. Since the calibration process is not straightforward and since a wrong calibration may lead to an over- or underestimation of the results, knowledge on the full calibration process and on the still remaining calibration problems is crucial when processing the data. One way to detect calibration problems is by comparing observed data with theoretical predictions of a whole sample of calibration sources. By using an iterative process in which improvements on both the calibration and the theoretical modelling are involved, we will demonstrate that a consistent theoretical data-set of infrared spectra has been constructed. This data-set has been used to derive the OLP10 flux calibration of the ISO-SWS. We will shown that the relative flux calibration accuracy of (high-flux) ISO-SWS observations has reached a 2% level of accuracy in band 1 and that the flux calibration accuracy in band 2 has improved significantly with the introduction of the memory-effect corrrection and the use of our synthetic spectra for the flux calibration derivation, reaching a 6% level of accuracy in band 2.

We present a detailed spectroscopic study of a sample of bright, mostly cool, stars observed with the Short-Wavelength Spectrometer (SWS) on board the Infrared Space Observatory (ISO), which enables the accurate determination of the stellar parameters of the cool giants, but also serves as a critical review of the ISO-SWS calibration.

A detailed spectroscopic study of the ISO-SWS data of the red giant Alpha Tau is presented, which enables not only the accurate determination of the stellar parameters of Alpha Tau, but also serves as a critical review of the ISO-SWS calibration. This study is situated in a broader context of an iterative process in which both accurate observations of stellar templates and cool star atmosphere models are involved to improve the ISO-SWS calibration process as well as the theoretical modelling of stellar atmospheres. Therefore a sample of cool stars, covering the whole A0 -- M8 spectral classification, has been observed in order to disentangle calibration problems and problems in generating the theoretical models and corresponding synthetic spectrum. By using stellar parameters found in the literature large discrepancies were seen between the ISO-SWS data and the generated synthetic spectrum of Alpha Tau. A study of the influence of various stellar parameters on the theoretical models and synthetic spectra, in conjunction with the Kolmogorov-Smirnov test to evaluate objectively the goodness-of-fit, enables us to pin down the stellar parameters with a high accuracy: Teff = 3850 +/- 70 K, log g = 1.50 +/- 0.15, M = 2.3 +/- 0.8 Msun, z = -0.15 +/- 0.20 dex, microturbulence = 1.7 +/- 0.3 km/s, 12C/13C= 10 +/- 1, abundance of C = 8.35 +/- 0.20 dex, abundance of N= 8.35 +/- 0.25 dex, abundance of O = 8.83 +/- 0.15 dex and the angular diameter is 20.77 +/- 0.83 mas. These atmospheric parameters were then compared with the results provided by other authors using other methods and/or spectra.

We present observations of the stretching mode of 13CO2 ice along 13 lines of sight in the Galaxy, using the Short Wavelength Spectrometer on board of the Infrared Space Observatory. Remarkable variations are seen in the absorption band profile in the different lines of sight. The main feature is attributed to 13CO2 mixed with polar molecules such as H2O, and CH3OH. Several high-mass protostars show an additional narrow substructure at 4.382 um, which we attribute to a polar CO2 ice, that experienced heating. This heating effect is sustained by a good correlation of the strength of the substructure with dust and CO gas temperatures along the line of sight, and anti-correlation with ice abundances. Thus, our main conclusion is that interstellar CO2 ices around luminous protostars are subjected to, and altered by, thermal processing and that it may reflect the evolutionary stage of the nearby protostar. In contrast, the ices around low mass protostars and in a quiescent cloud in our sample do not show signs of thermal processing. Furthermore, we determine for the first time the Galactic 12C/13C ratio from the solid state. The 12CO2/13CO2 ratio for the local ISM (69+/-15), as well as the dependence on Galacto-centric radius, are in good agreement with gas phase (C18O, H2CO) studies. For the few individual objects for which gas phase values are available, the 12C/13C ratios derived from CO2 tend to be higher compared to CO studies (albeit with 2.5 sigma significance only). We discuss the implications of this possible difference for the chemical origin of interstellar CO2.

We present infrared spectroscopy and imaging of AFGL~4106. The 2.4-5 micron ISO-SWS spectrum reveals the presence of a cool, luminous star (T_eff ~ 3750 K) in addition to an almost equally luminous F star (T_eff ~ 7250 K). The 5-195 micron SWS and LWS spectra are dominated by strong emission from circumstellar dust. We find that the dust consists of amorphous silicates, with a minor but significant contribution from crystalline silicates. The amorphous silicates consist of Fe-rich olivines. The presence of amorphous pyroxenes cannot be excluded but if present they contain much less Fe than the amorphous olivines. Comparison with laboratory data shows that the pure Mg-end members of the crystalline olivine and pyroxene solid solution series are present. In addition, we find strong evidence for simple oxides (FeO and Al2O3) as well as crystalline H2O ice. Several narrow emission features remain unidentified. Modelling of the dust emission using a dust radiation transfer code shows that large grains (~1 micron) must be present and that the abundance of the crystalline silicates is between 7 and 15% of the total dust mass, depending on the assumed enstatite to forsterite ratio, which is estimated to be between 1 and 3. The amorphous and crystalline dust components in the shell do not have the same temperature, implying that the different dust species are not thermally coupled. We find a dust mass of ~3.9 x 10^-2 M_sol expelled over a period of 4 x 10^3 years for a distance of 3.3 kpc. The F-star in the AFGL~4106 binary is likely a post-red-supergiant in transition to a blue supergiant or WR phase.

We present a detailed chemical analysis on the basis of high-resolution, high signal-to-noise optical spectra of two post-AGB objects IRAS 22223+4327 and IRAS 04296+3429. Both display the unidentified $21 \mu m$ feature in their IR-spectra. The spectroscopic indicators provide accurate atmospheric parameters of $T_{eff}$=6500 K, $log g=1.0$ and $\xi_t = 5.5 km/s$ for IRAS 2223+4327 and $T_{eff}$=7000 K, $log g=1.0$ and $\xi_t = 4.0 km/s$ for IRAS 04296+3429. Both photospheres are found to be metal-deficient with [Fe/H]= -0.4 and -0.7 respectively. C and N are found to be overabundant. The mean abundance of all the measured s-process-elements is [s/Fe]=+1.0 for IRAS 2223+4327 and +1.4 for IRAS 04296+3429. The distribution of the s-process elements can best be described as due to a distribution of neutron exposures with a low mean neutron exposure of $\tau_{0} = 0.2 mbarn^{-1}$. The 21 $\mu$m stars form an interesting sub-group in the total post-AGB sample of stars, not only for their IR characteristics, but also in a broader context of stellar (chemical) evolution theory. They show, in contrast to other post-AGB stars, that the 3rd dredge-up has been efficient during their AGB evolution. The mean neutron exposure is lower than expected for their metallicity.