Free Access
Issue
A&A
Volume 570, October 2014
Article Number A101
Number of page(s) 9
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201424892
Published online 30 October 2014

© ESO, 2014

1. Introduction

During the late stages of stellar evolution, asymptotic giant branch (AGB) stars build up large circumstellar envelopes (CSEs) of dust and gas. The CSEs are a direct consequence of the mass loss the stars experience on the AGB, carrying material processed inside the star to the interstellar medium (ISM). CSEs are generally assumed to be formed through an isotropic mass loss from the stellar surface, building up a spherically symmetric envelope (e.g., Castro-Carrizo et al. 2010). Through their relatively simple geometry, the CSEs around AGB stars are excellent astrochemical and astrophysical laboratories. However, variations in the mass-loss rate over time, and any directional dependence of the mass loss will influence the shape of the CSE.

Recent images of thermal dust emission in the far-infrared taken with the PACS instrument aboard the Herschel Space Telescope show a variety of envelope shapes on very large scales (Cox et al. 2012), including arc-like structures (Decin et al. 2012). Spiral-shaped CSEs have been found around a handful of sources, such as AFGL 3068 (Mauron & Huggins 2006), CIT 6 (Kim et al. 2013), and o Ceti (Mayer et al. 2011). The most detailed observations so far were taken with the Atacama Large Millimeter/submillimeter Array (ALMA) of the CSE around the carbon star R Scl in molecular line emission (Maercker et al. 2012, hereafter M2012). They show a spiral structure extending from a detached shell inwards to the central star and present-day mass loss, indicating the presence of a binary companion shaping the wind. The asymmetries observed in the CSEs around AGB stars may be due to several processes. Non-isotropic mass loss, binary interaction, and outflows can shape the envelopes around AGB stars into asymmetric shapes (Mastrodemos & Morris 1999; Mauron & Huggins 2006; Kim & Taam 2012). Temporal variations in the mass loss will create rings and/or arc-like structures in the form of density enhancements in the envelope. An example of the most extreme cases of this can be seen in the detached shells around a number of carbon AGB stars, which are believed to be caused by the increase in mass-loss rate and expansion velocity during a thermal pulse (Olofsson et al. 1990, 1993, 1996; Steffen et al. 1998; Steffen & Schönberner 2000; Schöier et al. 2005; Mattsson et al. 2007; Maercker et al. 2010; Olofsson et al. 2010; Maercker et al. 2012).

The duration of a thermal pulse is a few hundred years, with timescales between pulses being several 10 000 years (Karakas & Lattanzio 2007). The chances of observing a pulse directly are extremely low. Currently less than around ten detached shells of gas and dust formed during a thermal pulse are known around carbon stars. Considering the typical time scales between thermal pulses and the lifetime of a detached shell of CO gas due to photodissociation, it is not likely that more sources of this kind will be found (Olofsson et al. 1990). Detailed observations of these detached shell sources may thus be the only way for us to effectively constrain the evolution of the star throughout the thermal pulse cycle. The spiral structure observed around R Scl in the ALMA observations allowed the evolution of the mass-loss rate and expansion velocity since the last thermal pulse to be measured (M2012) and is the only case where direct determination has been possible.

Direct imaging of the envelopes is difficult, in particular of the extended, dusty CSE. Thermal dust emission in the far infrared from detached shells has been observed by the AKARI instrument (Izumiura et al. 2011) and the PACS instrument (Kerschbaum et al. 2010; Mečina et al. 2014). PACS, in particular, delivered comparatively high-resolution images of a large number of AGB CSEs (Cox et al. 2012).

Observations of dust scattered stellar light have proven to be an effective way to obtain high-resolution images of the dusty CSEs. Hubble Space Telescope images of R Scl give the highest-resolution images of the detached shell of dust around this object to date (Olofsson et al. 2010, hereafter O2010), although only approximately one-third of the shell could be imaged. Images of polarised, dust-scattered stellar light add a powerful tool for investigating the distribution of the dust in the plane of the sky. This technique was used to successfully image the detached shells around R Scl and U Ant (González Delgado et al. 2003, hereafter GD2003; Maercker et al. 2010). The detached shells around DR Ser and U Cam, as well as the circumstellar environment around the S-type AGB star W Aql, were imaged in polarised, dust-scattered stellar light using the PolCor instrument mounted on the Nordic Optical Telescope (Ramstedt et al. 2011). This was the first direct image of the detached shell around the carbon star DR Ser. Until then the shell had only been inferred indirectly through the line profiles in molecular line observations (Schöier et al. 2005, hereafter S2005). The polarised, dust-scattered light observations allowed us to directly measure the radius and thickness of the shell around DR Ser.

