Search

The spatial distribution of molecules in AGB circumstellar envelopes is regulated by different processes. In the outer layers all molecules are destroyed due to the interaction with interstellar ultraviolet photons. Here we aim to characterize in a coherent and uniform way the radial extent of three molecules (SiO, CS, and SiS) in envelopes around AGB stars of O- and C-rich character, and to study their dependence with mass loss rate. To that purpose, we used the Yebes 40m and IRAM 30m telescopes to observe 7 M-type and 7 C-type AGB envelopes covering a wide range of mass loss rates (1e-7 - 1e-5 Msun/yr) in lines of SiO, CS, and SiS spanning a range of upper level energies of 2-130 K. We carried out excitation and radiative transfer calculations over a wide parameter space to characterize the molecular abundance and radial extent. A chi2 analysis indicates that the abundance is well constrained while the radial extent is more difficult to constrain. The radial extent increases with increasing envelope density, in agreement with previous observational findings. At high envelope densities, Mdot/vexp > 1e-6 (Msun/yr)/(km/s), the radial extent of SiO, CS, and SiS are similar, while at low envelope densities, Mdot/vexp < 1e-7 (Msun/yr)/(km/s), the radial extent differ among the three molecules, in agreement with theoretical expectations based on destruction due to photodissociation. At low envelope densities we find a sequence of increasing radial extent, SiS -> CS -> SiO. We also find a tentative dependence of the radial extent with the chemical type (O- or C-rich) of the star for SiO and CS. Interferometric observations and further investigation of the photodissociation of SiO, CS, and SiS should allow to clarify the situation on the relative photodissociation radius of SiO, CS, and SiS in AGB envelopes and the dependence with envelope density and C/O ratio., Comment: Accepted for publication in A&A. Abstract abridged

Polarisation observations of masers in the circumstellar envelopes (CSEs) around Asymptotic Giant Branch (AGB) stars have revealed strong magnetic fields. However, masers probe only specific lines-of-sight through the CSE. Non-masing molecular line polarisation observation can more directly reveal the large scale magnetic field morphology and hence probe the effect of the magnetic field on AGB mass-loss and the shaping of the AGB wind. Observations and models of CSE molecular line polarisation can now be used to describe the magnetic field morphology and estimate its strength throughout the entire CSE. We use observations taken with ALMA of molecular line polarisation in the envelope of two AGB stars (CW~Leo and R~Leo). We model the observations using the multi-dimensional polarised radiative transfer tool PORTAL. We find linearly polarised emission, with maximum fractional polarisation on the order of a few percent, in several molecular lines towards both stars. We can explain the observed differences in polarisation structure between the different molecular lines by alignment of the molecules by a combination of the Goldreich-Kylafis effect and radiative alignment effects. We specifically show that the polarisation of CO traces the morphology of the magnetic field. Competition between the alignment mechanisms allows us to describe the behaviour of the magnetic field strength with radius throughout the circumstellar envelope of CW~Leo. The magnetic field strength derived using this method is inconsistent with the magnetic field strength derived using a structure function analysis of the CO polarisation and the strength previously derived using CN Zeeman observations. In contrast with CW~Leo, the magnetic field in the outer envelope of R~Leo appears to be advected outwards by the stellar wind. {abridged abstract}, Comment: 23 pages, 21 figures, accepted for publication in A&A

Red supergiant stars (RSGs, Minit = 10-40Msun) are known to eject large amounts of material, as much as half of their initial mass during this evolutionary phase. However, the processes powering the mass ejection in low- and intermediate-mass stars do not work for RSGs and the mechanism that drives the ejection remains unknown. Different mechanisms have been proposed as responsible for this mass ejection but so far little is known about the actual processes taking place in these objects. Here we present high angular resolution interferometric ALMA maps of VY CMa continuum and molecular emission, which resolve the structure of the ejecta with unprecedented detail. The study of the molecular emission from the ejecta around evolved stars has been shown to be an essential tool in determining the characteristics of the mass loss ejections. Our aim is thus to use the information provided by these observations to understand the ejections undergone by VY CMa and to determine their possible origins. We inspected the kinematics of molecular emission observed. We obtained position-velocity diagrams and reconstructed the 3D structure of the gas traced by the different species. It allowed us to study the morphology and kinematics of the gas traced by the different species surrounding VY CMa. Two types of ejecta are clearly observed: extended, irregular, and vast ejecta surrounding the star that are carved by localized fast outflows. The structure of the outflows is found to be particularly flat. We present a 3D reconstruction of these outflows and proof of the carving. This indicates that two different mass loss processes take place in this massive star. We tentatively propose the physical cause for the formation of both types of structures. These results provide essential information on the mass loss processes of RSGs and thus of their further evolution., Comment: 13 pages, 10 figures

