Search

Context. The chemical evolution of asymptotic giant branch (AGB) stars is driven by repeated thermal pulses (TPs). The duration of a TP is only a few hundred years, whereas an inter-pulse period lasts 104 − 105 yr. Direct observations of TPs are hence unlikely. However, the detached shells seen in CO line emission that are formed as a result of a TP provide indirect constraints on the changes experienced by the star during the pulse. Aims. We aim to resolve the spatial and kinematic sub-structures in five detached-shell sources to provide detailed constraints for hydrodynamic models that describe the formation and evolution of the shells. Methods. We used observations of the 12CO (1 − 0) emission towards five carbon-AGB stars with ALMA (Atacama Large Millimeter/submillimeter Array), including previously published observations of the carbon AGB star U Ant. The data have angular resolutions of 0″​​.3 to 1″ and a velocity resolution of 0.3 km s−1. This enabled us to quantify spatial and kinematic structures in the shells. Combining the ALMA data with single-dish observations of the 12CO (1 − 0) to 12CO (4 − 3) emission towards the sources, we used radiative transfer models to compare the observed structures with previous estimates of the shell masses and temperatures. Results. The observed emission is separated into two distinct components: a more coherent, bright outer shell and a more filamentary, fainter inner shell. The kinematic information shows that the inner sub-shells move at a higher velocity relative to the outer sub-shells. The observed sub-structures reveal a negative velocity gradient outwards across the detached shells, confirming the predictions from hydrodynamical models. However, the models do not predict a double-shell structure, and the CO emission likely only traces the inner and outer edges of the shell, implying a lack of CO in the middle layers of the detached shell. Previous estimates of the masses and temperatures are consistent with originating mainly from the brighter subshell, but the total shell masses are likely lower limits. Also, additional structures in the form of partial shells outside the detached shell around V644 Sco, arcs within the shell of R Scl, and a partially filled shell for DR Ser indicate a more complicated evolution of the shells and mass-loss process throughout the TP cycle than previously assumed. Conclusions. The observed spatial and kinematical splittings of the shells appear consistent with results from the hydrodynamical models, provided the CO emission does not trace the H2 density distribution in the shell but rather traces the edges of the shells. The hydrodynamical models predict very different density profiles depending on the evolution of the shells and the different physical processes involved in the wind-wind interaction (e.g. heating and cooling processes). It is therefore not possible to constrain the total shell mass based on the CO observations alone. Additional features outside and inside the shells complicate the interpretation of the data. Complementary observations of, for example, CI as a dissociation product of CO would be necessary to understand the distribution of CO compared to H2, in addition to new detailed hydrodynamical models of the pre-pulse, pulse, and post-pulse wind. Only a comprehensive combination of observations and models will allow us to constrain the evolution of the shells and the changes in the star during the thermal-pulse cycle. [ABSTRACT FROM AUTHOR]