We now present new observations with the PolCor instrument mounted on the ESO 3.6 m telescope in La Silla, Chile, of the dusty CSEs around two detached-shell AGB stars. The instrument provides images in polarised light combined with a coronographic mask. This allows us to block the direct, stellar light, making the faint, dust-scattered stellar light in the circumstellar medium visible. In Sect. 2 we present the observations, source properties, and data treatment, Sect. 3 gives an overview of the radiative transfer modelling, and Sects. 4 to 6 present and discuss the results, and give our conclusions.

2. Observations and data reduction

The observations were taken with the PolCor instrument mounted on the ESO 3.6 m telescope in La Silla, Chile, during 5 nights between October 18 and October 28 2011. We here present the observations of the carbon AGB stars R Scl taken on October 19 and V644 Sco taken on October 24, 26, and 27. They were part of a sample of 14 evolved stars to study the morphology of the dusty circumstellar environment around AGB stars. Observations of these sources were taken in the V and R bands centred at wavelengths of 0.55 μm and 0.64 μm, respectively. The pixel scale in the images is 0.̋114 per pixel. Typical values of the seeing were 1.̋3. All images were flat-fielded and dark-subtracted. PolCor reads out images at rates of 0.04  s and 0.1 s, allowing us to use the lucky imaging technique to shift-and-add the individual exposures, effectively improving the seeing in the image. Typically between 1020% of the images were rejected in the lucky imaging, resulting in a seeing of 0.̋9. In addition, the instrument takes images at four different polarisation angles (0, 45, 90, and 135 degrees) to determine the full set of Stokes parameters I, U, and Q, as well as the polarisation angle. The instrument is equipped with coronographic masks of varying sizes and densities to block the direct stellar light, making the faint scattered light visible. For the observations presented here a mask with a diameter of 3′′  was used. In order to more easily place the star in the centre of the mask, we used a mask that reduced the direct stellar light by 99%. Total on-source integration times (before image rejection in the lucky imaging) are listed in Table 1. The three nights of observations for V644 Sco were averaged. The PolCor instrument is described in detail in Ramstedt et al. (2011).

Table 1

PolCor observations with the 3.6 m telescope at La Silla observatory in October 2011.

2.1. Observed sources

2.1.1. R Scl

R Scl is a carbon-rich, semi-regular pulsating AGB star with a pulsation period of 370 days. The luminosity is estimated to 4300 L using a period-luminosity relation (Knapp et al. 2003), giving a distance of 290 pc. The parallax gives a distance of 266 pc. The detached shell around R Scl was first observed by Olofsson et al. (1990) in CO emission with the single-dish telescope SEST (Swedish-European Submillimeter Telescope). The gas has since been imaged in scattered stellar light in the resonance lines of Na and K (González Delgado et al. 2001, hereafter GD2001), and in high-spatial resolution interferometer maps in CO emission by ALMA (M2012). The dusty detached shell was observed in polarised, dust-scattered stellar light (GD2003), in dust-scattered stellar light with the Hubble Space Telescope (O2010), and in thermal dust emission with the PACS instrument aboard Herschel (Cox et al. 2012). All observations show a detached shell at 19′′20′′. The expansion velocity of the shell is measured to 14.5 km s-1(M2012), and the present-day expansion velocity to 10.5 km s-1  (S2005). The present-day mass-loss rate is estimated to be 3 × 10-7M yr-1 (S2005), while M2012 determine a mass-loss rate during the creation of the shell of 10-5M yr-1. The total estimated dust mass in the shell is very uncertain, and ranges between (2−3) × 10-6M (GD2003; O2010) to 3 × 10-5M (S2005). The gas mass is estimated to be 2.5 × 10-3M (S2005).

2.1.2. V644 Sco

V644 Sco is an irregularly pulsating carbon AGB star. Its distance is estimated to be 700 pc, based on SED modelling and assuming a luminosity of 4000 L (S2005). V644 Sco was first observed in CO emission using SEST by Olofsson et al. (1996), where a detached gas shell at less than 7′′ was found. The shells of dust and gas were modelled by S2005, giving estimates of the dust and gas masses of 1.4 × 10-4M  and 2.5 × 10-3M, respectively. They derive a radius for the detached gas shell of 10.̋5. Previous to this publication, the detached shell had not been observed directly.

2.2. Polarised dust-scattered stellar light

