Modeling the dust emission from PN IC418

We construct a detailed model for the IR dust emission from the PN IC 418. We succeed to reproduce the emission from 2 to 200$\mu$m. We can determine the amount of emitting dust as well as its composition, and compare to the depletion of elements determined for the photoionized region.Comment: Poster contribution (2 pages, 1 figure) to IAU Symposium 283: "Planetary Nebulae: An Eye to the Future" held in Puerto de la Cruz, Tenerife, Spain in July 25th-29th 2011. Few typos correcte

Modeling dust emission in PN IC 418

We investigated the infrared (IR) dust emission from PN IC 418, using a detailed model controlled by a previous determination of the stellar properties and the characteristics of the photoionized nebula, keeping as free parameters the dust types, amounts and distributions relative to the distance of the central star. The model includes the ionized region and the neutral region beyond the recombination front (Photodissociation region, or PDR), where the [OI] and [CII] IR lines are formed. We succeeded in reproducing the observed infrared emission from 2 to 200~\mm. The global energy budget is fitted by summing up contributions from big grains of amorphous carbon located in the neutral region and small graphite grains located in the ionized region (closer to the central star). Two emission features seen at 11.5 and 30~\mm are also reproduced by assuming them to be due to silicon carbide (SiC) and magnesium and iron sulfides (Mg$_x$Fe$_{1-x}$S), respectively. For this, we needed to consider ellipsoidal shapes for the grains to reproduce the wavelength distribution of the features. Some elements are depleted in the gaseous phase: Mg, Si, and S have sub-solar abundances (-0.5 dex below solar by mass), while the abundance of C+N+O+Ne by mass is close to solar. Adding the abundances of the elements present in the dusty and gaseous forms leads to values closer to but not higher than solar, confirming that the identification of the feature carriers is plausible. Iron is strongly depleted (3 dex below solar) and the small amount present in dust in our model is far from being enough to recover the solar value. A remaining feature is found as a residue of the fitting process, between 12 and 25~\mm, for which we do not have identification.Comment: Accepted for publication in Astronomy & Astrophysics. V2: adding reference

Morphological classification of post-AGB stars

