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.

Interstellar carbonaceous dust is mainly formed in the innermost regions of circumstellar envelopes around carbon-rich asymptotic giant branch (AGB) stars. In these highly chemically stratified regions, atomic and diatomic carbon, along with acetylene are the most abundant species after H$_2$ and CO. In a previous study, we addressed the chemistry of carbon (C and C$_2$) with H$_2$ showing that acetylene and aliphatic species form efficiently in the dust formation region of carbon-rich AGBs whereas aromatics do not. Still, acetylene is known to be a key ingredient in the formation of linear polyacetylenic chains, benzene and polycyclic aromatic hydrocarbons (PAHs), as shown by previous experiments. However, these experiments have not considered the chemistry of carbon (C and C$_2$) with C$_2$H$_2$. In this work, by employing a sufficient amount of acetylene, we investigate its gas-phase interaction with atomic and diatomic carbon. We show that the chemistry involved produces linear polyacetylenic chains, benzene and other PAHs, which are observed with high abundances in the early evolutionary phase of planetary nebulae. More importantly, we have found a non-negligible amount of pure and hydrogenated carbon clusters as well as aromatics with aliphatic substitutions, both being a direct consequence of the addition of atomic carbon. The incorporation of alkyl substituents into aromatics can be rationalized by a mechanism involving hydrogen abstraction followed by methyl addition. All the species detected in gas phase are incorporated into the nanometric sized dust analogues, which consist of a complex mixture of sp, sp$^2$ and sp$^3$ hydrocarbons with amorphous morphology.

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.

Aim: We aim to determine the abundances of CS, SiO, and SiS in a large sample of carbon star 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 the surroundings of the central AGB star. Methods: We surveyed a sample of 25 carbon-rich AGB stars in the $\lambda$ 2 mm band, using the IRAM 30 m 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 CSEs. Results: We detected CS in all 25 CSEs, SiO in 24 of them, and SiS in 17 sources. We found that CS and SiS have similar abundances in carbon star envelopes, while SiO is present with a lower abundance. We also found a strong correlation in which the denser the envelope, the less abundant are CS and SiO. The trend is however only tentatively seen for SiS in the range of high mass loss rates. Furthermore, we found a relation in which the integrated flux of the MgS dust feature at 30 um increases as the fractional abundance of CS decreases. Conclusions: The decline in the fractional abundance of CS with increasing density could be due to gas-phase chemistry in the inner envelope or to adsorption onto dust grains. The latter possibility is favored by a correlation between the CS fractional abundance and the 30 um feature, which suggests that CS is efficiently incorporated onto MgS dust around C-rich AGB stars. In the case of SiO, the observed abundance depletion with increasing density is most likely caused by an efficient incorporation onto dust grains. We conclude that CS, SiO (very likely), and SiS (tentatively) are good candidates to act as gas-phase precursors of dust in C-rich AGB envelopes.

AB Aur is a Herbig Ae star hosting a well-known transitional disk. Because of its proximity and low inclination angle, it is an excellent object to study planet formation. Our goal is to investigate the chemistry and dynamics of the molecular gas component in the AB Aur disk, and its relation with the prominent horseshoe shape observed in continuum mm emission. We used the NOEMA interferometer to map with high angular resolution the J = 3-2 lines of HCO+ and HCN. By combining both, we can gain insight into the AB Aur disk structure. Chemical segregation is observed in the AB Aur disk: HCO+ shows intense emission toward the star position, at least one bright molecular bridge within the dust cavity, and ring-like emission at larger radii, while HCN is only detected in an annular ring that is coincident with the dust ring and presents an intense peak close to the dust trap. We use HCO+ to investigate the gas dynamics inside the cavity. The observed bright HCO+ bridge connects the compact central source with the outer dusty ring. This bridge can be interpreted as an accretion flow from the outer ring to the inner disk/jet system proving gas accretion through the cavity.

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.

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

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.

