We performed a far-IR imaging survey of the circumstellar dust shells of 144 evolved stars as a mission programme of the AKARI infrared astronomical satellite using the Far-Infrared Surveyor (FIS) instrument. With this survey, we deliver far-IR surface brightness distributions of roughly 10' x 40' or 10' x 20' areas of the sky around the target evolved stars in the four FIS bands at 65, 90, 140, and 160 microns. Our objectives are to characterize the far-IR surface brightness distributions of the cold dust component in the circumstellar dust shells, from which we derive the amount of cold dust grains as low as 20 K and empirically establish the history of the early mass loss history. In this first installment of the series, we introduce the project and its aims, describe the observations, data reduction, and surface brightness correction process, and present the entire data set along with the results of integrated photometry measurements (i.e., the central source and circumstellar dust shell altogether). We find that (1) far-IR emission is detected from all but one object at the spatial resolution about 30" - 50" in the corresponding bands, (2) roughly 60 - 70 % of the target sources show some extension, (3) previously unresolved nearby objects in the far-IR are now resolved around 28 target sources, (4) the results of photometry measurements are reasonable with respect to the entries in the AKARI/FIS Bright Source Catalogue, despite the fact that the targets are assumed to be point-sources when catalogue flux densities were computed, and (5) an IR two-color diagram would place the target sources in a roughly linear distribution that may correlate with the age of the circumstellar dust shell and can potentially be used to identify which targets are more extended than others.

The spatial structure of the emission lines and continuum over the 50 arcsecond extent of the nearby, O-rich, planetary nebula NGC 7009 (Saturn Nebula) have been observed with the MUSE integral field spectrograph on the ESO Very Large Telescope. Science Verification data, in <0.6 arcsecond seeing, have been reduced and analysed as images over the wavelength range 4750-9350A. Emission line maps over the bright shells are presented, from neutral to the highest ionization available (He II and [Mn V]). For collisionally excited lines (CELs), maps of electron temperature (T_e from [N II] and [S III]) and electron density (N_e from [S II] and [Cl III]) are available and for optical recombination lines (ORLs) temperature (from the Paschen jump and ratio of He I lines) and density (from high Paschen lines). These estimates are compared: for the first time, maps of the differences in CEL and ORL T_e's have been derived, and correspondingly a map of t^2 between a CEL and ORL temperature, showing considerable detail. Total abundances of He and O were formed, the latter using three ionization correction factors. However the map of He/H is not flat, departing by ~2% from a constant value, with remnants corresponding to ionization structures. Ionization correction factor methods are compared for O abundance, but none delivers a flat map. An integrated spectrum over an area of 2340 square arcseconds was also formed and compared to 1D photoionization models. The spatial variation of a range of nebular parameters illustrates the complexity of the ionized media in NGC 7009. These MUSE data are very rich with detections of many lines over areas of hundreds of square arcseconds and follow-on studies are indicated. (Abridged)

We observed two proto-planetary nebulae, HD 56126 representing a source with an elliptical circumstellar shell and IRAS 16594-4656 representing a source with a bipolar circumstellar shell, with ALMA in the 12CO and 13CO J=3-2 lines and neighboring continuum to see how the morpho-kinematics of CO gas and dust emission properties in their circumstellar environments differ.

The AKARI Infrared Astronomical Satellite produced the all-sky survey (AFASS) maps in the far-IR at roughly arc-minute spatial resolution, enabling us to investigate the whole sky in the far-IR for objects having surface brightnesses greater than a few to a couple of dozen MJy/sr. While the AFASS maps are absolutely calibrated against large-scale diffuse emission, it was uncertain whether or not an additional flux correction for point sources was necessary. Here, we verify that calibration for point-source photometry in the AFASS maps is proper. With the aperture correction method based on the empirical point-spread-function templates derived directly from the AFASS maps, fluxes in the AKARI bright source catalogue (BSC) are reproduced. The AKARI BSC fluxes are also satisfactorily recovered with the 1 sigma aperture, which is the empirical equivalent of an infinite aperture. These results confirm that in the AFASS maps far-IR photometry can be properly performed by using the aperture correction method for point sources and by summing all pixel values within an appropriately defined aperture of the intended target (i.e., the aperture photometry method) for extended sources.

The bipolar nebula Menzel 3 (Mz 3) was observed as part of the \textit{Herschel} Planetary Nebula Survey (\textit{HerPlaNS}), which used the PACS and SPIRE instruments aboard the \textit{Herschel Space Observatory} to study a sample of planetary nebulae (PNe). In this paper, one of the series describing \textit{HerPlaNS} results, we report the detection of H I recombination lines (HRLs) in the spectrum of Mz 3. Inspection of the spectrum reveals the presence of 12 HRLs in the 55 to 680 $\mu$m range covered by the PACS and SPIRE instruments (H11$\alpha$ to H21$\alpha$ and H14$\beta$). The presence of HRLs in this range is unusual for PNe and has not been reported in Mz 3 before. Our analysis indicates that the HRLs we observed are enhanced by laser effect occurring in the core of Mz 3. Our arguments for this are: (i) the available Mz 3 optical to submillimetre HRL $\alpha$ line intensity ratios are not well reproduced by the spontaneous emission of optically thin ionized gas, as would be typical for nebular gas in PNe; (ii) the compact core of Mz 3 is responsible for a large fraction of the Herschel HRLs emission; (iii) the line intensity ratios for Mz 3 are very similar to those in the core emission of the well known star MWC 349A, where laser effect is responsible for the enhancement of HRLs in the Herschel wavelength range; (iv) the physical characteristics relevant to cause laser effect in the core of MWC 349A are very similar to those in the core of Mz 3.

We perform a comprehensive analysis of the planetary nebula (PN) NGC6781 to investigate the physical conditions of each of its ionized, atomic, and molecular gas and dust components and the object's evolution, based on panchromatic observational data ranging from UV to radio. Empirical nebular elemental abundances, compared with theoretical predictions via nucleosynthesis models of asymptotic giant branch (AGB) stars, indicate that the progenitor is a solar-metallicity, 2.25-3.0 Msun initial-mass star. We derive the best-fit distance of 0.46 kpc by fitting the stellar luminosity (as a function of the distance and effective temperature of the central star) with the adopted post-AGB evolutionary tracks. Our excitation energy diagram analysis indicate high excitation temperatures in the photodissociation region (PDR) beyond the ionized part of the nebula, suggesting extra heating by shock interactions between the slow AGB wind and the fast PN wind. Through iterative fitting using the Cloudy code with empirically-derived constraints, we find the best-fit dusty photoionization model of the object that would inclusively reproduce all of the adopted panchromatic observational data. The estimated total gas mass (0.41 Msun) corresponds to the mass ejected during the last AGB thermal pulse event predicted for a 2.5 Msun initial-mass star. A significant fraction of the total mass (about 70 percent) is found to exist in the PDR, demonstrating the critical importance of the PDR in PNe that are generally recognized as the hallmark of ionized/H+ regions.

We present a general surface brightness correction method for compact extended sources imaged in the slow-scan pointed observation mode of the Far-Infrared Surveyor (FIS) aboard the AKARI Infrared Astronomical Satellite. Our method recovers correct surface brightness distribution maps by re-scaling archived raw FIS maps using the surface-brightness-dependent inverse FIS response function. The flux of a target source is then automatically corrected for as the simple sum of surface brightnesses within the adopted contour encircling the perimeter of the target (i.e., contour photometry). This correction method is contrasted to the previous aperture photometry method for point sources, which directly corrects for the target flux with a flux-dependent scaling law. The new surface brightness correction scheme is applicable to objects of any shape from unresolved point sources to resolved extended objects, as long as the target is not deemed diffuse, i.e., the total extent of the target source does not exceed too much more than a single FIS scan width of 10 arcmin. The new correction method takes advantage of the well-defined shape (i.e., the scale invariance) of the point-spread function, which enables us to adopt a power-law FIS response function. We analyze the point-source photometric calibrator data using the FIS AKARI Slow-scan Tool (FAST) and constrained the parameters of the adopted power-law FIS response function. We conclude that the photometric accuracy of the new correction method is better than 10% error based on comparisons with the expected fluxes of the photometric calibrators and that resulting fluxes without the present correction method can lead up to 230% overestimates or down to 50% underestimates.