Context. The mass loss experienced on the asymptotic giant branch (AGB) at the end of the lives of low- and intermediate-mass stars is widely accepted to rely on radiation pressure acting on newly formed dust grains. Dust formation happens in the extended atmospheres of these stars, where the density, velocity, and temperature distributions are strongly affected by convection, stellar pulsation, and heating and cooling processes. The interaction between these processes and how that affects dust formation and growth is complex. Hence, characterising the extended atmospheres empirically is paramount to advance our understanding of the dust formation and wind-driving processes. Aims. We aim to determine the density, temperature, and velocity distributions of the gas in the extended atmosphere of the AGB star R Dor. Methods. We acquired observations using ALMA towards R Dor to study the gas through molecular line absorption and emission. We modelled the observed 12CO v = 0,  J = 2 − 1, v = 1,  J = 2 − 1, and 3 − 2 and 13CO v = 0,  J = 3 − 2 lines using the 3D radiative transfer code LIME to determine the density, temperature, and velocity distributions up to a distance of four times the radius of the star at sub-millimetre wavelengths. Results. The high angular resolution of the sub-millimetre maps allows for even the stellar photosphere to be spatially resolved. By analysing the absorption against the star, we infer that the innermost layer in the near-side hemisphere is mostly falling towards the star, while gas in the layer above that seems to be mostly outflowing. Interestingly, the high angular resolution of the ALMA Band 7 observations reveal that the velocity field of the gas seen against the star is not homogenous across the stellar disc. The gas temperature and density distributions have to be very steep close to the star to fit the observed emission and absorption, but they become shallower for radii larger than ∼1.6 times the stellar sub-millimetre radius. While the gas mass in the innermost region is hundreds of times larger than the mass lost on average by R Dor per pulsation cycle, the gas densities just above this region are consistent with those expected based on the mass-loss rate and expansion velocity of the large-scale outflow. Our fits to the line profiles require the velocity distribution on the far side of the envelope to be mirrored, on average, with respect to that on the near side. Using a sharp absorption feature seen in the CO v = 0,  J = 2 − 1 line, we constrained the standard deviation of the stochastic velocity distribution in the large-scale outflow to be ≲0.4 km s−1. We characterised two blobs detected in the CO v = 0,  J = 2 − 1 line and found densities substantially larger than those of the surrounding gas. The two blobs also display expansion velocities that are high relative to that of the large-scale outflow. Monitoring the evolution of these blobs will lead to a better understanding of the role of these structures in the mass-loss process of R Dor. [ABSTRACT FROM AUTHOR]

Context. Low- and intermediate-mass asymptotic giant stars and massive red supergiant stars are important contributors to the chemical enrichment of the Universe. They are among the most efficient dust factories of the Galaxy, harboring chemically rich circumstellar environments. Yet, the processes that lead to dust formation or the large-scale shaping of the mass loss still escape attempts at modeling. Aims. Through the ATOMIUM project, we aim to present a consistent view of a sample of 17 nearby cool evolved stars. Our goals are to unveil the dust-nucleation sites and morphologies of the circumstellar envelope of such stars and to probe ambient environments with various conditions. This will further enhance our understanding of the roles of stellar convection and pulsations, and that of companions in shaping the dusty circumstellar medium. Methods. Here we present and analyze VLT/SPHERE-ZIMPOL polarimetric maps obtained in the visible (645–820 nm) of 14 out of the 17 ATOMIUM sources. They were obtained contemporaneously with the ALMA high spatial resolution data. To help interpret the polarized signal, we produced synthetic maps of light scattering by dust, through 3D radiative transfer simulations with the RADMC3D code. Results. The degree of linear polarization (DoLP) observed by ZIMPOL spreads across several optical filters. We infer that it primarily probes dust located just outside of the point spread function of the central source, and in or near the plane of the sky. The polarized signal is mainly produced by structures with a total optical depth close to unity in the line of sight, and it represents only a fraction of the total circumstellar dust. The maximum DoLP ranges from 0.03–0.38 depending on the source, fractions that can be reproduced by our 3D pilot models for grains composed of olivine, melilite, corundum, enstatite, or forsterite. The spatial structure of the DoLP shows a diverse set of shapes, including clumps, arcs, and full envelopes. Only for three sources do we note a correlation between the ALMA CO υ = 0, J = 2−1 and SiO υ = 0, J = 5−4 lines, which trace the gas density, and the DoLP, which traces the dust. Conclusions. The clumpiness of the DoLP and the lack of a consistent correlation between the gas and the dust location show that, in the inner environment, dust formation occurs at very specific sites. This has potential consequences for the derived mass-loss rates and dust-to-gas ratio in the inner region of the circumstellar environment. Except for π1 Gru and perhaps GY Aql, we do not detect interactions between the circumstellar wind and the hypothesized companions that shape the wind at larger scales. This suggests that the orbits of any other companions are tilted out of the plane of the sky. [ABSTRACT FROM AUTHOR]