The presence of SiS in space seems to be restricted to a few selected types of astronomical environments. It is long known to be present in circumstellar envelopes around evolved stars and it has also been detected in a handful of star-forming regions with evidence of outflows, like Sgr\,B2, Orion\,KL and more recently L1157-B1. The kinetics of reactions involving SiS is very poorly known and here we revisit the chemistry of SiS in space by studying some potentially important reactions of formation and destruction of this molecule. We calculated {\it ab initio} potential energy surfaces of the SiOS system and computed rate coefficients in the temperature range 50-2500 K for the reaction of destruction of SiS, in collisions with atomic O, and of its formation, through the reaction between Si and SO. We find that both reactions are rapid, with rate coefficients of a few times 10$^{-10}$ cm$^3$ s$^{-1}$, almost independent of temperature. In the reaction between Si and SO, SiO production is 5-7 times more efficient than SiS formation. The reaction of SiS with O atoms can play an important role in destroying SiS in envelopes around evolved stars. We built a simple chemical model of a postshock gas to study the chemistry of SiS in protostellar outflows and we found that SiS forms with a lower abundance and later than SiO, that SiS is efficiently destroyed through reaction with O, and that the main SiS-forming reactions are Si + SO and Si + SO$_2$.

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.

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.

Recent Herschel and ALMA observations of Photodissociation Regions (PDRs) have revealed the presence of a high thermal pressure (P ~ 10^7-10^8 K cm-3) thin compressed layer at the PDR surface where warm molecular tracer emission (e.g. CH+, SH+, high-J CO, H2,...) originate. These high pressures (unbalanced by the surrounding environment) and a correlation between pressure and incident FUV field (G0) seem to indicate a dynamical origin with the radiation field playing an important role in driving the dynamics. We investigate whether photoevaporation of the illuminated edge of a molecular cloud could explain these high pressures and pressure-UV field correlation. We developed a 1D hydrodynamical PDR code coupling hydrodynamics, EUV and FUV radiative transfer and time-dependent thermo-chemical evolution. We applied it to a 1D plane-parallel photoevaporation scenario where a UV-illuminated molecular cloud can freely evaporate in a surrounding low-pressure medium. We find that photoevaporation can produce high thermal pressures and the observed P-G0 correlation, almost independently from the initial gas density. In addition, we find that constant-pressure PDR models are a better approximation to the structure of photoevaporating PDRs than constant-density PDR models, although moderate pressure gradients are present. Strong density gradients from the molecular to the neutral atomic region are found, which naturally explain the large density contrasts (1-2 orders of magnitude) derived from observations of different tracers. The photoevaporating PDR is preceded by a low velocity shock (a few km/s) propagating into the molecular cloud. Photoevaporating PDR models offer a promising explanation to the recent observational evidence of dynamical effects in PDRs.

We present the detection of C4H2 for first time in the envelope of the C-rich AGB star IRC+10216 based on high spectral resolution mid-IR observations carried out with the Texas Echelon-cross-Echelle Spectrograph (TEXES) mounted on the Infrared Telescope Facility (IRTF). The obtained spectrum contains 24 narrow absorption features above the detection limit identified as lines of the ro-vibrational C4H2 band nu6+nu8(sigma_u^+). The analysis of these lines through a ro-vibrational diagram indicates that the column density of C4H2 is 2.4(1.5)E+16 cm^(-2). Diacetylene is distributed in two excitation populations accounting for 20 and 80% of the total column density and with rotational temperatures of 47(7) and 420(120) K, respectively. This two-folded rotational temperature suggests that the absorbing gas is located beyond ~0.4"~20R* from the star with a noticeable cold contribution outwards from ~10"~500R*. This outer shell matches up with the place where cyanoacetylenes and carbon chains are known to form due to the action of the Galactic dissociating radiation field on the neutral gas coming from the inner layers of the envelope.

We present a proof of concept on the coupling of radio astronomical receivers and spectrometers with chemical reactorsand the performances of the resulting setup for spectroscopy and chemical simulations in laboratory astrophysics. Several experiments including cold plasma generation and UV photochemistry were performed in a 40\,cm long gas cell placed in the beam path of the Aries 40\,m radio telescope receivers operating in the 41-49 GHz frequency range interfaced with fast Fourier transform spectrometers providing 2 GHz bandwidth and 38 kHz resolution. The impedance matching of the cell windows has been studied using different materials. The choice of the material and its thickness was critical to obtain a sensitivity identical to that of standard radio astronomical observations. Spectroscopic signals arising from very low partial pressures of CH3OH, CH3CH2OH, HCOOH, OCS,CS, SO2 (<1E-03 mbar) were detected in a few seconds. Fast data acquisition was achieved allowing for kinetic measurements in fragmentation experiments using electron impact or UV irradiation. Time evolution of chemical reactions involving OCS, O2 and CS2 was also observed demonstrating that reactive species, such as CS, can be maintained with high abundance in the gas phase during these experiments.

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.

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.