We present ALMA continuum and molecular line emission maps at $\sim$1 mm of OH 231.8, a well studied bipolar nebula around an AGB star. The excellent angular resolution of our maps ($\sim$20 mas) allows us to scrutinise the central nebular regions of OH 231.8, which hold the clues to unravel how this iconic object assembled its complex nebular architecture. We report, for the first time in this object and others of its kind, the discovery of a rotating circumbinary disk selectively traced by NaCl, KCl, and H$_2$O emission lines. This represents the first detection of KCl in an oxygen-rich AGB circumstellar envelope. The rotating disk, of radius $\sim$30 au, lies at the base of a young bipolar wind traced by SiO and SiS emission, which also presents signs of rotation at its base. The NaCl equatorial structure is characterised by a mean rotation velocity of $\sim$4 km s$^{-1}$ and extremely low expansion speeds, $\sim$3 km s$^{-1}$. The outflow has a predominantly expansive kinematics characterised by a constant radial velocity gradient of $\sim$65 km s$^{-1}$ arcsec$^{-1}$ at its base. Beyond $r$$\sim$350 au, the gas in the outflow continues radially flowing at a constant terminal speed of $\sim$16 km s$^{-1}$. Our continuum maps reveal a spatially resolved dust disk-like structure perpendicular to the outflow, with the NaCl, KCl and H$_2$O emission arising from the disk's surface layers. Within the disk, we also identify an unresolved point continuum source, which likely represents the central Mira-type star QX Pup enshrouded by a $\sim$3 $R_{\star}$ component of hot ($\sim$1400 K) freshly formed dust. The point source is slightly off-centered from the disk centroid, enabling us for the first time to place constraints to the orbital separation and period of the central binary system, $a$$\sim$20 au and $P_{\rm orb}$$\sim$55 yr, respectively. (abridged)., Comment: 30 pages, 18 figures. Accepted for publication in Astronomy & Astrophysics (on 19/06/2022)

The nucleosynthesis production of fluorine (F) is still a matter of debate. Asymptotic giant branch (AGB) stars are one of the main candidates for F production. However, their contribution to the total F budget is not fully known due to the lack of observations. In this paper, we report the detection of AlF line emission, one of the two main carriers of F in the gas-phase in the outflow of evolved stars, towards five nearby oxygen-rich AGB stars, $o$ Ceti, R Leo, IK Tau, R Dor, and W Hya. From spatially resolved observations, we estimated the AlF emitting region with a radius $\sim11R_{\star}$ for $o$ Ceti and $\sim9R_{\star}$ for R Leo. From population diagram analysis, we report the AlF column densities of $\sim 5.8\times10^{15}$ cm$^{-2}$ and $\sim 3\times10^{15}$ cm$^{-2}$ for $o$ Ceti and R Leo, respectively, within these regions. For $o$ Ceti, we used the C$^{18}$O ($v=0$, $J=3-2$) observations to estimate the H$_2$ column density of the emitting region. We found a fractional abundance of $f_{\rm AlF/H_2}\sim(2.5\pm1.7)\times10^{-8}$. This gives a lower limit on the F budget in $o$ Ceti and is compatible with the solar F budget $f_{\rm F/H_2}=(5\pm2)\times10^{-8}$. For R Leo, a fractional abundance $f_{\rm AlF/H_2}=(1.2\pm0.5)\times10^{-8}$ is estimated. For other sources, we cannot precisely determine the emitting region based on the available data. Assuming an emitting region with a radius of $\sim 11R_{\star}$ and the rotational temperatures derived for $o$ Ceti and R Leo, we crudely approximated the AlF column density to be $\sim(1.2-1.5)\times10^{15}$ cm$^{-2}$ in W Hya, $\sim(2.5-3.0)\times10^{14}$ cm$^{-2}$ in R Dor, and $\sim(0.6-1.0)\times10^{16}$ cm$^{-2}$ in IK Tau. These result in fractional abundances within a range of $f_{\rm AlF/H_2}\sim(0.1-4)\times10^{-8}$ in W Hya, R Dor, and IK Tau., Comment: 11 pages, 9 figures