We present a complete study of the morphology of post-Asymptotic Giant Branch (AGB) stars. Post-AGB is a very short evolutionary phase between the end of the AGB and the beginning of the Planetary Nebula (PN) stage (between 100 and 10,000 yrs). We have defined the end of the post-AGB phase and the beginning of the PN phase when the star is hot enough to fully ionize the hydrogen envelope. Post-AGB stars have a circumstellar shell that is illuminated by the central stars or partially ionized. However, this circumstellar shell is too small to be resolved from ground-based observations. Thus, we have used data from the Hubble Space Telescope (HST) database to resolve these shells. About 150 post-AGB were found in this database. Here we present the preliminary results on their morphological classification and the correlation with several parameters such as galactic latitude and IRAS fluxes. Our preliminary results show that 40% of the sample are stellar-like (S), 33 % bipolar (B), 12 % multi-polar (M) and 15 % elliptical (E).Comment: proceedings of the conference "Why Galaxies Care About AGB Stars II", Vienna 2010, eds. Franz Kerschbaum, Thomas Lebzelter and Bob Wing, ASP Conf.Ser (in press

New groups of planetary nebulae with peculiar dust chemistry towards the Galactic bulge

We investigate Galactic bulge planetary nebulae without emission-line central stars for which peculiar infrared spectra have been obtained with the Spitzer Space Telescope, including the simultaneous signs of oxygen and carbon based dust. Three separate sub-groups can be defined characterized by the different chemical composition of the dust and the presence of crystalline and amorphous silicates. We find that the classification based on the dust properties is reflected in the more general properties of these planetary nebulae. However, some observed properties are difficult to relate to the common view of planetary nebulae. In particular, it is challenging to interpret the peculiar gas chemical composition of many analyzed objects in the standard picture of the evolution of planetary nebulae progenitors. We confirm that the dual-dust chemistry phenomenon is not limited to planetary nebulae with emission-line central stars.Comment: 17 pages, 13 figure

Detection of a Far-Infrared Bow-Shock Nebula Around R Hya: the First MIRIAD Results

We present the first results of the MIRIAD (MIPS [Multiband Imaging Photometer for Spitzer] Infra-Red Imaging of AGB [asymptotic giant branch] Dustshells) project using the Spitzer Space Telescope. The primary aim of the project is to probe the material distribution in the extended circumstellar envelopes (CSE) of evolved stars and recover the fossil record of their mass loss history. Hence, we must map the whole of the CSEs plus the surrounding sky for background subtraction, while avoiding the central star that is brighter than the detector saturation limit. With our unique mapping strategy, we have achieved better than one MJy/sr sensitivity in three hours of integration and successfully detected a faint (< 5 MJy/sr), extended (~400 arcsec) far-infrared nebula around the AGB star R Hya. Based on the parabolic structure of the nebula, the direction of the space motion of the star with respect to the nebula shape, and the presence of extended H alpha emission co-spatial to the nebula, we suggest that the detected far-IR nebula is due to a bow shock at the interface of the interstellar medium and the AGB wind of this moving star. This is the first detection of the stellar-wind bow-shock interaction for an AGB star and exemplifies the potential of Spitzer as a tool to examine the detailed structure of extended far-IR nebulae around bright central sources. \Comment: 10 pages, 2 figures, accepted for publication in ApJ

Radiative transfer models of non-spherical prestellar cores

We present 2D Monte Carlo radiative transfer simulations of prestellar cores. We consider two types of asymmetry: disk-like asymmetry, in which the core is denser towards the equatorial plane than towards the poles; and axial asymmetry, in which the core is denser towards the south pole than the north pole. We limit our treatment to cores with mild asymmetries, which are exposed directly to the interstellar radiation field or are embedded inside molecular clouds. The isophotal maps of a core depend strongly on the viewing angle. Maps at wavelengths longer than the peak of the SED (e.g. 850 micron) essentially trace the column-density. Thus, for instance, cores with disk-like asymmetry appear elongated when mapped at 850 micron from close to the equatorial plane. However, at wavelengths near the peak of the SED (e.g. 200 micron), the emissivity is more strongly dependent on the temperature, and therefore, at particular viewing angles, there are characteristic features which reflect a more complicated convolution of the density and temperature fields within the core. These characteristic features are on scales 1/5 to 1/3 of the overall core size, and so high resolution observations are needed to observe them. They are also weaker if the core is embedded in a molecular cloud (because the range of temperature within the core is then smaller), and so high sensitivity is needed to detect them. Herschel, to be launched in 2007, will in principle provide the necessary resolution and sensitivity at 170 to 250 micron.Comment: 16 pages, 22 figures, accepted by A&A, also available (with high resolution figures) at http://www.astro.cf.ac.uk/pub/Dimitrios.Stamatellos/publications

On the difference between type E and A OH/IR stars

The observed SEDs of a sample of 60 OH/IR stars are fitted using a radiative transfer model of a dusty envelope. Among the whole sample, 21 stars have reliable phase-lag distances while the others have less accurate distances. L*-P,Mlr-P and Mlr-L* relations have been plotted for these stars. It is found that type E (with emission feature at 10um and type A (with absorption feature at 10um) OH/IR stars have different L*-P and Mlr-L* relations while both of them follow a single Mlr-P relation. The type E stars are proven to be located in the area without large scale dense interstellar medium while the type A stars are located probably in dense interstellar medium. It is argued here that this may indicate the two types of OH/IR stars have different chemical composition or zero age main sequence mass and so evolve in different ways. This conclusion has reinforced the argument by Chen et al.(2001) who reached a similar conclusion from the galactic distribution of about 1000 OH/IR stars with the IRAS low-resolution spectra (LRS).Comment: 6 pages, 9 figures, 2 table

The wind of W Hya as seen by Herschel. II. The molecular envelope of W Hya

The evolution of low- and intermediate-mass stars on the asymptotic giant branch (AGB) is mainly controlled by the rate at which these stars lose mass in a stellar wind. Understanding the driving mechanism and strength of the stellar winds of AGB stars and the processes enriching their surfaces with products of nucleosynthesis are paramount to constraining AGB evolution and predicting the chemical evolution of galaxies. In a previous paper we have constrained the structure of the outflowing envelope of W Hya using spectral lines of the $^{12}$CO molecule. Here we broaden this study by modelling an extensive set of H$_{2}$O and $^{28}$SiO lines observed by the three instruments on board Herschel using a state-of-the-art molecular excitation and radiative transfer code. The oxygen isotopic ratios and the $^{28}$SiO abundance profile can be connected to the initial stellar mass and to crucial aspects of dust formation at the base of the stellar wind, respectively. The modelling of H$_{2}$O and $^{28}$SiO confirms the properties of the envelope model of W Hya derived from $^{12}$CO lines. We find an H$_2$O ortho-to-para ratio of 2.5\,$^{+2.5}_{-1.0}$, consistent with what is expected for an AGB wind. The O$^{16}$/O$^{17}$ ratio indicates that W Hya has an initial mass of about 1.5 M$_\odot$. Although the ortho- and para-H$_{2}$O lines observed by HIFI appear to trace gas of slightly different physical properties, a turbulence velocity of $0.7\pm0.1$ km s$^{-1}$ fits the HIFI lines of both spin isomers and those of $^{28}$SiO well. The ortho- and para-H$_2^{16}$O and $^{28}$SiO abundances relative to H$_{2}$ are $(6^{+3}_{-2}) \times 10^{-4}$, $(3^{+2}_{-1}) \times 10^{-4}$, and $(3.3\pm 0.8)\times 10^{-5}$, respectively. Assuming a solar silicon-to-carbon ratio, the $^{28}$SiO line emission model is consistent with about one-third of the silicon atoms being locked up in dust particles