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

The rotational spectrum of silyl isocyanide (SiH$_3$NC), an isomer of the well studied silyl cyanide (SiH$_3$CN), has been detected in the laboratory in a supersonic molecular beam, and the identification was confirmed by observations of the corresponding rotational transitions in the rare isotopic species SiH$_3$$^{15}$NC and SiH$_3$N$^{13}$C. Spectroscopic constants derived from 19 transitions between $11 - 35$~GHz in the three lowest harmonically related rotational transitions in the $K = 0 ~{\rm{and}}~1$ ladders of the normal isotopic species including the nitrogen nuclear quadrupole hyperfine constant, allow the principal astronomical transitions of SiH$_3$NC to be calculated to an uncertainty of about 4~km~s$^{-1}$ in equivalent radial velocity, or within the FWHM of narrow spectral features in the inner region of IRC+10216 near 200~GHz. The concentration of SiH$_3$NC in our molecular beam is three times less than SiH$_3$CN, or about the same as the corresponding ratio of the isomeric pair SiNC and SiCN produced under similar conditions. Silyl isocyanide is an excellent candidate for astronomical detection, because the spectroscopic and chemical properties are very similar to SiH$_3$CN which was recently identified in the circumstellar envelope of IRC+10216 by \citet{cernicharo_discovery_2017} and of SiNC and SiCN in the same source.