The large field and wavelength range of MUSE is well suited to mapping Galactic planetary nebulae (PN). The bright PN NGC 7009 was observed with MUSE on the VLT during the Science Verification of the instrument in seeing of 0.6". Emission line maps in hydrogen Balmer and Paschen lines were formed from analysis of the MUSE cubes. The measured electron temperature and density from the MUSE cube were employed to predict the theoretical hydrogen line ratios and map the extinction distribution across the nebula. After correction for the interstellar extinction to NGC 7009, the internal dust-to-gas ratio (A_V/N_H) has been mapped for the first time in a PN. The extinction map of NGC 7009 has considerable structure, broadly corresponding to the morphological features of the nebula. A large-scale feature in the extinction map, consisting of a crest and trough, occurs at the rim of the inner shell. The nature of this feature was investigated and instrumental and physical causes considered; no convincing mechanisms were identified to produce this feature, other than mass loss variations in the earlier asymptotic giant branch phase. The dust-to-gas ratio A_V/N_H increases from 0.7 times the interstellar value to >5 times from the centre towards the periphery of the ionized nebula. The integrated A_V/N_H is about 2 times the mean ISM value. It is demonstrated that extinction mapping with MUSE provides a powerful tool for studying the distribution of PN internal dust and the dust-to-gas ratio. (Abridged.)

In our series of papers presenting the Herschel imaging of evolved planetary nebulae, we present images of the dust distribution in the Helix nebula (NGC 7293). Images at 70, 160, 250, 350, and 500 micron were obtained with the PACS and SPIRE instruments on board the Herschel satellite. The broadband maps show the dust distribution over the main Helix nebula to be clumpy and predominantly present in the barrel wall. We determined the spectral energy distribution of the main nebula in a consistent way using Herschel, IRAS, and Planck flux values. The emissivity index of 0.99 +/- 0.09, in combination with the carbon rich molecular chemistry of the nebula, indicates that the dust consists mainly of amorphous carbon. The dust excess emission from the central star disk is detected at 70 micron and the flux measurement agree with previous measurement. We present the temperature and dust column density maps. The total dust mass across the Helix nebula (without its halo) is determined to be 0.0035 solar mass at a distance of 216 pc. The temperature map shows dust temperatures between 22 and 42 K, which is similar to the kinetic temperature of the molecular gas, strengthening the fact that the dust and gas co-exist in high density clumps. Archived images are used to compare the location of the dust emission in the far infrared (Herschel) with the ionized (GALEX, Hbeta) and molecular hydrogen component. The different emission components are consistent with the Helix consisting of a thick walled barrel-like structure inclined to the line of sight. The radiation field decreases rapidly through the barrel wall.

We present X-ray spectral analysis of 20 point-like X-ray sources detected in Chandra Planetary Nebula Survey (ChanPlaNS) observations of 59 planetary nebulae (PNe) in the solar neighborhood. Most of these 20 detections are associated with luminous central stars within relatively young, compact nebulae. The vast majority of these point-like X-ray-emitting sources at PN cores display relatively "hard" ($\geq0.5$~keV) X-ray emission components that are unlikely to be due to photospheric emission from the hot central stars (CSPN). Instead, we demonstrate that these sources are well modeled by optically-thin thermal plasmas. From the plasma properties, we identify two classes of CSPN X-ray emission: (1) high-temperature plasmas with X-ray luminosities, $L_{\rm X}$, that appear uncorrelated with the CSPN bolometric luminosity, $L_{\rm bol}$; and (2) lower-temperature plasmas with $L_{\rm X}/L_{\rm bol}\sim10^{-7}$. We suggest these two classes correspond to the physical processes of magnetically active binary companions and self-shocking stellar winds, respectively. In many cases this conclusion is supported by corroborative multiwavelength evidence for the wind and binary properties of the PN central stars. By thus honing in on the origins of X-ray emission from PN central stars, we enhance the ability of CSPN X-ray sources to constrain models of PN shaping that invoke wind interactions and binarity.

We present results from the most recent set of observations obtained as part of the Chandra X-ray observatory Planetary Nebula Survey (ChanPlaNS), the first comprehensive X-ray survey of planetary nebulae (PNe) in the solar neighborhood (i.e., within ~1.5 kpc of the Sun). The survey is designed to place constraints on the frequency of appearance and range of X-ray spectral characteristics of X-ray-emitting PN central stars and the evolutionary timescales of wind-shock-heated bubbles within PNe. ChanPlaNS began with a combined Cycle 12 and archive Chandra survey of 35 PNe. ChanPlaNS continued via a Chandra Cycle 14 Large Program which targeted all (24) remaining known compact (R_neb <~ 0.4 pc), young PNe that lie within ~1.5 kpc. Results from these Cycle 14 observations include first-time X-ray detections of hot bubbles within NGC 1501, 3918, 6153, and 6369, and point sources in HbDs 1, NGC 6337, and Sp 1. The addition of the Cycle 14 results brings the overall ChanPlaNS diffuse X-ray detection rate to ~27% and the point source detection rate to ~36%. It has become clearer that diffuse X-ray emission is associated with young (<~5x10^3 yr), and likewise compact (R_neb<~0.15 pc), PNe with closed structures and high central electron densities (n_e>~1000 cm^-3), and rarely associated with PNe that show H_2 emission and/or pronounced butterfly structures. Hb 5 is one such exception of a PN with a butterfly structure that hosts diffuse X-ray emission. Additionally, of the five new diffuse X-ray detections, two host [WR]-type CSPNe, NGC 1501 and NGC 6369, supporting the hypothesis that PNe with central stars of [WR]-type are likely to display diffuse X-ray emission.

We report the first detections of OH$^+$ emission in planetary nebulae (PNe). As part of an imaging and spectroscopy survey of 11 PNe in the far-IR using the PACS and SPIRE instruments aboard the Herschel Space Observatory, we performed a line survey in these PNe over the entire spectral range between 51 and 672$\mu$m to look for new detections. OH$^+$ rotational emission lines at 152.99, 290.20, 308.48, and 329.77$\mu$m were detected in the spectra of three planetary nebulae: NGC 6445, NGC 6720, and NGC 6781. Excitation temperatures and column densities derived from these lines are in the range of 27 to 47 K and 2$\times$10$^{10}$ to 4 $\times$10$^{11}$ cm$^{-2}$, respectively. In PNe, the OH+ rotational line emission appears to be produced in the photodissociation region (PDR) in these objects. The emission of OH+ is observed only in PNe with hot central stars (T$_{eff}$ > 100000 K), suggesting that high-energy photons may play a role in the OH+ formation and its line excitation in these objects, as it seems to be the case for ultraluminous galaxies.

The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyze the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. For the first time, we have detected extended OH+ emission in a planetary nebula. The spectra towards the Helix Nebula also show CO emission lines (from J= 4 to 8), [NII] at 1461 GHz from ionized gas, and [CI] (2-1), which together with the OH+ lines, trace extended CO photodissociation regions along the rims. The estimated OH+ column density is (1-10)x1e12 cm-2. The CH+ (1-0) line was not detected at the sensitivity of our observations. Non-LTE models of the CO excitation were used to constrain the average gas density (n(H2)=(1-5)x1e5 cm-3) and the gas temperature (Tk= 20-40 K). The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix Nebula whereas OH+ and [CI] lines trace the diffuse gas and the UV and X-ray illuminated clumps surface where molecules reform after CO photodissociation. [NII] traces a more diffuse ionized gas component in the interclump medium.

We present 48 Herschel/PACS spectra of evolved stars in the wavelength range of 67-72 $\mu$m. This wavelength range covers the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{2-2x}\text{Fe}_{(2x)}\text{SiO}_{4}$). The width and wavelength position of this band are sensitive to the temperature and composition of the crystalline olivine. Our sample covers a wide range of objects: from high mass-loss rate AGB stars (OH/IR stars, $\dot M \ge 10^{-5}$ M$_\odot$/yr), through post-AGB stars with and without circumbinary disks, to planetary nebulae and even a few massive evolved stars. The goal of this study is to exploit the spectral properties of the 69 $\mu$m band to determine the composition and temperature of the crystalline olivine. Since the objects cover a range of evolutionary phases, we study the physical and chemical properties in this range of physical environments. We fit the 69 $\mu$m band and use its width and position to probe the composition and temperature of the crystalline olivine. For 27 sources in the sample, we detected the 69 $\mu$m band of crystalline olivine ($\text{Mg}_{(2-2x)}\text{Fe}_{(2x)}\text{SiO}_{4}$). The 69 $\mu$m band shows that all the sources produce pure forsterite grains containing no iron in their lattice structure. The temperature of the crystalline olivine as indicated by the 69 $\mu$m band, shows that on average the temperature of the crystalline olivine is highest in the group of OH/IR stars and the post-AGB stars with confirmed Keplerian disks. The temperature is lower for the other post-AGB stars and lowest for the planetary nebulae. A couple of the detected 69 $\mu$m bands are broader than those of pure magnesium-rich crystalline olivine, which we show can be due to a temperature gradient in the circumstellar environment of these stars. continued...