Stellar light is unpolarised unless it is intercepted by a medium and scattered. For dust grains the scattering of stellar light occurs in all directions (however, not necessarily isotropically), while the highest degree of polarisation is reached if the angle between the incoming light and the outgoing scattered light is at 90 degrees (e.g., Zubko & Laor 2000). For stellar light scattered by a dusty CSE surrounding the star, the polarised light hence shows the distribution of the dust in the plane of the sky. Observations at four different polarisation angles can be used to derive the full set of Stokes parameters I, U and Q, as well as the polarised intensity, P, and the polarisation degree, p (see, e.g., Maercker et al. 2010; Ramstedt et al. 2011).

We assume that interstellar space and Earth’s atmosphere have a negligible effect on the detected light, and that the telescope does not introduce a significant degree of polarisation (e.g., GD2003; Ramstedt et al. 2011). The measured images then contain contributions from only the star and the CSE: (1)(2)(3)(4)The direct stellar light is unpolarised, and the measured Stokes U and Q parameters, and hence also the polarised intensity, are equal to the contribution from the circumstellar medium only. Any residual polarisation added by the direct stellar light is generally behind and close to the edges of the mask where interpretation of the data already is difficult, and hence does not affect the analysis significantly.

2.3. PSF subtraction

In order to correct for the direct, stellar light in the total intensity images, template stars of similar magnitude, but without any (known) circumstellar emission were observed in the same setup as the science targets. The direct stellar contribution to the total intensity can then be subtracted by scaling the brightness of the template star to the science target (see also GD2001; Maercker et al. 2010; Ramstedt et al. 2011).

For several reasons the template subtraction is more complicated than anticipated: 1) the shape of the psf depends on the position of the star behind the mask; 2) the scaling of the template star to the science target is not trivial, as the CSE may contribute to the total intensity all the way to the edge of the mask; and 3) the shape of the psf can differ between science target and template star. Incorrect subtraction of the template star may result in significant artefacts in the final images. The observations of R Scl and V644 Sco were part of a sample of 14 sources to observe the dusty circumstellar structure around AGB stars. As such, the time spent on each object was limited, making a good psf subtraction difficult.

Not subtracting a template psf leads to a lower limit when calculating the polarisation degree of the circumstellar envelope. For V644 Sco template subtraction was not possible, and in all cases we present the polarisation degree as a lower limit. Due to the larger angular size of the detached shell around R Scl, template subtraction was somewhat simpler. In this case, in addition to calculating the lower limit, we also scaled the observed template star to the total intensity images and subtracted the psf.

Observations of flux template stars in principle allows absolute flux-calibration of the images. However, the focus of this study is the morphological information contained in the polarised light, and absolute flux calibration is therefore not necessary.

3. Dust radiative transfer analysis of scattered stellar light

We used the dust radiative transfer code radmc-3d (Dullemond 2012) in order to model the polarised, circumstellar emission towards R Scl and V644 Sco. Radmc-3d allows us to derive the full set of Stokes vectors for stellar light scattered by dust grains of any type of dust. We used amorphous carbon grains with optical constants from Suh (2000) with a grain density of 2 g/cm3 and a grain radius of 0.1 μm.

We modelled the polarised, stellar scattered light assuming a constant present-day mass-loss rate expanding at a constant velocity, giving a r-2 density profile. The detached shell is assumed to follow a Gaussian density profile in the radial direction with a radius R and a full width at half maximum (FWHM) of ΔR. The parameters for the present-day mass-loss rate and detached shell are taken from S2005. Due to large uncertainties in the dust properties (composition, size, shape), a detailed modelling of the dust to determine new dust masses and dust mass-loss-rates will not improve on previous results, and we instead concentrate on the detailed morphological information of the detached shells.

4. Results

thumbnail Fig. 1

Observations of R Scl in the R-band. Top left to right: total intensity, total intensity with template star subtracted, polarised intensity, and the lower limit of the polarisation degree derived using the measured total intensity (without psf subtraction). The images are smoothed by a Gaussian kernel with a FWHM of two pixels. The colour scale is given in counts/s. The bottom row shows the corresponding radial profiles. The total intensity radial profile shows R Scl (solid line) and the scaled template star (dashed line). The radial profile of the polarisation degree shows the polarisation degree using the measured total intensity (solid line) and the template-subtracted total intensity (dashed line).

thumbnail Fig. 2

Same as Fig. 1 but in the V-band.

Table 2

Measured shell radii (R) and shell widths (ΔR – FWHM of Gaussian distribution) for R Scl and V644 Sco.