The extremely young Class 0 object B1b-S and the first hydrostatic core (FSHC) candidate, B1b-N, provide a unique opportunity to study the chemical changes produced in the elusive transition from the prestellar core to the protostellar phase. We present 40"x70" images of Barnard 1b in the 13CO 1->0, C18O 1->0, NH2D 1_{1,1}a->1_{0,1}s, and SO 3_2->2_1 lines obtained with the NOEMA interferometer. The observed chemical segregation allows us to unveil the physical structure of this young protostellar system down to scales of ~500au. The two protostellar objects are embedded in an elongated condensation, with a velocity gradient of ~0.2-0.4 m s^{-1} au^{-1} in the east-west direction, reminiscent of an axial collapse. The NH2D data reveal cold and dense pseudo-disks (R~500-1000 au) around each protostar. Moreover, we observe evidence of pseudo-disk rotation around B1b-S. We do not see any signature of the bipolar outflows associated with B1b-N and B1b-S, which were previously detected in H2CO and CH3OH, in any of the imaged species. The non-detection of SO constrains the SO/CH3OH abundance ratio in the high-velocity gas.

One serious challenge for planet formation is the rapid inward drift of pebble-sized dust particles in protoplanetary disks. Dust trapping at local maxima in the disk gas pressure has received much theoretical attention but still lacks observational support. The cold dust emission in the AB Aur disk forms an asymmetric ring at a radius of about 120 au, which is suggestive of dust trapping in a gas vortex. We present high spatial resolution (0".58x0".78 ~ 80x110 au) NOEMA observations of the 1.12 mm and 2.22 mm dust continuum emission from the AB Aur disk. Significant azimuthal variations of the flux ratio at both wavelengths indicate a size segregation of the large dust particles along the ring. Our continuum images also show that the intensity variations along the ring are smaller at 2.22 mm than at 1.12 mm, contrary to what dust trapping models with a gas vortex have predicted. Our two-fluid (gas+dust) hydrodynamical simulations demonstrate that this feature is well explained if the gas vortex has started to decay due to turbulent diffusion, and dust particles are thus losing the azimuthal trapping on different timescales depending on their size. The comparison between our observations and simulations allows us to constrain the size distribution and the total mass of solid particles in the ring, which we find to be of the order of 30 Earth masses, enough to form future rocky planets.

Silicon carbide dust grains are 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. To date, only three molecules containing an Si--C bond have been identified to have significant abundances in C-rich AGB stars: SiC$_2$, SiC, and Si$_2$C. The ring molecule SiC$_2$ has been observed in a handful of evolved stars, while SiC and Si$_2$C have only been detected in the C-star envelope IRC +10216. We aim to study how widespread and abundant SiC$_2$, SiC, and Si$_2$C 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 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.

Our knowledge of the chemical properties of the circumstellar ejecta of the most massive evolved stars is particularly poor. We aim to study the chemical characteristics of the prototypical yellow hypergiant star, IRC +10420. For this purpose, we obtained full line surveys at 1 and 3mm atmospheric windows. We have identified 106 molecular emission lines from 22 molecular species. Approximately half of the molecules detected are N-bearing species, in particular HCN, HNC, CN, NO, NS, PN, and N2H+. We used rotational diagrams to derive the density and rotational temperature of the different molecular species detected. We introduced an iterative method that allows us to take moderate line opacities into account. We have found that IRC +10420 presents high abundances of the N-bearing molecules compared with O-rich evolved stars. This result supports the presence of a N-rich chemistry, expected for massive stars. Our analysis also suggests a decrease of the 12C/13C ratio from \gtrsim 7 to \sim 3.7 in the last 3800 years, which can be directly related to the nitrogen enrichment observed. In addition, we found that SiO emission presents a significant intensity decrease for high-J lines when compared with older observations. Radiative transfer modeling shows that this variation can be explained by a decrease in the infrared (IR) flux of the dust. The origin of this decrease might be an expansion of the dust shell or a lower stellar temperature due to the pulsation of the star.