This is the first of a series of investigations into far-IR characteristics of 11 planetary nebulae (PNs) under the Herschel Space Observatory Open Time 1 program, Herschel Planetary Nebula Survey (HerPlaNS). Using the HerPlaNS data set, we look into the PN energetics and variations of the physical conditions within the target nebulae. In the present work, we provide an overview of the survey, data acquisition and processing, and resulting data products. We perform (1) PACS/SPIRE broadband imaging to determine the spatial distribution of the cold dust component in the target PNs and (2) PACS/SPIRE spectral-energy-distribution (SED) and line spectroscopy to determine the spatial distribution of the gas component in the target PNs. For the case of NGC 6781, the broadband maps confirm the nearly pole-on barrel structure of the amorphous carbon-richdust shell and the surrounding halo having temperatures of 26-40 K. The PACS/SPIRE multi-position spectra show spatial variations of far-IR lines that reflect the physical stratification of the nebula. We demonstrate that spatially-resolved far-IR line diagnostics yield the (T_e, n_e) profiles, from which distributions of ionized, atomic, and molecular gases can be determined. Direct comparison of the dust and gas column mass maps constrained by the HerPlaNS data allows to construct an empirical gas-to-dust mass ratio map, which shows a range of ratios with the median of 195+-110. The present analysis yields estimates of the total mass of the shell to be 0.86 M_sun, consisting of 0.54 M_sun of ionized gas, 0.12 M_sun of atomic gas, 0.2 M_sun of molecular gas, and 4 x 10^-3 M_sun of dust grains. These estimates also suggest that the central star of about 1.5 M_sun initial mass is terminating its PN evolution onto the white dwarf cooling track.