Aims: We aim to determine the abundances of SiO, CS, SiS, SO, and SO$_2$ in a large sample of oxygen-rich AGB envelopes covering a wide range of mass loss rates to investigate the potential role that these molecules could play in the formation of dust in these environments. Methods:We surveyed a sample of 30 oxygen-rich AGB stars in the $\lambda$ 2 mm band using the IRAM 30m telescope. We performed excitation and radiative transfer calculations based on the LVG method to model the observed lines of the molecules and to derive their fractional abundances in the observed envelopes. Results:We detected SiO in all 30 targeted envelopes, as well as CS, SiS, SO, and SO$_2$ in 18, 13, 26, and 19 sources, respectively. Remarkably, SiS is not detected in any envelope with a mass loss rate below $10^{-6}$ M$_{\odot}$ yr$^{-1}$, whereas it is detected in all envelopes with mass loss rates above that threshold. From a comparison with a previous, similar study on C-rich sources, it becomes evident that the fractional abundances of CS and SiS show a marked differentiation between C-rich and O-rich sources, being two orders of magnitude and one order of magnitude more abundant in C-rich sources, respectively, while the fractional abundance of SiO turns out to be insensitive to the C/O ratio. The abundance of SiO in O-rich envelopes behaves similarly to C-rich sources, that is, the denser the envelope the lower its abundance. A similar trend, albeit less clear than for SiO, is observed for SO in O-rich sources. Conclusions: The marked dependence of CS and SiS abundances on the C/O ratio indicates that these two molecules form more efficiently in C- than O-rich envelopes. The decline in the abundance of SiO with increasing envelope density and the tentative one for SO indicate that SiO and possibly SO act as gas-phase precursors of dust in circumstellar envelopes around O-rich AGB stars., Comment: 22 pages, 3 figures, Accepted for publications in Astronomy & Astrophysics

The study of the gas in the envelopes surrounding asymptotic giant branch (AGB) stars through observations in the millimetre wavelength range provides information about the history and nature of these molecular factories. Here we present ALMA observations at subarsecond resolution, complemented with IRAM-30m data, of several lines of SiO, SiS, and CS towards the best-studied AGB circumstellar envelope, IRC+10216. We aim to characterise their spatial distribution and determine their fractional abundances mainly through radiative transfer and chemical modelling. The three species display extended emission with several enhanced emission shells. CS displays the most extended distribution reaching distances up to approximately 20''. SiS and SiO emission have similar sizes of approximately 11'', but SiS emission is slightly more compact. We have estimated fractional abundances relative to H$_2$, which on average are equal to f(SiO)$\sim$10$^{-7}$, f(SiS)$\sim$10$^{-6}$, and f(CS)$\sim$10$^{-6}$ up to the photo-dissociation region. The observations and analysis presented here show evidence that the circumstellar material displays clear deviations from an homogeneous spherical wind, with clumps and low density shells that may allow UV photons from the interstellar medium (ISM) to penetrate deep into the envelope, shifting the photo-dissociation radius inwards. Our chemical model predicts photo-dissociation radii compatible with those derived from the observations, although it is unable to predict abundance variations from the starting radius of the calculations ($\sim$10$R_{*}$), which may reflect the simplicity of the model. We conclude that the spatial distribution of the gas proves the episodic and variable nature of the mass loss mechanism of IRC+10216, on timescales of hundreds of years., Comment: 21 pages, 24 figures. Accepted for publication in A&A