We present a microwave spectral taxonomy study of several hydrocarbon/CS$_2$ discharge mixtures in which more than 60 distinct chemical species, their more abundant isotopic species, and/or their vibrationally excited states were detected using chirped-pulse and cavity Fourier-transform microwave spectroscopies. Taken together, in excess of 85 unique variants were detected, including several new isotopic species and more than 25 new vibrationally excited states of C$_2$S, C$_3$S, and C$_4$S, which have been assigned on the basis of published vibration-rotation interaction constants for C$_3$S, or newly calculated ones for C$_2$S and C$_4$S. On the basis of these precise, low-frequency measurements, several vibrationally exited states of C$_2$S and C$_3$S were subsequently identified in archival millimeter-wave data in the 253--280 GHz frequency range, ultimately providing highly accurate catalogs for astronomical searches. As part of this work, formation pathways of the two smaller carbon-sulfur chains were investigated using $^{13}$C isotopic spectroscopy, as was their vibrational excitation. The present study illustrates the utility of microwave spectral taxonomy as a tool for complex mixture analysis, and as a powerful and convenient `stepping stone' to higher frequency measurements in the millimeter and submillimeter bands.

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

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.

The formation of silicate dust in oxygen-rich envelopes of evolved stars is thought to be initiated by formation of seed particles that can withstand the high temperatures close to the stellar photosphere and act as condensation cores farther away from the star. Candidate species considered as first condensates are TiO and TiO$_2$. We aim to identify and characterize the circumstellar gas-phase chemistry of titanium that leads to the formation of solid titanium compounds in the envelope of $o$ Cet, the prototypical Mira, and seek an observational verification of whether titanium oxides play a major role in the onset of dust formation in M-type AGB stars. We present high angular-resolution ALMA observations at submillimeter (submm) wavelengths supplemented by APEX and Herschel spectra of the rotational features of TiO and TiO$_2$. In addition, circumstellar features of TiO and TiI are identified in optical spectra which cover multiple pulsation cycles of $o$ Cet. The submm ALMA data reveal TiO and TiO$_2$ bearing gas within the extended atmosphere of Mira. While TiO is traceable up to a radius (FWHM/2) of 4.0 R$_{\star}$, TiO$_2$ extends as far as 5.5 R$_{\star}$ and unlike TiO appears to be anisotropically distributed. Optical spectra display variable emission of TiI and TiO from inner parts of the extended atmosphere (<3 R$_{\star}$). Chemical models which include shocks are in general agreement with the observations of gas-phase titanium-bearing molecules. It is unlikely that substantial amounts of titanium is locked up in solids because the abundance of the gaseous titanium species is very high. In particular, formation of hot titanium-rich condensates is very improbable because we find no traces of their hot precursor species in the gas phase. It therefore appears unlikely that the formation of dust in Mira, and possibly other M-type AGB stars, is initiated by titanium oxides.

Context: Dust is efficiently produced by cool giant stars, but the condensation of inorganic dust is poorly understood. Aims: Identify and characterize aluminum bearing species in the circumstellar gas of Mira ($o$ Ceti) in order to elucidate their role in the production of Al$_2$O$_3$ dust. Methods: Multiepoch spectral line observations at (sub-)millimeter, far-infrared, and optical wavelengths including: maps with ALMA which probe the gas distribution in the immediate vicinity of the star at ~30 mas; observations with ALMA, APEX, and Herschel in 2013-2015 for studying cycle and inter-cycle variability of the rotational lines of Al bearing molecules; optical records as far back as 1965 to examine variations in electronic transitions over time spans of days to decades; and velocity measurements and excitation analysis of the spectral features which constrain the physical parameters of the gas. Results: Three diatomic molecules AlO, AlOH, and AlH, and atomic Al I are the main observable aluminum species in Mira, although a significant fraction of aluminum might reside in other species that have not yet been identified. Strong irregular variability in the (sub-)millimeter and optical features of AlO (possibly the direct precursor of Al$_2$O$_3$) indicates substantial changes in the excitation conditions, or varying abundance that is likely related to shocks in the star. The inhomogeneous distribution of AlO might influence the spatial and temporal characteristics of dust production. Conclusions: We are unable to quantitatively trace aluminum depletion from the gas, but the rich observational material constrains time dependent chemical networks. Future improvements should include spectroscopic characterization of higher aluminum oxides, coordinated observations of dust and gas species at different variability phases, and tools to derive abundances in shock excited gas.

The partially deuterated linear isomer HDCCC of the ubiquitous cyclic carbene ($c$-C$_3$H$_2$) was observed in the starless cores TMC-1C and L1544 at 96.9 GHz, and a confirming line was observed in TMC-1 at 19.38 GHz. To aid the identification in these narrow line sources, four centimetre-wave rotational transitions (two in the previously reported $K_a =0$ ladder, and two new ones in the $K_a =1$ ladder), and 23 transitions in the millimetre band between 96 and 272 GHz were measured in high-resolution laboratory spectra. Nine spectroscopic constants in a standard asymmetric top Hamiltonian allow the principal transitions of astronomical interest in the $K_a \le 3$ rotational ladders to be calculated to within 0.1 km s$^{-1}$ in radial velocity up to 400 GHz. Conclusive evidence for the identification of the two astronomical lines of HDCCC was provided by the $V_{\rm{LSR}}$ which is the same as that of the normal isotopic species (H$_2$CCC) in the three narrow line sources. In these sources, deuterium fractionation in singly substituted H$_2$CCC (HDCCC/H$_2$CCC $\sim4\%\text{-}19\%$) is comparable to that in $c$-C$_3$H$_2$ ($c$-C$_3$H$_2$/$c$-C$_3$HD $\sim5\%\text{-}17\%$), and similarly in doubly deuterated $c$-C$_3$H$_2$ ($c$-C$_3$D$_2$/$c$-C$_3$HD $\sim3\%\text{-}17\%$), implying that the efficiency of the deuteration processes in the H$_2$CCC and $c$-C$_3$H$_2$ isomers are comparable in dark clouds.

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.

Rotational spectra in four new excited vibrational levels of the linear carbon chain radical C$_4$H radical were observed in the millimeter band between 69 and 364 GHz in a low pressure glow discharge, and two of these were observed in a supersonic molecular beam between 19 and 38 GHz. All have rotational constants within 0.4% of the $^2\Sigma^+$ ground vibrational state of C$_4$H and were assigned to new bending vibrational levels, two each with $^2\Sigma$ and $^2\Pi$ vibrational symmetry. The new levels are tentatively assigned to the $1\nu_6$ and $1\nu_5$ bending vibrational modes (both with $^2\Pi$ symmetry), and the $1\nu_6 + 1\nu_7$ and $1\nu_5 + 1\nu_6$ combination levels ($^2\Sigma$ symmetry) on the basis of the derived spectroscopic constants, relative intensities in our discharge source, and published laser spectroscopic and quantum calculations. Prior spectroscopic constants in the $1\nu_7$ and $2\nu_7$ levels were refined. Also presented are interferometric maps of the ground state and the $1\nu_7$ level obtained with the SMA near 257 GHz which show that C$_4$H is present near the central star in IRC+10216. We found no evidence with the SMA for the new vibrationally excited levels of C$_4$H at a peak flux density averaged over a $3^{\prime\prime}$ synthesized beam of $\ge 0.15$ Jy/beam in the 294-296 and 304-306 GHz range, but it is anticipated that rotational lines in the new levels might be observed in IRC+10216 when ALMA attains its full design capability.

A spectral line survey of the oxygen-rich red supergiant VY Canis Majoris was made between 279 and 355 GHz with the Submillimeter Array. Two hundred twenty three spectral features from 19 molecules (not counting isotopic species of some of them) were observed, including the rotational spectra of TiO, TiO2, and AlCl for the first time in this source. The parameters and an atlas of all spectral features is presented. Observations of each line with a synthesized beam of ~0.9 arcsec, reveal the complex kinematics and morphology of the nebula surrounding VY CMa. Many of the molecules are observed in high lying rotational levels or in excited vibrational levels. From these, it was established that the main source of the submillimeter-wave continuum (dust) and the high excitation molecular gas (the star) are separated by about 0.15 arcsec. Apparent coincidences between the molecular gas observed with the SMA, and some of the arcs and knots observed at infrared wavelengths and in the optical scattered light by the Hubble Space Telescope are identified. The observations presented here provide important constraints on the molecular chemistry in oxygen-dominated circumstellar environments and a deeper picture of the complex circumstellar environment of VY CMa.

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

The formation of inorganic dust in circumstellar environments of evolved stars is poorly understood. Spectra of molecules thought to be most important for the nucleation, i.e. AlO, TiO, and TiO2, have been recently detected in the red supergiant VY CMa. These molecules are effectively formed in VY CMa and the observations suggest that non-equilibrium chemistry must be involved in their formation and nucleation into dust. In addition to exploring the recent observations of VY CMa, we briefly discuss the possibility of detecting these molecules in the dust-poor circumstellar environment of Betelgeuse.

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

We report a spectral-line survey of the extreme carbon star IRC+10216 carried out between 293.9 and 354.8 GHz with the Submillimeter Array. A total of 442 lines were detected, more than 200 for the first time; 149 are unassigned. Maps at an angular resolution of ~3" were obtained for each line. A substantial new population of narrow lines with an expansion velocity of ~4 km/s (i.e. ~ 30% of the terminal velocity) was detected. Most of these are attributed to rotational transitions within vibrationally excited states, emitted from energy levels above the v=0, J=0 ground state with excitation energy of 1000-3000 K. Emission from these lines appears to be centered on the star with an angular extent of <1". We use multiple transitions detected in several molecules to derive physical conditions in this inner envelope of IRC+10216.

We present the first astronomical detection of a diatomic negative ion, the cyanide anion CN-, as well as quantum mechanical calculations of the excitation of this anion through collisions with para-H2. CN- is identified through the observation of the J = 2-1 and J = 3-2 rotational transitions in the C-star envelope IRC +10216 with the IRAM 30-m telescope. The U-shaped line profiles indicate that CN-, like the large anion C6H-, is formed in the outer regions of the envelope. Chemical and excitation model calculations suggest that this species forms from the reaction of large carbon anions with N atoms, rather than from the radiative attachment of an electron to CN, as is the case for large molecular anions. The unexpectedly large abundance derived for CN-, 0.25 % relative to CN, makes likely its detection in other astronomical sources. A parallel search for the small anion C2H- remains so far unconclusive, despite the previous tentative identification of the J = 1-0 rotational transition. The abundance of C2H- in IRC +10216 is found to be vanishingly small, < 0.0014 % relative to C2H.

A spectral-line survey of IRC+10216 in the 345 GHz band has been undertaken with the Submillimeter Array. Although not yet completed, it has already yielded a fairly large sample of narrow molecular emission lines with line-widths indicating expansion velocities of ~4 km/s, less than 3 times the well-known value of the terminal expansion velocity (14.5 km/s) of the outer envelope. Five of these narrow lines have now been identified as rotational transitions in vibrationally excited states of previously detected molecules: the v=1, J=17--16 and J=19--18 lines of Si34S and 29SiS and the v=2, J=7--6 line of CS. Maps of these lines show that the emission is confined to a region within ~60 AU of the star, indicating that the narrow-line emission is probing the region of dust-formation where the stellar wind is still being accelerated.

Pure rotational spectra of three polycyclic aromatic hydrocarbons - acenaphthene, acenaphthylene and fluorene - have been obtained by Fourier transform microwave spectroscopy of a molecular beam and subsequently by millimeter wave absorption spectroscopy for acenaphthene and fluorene. The data presented here will be useful for deep radio astronomical searches for PAHs employing large radio telecopes.

Two consecutive rotational transitions of the long cyanopolyyne HC11N, J=39-38, and J=38-37, have been detected in the cold dust cloud TMC-1 at the frequencies expected from recent laboratory measurements by Travers et al. (1996), and at about the expected intensities. The astronomical lines have a mean radial velocity of 5.8(1) km/s, in good agreement with the shorter cyanopolyynes HC7N and HC9N observed in this very sharp-lined source [5.82(5) and 5.83(5) km/s, respectively]. The column density of HC11N is calculated to be 2.8x10^(11) cm^(-2). The abundance of the cyanopolyynes decreases smoothly with length to HC11N, the decrement from one to the next being about 6 for the longer carbon chains.

The cumulene carbenes are important components of hydrocarbon chemistry in low mass star forming cores. Here we report the first astronomical detection of the long chain cumulene carbene H_2C_6 in the interstellar cloud TMC1, from observations of two of its rotational transitions: $J_{K,K'} = 7_{1,7} \rightarrow 6_{1,6}$ at 18.8 GHz and $8_{1,8} \rightarrow 7_{1,7}$ at 21.5 GHz, using NASA's Deep Space Network 70 m antenna at Goldstone, California. In addition we also observed the shorter cumulene carbene, H_2C_4 at the same position. The fractional abundance of H_2C_6 relative to H_2 is about $4.7 \times 10^{-11}$ and H_2C_4 is about $1.1 \times 10^{-9}$. The abundance of H_2C_6 is in fairly good agreement with gas phase chemical models for young molecular cloud cores, but the abundance of H_2C_4 is significantly larger than predicted.