Atmospheric molecular blobs shape up circumstellar envelopes of AGB stars

During their thermally pulsing phase, asymptotic giant branch (AGB) stars eject material that forms extended dusty envelopes1. Visible polarimetric imaging found clumpy dust clouds within two stellar radii of several oxygen-rich stars2-6. Inhomogeneous molecular gas has also been observed in multiple emission lines within several stellar radii of different oxygen-rich stars, including W Hya and Mira7-10. At the stellar surface level, infrared images have shown intricate structures around the carbon semiregular variable R Scl and in the S-type star π1 Gru11,12. Infrared images have also shown clumpy dust structures within a few stellar radii of the prototypical carbon AGB star IRC+10°216 (refs. 13,14), and studies of molecular gas distribution beyond the dust formation zone have also shown complex circumstellar structures15. Because of the lack of sufficient spatial resolution, however, the distribution of molecular gas in the stellar atmosphere and the dust formation zone of AGB carbon stars is not known, nor is how it is subsequently expelled. Here we report observations with a resolution of one stellar radius of the recently formed dust and molecular gas in the atmosphere of IRC+10°216. Lines of HCN, SiS and SiC2 appear at different radii and in different clumps, which we interpret as large convective cells in the photosphere, as seen in Betelgeuse16. The convective cells coalesce with pulsation, causing anisotropies that, together with companions17,18, shape its circumstellar envelope.

Detection of vibrationally excited HC7N and HC9N in IRC+10216 ⋆

Observations of IRC +10216 with the Yebes 40m telescope between 31 and 50 GHz have revealed more than 150 unidentified lines. Some of them can be grouped into a new series of 26 doublets, harmonically related with integer quantum numbers ranging from J up=54 to 80. The separation of the doublets increases systematically with J, i.e., as expected for a linear species in one of its bending modes. The rotational parameters resulting from the fit to these data are B = 290.8844 ± 0.0004 MHz, D = 0.88 ± 0.04 Hz, q = 0.1463 ± 0.0001 MHz. The rotational constant is very close to that of the ground state of HC9N. Ab initio calculations show an excellent agreement between these parameters and those predicted for the lowest energy vibrationally excited state, ν 19=1, of HC9N. This is the first detection, and complete characterization in space, of vibrationally excited HC9N. An energy of 41.5 cm-1 is estimated for the ν 19 state. In addition, 17 doublets of HC7N in the ν 15=1 state, for which laboratory spectroscopy is available, have been detected for the first time in IRC+10216. Several doublets of HC5N in its ν 11=1 state have been also observed. The column density ratio between the ground and the lowest excited vibrational states are ≈127, 9.5, and 1.5 for HC5N, HC7N, and HC9N, respectively. We find that these lowest-lying vibrational states are most probably populated via infrared pumping to vibrationally excited states lying at ≈600 cm-1. The lowest vibrationally excited states thus need to be taken into account to precisely determine absolute abundances and abundanceratios for long carbon chains. The abundance ratios N(HC5N)/N(HC7N) and N(HC7N)/N(HC9N) are 2.4 and 7.7 respectively.

New molecular species at redshift z = 0.89

We present the first detections of CH₃SH, C₃H⁺, C₃N, HCOOH, CH₂CHCN, and H₂CN in an extragalactic source. Namely the spiral arm of a galaxy located at z = 0.89 on the line of sight to the radio-loud quasar PKS 1830-211. OCS, SO₂, and NH₂CN were also detected, raising the total number of molecular species identified in that early time galaxy to 54, not counting isotopologues. The detections were made in absorption against the SW quasar image, at 2 kpc from the galaxy centre, over the course of a Q band spectral line survey made with the Yebes 40 m telescope (rest-frame frequencies: 58.7-93.5 GHz). We derived the rotational temperatures and column densities of those species, which are found to be subthermally excited. The molecular abundances, and in particular the large abundances of C₃H⁺ and of several previously reported cations, are characteristic of diffuse or translucent clouds with enhanced UV radiation or strong shocks.

Discovery of two new magnesium-bearing species in IRC+10216: MgC3N and MgC4H