We report on the detection of calcium isocyanide, CaNC, in the carbon-rich evolved star IRC+10216. We derived a column density for this species of (2.0$\pm$0.5)$\times$10$^{11}$ cm$^{-2}$. Based on the observed line profiles and the modelling of its emission through the envelope, the molecule has to be produced in the intermediate and outer layers of the circumstellar envelope where other metal-isocyanides have previously been found in this source. The abundance ratio of CaNC relative to MgNC and FeCN is $\simeq$1/60 and $\simeq$1, respectively. We searched for the species CaF, CaCl, CaC, CaCCH, and CaCH$_3$ for which accurate frequency predictions are available. Only upper limits have been obtained for these molecules., Comment: 5 pages, 2 figures, accepted for publication in A&A

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., Comment: Submitted to Astro2020 decadal survey; 8 pages, 2 figures

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., Comment: Published in Astronomy and Astrophysics Letters

We present new interferometer molecular observations of R Leo taken at 1.2 mm with the Atacama Large Millimeter Array with an angular resolution up to ~0.026 arcsec. These observations permit us to resolve the innermost envelope of this star revealing the existence of a complex structure that involves extended continuum emission and molecular emission showing a non-radial gas velocity distribution. This molecular emission displays prominent red-shifted absorptions located right in front to the star typical of material infall and lateral gas motions compatible with the presence of a torus-like structure., Comment: 5 pages, 3 figures, 1 table; accepted for publication in A&A

We present continuum and molecular line emission ALMA observations of OH 231.8+4.2, a well studied bipolar nebula around an asymptotic giant branch (AGB) star. The high angular resolution (~0.2-0.3 arcsec) and sensitivity of our ALMA maps provide the most detailed and accurate description of the overall nebular structure and kinematics of this object to date. We have identified a number of outflow components previously unknown. Species studied in this work include 12CO, 13CO, CS, SO, SO2, OCS, SiO, SiS, H3O+, Na37Cl, and CH3OH. The molecules Na37Cl and CH3OH are first detections in OH 231.8+4.2, with CH3OH being also a first detection in an AGB star. Our ALMA maps bring to light the totally unexpected position of the mass-losing AGB star (QX Pup) relative to the large-scale outflow. QX Pup is enshrouded within a compact (<60 AU) parcel of dust and gas (clump S) in expansion (V~5-7 km/s) that is displaced by 0.6arcsec to the south of the dense equatorial region (or waist) where the bipolar lobes join. Our SiO maps disclose a compact bipolar outflow that emerges from QX Pup's vicinity. This outflow is oriented similarly to the large-scale nebula but the expansion velocities are about ten times lower (~35 km/s). We deduce short kinematical ages for the SiO outflow, ranging from ~50-80 yr, in regions within ~150 AU, to ~400-500 yr at the lobe tips (~3500 AU). Adjacent to the SiO outflow, we identify a small-scale hourglass-shaped structure (mini-hourglass) that is probably made of compressed ambient material formed as the SiO outflow penetrates the dense, central regions of the nebula. The lobes and the equatorial waist of the mini-hourglass are both radially expanding with a constant velocity gradient. The mini-waist is characterized by extremely low velocities, down to ~1 km/s at ~150 AU, which tentatively suggest the presence of a stable structure. (abridged)

We present new high angular resolution interferometer observations of the v=0 J=14-13 and 15-14 SiS lines towards IRC+10216, carried out with CARMA and ALMA. The maps, with angular resolutions of ~0.25"and 0.55", reveal (1) an extended, roughly uniform, and weak emission with a size of ~0.5", (2) a component elongated approximately along the East-West direction peaking at ~0.13" and 0.17" at both sides of the central star, and (3) two blue- and red-shifted compact components peaking around 0.07" to the NW of the star. We have modeled the emission with a 3D radiation transfer code finding that the observations cannot be explained only by thermal emission. Several maser clumps and one arc-shaped maser feature arranged from 5 to 20R* from the central star, in addition to a thin shell-like maser structure at ~13R* are required to explain the observations. This maser emitting set of structures accounts for 75% of the total emission while the other 25% is produced by thermally excited molecules. About 60% of the maser emission comes from the extended emission and the rest from the set of clumps and the arc. The analysis of a time monitoring of these and other SiS and 29SiS lines carried out with the IRAM 30m telescope from 2015 to present suggests that the intensity of some spectral components of the maser emission strongly depends on the stellar pulsation while other components show a mild variability. This monitoring evidences a significant phase lag of ~0.2 between the maser and NIR light-curves., Comment: 19 pages, 11 figures, 3 tables; accepted for publication in The ApJ