As part of the Herschel Guaranteed Time Key Project MESS (Mass loss of Evolved StarS) we have imaged a sample of planetary nebulae. In this paper we present the PACS and SPIRE images of the classical bipolar planetary nebula NGC 650. We used these images to derive a temperature map of the dust. We also constructed a photoionization and dust radiative transfer model using the spectral synthesis code Cloudy. To constrain this model, we used the PACS and SPIRE fluxes and combined these with hitherto unpublished IUE and Spitzer IRS spectra as well as various other data from the literature. The temperature map combined with the photoionization model were used to study various aspects of the central star, the nebula, and in particular the dust grains in the nebula. The central star parameters are determined to be T_eff = 208 kK and L = 261 L_sol assuming a distance of 1200 pc. The stellar temperature is much higher than previously published values. We confirm that the nebula is carbon-rich with a C/O ratio of 2.1. The nebular abundances are typical for a type IIa planetary nebula. With the photoionization model we determined that the grains in the ionized nebula are large (assuming single-sized grains, they would have a radius of 0.15 micron. Most likely these large grains were inherited from the asymptotic giant branch phase. The PACS 70/160 micron temperature map shows evidence for two radiation components heating the grains. The first component is direct emission from the central star, while the second component is diffuse emission from the ionized gas (mainly Ly alpha). We show that previous suggestions that there is a photo-dissociation region surrounding the ionized region are incorrect. The neutral material resides in dense clumps inside the ionized region. These may also harbor stochastically heated very small grains in addition to the large grains.

We present the results of differential proper-motion analyses of the Egg Nebula (RAFGL 2688, V1610 Cyg) based on the archived two-epoch optical data taken with the Hubble Space Telescope. First, we determined that the polarization characteristics of the Egg Nebula is influenced by the higher optical depth of the central regions of the nebula (i.e., the "dustsphere" of about 1000 AU radius), causing the nebula illuminated in two steps -- the direct starlight is first channeled into bipolar cavities and then scattered off to the rest of the nebula. We then measured the amount of motion of local structures and the signature concentric arcs by determining their relative shifts over the 7.25 yr interval. Based on our analysis, which does not rely on the single-scattering assumption, we concluded that the lobes have been excavated by a linear expansion along the bipolar axis for the past 400 yr, while the concentric arcs have been generated continuously and moving out radially at about 10 km/s for the past 5,500 yr, and there appears to be a colatitudinally-increasing trend in the radial expansion velocity field of the concentric arcs. There exist numerical investigations into the mass-loss modulation by the central binary system, which predict such a colatitudinally-increasing expansion velocity field in the spiral-shock trails of the mass-loss ejecta. Therefore, the Egg Nebula may represent a rare edge-on case of the binary-modulated circumstellar environs, corroborating the previous theoretical predictions.

CW Leo has been observed six times between October 2009 and June 2012 with the SPIRE instrument on board the Herschel satellite. Variability has been detected in the flux emitted by the central star with a period of 639 \pm 4 days, in good agreement with determinations in the literature. Variability is also detected in the bow shock around CW Leo that had previously been detected in the ultraviolet and Herschel PACS/SPIRE data. Although difficult to prove directly, our working hypothesis is that this variability is directly related to that of the central star. In this case, fitting a sine curve with the period fixed to 639 days results in a time-lag in the variability between bow shock and the central star of 402 \pm 37 days. The orientation of the bow shock relative to the plane of the sky is unknown (but see below). For an inclination angle of zero degrees, the observed time-lag translates into a distance to CW Leo of 130 \pm 13 pc, and for non-zero inclination angles the distance is smaller. Fitting the shape of the bow shock with an analytical model (Wilkin 1996), the effect of the inclination angle on the distance may be estimated. Making the additional assumption that the relative peculiar velocity between the interstellar medium (ISM) and CW Leo is determined entirely by the star space velocity with respect to the local standard of rest (i.e. a stationary ISM), the inclination angle is found to be (-33.3 \pm 0.8) degrees based on the observed proper motion and radial velocity. Using the Wilkin model, our current best estimate of the distance to CW Leo is 123 \pm 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 \pm 150 Lsol (internal error).

We present an overview of the initial results from the Chandra Planetary Nebula Survey (ChanPlaNS), the first systematic (volume-limited) Chandra X-ray Observatory survey of planetary nebulae (PNe) in the solar neighborhood. The first phase of ChanPlaNS targeted 21 mostly high-excitation PNe within ~1.5 kpc of Earth, yielding 4 detections of diffuse X-ray emission and 9 detections of X-ray-luminous point sources at the central stars (CSPNe) of these objects. Combining these results with those obtained from Chandra archival data for all (14) other PNe within ~1.5 kpc that have been observed to date, we find an overall X-ray detection rate of ~70%. Roughly 50% of the PNe observed by Chandra harbor X-ray-luminous CSPNe, while soft, diffuse X-ray emission tracing shocks formed by energetic wind collisions is detected in ~30%; five objects display both diffuse and point-like emission components. The presence of X-ray sources appears correlated with PN density structure, in that molecule-poor, elliptical nebulae are more likely to display X-ray emission (either point-like or diffuse) than molecule-rich, bipolar or Ring-like nebulae. All but one of the X-ray point sources detected at CSPNe display X-ray spectra that are harder than expected from hot (~100 kK) central star photospheres, possibly indicating a high frequency of binary companions to CSPNe. Other potential explanations include self-shocking winds or PN mass fallback. Most PNe detected as diffuse X-ray sources are elliptical nebulae that display a nested shell/halo structure and bright ansae; the diffuse X-ray emission regions are confined within inner, sharp-rimmed shells. All sample PNe that display diffuse X-ray emission have inner shell dynamical ages <~5x10^3 yr, placing firm constraints on the timescale for strong shocks due to wind interactions in PNe.

Far-infrared Herschel/PACS images at 70 and 160 micron of a sample of 78 Galactic evolved stars are used to study the (dust) emission structures, originating from stellar wind-ISM interaction. In addition, two-fluid hydrodynamical simulations of the coupled gas and dust in wind-ISM interactions are used to compare with the observations. Four distinct classes of wind-ISM interaction (i.e. "fermata", "eyes", "irregular", and "rings") are identified and basic parameters affecting the morphology are discussed. We detect bow shocks for ~40% of the sample and detached rings for ~20%. De-projected stand-off distances (R_0) -- defined as the distance between the central star and the nearest point of the interaction region -- of the detected bow shocks ("fermata" and "eyes") are derived from the PACS images and compared to previous results, model predictions and the simulations. All observed bow shocks have stand-off distances smaller than 1 pc. Observed and theoretical stand-off distances are used together to independently derive the local ISM density. Both theoretical (analytical) models and hydrodynamical simulations give stand-off distances for adopted stellar properties that are in good agreement with the measured de-projected stand-off distance of wind-ISM bow shocks. The possible detection of a bow shock -- for the distance limited sample -- appears to be governed by its physical size as set roughly by the stand-off distance. In particular the star's peculiar space velocity and the density of the ISM appear decisive in detecting emission from bow shocks or detached rings. Tentatively, the "eyes" class objects are associated to (visual) binaries, while the "rings" generally appear not to occur for M-type stars, only for C or S-type objects that have experienced a thermal pulse.

In this paper we give a progress report on the Herschel imaging and spectroscopic observations of planetary nebulae that are carried out as part of the MESS guaranteed time key program. We present and discuss imaging and temperature maps of NGC 6720, NGC 650, and NGC 6853, as well as PACS and SPIRE spectroscopy of NGC 7027.

We present Herschel PACS and SPIRE images of the dust shells around the planetary nebulae NGC 650, NGC 6853, and NGC 6720, as well as images showing the dust temperature in their shells. The latter shows a rich structure, which indicates that internal extinction in the UV is important despite the highly evolved status of the nebulae.

The asymptotic giant branch (AGB) stars X Her and TX Psc have been imaged at 70 and 160 microns with the PACS instrument onboard the Herschel satellite, as part of the large MESS (Mass loss of Evolved StarS) Guaranteed Time Key Program. The images reveal an axisymmetric extended structure with its axis oriented along the space motion of the stars. This extended structure is very likely to be shaped by the interaction of the wind ejected by the AGB star with the surrounding interstellar medium (ISM). As predicted by numerical simulations, the detailed structure of the wind-ISM interface depends upon the relative velocity between star+wind and the ISM, which is large for these two stars (108 and 55 km/s for X Her and TX Psc, respectively). In both cases, there is a compact blob upstream whose origin is not fully elucidated, but that could be the signature of some instability in the wind-ISM shock. Deconvolved images of X Her and TX Psc reveal several discrete structures along the outermost filaments, which could be Kelvin-Helmholtz vortices. Finally, TX Psc is surrounded by an almost circular ring (the signature of the termination shock?) that contrasts with the outer, more structured filaments. A similar inner circular structure seems to be present in X Her as well, albeit less clearly.

The detected variety in chemistry and circumstellar shell morphology of the limited sample of Galactic post-AGB stars is so large that there is no consensus yet on how the different objects are linked by evolutionary channels. The evaluation is complicated by the fact that their distances and hence luminosities remain largely unknown. Via cross-correlation of the Spitzer SAGE catalogue with optical catalogues we selected a sample of LMC post-AGB candidates based on their [8]-[24] colour index and estimated luminosity. We determined the fundamental properties of the central stars of 105 of these objects using low-resolution, optical spectra that we obtained at Siding Spring Observatory and SAAO, and constructed a catalogue of 70 high probability and 1337 candidate post-AGB stars that is available at the CDS. The sample forms an ideal testbed for stellar evolution theory predictions of the final phase of low- and intermediate-mass stars, because the distance and hence luminosity and also the current and initial mass of these objects is well constrained. About half of the objects in our sample of post-AGB candidates show a spectral energy distribution (SED) that is indicative of a disc rather than an expanding and cooling AGB remnant. Like in the Galaxy, the disc sources are likely associated with binary evolution. Important side products of this research are catalogues of candidate young stellar objects, candidate supergiants with circumstellar dust, and discarded objects for which a spectrum was obtained. These too are available at the CDS.

We report mid-infrared (MIR) imaging observations of the carbon star U Ant made with the Infrared Camera (IRC) on board AKARI. Subtraction of the artifacts and extended PSF of the central star reveals the detached dust shell around the carbon star at MIR wavelengths (15 and 24{\mu}m) for the first time. The observed radial brightness profiles of the MIR emission are well explained by two shells at 43" and 50" from the central star detected in optical scattered light observations. Combining Herschel/PACS, AKARI/FIS, and AKARI/IRC data, we obtain the infrared spectral energy distribution (SED) of the thermal emission from the detached shell of U Ant in a wide infrared spectral range of 15- 160 {\mu}m. Thermal emission of amorphous carbon grains with a single temperature cannot account for the observed SED from 15 to 160 {\mu}m: it underestimates the emission at 15 {\mu}m. Alternatively, the observed SED is fitted by the model that amorphous carbon grains in the two shells have different temperatures of 60 and 104 K, which allocates most dust mass in the shell at 50". This supports previous suggestion that the 43" shell is gas-rich and the 50" one is dust-rich. We suggest a possibility that the segregation of the gas and dust resulting from the drift motion of submicron-sized dust grains relative to the gas and that the hot dust component associated with the gas-rich shell is composed of very small grains that are strongly coupled with the gas.

MESS (Mass-loss of Evolved StarS) is a Guaranteed Time Key Program that uses the PACS and SPIRE instruments on board the Herschel Space Observatory to observe a representative sample of evolved stars, that include asymptotic giant branch (AGB) and post-AGB stars, planetary nebulae and red supergiants, as well as luminous blue variables, Wolf-Rayet stars and supernova remnants. In total, of order 150 objects are observed in imaging and about 50 objects in spectroscopy. This paper describes the target selection and target list, and the observing strategy. Key science projects are described, and illustrated using results obtained during Herschel's science demonstration phase. Aperture photometry is given for the 70 AGB and post-AGB stars observed up to October 17, 2010, which constitutes the largest single uniform database of far-IR and sub-mm fluxes for late-type stars.

In this article we present the detection of the 69 {\mu}m band of the crystalline olivine forsterite within the MESS key program of Herschel. We determine the temperature of the forsterite grains by fitting the 69 {\mu}m band.

In this paper we will discuss the images of Planetary Nebulae that have recently been obtained with PACS and SPIRE on board the Herschel satellite. This comprises results for NGC 650 (the little Dumbbell nebula), NGC 6853 (the Dumbbell nebula), and NGC 7293 (the Helix nebula).

We present far-infrared and submillimetre spectra of three carbon-rich evolved objects, AFGL 2688, AFGL 618 and NGC 7027. The spectra were obtained with the SPIRE Fourier transform spectrometer on board the Herschel Space Observatory, and cover wavelengths from 195-670 um, a region of the electromagnetic spectrum hitherto difficult to study in detail. The far infrared spectra of these objects are rich and complex, and we measure over 150 lines in each object. Lines due to 18 different species are detected. We determine physical conditions from observations of the rotational lines of several molecules, and present initial large velocity gradient models for AFGL 618. We detect water in AFGL 2688 for the first time, and confirm its presence in AFGL 618 in both ortho and para forms. In addition, we report the detection of the J=1-0 line of CH+ in NGC 7027.

Herschel PACS and SPIRE images have been obtained over a 30'x30' area around the well-known carbon star CW Leo (IRC +10 216). An extended structure is found in an incomplete arc of ~22' diameter, which is cospatial with the termination shock due to interaction with the interstellar medium (ISM) as defined by Sahai & Chronopoulos from ultraviolet GALEX images. Fluxes are derived in the 70, 160, 250, 350, and 550 um bands in the region where the interaction with the ISM takes place, and this can be fitted with a modified black body with a temperature of 25+-3 K. Using the published proper motion and radial velocity for the star, we derive a heliocentric space motion of 25.1 km/s. Using the PACS and SPIRE data and the analytical formula of the bow shock structure, we infer a de-projected standoff distance of the bow shock of R0 = (8.0+-0.3)x10^17 cm. We also derive a relative velocity of the star with respect to the ISM of (106.6+-8.7)/sqrt(n_ISM) km/s, where n_ISM is the number density of the local ISM.

Herschel PACS and SPIRE images have been obtained of NGC 6720 (the Ring Nebula). This is an evolved planetary nebula with a central star that is currently on the cooling track, due to which the outer parts of the nebula are recombining. From the PACS and SPIRE images we conclude that there is a striking resemblance between the dust distribution and the H2 emission, which appears to be observational evidence that H2 forms on grain surfaces. We have developed a photoionization model of the nebula with the Cloudy code which we used to determine the physical conditions of the dust and investigate possible formation scenarios for the H2. We conclude that the most plausible scenario is that the H2 resides in high density knots which were formed after the recombination of the gas started when the central star entered the cooling track. Hydrodynamical instabilities due to the unusually low temperature of the recombining gas are proposed as a mechanism for forming the knots. H2 formation in the knots is expected to be substantial after the central star underwent a strong drop in luminosity about one to two thousand years ago, and may still be ongoing at this moment, depending on the density of the knots and the properties of the grains in the knots.

HD179821 is an enigmatic evolved star that possesses characteristics of both a post-asymptotic giant branch star and a yellow hyper-giant, and there has been no evidence that unambiguously defines its nature. These two hypotheses are products of an indeterminate distance, presumed to be 1 kpc or 6 kpc. We have obtained the two-epoch Hubble Space Telescope WFPC2 data of its circumstellar shell, which shows multiple concentric arcs extending out to about 8 arcsec. We have performed differential proper-motion measurements on distinct structures within the circumstellar shell of this mysterious star in hopes of determining the distance to the object, and thereby distinguishing the nature of this enigmatic stellar source. Upon investigation, rather than azimuthal radially symmetric expansion, we discovered a bulk motion of the circumstellar shell of (2.41+-0.43, 2.97+-0.32) mas/yr. This corresponded to a translational ISM flow of (1.28+-0.95, 7.27+-0.75) mas/yr local to the star. This finding implies that the distance to HD 179821 should be rather small in order for its circumstellar shell to preserve its highly intact spherical structure in the presence of the distorting ISM flow, therefore favoring the proposition that HD 179821 is a post-AGB object.

We carried out a near- to mid-infrared imaging and spectroscopic observations of the patchy areas in the Small Magellanic Cloud using the Infrared Camera on board AKARI. Two 100 arcmin2 areas were imaged in 3.2, 4.1, 7, 11, 15, and 24 um and also spectroscopically observed in the wavelength range continuously from 2.5 to 13.4 um. The spectral resolving power (lambda/Delta lambda) is about 20, 50, and 50 at 3.5, 6.6 and 10.6 um, respectively. Other than the two 100 arcmin2 areas, some patchy areas were imaged and/or spectroscopically observed as well. In this paper, we overview the observations and present a list of near- to mid-infrared photometric results, which lists ~ 12,000 near-infrared and ~ 1,800 mid-infrared bright point sources detected in the observed areas. The 10 sigma limits are 16.50, 16.12, 13.28, 11.26, 9.62, and 8.76 in Vega magnitudes at 3.2, 4.1, 7, 11, 15, and 24 um bands, respectively.

We investigate the dust size and dust shell structure of the bipolar proto-planetary nebula M 1--92 by means of radiative transfer modeling. Our models consists of a disk and bipolar lobes that are surrounded by an AGB shell, each component having different dust characteristics. The upper limit of the grain size $a_\mathrm{max}$ in the lobes is estimated to be $0.5 \mu$m from the polarization value in the bipolar lobe. The $a_\mathrm{max}$ value of the disk is constrained with the disk mass (0.2 $M_{\sun}$), which was estimated from a previous CO emission line observation. We find a good model with $a_\mathrm{max}=1000.0 \mu$m, which provides an approximated disk mass of 0.15 $M_{\sun}$. Even taking into account uncertainties such as the gas-to-dust mass ratio, a significantly larger dust of $a_\mathrm{max}>100.0 \mu$m, comparing to the dust in the lobe, is expected. We also estimated the disk inner radius, the disk outer radius, and the envelope mass to be 30 $R_\star$(=9 AU), 4500 AU, and 4 $M_{\sun}$, respectively, where $v_\mathrm{exp}$ is the expansion velocity. If the dust existing in the lobes in large separations from the central star undergoes little dust processing, the dust sizes preserves the ones in the dust formation. Submicron-sized grains are found in many objects besides M 1--92, suggesting that the size does not depend much on the object properties, such as initial mass of the central star and chemical composition of the stellar system. On the other hand, the grain sizes in the disk do. Evidence of large grains has been reported in many bipolar PPNs, including M 1--92. This result suggests that disks play an important role in grain growth.

The circumstellar dust shells of intermediate initial-mass (about 1 to 8 solar masses) evolved stars are generated by copious mass loss during the asymptotic giant branch phase. The density structure of their circumstellar shell is the direct evidence of mass loss processes, from which we can investigate the nature of mass loss. We used the AKARI Infrared Astronomy Satellite and the Spitzer Space Telescope to obtain the surface brightness maps of an evolved star R Cas at far-infrared wavelengths, since the temperature of dust decreases as the distance from the star increases and one needs to probe dust at lower temperatures, i.e., at longer wavelengths. The observed shell structure and the star's known proper motion suggest that the structure represents the interface regions between the dusty wind and the interstellar medium. The deconvolved structures are fitted with the analytic bow shock structure to determine the inclination angle of the bow shock cone. Our data show that (1) the bow shock cone of 1 - 5 x 10^-5 solar masses (dust mass) is inclined at 68 degrees with respect to the plane of the sky, and (2) the dust temperature in the bow shock cone is raised to more than 20 K by collisional shock interaction in addition to the ambient interstellar radiation field. By comparison between the apex vector of the bow shock and space motion vector of the star we infer that there is a flow of interstellar medium local to R Cas whose flow velocity is at least 55.6 km/s, consistent with an environment conducive to dust heating by shock interactions.

We report here on the discovery of an extended far-infrared shell around the AGB star, R Cassiopeia, made by AKARI and Spitzer. The extended, cold circumstellar shell of R Cas spans nearly 3 arcmin and is probably shaped by interaction with the interstellar medium. This report is one of several studies of well-resolved mass loss histories of AGB stars under AKARI and Spitzer observing programs labeled "Excavating Mass Loss History in Extended Dust Shells of Evolved Stars (MLHES)".

Recent far-infrared mapping of mass-losing stars by the AKARI Infrared Astronomy Satellite and Spitzer Space Telescope have suggested that far-infrared bow shock structures are probably ubiquitous around these mass-losing stars, especially when these stars have high proper motion. Higher spatial resolution data of such far-infrared bow shocks now allow detailed fitting to yield the orientation of the bow shock cone with respect to the heliocentric space motion vector of the central star, using the analytical solution for these bow shocks under the assumption of momentum conservation across a physically thin interface between the stellar winds and interstellar medium (ISM). This fitting analysis of the observed bow shock structure would enable determination of the ambient ISM flow vector, founding a new technique to probe the 3-D ISM dynamics that are local to these interacting systems. In this review, we will demonstrate this new technique for three particular cases, Betelgeuse, R Hydrae, and R Cassiopeiae.

We report the observations of 12 globular clusters with the AKARI/FIS. Our goal is to search for emission from the cold dust within clusters. We detect diffuse emissions toward NGC 6402 and 2808, but the IRAS 100-micron maps show the presence of strong background radiation. They are likely emitted from the galactic cirrus, while we cannot rule out the possible association of a bump of emission with the cluster in the case of NGC 6402. We also detect 28 point-like sources mainly in the WIDE-S images (90 micron). At least several of them are not associated with the clusters but background galaxies based on some external catalogs. We present the SEDs by combining the near-and-mid infrared data obtained with the IRC if possible. The SEDs suggest that most of the point sources are background galaxies. We find one candidate of the intracluster dust which has no mid-infrared counterpart unlike the other point-like sources, although some features such as its point-like appearance should be explained before we conclude its intracluster origin. For most of the other clusters, we have confirmed the lack of the intracluster dust. We evaluate upper limits of the intracluster dust mass to be between 1.0E-05 and 1.0E-03 solar mass depending on the dust temperature. The lifetime of the intracluster dust inferred from the upper limits is shorter than 5 Myr (T=70K) or 50 Myr (35K). Such short lifetime indicates some mechanism(s) are at work to remove the intracluster dust. We also discuss its impact on the chemical evolution of globular clusters.

Evolved mass-losing stars such as Mira enrich the interstellar medium (ISM) significantly by their dust-rich molecular wind. When these stars move fast enough relative to the ISM, the interaction between the wind and ISM generates the structure known as the astropause (a stellar analog of the heliopause), which is a cometary stellar wind cavity bounded by the contact discontinuity surface between the wind and ISM. Far-infrared observations of Mira spatially resolve the structure of its astropause for the first time, distinguishing the contact surface between Mira's wind and the ISM and the termination shock due to Mira's wind colliding with the ISM. The physical size of the astropause and the estimated speed of the termination shock suggest the age of the astropause to be about 40,000 yr, confirming a theoretical prediction of the shock re-establishment time after Mira has entered the Local Bubble.

We present 10' x 50' scan maps around an M supergiant Alpha Ori at 65, 90, 140 and 160 microns obtained with the AKARI Infrared Astronomy Satellite. Higher spatial resolution data with the exact analytic solution permit us to fit the de-projected shape of the stellar wind bow shock around Alpha Ori to have the stand-off distance of 4.8', position angle of 55 degrees and inclination angle of 56 degrees. The shape of the bow shock suggests that the peculiar velocity of Alpha Ori with respect to the local medium is v_* = 40 (n_H)^(-1/2), where n_H is the hydrogen nucleus density at Alpha Ori. We find that the local medium is of n_H = 1.5 to 1.9 cm^(-3) and the velocity of the local flow is at 11 km s^(-1) by using the most recent astrometric solutions for Alpha Ori under the assumption that the local medium is moving away from the Orion OB 1 association. AKARI images may also reveal a vortex ring due to instabilities on the surface of the bow shock as demonstrated by numerical models. This research exemplifies the potential of AKARI All-Sky data as well as follow-up observations with Herschel Space Telescope and Stratospheric Observatory for Infrared Astronomy for this avenue of research in revealing the nature of interaction between the stellar wind and interstellar medium.

The results from a Hubble Space Telescope (HST) snapshot survey of post-AGB objects are shown. The aim of the survey is to complement existing HST images of PPN and to connect various types of nebulosities with physical and chemical properties of their central stars. Nebulosities are detected in 15 of 33 sources. Images and photometric and geometric measurements are presented. For sources with nebulosities we see a morphological bifurcation into two groups, DUPLEX and SOLE, as previous studies have found. We find further support to the previous results suggesting that this dichotomy is caused by a difference in optical thickness of the dust shell. The remaining 18 sources are classified as stellar post-AGB objects, because our observations indicate a lack of nebulosity. We show that some stellar sources may in fact be DUPLEX or SOLE based on their infrared colors. The cause of the differences among the groups are investigated. We discuss some evidence suggesting that high progenitor-mass AGB stars tend to become DUPLEX post-AGB objects. Intermediate progenitor-mass AGB stars tend to be SOLE post-AGB objects. Most of the stellar sources probably have low mass progenitors and do not seem to develop nebulosities during the post-AGB phase and therefore do not become planetary nebulae.

context: The spectrum of IRAS 16594-4656 shows shock excited H_2 emission and collisionally excited emission lines such as[O I],[C I],and [Fe II]. aim: The goal is to determine the location of the H_2 and [Fe II] shock emission, to determine the shock velocities,and constrain the physical properties in the shock. methods: High resolution spectra of the H_2 1-0 S(1),H_2 2-1 S(1), [Fe II], and Pa$\beta$ emission lines were obtained with the near infrared spectrograph Phoenix on Gemini South. results: The position-velocity diagrams of H_2 1-0 S(1), H_2 2-1 S(1), and [Fe II] are presented. The H_2 and [Fe II] emission is spatially extended. The collisionally excited [O I] and [C I] optical emission lines have a similar double peaked profile compared to the extracted H_2 profile and appear to be produced in the same shock. They all indicate an expansion velocity of ~8 km/s and the presence of a neutral, very high density region with $n_{\rm e}$ about 3 x 10^6 to 5 x10^7 cm$^{-3}$. The [Fe II] emission however is single peaked. It has a gaussian FWHM of 30 km/s and a total width of 62 km/s at 1% of the peak. The Pa$\beta$ profile is even wider with a gaussian FWHM of 48 km/s and a total width of 75 km/s at 1% of the peak. conclusions: The H$_2$ emission is excited in a slow 5 to 20 km/s shock into dense material at the edge of the lobes, caused by the interaction of the AGB ejecta and the post-AGB wind. The 3D representation of the H_2 data shows a hollow structure with less H_2 emission in the equatorial region. The [Fe II] emission is not present in the lobes, but originates close to the central star in fast shocks in the post-AGB wind or in a disk. The Pa$\beta$ emission also appears to originate close to the star.

The structure of the dusty circumstellar envelopes (CSEs) of proto-planetary nebulae (PPNs) reveals the mass-loss history of these sources and how such histories may differ for elliptical (SOLE) and bipolar (DUPLEX) PPNs. To study the PPN structures via dust-scattered linearly polarized starlight, we have compiled the imaging-polarimetric data for all 18 evolved stars that have been obtained to date with NICMOS on-board the Hubble Space Telescope (HST). This alternative imaging technique provides a unique way to probe the distribution of dust grains that scatter light around evolved stars. The new perspective gained from the imaging-polarimetric data has revealed several new aspects to the structures of PPNs. Point-symmetry is a prevalent imaging-polarimetric characteristic resulting from the azimuthal density gradient in the CSEs. Among these point-symmetric nebulae, three detailed morphological types can be differentiated by their polarized intensity, I_pol, and polarization strength, P. While the azimuthal density gradient is reversed above and below the equatorial plane in optically thicker bipolar nebulae, there is no gradient reversal in optically thinner elliptical nebulae. The equatorial plane of the system defined by the integrated angle of polarization is not necessarily orthogonal to the axis of the apparent bipolar structure in the total intensity data.

Color-magnitude diagrams (CMDs) are presented for the Spitzer SAGE (Surveying the Agents of a Galaxy's Evolution) survey of the Large Magellanic Cloud (LMC). IRAC and MIPS 24 um epoch one data are presented. These data represent the deepest, widest mid-infrared CMDs of their kind ever produced in the LMC. Combined with the 2MASS survey, the diagrams are used to delineate the evolved stellar populations in the Large Magellanic Cloud as well as Galactic foreground and extragalactic background populations. Some 32000 evolved stars brighter than the tip of the red giant branch are identified. Of these, approximately 17500 are classified as oxygen-rich, 7000 carbon-rich, and another 1200 as ``extreme'' asymptotic giant branch (AGB) stars. Brighter members of the latter group have been called ``obscured'' AGB stars in the literature owing to their dusty circumstellar envelopes. A large number (1200) of luminous oxygen--rich AGB stars/M supergiants are also identified. Finally, there is strong evidence from the 24 um MIPS channel that previously unexplored, lower luminosity oxygen-rich AGB stars contribute significantly to the mass loss budget of the LMC (1200 such sources are identified).

New Spitzer imaging observations have revealed the structure around the Mira variable star R Hya to be a one-sided parabolic arc 100 arcsec to the West stretching from North to South. We successfully model R Hya and its surroundings in terms of an interaction of the stellar wind from an asymptotic giant branch (AGB) star with the interstellar medium (ISM) the star moves through. Our three-dimensional hydrodynamic simulation reproduces the structure as a bow shock into the oncoming ISM. We propose this as another explanation of detached shells around such stars which should be considered alongside current theories of internal origin. The simulation predicts the existence of a tail of ram-pressure-stripped AGB material stretching downstream. Indications for such a tail behind R Hya are seen in IRAS maps.

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. \

We present the far-infrared (IR) maps of a bipolar planetary nebula (PN), NGC 650, at 24, 70, and 160 micron taken with the Multiband Imaging Photometer for Spitzer (MIPS) on-board the Spitzer Space Telescope. While the two-peak emission structure seen in all MIPS bands suggests the presence of a near edge-on dusty torus, the distinct emission structure between the 24 micron map and the 70/160 micron maps indicates the presence of two distinct emission components in the central torus. Based on the spatial correlation of these two far-IR emission components with respect to various optical line emission, we conclude that the 24 micron emission is largely due to the [O IV] line at 25.9 micron arising from highly ionized regions behind the ionization front, whereas the 70 and 160 micron emission is due to dust continuum arising from low-temperature dust in the remnant asymptotic giant branch (AGB) wind shell. The far-IR nebula structure also suggests that the enhancement of mass loss at the end of the AGB phase has occurred isotropically, but has ensued only in the equatorial directions while ceasing in the polar directions. The present data also show evidence for the prolate spheroidal distribution of matter in this bipolar PN. The AGB mass loss history reconstructed in this PN is thus consistent with what has been previously proposed based on the past optical and mid-IR imaging surveys of the post-AGB shells.

We are performing a uniform and unbiased, ~7x7 degrees imaging survey of the Large Magellanic Cloud (LMC), using the IRAC and MIPS instruments on board the Spitzer Space Telescope in order to survey the agents of a galaxy's evolution (SAGE), the interstellar medium (ISM) and stars in the LMC. The detection of diffuse ISM with column densities >1.2x10^21 H cm^-2 permits detailed studies of dust processes in the ISM. SAGE's point source sensitivity enables a complete census of newly formed stars with masses >3 solar masses that will determine the current star formation rate in the LMC. SAGE's detection of evolved stars with mass loss rates >1x10^-8 solar masses per year will quantify the rate at which evolved stars inject mass into the ISM of the LMC. The observing strategy includes two epochs in 2005, separated by three months, that both mitigate instrumental artifacts and constrain source variability. The SAGE data are non-proprietary. The data processing includes IRAC and MIPS pipelines and a database for mining the point source catalogs, which will be released to the community in support of Spitzer proposal cycles 4 and 5. We present initial results on the epoch 1 data with a special focus on the N79 and N83 region. The SAGE epoch 1 point source catalog has ~4 million sources. The point source counts are highest for the IRAC 3.6 microns band and decrease dramatically towards longer wavelengths consistent with the fact that stars dominate the point source catalogs and that the dusty objects, e.g. young stellar objects and dusty evolved stars that detected at the longer wavelengths, are rare in comparison. We outline a strategy for identifying foreground MW stars, that may comprise as much as 18% of the source list, and background galaxies, that may comprise ~12% of the source list.

We present the results of proper-motion measurements of the dust shell structure in the Egg Nebula (AFGL 2688, CRL 2688, V1610 Cyg), based on the archived two-epoch imaging-polarimetry data in the 2 micron band taken with NICMOS on-board the Hubble Space Telescope. We measured the amount of motions of local structures in the nebula by determining their relative shifts over an interval of 5.5 years. The dynamical age of the nebula is found to be roughly 350 years based on the overall motion of the nebula that exhibits a Hubble-law-esque linear relation between the measured proper motion and the projected radial distance from the origin of the expansion. By adopting the de-projected velocity of 45 km/s at the tips of the bipolar lobes, our proper-motion measurements indicate that the distance to the Egg Nebula is about 420 pc and that the lobes are inclined at 7.7 degrees with respect to the plane of the sky. The refined distance estimate yields the luminosity of the central star of 3.3 x 10^3 L_sun, the total shell mass of 1.2 M_sun, and the mass loss rate (the upper limit) of 3.6 x 10^(-3) M_sun/yr. Assuming 0.6 M_sun central post-AGB stellar mass, the initial mass of the Egg is 1.8 M_sun. Upon analysis, we also discovered that the central star of the Egg Nebula has proper motion of its own at the rate of (14, -10) mas/yr.

Using NICMOS on HST, we have performed imaging polarimetry of proto-planetary nebulae. Our objective is to study the structure of optically thin circumstellar shells of post-asymptotic giant branch stars by separating dust-scattered, linearly polarized star light from unpolarized direct star light. This unique technique allows us to probe faint reflection nebulae around the bright central star, which can be buried under the point-spread-function of the central star in conventional imaging. Our observations and archival search have yielded polarimetric images for five sources: IRAS 07134+1005 (HD 56126), IRAS 06530-0213, IRAS 04296+3429, IRAS (Z)02229+6208, and IRAS 16594-4656. These images have revealed the circumstellar dust distribution in an unprecedented detail via polarized intensity maps, providing a basis to understand the 3-D structure of these dust shells. We have observationally confirmed the presence of the inner cavity caused by the cessation of AGB mass loss and the internal shell structures which is strongly tied to the progenitor star's mass loss history on the AGB. We have also found that equatorial enhancement in these circumstellar shells comes with various degrees of contrast, suggesting a range of optical depths in these optically thin shells. Our data support the interpretation that the dichotomy of PPN morphologies is due primarily to differences in optical depth and secondary to the inclination effect. The polarization maps reveal a range of inclination angles for these optically thin reflection nebulae, dispelling the notion that elliptical nebulae are pole-on bipolar nebulae.

We present millimeter interferometry images of the CO J=1-0 line emission arising in the circumstellar envelope of HD 56126 (a.k.a. IRAS 07134+1005), which is one of the best studied 21-micron proto-planetary nebulae (PPNs). The CO emission extends from 1.2" to 7" in radius from the central star and appears consistent with a simple expanding envelope, as expected for a post-AGB star. We quantitatively model the molecular envelope using a radiative transfer code that we have modified for detached shells. Our best fit model reveals that two sequential winds created the circumstellar envelope of HD 56126: an AGB wind that lasted 6500 years with a mass-loss rate of 5.1x10^{-6} M_odot yr^{-1} and a more intense superwind that lasted 840 years with a mass-loss rate of 3x10^{-5} M_odot yr^{-1} and that ended the star's life on the AGB 1240 years ago. The total mass of this envelope is 0.059 M_odot which indicates a lower limit progenitor mass for the system of 0.66 M_odot, quite reasonable for this low-metallicity star which probably resides in the thick disk of the Galaxy. Comparison with images of the dust emission reveal a similar structure with the gas in the inner regions. Using 2-DUST, we model the dust emission of this source so that the model is consistent with the CO emission model and find an average gas-to-dust mass ratio of 75.

We present the preliminary results from our HST/NICMOS imaging polarimetry on optically thin PPNs. The data show the dust distribution in an unprecedented detail. The structure of these faint PPNs is revealed via polarization properties of the dust-scattered light without much interference from the bright central star, proving the unique advantages of the technique over simple imaging.

Past imaging observations of a post-AGB star, HD 161796, have probed different parts of its circumstellar shell which correspond to different epochs of the star's AGB mass loss history. While the overall structure of the shell can be described by an axisymmetric model consisting of three layers of characteristic structure, the mass distribution in the shell appears to rotate its axis of symmetry in a continuous manner. It is thus interesting to observe the halo region of the shell using from far-IR to sub-mm wavelengths and increase the ``time resolution'' during this critical epoch of the star's mass loss history.

We present JHK' photometric measurements of 78 objects mostly consisting of proto-planetary nebula candidates. Photometric magnitudes are determined by means of imaging and aperture photometry. Unlike the observations with a photometer with a fixed-sized beam, the method of imaging photometry permits accurate derivation of photometric values because the target sources can be correctly identified and confusion with neighboring sources can be easily avoided. Of the 78 sources observed, we report 10 cases in which the source seems to have been misidentified or confused by nearby bright sources. We also present nearly two dozen cases in which the source seems to have indicated a variability which prompts a follow-up monitoring. There are also a few sources that show previously unreported extendedness. In addition, we present H band finding charts of the target sources.

We have developed a general purpose dust radiative transfer code for an axisymmetric system, 2-Dust, motivated by the recent increasing availability of high-resolution images of circumstellar dust shells at various wavelengths. This code solves the equation of radiative transfer following the principle of long characteristic in a 2-D polar grid while considering a 3-D radiation field at each grid point. A solution is sought through an iterative scheme in which self-consistency of the solution is achieved by requiring a global luminosity constancy throughout the shell. The dust opacities are calculated through Mie theory from the given size distribution and optical properties of the dust grains. The main focus of the code is to obtain insights on (1) the global energetics of dust grains in the shell (2) the 2-D projected morphologies that are strongly dependent on the mixed effects of the axisymmetric dust distribution and inclination angle of the shell. Here, test models are presented with discussion of the results. The code can be supplied with a user-defined density distribution function, and thus, is applicable to a variety of dusty astronomical objects possessing the axisymmetric geometry.

We use detailed radiative transfer models to investigate the differences between the Star-Obvious Low-level-Elongated proto-planetary nebulae (SOLE PPNs) and DUst-Prominent Longitudinally-EXtended proto-planetary nebulae (DUPLEX PPNs) which are two subclasses of PPNs suggested by Ueta, Meixner, & Bobrowsky (2000). We select one SOLE PPN, HD 161796, and one DUPLEX PPN, IRAS 17150$-$3224, both of which are well studied and representative of their PPN classes. Using an axisymmetric dust shell radiative transfer code, we model these two sources in detail and constrain their mass-loss histories, inclination angles and dust composition. The physical parameters derived for HD 161796 and IRAS 17150$-$3224 demonstrate that they are physically quite different and that their observed differences cannot be attributed to inclination angle effects. If these calculations reflect a more general truth about SOLE vs. DUPLEX PPNs, then these two subclasses of PPNs are physically distinct with the SOLE PPNs derived from low mass progenitors and DUPLEX PPNs derived from high mass progenitors.

We present high-sensitivity near-IR images of a carbon-rich proto-planetary nebula, AFGL 618, obtained with the Subaru Telescope. These images have revealed ``bullets'' and ``horns'' extending farther out from the edges of the previously known bipolar lobes. The spatial coincidence between these near-IR micro-structures and the optical collimated outflow structure, together with the detection of shock-excited, forbidden IR lines of atomic species, strongly suggests that these bullets and horns represent the locations from which [\ion{Fe}{2}] IR lines arise. We have also discovered CO clumps moving at $> 200$ km s$^{-1}$ at the positions of the near-IR bullets by re-analyzing the existing $^{12}$CO $J=1-0$ interferometry data. These findings indicate that the near-IR micro-structures represent the positions of shocked surfaces at which fast-moving molecular clumps interface with the ambient circumstellar shell.

We report sub-arcsecond imaging of extended mid-infrared emission from a proto-planetary nebula (PPN), \iras 22272+5435, performed at the MMT observatory with its newly upgraded 6.5 m aperture telescope and at the Keck observatory. The mid-infrared emission structure is resolved into two emission peaks separated by $0\arcsec.5 - 0\arcsec.6$ in the MMT 11.7 $\um$ image and in the Keck 7.9, 9.7, and 12.5 $\um$ images, corroborating the predictions based on previous multi-wavelength morphological studies and radiative transfer calculations. The resolved images show that the PPN dust shell has a toroidal structure with the $0\arcsec.5$ inner radius. In addition, an unresolved mid-IR excess appears at the nebula center. Radiative transfer model calculations suggest that the highly equatorially-enhanced ($\rho_{\rm eq}/\rho_{\rm pole} = 9$) structure of the PPN shell was generated by an axisymmetric superwind with ${\dot M}_{\rm sw} = 4 \times 10^{-6} M_{\odot}$ yr$^{-1}$, which was preceded by a spherical asymptotic giant branch (AGB) wind with ${\dot M}_{\rm AGB} = 8 \times 10^{-7} M_{\odot}$ yr$^{-1}$. These model calculations also indicate that the superwind shell contains larger dust grains than the AGB wind shell. The unresolved mid-infrared excess may have been produced by a post-AGB mass loss at a rate of $2 - 6 \times 10^{-7} M_{\odot}$ yr$^{-1}$. While the central star left the AGB about 380 years ago after the termination of the superwind, the star seems to have been experiencing an ambient post-AGB mass loss with a sudden, increased mass ejection about 10 years ago.

We will review our mid-infrared and HST imaging surveys of the circumstellar dust shells of proto-planetary nebulae. While optical imaging indirectly probes the dust distribution via dust-scattered starlight, mid-IR imaging directly maps the distribution of warm dust grains. Both imaging surveys revealed preferencially axisymmetric nature of PPN dust shells, suggesting that axisymmetry in planetary nebulae sets in by the end of the asymptotic giant branch phase, most likely by axisymmetric superwind mass loss. Moreover, both surveys yielded two morphological classes which have one-to-one correspondence between the two surveys, indicating that the optical depth of circumstellar dust shells plays an equally important role as the inclination angle in determining the morphology of the PPN shells.

We investigate spatial distributions of specific dust components in the circumstellar shells of a proto-planetary nebula candidate, HD 179821, and a planetary nebula, BD$+30^{\circ}$3639, by means of spectral imaging. With high-resolution ground-based images and ISO spectra in the mid-infrared, we can derive ``dust feature only'' maps by subtracting synthesized continuum maps from the observed images at the feature wavelength. Such spatially detailed information will help to develop models for the evolution of dust grains around evolved stars.

We report detection of an extended mid-infrared emission from \iras 18576+0341 (AFGL 2298). The object shows a dusty circumstellar shell that has diameter of $\age 7\arcsec$ at 10.3 and 18.0 $\um$. The dust nebula shows two emission peaks concentrically elongated and symmetrically oriented on the opposite sides of the third, central peak, which appears to be the central star of the system. The observed mid-infrared morphology indicates that the circumstellar dust shell has an equatorially-enhanced material distribution, which is a common signature of stellar objects that have experienced mass loss. Radiative transfer model calculations suggest that the central star is an extremely bright ($\lstar = 10^{6.4} \lsun$) star at a distance of about 10 kpc: this object is best described as a new luminous blue variable candidate. The circumstellar dust shell seems to have been generated by an equatorially-enhanced mass loss process with ${\dot M} \ge 6.8 \times 10^{-6} \msun$ yr$^{-1}$ and ${\dot M}_{\rm pole}/{\dot M}_{\rm eq} \sim 0.5$.

We summarize results of our {\it HST} imaging survey of proto-planetary nebula (PPN) candidates, in which we discovered two types of axisymmetric reflection nebulosities. The Star-Obvious Low-level-Elongated (SOLE) nebulae show a bright central star embedded in a smooth, faint extended nebulosity, whereas the DUst-Prominent Longitudinally-EXtended (DUPLEX) nebulae have remarkable bipolar structure with a completely or partially obscured central star. The intrinsic axisymmetry of the PPN reflection nebulosities demonstrates that the axisymmetry often found in planetary nebulae predates the PPN phase. We suggest that the major cause for the apparent dual morphology is the optical depth difference in the circumstellar dust shell rather than the inclination angle effect.

I-band images of the OH maser source OH44.8-2.3 are presented. Having a well-determined distance, this object is shown to be very young and physically small, with an asymmetric compact bipolar nebula already present. In the morphological terminology of Meixner et al. (2000) and Ueta et al. (2000), this object would be classified as a DUPLEX nebula because of its high optical depth and pole-to-equator density contrast.

We report the results from an optical imaging survey of proto-planetary nebula candidates using the HST. We exploited the high resolving power and wide dynamic range of HST and detected nebulosities in 21 of 27 sources. All detected reflection nebulosities show elongation, and the nebula morphology bifurcates depending on the degree of the central star obscuration. The Star-Obvious Low-level-Elongated (SOLE) nebulae show a bright central star embedded in a faint, extended nebulosity, whereas the DUst-Prominent Longitudinally-EXtended (DUPLEX) nebulae have remarkable bipolar structure with a completely or partially obscured central star. The intrinsic axisymmetry of these proto-planetary nebula reflection nebulosities demonstrates that the axisymmetry frequently found in planetary nebulae predates the proto-planetary nebula phase, confirming previous independent results. We suggest that axisymmetry in proto-planetary nebulae is created by an equatorially enhanced superwind at the end of the asymptotic giant branch phase. We discuss that the apparent morphological dichotomy is caused by a difference in the optical thickness of the circumstellar dust/gas shell with a differing equator-to-pole density contrast. Moreover, we show that SOLE and DUPLEX nebulae are physically distinct types of proto-planetary nebulae, with a suggestion that higher mass progenitor AGB stars are more likely to become DUPLEX proto-planetary nebulae.

We present the data from a mid-infrared imaging survey of 66 proto-planetary nebula candidates using two mid-IR cameras (MIRAC2 and Berkcam) at the NASA Infrared Telescope Facility and the United Kingdom Infrared Telescope. The goal of this survey is to determine the size, flux, and morphology of the mid-IR emission regions, which sample the inner regions of the circumstellar dust shells of proto-planetary nebulae. We imaged these proto-planetary nebulae with narrow-band filters ($\Delta\lambda / \lambda \sim 10%$) at wavelengths of notable dust features. With our typical angular resolution of 1\arcsec, we resolve 17 sources, find 48 objects unresolved, and do not detect 1 source. For several sources, we checked optical and infrared associations and positions of the sources. In table format, we list the size and flux measurements for all the detected objects and show figures of all the resolved sources. Images for all the detected objects are available on line in FITS format from the Astronomy Digital Image Library at the National Center for Supercomputing Application. The proto-planetary nebula candidate sample includes, in addition to the predominant proto-planetary nebulae, extreme asymptotic giant branch stars, young planetary nebulae, a supergiant, and a luminous blue variable. We find that dust shells which are cooler ($\rm T \sim 150$ K) and brighter in the infrared are more easily resolved. Eleven of the seventeen resolved sources are extended and fall into one of two types of mid-IR morphological classes: core/elliptical or toroidal. Core/elliptical structures show unresolved cores with lower surface brightness elliptical nebulae. Toroidal structures show limb-brightened peaks suggesting equatorial density enhancements. We argue that core/ellipticals have denser dust shells than toroidals.