Figures 1 to 4 show the PolCor observations for R Scl and V644 Sco in the R-band and V-band. For both sources the images of the observed total intensity, the total intensity after psf subtraction, the polarised intensity, and the lower limit of the polarisation degree are shown. For each of these we also present azimuthally averaged radial profiles (AARP). The profiles are averaged over regions that avoid artefacts introduced by the telescope spider. The detached shells are detected in the images of polarised emission at signal-to-noise ratios of S/N = 10−20. Due to the difficulty of the template star subtraction (Sect. 2.3), only a lower limit of the polarisation degree is estimated. Table 2 shows the determined shell radii and widths (defined as FWHM of a Gaussian distribution) for R Scl and V644 Sco.

The sizes and widths of the shells are determined by creating AARPs of the polarised intensity in the radmc-3d model images. The peak of the intensity at the position of the shell in the model AARP is then scaled to the peak in the observed AARP at the position of the shell. The radius of the shell and the width are varied until a satisfactory fit to the observed profile is achieved (Fig. 5). The radius is constrained by fitting the position of the peak, while the width is determined by fitting the model AARP to the outer edge of the detached shell in the observed AARP. Both the radii and widths can be determined to within ±0.̋5.

In principle, the increase in intensity towards the inner regions could be used to constrain the density contrast between the present-day mass-loss rate and the detached shells. However, effects due to the coronographic mask and a potential polarisation of the central star make any such estimate extremely uncertain.

4.1. R Scl

A shell radius for R Scl is determined from the dust-scattered polarised light images to 19.̋7 and a FWHM width of the Gaussian shell to 3.̋2. The polarisation degree reaches values of >20% at the position of the shell, showing that the dust is indeed located in a thin shell around the star. A comparison with determined shell radii and widths using different probes for the dust and gas is presented in Table 2. Our results are consistent with the radius derived by (GD2003). The determined radius also fits well with the radius determined in the HST images in O2010. We derive a larger FWHM for the shell than O2010. However, the HST images only cover 1 / 3 of the total shell. The region imaged in the HST data is dominated by a bright arc along the limb-brightened shell, while the radmc-3d models in this work are fit to the AARP of the entire shell, which may lead to a broader average width. In particular the apparent flattening of the southern part of the shell will lead to a slightly broader shell. The shell radius derived in the dust models in S2005 has a large uncertainty (± 14′′). These models depend critically on the assumed grain properties. If the grains are too small for example, the radius of the shell will be overestimated, since the grains have to be moved to larger distances to be cool enough to achieve the right temperature and give the correct infrared flux.

Imaging of the detached gas shell around R Scl in stellar light scattered in resonance lines of NaD and K (GD2001) gives results that are consistent with ours. The radius of the shell in the CO data from ALMA (M2012) is determined by measuring the peak of the emission in the AARP. This may lead to underestimating the size, due to limb-brightening along the inner edge of the shell, slightly moving the observed peak inwards compared to the density peak of the shell. This was already observed in M2010 and O2010. The width of the shell determined in ALMA data is limited by the shell being marginally resolved (with a beam-size of 1.̋2). Hence, the results presented here for the radius and the width of the dust shell are consistent with the measurements of the CO(3−2) emission line in the ALMA data.

In Fig. 6 (left) we compare the polarised intensity PolCor image in the R-band with the ALMA image of the CO(3−2) emission centred on the vlsr in a 0.5 km s-1 bin. The image of the polarised intensity shows the morphology of the dust in the detached shell in the plane of the sky, while the image at the stellar vlsr from ALMA shows the CO(3−2) emission in the plane of the sky. Assuming that the gas emission (roughly) traces the density distribution of the gas, the distribution of the dust and gas can be compared. The ALMA contours trace the PolCor image closely, showing that the spatial distribution of the dust and gas coincides. In particular, the flattening of the shell in the south is present in both images, as well as an apparent gap in the south-western part of the shell. The peaks of the emission also generally coincide between the dust and the gas (Fig. 7), although there are a few exceptions, most notably the emission peak in the molecular gas at a position angle PA = 120° where the dust has a clear gap. Nevertheless, the comparison between the PolCor and ALMA data suggests that the distribution of the dust and gas, including different clumps, is the same.