Characterizing the atmospheres of extrasolar planets is the new frontier in exoplanetary science. The last two decades of exoplanet discoveries have revealed that exoplanets are very common and extremely diverse in their orbital and bulk properties. We now enter a new era as we begin to investigate the chemical diversity of exoplanets, their atmospheric and interior processes, and their formation conditions. Recent developments in the field have led to unprecedented advancements in our understanding of atmospheric chemistry of exoplanets and the implications for their formation conditions. We review these developments in the present work. We review in detail the theory of atmospheric chemistry in all classes of exoplanets discovered to date, from highly irradiated gas giants, ice giants, and super-Earths, to directly imaged giant planets at large orbital separations. We then review the observational detections of chemical species in exoplanetary atmospheres of these various types using different methods, including transit spectroscopy, doppler spectroscopy, and direct imaging. In addition to chemical detections, we discuss the advances in determining chemical abundances in these atmospheres and how such abundances are being used to constrain exoplanetary formation conditions and migration mechanisms. Finally, we review recent theoretical work on the atmospheres of habitable exoplanets, followed by a discussion of future outlook of the field.

Transitional disks are structures of dust and gas around young stars with large inner cavities in which planet formation may occur. Lopsided dust distributions are observed in the dust continuum emission at millimeter wavelengths. These asymmetrical structures can be explained as the result of an enhanced gas density vortex where the dust is trapped potentially promoting the rapid growth to the planetesimal scale. AB Aur hosts a transitional disk with a clear horseshoe morphology which strongly suggests the presence of a dust trap. Our goal is to investigate its formation and the possible effects on the gas chemistry. We used the NOEMA interferometer to image the 1mm continuum dust emission and the 13CO J=2->1, C18O J=2->1, SO J=56->45 and H2CO J=303->202 rotational lines. Line integrated intensity ratio images are built to investigate the chemical changes within the disk. We have used a single point (n,T) chemical model to investigate the lifetime of gaseous CO, H2CO and SO in the dust trap. Our model shows that for densities >10^7 cm^-3, the SO molecules are depleted (directly frozen or converted into SO2 and then frozen out) in less than 0.1 Myr. The lower SO abundance towards the dust trap could indicate that a larger fraction of the gas is in a high density environment. Gas dynamics, grain growth and gas chemistry are coupled in the planet formation process. Because of the strong dependence of SO abundance on the gas density, the sulfur chemistry can be used as a chemical diagnostic to detect the birthsites of future planets. However, the large uncertainties inherent to chemical models and the limited knowledge of the disk physical structure and initial conditions are important drawbacks.

Context. The recent detection of warm H$_2$O vapor emission from the outflows of carbon-rich asymptotic giant branch (AGB) stars challenges the current understanding of circumstellar chemistry. Two mechanisms have been invoked to explain warm H$_2$O vapor formation. In the first, periodic shocks passing through the medium immediately above the stellar surface lead to H$_2$O formation. In the second, penetration of ultraviolet interstellar radiation through a clumpy circumstellar medium leads to the formation of H$_2$O molecules in the intermediate wind. Aims. We aim to determine the properties of H$_2$O emission for a sample of 18 carbon-rich AGB stars and subsequently constrain which of the above mechanisms provides the most likely warm H$_2$O formation pathway. Methods, Results, and Conclusions. See paper.

We present the detection of two H2C3O isomers, propynal and cyclopropenone, toward various starless cores and molecular clouds, together with upper limits for the third isomer propadienone. We review the processes controlling the abundances of H2C3O isomers in interstellar media showing that the reactions involved are gas-phase ones. We show that the abundances of these species are controlled by kinetic rather than thermodynamic effects.

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.

