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

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.

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.

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.

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.

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.

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 millimeter IRAM-30 m line survey toward IK Tau

AIMS

We aim to investigate the physical and chemical properties of the molecular envelope of the oxygen-rich AGB star IK Tau.

METHODS

We carried out a millimeter wavelength line survey between ~79 and 356 GHz with the IRAM-30 m telescope. We analysed the molecular lines detected in IK Tau using the population diagram technique to derive rotational temperatures and column densities. We conducted a radiative transfer analysis of the SO2 lines, which also helped us to verify the validity of the approximated method of the population diagram for the rest of the molecules.

RESULTS

For the first time in this source we detected rotational lines in the ground vibrational state of HCO+, NS, NO, and H2CO, as well as several isotopologues of molecules previously identified, namely, C18O, Si17O, Si18O, 29SiS, 30SiS, Si34S, H13CN, 13CS, C34S, H234S, 34SO, and 34SO2. We also detected several rotational lines in vibrationally excited states of SiS and SiO isotopologues, as well as rotational lines of H2O in the vibrationally excited state ν2=2. We have also increased the number of rotational lines detected of molecules that were previously identified toward IK Tau, including vibrationally excited states, enabling a detailed study of the molecular abundances and excitation temperatures. In particular, we highlight the detection of NS and H2CO with fractional abundances of f(NS)~10-8 and f(H2CO)~[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 SO2 a warm component with Trot~290 K is also detected.

CONCLUSIONS

With a total of ~350 lines detected of 34 different molecular species (including different isotopologues), IK Tau displays a rich chemistry for an oxygen-rich circumstellar envelope. The detection of carbon bearing molecules like H2CO, as well as the discrepancies found between our derived abundances and the predictions from chemical models for some molecules, highlight the need for a revision of standard chemical models. We were able to identify at least two different emission components in terms of rotational temperatures. The warm component, which is mainly traced out by SO2, is probably arising from the inner regions of the envelope (at ≲8R∗) where SO2 has a fractional abundance of f(SO2)~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.

High-Resolution Rotational Spectrum, Dunham Coefficients, and Potential Energy Function of NaCl

We report laboratory spectroscopy for the first time of the J = 1-0 and J = 2-1 lines of Na35Cl and Na37Cl in several vibrational states. The hyperfine structure has been resolved in both transitions for all vibrational levels, which permit us to predict with high accuracy the hyperfine splitting of the rotational transitions of the two isotopologues at higher frequencies. The new data have been merged with all previous works at microwave, millimeter, and infrared wavelengths and fitted to a series of mass-independent Dunham parameters and to a potential energy function. The obtained parameters have been used to compute a new dipole moment function, from which the dipole moment for infrared transitions up to Δv = 8 has been derived. Frequency and intensity predictions are provided for all rovibrational transitions up to J = 150 and v = 8, from which the ALMA data of evolved stars can be modeled and interpreted.