Silicon carbide dust is ubiquitous in circumstellar envelopes around C-rich AGB stars. However, the main gas-phase precursors leading to the formation of SiC dust have not yet been identified. The most obvious candidates among the molecules containing an Si--C bond detected in C-rich AGB stars are SiC2, SiC, and Si2C. We aim to study how widespread and abundant SiC2, SiC, and Si2C are in envelopes around C-rich AGB stars and whether or not these species play an active role as gas-phase precursors of silicon carbide dust in the ejecta of carbon stars. We carried out sensitive observations with the IRAM 30m telescope of a sample of 25 C-rich AGB stars to search for emission lines of SiC2, SiC, and Si2C in the 2 mm band. We performed non-LTE excitation and radiative transfer calculations based on the LVG method to model the observed lines of SiC2 and to derive SiC2 fractional abundances in the observed envelopes. We detect SiC2 in most of the sources, SiC in about half of them, and do not detect Si2C in any source, at the exception of IRC +10216. Most of these detections are reported for the first time in this work. We find a positive correlation between the SiC and SiC2 line emission, which suggests that both species are chemically linked, the SiC radical probably being the photodissociation product of SiC2 in the external layer of the envelope. We find a clear trend in which the denser the envelope, the less abundant SiC2 is. The observed trend is interpreted as an evidence of efficient incorporation of SiC2 onto dust grains, a process which is favored at high densities owing to the higher rate at which collisions between particles take place. The observed behavior of a decline in the SiC2 abundance with increasing density strongly suggests that SiC2 is an important gas-phase precursor of SiC dust in envelopes around carbon stars., Comment: Published in A&A. 16 pages and 10 figures

ALMA is providing us essential information on where certain molecules form. Observing where these molecules emission arises from, the physical conditions of the gas, and how this relates with the presence of other species allows us to understand the formation of many species, and to significantly improve our knowledge of the chemistry that occurs in the space. We studied the molecular distribution of NaCN around IRC +10216, a molecule detected previously, but whose origin is not clear. High angular resolution maps allow us to model the abundance distribution of this molecule and check suggested formation paths. We modeled the emission of NaCN assuming local thermal equilibrium (LTE) conditions. These profiles were fitted to azimuthal averaged intensity profiles to obtain an abundance distribution of NaCN. We found that the presence of NaCN seems compatible with the presence of CN, probably as a result of the photodissociation of HCN, in the inner layers of the ejecta of IRC +10216. However, similar as for CH 3 CN, current photochemical models fail to reproduce this CN reservoir. We also found that the abundance peak of NaCN appears at a radius of 3 x 10 15 cm, approximately where the abundance of NaCl, suggested to be the parent species, starts to decay. However, the abundance ratio shows that the NaCl abundance is lower than that obtained for NaCN. We expect that the LTE assumption might result in NaCN abundances higher than the real ones. Updated photochemical models, collisional rates, and reaction rates are essential to determine the possible paths of the NaCN formation., Comment: 7 pages, 10 figures. Accepted for publication in A&A letters

