Search

Despite their large impact on stellar and galactic evolution, the properties of outflows from red supergiants are not well characterised. We used the Onsala 20m telescope to perform a spectral survey at 3mm and 4mm (68 - 116 GHz) of the red supergiant NML Cyg, alongside the yellow hypergiant IRC+10420. Our observations of NML Cyg were combined with complementary archival data to enable a search for signatures of morphological complexity in the circumstellar environment, using emission lines from 15 molecular species. The recovered parameters imply the presence of three distinct, coherent and persistent components, comprised of blue-shifted and red-shifted components, in addition to an underlying outflow centred at the stellar systemic velocity. Furthermore, to reproduce CO emission with three-dimensional radiative transfer models required a spherical outflow with three superposed conical outflows, one towards and one away from the observer, and one in the plane of the sky. These components are higher in density than the spherical outflow by up to an order of magnitude. We hence propose that NML Cyg's circumstellar environment consists of a small number of high-density large-scale coherent outflows embedded in a spherical wind. This would make the mass-loss history similar to that of VY CMa, and distinct from mu Cep, where the outflow contains many randomly distributed smaller clumps. A possible correlation between stellar properties, outflow structures and content is critical in understanding the evolution of massive stars and their environmental impact., 17 pages, 8 figures, 3 appendices, (37 figures and 47 pages in total), accepted for publication in MNRAS