Aims. Our goal is to determine the molecular composition of the circumstellar disk around AB Aurigae (hereafter, AB Aur). AB Aur is a prototypical Herbig Ae star and the understanding of its disk chemistry is of paramount importance to understand the chemical evolution of the gas in warm disks. Methods. We used the IRAM 30-m telescope to perform a sensitive search for molecular lines in AB Aur as part of the IRAM Large program ASAI (A Chemical Survey of Sun-like Star-forming Regions). These data were complemented with interferometric observations of the HCO+ 1-0 and C17O 1-0 lines using the IRAM Plateau de Bure Interferometer (PdBI). Single-dish and interferometric data were used to constrain chemical models. Results. Throughout the survey, several lines of CO and its isotopologues, HCO+, H2CO, HCN, CN and CS, were detected. In addition, we detected the SO 54-33 and 56-45 lines, confirming the previous tentative detection. Comparing to other T Tauri's and Herbig Ae disks, AB Aur presents low HCN 3-2/HCO+ 3-2 and CN 2-1/HCN 3-2 line intensity ratios, similar to other transition disks. AB Aur is the only protoplanetary disk detected in SO thus far. Conclusions. We modeled the line profiles using a chemical model and a radiative transfer 3D code. Our model assumes a flared disk in hydrostatic equilibrium. The best agreement with observations was obtained for a disk with a mass of 0.01 Msun , Rin=110 AU, Rout=550 AU, a surface density radial index of 1.5 and an inclination of 27 deg. The intensities and line profiles were reproduced within a factor of 2 for most lines. This agreement is reasonable taking into account the simplicity of our model that neglects any structure within the disk. However, the HCN 3-2 and CN 2-1 line intensities were predicted more intense by a factor of >10. We discuss several scenarios to explain this discrepancy.

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.

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

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

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

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.

Hydrides represent the first steps of interstellar chemistry. Sulfanylium (SH$^+$), in particular, is a key tracer of energetic processes. We used ALMA and the IRAM 30 m telescope to search for the lowest frequency rotational lines of SH$^+$ toward the Orion Bar, the prototypical photo-dissociation region illuminated by a strong UV radiation field. On the basis of previous $Herschel$/HIFI observations of SH$^+$, we expected to detect emission of the two SH$^+$ hyperfine structure (HFS) components of the $N_J = 1_0 - 0_1$ fine structure (FS) component near 346 GHz. While we did not observe any lines at the frequencies predicted from laboratory data, we detected two emission lines, each $\sim$15 MHz above the SH$^+$ predictions and with relative intensities and HFS splitting expected for SH$^+$. The rest frequencies of the two newly detected lines are more compatible with the remainder of the SH$^+$ laboratory data than the single line measured in the laboratory near 346 GHz and previously attributed to SH$^+$. Therefore, we assign these new features to the two SH$^+$ HFS components of the $N_J = 1_0 - 0_1$ FS component and re-determine its spectroscopic parameters, which will be useful for future observations of SH$^+$, in particular if its lowest frequency FS components are studied. Our observations demonstrate the suitability of these lines for SH$^+$ searches at frequencies easily accessible from the ground.

The C-star envelope IRC+10216 harbors a rich variety of molecules, with more than 80 detected to date. During the course of a $\lambda$ 3 mm survey of IRC+10216 carried out with the IRAM 30-m telescope we have detected various weak lines, with antenna temperatures of a few mK, that we assign to rotational transitions of four new molecules. The observation of three lines of C$_5$S confirms a previous tentative identification of this molecule by Bell et al. (1993) based on a line at 24.0 GHz. We also report the tentative identification of three molecules not yet observed in space: MgCCH, the first metal acetylide detected in space, and NCCP and SiH$_3$CN, the phosphorus and silicon analogs of cyanogen (NCCN) and methyl cyanide (CH$_3$CN). We derive the following column densities: $N$(C$_5$S) = (2-14) $\times$ 10$^{12}$ cm$^{-2}$ (depending on the rotational temperature adopted), $N$(MgCCH) = 2 $\times$ 10$^{12}$ cm$^{-2}$, $N$(NCCP) = 7 $\times$ 10$^{11}$ cm$^{-2}$, and $N$(SiH$_3$CN) = 10$^{12}$ cm$^{-2}$. The S-bearing carbon chain C$_5$S is less abundant than C$_3$S, while MgCCH has an abundance in between that of MgNC and those of MgCN and HMgNC. On the other hand, NCCP and SiH$_3$CN are the least abundant P- and Si-bearing molecules observed to date in IRC+10216. Based on the behavior of similar molecules it is likely that these four species are formed in the outer circumstellar layers of IRC+10216. We discuss possible gas-phase formation routes.