During their late pulsating phase, AGB stars expel most of their mass in the form of massive dusty envelopes, an event that largely controls the composition of interstellar matter. The envelopes, however, are distant and opaque to visible and NIR radiation: Their structure remains poorly known and the mass-loss process poorly understood. Millimeter-wave interferometry is the optimal investigative tool for this purpose. The circumstellar envelope IRC +10 216 and its central star, the C-rich TP-AGB star closest to the Sun, are the best objects for such an investigation. Two years ago, we reported on IRAM 30-m telescope CO(2-1) line emission observations in that envelope (HPBW 11"). We now report much higher angular resolution observations of CO(2-1), CO(1-0), CN(2-1) and C$_4$H(24-23) made with the SMA, PdB and ALMA interferometers (with synthesized half-power beamwidths of 3", 1" and 0.3", respectively). Although the envelope appears more intricate at high resolution, its prevailing structure remains a pattern of thin, nearly concentric shells. Outside the small (r<0.3") dust formation zone, the gas appears to expand radially at a constant velocity, 14.5 km/s, with small turbulent motions. Based on that property, we have reconstructed the 3-D structure of the outer envelope and have derived the gas temperature and density radial profiles in the inner (r<25") envelope. The over-dense shells have spherical or slightly oblate shapes and typically extend over a few steradians, implying isotropic mass loss. The regular spacing of shells in the outer envelope supports the model of a binary star system with a period of 700 years and a near face-on elliptical orbit. The companion fly-by triggers enhanced episodes of mass loss near periastron. The densification of the shell pattern observed in the central part of the envelope suggests a more complex scenario for the last few thousand years., Comment: 22 pages, 28 figures. Accepted for publication in A&A

Linear carbon chains are common in various types of astronomical molecular sources. Possible formation mechanisms involve both bottom-up and top-down routes. We have carried out a combined observational and modeling study of the formation of carbon chains in the C-star envelope IRC+10216, where the polymerization of acetylene and hydrogen cyanide induced by ultraviolet photons can drive the formation of linear carbon chains of increasing length. We have used ALMA to map the emission of 3 mm rotational lines of the hydrocarbon radicals C2H, C4H, and C6H, and the CN-containing species CN, C3N, HC3N, and HC5N with an angular resolution of 1". The spatial distribution of all these species is a hollow, 5-10" wide, spherical shell located at a radius of 10-20" from the star, with no appreciable emission close to the star. Our observations resolve the broad shell of carbon chains into thinner sub-shells which are 1-2" wide and not fully concentric, indicating that the mass loss process has been discontinuous and not fully isotropic. The radial distributions of the species mapped reveal subtle differences: while the hydrocarbon radicals have very similar radial distributions, the CN-containing species show more diverse distributions, with HC3N appearing earlier in the expansion and the radical CN extending later than the rest of the species. The observed morphology can be rationalized by a chemical model in which the growth of polyynes is mainly produced by rapid gas-phase chemical reactions of C2H and C4H radicals with unsaturated hydrocarbons, while cyanopolyynes are mainly formed from polyynes in gas-phase reactions with CN and C3N radicals., Comment: Accepted for publication in A&A

We carried out a millimeter wavelength line survey between 79 and 356 GHz with the IRAM-30m telescope to investigate the physical and chemical properties of the molecular envelope of the oxygen-rich AGB star IK Tau. We analysed the molecular lines detected using the population diagram technique to derive rotational temperatures and column densities. Additionally, we conducted a radiative transfer analysis of the SO$_2$ lines detected. For the first time in this source we detected rotational lines in the ground vibrational state of HCO$^+$, NS, NO, and H$_2$CO, as well as several isotopologues of molecules previously identified. We also detected several rotational lines in vibrationally excited states of SiS and SiO isotopologues, and rotational lines of H$_2$O ($\nu_{\mathrm{2}}$=2). We have also increased the number of rotational lines detected of molecules that were previously identified, enabling a detailed study of the molecular abundances and excitation temperatures. IK Tau displays a rich chemistry for an oxygen-rich circumstellar envelope. We highlight the detection of NS and H$_2$CO with fractional abundances of f(NS)$\sim$10$^{-8}$ and f(H$_2$CO)$\sim$[10$^{-7}$--10$^{-8}$]. Most of the molecules display rotational temperatures between 15 and 40 K. NaCl and SiS isotopologues display rotational temperatures higher than the average (~65 K). In the case of SO$_2$ a warm component with $T_{\mathrm{rot}}$$\sim$290 K is also detected. This SO$_2$ warm component is probably arising from the inner regions of the envelope (at $\sim$8$R_{*}$) where SO$_2$ has a fractional abundance of f(SO$_2$)$\sim$10$^{-6}$. This result should be considered for future investigation of the main formation channels of this, and other, parent species in the inner winds of O-rich AGB stars, which at present are not well reproduced by current chemistry models., Comment: 47 pages, 26 figures. Accepted for Astronomy and Astrophysics