The total intensity PolCor image on the other hand should be directly comparable to the HST image observed by O2010. Figure 6 (right) shows the contours of the HST image in the f814 filter on top of the PolCor R-band total intensity image where the stellar psf has been subtracted. The psf subtraction in the HST data is of much higher quality, and structures observed close to the star are more reliable than for the PolCor data. The comparison between the two images shows that in particular for the shell the PolCor image shows the same features as the HST observations. The advantage of the PolCor data is that it covers the entire shell, while the limited field of view in the HST data misses potentially interesting features, such as the flattening of the shell in the southern part and an apparent gap to the south-west.

Except for the uncertain estimate from the dust modelling in S2005, all estimates of the size and width of the detached shell of dust and gas around R Scl give consistent results. The average value for the radius of the shell (taking all measurements into account) is 19.̋5±0.̋7. The width of the shell varies more. This however depends on the type of observation and measurement method. An average width of 2′′, with 1′′and 3′′as lower and upper limits for a varying width, appears reasonable. The observations of the detached shell presented here are the only images with high-spatial resolution of the entire detached shell of dust around R Scl, allowing us to measure the size of the shell with high accuracy.

thumbnail Fig. 3

Observations of V644 Sco in the R-band. Top left to right: total intensity, total intensity with template star subtracted, polarised intensity, and the lower limit of the polarisation degree derived using the measured total intensity (without psf subtraction). The images are smoothed by a Gaussian kernel with a FWHM of two pixels. The colour scale is given in counts/s. The bottom row shows the corresponding radial profiles. The total intensity radial profile shows V644 Sco (solid line) and the scaled template star (dashed line). The radial profile of the polarisation degree shows the polarisation degree using the measured total intensity only (solid line).

thumbnail Fig. 4

Same as Fig. 3 but in the V-band. The images are smoothed by a Gaussian kernel with a FWHM of four pixels.

4.2. V644 Sco

V644 Sco was observed in CO emission for the first time in Olofsson et al. (1996), showing the typical double-peaked line profile indicative of an expanding detached shell. They mapped the shell in the CO(3−2) transition in a 49-point map with a 16′′  beam and a 7.̋5 spacing. Although the emission appears to be partially resolved at velocities close to the vlsr, they conclude that the detached shell must have a radius of less that 7′′. S2005 determine the detached shell to have a radius of 10.̋5 (Table 2) based on molecular line modelling, and a radius of 17.̋0±7′′ based on the dust modelling. Also here the dust modelling suffers from large uncertainties. In the PolCor data (Figs. 3 and 4) the shell around V644 Sco can clearly be seen at a radius of 9.̋4. The width of a Gaussian shell fit to the AARP is 2.̋0. This is the first time the radius and width of the shell around V644 Sco have been measured directly, instead of being inferred from radiative transfer modelling. The radius derived from the molecular line modelling is remarkably close to the observed value in the PolCor data, lending strength to the results of the CO models in S2005.

The polarised data clearly show the detached shell, as well as structure/inhomogeneities within the shell. As with similar observations of other detached shell sources, V644 Sco shows a remarkable overall spherical symmetry. We only determine a lower limit of the polarisation degree due to the difficulties with the psf subtraction (Sect. 2.3). However, in all images the shell clearly appears with the degree of polarisation peaking at the position of the shell seen in the polarised images. As with R Scl, this implies that the dust is distributed in a thin shell with the determined radius and width. This indicates that the shell originated in a spherical symmetric mass loss event from the star, and evolved as an expanding density enhancement. The density enhancement may be retained through interaction with a surrounding spherically symmetric medium as described by Mattsson et al. (2007).

thumbnail Fig. 5

Radial profiles of the polarised intensity compared to the radiative transfer models for R Scl (top) and V644 Sco (bottom). The R-band and V-band data are plotted as blue and green lines, respectively, the model profiles as black lines. The profiles are normalised to 1.

5. Discussion

5.1. The detached shells around R Scl and V644 Sco

The PolCor images of R Scl allow a direct comparison with the ALMA observations of the CO(3−2) line emission (Fig. 6, left, and Fig. 7). The good correlation between the PolCor data and the HST and ALMA observations confirms that the observations reliably measure the radius and width of the observed detached shells. These are the most detailed observations of the shell of dust around R Scl to date. For V644 Sco we for the first time set reliable and direct constraints on the size and width of the detached shell of dust.

