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.

History of two mass loss processes in VY CMa

Red Supergiant stars (RSGs) are known to eject large amounts of material during this evolutionary phase. However, the processes powering the mass ejection in low- and intermediate-mass stars do not work for RSGs and the mechanism that drives the ejection remains unknown. Different mechanisms have been proposed as responsible for this mass ejection including Alfvén waves, large convective cells, and magnetohydrodynamical (MHD) disturbances at the photosphere, but so far little is known about the actual processes taking place in these objects. Here we present high angular resolution interferometric ALMA maps of VY CMa continuum and molecular emission, which resolve the structure of the ejecta with unprecedented detail. We reconstructed the 3D structure of the gas traced by the different species. It allowed us to study the morphology and kinematics of the gas traced by the different species surrounding VY CMa. Two types of ejecta are clearly observed: extended, irregular, and vast ejecta surrounding the star that are carved by localized fast outflows. The structure of the outflows is found to be particularly flat. We present a 3D reconstruction of these outflows and proof of the carving. This indicates that two different mass loss processes take place in this massive star. We tentatively propose the physical cause for the formation of both types of structures. These results provide essential information on the mass loss processes of RSGs and thus of their further evolution.

Detection of infrared fluorescence of carbon dioxide in R Leonis with SOFIA/EXES

We report on the detection of hot CO2 in the O-rich AGB star R Leo based on high spectral resolution observations in the range 12.8 - 14.3 μm carried out with the Echelon-cross-Echelle Spectrograph (EXES) mounted on the Stratospheric Observatory for Infrared Astronomy (SOFIA). We have found ≃ 240 CO2 emission lines in several vibrational bands. These detections were possible thanks to a favorable Doppler shift that allowed us to avoid contamination by telluric CO2 features. The highest excitation lines involve levels at an energy of ≃ 7000 K. The detected lines are narrow (average deconvolved width ≃ 2.5 km s-1) and weak (usually ≲ 10% the continuum). A ro-vibrational diagram shows that there are three different populations, warm, hot, and very hot, with rotational temperatures of ≃ 550, 1150, and 1600 K, respectively. From this diagram, we derive a lower limit for the column density of ≃ 2.2 × 1016 cm-2. Further calculations based on a model of the R Leo envelope suggest that the total column density can be as large as 7 × 1017 cm -2 and the abundance with respect to H2 - 2.5 × 10-5. The detected lines are probably formed due to de-excitation of CO2 molecules from high energy vibrational states, which are essentially populated by the strong R Leo continuum at 2.7 and 4.2 μm.

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.

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.

Gas infall and possible circumstellar rotation in R Leo

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

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.

The Maser Emitting Structure and Time Variability of the SIS Lines J = 14 - 13 and 15 - 14 in IRC+10216

We present new high angular resolution interferometer observations of the υ = 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 30 m 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.

The Abundance of SiC2 in Carbon Star Envelopes: Evidence that SiC2 is a gas-phase precursor of SiC dust

CONTEXT

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. To date, the ring molecule SiC2 has been observed in a handful of evolved stars, while SiC and Si2C have only been detected in the C-star envelope IRC +10216.

AIMS

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.

METHODS

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.

RESULTS

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.

CONCLUSIONS

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.

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.

Carbon Chemistry in IRC+10216: Infrared Detection of Diacetylene

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 ν 6 + ν 8 ( σ u + ) . The analysis of these lines through a ro-vibrational diagram indicates that the column density of C4H2 is (2.4 ± 1.5) × 1016 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 ≃ 20 R ⋆ 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.

Clues to NaCN formation

CONTEXT

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.

AIMS

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.

METHODS

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.

RESULTS

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 CH3CN, current photochemical models fail to reproduce this CN reservoir. We also found that the abundance peak of NaCN appears at a radius of 3 × 1015cm, 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.

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 Abundance of C2H4 in the Circumstellar Envelope of IRC+10216

High spectral resolution mid-IR observations of ethylene (C2H4) towards the AGB star IRC+10216 were obtained using the Texas Echelon Cross Echelle Spectrograph (TEXES) at the NASA Infrared Telescope Facility (IRTF). Eighty ro-vibrational lines from the 10.5 µm vibrational mode ν7 with J ≲ 30 were detected in absorption. The observed lines are divided into two groups with rotational temperatures of 105 and 400 K (warm and hot lines). The warm lines peak at ≃ -14 km s-1 with respect to the systemic velocity, suggesting that they are mostly formed outwards from ≃ 20R⋆. The hot lines are centered at -10 km s-1 indicating that they come from a shell between 10 and 20R⋆. 35% of the observed lines are unblended and can be fitted with a code developed to model the emission of a spherically symmetric circumstellar envelope. The analysis of several scenarios reveal that the C2H4 abundance relative to H2 in the range 5 - 20R⋆ is 6.9 × 10-8 in average and it could be as high as 1.1 × 10-7. Beyond 20R⋆, it is 8.2 × 10-8. The total column density is (6.5 ± 3.0) × 1015 cm-2. C2H4 is found to be rotationally under local thermodynamical equilibrium (LTE) and vibrationally out of LTE. One of the scenarios that best reproduce the observations suggests that up to 25% of the C2H4 molecules at 20R⋆ could condense onto dust grains. This possible depletion would not influence significantly the gas acceleration although it could play a role in the surface chemistry on the dust grains.

The abundance of 28Si32S, 29Si32S, 28Si34S, and 30Si32S in the inner layers of the envelope of IRC+10216

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