We present high spectral resolution mid-IR observations of SiS towards the C-rich AGB star IRC+10216 carried out with the Texas Echelon-cross-Echelle Spectrograph mounted on the NASA Infrared Telescope Facility. We have identified 204 ro-vibrational lines of 28Si32S, 26 of 29Si32S, 20 of 28Si34S, and 15 of 30Si32S in the frequency range 720-790 cm-1. These lines belong to bands v=1-0, 2-1, 3-2, 4-3, and 5-4, and involve rotational levels with Jlow<90. About 30 per cent of these lines are unblended or weakly blended and can be partially or entirely fitted with a code developed to model the mid-IR emission of a spherically symmetric circumstellar envelope composed of expanding gas and dust. The observed lines trace the envelope at distances to the star <35R* (~0.7 arcsec). The fits are compatible with an expansion velocity of 1+2.5(r/R*-1) km/s between 1 and 5R*, 11 km/s between 5 and 20R*, and 14.5 km/s outwards. The derived abundance profile of 28Si32S with respect to H2 is 4.9e-6 between the stellar photosphere and 5R*, decreasing linearly to 1.6e-6 at 20R* and to 1.3e-6 at 50R*. 28Si32S seems to be rotationally under LTE in the region of the envelope probed with our observations and vibrationally out of LTE in most of it. There is a red-shifted emission excess in the 28Si32S lines of band v=1-0 that cannot be found in the lines of bands v=2-1, 3-2, 4-3, and 5-4. This excess could be explained by an enhancement of the vibrational temperature around 20R* behind the star. The derived isotopic ratios 28Si/29Si, and 32S/34S are 17 and 14, compatible with previous estimates., Comment: 11 pages, 5 figures, and 4 tables. Accepted for publication in MNRAS

IRC+10216 is a circumstellar envelope around a carbon-rich evolved star which contains a large variety of molecules. According to interferometric observations, molecules are distributed either concentrated around the central star or as a hollow shell with a radius of 15". We present ALMA Cycle 0 band 6 observations of the J=14-13 rotational transition of CH3CN in IRC+10216, obtained with an angular resolution of 0.76x0.61". The bulk of the emission is distributed as a hollow shell located at just 2" from the star, with a void of emission in the central region up to a radius of 1". This spatial distribution is markedly different from those found to date in this source for other molecules. Our analysis indicate that methyl cyanide is not formed neither in the stellar photosphere nor far in the outer envelope, but at radial distances as short as 1-2", reaching a maximum abundance of 0.02 molecules cm-3 at 2" from the star. Standard chemical models of IRC+10216 predict that the bulk of CH3CN molecules should be present at a radius of 15", where other species such as polyyne radicals and cyanopolyynes are observed, with an additional inner component within 1" from the star. The non-uniform structure of the circumstellar envelope and grain surface processes are discussed as possible causes of the peculiar distribution of methyl cyanide in IRC+10216., Comment: Accepted for publication in ApJ. Revised version includes updates following the recommendations of our referee

We report the discovery in space of a disilicon species, SiCSi, from observations between 80 and 350 GHz with the IRAM 30m radio telescope. Owing to the close coordination between laboratory experiments and astrophysics, 112 lines have now been detected in the carbon-rich star CWLeo. The derived frequencies yield improved rotational and centrifugal distortion constants up to sixth order. From the line profiles and interferometric maps with the Submillimeter Array, the bulk of the SiCSi emis- sion arises from a region of 6 arcseconds in radius. The derived abundance is comparable to that of SiC2. As expected from chemical equilibrium calculations, SiCSi and SiC2 are the most abundant species harboring a SiC bond in the dust formation zone and certainly both play a key role in the formation of SiC dust grains., Comment: To be published in the Astrophysical Journal Letters; Accepted May 6 2015