For U Ant, Maercker et al. (2010) find that the dust has (mostly) separated from the gas, with a detached shell dominated by dust outside a detached shell dominated by gas. This implies a drift velocity between the dust and gas, as is expected in a dust-driven wind. However, the comparison between the PolCor data and ALMA data for R Scl show a remarkable correlation between the dust and gas. The overall shape of the shell is reproduced in both observations, including the flatter southern part of the shell, while even individual clumps of dust and gas coincide. This implies a common evolution of the dust and gas around R Scl. Also for V644 Sco, U Cam, and DR Ser (O2010; Ramstedt et al. 2011) the derived sizes for the detached dust shells coincide with the sizes determined for the gaseous detached shell. The shells around R Scl, V644 Sco, U Cam, and DR Ser are 10002000 years younger than the shell around U Ant (S2005). Maercker et al. (2010) suggest a common origin of the dust and gas in the shells around U Ant during a thermal pulse, with a drift velocity of a few km s-1. Kerschbaum et al. (2010) only detect the inner detached shell in thermal dust emission in far-infrared PACS observations, while Arimatsu et al. (2011) detect both shells in mid-infrared observations. This implies the presence of small, hot grains in the inner shell, and larger grains in the outer shell which have drifted away from the gas, confirming the conclusions by Maercker et al. (2010). It is possible that the dust has not yet separated significantly from the younger detached shell sources. The separation of the dust and gas may also be affected by the surrounding medium, and hence the pre-pulse mass-loss rate.

The observed structures support the hypothesis in which the detached shells originate in the high mass-loss-rate phase during a thermal pulse. The shells then expand away from the star, while the declining post-pulse mass-loss rate causes the shell to appear as detached. Comparing the results of dust-scattered light one can conclude that the radii of the shells increase as a function of age as expected. Furthermore, with the exception of U Cam, the widths of the shells are similar (1000−2000 AU – note that while the angular sizes of the shells can be relatively well determined, the physical sizes are very dependent on the estimated distance to the sources). This is consistent with the creation of a detached shell during a mass-loss eruption and a subsequent interaction with a surrounding medium, keeping the absolute width of the shell constant (Steffen & Schönberner 2000).

5.2. The importance of constraining the shells

thumbnail Fig. 6

Comparison of the PolCor polarised intensity image for R Scl in the R-band with the ALMA band 7 data at the stellar vlsr (black contours, left) and the PolCor total intensity image with the Hubble Space Telescope data at 0.8 μm (black contours, right). The colour scale shows the PolCor images in counts/s.

thumbnail Fig. 7

Position angle (PA) vs. radius in the range from 15′′ to 11′′ showing the detached shell around R Scl. The colour image shows the polarised intensity from PolCor in the R-band, the contours the ALMA CO(3−2) emission at the stellar vlsr. The PA starts north and moves counter-clockwise.

Assuming distances of 290 pc to R Scl and 700 pc to V644 Sco, we obtain physical radii of the detached shells of 8.4 × 1016 cm and 9.8 × 1016 cm, respectively. The widths of the shells are 8.7 × 1016 cm for R Scl and 2.1 × 1016 cm for V644 Sco. Constraining the physical sizes of the shells is important to constrain hydrodynamical models. The age of the shell can be constrained by combining the velocity information in the molecular data with the information on the spatial correlation of the dust and gas. The radii of the shells depend on the age of the shell and the expansion velocity as a function of time. The width of the shell depends on the formation conditions during the thermal pulse and on the subsequent interaction with the circumstellar medium (Steffen & Schönberner 2000; Mattsson et al. 2007). Hydrodynamical models have to model the pre-pulse mass loss, the formation of the detached shell during a thermal pulse, and the evolution of the shell. The models are ultimately constrained by matching the observed parameters of the shells. The presented observations hence constrain the timescales and conditions of the pre-pulse, pulse, and post-pulse evolution (M2012).

The total mass in the detached shell likewise depends on the formation process and subsequent evolution of the shell. However, accurate estimates of the dust masses require high-quality observations at optical wavelengths in different bands (to constrain the grain size distribution of small grains), as well as observations at far-infrared and radio wavelengths to also constrain the population of large dust grains. The shapes and composition of the grains are also large uncertainties, and they have to be constrained by spectral observations. For some detached shell objects, far-infrared observations exist (Kerschbaum et al. 2010; Izumiura et al. 2011; Mečina et al. 2014). R Scl has been observed with the Large Bolometer Camera (LABOCA) on the Atacama Pathfinder Experiment (APEX) at 850 μm, and we are currently performing observations with LABOCA of three additional detached shell sources (V644 Sco, U Ant, and DR Ser) to probe the contribution of large dust grains to the total mass.

The drawback of observations in the far-infrared and radio is that the spatial resolution is relatively low. The optical observations in polarised light obtained here can be used to effectively constrain the size and width of the detached shells in the modelling of the thermal dust emission, leading to more accurate dust-mass estimates. In particular, a comparison to the PACS (Cox et al. 2012) and LABOCA observations of R Scl will determine the properties of the detached dust shell with similar detail as the observations of the detached gas shell.