Context. The asymptotic giant branch (AGB) marks the final evolutionary stage of stars with initial masses between ~0.8 and 8 M®. During this phase, stars undergo copious mass loss, which contributes significantly to the enrichment of the interstellar medium. The well-accepted mass-loss mechanism requires radiation pressure acting on dust grains that form in the density-enhanced and extended AGB stellar atmospheres. The details of the mass-loss process are not yet well understood, however. For oxygen-rich AGB stars, which are the focus of this study, the dust grains that drive the wind are expected to scatter visible light very efficiently because their sizes are relative large. Aims. We study the distribution of dust in the inner wind of oxygen-rich AGB stars to advance our understanding of the wind-driving process. Methods. We observed light scattered off dust grains that form around three oxygen-rich AGB stars (W Hya, SWVir, and RCrt) with mass-loss rates between 10-7 and 10-6 M® yr-1 using the extreme-adaptive-optics imager and polarimeter SPHERE/ZIMPOL with three filters centred at 0.65, 0.75 and 0.82 μm. We compared the observed morphologies and the spectral dependence of the scattered light between the three sources and determined the radial profile, per image octant, of the dust density distribution around the closest target, WHya. Results. We find the distribution of dust to be asymmetric for the three targets. A biconical morphology is seen for R Crt, with a position angle that is very similar to those inferred from interferometric observations of maser emission and of mid-infrared continuum emission. The cause of the biconical outflow cannot be inferred from the ZIMPOL data, but we speculate that it might be the consequence of a circumstellar disc or of the action of strong magnetic fields. The dust grains polarise light more efficiently at 0.65 μm for R Crt and SWVir and at 0.82 μm for WHya. This indicates that at the time of the observations, the grains around SWVir and R Crt had sizes <0.1 μm, while those aroundWHya were larger, with sizes &0.1 μm. The asymmetric distribution of dust around R Crt makes the interpretation more uncertain for this star, however. We find that polarised light is produced already from within the visible photosphere of W Hya, which we reproduce using models with an inner dust shell that is optically thick to scattering. We fit radiative transfer models to the radial profile of the polarised light observed around WHya and find a steep dust density profile, with steepness varying considerably with direction. We find the wind-acceleration region of WHya to extend to at least ~7 R. This is in agreement with theoretical predictions of wind acceleration up to ~12 R, and highlights that ZIMPOL observations probe the crucial region around AGB stars where dust forms and is accelerated. [ABSTRACT FROM AUTHOR]

Context. On the asymptotic giant branch (AGB), Sun-like stars lose a large portion of their mass in an intensive wind and enrich the surrounding interstellar medium with nuclear processed stellar material in the form of molecular gas and dust. For a number of carbon-rich AGB stars, thin detached shells of gas and dust have been observed. These shells are formed during brief periods of increased mass loss and expansion velocity during a thermal pulse, and open up the possibility to study the mass-loss history of thermally pulsing AGB stars. Aims. We study the properties of dust grains in the detached shell around the carbon AGB star R Scl and aim to quantify the influence of the dust grain properties on the shape of the spectral energy distribution (SED) and the derived dust shell mass. Methods. We modelled the SED of the circumstellar dust emission and compared the models to observations, including new observations of Herschel/PACS and SPIRE (infrared) and APEX/LABOCA (sub-millimeter). We derived present-day mass-loss rates and detached shell masses for a variation of dust grain properties (opacities, chemical composition, grain size, and grain geometry) to quantify the influence of changing dust properties to the derived shell mass. Results. The best-fitting mass-loss parameters are a present-day dust mass-loss rate of 2 × 10−10M⊙ yr−1 and a detached shell dust mass of (2.9 ± 0.3) × 10−5M⊙. Compared to similar studies, the uncertainty on the dust mass is reduced by a factor of 4. We find that the size of the grains dominates the shape of the SED, while the estimated dust shell mass is most strongly affected by the geometry of the dust grains. Additionally, we find a significant sub-millimeter excess that cannot be reproduced by any of the models, but is most likely not of thermal origin. [ABSTRACT FROM AUTHOR]