We report on the detection of two series of harmonically related doublets in IRC +10216. From the observed frequencies, the rotational constant of the first series is B = 1380.888 MHz and that of the second series is B = 1381.512 MHz. The two series correspond to two species with a 2Σ electronic ground state. After considering all possible candidates, and based on quantum chemical calculations, the first series is assigned to MgC3N and the second to MgC4H. For the latter species, optical spectroscopy measurements support its identification. Unlike diatomic metal-containing molecules, the line profiles of the two new molecules indicate that they are formed in the outer layers of the envelope, as occurs for MgNC and other polyatomic metal-cyanides. We also confirm the detection of MgCCH that was previously reported from the observation of two doublets. The relative abundance of MgC3N with respect to MgNC is close to one while that of MgC4H relative to MgCCH is about ten. The synthesis of these magnesium cyanides and acetylides in IRC +10216 can be explained in terms of a two-step process initiated by the radiative association of Mg+ with large cyanopolyynes and polyynes followed by the dissociative recombination of the ionic complexes.

Discovery of the first Ca-bearing molecule in space: CaNC

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.5)×1011 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 ≃1/60 and ≃1, respectively. We searched for the species CaF, CaCl, CaC, CaCCH, and CaCH3 for which accurate frequency predictions are available. Only upper limits have been obtained for these molecules.

IRC +10216 as a spectroscopic laboratory: improved rotational constants for SiC2, its isotopologues, and Si2C

This work presents a detailed analysis of the laboratory and astrophysical spectral data available for 28SiC2, 29SiC2,30SiC2, Si13CC, and Si2C. New data on the rotational lines of these species between 70 and 350 GHz have been obtained with high spectral resolution (195 kHz) with the IRAM 30m telescope in the direction of the circumstellar envelope IRC +10216. Frequency measurements can reach an accuracy of 50 kHz for features observed with a good signal to noise ratio. From the observed astrophysical lines and the available laboratory data new rotational and centrifugal distortion constants have been derived for all the isotopologues of SiC2, allowing to predict their spectrum with high accuracy in the millimeter and submillimeter domains. Improved rotational and centrifugal distortion constants have also been obtained for disilicon carbide, Si2C. This work shows that observations of IRC +10216 taken with the IRAM 30m telescope, with a spectral resolution of 195 kHz, can be used for any molecular species detected in this source to derive, or improve, its rotational constants. Hence, IRC +10216 in addition to be one the richest sources in molecular species in the sky, can also be used as a state-of-the-art spectroscopy laboratory in the millimeter and submillimeter domains.

Time-dependent molecular emission in IRC+10216

CONTEXT

The variability in IRC+10216, the envelope of the asymptotic giant branch (AGB) star CW Leo, has attracted increasing attention in recent years. Studying the details of this variability in the molecular emission required a systematic observation program.

AIMS

We aim to reveal and characterize the periodical variability of the rotational lines from several molecules and radicals in IRC+10216, and to compare it with previously reported IR variability.

METHODS

We carried out systematic monitoring within the ~80 to 116 GHz frequency range with the IRAM 30m telescope.

RESULTS

We report on the periodical variability in IRC+10216 of several rotational lines from the following molecules and radicals: HC3N, HC5N, CCH, C4H, C5H, and CN. The analysis of the variable molecular lines provides periods that are consistent with previously reported IR variability, and interesting phase lags are revealed that point toward radiative transfer and pumping, rather than chemical effects.

CONCLUSIONS

This study indicates that observations of several lines of a given molecule have to be performed simultaneously or at least at the same phase in order to avoid erroneous interpretation of the data. In particular, merging ALMA data from different epochs may prove to be difficult, as shown by the example of the variability we studied here. Moreover, radiative transfer codes have to incorporate the effect of population variability in the rotational levels in CW Leo.