6. Conclusions

The observations of polarised, dust-scattered stellar light clearly show the detached dust shells around the carbon AGB stars R Scl and V644 Sco. The data presented here are the most detailed images of the entire dust shell around R Scl to date. They allow us to measure the size and width of the detached shell of dust and determine the morphology of the shell with high precision. For V644 Sco this is the first time that the detached shell of dust was constrained directly.

For R Scl a direct comparison of the dust-scattered optical emission observed here with ALMA observations of CO(3−2) (M2012) and optical Hubble Space Telescope images (O2010) shows that the dust and gas in the detached shell coincide spatially to a high degree. This implies that the dust and gas have evolved together since the detached shell was created.

For V644 Sco this is the first time that the radius and width of the detached dust shell were measured directly. Comparison with radiative transfer models of the molecular emission (S2005) also indicates a close correlation between the dust and gas.

The creation of detached shells of dust and gas by thermal pulses has been described previously, both observationally and theoretically (e.g., Schöier et al. 2005; Mattsson et al. 2007). The results found here are in line with a scenario where the increase in mass-loss rate and expansion velocity during a thermal pulse leads to the formation of a detached shell of dust and gas. The derived radii and widths of the shells will help constrain more detailed hydrodynamical and radiative transfer models of the detached shells, and hence help constrain the formation and evolution of the shells since formation during the thermal pulse.

The results demonstrate that observations of polarised, dust-scattered stellar light can trace intricate details in the circumstellar envelopes around AGB stars. They are particularly useful to constrain the structures in the detached shells. These observations were part of a pilot sample, and hence only have a limited amount of time spent per object. Already they give new constraints on the detached shells around R Scl and, in particular, V644 Sco. Dedicated observations with PolCor and similar instruments of the known detached shell sources will help constrain the evolution of the star during and after a thermal pulse. With ALMA, correlating the evolution of the dust with the evolution of the molecular gas will give the best observational description of the thermal pulse phenomenon to date.

Acknowledgments

The research leading to these results has received funding from the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme (FP7/2007-2013) under REA grant agreement No. 623898.11. M.L.L.-F. was supported by the Deutscher Akademischer Austausch Dienst (DAAD). The authors would like to thank W.H.T Vlemmings, E. De Beck, and H. Olofsson for their contribution in discussing the results.

References

  1. Arimatsu, K., Izumiura, H., Ueta, T., Yamamura, I., & Onaka, T. 2011, ApJ, 729, L19 [NASA ADS] [CrossRef] [Google Scholar]
  2. Castro-Carrizo, A., Quintana-Lacaci, G., Neri, R., et al. 2010, A&A, 523, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Cox, N. L. J., Kerschbaum, F., van Marle, A.-J., et al. 2012, A&A, 537, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Decin, L., Cox, N. L. J., Royer, P., et al. 2012, A&A, 548, A113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Dullemond, C. P. 2012, RADMC-3D: A multi-purpose radiative transfer tool, astrophysics Source Code Library [Google Scholar]
  6. González Delgado, D., Olofsson, H., Schwarz, H. E., Eriksson, K., & Gustafsson, B. 2001, A&A, 372, 885 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. González Delgado, D., Olofsson, H., Schwarz, H. E., et al. 2003, A&A, 399, 1021 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Izumiura, H., Ueta, T., Yamamura, I., et al. 2011, A&A, 528, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Karakas, A., & Lattanzio, J. C. 2007, PASA, 24, 103 [NASA ADS] [CrossRef] [Google Scholar]
  10. Kerschbaum, F., Ladjal, D., Ottensamer, R., et al. 2010, A&A, 518, L140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Kim, H., & Taam, R. E. 2012, ApJ, 744, 136 [NASA ADS] [CrossRef] [Google Scholar]
  12. Kim, H., Hsieh, I.-T., Liu, S.-Y., & Taam, R. E. 2013, ApJ, 776, 86 [NASA ADS] [CrossRef] [Google Scholar]
  13. Knapp, G. R., Pourbaix, D., Platais, I., & Jorissen, A. 2003, A&A, 403, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Maercker, M., Olofsson, H., Eriksson, K., Gustafsson, B., & Schöier, F. L. 2010, A&A, 511, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Maercker, M., Mohamed, S., Vlemmings, W. H. T., et al. 2012, Nature, 490, 232 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  16. Mastrodemos, N., & Morris, M. 1999, ApJ, 523, 357 [NASA ADS] [CrossRef] [Google Scholar]
  17. Mattsson, L., Höfner, S., & Herwig, F. 2007, A&A, 470, 339 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Mauron, N., & Huggins, P. J. 2006, A&A, 452, 257 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Mayer, A., Jorissen, A., Kerschbaum, F., et al. 2011, A&A, 531, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Mečina, M., Kerschbaum, F., Groenewegen, M. A. T., et al. 2014, A&A, 566, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Olofsson, H., Carlstrom, U., Eriksson, K., Gustafsson, B., & Willson, L. A. 1990, A&A, 230, L13 [NASA ADS] [Google Scholar]
  22. Olofsson, H., Eriksson, K., Gustafsson, B., & Carlstrom, U. 1993, ApJS, 87, 267 [NASA ADS] [CrossRef] [Google Scholar]
  23. Olofsson, H., Bergman, P., Eriksson, K., & Gustafsson, B. 1996, A&A, 311, 587 [NASA ADS] [Google Scholar]
  24. Olofsson, H., Maercker, M., Eriksson, K., Gustafsson, B., & Schöier, F. 2010, A&A, 515, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Ramstedt, S., Maercker, M., Olofsson, G., Olofsson, H., & Schöier, F. L. 2011, A&A, 531, A148 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Schöier, F. L., Lindqvist, M., & Olofsson, H. 2005, A&A, 436, 633 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Steffen, M., & Schönberner, D. 2000, A&A, 357, 180 [NASA ADS] [Google Scholar]
  28. Steffen, M., Szczerba, R., & Schoenberner, D. 1998, A&A, 337, 149 [NASA ADS] [Google Scholar]
  29. Suh, K.-W. 2000, MNRAS, 315, 740 [NASA ADS] [CrossRef] [Google Scholar]
  30. Zubko, V. G., & Laor, A. 2000, ApJS, 128, 245 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

