ABSTRACT We present and analyse SO and SO2, recently observed with high angular resolution and sensitivity in a spectral line survey with ALMA, for two oxygen-rich AGB stars: the low mass-loss rate R ...Dor and high mass-loss rate IK Tau. We analyse 8 lines of SO detected towards both stars, 78 lines of SO2 detected towards R Dor, and 52 lines of SO2 detected towards IK Tau. We detect several lines of 34SO, 33SO, and 34SO2 towards both stars, and tentatively S18O towards R Dor, and hence derive isotopic ratios for these species. The spatially resolved observations show us that the two sulphur oxides are co-located towards R Dor and trace out the same wind structures in the circumstellar envelope. Much of the emission is well reproduced with a Gaussian abundance distribution spatially centred on the star. Emission from the higher energy levels of SO and SO2 towards R Dor provides evidence in support of a rotating inner region of gas identified in earlier work. The new observations allow us to refine the abundance distribution of SO in IK Tau derived from prior observations with single antennas, and confirm that the distribution is shell like with the peak in the fractional abundance not centred on the star. The confirmation of different types of SO abundance distributions will help fine-tune chemical models and allows for an additional method to discriminate between low and high mass-loss rates for oxygen-rich AGB stars.
Context. At the end of their lives, most stars lose a significant amount of mass through a stellar wind. The specific physical and chemical circumstances that lead to the onset of the stellar wind ...for cool luminous stars are not yet understood. Complex geometrical morphologies in the circumstellar envelopes prove that various dynamical and chemical processes are interlocked and that their relative contributions are not easy to disentangle. Aims. We aim to study the inner-wind structure (R< 250 R⋆) of the well-known red supergiant VY CMa, the archetype for the class of luminous red supergiant stars experiencing high mass loss. Specifically, the objective is to unravel the density structure in the inner envelope and to examine the chemical interaction between gas and dust species. Methods. We analyse high spatial resolution (~0.̋24×0.̋13) ALMA science verification (SV) data in band 7, in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. Results. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50° measured from north to east. However, this picture cannot capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process that prevents all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form while, for lower values of this ratio, thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons does not seem to be a viable explanation for the detection of gaseous NaCl at 220 R⋆ from the central star, so instead, we propose shock-induced sputtering driven by localized mass ejection events as an alternative. Conclusions. The analysis of the NaCl lines demonstrates the capabilities of ALMA to decode the geometric morphologies and chemical pathways prevailing in the winds of evolved stars. These early ALMA results prove that the envelopes surrounding evolved stars are far from homogeneous, and that a variety of dynamical and chemical processes dictate the wind structure.
Context. In the inner regions of asymptotic giant branch (AGB) outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium chemical ...models. The presence of the majority of these species can be explained by shock-induced non-equilibrium chemical models, where shocks caused by the pulsating star take the chemistry out of equilibrium in the inner region. Moreover, a non-uniform density structure has been detected in several AGB outflows. Both large-scale structures, such as spirals and disks, and small-scale density inhomogeneities or clumps have been observed. These structures may also have a considerable impact on the circumstellar chemistry. A detailed parameter study on the quantitative effects of a non-homogeneous outflow has so far not been performed. Aims. We examine the effects of a non-uniform density distribution within an AGB outflow on its chemistry by considering a stochastic, clumpy density structure. Methods. We implement a porosity formalism for treating the increased leakage of light associated with radiation transport through a clumpy, porous medium. We then use this method to examine the effects from the altered UV radiation field penetration on the chemistry, accounting also for the increased reaction rates of two-body processes in the overdense clumps. The specific clumpiness is determined by three parameters: the characteristic length scale of the clumps at the stellar surface, the clump volume filling factor, and the inter-clump density contrast. In this paper, the clumps are assumed to have a spatially constant volume filling factor, which implies that they expand as they move outward in the wind. Results. We present a parameter study of the effect of clumping and porosity on the chemistry throughout the outflow. Both the higher density within the clumps and the increased UV radiation field penetration have an important impact on the chemistry, as they both alter the chemical pathways throughout the outflow. The increased amount of UV radiation in the inner region leads to photodissociation of parent species, releasing the otherwise deficient elements. We find an increased abundance in the inner region of all species not expected to be present assuming thermodynamic equilibrium chemistry, such as HCN in O-rich outflows, H2O in C-rich outflows, and NH3 in both. Conclusions. A non-uniform density distribution directly influences the chemistry throughout the AGB outflow, both through the density structure itself and through its effect on the UV radiation field. Species not expected to be present in the inner region of the outflow assuming thermodynamic equilibrium chemistry are now formed in this region, including species that are not formed in greater abundance by shock-induced non-equilibrium chemistry models. Outflows whose clumps have a large overdensity and that are very porous to the interstellar UV radiation field yield abundances comparable to those observed in O-rich and C-rich outflows for most of the unexpected species investigated. The inner wind abundances of H2O in C-rich outflows and of NH3 in O-rich and C-rich outflows are however underpredicted.
Sulphur-bearing molecules in AGB stars Danilovich, T.; Ramstedt, S.; Gobrecht, D. ...
Astronomy and astrophysics (Berlin),
10/2018, Letnik:
617
Journal Article
Recenzirano
Odprti dostop
Context. Sulphur has long been known to form different molecules depending on the chemical composition of its environment. More recently, the sulphur-bearing molecules SO and H2S have been shown to ...behave differently in oxygen-rich asymptotic giant branch (AGB) circumstellar envelopes of different densities. Aims. By surveying a diverse sample of AGB stars for CS and SiS emission, we aim to determine in which environments these sulphur-bearing molecules most readily occur. We include sources with a range of mass-loss rates and carbon-rich, oxygen-rich, and mixed S-type chemistries. Where these molecules are detected, we aim to determine their CS and SiS abundances. Methods. We surveyed 20 AGB stars of different chemical types using the APEX telescope, and combined this with an IRAM 30 m and APEX survey of CS and SiS emission towards over 30 S-type stars. For those stars with detections, we performed radiative transfer modelling to determine abundances and abundance distributions. Results. We detect CS towards all the surveyed carbon stars, some S-type stars, and the highest mass-loss rate oxygen-rich stars, (Ṁ ≥ 5 × 10−6 M⊙ yr−1). SiS is detected towards the highest mass-loss rate sources of all chemical types (Ṁ ≥ 8 × 10−7 M⊙ yr−1). We find CS peak fractional abundances ranging from ~4 × 10−7 to ~2 × 10−5 for the carbon stars, from ~3 × 10−8 to ~1 × 10−7 for the oxygen-rich stars, and from ~1 × 10−7 to ~8 × 10−6 for the S-type stars. We find SiS peak fractional abundances ranging from ~9 × 10−6 to ~2 × 10−5 for the carbon stars, from ~5 × 10−7 to ~2 × 10−6 for the oxygen-rich stars, and from ~2 × 10−7 to ~2 × 10−6 for the S-type stars. Conclusions. Overall, we find that wind density plays an important role in determining the chemical composition of AGB circumstellar envelopes. It is seen that for oxygen-rich AGB stars both CS and SiS are detected only in the highest density circumstellar envelopes and their abundances are generally lower than for carbon-rich AGB stars by around an order of magnitude. For carbon-rich and S-type stars SiS was also only detected in the highest density circumstellar envelopes, while CS was detected consistently in all surveyed carbon stars and sporadically among the S-type stars.
The mechanisms responsible for heating the extended atmospheres of early-M spectral-type supergiants are poorly understood. So too is the subsequent role these mechanisms play in driving the large ...mass-loss rates of these stars. Here we present ALMA long (i.e., ~16 km) baseline 338 GHz (0.89 mm) continuum observations of the free-free emission in the extended atmosphere of the M2 spectral-type supergiant Betelgeuse. The spatial resolution of 14 mas exquisitely resolves the atmosphere, revealing it to have a mean temperature of 2760 K at ~1.3 R⋆, which is below both the photospheric effective temperature (Teff = 3690 K) and the temperatures at ~2 R⋆. This is unambiguous proof for the existence of an inversion of the mean temperature in the atmosphere of a red supergiant. The emission is clearly not spherically symmetric with two notable deviations from a uniform disk detected in both the images and visibilities. The most prominent asymmetry is located in the north-east quadrant of the disk and is spatially resolved showing it to be highly elongated with an axis-ratio of 2.4 and occupying ~5% of the disk projected area. Its temperature is approximately 1000 K above the measured mean temperature at 1.3 R⋆. The other main asymmetry is located on the disk limb almost due east of the disk center and occupies ~3% of the disk projected area. Both emission asymmetries are clear evidence for localized heating taking place in the atmosphere of Betelgeuse. We suggest that the detected localized heating is related to magnetic activity generated by large-scale photospheric convection.
Abstract
Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss, which spreads heavy elements into the interstellar medium. The structure and ...dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO J = 2 − 1 line from the red supergiant star $\mu$ Cep with the NOEMA interferometer. In the line the synthesized beam was 0.92 × 0.72 arcsec (590 × 462 au at 641 pc). The continuum map shows only the unresolved contribution of the free–free emission of the star chromosphere. The continuum-subtracted channel maps reveal a very inhomogeneous and clumpy circumstellar environment. In particular, we detected a bright CO clump, as bright as the central source in the line, at 1.80 arcsec south-west from the star, in the blue channel maps. After a deprojection of the radial velocity assuming two different constant wind velocities, the observations were modelled using the 3D radiative transfer code lime to derive the characteristics of the different structures. We determine that the gaseous clumps observed around $\mu$ Cep are responsible for a mass-loss rate of $(4.9 \pm 1.0) \times 10^{-7}\, {\rm M}_\odot \, {\rm yr}^{-1}$, in addition to a spatially unresolved wind component with an estimated mass-loss rate of $2.0 \times 10^{-6}\, {\rm M}_\odot \, {\rm yr}^{-1}$. Therefore, the clumps have a significant role in $\mu$ Cep’s mass loss ($\ge 25 {{\ \rm per\ cent}}$). We cannot exclude that the unresolved central outflow may be made of smaller unresolved clumps.
Context. Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important ...in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims. This study attempts to spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud (LMC). Methods. The Atacama Large Millimeter Array was used to observe two OH/IR stars and four carbon stars in the LMC in the CO J = 2−1 line. Results. We present the first measurement of expansion velocities in extragalactic carbon stars. All four C stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. The derived gas-to-dust ratios allow the predicted velocities to agree with the observed gas-to-dust ratios. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with extreme C stars with properties that are more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than those derived from dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions. Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances (hence luminosities), it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear to be similar. This is in agreement with dynamical dust-driven wind models.
The processes leading to dust formation and the subsequent role it plays in driving mass loss in cool evolved stars is an area of intense study. Here we present high resolution ALMA Science ...Verification data of the continuum emission around the highly evolved oxygen-rich red supergiant VY CMa. These data enable us to study the dust in its inner circumstellar environment at a spatial resolution of 129 mas at 321 GHz and 59 mas at 658 GHz, thus allowing us to trace dust on spatial scales down to 11 R⋆ (71 AU). Two prominent dust components are detected and resolved. The brightest dust component, C, is located 334 mas (61 R⋆) southeast of the star and has a dust mass of at least 2.5 × 10-4 M⊙. It has a dust emissivity spectral index of β = −0.1 at its peak, implying that it is optically thick at these frequencies with a cool core of Td ≲ 100 K. Interestingly, not a single molecule in the ALMA data has emission close to the peak of this massive dust clump. The other main dust component, VY, is located at the position of the star and contains a total dust mass of 4.0 × 10-5 M⊙. It also contains a weaker dust feature extending over 60 R⋆ to the north with the total component having a typical dust emissivity spectral index of β = 0.7. We find that at least 17% of the dust mass around VY CMa is located in clumps ejected within a more quiescent roughly spherical stellar wind, with a quiescent dust mass loss rate of 5 × 10-6 M⊙yr-1. The anisotropic morphology of the dust indicates a continuous, directed mass loss over a few decades, suggesting that this mass loss cannot be driven by large convection cells alone.
Aims. In this study we intend to examine rotational emission lines of two isotopologues of water: H217O and H218O. By determining the abundances of these molecules, we aim to use the derived ...isotopologue – and hence oxygen isotope – ratios to put constraints on the masses of a sample of M-type AGB stars that have not been classified as OH/IR stars. Methods. We have used detailed radiative transfer analysis based on the accelerated lambda iteration method to model the circumstellar molecular line emission of H217O and H218O for IK Tau, R Dor, W Hya, and R Cas. The emission lines used to constrain our models came from Herschel/HIFI and Herschel/PACS observations and are all optically thick, meaning that full radiative transfer analysis is the only viable method of estimating molecular abundance ratios. Results. We find generally low values of the 17O/18O ratio for our sample, ranging from 0.15 to 0.69. This correlates with relatively low initial masses, in the range ~ 1.0 to 1.5 M⊙ for each source, based on stellar evolutionary models. We also find ortho-to-para ratios close to 3, which are expected from warm formation predictions. Conclusions. The 17O/18O ratios found for this sample are at the lower end of the range predicted by stellar evolutionary models, indicating that the sample chosen had relatively low initial masses.