IRC +10 216 in 3-D: morphology of a TP-AGB star envelope

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, which combines the advantages of longer wavelength, high angular resolution and very high spectral resolution is the optimal investigative tool for this purpose. Mm waves pass through dust with almost no attenuation. Their spectrum is rich in molecular lines and hosts the fundamental lines of the ubiquitous CO molecule, allowing a tomographic reconstruction of the envelope structure. 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 the first detailed study of the CO(2-1) line emission in that envelope, made with the IRAM 30-m telescope. It revealed a series of dense gas shells, expanding at a uniform radial velocity. The limited resolution of the telescope (HPBW 11″) did not allow us to resolve the shell structure. We now report much higher angular resolution observations of CO(2-1), CO(1-0), CN(2-1) and C4H(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 much more intricate at high resolution than with an 11″ beam, its prevailing structure remains a pattern of thin, nearly concentric shells. The average separation between the brightest CO shells is 16″ in the outer envelope, where it appears remarkably constant. Closer to the star (< 40″), the shell pattern is denser and less regular, showing intermediary arcs. Outside the small (r < 0.3″) dust formation zone, the gas appears to expand radially at a constant velocity, 14.5 km s-1, 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 shell-intershell density contrast is found to be typically 3. 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.

Discovery of methyl silane and confirmation of silyl cyanide in IRC +10216

We report the discovery in space of methyl silane, CH3SiH3, from observations of ten rotational transitions between 80 and 350 GHz (Ju from 4 to 16) with the IRAM 30 m radio telescope. The molecule was observed in the envelope of the C-star IRC +10216. The observed profiles and our models for the expected emission of methyl silane suggest that the it is formed in the inner zones of the circumstellar envelope, 1-40 R*, with an abundance of (0.5-1) × 10-8 relative to H2. We also observed several rotational transitions of silyl cyanide (SiH3CN), confirming its presence in IRC +10216 in particular, and in space in general. Our models indicate that silyl cyanide is also formed in the inner regions of the envelope, around 20 R*, with an abundance relative to H2 of 6×10-10. The possible formation mechanisms of both species are discussed. We also searched for related chemical species but only upper limits could be obtained.

The growth of carbon chains in IRC +10216 mapped with ALMA

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.

Hints of a rotating spiral structure in the innermost regions around IRC +10216

The Atacama Large Millimeter/submillimeter Array (ALMA) is allowing us to study the innermost regions of the circumstellar envelopes of evolved stars with un-precedented precision and sensitivity. Key processes in the ejection of matter and dust from these objects occur in their inner zones. In this work, we present sub-arcsecond interferometric maps of transitions of metal-bearing molecules towards the prototypical C-rich evolved star IRC +10216. While Al-bearing molecules seem to be present as a roughly spherical shell, the molecular emission from the salts NaCl and KCl presents an elongation in the inner regions, with a central minimum. In order to accurately analyze the emission from the NaCl rotational lines, we present new calculations of the collisional rates for this molecule based on new spectroscopic constants. The most plausible interpretation for the spatial distribution of the salts is a spiral with a NaCl mass of 0.08M☉. Alternatively, a torus of gas and dust would result in similar structures as those observed. From the torus scenario we derive a mass of ~ 1.1 × 10-4M☉. In both cases, the spiral and the torus, the NaCl structure presents an inner minimum of 27 AU. In the case of the torus, the outer radius is 73 AU. The kinematics of both the spiral and the torus suggests that they are slowly expanding and rotating. Alternative explanations for the presence of the elongation are explored. The presence of these features only in KCl and NaCl might be a result of their comparatively high dipole moment with respect to the Al-bearing species.

THE PECULIAR DISTRIBUTION OF CH3CN IN IRC +10216 SEEN BY ALMA

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 [Formula: see text]. 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.

Probing non polar interstellar molecules through their protonated form: Detection of protonated cyanogen (NCCNH+)

Cyanogen (NCCN) is the simplest member of the series of dicyanopolyynes. It has been hypothesized that this family of molecules can be important constituents of interstellar and circumstellar media, although the lack of a permanent electric dipole moment prevents its detection through radioastronomical techniques. Here we present the first solid evidence of the presence of cyanogen in interstellar clouds through the detection of its protonated form toward the cold dark clouds TMC-1 and L483. Protonated cyanogen (NCCNH+) has been identified through the J = 5 - 4 and J = 10 - 9 rotational transitions using the 40m radiotelescope of Yebes and the IRAM 30m telescope. We derive beam averaged column densities for NCCNH+ of (8.6 ± 4.4) × 1010 cm-2 in TMC-1 and (3.9 ± 1.8) × 1010 cm-2 in L483, which translate to fairly low fractional abundances relative to H2, in the range (1-10) × 10-12. The chemistry of protonated molecules in dark clouds is discussed, and it is found that, in general terms, the abundance ratio between the protonated and non protonated forms of a molecule increases with increasing proton affinity. Our chemical model predicts an abundance ratio NCCNH+/NCCN of ~ 10-4, which implies that the abundance of cyanogen in dark clouds could be as high as (1-10) × 10-8 relative to H2, i.e., comparable to that of other abundant nitriles such as HCN, HNC, and HC3N.

Discovery of SiCSi in IRC +10216: A missing link between gas and dust carriers of Si-C bonds

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 CW Leo. 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 emission arises from a region of 6″ 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 Si-C bond in the dust formation zone and certainly both play a key role in the formation of SiC dust grains.

SI-BEARING MOLECULES TOWARD IRC+10216: ALMA UNVEILS THE MOLECULAR ENVELOPE OF CWLEO

We report the detection of SiS rotational lines in high-vibrational states as well as SiO and SiC2 lines in their ground vibrational state toward IRC+10216 during the Atacama Large Millimeter Array Cycle 0. The spatial distribution of these molecules shows compact emission for SiS and a more extended emission for SiO and SiC2, and also proves the existence of an increase in the SiC2 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 use 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 be 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.

Molecular shells in IRC+10216: tracing the mass loss history,

Thermally-pulsating AGB stars provide three-fourths of the matter returned to the interstellar medium. The mass and chemical composition of their ejecta largely control the chemical evolution of galaxies. Yet, both the mass loss process and the gas chemical composition remain poorly understood. We present maps of the extended ¹²CO and ¹³CO emissions in IRC+10216, the envelope of CW Leo, the high mass loss star the closest to the Sun. IRC+10216 is nearly spherical and expands radially with a velocity of 14.5 km s^-1. The observations were made On-the-Fly with the IRAM 30 m telescope; their sensibility, calibration, and angular resolution are far higher than all previous studies. The telescope resolution at λ = 1.3 mm (11″ HPBW) corresponds to an expansion time of 500 yr. The CO emission consists of a centrally peaked pedestal and a series of bright, nearly spherical shells. It peaks on CW Leo and remains relatively strong up to r_phot = 180″. Further out the emission becomes very weak and vanishes as CO gets photodissociated. As CO is the best tracer of the gas up to r_phot, the maps show the mass loss history in the last 8000 yr. The bright CO shells denote over-dense regions. They show that the mass loss process is highly variable on timescales of hundreds of years. The new data, however, do not support previous claims of a strong decrease of the average mass loss in the last few thousand years. The over-dense shells are not perfectly concentric and extend farther to the N-NW. The typical shell separation is 800-1000 yr in the middle of the envelope, but seems to increase outwards. The shell-intershell brightness contrast is ≥3. All those key features can be accounted for if CW Leo has a companion star with a period ≃800 yr that increases the mass loss rate when it comes close to periastron. Higher angular resolution observations are needed to fully resolve the dense shells and measure the density contrast. The latter plays an essential role in our understanding of the envelope chemistry.

Discovery of Time Variation of the Intensity of Molecular Lines in IRC+10216 in The Submillimeter and Far Infrared Domains

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 Herschel and with the IRAM 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.

Warm water vapour in the sooty outflow from a luminous carbon star

The detection of circumstellar water vapour around the ageing carbon star IRC +10216 challenged the current understanding of chemistry in old stars, because water was predicted to be almost absent in carbon-rich stars. Several explanations for the water were postulated, including the vaporization of icy bodies (comets or dwarf planets) in orbit around the star, grain surface reactions, and photochemistry in the outer circumstellar envelope. With a single water line detected so far from this one carbon-rich evolved star, it is difficult to discriminate between the different mechanisms proposed. Here we report the detection of dozens of water vapour lines in the far-infrared and sub-millimetre spectrum of IRC +10216 using the Herschel satellite. This includes some high-excitation lines with energies corresponding to approximately 1,000 K, which can be explained only if water is present in the warm inner sooty region of the envelope. A plausible explanation for the warm water appears to be the penetration of ultraviolet photons deep into a clumpy circumstellar envelope. This mechanism also triggers the formation of other molecules, such as ammonia, whose observed abundances are much higher than hitherto predicted.