We report the detection during the ALMA Cycle 0 of SiS rotational lines in high-vibrational states as well as SiO and SiC$_2$ lines in their ground vibrational state, towards IRC+10216. The spatial distribution of these molecules shows compact emission for SiS and a more extended emission for SiO and SiC$_2$, and also proves the existence of an increase in the SiC$_2$ emission at the outer shells of the circumstellar envelope. We analyze the excitation conditions of the vibrationally excited SiS using the population diagram technique and we used a large velocity gradient model to compare with the observations. We found moderate discrepancies between the observations and the models that could be explained if SiS lines detected are optically thick. Additionally, the line profiles of the detected rotational lines in the high energy vibrational states show a decreasing linewidth with increasing energy levels. This may evidence that these lines could be excited only in the inner shells, i.e. the densest and hottest, of the circumstellar envelope of IRC+10216., Comment: 20 pages, 3 figures. Published in ApJL

OH 231.8+4.2, a bipolar outflow around a Mira-type variable star, displays a unique molecular richness amongst circumstellar envelopes (CSEs) around O-rich AGB and post-AGB stars. We report line observations of the HCO+ and H13CO+ molecular ions and the first detection of SO+, N2H+, and (tentatively) H3O+ in this source. SO+ and H3O+ have not been detected before in CSEs around evolved stars. These data have been obtained as part of a full mm-wave and far-IR spectral line survey carried out with the IRAM 30 m radio telescope and with Herschel/HIFI. Except for H3O+, all the molecular ions detected in this work display emission lines with broad profiles (FWHM 50-90 km/s), which indicates that these ions are abundant in the fast bipolar outflow of OH 231.8. The narrow profile (FWHM 14 km/s) and high critical densities (>1e6cm-3 ) of the H3O+ transitions observed are consistent with this ion arising from denser, inner (and presumably warmer) layers of the fossil remnant of the slow AGB CSE at the core of the nebula. From rotational diagram analysis, we deduce excitation temperatures of Tex 10-20 K for all ions except for H3O+, which is most consistent with Tex 100 K. Although uncertain, the higher excitation temperature suspected for H3O+ is similar to that recently found for H2O and a few other molecules, which selectively trace a previously unidentified, warm nebular component.The column densities of the molecular ions reported here are in the range Ntot [1-8]x1e13 cm-2, leading to beam-averaged fractional abundances relative to H2 of X(HCO+) 1e-8, X(H13CO+) 2e-9, X(SO+) 4e-9, X(N2H+) 2e-9, and X(H3O+) 7e-9 cm-2. We have performed chemical kinetics models to investigate the formation of these ions in OH 231.8 as the result of standard gas phase reactions initiated by cosmic-ray and UV-photon ionization. (abridged)., Comment: 24 pages, 11 figures

OH231.8+4.2 is a well studied O-rich CSE around an intermediate-mass evolved star that displays bipolar molecular outflows accelerated up to 400 km/s. OH231 also presents an exceptional molecular richness probably due to shock-induced chemical processes. We report the first detection in this source of HNCO, HNCS, HC3N and NO, with the IRAM-30m telescope in a sensitive mm-wavelength survey towards this target. HNCO and HNCS are also first detections in CSEs. The observed line profiles show that the emission arises in the massive central component of the envelope and at the base of the fast lobes. The NO profiles are broader than those of HNCO, HNCS, and HC3N, and most importantly, broader than the line profiles of 13CO. This indicates that the NO abundance is enhanced in the fast lobes relative to the slow, central parts. From LTE and non-LTE excitation analysis, we estimate beam-average rotational temperatures of 15-30 K (and, maybe, up to 55 K for HC3N) and fractional abundances of X(HNCO)=[0.8-1]E-7, X(HNCS)=[0.9-1]E-8, X(HC3N)=[5-7]E-9 and X(NO)=[1-2]E-6. NO is, therefore, amongst the most abundant N-bearing species in OH231. We have performed thermodynamical chemical equilibrium and chemical kinetics models to investigate the formation of these N-bearing species in OH231. The model underestimates the observed abundances for HNCO, HNCS, and HC3N by several orders of magnitude, which indicates that these molecules can hardly be products of standard UV-photon and/or cosmic-ray induced chemistry in OH231, and that other processes (e.g. shocks) play a major role in their formation. For NO, the model abundance is compatible with the observed average value. The new detections presented in this work corroborate the particularly rich chemistry of OH231, which is likely profoundly influenced by shock-induced processes, as proposed in earlier works., Comment: 22 pages, 14 figures. Published in A&A

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.