Estimating the distances to asymptotic giant branch (AGB) stars using optical measurements of their parallaxes is not straightforward because of the large uncertainties introduced by their dusty envelopes, their large angular sizes, and their surface brightness variability. We compared the parallaxes from Gaia DR3 with parallaxes measured with maser observations with very long baseline interferometry (VLBI) to determine a statistical correction factor for the DR3 parallaxes using a sub-sample of 33 maser-emitting oxygen-rich nearby AGB stars. We found that the nominal errors on the Gaia DR3 parallaxes are underestimated by a factor of 5.44 for the brightest sources ($G, 15 pages, 12 figures, 3 appendices, accepted for publication in A&A, corrected typos, distance catalogue available on the CDS (and any updates that may occur)

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

This overview paper presentsATOMIUM, a Large Programme in Cycle 6 with the Atacama Large Millimeter/submillimeter Array (ALMA). The goal ofATOMIUMis to understand the dynamics and the gas phase and dust formation chemistry in the winds of evolved asymptotic giant branch (AGB) and red supergiant (RSG) stars. A more general aim is to identify chemical processes applicable to other astrophysical environments. Seventeen oxygen-rich AGB and RSG stars spanning a range in (circum)stellar parameters and evolutionary phases were observed in a homogeneous observing strategy allowing for an unambiguous comparison. Data were obtained between 213.83 and 269.71 GHz at high (∼0″​​.025–0″​​.050), medium (∼0″​​.13–0″​​.24), and low (∼1″) angular resolution. The sensitivity per ∼1.3 km s−1channel was 1.5–5 mJy beam−1, and the line-free channels were used to image the millimetre wave continuum. Our primary molecules for studying the gas dynamics and dust formation are CO, SiO, AlO, AlOH, TiO, TiO2, and HCN; secondary molecules include SO, SO2, SiS, CS, H2O, and NaCl. The scientific motivation, survey design, sample properties, data reduction, and an overview of the data products are described. In addition, we highlight one scientific result – the wind kinematics of theATOMIUMsources. Our analysis suggests that theATOMIUMsources often have a slow wind acceleration, and a fraction of the gas reaches a velocity which can be up to a factor of two times larger than previously reported terminal velocities assuming isotropic expansion. Moreover, the wind kinematic profiles establish that the radial velocity described by the momentum equation for a spherical wind structure cannot capture the complexity of the velocity field. In fifteen sources, some molecular transitions other than12COv = 0J = 2 − 1 reach a higher outflow velocity, with a spatial emission zone that is often greater than 30 stellar radii, but much less than the extent of CO. We propose that a binary interaction with a (sub)stellar companion may (partly) explain the non-monotonic behaviour of the projected velocity field. TheATOMIUMdata hence provide a crucial benchmark for the wind dynamics of evolved stars in single and binary star models.

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

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

Aims.In this paper, we aim to constrain the dust mass and grain sizes in the interaction regions between the stellar winds and the interstellar medium (ISM) around asymptotic giant branch (AGB) stars. By describing the dust in these regions, we aim to shed light on the role of evolved low-mass stars in the origin of dust in galaxies.Methods.We use images in the far-infrared (FIR) at 70 and 160 µm to derive dust temperatures and dust masses in the wind-ISM interaction regions around a sample of carbon-rich and oxygen-rich AGB stars. The dust temperature and mass are determined in two ways: first, directly from the data using the ratio of the measured fluxes and assuming opacities for dust with a constant grain size of 0.1 µm, and then using three-dimensional dust-radiative transfer models spatially constrained by the observations. Each of the radiative transfer models contains one constant grain size, which is varied between 0.01 and 5.0 µm.Results.We find that the observed dust mass in the wind-ISM interaction regions is consistent with mass accumulated from the stellar winds. For the carbon-rich sources, adding the spatial constraints in the radiative transfer models results in preferentially larger grain sizes (≈2 µm). For the oxygen-rich sources, the spatial constraints result in overly high temperatures in the models, making it impossible to fit the observed FIR ratio irrespective of the grain size used, indicating a more complex interplay of grain properties and the stellar radiation field.Conclusions.Our results have implications for how likely it is for the grains to survive the transition into the ISM, and the properties of dust particles that later act as seeds for grain growth in the ISM. However, our results for the oxygen-rich sources show that the derivation of dust properties is not straight forward, requiring more complex modelling.

The Origins Space Telescope (Origins) is one of four science and technology definition studies selected by the National Aeronautics and Space Administration (NASA) in preparation of the 2020 Astronomy and Astrophysics Decadal survey in the US. Origins will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. It is designed to answer three major science questions: How do galaxies form stars, make metals, and grow their central supermassive black holes from reionization? How do the conditions for habitability develop during the process of planet formation? Do planets orbiting M-dwarf stars support life? Origins operates at mid- to far-infrared wavelengths from ~ 2.8 μm to 588 μm, and is more than 1000 times more sensitive than prior far-IR missions due to its cold (~ 4.5 K) aperture and state-of-the-art instruments.

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

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$., Comment: 16 pages, 9 figures, 3 tables

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

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

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

The Origins Space Telescope will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. How did galaxies evolve from the earliest galactic systems to those found in the universe today? How do habitable planets form? How common are life-bearing worlds? To answer these alluring questions, Origins will operate at mid- and far-infrared wavelengths and offer powerful spectroscopic instruments and sensitivity three orders of magnitude better than that of Herschel, the largest telescope flown in space to date. After a 3 ½ year study, the Origins Science and Technology Definition Team will recommend to the Decadal Survey a concept for Origins with a 5.9-m diameter telescope cryocooled to 4.5 K and equipped with three scientific instruments. A mid-infrared instrument (MISC-T) will measure the spectra of transiting exoplanets in the 2.8 – 20 μm wavelength range and offer unprecedented sensitivity, enabling definitive biosignature detections. The Far-IR Imager Polarimeter (FIP) will be able to survey thousands of square degrees with broadband imaging at 50 and 250 μm. The Origins Survey Spectrometer (OSS) will cover wavelengths from 25 – 588 μm, make wide-area and deep spectroscopic surveys with spectral resolving power R ~ 300, and pointed observations at R ~ 40,000 and 300,000 with selectable instrument modes. Origins was designed to minimize complexity. The telescope has a Spitzer-like architecture and requires very few deployments after launch. The cryo-thermal system design leverages JWST technology and experience. A combination of current-state-of-the-art cryocoolers and next-generation detector technology will enable Origins’ natural backgroundlimited sensitivity.

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

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

Eleven nearby (~ 0.55 mag; AK Hya, V744 Cen, RU Crt, alpha Her). Stars below this limit were not detected (BQ Gem, eps Oct, NU Pav, II Hya, CL Hyi, ET Vir, SX Pav). Ks-[22] colour is found to trace mass-loss rate to well within an order of magnitude. This confirms existing results, indicating a factor of 100 increase in AGB-star mass-loss rates at a pulsation period of ~60 days, similar to the known "superwind" trigger at ~300 days. Between ~60 and ~300 days, an approximately constant mass-loss rate and wind velocity of ~3.7 x 10^-7 solar masses per year and ~8 km/s is found. While this has not been corrected for observational biases, this rapid increase in mass-loss rate suggests a need to recalibrate the treatment of AGB mass loss in stellar evolution models. The comparative lack of correlation between mass-loss rate and luminosity (for L, Comment: 18 pages, 7 figures, accepted for publication in MNRAS

The Origins Space Telescope (OST) will trace the history of our origins from the time dust and heavy elements permanently altered the cosmic landscape to present-day life. How did the universe evolve in response to its changing ingredients? How common are life-bearing planets? To accomplish its scientific objectives, OST will operate at mid- and far-infrared wavelengths and offer superlative sensitivity and new spectroscopic capabilities. The OST study team will present a scientifically compelling, executable mission concept to the 2020 Decadal Survey in Astrophysics. To understand the concept solution space, our team studied two alternative mission concepts. We report on the study approach and describe both of these concepts, give the rationale for major design decisions, and briefly describe the mission-enabling technology.

© ESO 2018. Context. With the outstanding spatial resolution and sensitivity of the Atacama Large Millimeter/sub-millimeter Array (ALMA), molecular gas other than the abundant CO can be observed and resolved in circumstellar envelopes (CSEs) around evolved stars, such as the binary S-type asymptotic giant branch (AGB) star W Aquilae. Aims. We aim to constrain the chemical composition of the CSE and determine the radial abundance distribution, the photospheric peak abundance, and isotopic ratios of a selection of chemically important molecular species in the innermost CSE of W Aql. The derived parameters are put into the context of the chemical evolution of AGB stars and are compared with theoretical models. Methods. We employ one-dimensional radiative transfer modeling - with the accelerated lambda iteration (ALI) radiative transfer code-of the radial abundance distribution of a total of five molecular species (CS, SiS, 30SiS, 29SiO and H13CN) and determine the best fitting model parameters based on high-resolution ALMA observations as well as archival single-dish observations. The additional advantage of the spatially resolved ALMA observations is that we can directly constrain the radial profile of the observed line transitions from the observations. Results. We derive abundances and e-folding radii for CS, SiS, 30SiS, 29SiO and H13CN and compare them to previous studies, which are based only on unresolved single-dish spectra. Our results are in line with previous results and are more accurate due to resolution of the emission regions. adsnote: Provided by the SAO/NASA Astrophysics Data System adsurl: http://cdsads.u-strasbg.fr/abs/2018A%26A...617A..23B archiveprefix: arXiv date-added: 2018-09-27 16:25:26 +0200 date-modified: 2018-09-27 16:25:35 +0200 eid: A23 keywords: stars: abundances, stars: AGB and post-AGB, circumstellar matter, stars: mass-loss, stars: winds, outflows, submillimeter: stars primaryclass: astro-ph.SR bdsk-url-1: https://doi.org/10.1051/0004-6361/201832724 ispartof: ASTRONOMY & ASTROPHYSICS vol:617 status: published

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

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

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

Context.Sulphur has long been known to form different molecules depending on the chemical composition of its environment. More recently, the sulphur-bearing molecules SO and H2S have been shown to behave differently in oxygen-rich asymptotic giant branch (AGB) circumstellar envelopes of different densities.Aims.By surveying a diverse sample of AGB stars for CS and SiS emission, we aim to determine in which environments these sulphur-bearing molecules most readily occur. We include sources with a range of mass-loss rates and carbon-rich, oxygen-rich, and mixed S-type chemistries. Where these molecules are detected, we aim to determine their CS and SiS abundances.Methods.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.Results.We detect CS towards all the surveyed carbon stars, some S-type stars, and the highest mass-loss rate oxygen-rich stars, (Ṁ≥ 5 × 10−6M⊙ yr−1). SiS is detected towards the highest mass-loss rate sources of all chemical types (Ṁ≥ 8 × 10−7M⊙ yr−1). We find CS peak fractional abundances ranging from ~4 × 10−7to ~2 × 10−5for the carbon stars, from ~3 × 10−8to ~1 × 10−7for the oxygen-rich stars, and from ~1 × 10−7to ~8 × 10−6for the S-type stars. We find SiS peak fractional abundances ranging from ~9 × 10−6to ~2 × 10−5for the carbon stars, from ~5 × 10−7to ~2 × 10−6for the oxygen-rich stars, and from ~2 × 10−7to ~2 × 10−6for the S-type stars.Conclusions.Overall, we find that wind density plays an important role in determining the chemical composition of AGB circumstellar envelopes. It is seen that for oxygen-rich AGB stars both CS and SiS are detected only in the highest density circumstellar envelopes and their abundances are generally lower than for carbon-rich AGB stars by around an order of magnitude. For carbon-rich and S-type stars SiS was also only detected in the highest density circumstellar envelopes, while CS was detected consistently in all surveyed carbon stars and sporadically among the S-type stars.

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

Context. The S-type asymptotic giant branch (AGB) star ?1 Gru has a known companion at a separation of 2? 7 (? 400 AU). Previous observations of the circumstellar envelope (CSE) show strong deviations from spherical symmetry. The envelope structure, including an equatorial torus and a fast bipolar outflow, is rarely seen in the AGB phase and is particularly unexpected in such a wide binary system. Therefore a second, closer companion has been suggested, but the evidence is not conclusive. Aims. The aim is to make a 3D model of the CSE and to constrain the density and temperature distribution using new spatially resolved observations of the CO rotational lines. Methods. We have observed the J = 3-2 line emission from 12CO and 13CO using the compact arrays of the Atacama Large Millimeter/submillimeter Array (ALMA). The new ALMA data, together with previously published 12CO J = 2-1 data from the Submillimeter Array (SMA), and the 12CO J = 5-4 and J = 9-8 lines observed with Herschel/Heterodyne Instrument for the Far-Infrared (HIFI), is modeled with the 3D non-LTE radiative transfer code SHAPEMOL. Results. The data analysis clearly confirms the torus-bipolar structure. The 3D model of the CSE that satisfactorily reproduces the data consists of three kinematic components: a radially expanding torus with velocity slowly increasing from 8 to 13 km s-1 along the equator plane; a radially expanding component at the center with a constant velocity of 14 km s-1; and a fast, bipolar outflow with velocity proportionally increasing from 14 km s-1 at the base up to 100 km s-1 at the tip, following a linear radial dependence. The results are used to estimate an average mass-loss rate during the creation of the torus of 7.7 × 10-7 M? yr-1. The total mass and linear momentum of the fast outflow are estimated at 7.3 × 10-4 M? and 9.6 × 1037 g cm s-1, respectively. The momentum of the outflow is in excess (by a factor of about 20) of what could be generated by radiation pressure alone, in agreement with recent findings for more evolved sources. The best-fit model also suggests a 12CO/13CO abundance ratio of 50. Possible shaping scenarios for the gas envelope are discussed.

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., Comment: 9 pages, 9 figures, Accepted for publication in Astronomy & Astrophysics

Context. Born-Again stars provide a unique possibility to study the evolution of the circumstellar envelope of evolved stars in human timescales. Up until now, most of the observations of the circumstellar material in these stars have been limited to studying the relatively hot gas and dust. In other evolved stars, the emission from rotational transitions of molecules, such as CO, is commonly used to study the cool component of their circumstellar envelopes. Thus, the detection and study of molecular gas in born-Again stars is of great importance when attempting to understand their composition and chemical evolution. In addition, the molecular emission is an invaluable tool for exploring the physical conditions, kinematics, and formation of asymmetric structures in the circumstellar envelopes of these evolved stars. However, up until now, all attempts to detect molecular emission from the cool material around born-Again stars have failed. Aims. We searched for emission from rotational transitions of molecules in the hydrogen-deficient circumstellar envelopes of bornagain stars to explore the chemical composition, kinematics, and physical parameters of the relatively cool gas. Methods. We carried out observations using the APEX and IRAM 30 m telescopes to search for molecular emission toward four well-studied born-Again stars, V4334 Sgr, V605 Aql, A30, and A78, that are thought to represent an evolutionary sequence. Results. For the first time, we detected emission from HCN and HCN molecules toward V4334 Sgr, and CO emission in V605 Aql. No molecular emission was detected above the noise level toward A30 and A78. The detected lines exhibit broad linewidths 150 km s, which indicates that the emission comes from gas ejected during the born-Again event, rather than from the old planetary nebula. A first estimate of the HCN/HCN abundance ratio in the circumstellar environment of V4334 Sgr is ≈ 3, which is similar to the value of the 12C/13C ratio measured from other observations. We derived a rotational temperature of Trot = 13 ± 1 K, and a total column density of NHCN = 1:6 ± 0:1 × 10 cm for V4334 Sgr. This result sets a lower limit on the amount of hydrogen that was ejected into the wind during the born-Again event of this source. For V605 Aql, we obtained a lower limit for the integrated line intensities I12CO/I13CO > 4.©ESO 2017., D.T. and W.V. acknowledge support from ERC consolidator grant 614264. A.A.Z. acknowledges funding from the UK STFC research council. A. S.-M. thanks the Deutsche Forschungsgemeinschaft (DFG) for funding support via the collaborative research grant SFB 956, project A6. S.P.T.M. thanks the Spanish MINECO for funding support from grants AYA2012-32032, CSD2009-00038, FIS2012-32096, and ERC under ERC-2013-SyG, G. A. 610256 NANOCOSMOS

[Context] Recent observations at subarcsecond resolution, now possible also at submillimeter wavelengths, have shown intricate circumstellar structures around asymptotic giant branch (AGB) stars, mostly attributed to binary interaction. The results presented here are part of a larger project aimed at investigating the effects of a binary companion on the morphology of circumstellar envelopes (CSEs) of AGB stars., [Aims] AGB stars are characterized by intense stellar winds that build CSEs around the stars. Here, the CO(J = 3 → 2) emission from the CSE of the binary S-Type AGB star W Aql has been observed at subarcsecond resolution using ALMA. The aim of this paper is to investigate the wind properties of the AGB star and to analyse how the known companion has shaped the CSE., [Methods] The average mass-loss rate during the creation of the detected CSE is estimated through modelling, using the ALMA brightness distribution and previously published single-dish measurements as observational constraints. The ALMA observations are presented and compared to the results from a 3D smoothed particle hydrodynamics (SPH) binary interaction model with the same properties as the W Aql system and with two different orbital eccentricities. Three-dimensional radiative transfer modelling is performed and the response of the interferometer is modelled and discussed., [Results] The estimated average mass-loss rate of W Aql is = 3.0 × 10M yr and agrees with previous results based on single-dish CO line emission observations. The size of the emitting region is consistent with photodissociation models. The inner 10′′ of the CSE is asymmetric with arc-like structures at separations of 2-3′′ scattered across the denser sections. Further out, weaker spiral structures at greater separations are found, but this is at the limit of the sensitivity and field of view of the ALMA observations., [Conclusions] The CO(J = 3 → 2) emission is dominated by a smooth component overlayed with two weak arc patterns with different separations. The larger pattern is predicted by the binary interaction model with separations of ~10′′ and therefore likely due to the known companion. It is consistent with a binary orbit with low eccentricity. The smaller separation pattern is asymmetric and coincides with the dust distribution, but the separation timescale (200 yr) is not consistent with any known process of the system. The separation of the known companions of the system is large enough to not have a very strong effect on the circumstellar morphology. The density contrast across the envelope of a binary with an even larger separation will not be easily detectable, even with ALMA, unless the orbit is strongly asymmetric or the AGB star has a much larger mass-loss rate., S.M. is grateful to the South African National Research Foundation (NRF) for a research grant. W.V. acknowledges support from ERC consolidator grant 614264. T.D. acknowledges support from the ERC consolidator grant 646758 AEROSOL and the FWO Research Project grant G024112N. G.Q.L. acknowledges support from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007−2013)/ERC Grant Agreement No. 610256 (NANOCOSMOS).

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

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

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

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

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

(abreviated) We aim to study the inner-wind structure (R, 17 pages, Astronomy & Astrophysics

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

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

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

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

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., Accepted for publication in Astronomy & Astrophysics, 25 pages, 20 figures

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

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

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., Comment: Letter 4 pages, 5 figures plus appendix with 3 figures. Accepted by Astronomy and Astrophysics Letters

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.

© 2016 ESO. Aims. We aim to probe the close and distant circumstellar environments of the stellar outburst object V838 Mon. Methods. Herschel far-infrared 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 dust immediately around V838 Mon, and in the surrounding infrared-bright region. These maps were fitted in 1d and 2d to measure the temperature, mass, and β of the two dust sources. PACS and SPIRE spectra were used to detect emission lines from the extended atmosphere of the star, which were then modelled to study the physical conditions in the emitting material. HIFI spectra were taken to measure the kinematics of the extended atmosphere but unfortunately yielded no detections. Results. Fitting of the far-infrared imaging of V838 Mon reveals 0.5-0.6 M⊙ of ≈ 19 K dust in the environs (≈ 2.7 pc) surrounding V838 Mon. The surface-integrated infrared flux (signifying the thermal light echo), and derived dust properties do not vary significantly between the different epochs. We measured the photometry of the point source. As the peak of the SED (Spectral Energy Distribution) lies outside the Herschel spectral range, it is only by incorporating data from other observatories and previous epochs that we can usefully fit the SED; with this we explicitly assume no evolution of the point source between the epochs. We find that warm dust with a temperature ~ 300 K distributed over a radius of 150-200 AU. 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 ≠18 K) that could be associated with the ISM or possibly with a cold layer in the outermost part of the shell, and a warm (Tkin ≈ 400 K) zone that is associated with the extended environment of V838 Mon within a region of radius of ≈ 210 AU. 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. ispartof: Astronomy & Astrophysics vol:596 pages:23-23 status: published

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 (, Accepted for publication in Astronomy & Astrophysics, 10 pages, 8 figures

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., Comment: 25 pags. 18 figs

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