The asymptotic-giant-branch star R Sculptoris is surrounded by a detached shell of dust and gas. The shell originates from a thermal pulse during which the star underwent a brief period of increased ...mass loss. It has hitherto been impossible to constrain observationally the timescales and mass-loss properties during and after a thermal pulse--parameters that determine the lifetime of the asymptotic giant branch and the amount of elements returned by the star. Here we report observations of CO emission from the circumstellar envelope and shell around R Sculptoris with an angular resolution of 1.3″. What was previously thought to be only a thin, spherical shell with a clumpy structure is revealed to also contain a spiral structure. Spiral structures associated with circumstellar envelopes have been previously seen, leading to the conclusion that the systems must be binaries. Combining the observational data with hydrodynamic simulations, we conclude that R Sculptoris is a binary system that underwent a thermal pulse about 1,800 years ago, lasting approximately 200 years. About 3 × 10(-3) solar masses of material were ejected at a velocity of 14.3 km s(-1) and at a rate around 30 times higher than the pre-pulse mass-loss rate. This shows that about three times more mass was returned to the interstellar medium during and immediately after the pulse than previously thought.
Context. Asymptotic giant branch (AGB) stars experience strong mass loss driven by dust particles formed in the upper atmospheres. The dust is released into the interstellar medium, and replenishes ...galaxies with synthesised material from the star. The dust grains further act as seeds for continued dust growth in the diffuse medium of galaxies. As such, understanding the properties of dust produced during the asymptotic giant branch phase of stellar evolution is important for understanding the evolution of stars and galaxies. Recent observations of the carbon AGB star R Scl have shown that observations at far-infrared and submillimetre wavelengths can effectively constrain the grain sizes in the shell, while the total mass depends on the structure of the grains (solid vs. hollow or fluffy). Aims. We aim to constrain the properties of the dust observed in the submillimetre in the detached shells around the three carbon AGB stars U Ant, DR Ser, and V644 Sco, and to investigate the constraints on the dust masses and grain sizes provided by far-infrared and submm observations. Methods. We observed the carbon AGB stars U Ant, DR Ser, and V644 Sco at 870 μm using LABOCA on APEX. Combined with observations from the optical to far-infrared, we produced dust radiative transfer models of the spectral energy distributions (SEDs) with contributions from the stars, present-day mass-loss and detached shells. We assume spherical, solid dust grains, and test the effect of different total dust masses and grain sizes on the SED, and attempted to consistently reproduce the SEDs from the optical to the submm. Results. We derive dust masses in the shells of a few 10−5 M ⊙. The best-fit grain radii are comparatively large, and indicate the presence of grains between 0.1 μm and 2 μm. The LABOCA observations suffer from contamination from 12CO (3 − 2), and hence gives fluxes that are higher than the predicted dust emission at submm wavelengths. We investigate the effect on the best-fitting models by assuming different degrees of contamination and show that far-infrared and submillimetre observations are important to constrain the dust mass and grain sizes in the shells. Conclusions. Spatially resolved observations of the detached shells in the far-infrared and submillimetre effectively constrain the temperatures in the shells, and hence the grain sizes. The dust mass is also constrained by the observations, but additional observations are needed to constrain the structure of the grains.
Context. On the asymptotic giant branch (AGB) low- and intermediate-mass stars eject a large fraction of their envelope, but the mechanism driving these outflows is still poorly understood. For ...oxygen-rich AGB stars, the wind is thought to be driven by radiation pressure caused by scattering of radiation off dust grains. Aims. We study the photosphere, the warm molecular layer, and the inner wind of the close-by oxygen-rich AGB star R Doradus. We focus on investigating the spatial distribution of the dust grains that scatter light and whether these grains can be responsible for driving the outflow of this star. Methods. We use high-angular-resolution images obtained with SPHERE/ZIMPOL to study R Dor and its inner envelope in a novel way. We present observations in filters V, cntHα, and cnt820 and investigate the surface brightness distribution of the star and of the polarised light produced in the inner envelope. Thanks to second-epoch observations in cntHα, we are able to see variability on the stellar photosphere. We study the polarised-light data using a continuum-radiative-transfer code that accounts for direction-dependent scattering of photons off dust grains. Results. We find that in the first epoch the surface brightness of R Dor is asymmetric in V and cntHα, the filters where molecular opacity is stronger, while in cnt820 the surface brightness is closer to being axisymmetric. The second-epoch observations in cntHα show that the morphology of R Dor has changed completely in a timespan of 48 days to a more axisymmetric and compact configuration. This variable morphology is probably linked to changes in the opacity provided by TiO molecules in the extended atmosphere. The observations show polarised light coming from a region around the central star. The inner radius of the region from where polarised light is seen varies only by a small amount with azimuth. The value of the polarised intensity, however, varies by between a factor of 2.3 and 3.7 with azimuth for the different images. We fit the radial profile of the polarised intensity using a spherically symmetric model and a parametric description of the dust density profile, ρ(r) = ρ°r− n. On average, we find exponents of − 4.5 ± 0.5 that correspond to a much steeper density profile than that of a wind expanding at constant velocity. The dust densities we derive imply an upper limit for the dust-to-gas ratio of ~2 × 10-4 at 5.0 R⋆. Considering all the uncertainties in observations and models, this value is consistent with the minimum values required by wind-driving models for the onset of a wind, of ~3.3 × 10-4. However, if the steep density profile we find extends to larger distances from the star, the dust-to-gas ratio will quickly become too small for the wind of R Dor to be driven by the grains that produce the scattered light.
Context. The outflows of oxygen-rich asymptotic giant branch (AGB) stars are thought to be driven by radiation pressure due to the scattering of photons on relatively large grains, with sizes of ...tenths of microns. The details of the formation of dust in the extended atmospheres of these stars and, therefore, the mass-loss process, is still not well understood. Aims. We aim to constrain the distribution of the gas and the composition and properties of the dust grains that form in the inner circumstellar environment of the archetypal Mira variable o Cet. Methods. We obtained quasi-simultaneous observations using ALMA and SPHERE/ZIMPOL on the Very Large Telescope (VLT) to probe the distribution of gas and large dust grains, respectively. Results. The polarized light images show dust grains around Mira A, but also around the companion, Mira B, and a dust trail that connects the two sources. The ALMA observations show that dust around Mira A is contained in a high-gas-density region with a significant fraction of the grains that produce the polarized light located at the edge of this region. Hydrodynamical and wind-driving models show that dust grains form efficiently behind shock fronts caused by stellar pulsation or convective motions. The distance at which we observe the density decline (a few tens of au) is, however, significantly larger than expected for stellar-pulsation-induced shocks. Other possibilities for creating the high-gas-density region are a recent change in the mass-loss rate of Mira A or interactions with Mira B. We are not able to determine which of these scenarios is correct. We constrained the gas density, temperature, and velocity within a few stellar radii from the star by modelling the CO v = 1, J = 3−2 line. We find a mass (~3.8 ± 1.3) × 10−4 M⊙ to be contained between the stellar millimetre photosphere, R⋆338 GHz $R^{\textrm{338~GHz}}_{\star}$ R⋆338 GHz , and 4 R⋆338 GHz $4~R^{\textrm{338~GHz}}_{\star}$ 4 R⋆338 GHz . Our best-fit models with lower masses also reproduce the 13CO v = 0, J = 3−2 line emission from this region well. We find TiO2 and AlO abundances corresponding to 4.5% and <0.1% of the total titanium and aluminium expected for a gas with solar composition. The low abundance of AlO allows for a scenario in which Al depletion into dust happens already very close to the star, as expected from thermal dust emission observations and theoretical calculations of Mira variables. The relatively large abundance of aluminium for a gas with solar composition allows us to constrain the presence of aluminium oxide grains based on the scattered light observations and on the gas densities we obtain. These models imply that aluminium oxide grains could account for a significant fraction of the total aluminium atoms in this region only if the grains have sizes ≲0.02 μm. This is an order of magnitude smaller than the maximum sizes predicted by dust-formation and wind-driving models. Conclusions. The study we present highlights the importance of coordinated observations using different instruments to advance our understanding of dust nucleation, dust growth, and wind driving in AGB stars.
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 10 4 − 10 5 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 12 CO (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 12 CO (1 − 0) to 12 CO (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 H 2 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 H 2 , 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.
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.
Recent observations at subarcsecond resolution, now possible also at submillimeter wavelengths, have shown intricate circumstellar structures around asymptotic giant branch (AGB) stars, mostly ...attributed to binary interaction. The results presented here are part of a larger project aimed at investigating the effects of a binary companion on the morphology of circumstellar envelopes (CSEs) of AGB stars.
AGB stars are characterized by intense stellar winds that build CSEs around the stars. Here, the CO(
= 3→2) emission from the CSE of the binary S-type AGB star W Aql has been observed at subarcsecond resolution using ALMA. The aim of this paper is to investigate the wind properties of the AGB star and to analyse how the known companion has shaped the CSE.
The average mass-loss rate during the creation of the detected CSE is estimated through modelling, using the ALMA brightness distribution and previously published single-dish measurements as observational constraints. The ALMA observations are presented and compared to the results from a 3D smoothed particle hydrodynamics (SPH) binary interaction model with the same properties as the W Aql system and with two different orbital eccentricities. Three-dimensional radiative transfer modelling is performed and the response of the interferometer is modelled and discussed.
The estimated average mass-loss rate of W Aql is
= 3.0×10
M
yr
and agrees with previous results based on single-dish CO line emission observations. The size of the emitting region is consistent with photodissociation models. The inner 10″ of the CSE is asymmetric with arc-like structures at separations of 2-3″ scattered across the denser sections. Further out, weaker spiral structures at greater separations are found, but this is at the limit of the sensitivity and field of view of the ALMA observations.
The CO(
= 3→2) emission is dominated by a smooth component overlayed with two weak arc patterns with different separations. The larger pattern is predicted by the binary interaction model with separations of ~10″ and therefore likely due to the known companion. It is consistent with a binary orbit with low eccentricity. The smaller separation pattern is asymmetric and coincides with the dust distribution, but the separation timescale (200 yrs) is not consistent with any known process of the system. The separation of the known companions of the system is large enough to not have a very strong effect on the circumstellar morphology. The density contrast across the envelope of a binary with an even larger separation will not be easily detectable, even with ALMA, unless the orbit is strongly asymmetric or the AGB star has a much larger mass-loss rate.
High-resolution observations of the extended atmospheres of asymptotic giant branch (AGB) stars can now directly be compared to the theories that describe stellar mass loss. Using Atacama Large ...Millimeter/submillimeter Array (ALMA) high angular resolution (30 × 42 mas) observations, we have for the first time resolved stellar rotation of an AGB star, R Dor. We measure an angular rotation velocity of ωR sin i = (3.5 ± 0.3) × 10−9 rad s−1, which indicates a rotational velocity of |υrot sin i| = 1.0 ± 0.1 km s−1 at the stellar surface (R* = 31.2 mas at 214 GHz). The rotation axis projected on the plane of the sky has a position angle Φ = 7 ± 6°. We find that the rotation of R Dor is two orders of magnitude faster than expected for a solitary AGB star that will have lost most of its angular momentum. Its rotational velocity is consistent with angular momentum transfer from a close companion. As a companion has not been directly detected, we suggest R Dor has a low-mass, close-in companion. The rotational velocity approaches the critical velocity, set by the local sound speed in the extended envelope, and is thus expected to affect the mass-loss characteristics of R Dor.
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 12 CO v = 0, J = 2 − 1, v = 1, J = 2 − 1, and 3 − 2 and 13 CO 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.
Aims. We study the circumstellar evolution of the binary HD 101584, consisting of a post-AGB star and a low-mass companion, which is most likely a post-common-envelope-evolution system. Methods. We ...used ALMA observations of the 12CO, 13CO, and C18O J = 2−1 lines and the 1.3 mm continuum to determine the morphology, kinematics, masses, and energetics of the circumstellar environment. Results. The circumstellar medium has a bipolar hour-glass structure, seen almost pole-on, formed by an energetic jet, ≈150 km s-1. We conjecture that the circumstellar morphology is related to an event that took place ≈500 yr ago, possibly a capture event where the companion spiraled in towards the AGB star. However, the kinetic energy of the accelerated gas exceeds the released orbital energy, and, taking into account the expected energy transfer efficiency of the process, the observed phenomenon does not match current common-envelope scenarios. This suggests that another process must augment, or even dominate, the ejection process. A significant amount of material resides in an unresolved region, presumably in the equatorial plane of the binary system.