PolCor observations with the 3.6 m telescope at La Silla observatory in October 2011.

Table 2

Measured shell radii (R) and shell widths (ΔR – FWHM of Gaussian distribution) for R Scl and V644 Sco.

All Figures

thumbnail Fig. 1

Observations of R Scl in the R-band. Top left to right: total intensity, total intensity with template star subtracted, polarised intensity, and the lower limit of the polarisation degree derived using the measured total intensity (without psf subtraction). The images are smoothed by a Gaussian kernel with a FWHM of two pixels. The colour scale is given in counts/s. The bottom row shows the corresponding radial profiles. The total intensity radial profile shows R Scl (solid line) and the scaled template star (dashed line). The radial profile of the polarisation degree shows the polarisation degree using the measured total intensity (solid line) and the template-subtracted total intensity (dashed line).

In the text
thumbnail Fig. 2

Same as Fig. 1 but in the V-band.

In the text
thumbnail Fig. 3

Observations of V644 Sco in the R-band. Top left to right: total intensity, total intensity with template star subtracted, polarised intensity, and the lower limit of the polarisation degree derived using the measured total intensity (without psf subtraction). The images are smoothed by a Gaussian kernel with a FWHM of two pixels. The colour scale is given in counts/s. The bottom row shows the corresponding radial profiles. The total intensity radial profile shows V644 Sco (solid line) and the scaled template star (dashed line). The radial profile of the polarisation degree shows the polarisation degree using the measured total intensity only (solid line).

In the text
thumbnail Fig. 4

Same as Fig. 3 but in the V-band. The images are smoothed by a Gaussian kernel with a FWHM of four pixels.

In the text
thumbnail Fig. 5

Radial profiles of the polarised intensity compared to the radiative transfer models for R Scl (top) and V644 Sco (bottom). The R-band and V-band data are plotted as blue and green lines, respectively, the model profiles as black lines. The profiles are normalised to 1.

In the text
thumbnail Fig. 6

Comparison of the PolCor polarised intensity image for R Scl in the R-band with the ALMA band 7 data at the stellar vlsr (black contours, left) and the PolCor total intensity image with the Hubble Space Telescope data at 0.8 μm (black contours, right). The colour scale shows the PolCor images in counts/s.

In the text
thumbnail Fig. 7

Position angle (PA) vs. radius in the range from 15′′ to 11′′ showing the detached shell around R Scl. The colour image shows the polarised intensity from PolCor in the R-band, the contours the ALMA CO(3−2) emission at the stellar vlsr. The PA starts north and moves counter-clockwise.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.