The past twenty years have revealed the diversity of planets that exist in the Universe. It turned out that most of exoplanets are different from the planets of our Solar System and thus, everything about them needs to be explored. Thanks to current observational technologies, we are able to determine some information about the atmospheric composition, the thermal structure and the dynamics of these exoplanets, but many questions remain still unanswered. To improve our knowledge about exoplanetary systems, more accurate observations are needed and that is why the Exoplanet Characterisation Observatory (EChO) is an essential space mission. Thanks to its large spectral coverage and high spectral resolution, EChO will provide exoplanetary spectra with an unprecedented accuracy, allowing to improve our understanding of exoplanets. In this work, we review what has been done to date concerning the chemical modeling of exoplanet atmospheres and what are the main characteristics of warm exoplanet atmospheres, which are one of the main targets of EChO. Finally we will present the ongoing developments that are necessary for the chemical modeling of exoplanet atmospheres.

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

The Plateau de Bure Interferometer has been used to map the continuum emission at 3.4 mm and 1.1 mm together with the J=1->0 and J=3->2 lines of HCN and HCO+ towards the binary star GV Tau. The 3.4 mm observations did not resolve the binary components and the HCN J=1->0 and HCO+ J=1->0 line emissions trace the circumbinary disk and the flattened envelope. However, the 1.1 mm observations resolved the individual disks of GV Tau N and GV Tau S and allowed us to study their chemistry. We detected the HCN 3->2 line only towards the individual disk of GV Tau N, and the emission of the HCO+ 3->2 line towards GV Tau S. Simple calculations indicate that the 3->2 line of HCN is formed in the inner R<12 AU of the disk around GV Tau N where the HCN/HCO+ abundance ratio is >300. On the contrary, this ratio is <1.6 in the disk around GV Tau S. The high HCN abundance measured in GV Tau N is well explained by photo-chemical processes in the warm (>400K) and dense disk surface.

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

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

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

We report the detection in the envelope of the C-rich star IRC +10216 of four series of lines with harmonically related frequencies: B1389, B1390, B1394 and B1401. The four series must arise from linear molecules with mass and size close to those of C6H and C5N. Three of the series have half-integer rotational quantum numbers; we assign them to the 2Delta and 2Sigma vibronic states of C6H in its lowest (v_11) bending mode. The fourth series, B1389, has integer J with no evidence of fine or hyperfine structure; it has a rotational constant of 1388.860(2) MHz and a centrifugal distortion constant of 33(1) Hz; it is almost certainly the C5N- anion.

The J,K = 1,0-0,0 rotational transition of phosphine (PH3) at 267 GHz has been tentatively identified with a T_MB = 40 mK spectral line observed with the IRAM 30-m telescope in the C-star envelope IRC+10216. A radiative transfer model has been used to fit the observed line profile. The derived PH3 abundance relative to H2 is 6 x 10^(-9), although it may have a large uncertainty due to the lack of knowledge about the spatial distribution of this species. If our identification is correct, it implies that PH3 has a similar abundance to that reported for HCP in this source, and that these two molecules (HCP and PH3) together take up about 5 % of phosphorus in IRC+10216. The abundance of PH3, as that of other hydrides in this source, is not well explained by conventional gas phase LTE and non-LTE chemical models, and may imply formation on grain surfaces.

We report on the first detection of maser emission in the J=11-10, J=14-13 and J=15-14 transitions of the v=0 vibrational state of SiS toward the C-rich star IRC+10216. These masers seem to be produced in the very inhomogeneous region between the star and the inner dust formation zone, placed at 5-7 R*, with expansion velocities below 10 km/s. We interpret the pumping mechanism as due to overlaps between v=1-0 ro-vibrational lines of SiS and mid-IR lines of C2H2, HCN and their 13C isotopologues. The large number of overlaps found suggests the existence of strong masers for high-J v=0 and v=1 SiS transitions, located in the submillimeter range. In addition, it could be possible to find several rotational lines of the SiS isotopologues displaying maser emission.