Context. The stellar outflows of low- to intermediate-mass stars are characterised by a rich chemistry. Condensation of molecular gas species into dust grains is a key component in a chain of physical processes that leads to the onset of a stellar wind. In order to improve our understanding of the coupling between the micro-scale chemistry and macro-scale dynamics, we need to retrieve the abundance of molecules throughout the outflow. Aims. Our aim is to determine the radial abundance profile of SiO and HCN throughout the stellar outflow of R Dor, an oxygen-rich AGB star with a low mass-loss rate. SiO is thought to play an essential role in the dust-formation process of oxygen-rich AGB stars. The presence of HCN in an oxygen-rich environment is thought to be due to non-equilibrium chemistry in the inner wind. Methods. We analysed molecular transitions of CO, SiO, and HCN measured with the APEX telescope and all three instruments on the Herschel Space Observatory, together with data available in the literature. Photometric data and the infrared spectrum measured by ISO-SWS were used to constrain the dust component of the outflow. Using both continuum and line radiative transfer methods, a physical envelope model of both gas and dust was established. We performed an analysis of the SiO and HCN molecular transitions in order to calculate their abundances. Results. We have obtained an envelope model that describes the dust and the gas in the outflow, and determined the abundance of SiO and HCN throughout the region of the stellar outflow probed by our molecular data. For SiO, we find that the initial abundance lies between 5.5 × 10-5 and 6.0 × 10-5 with respect to H2. The abundance profile is constant up to 60 ± 10 R*, after which it declines following a Gaussian profile with an e-folding radius of 3.5 ± 0.5 × 1013 cm or 1.4 ± 0.2 R*. For HCN, we find an initial abundance of 5.0 × 10-7 with respect to H2. The Gaussian profile that describes the decline starts at the stellar surface and has an e-folding radius re of 1.85 ± 0.05 × 1015 cm or 74 ± 2 R*. Conclusions. We cannot unambiguously identify the mechanism by which SiO is destroyed at 60 ± 10 R*. The initial abundances found are higher than previously determined (except for one previous study on SiO), which might be due to the inclusion of higher-J transitions. The difference in abundance for SiO and HCN compared to high mass-loss rate Mira star IK Tau might be due to different pulsation characteristics of the central star and/or a difference in dust condensation physics. [ABSTRACT FROM AUTHOR]

Context. Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. Aims. We characterize the CO emission from the wind of the low mass-loss rate oxygen-richAGB star WHya using data obtained by the HIFI, PACS, and SPIRE instruments on board the Herschel Space Observatory and ground-based telescopes. 12CO and 13CO lines are used to constrain the intrinsic 12C/13C ratio from resolved HIFI lines. Methods. We combined a state-of-the-art molecular line emission code and a dust continuum radiative transfer code to model the CO lines and the thermal dust continuum. Results. The acceleration of the outflow up to about 5.5 kms-1 is quite slow and can be represented by a μ-type velocity law with index μ = 5. Beyond this point, acceleration up the terminal velocity of 7 kms-1 is faster. Using the J =10-9, 9-8, and 6-5 transitions, we find an intrinsic 12C/13C ratio of 18 ± 10 for WHya, where the error bar is mostly due to uncertainties in the 12CO abundance and the stellar flux around 4.6 µm. To match the low-excitation CO lines, these molecules need to be photo-dissociated at ~500 stellar radii. The radialdust emission intensity profile of our stellar wind model matches PACS images at 70 μm out to 20" (or 800 stellar radii). For larger radii the observed emission is substantially stronger than our model predicts, indicating that at these locations there is extra material present. Conclusions. The initial slow acceleration of the wind may imply inefficient dust formation or dust driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for WHya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M☉ or more. However, the uncertainty in the isotopologic ratio is large, which makes it difficult to set reliable limits on WHya's main-sequence mass. [ABSTRACT FROM AUTHOR]