In a continuing effort to investigate the role of magnetic fields in evolved low and intermediate mass stars (principally regarding the shaping of their envelopes), we present new ALMA high resolution polarization data obtained for the nebula OH 231.8+4.2. We found that the polarized emission likely arises from aligned grains in the presence of magnetic fields rather than radiative alignment and self scattering. The ALMA data show well organized electric field orientations in most of the nebula and the inferred magnetic field vectors (rotated by 90 degrees) trace an hourglass morphology centred on the central system of the nebula. One region in the southern part of OH 231.8+4.2 shows a less organized distribution probably due to the shocked environment. These findings, in conjunction with earlier investigations (maser studies and dust emission analysis at other scales and wavelengths) suggest an overall magnetic hourglass located inside a toroidal field. We propose the idea that the magnetic field structure is closely related to the architecture of a magnetic tower and that the outflows were therefore magnetically launched. While the current dynamical effect of the fields might be weak in the equatorial plane principally due to the evolution of the envelope, it would still be affecting the outflows. In that regard, the measurement of the magnetic field at the stellar surface, which is still missing, combined with a full MHD treatment are required to better understand and constrain the events occurring in OH 231.8+4.2.

We have mapped the pre-planetary nebula IRAS 17150-3224 in 350 GHz continuum and CO J=3-2 line at an angular resolution of about 0".09 using the Atacama Large Millimeter/submillimeter Array. A continuum source is detected at the center of the nebula, elongated along the equatorial plane, likely tracing a dusty torus around the central source. Continuum emission is also detected on both sides of the central continuum source in the equatorial plane, probably resulting from interactions of collimated fast winds with envelope material in the equator. CO emission is detected along the optical lobe. Although the optical lobe appears as bipolar, CO map shows that it is actually a quadrupolar outflow consisting of two overlapping bipolar outflows. Two additional younger bipolar outflows are also detected in CO, one at lower latitude and the other along the equatorial plane. In the CO position-velocity maps, blueshifted absorption stripes are detected in the outflow emissions, due to absorption by a series of shells produced by a series of asymptotic giant branch (AGB) wind ejections. By modeling the morphology and kinematics of the AGB wind and outflows, we find that the AGB wind could end about 1300 yr ago, the quadrupolar outflow was ejected about 350 yr ago, and the two additional bipolar outflows were ejected about 280 and 200 ago, respectively. The outflows could be produced either by bullets coming from an explosion, or by a precessing collimated fast wind with a time-dependent ejection velocity.

In this work we aimed to describe the three-dimensional morphology and kinematics of the molecular gas of the water-fountain nebula IRAS 16342-3814. In order to do this, we retrieved data from the ALMA archive to analyse it using a simple spatio-kinematical model. We used the software SHAPE to construct a three-dimensional spatio-kinematical model of the molecular gas in IRAS 16342-3814. By reproducing the intensity distribution and position-velocity diagram of the CO emission from the ALMA observations we derived the morphology and velocity field of the gas. We used CO(1-0) data to support the physical interpretation of the model. A spatio-kinematical model that includes a high-velocity collimated outflow embedded within material expanding at relatively lower velocity reproduces the images and position-velocity diagrams from the observations. The high-velocity collimated outflow exhibits deceleration across its length, while the velocity of the surrounding component increases with distance. The morphology of the emitting region; the velocity field and the mass of the gas as function of velocity are in excellent agreement with the properties predicted for a molecular outflow driven by a jet. The timescale of the molecular outflow is estimated to be ~70-100 years. An oscillating pattern was found associated to the high-velocity collimated outflow. The oscillation period of the pattern is T~60-90 years and its opening angle is ~2 degrees. The CO (3-2) emission in IRAS 16342-3814 is interpreted in terms of a jet-driven molecular outflow expanding along an elongated region. The position-velocity diagram and the mass spectrum reveal a feature due to entrained material that is associated to the driving jet. It is likely that the jet in those objects has already disappeared since it is expected to last only for a couple of hundred years.

The carbon star V Hydrae (V Hya) provides new insight into the nature of the launching mechanism of jet-like outflows that are believed to be the cause of the poorly understood transition phase of AGB stars into aspherical planetary nebulae. V Hya has been shown to periodically eject collimated gas blobs at high velocities ("bullets"). By analyzing data from HST/STIS 2-D spectra, obtained at six epochs spaced over a decade that show 4 successively ejected bullets with a spacing of ~8.5 years, we have created kinematic models of the dynamical evolution of a specific bullet (#1) for the first three observed epochs (2002, 2003, 2004) using a 3D spatio-kinematic code, SHAPE. Using these models, we fit the observed morphology, line-of-sight velocity, proper motion and intensity for the extended, gaseous bullet as a function of time over a period of 2 years, in order to constrain its 3D movement and the evolution of its physical properties over this period. Our results suggest that although bullet #1's motion is predominantly ballistic, there are small but significant changes in the position angle and inclination angle of the long (symmetry) axis of the bullet that tilt it progressively towards the symmetry axis of the bipolar molecular nebula around V Hya. In contrast, bullet #3 shows strong acceleration soon after ejection. We discuss the possibilities that bullet acceleration is caused by either a non-radial magnetic field and/or by hydrodynamic interaction with the ambient gas through which the bullet is traveling.

Understanding strong binary interactions is of wide astrophysical importance, and the deaths of most stars in the Universe that evolve in a Hubble time could be fundamentally affected by such interactions. These stars end their lives, evolving from Asymptotic Giant Branch stars with extensive mass-loss into planetary nebulae with a spectacular array of morphologies. Binarity, and the associated formation of accretion disks (that drive collimated, fast jets) during the very late AGB or early post-AGB phase is believed to produce this dramatic morphological transformation. But the evidence for binarity and accretion during the AGB phase has been hard to obtain due to observational limitations. However, recent observations at UV and X-ray wavelengths have broken thru the observational barrier -- our studies using GALEX reveal a candidate population of AGB stars, generally with strongly-variable far-ultraviolet (FUV) emission (fuvAGB stars), and our follow-up studies with XMM-Newton, Chandra, and HST of a few key objects supports our hypothesis that these objects have companions that are actively accreting material from the primary. The most prominent fuvAGB star has been detected with the VLA, showing the presence of variable non-thermal emission. The ngVLA, with its unprecedented sensitivity, is needed to survey a statistical sample of fuvAGB stars over the ~3-90 Ghz range to search for and characterize the nature of the radio emission from fuvAGB stars and test our binarity+accretion hypothesis. Such a survey will distinguish between binarity/accretion-related radio emission that is expected to have both thermal and non-thermal components and display significantly time-variable on short time-scales (minutes to weeks), and single-star chromospheric emission from the primary that is expected to be thermal, possibly with time-variability, but only on long time-scales (many months to a year).

Planetary nebulae (PNe) represent the near endpoints of evolution for stars of initial mass $\sim$1-8 $M_\odot$, wherein the envelope of an asymptotic giant branch (AGB) star becomes photodissociated and ionized by high-energy radiation from a newly emerging white dwarf that was the progenitor star's core. It is increasingly evident that most PNe are descended from binary systems. PNe hence provide unique insight into a diverse range of astrophysical phenomena, including the influence of companion stars on the late stages of stellar evolution; stellar wind interactions and shocks; the physics and chemistry of photoionized plasmas and photon-dominated regions (PDRs); and enrichment of the ISM in the products of intermediate-mass stellar nucleosynthesis. We describe specific examples of the potential impact of the ngVLA in each of these areas.

Binarity is believed to dramatically affect the history and geometry of mass loss in AGB and post-AGB stars, but observational evidence of binarity is sorely lacking. As part of a project to search for hot binary companions to cool AGB stars using the GALEX archive, we discovered a late-M star, Y Gem, to be a source of strong and variable UV and X-ray emission. Here we report UV spectroscopic observations of Y Gem obtained with the Hubble Space Telescope that show strong flickering in the UV continuum on time-scales of <~20 s, characteristic of an active accretion disk. Several UV lines with P-Cygni-type profiles from species such as Si IV and C IV are also observed, with emission and absorption features that are red- and blue- shifted by velocities of ~500 km/s from the systemic velocity. Our model for these (and previous) observations is that material from the primary star is gravitationally captured by a companion, producing a hot accretion disk. The latter powers a fast outflow that produces blue-shifted features due to absorption of UV continuum emitted by the disk, whereas the red-shifted emission features arise in heated infalling material from the primary. The outflow velocities support a previous inference by Sahai et al. (2015) that Y Gem's companion is a low-mass main-sequence star. Blackbody fitting of the UV continuum implies an accretion luminosity of about 13 Lsun and thus a mass-accretion rate >5e-7 Msun/yr; we infer that Roche lobe overflow is the most likely binary accretion mode for Y Gem.

We present interferometric and single-dish molecular line observations of the interstellar bullet-outflow source IRAS 05506+2414, whose wide-angle bullet spray is similar to the Orion BN/KL explosive outflow and likely arises from an entirely different mechanism than the classical accretion-disk-driven bipolar flows in young stellar objects. The bullet-outflow source is associated with a large pseudo-disk and three molecular outflows -- a high-velocity outflow (HVO), a medium-velocity outflow (MVO), and a slow, extended outflow (SEO). The size (mass) of the pseudo-disk is 10,350 AU x 6,400 AU (0.64-0.17 Msun); from a model-fit assuming infall and rotation we derive a central stellar mass of 8--19 Msun. The HVO (MVO) has an angular size ~5180 (~3330) AU, and a projected outflow velocity of ~140 km/s (~30 km/s). The SEO size (outflow speed) is ~0.9 pc (~6 km/s). The HVO's axis is aligned with (orthogonal to) that of the SEO (pseudo-disk). The velocity structure of the MVO is unresolved. The scalar momenta in the HVO and SEO are very similar, suggesting that the SEO has resulted from the HVO interacting with ambient cloud material. The bullet spray shares a common axis with the pseudo-disk, and has an age comparable to that of MVO (few hundred years), suggesting that these three structures are intimately linked together. We discuss several models for the outflows in IRAS 05506+2414 (including dynamical decay of a stellar cluster, chance encounter of a runaway star with a dense cloud, and close passage of two protostars), and conclude that 2nd-epoch imaging to derive proper motions of the bullets and nearby stars can help to discriminate between them.

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.

Our Cycle 0 ALMA observations confirmed that the Boomerang Nebula is the coldest known object in the Universe, with a massive high-speed outflow that has cooled significantly below the cosmic background temperature. Our new CO 1-0 data reveal heretofore unseen distant regions of this ultra-cold outflow, out to $\gtrsim120,000$ AU. We find that in the ultra-cold outflow, the mass-loss rate (dM/dt) increases with radius, similar to its expansion velocity ($V$) - taking $V\propto r$, we find $dM/dt \propto r^{0.9-2.2}$. The mass in the ultra-cold outflow is $\gtrsim3.3$ Msun, and the Boomerang's main-sequence progenitor mass is $\gtrsim4$ Msun. Our high angular resolution ($\sim$0".3) CO J=3-2 map shows the inner bipolar nebula's precise, highly-collimated shape, and a dense central waist of size (FWHM) $\sim$1740 AU$\times275$ AU. The molecular gas and the dust as seen in scattered light via optical HST imaging show a detailed correspondence. The waist shows a compact core in thermal dust emission at 0.87-3.3 mm, which harbors $(4-7)\times10^{-4}$ Msun~of very large ($\sim$mm-to-cm sized), cold ($\sim20-30$ K) grains. The central waist (assuming its outer regions to be expanding) and fast bipolar outflow have expansion ages of $\lesssim1925$ yr and $\le1050$ yr: the "jet-lag" (i.e., torus age minus the fast-outflow age) in the Boomerang supports models in which the primary star interacts directly with a binary companion. We argue that this interaction resulted in a common-envelope configuration while the Boomerang's primary was an RGB or early-AGB star, with the companion finally merging into the primary's core, and ejecting the primary's envelope that now forms the ultra-cold outflow.

Preplanetary nebulae and planetary nebulae are evolved, mass-losing stellar objects that show a wide variety of morphologies. Many of these nebulae consist of outer structures that are nearly spherical (spiral/shell/arc/halo) and inner structures that are highly asymmetric (bipolar/multipolar). The coexistence of such geometrically distinct structures is enigmatic because it hints at the simultaneous presence of both wide and close binary interactions, a phenomenon that has been attributed to stellar binary systems with eccentric orbits. Here we report new high-resolution molecular-line observations of the circumstellar spiral-shell pattern of AFGL 3068, an asymptotic giant branch star transitioning to the preplanetary nebula phase. The observations clearly reveal that the dynamics of the mass loss is influenced by the presence of an eccentric-orbit binary. This quintessential object opens a new window on the nature of deeply embedded binary stars through the circumstellar spiral-shell patterns that reside at distances of several thousand Astronomical Units from the stars.

Extremely high velocity emission likely related to jets is known to occur in some proto-Planetary Nebulae. However, the molecular complexity of this kinematic component is largely unknown. We observed the known extreme outflow from the proto-Planetary Nebula IRAS 16342-3814, a prototype water fountain, in the full frequency range from 73 to 111 GHz with the RSR receiver on the Large Millimetre Telescope. We detected the molecules SiO, HCN, SO, and $^{13}$CO. All molecular transitions, with the exception of the latter are detected for the first time in this source, and all present emission with velocities up to a few hundred km s$^{-1}$. IRAS 16342-3814 is therefore the only source of this kind presenting extreme outflow activity simultaneously in all these molecules, with SO and SiO emission showing the highest velocities found of these species in proto-Planetary Nebulae. To be confirmed is a tentative weak SO component with a FWHM $\sim$ 700 km s$^{-1}$. The extreme outflow gas consists of dense gas (n$_{\rm H_2} >$ 10$^{4.8}$--10$^{5.7}$ cm$^{-3}$), with a mass larger than $\sim$ 0.02--0.15 M$_{\odot}$. The relatively high abundances of SiO and SO may be an indication of an oxygen-rich extreme high velocity gas.

We have mapped 12CO J=3-2 and other molecular lines from the "water-fountain" bipolar pre-planetary nebula (PPN) IRAS 16342-3814 with ~0."35 resolution using ALMA. We find (i) two very high-speed knotty, jet-like molecular outflows, (ii) a central high-density (> few x 10^6 cm^{-3}), expanding torus of diameter 1300 AU, and (iii) the circumstellar envelope of the progenitor AGB, generated by a sudden, very large increase in the mass-loss rate to >3.5 x 10^{-4} Msun/yr in the past ~455 yr. Strong continuum emission at 0.89 mm from a central source (690 mJy), if due to thermally-emitting dust, implies a substantial mass (0.017 Msun) of very large (~mm-sized) grains. The measured expansion ages of the above structural components imply that the torus (age~160 yr) and the younger high-velocity outflow (age~110 yr) were formed soon after the sharp increase in the AGB mass-loss rate. Assuming a binary model for the jets in IRAS 16342, the high momentum rate for the dominant jet-outflow in IRAS 16342 implies a high minimum accretion rate, ruling out standard Bondi-Hoyle-Lyttleton wind accretion and wind Roche lobe overflow (RLOF) models with white-dwarf or main-sequence companions. Most likely, enhanced RLOF from the primary or accretion modes operating within common envelope evolution are needed.

Context: Low- and intermediate-mass stars lose most of their stellar mass at the end of their lives on the asymptotic giant branch (AGB). Determining gas and dust mass-loss rates (MLRs) is important in quantifying the contribution of evolved stars to the enrichment of the interstellar medium. Aims: Attempt to, for the first time, spectrally resolve CO thermal line emission in a small sample of AGB stars in the Large Magellanic Cloud. Methods: ALMA was used to observe 2 OH/IR stars and 4 carbon stars in the LMC in the CO J= 2-1 line. Results: We present the first measurement of expansion velocities in extragalactic carbon stars. All four C-stars are detected and wind expansion velocities and stellar velocities are directly measured. Mass-loss rates are derived from modelling the spectral energy distribution and Spitzer/IRS spectrum with the DUSTY code. Gas-to-dust ratios are derived that make the predicted velocities agree with the observed ones. The expansion velocities and MLRs are compared to a Galactic sample of well-studied relatively low MLRs stars supplemented with "extreme" C-stars that have properties more similar to the LMC targets. Gas MLRs derived from a simple formula are significantly smaller than derived from the dust modelling, indicating an order of magnitude underestimate of the estimated CO abundance, time-variable mass loss, or that the CO intensities in LMC stars are lower than predicted by the formula derived for Galactic objects. This could be related to a stronger interstellar radiation field in the LMC. Conclusions: Although the LMC sample is small and the comparison to Galactic stars is non-trivial because of uncertainties in their distances it appears that for C stars the wind expansion velocities in the LMC are lower than in the solar neighbourhood, while the MLRs appear similar. This is in agreement with dynamical dust-driven wind models.

We report spatially resolved (FWHM$\sim3.8-4.6"$) mid-IR imaging observations of the planetary nebula (PN) NGC 7027 taken with the 2.5-m telescope aboard the Stratospheric Observatory for Infrared Astronomy (SOFIA). Images of NGC 7027 were acquired at 6.3, 6.6, 11.1, 19.7, 24.2, 33.6, and 37.1 $\mu\mathrm{m}$ using the Faint Object Infrared Camera for the SOFIA Telescope (FORCAST).The observations reveal emission from Polycyclic Aromatic Hydrocarbon (PAH) and warm dust ($T_D\sim90$ K) from the illuminated inner edge of the molecular envelope surrounding the ionized gas and central star. The DustEM code was used to fit the spectral energy distribution of fluxes obtained by FORCAST and the archival infrared spectrum of NGC 7027 acquired by the Short Wavelength Spectrometer (SWS) on the Infrared Space Observatory (ISO). Best-fit dust models provide a total dust mass of $5.8^{+2.3}_{-2.6}\times10^{-3}$ $\mathrm{M}_\odot$, where carbonaceous large ($a=1.5$ $\mu$m) and very small ($a \sim12\AA$) grains, and PAHs ($3.1\AA<a<12\AA$) compose 96.5, 2.2, and 1.3 $\%$ of the dust by mass, respectively. The 37 $\mu$m optical depth map shows minima in the dust column density at regions in the envelope that are coincident with a previously identified collimated outflow from the central star. The optical depth minima are also spatially coincident with enhancements in the 6.2 $\mu$m PAH feature, which is derived from the 6.3 and 6.6 $\mu$m maps. We interpret the spatial anti-correlation of the dust optical depth and PAH 6.2 $\mu$m feature strength and their alignment with the outflow from the central star as evidence of dust processing and rapid PAH formation via grain-grain collisions in the post-shock environment of the dense ($n_H\sim10^5\,\mathrm{cm}^{-3}$) photo-dissociation region (PDR) and molecular envelope.

The well-studied carbon star, V Hya, showing evidence for high-speed, collimated outflows and dense equatorial structures, is a key object in the study of the poorly understood transition of AGB stars into aspherical planetary nebulae. Using the STIS instrument onboard HST, we have obtained high spatial-resolution long-slit optical spectra of V Hya that show high-velocity emission in [SII] and [FeII] lines. Our dataset, spanning three epochs spaced apart by a year during each of two periods (in 2002-2004 and 2011-2013), shows that V Hya ejects high-speed (about 200-250 km/s) bullets once every ~8.5 yr. The ejection axis flip-flops around a roughly eastern direction, both in and perpendicular to the sky-plane, and the radial velocities of the ejecta also vary in concert between low and high values. We propose a model in which the bullet ejection is associated with the periastron passage of a binary companion in an eccentric orbit around V Hya with an orbital period of ~8.5 yr. The flip-flop phenomenon is likely the result of collimated ejection from an accretion disk (produced by gravitational capture of material from the primary) that is warped and precessing, and/or that has a magnetic field that is misaligned with that of the companion or the primary star. We show how a previously observed 17-year period in V Hya's light-cycle can also be explained in our model. Additionally, we describe how the model proposed here can be extended to account for multipolar nebulae.

The Infrared Spectrograph on the Spitzer Space Telescope observed 184 carbon stars in the Magellanic Clouds. This sample reveals that the dust-production rate (DPR) from carbon stars generally increases with the pulsation period of the star. The composition of the dust grains follows two condensation sequences, with more SiC condensing before amorphous carbon in metal-rich stars, and the order reversed in metal-poor stars. MgS dust condenses in optically thicker dust shells, and its condensation is delayed in more metal-poor stars. Metal-poor carbon stars also tend to have stronger absorption from C2H2 at 7.5 um. The relation between DPR and pulsation period shows significant apparent scatter, which results from the initial mass of the star, with more massive stars occupying a sequence parallel to lower-mass stars, but shifted to longer periods. Accounting for differences in the mass distribution between the carbon stars observed in the Small and Large Magellanic Clouds reveals a hint of a subtle decrease in the DPR at lower metallicities, but it is not statistically significant. The most deeply embedded carbon stars have lower variability amplitudes and show SiC in absorption. In some cases they have bluer colors at shorter wavelengths, suggesting that the central star is becoming visible. These deeply embedded stars may be evolving off of the AGB and/or they may have non-spherical dust geometries.

We report interferometric mapping of the bipolar pre-planetary nebula IRAS 08005-2356 with an angular-resolution of ~1"-5", using the Submillimeter Array (SMA), in the 12CO J=2-1, 3-2, 13CO J=2-1 and SiO J=5-4 (v=0) lines. Single-dish observations, using the SMT 10-m, were made in these lines as well as in the CO J=4-3 and SiO J-6-5 (v=0) lines. The lines profiles are very broad, showing the presence of a massive (>0.1 Msun), extreme high-velocity outflow (V~200 km/s) directed along the nebular symmetry axis derived from the HST imaging of this object. The outflow's scalar momentum far exceeds that available from radiation pressure of the central post-AGB star, and it may be launched from an accretion disk around a main-sequence companion. We provide indirect evidence for such a disk from its previously published, broad H-alpha emission profile, which we propose results from Ly-beta emission generated in the disk followed by Raman-scattering in the innermost regions of a fast, neutral wind.

We report the results of a pilot survey for X-ray emission from a newly discovered class of AGB stars with far-ultraviolet excesses (fuvAGB stars) using XMM-Newton and Chandra. We detected X-ray emission in 3 of 6 fuvAGB stars observed -- the X-ray fluxes are found to vary in a stochastic or quasi-periodic manner on roughly hour-long times-scales, and simultaneous UV observations using the Optical Monitor on XMM for these sources show similar variations in the UV flux. These data, together with previous studies, show that X-ray emission is found only in fuvAGB stars. From modeling the spectra, we find that the observed X-ray luminosities are ~(0.002-0.2 ) Lsun, and the X-ray emitting plasma temperatures are ~(35-160) x 10^6 K. The high X-ray temperatures argue against the emission arising in stellar coronae, or directly in an accretion shock, unless it occurs on a WD companion. However, none of the detected objects is a known WD-symbiotic star, suggesting that if WD companions are present, they are relatively cool (<20,000 K). In addition, the high X-ray luminosities specifically argue against emission originating in the coronae of main-sequence companions. We discuss several models for the X-ray emission and its variability and find that the most likely scenario for the origin of the X-ray (and FUV) emission involves accretion activity around a companion star, with confinement by strong magnetic fields associated with the companion and/or an accretion disk around it.

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 have discovered two extended half-ring structures in a far-ultraviolet image taken with the GALEX satellite of the well-known mass-losing carbon star CIT 6 (RW LMi). The northern (southern) ring is brighter (fainter) with a diameter of about 15 arcmin (18 arcmin). These structures most likely represent the astrosphere resulting from the shock interaction of CIT 6's molecular wind with the Warm Interstellar Medium, as it moves through the latter. These data provide a direct estimate of the size of CIT 6's circumstellar envelope that is a factor ~20 larger than previous estimates based on CO millimeter-wave line data. We find that CIT 6 has been undergoing heavy mass-loss for at least 93,000 yr and the total envelope mass is 0.29 Msun or larger, assuming a constant mass-loss rate of 3.2 x 10^(-6) Msun/yr. Assuming that the shock front has reached a steady-state and CIT 6's motion relative to the ISM is in the sky-plane, we measure the termination-shock standoff distance directly from the image and find that CIT 6 is moving at a speed of about 39 (0.17 cm^(-3)/n_ISM)^(1/2) km/s through the interstellar medium around it. However, comparisons with published numerical simulations and analytical modelling shows that CIT 6's forward shock (the northern ring) departs from the parabolic shape expected in steady-state. We discuss several possible explanations for this departure.

Infrared spectra of carbon-rich objects which have evolved off the asymptotic giant branch reveal a range of dust properties, including fullerenes, polycyclic aromatic hydrocarbons (PAHs), aliphatic hydrocarbons, and several unidentified features, including the 21 um emission feature. To test for the presence of fullerenes, we used the position and width of the feature at 18.7-18.9 um and examined other features at 17.4 and 6-9 um. This method adds three new fullerene sources to the known sample, but it also calls into question three previous identifications. We confirm that the strong 11 um features seen in some sources arise primarily from SiC, which may exist as a coating around carbonaceous cores and result from photo-processing. Spectra showing the 21 um feature usually show the newly defined Class D PAH profile at 7-9 um. These spectra exhibit unusual PAH profiles at 11-14 um, with weak contributions at 12.7 um, which we define as Class D1, or show features shifted to ~11.4, 12.4, and 13.2 um, which we define as Class D2. Alkyne hydrocarbons match the 15.8 um feature associated with 21 um emission. Sources showing fullerene emission but no PAHs have blue colors in the optical, suggesting a clear line of sight to the central source. Spectra with 21 um features and Class D2 PAH emission also show photometric evidence for a relatively clear line of sight to the central source. The multiple associations of the 21 um feature to aliphatic hydrocarbons suggest that the carrier is related to this material in some way.

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.

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.

This paper reports variations of polycyclic aromatic hydrocarbons (PAHs) features that were found in Spitzer Space Telescope spectra of carbon-rich post-asymptotic giant branch (post-AGB) stars in the Large Magellanic Cloud (LMC). The paper consists of two parts. The first part describes our Spitzer spectral observing programme of 24 stars including post-AGB candidates. The latter half of this paper presents the analysis of PAH features in 20 carbon-rich post-AGB stars in the LMC, assembled from the Spitzer archive as well as from our own programme. We found that five post-AGB stars showed a broad feature with a peak at 7.7 micron, that had not been classified before. Further, the 10--13 micron PAH spectra were classified into four classes, one of which has three broad peaks at 11.3, 12.3 and 13.3 micron rather than two distinct sharp peaks at 11.3 and 12.7 micron, as commonly found in HII regions. Our studies suggest that PAHs are gradually processed while the central stars evolve from post-AGB phase to PNe, changing their composition before PAHs are incorporated into the interstellar medium. Although some metallicity dependence of PAH spectra exists, the evolutionary state of an object is more significant than its metallicity in determining the spectral characteristics of PAHs for LMC and Galactic post-AGB stars.

We have imaged the bipolar planetary nebula M2-9 using SOFIA's FORCAST instrument in six wavelength bands between 6.6 and 37.1 $\mu m$. A bright central point source, unresolved with SOFIA's $\sim$ 4${''}$-to-5${''}$ beam, is seen at each wavelength, and the extended bipolar lobes are clearly seen at 19.7 $\mu m$ and beyond. The photometry between 10 and 25 $\mu m$ is well fit by the emission predicted from a stratified disk seen at large inclination, as has been proposed for this source by Lykou et al and by Smith and Gehrz. The principal new results in this paper relate to the distribution and properties of the dust that emits the infrared radiation. In particular, a considerable fraction of this material is spread uniformly through the lobes, although the dust density does increase at the sharp outer edge seen in higher resolution optical images of M2-9. The dust grain population in the lobes shows that small ($<$ 0.1 $\mu m$) and large ($>$ 1 $\mu m$) particles appear to be present in roughly equal amounts by mass. We suggest that collisional processing within the bipolar outflow plays an important role in establishing the particle size distribution.

CRL 618 is a well-studied pre-planetary nebula. It has multiple highly collimated optical lobes, fast molecular outflows along the optical lobes, and an extended molecular envelope that consists of a dense torus in the equator and a tenuous round halo. Here we present our observations of this source in CO J=3-2 and HCN J=4-3 obtained with the Submillimeter Array at up to ~ 0.3" resolutions. We spatially resolve the fast molecular-outflow region previously detected in CO near the central star and find it to be composed of multiple outflows that have similar dynamical ages, and are oriented along the different optical lobes. We also detect fast molecular outflows further away from the central star near the tips of the extended optical lobes and a pair of equatorial outflows inside the dense torus. We find that two episodes of bullet ejections in different directions are needed, one producing the fast molecular outflows near the central star, and one producing the fast molecular outflows near the tips of the extended optical lobes. One possibility to launch these bullets is the magneto-rotational explosion of the stellar envelope.

The Boomerang Nebula is the coldest known object in the Universe, and an extreme member of the class of Pre-Planetary Nebulae, objects which represent a short-lived transitional phase between the AGB and Planetary Nebula evolutionary stages. Previous single-dish CO (J=1-0) observations (with a 45 arcsec beam) showed that the high-speed outflow in this object has cooled to a temperature significantly below the temperature of the cosmic background radiation. Here we report the first observations of the Boomerang with ALMA in the CO J=2-1 and J=1-0 lines to resolve the structure of this ultra-cold nebula. We find a central hourglass-shaped nebula surrounded by a patchy, but roughly round, cold high-velocity outflow. We compare the ALMA data with visible-light images obtained with HST and confirm that the limb-brightened bipolar lobes seen in these data represent hollow cavities with dense walls of molecular gas and dust producing both the molecular-emission-line and scattered-light structures seen at millimeter and visible wavelengths. The large diffuse biconical shape of the nebula seen in the visible is likely due to preferential illumination of the cold high-velocity outflow. We find a compact source of millimeter-wave continuum in the nebular waist -- these data, together with sensitive upper limits on the radio continuum using observations with ATCA, indicate the presence of a substantial mass of very large (mm-sized) grains in the waist of the nebula. Another unanticipated result is the detection of CO emission regions beyond the ultracold region which indicate the re-warming of the cold gas, most likely due to photoelectric grain heating.

CRL 618 is a well-studied pre-planetary nebula. We have mapped its central region in continuum and molecular lines with the Submillimeter Array at 350 GHz at ~ 0.3" to 0.5" resolutions. Two components are seen in 350 GHz continuum: (1) a compact emission at the center tracing the dense inner part of the H II region previously detected in 23 GHz continuum and it may trace a fast ionized wind at the base, and (2) an extended thermal dust emission surrounding the H II region, tracing the dense core previously detected in HC3N at the center of the circumstellar envelope. The dense core is dusty and may contain mm-sized dust grains. It may have a density enhancement in the equatorial plane. It is also detected in carbon chain molecules HC3N and HCN, and their isotopologues, with higher excitation lines tracing closer to the central star. It is also detected in C2H3CN toward the innermost part. Most of the emission detected here arises within ~ 630 AU (~0.7") from the central star. A simple radiative transfer model is used to derive the kinematics, physical conditions, and the chemical abundances in the dense core. The dense core is expanding and accelerating, with the velocity increasing roughly linearly from ~ 3 km/s in the innermost part to ~ 16 km/s at 630 AU. The mass-loss rate in the dense core is extremely high with a value of ~ 1.15x10^{-3} solar mass per year. The dense core has a mass of ~ 0.47 solar mass and a dynamical age of ~ 400 yrs. It could result from a recent enhanced heavy mass-loss episode that ends the AGB phase. The isotopic ratios of 12C/13C and 14N/15N are 9$\pm$4 and 150$\pm50$, respectively, both lower than the solar values.

We report the detection of strong and compact molecular line emission (in the CO J=3-2, 4-3, 6-5, 7-6, 13CO J=3-2, HCN and HCO+ J=4-3 transitions) from a cometary-shaped object (Carina-frEGG1) in the Carina star-forming region (SFR) previously classified as a photoevaporating protoplanetary disk (proplyd). We derive a molecular mass of 0.35 Msun for Carina-frEGG1, which shows that it is not a proplyd, but belongs to a class of free-floating evaporating gas globules (frEGGs) recently found in the Cygnus SFR by Sahai, Morris & Claussen (2012). Archival Adaptive Optics near-IR (Ks) images show a central hourglass-shaped nebula. The derived source luminosity (about 8-18 Lsun), the hourglass morphology, and the presence of collimated jets seen in HST images, imply the presence of a jet-driving, young, low-mass star deeply embedded in the dust inside Carina-frEGG1. Our results suggest that the true nature of many or most such cometary-shaped objects seen in massive SFRs and previously labelled as proplyds has been misunderstood, and that these are really frEGGs.

The theory of massive star formation currently suffers from a scarce observational base of massive young stellar objects to compare with. In this paper, we present OVRO CO(1-0), 13CO(1-0), C18O(1-0), and 2.6 mm continuum images of the infrared source IRAS 19520+2759 together with complementary single-dish observations of CS(1-0), obtained with the 34m antenna DSS-54 at the Madrid Deep Space Communications Complex, as well as archive images at different wavelengths. As a result from our work, IRAS 19520+2759, with a controversial nature in the past, is firmly established as a massive young stellar object associated with a strong and compact millimetre source and driving a collimated outflow. In addition, a second fainter millimetre source is discovered about 4 arcsec to the south, which is also driving an outflow. Furthermore, the two millimetre sources are associated with C18O clumps elongated perpendicularly to the outflows, which may be related to rotating toroids. The masses of gas and dust of the millimetre sources are estimated to be around 100 and 50 Msun. MM1, the dominant source at all wavelengths, with a total luminosity of (1-2)x10^5 Lsun at 9 kpc, is however not associated with 6 cm emission down to a rms noise level of 0.1 mJy. We propose that IRAS 19520+2759 could be an example of the recent theoretical prediction of `bloated' or `swollen' star, i.e., a massive young stellar object whose radius has increased due to effects of accretion at a high-mass accretion rate.

The majority of Proto-planetary nebulae (PPN) are observed to have bipolar morphologies. The majority of mature PN are observed to have elliptical shapes. In this paper we address the evolution of PPN/PN morphologies attempting to understand if a transition from strongly bipolar to elliptical shape can be driven by changes in the parameters of the mass loss process. To this end we present 2.5D hydrodynamical simulations of mass loss at the end stages of stellar evolution for intermediate mass stars. We track changes in wind velocity, mass loss rate and mass loss geometry. In particular we focus on the transition from mass loss dominated by a short duration jet flow (driven during the PPN phase) to mass loss driven by a spherical fast wind (produced by the central star of the PN). We address how such changes in outflow characteristics can change the nebula from a bipolar to an elliptical morphology. Our results show that including a period of jet formation in the temporal sequence of PPN to PN produces realistic nebular synthetic emission geometries. More importantly such a sequence provides insight, in principle, into the apparent difference in morphology statistics characterizing PPN and PN systems. In particular we find that while jet driven PPN can be expected to be dominated by bipolar morphologies, systems that begin with a jet but are followed by a spherical fast wind will evolve into elliptical nebulae. Furthermore, we find that spherical nebulae are highly unlikely to ever derive from either bipolar PPN or elliptical PN.

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.

We have obtained new velocity-resolved spectra of the [CII] 158 micron line towards the Ring Nebula in Lyra (NGC 6720), one of the best-studied planetary nebulae, in order to probe its controversial 3-dimensional structure and to estimate the mass of circumstellar material in this object. We used the Terahertz receiver GREAT aboard the SOFIA airborne telescope to obtain the [CII] spectra at eight locations within and outside the bright optical ring of NGC 6720. Emission was detected at all positions except for the most distant position along the nebula's minor axis, and generally covers a broad velocity range, ~50 km/s (FWZI), except at a position along the major axis located just outside the optical ring, where it is significantly narrower (~25 km/s). The one narrow spectrum appears to be probing circumstellar material lying outside the main nebular shell that has not been accelerated by past fast wind episodes from the central star, and therefore most likely comes from equatorial and/or low-latitude regions of this bipolar nebula. Along lines-of-sight passing within about 10 arcsec of the nebular center, the CII column density is a factor 46 higher than the CO column density. The total mass of gas associated with the [CII] emission inside a circular region of diameter 87.5 arcsec is at least 0.11 solar masses. A significant amount of [CII] flux arises from a photodissociation region immediately outside the bright optical ring, where we find a CII to CO ratio of 6.5, lower than that seen towards the central region. Comparing our data with lower-quality CI spectra, which indicate similarly large CI/CO ratios in NGC 6720, we conclude that the bulk of elemental carbon in NGC 6720 is divided roughly equally between CII and CI, and that the emissions from these species are far more robust tracers of circumstellar material than CO in this object and other evolved planetary nebulae.

We report multi-wavelength observations of the far-infrared source IRAS 20324+4057, including high-resolution optical imaging with HST, and ground-based near-infrared, millimeter-wave and radio observations. These data show an extended, limb-brightened, tadpole-shaped nebula with a bright, compact, cometary nebula located inside the tadpole head. Our molecular line observations indicate that the Tadpole is predominantly molecular, with a total gas mass exceeding 3.7 Msun. Our radio continuum imaging, and archival Spitzer IRAC images, show the presence of additional tadpole-shaped objects in the vicinity of IRAS 20324+4057 that share a common E-W head-tail orientation: we propose that these structures are small, dense molecular cores that originated in the Cygnus cloud and are now being (i) photoevaporated by the ultraviolet radiation field of the Cyg OB2 No. 8 cluster located to the North-West, and (ii) shaped by ram pressure of a distant wind source or sources located to the West, blowing ablated and photoevaporated material from their heads eastwards. The ripples in the tail of the Tadpole are interpreted in terms of instabilities at the interface between the ambient wind and the dense medium of the former.

S-type AGB stars are thought to be in the transitional phase between M-type and C-type AGB stars. Because of their peculiar chemical composition, one may expect a strong influence of the stellar C/O ratio on the molecular chemistry and the mineralogy of the circumstellar dust. In this paper, we present a large sample of 87 intrinsic galactic S-type AGB stars, observed at infrared wavelengths with the Spitzer Space Telescope, and supplemented with ground-based optical data. On the one hand, we derive the stellar parameters from the optical spectroscopy and photometry, using a grid of model atmospheres. On the other, we decompose the infrared spectra to quantify the flux-contributions from the different dust species. Finally, we compare the independently determined stellar parameters and dust properties. For the stars without significant dust emission, we detect a strict relation between the presence of SiS absorption in the Spitzer spectra and the C/O ratio of the stellar atmosphere. These absorption bands can thus be used as an additional diagnostic for the C/O ratio. For stars with significant dust emission, we define three groups, based on the relative contribution of certain dust species to the infrared flux. We find a strong link between group-membership and C/O ratio. We show that these groups can be explained by assuming that the dust-condensation can be cut short before silicates are produced, while the remaining free atoms and molecules can then form the observed magnesium sulfides or the carriers of the unidentified 13 and 20 micron features. Finally, we present the detection of emission features attributed to molecules and dust characteristic to C-type stars, such as molecular SiS, hydrocarbons and magnesium sulfide grains. We show that we often detect magnesium sulfides together with molecular SiS and we propose that it is formed by a reaction of SiS molecules with Mg.

Binarity is believed to dramatically affect the history and geometry of mass loss in AGB and post-AGB stars, but observational evidence of binarity is sorely lacking. As part of a project to look for hot binary companions to cool AGB stars using the GALEX archive, we have discovered a late-M star, Y Gem, to be a source of strong and variable UV emission. Y Gem is a prime example of the success of our technique of UV imaging of AGB stars in order to search for binary companions. Y Gem's large and variable UV flux makes it one of the most prominent examples of a late AGB star with a mass accreting binary companion. The UV emission is most likely due to emission associated with accretion activity and a disk around a main-sequence companion star. The physical mechanism generating the UV emission is extremely energetic, with an integrated luminosity of a few L(sun) at its peak. We also find weak CO J=2-1 emission from Y Gem with a very narrow line profile (FWHM of 3.4 km/s). Such a narrow line is unlikely to arise in an outflow, and is consistent with emission from an orbiting, molecular reservoir of radius 300 AU. Y Gem may be the progenitor of the class of post-AGB stars which are binaries and possess disks but no outflows.

We report the results of a pilot multiwavelength survey in the radio continuum (X, Ka and Q bands, i.e., from 3.6 cm to 7 mm) carried out with the EVLA in order to confirm the presence of very large dust grains in dusty disks and torii around the central stars in a small sample of post-AGB objects, as inferred from millimeter and sub-millimeter observations. Supporting mm-wave observations were also obtained with CARMA towards three of our sources. Our EVLA survey has resulted in a robust detection of our most prominent submm emission source, the pre-planetary nebula IRAS 22036+5306, in all three bands, and the disk-prominent post-AGB object, RV Tau, in one band. The observed fluxes are consistent with optically-thin free-free emission, and since they are insignificant compared to their submm/mm fluxes, we conclude that the latter must come from substantial masses of cool, large (mm-sized) grains. We find that the power-law emissivity in the cm-to-submm range for the large grains in IRAS22036 is nu^{beta}, with beta=1-1.3. Furthermore, the value of beta in the 3 to 0.85 mm range for the three disk-prominent post-AGB sources (beta < or = 0.4) is significantly lower than that of IRAS22036, suggesting that the grains in post-AGB objects with circumbinary disks are likely larger than those in the dusty waists of pre-planetary nebulae.

Using Hubble Space Telescope images of 119 young planetary nebulae, most of which have not previously been published, we have devised a comprehensive morphological classification system for these objects. This system generalizes a recently devised system for pre-planetary nebulae, which are the immediate progenitors of planetary nebulae (PNs). Unlike previous classification studies, we have focussed primarily on young PNs rather than all PNs, because the former best show the influences or symmetries imposed on them by the dominant physical processes operating at the first and primary stage of the shaping process. Older PNs develop instabilities, interact with the ambient interstellar medium, and are subject to the passage of photoionization fronts, all of which obscure the underlying symmetries and geometries imposed early on. Our classification system is designed to suffer minimal prejudice regarding the underlying physical causes of the different shapes and structures seen in our PN sample, however, in many cases, physical causes are readily suggested by the geometry, along with the kinematics that have been measured in some systems. Secondary characteristics in our system such as ansae indicate the impact of a jet upon a slower-moving, prior wind; a waist is the signature of a strong equatorial concentration of matter, whether it be outflowing or in a bound Keplerian disk, and point symmetry indicates a secular trend, presumably precession, in the orientation of the central driver of a rapid, collimated outflow.

Proto-planetary nebulae (pPN) and planetary nebulae (PN) seem to be formed by interacting winds from asymptotic giant branch (AGB) stars. The observational issue that most pPN are bipolar but most older PN are elliptical is addressed. We present 2.5D hydrodynamical numerical simulations of episodic cooling interacting winds to investigate the long term evolution of PN morphologies. We track wind acceleration, decrease in mass-loss and episodic change in wind geometry from spherical (AGB) to collimated (pPN) and back to spherical again (PN). This outflow sequence is found to produce realistic PN dynamics and morphological histories. Effects from different AGB distributions and jet duty cycles are also investigated.

We report on the results of two epochs of Very Long Baseline Array (VLBA) observations of the 22 GHz water masers toward IRAS 19190+1102. The water maser emission from this object shows two main arc-shaped formations perpendicular to their NE-SW separation axis. The arcs are separated by ~280 mas in position, and are expanding outwards at an angular rate of 2.35 mas/yr. We detect maser emission at velocities between -53.3 km/s to +78.0 km/s and there is a distinct velocity pattern where the NE masers are blueshifted and the SW masers are redshifted. The outflow has a three-dimensional outflow velocity of 99.8 km/s and a dynamical age of about 59 yr. A group of blueshifted masers not located along the arcs shows a change in velocity of more than 35 km/s between epochs, and may be indicative of the formation of a new lobe. These observations show that IRAS 19190+1102 is a member of the class of "water fountain"' pre-planetary nebulae displaying bipolar structure

Aims: We aim at measuring mass-loss rates and the luminosities of a statistically large sample of Galactic bulge stars at several galactocentric radii. The sensitivity of previous infrared surveys of the bulge has been rather limited, thus fundamental questions for late stellar evolution, such as the stage at which substantial mass-loss begins on the red giant branch and its dependence on fundamental stellar properties, remain unanswered. We aim at providing evidence and answers to these questions. Methods: To this end, we observed seven 15 times 15 arcmin^2 fields in the nuclear bulge and its vicinity with unprecedented sensitivity using the IRAC and MIPS imaging instruments on-board the Spitzer Space Telescope. In each of the fields, tens of thousands of point sources were detected. Results: In the first paper based on this data set, we present the observations, data reduction, the final catalogue of sources, and a detailed comparison to previous mid-IR surveys of the Galactic bulge, as well as to theoretical isochrones. We find in general good agreement with other surveys and the isochrones, supporting the high quality of our catalogue.

We have discovered a very extended shock structure (i.e., with a diameter of about 24 arcmin) surrounding the well-known carbon star IRC+10216 in ultraviolet images taken with the GALEX satellite. We conclude that this structure results from the interaction of IRC+10216's molecular wind with the interstellar medium (ISM), as it moves through the latter. All important structural features expected from theoretical models of such interactions are identified: the termination shock, the astrosheath, the astropause, the bowshock, and an astrotail (with vortices). The extent of the astropause provides new lower limits to the envelope age (69,000 years) and mass (1.4 solar masses, for a mass-loss rate of 2 x 10^(-5) solar masses/year). From the termination-shock standoff distance, we find that IRC+10216 is moving at a speed of about > 91 km/s [1 cm^(-3)/n_(ISM)]^1/2 through the local ISM.

We present archival Spitzer IRS spectra of 19 luminous 8 micron selected sources in the Large Magellanic Cloud (LMC). The object classes derived from these spectra and from an additional 24 spectra in the literature are compared with classifications based on 2MASS/MSX (J, H, K, and 8 micron) colors in order to test the "JHK8" classification scheme (Kastner et al. 2008). The IRS spectra confirm the classifications of 22 of the 31 sources that can be classified under the JHK8 system. The spectroscopic classification of 12 objects that were unclassifiable in the JHK8 scheme allow us to characterize regions of the color-color diagrams that previously lacked spectroscopic verification, enabling refinements to the JHK8 classification system. The results of these new classifications are consistent with previous results concerning the identification of the most infrared-luminous objects in the LMC. In particular, while the IRS spectra reveal several new examples of asymptotic giant branch (AGB) stars with O-rich envelopes, such objects are still far outnumbered by carbon stars (C-rich AGB stars). We show that Spitzer IRAC/MIPS color-color diagrams provide improved discrimination between red supergiants and oxygen-rich and carbon-rich asymptotic giant branch stars relative to those based on 2MASS/MSX colors. These diagrams will enable the most luminous IR sources in Local Group galaxies to be classified with high confidence based on their Spitzer colors. Such characterizations of stellar populations will continue to be possible during Spitzer's warm mission, through the use of IRAC [3.6]-[4.5] and 2MASS colors.

We have carried out high-resolution spectroscopic observations of the carbon star V Hya, covering the 4.6 micron band of CO. These data, taken over 7 epochs, show that the circumstellar environment of V Hya consists of a complex high-velocity (HV) outflow containing at least six kinematic components with expansion velocities ranging between 70 and 120 km/s, together with a slow-moving normal outflow at about 10 km/s. Physical changes occur in the HV outflow regions on a time-scale as short as two days, limiting their extent to < ~ 10^{16} cm. The intrinsic line-width for each HV component is quite large (6-8 km/s) compared to the typical values (~1 km/s) appropriate for normal AGB circumstellar envelopes (CSEs), due to excess turbulence and/or large velocity gradients resulting from the energetic interaction of the HV outflow with the V Hya CSE. We have modelled the absorption features to set constraints on the temperature distribution in, and the mass ejection-rates for gas in the main HV components.

Almost all stars in the 1-8 Msun range evolve through the Asymptotic Giant Branch (AGB), preplanetary nebula (PPN) and planetary nebula (PN) evolutionary phases. Most stars that leave the main sequence in a Hubble time will end their lives in this way. The heavy mass loss which occurs during the AGB phase is important across astrophysics, and the particulate matter crucial for the birth of new solar systems is made and ejected by AGB stars. Yet stellar evolution from the beginning of the AGB phase to the PN phase remains poorly understood. We do not understand how the mass-loss (rate, geometry, temporal history) depends on fundamental stellar parameters or the presence of a binary companion. While the study of evolved non-massive stars has maintained a relatively modest profile in recent decades, we are nonetheless in the midst of a quiet but exciting revolution in this area, driven by new observational results, such as the discovery of jets and disks in stellar environments where these were never expected, and by the recognition of new symmetries such as multipolarity and point-symmetry occuring frequently in the nebulae resulting from the outflows. In this paper we summarise the major unsolved problems in this field, and specify the areas where allocation of effort and resources is most likely to help make significant progress.

Collimated fast winds (CFWs) have been proposed to operate during the post-AGB evolutionary phase (and even earlier during the late AGB phase), responsible for the shaping of pre-planetary nebulae (PPNs) and young planetary nebulae (PNs). This paper is a follow-up to our previous study of CFW models for the well-studied PPN CRL 618. Previously, we compared our CFW models with optical observations of CRL 618 in atomic and ionic lines and found that a CFW with a small opening angle can readily reproduce the highly collimated shape of the northwestern (W1) lobe of CRL 618 and the bow-like structure seen at its tip. In this paper, we compare our CFW models with recent observations of CRL 618 in CO J=2-1, J=6-5, and H2 1-0 S(1). In our models, limb-brightened shell structures are seen in CO and H2 at low velocity arising from the shocked AGB wind in the shell, and can be identified as the low-velocity (LV) components in the observations. However, the shell structure in CO J=2-1 is significantly less extended than that seen in the observations. None of our models can properly reproduce the observed high-velocity (HV) molecular emission near the source along the body of the lobe. In order to reproduce the HV molecular emission in CRL 618, the CFW is required to have a different structure. One possible CFW structure is the cylindrical jet, with the fast wind material confined to a small cross section and collimated to the same direction along the outflow axis.

We report the detection of X-ray emission from the jet-driving symbiotic star MWC 560. We observed MWC 560 with XMM-Newton for 36 ks. We fitted the spectra from the EPIC pn, MOS1 and MOS2 instruments with XSPEC and examined the light curves with the package XRONOS. The spectrum can be fitted with a highly absorbed hard X-ray component from an optically-thin hot plasma, a Gaussian emission line with an energy of 6.1 keV and a less absorbed soft thermal component. The best fit is obtained with a model in which the hot component is produced by optically thin thermal emission from an isobaric cooling flow with a maximum temperature of 61 keV, which might be created inside an optically-thin boundary layer on the surface of the accreting with dwarf. The derived parameters of the hard component detected in MWC 560 are in good agreement with similar objects as CH Cyg, SS7317, RT Cru and T CrB, which all form a new sub-class of symbiotic stars emitting hard X-rays. Our previous numerical simulations of the jet in MWC 560 showed that it should produce detectable soft X-ray emission. We infer a temperature of 0.17 keV for the observed soft component, i.e. less than expected from our models. The total soft X-ray flux (i.e. at < 3 keV) is more than a factor 100 less than predicted for the propagating jet soon after its birth (<0.3 yr), but consistent with the value expected due its decrease with age. The ROSAT upper limit is also consistent with such a decrease. We find aperiodic or quasi-periodic variability on timescales of minutes and hours, but no periodic rapid variability. All results are consistent with an accreting white dwarf powering the X-ray emission and the existence of an optically-thin boundary layer around it.

Asymptotic Giant Branch variables are found to obey period-luminosity relations in the mid-IR similar to those seen at K_S (2.14 microns), even at 24 microns where emission from circumstellar dust is expected to be dominant. Their loci in the M, logP diagrams are essentially the same for the LMC and for NGC6522 in spite of different ages and metallicities. There is no systematic trend of slope with wavelength. The offsets of the apparent magnitude vs. logP relations imply a difference between the two fields of 3.8 in distance modulus. The colours of the variables confirm that a principal period with log P > 1.75 is a necessary condition for detectable mass-loss. At the longest observed wavelength, 24 microns, many semi-regular variables have dust shells comparable in luminosity to those around Miras. There is a clear bifurcation in LMC colour-magnitude diagrams involving 24 micron magnitudes.

We present high angular resolution observations, using the Very Long Baseline Array (VLBA) of the NRAO, of the high-velocity water masers toward the "water-fountain" pre-planetary nebula, IRAS 16342-3814. The detailed structure of the water masers appears to be that of bow shocks on either side of a highly collimated jet. The proper motions of the water masers are approximately equal to the radial velocities; the three-dimensional velocities are approximately +/-180 km/s, which leads to a very short dynamical time-scale of ~100 years. Although we do not find direct evidence for precession of the fast collimated jet, there may be indirect evidence for such precession.

We have carried out a spatial-kinematical study of three proto-planetary nebulae, IRAS 16594-4656, Hen 3-401, and Rob 22. High-resolution H2 images were obtained with NICMOS on the HST and high-resolution spectra were obtained with the Phoenix spectrograph on Gemini-South. IRAS 16594-4656 shows a "peanut-shaped" bipolar structure with H2 emission from the walls and from two pairs of more distant, point-symmetric faint blobs. The velocity structure shows the polar axis to be in the plane of the sky, contrary to the impression given by the more complex visual image and the visibility of the central star, with an ellipsoidal velocity structure. Hen 3-401 shows the H2 emission coming from the walls of the very elongated, open-ended lobes seen in visible light, along with a possible small disk around the star. The bipolar lobes appear to be tilted 10-15 deg with respect to the plane of the sky and their kinematics display a Hubble-like flow. In Rob 22, the H2 appears in the form of an "S" shape, approximately tracing out the similar pattern seen in the visible. H2 is especially seen at the ends of the lobes and at two opposite regions close to the unseen central star. The axis of the lobes is nearly in the plane of the sky. Expansion ages of the lobes are calculated to be approximately 1600 yr (IRAS 16594-4656), 1100 yr (Hen 3-401), and 640 yr (Rob 22), based upon approximate distances.

The search for binarity in AGB stars is of critical importance for our understanding of how planetary nebulae acquire the dazzling variety of aspherical shapes which characterises this class. However, detecting binary companions in such stars has been severely hampered due to their extreme luminosities and pulsations. We have carried out a small imaging survey of AGB stars in ultraviolet light (using GALEX) where these cool objects are very faint, in order to search for hotter companions. We report the discovery of significant far-ultraviolet excesses towards nine of these stars. The far-ultraviolet excess most likely results either directly from the presence of a hot binary companion, or indirectly from a hot accretion disk around the companion.

We present echelle long-slit optical spectra of a sample of objects evolving off the AGB, most of them in the pre-planetary nebula (pPN) phase, obtained with the ESI and MIKE spectrographs at Keck-II and Magellan-I, respectively. The total wavelength range covered with ESI (MIKE) is ~3900 to 10900 A (~3600 to 7200A). In this paper, we focus our analysis mainly on the Halpha profiles. Prominent Halpha emission is detected in half of the objects, most of which show broad Halpha wings (up to ~4000 km/s). In the majority of the Halpha-emission sources, fast, post-AGB winds are revealed by P-Cygni profiles. In ~37% of the objects Halpha is observed in absorption. In almost all cases, the absorption profile is partially filled with emission, leading to complex, structured profiles that are interpreted as an indication of incipient post-AGB mass-loss. All sources in which Halpha is seen mainly in absorption have F-G type central stars, whereas sources with intense Halpha emission span a larger range of spectral types from O to G. Shocks may be an important excitation agent of the close stellar surroundings for objects with late type central stars. Sources with pure emission or P Cygni Halpha profiles have larger J-K color excess than objects with Halpha mainly in absorption, which suggests the presence of warm dust near the star in the former. The two classes of profile sources also segregate in the IRAS color-color diagram in a way that intense Halpha-emitters have dust grains with a larger range of temperatures. (abridged)

We have made a serendipitous discovery of an enigmatic outflow source, IRAS 05506+2414 (hereafter IRAS 05506), as part of a multi-wavelength survey of pre-planetary nebulae (PPNs). The HST optical and near-infrared images show a bright compact central source with a jet-like extension, and a fan-like spray of high-velocity (with radial velocities upto 350 kms/s) elongated knots which appear to emanate from it. These structures are possibly analogous to the near-IR "bullets" seen in the Orion nebula. Interferometric observations at 2.6 mm show the presence of a high-velocity CO outflow and a continuum source also with a faint extension, both of which are aligned with the optical jet structure. IRAS 05506 is most likely not a PPN. We find extended NH3 (1,1) emission towards IRAS 05506; these data together with the combined presence of far-IR emission, H2O and OH masers, and CO and CS J=2-1 emission, strongly argue for a dense, dusty star-forming core associated with IRAS 05506. IRAS 05506 is probably an intermediate-mass or massive protostar, and the very short time-scale (200 yr) of its outflows indicates that it is very young. If IRAS 05506 is a massive star, then the lack of radio continuum and the late G -- early K spectral type we find from our optical spectra implies that in this object we are witnessing the earliest stages of its life, while its temperature is still too low to provide sufficient UV flux for ionisation.

To ascertain the nature of the brightest compact mid-infrared sources in the Large Magellanic Cloud (LMC), we have applied an updated version of the Buchanan et al. (2006) 2MASS-MSX color classification system, which is based on the results of Spitzer Space Telescope spectroscopy, to a mid-infrared (8 micron) flux-limited sample of 250 LMC objects for which 2MASS and MSX photometry is available. The resulting 2MASS-MSX ("JHK8") color-based classifications of these sources, which constitute the most mid-IR-luminous objects in the LMC, were augmented, cross-checked, and corrected where necessary via a variety of independent means, such that only 47 sources retain tentative classifications and only 10 sources cannot be classified at all. The sample is found to consist primarily of carbon-rich AGB stars (~35%), red supergiants (~18%), and compact H II regions (~30%), with additional, small populations of oxygen-rich AGB stars (~4%), dusty, early-type emission-line stars (~3%), and foreground, O-rich AGB stars in the Milky Way (~3%). The very large ratio of C-rich to O-rich objects among the luminous and heavily dust-enshrouded AGB stars in our LMC IR source sample is consistent with the hypothesis that carbon stars form easily in lower metallicity environments. We demonstrate that very luminous C-rich and O-rich AGB stars and red supergiants, identified here primarily on the basis of their JHK8 colors, also appear as distinct clusters in Spitzer IRAC/MIPS color-color diagrams. Thus, in principle, the IRS-based IR photometric classification techniques applied here to the LMC can be applied to any external galaxy whose most luminous IR point sources are detectable and resolvable by 2MASS and Spitzer.

The Egg Nebula has been regarded as the archetype of bipolar proto-planetary nebulae, yet we lack a coherent model that can explain the morphology and kinematics of the nebular and dusty components observed at high-spatial and spectral resolution. Here, we report on two sets of observations obtained with the Keck Adaptive Optics Laser Guide Star: H to M-band NIRC2 imaging, and narrow bandpath K-band OSIRIS 3-D imaging-spectroscopy (through the H2 2.121micron emission line). While the central star or engine remains un-detected at all bands, we clearly resolve the dusty components in the central region and confirm that peak A is not a companion star. The spatially-resolved spectral analysis provide kinematic information of the H_2 emission regions in the eastern and central parts of the nebula and show projected velocities for the H_2 emission higher than 100 km/s. We discuss these observations against a possible formation scenario for the nebular components.

Using the Hubble Space Telescope (HST), we have carried out a survey of candidate preplanetary nebulae (PPNs). We report here our discoveries of objects having well-resolved geometrical structures, and use the large sample of PPNs now imaged with HST (including previously studied objects in this class) to devise a comprehensive morphological classification system for this category of objects. The wide variety of aspherical morphologies which we have found for PPNs are qualitatively similar to those found for young planetary nebulae in previous surveys. We also find prominent halos surrounding the central aspherical shapes in many of our objects -- these are direct signatures of the undisturbed circumstellar envelopes of the progenitor AGB stars. Although the majority of these have surface-brightness distributions consistent with a constant mass-loss rate with a constant expansion velocity, there are also examples of objects with varying mass-loss rates. As in our surveys of young planetary nebulae (PNs), we find no round PPNs. The similarities in morphologies between our survey objects and young PNs supports the view that the former are the progenitors of aspherical planetary nebulae, and that the onset of aspherical structure begins during the PPN phase (or earlier). Thus, the primary shaping of a PN clearly does not occur during the PN phase via the fast radiative wind of the hot central star, but significantly earlier in its evolution.

Using the VLBA, we have observed H2O maser emission in the pre-planetary nebula IRAS 19134+2131(I19134), in which the H2O maser spectrum has two groups of emission features separated in radial velocity by ~100 km/s. We also obtained optical images of I19134 with the HST to locate the bipolar reflection nebula in this source for the first time. The spatio-kinematical structure of the H2O masers indicates the existence of a fast, collimated (precessing) flow having a projected extent of ~140 mas and an expansion rate of ~1.9 mas/yr on the sky plane, which gives a dynamical age of only ~40 yr. The two detected optical lobes are also separated by ~150 mas in almost the same direction as that of the collimated flow. The good agreement between the extent and orientation of the H2O maser outflow and optical lobes suggests that the lobes have been recently formed along the collimated fast flow. The positions of all of the detected maser features have been measured with respect to the reference source J1925+2106 over one year. Therefore we analyzed maser feature motions that consist of the combination of an annual parallax, a secular motion following Galactic rotation, and the intrinsic motions within the flow. We obtain an annual-parallax distance to I19134 of D~8 kpc kpc and estimate its location in the Galaxy to be (R, theta, z)=(7.4 kpc, 62 deg, 0.65 kpc). From the mean motion of the blue-shifted and red-shifted clusters of maser features, we estimate the 3-D secular motion of I19134 to be (V_{R}, V_{theta}, V_{z})=(3, 125, 8) [km/s]. From the height from the Galactic plane, z, and the velocity component perpendicular to the Galactic plane, V_{z}, we estimate a rough upper limit of ~9 M_{sun} to the stellar mass of I19134's progenitor.

We present the results of far-infrared imaging of extended regions around three bipolar pre-planetary nebulae, AFGL 2688, OH 231.8+4.2, and IRAS 16342$-$3814, at 70 and 160 $\mu$m with the MIPS instrument on the Spitzer Space Telescope. After a careful subtraction of the point spread function of the central star from these images, we place constraints on the existence of extended shells and thus on the mass outflow rates as a function of radial distance from these stars. We find no apparent extended emission in AFGL 2688 and OH 231.8+4.2 beyond 100 arcseconds from the central source. In the case of AFGL 2688, this result is inconsistent with a previous report of two extended dust shells made on the basis of ISO observations. We derive an upper limit of $2.1\times10^{-7}$ M$_\odot$ yr$^{-1}$ and $1.0\times10^{-7}$ M$_\odot$ yr$^{-1}$ for the dust mass loss rate of AFGL 2688 and OH 231.8, respectively, at 200 arcseconds from each source. In contrast to these two sources, IRAS 16342$-$3814 does show extended emission at both wavelengths, which can be interpreted as a very large dust shell with a radius of $\sim$ 400 arcseconds and a thickness of $\sim$ 100 arcseconds, corresponding to 4 pc and 1 pc, respectively, at a distance of 2 kpc. However, this enhanced emission may also be galactic cirrus; better azimuthal coverage is necessary for confirmation of a shell. If the extended emission is a shell, it can be modeled as enhanced mass outflow at a dust mass outflow rate of $1.5\times10^{-6}$ M$_\odot$ yr$^{-1}$ superimposed on a steady outflow with a dust mass outflow rate of $1.5\times10^{-7}$ M$_\odot$ yr$^{-1}$. It is likely that this shell has swept up a substantial mass of interstellar gas during its expansion, so these estimates are upper limits to the stellar mass loss rate.

In papers I and II in this series, we presented hydrodynamical simulations of jet models with parameters representative of the symbiotic system MWC 560. These were simulations of a pulsed, initially underdense jet in a high density ambient medium. Since the pulsed emission of the jet creates internal shocks and since the jet velocity is very high, the jet bow shock and the internal shocks are heated to high temperatures and should therefore emit X-ray radiation. In this paper, we investigate in detail the X-ray properties of the jets in our models. We have focused our study on the total X-ray luminosity and its temporal variability, the resulting spectra and the spatial distribution of the emission. Temperature and density maps from our hydrodynamical simulations with radiative cooling presented in the second paper are used together with emissivities calculated with the atomic database ATOMDB. The jets in our models show extended and variable X-ray emission which can be characterized as a sum of hot and warm components with temperatures that are consistent with observations of CH Cyg and R Aqr. The X-ray spectra of our model jets show emission line features which correspond to observed features in the spectra of CH Cyg. The innermost parts of our pulsed jets show iron line emission in the 6.4 - 6.7 keV range which may explain such emission from the central source in R Aqr. We conclude that MWC 560 should be detectable with Chandra or XMM-Newton, and such X-ray observations will provide crucial for understanding jets in symbiotic stars.

We present high-angular resolution images of the post-AGB nebula IRAS18276-1431 (also known as OH17.7-2.0) obtained with the Keck II Adaptive Optics (AO) system in its Natural Guide Star (NGS) mode in the Kp, Lp, and Ms near-infrared bands. We also present supporting optical F606W and F814W HST images as well as interferometric observations of the 12CO(J=1-0), 13CO(J=1-0), and 2.6mm continuum emission with OVRO. The envelope of IRAS18276-1431 displays a clear bipolar morphology in our optical and NIR images with two lobes separated by a dark waist and surrounded by a faint 4.5"x3.4" halo. Our Kp-band image reveals two pairs of radial ``searchlight beams'' emerging from the nebula center and several intersecting, arc-like features. From our CO data we derive a mass of M>0.38[D/3kpc]^2 Msun and an expansion velocity v_exp=17km/s for the molecular envelope. The density in the halo follows a radial power-law proportional to r^-3, which is consistent with a mass-loss rate increasing with time. Analysis of the NIR colors indicates the presence of a compact central source of ~300-500K dust illuminating the nebula in addition to the central star. Modeling of the thermal IR suggests a two-shell structure in the dust envelope: 1) an outer shell with inner and outer radius R_in~1.6E16cm and R_out>~1.25E17cm, dust temperature T_d~105-50K, and a mean mass-loss rate of Mdot~1E-3Msun/yr; and 2) an inner shell with R_in~6.3E14cm, T_dust~500-105K, and Mdot~3E-5Msun/yr. An additional population of big dust grains (radius a>~0.4mm) with T_dust=150-20K and mass M_dust=(0.16-1.6)E-3 [D/3kpc]^2 Msun can account for the observed sub-mm and mm flux excess. The mass of the envelope enclosed within R_out=1.25E17cm derived from SED modeling is ~1[D/3kpc]^2 Msun.

We report high angular-resolution (~1") CO J=3--2 interferometric mapping, using the Submillimeter Array (SMA), of IRAS22036+5306 (I22036), a bipolar pre-planetary nebula (PPN) with knotty jets discovered in our HST SNAPshot survey of young PPNs. In addition, we have obtained supporting lower-resolution (~10") CO and 13CO J=1-0 observations with the Owens Valley Radio Observatory (OVRO) interferometer, as well as optical long-slit echelle spectra at the Palomar Observatory. The CO J=3-2 observations show the presence of a very fast (~220 km/s), highly collimated, massive (0.03 Msun) bipolar outflow with a very large scalar momentum (about 10^{39} g cm s^{-1}), and the characteristic spatio-kinematic structure of bow-shocks at the tips of this outflow. The Halpha line shows an absorption feature blue-shifted from the systemic velocity by ~100 km/s, which most likely arises in neutral interface material between the fast outflow and the dense walls of the bipolar lobes at low latitudes. The fast outflow in I22036, as in most PPNs, cannot be driven by radiation pressure. We find an unresolved source of submillimeter (and millimeter-wave) continuum emission in I22036, implying a very substantial mass (0.02-0.04 Msun) of large (radius >~1 mm), cold (< ~50 K) dust grains associated with I22036's toroidal waist. We also find that the 13C/12C ratio in I22036 is very high (0.16), close to the maximum value achieved in equilibrium CNO-nucleosynthesis (0.33). The combination of the high circumstellar mass (i.e., in the extended dust shell and the torus) and the high 13C/12C ratio in I22036 provides strong support for this object having evolved from a massive (>~4 Msun) progenitor in which hot-bottom-burning has occurred.

(abridged) The interaction of a fast wind with a spherical Asymptotic Giant Branch (AGB) wind is thought to be the basic mechanism for shaping Pre-Planetary Nebulae (PPN) and later Planetary Nebulae (PN). Due to the large speed of the fast wind, one expects extended X-ray emission from these objects, but X-ray emission has only been detected in a small fraction of PNs and only in one PPN. Using numerical simulations we investigate the constraints that can be set on the physical properties of the fast wind (speed, mass-flux, opening angle) in order to produce the observed X-ray emission properties of PPNs and PNs. We combine numerical hydrodynamical simulations including radiative cooling using the code FLASH with calculations of the X-ray properties of the resulting expanding hot bubble using the atomic database ATOMDB. In this first study, we compute X-ray fluxes and spectra using one-dimensional models. Comparing our results with analytical solutions, we find some agreements and many disagreements. In particular, we test the effect of different time histories of the fast wind on the X-ray emission and find that it is determined by the final stage of the time history during which the fast wind velocity has its largest value. The disagreements which are both qualitative and quantitative in nature argue for the necessity of using numerical simulations for understanding the X-ray properties of PNs.

We present new, sensitive, near-infrared images of the Class I protostar, Elias 29, in the Ophiuchus cloud core. To explore the relationship between the infall envelope and the outflow, narrowband H2 1-0 S(1), Br-gamma, and narrowband K-continuum filters were used to image the source with the Wide-Field Infrared Camera on the Hale 5m telescope and with Persson's Auxiliary Nasmyth Infrared Camera on the Baade 6.5 m telescope. The source appears as a bipolar, scattered light nebula, with a wide opening angle in all filters, as is typical for late-stage protostars. However, the pure H2 emission-line images point to the presence of a heretofore undetected precessing jet. It is argued that high-velocity, narrow, precessing jets provide the mechanism for creating the observed wide-angled outflow cavity in this source.

We present an atlas of Spitzer Space Telescope Infrared Spectrograph (IRS) spectra of highly luminous, compact mid-infrared sources in the Large Magellanic Cloud. Sources were selected on the basis of infrared colors and 8 micron (MSX) fluxes indicative of highly evolved, intermediate- to high-mass stars with current or recent mass loss at large rates. We determine the chemistry of the circumstellar envelope from the mid-IR continuum and spectral features and classify the spectral types of the stars. In the sample of 60 sources, we find 21 Red Supergiants (RSGs), 16 C-rich Asymptotic Giant Branch (AGB) stars, 11 HII regions, 4 likely O-rich AGB stars, 4 Galactic O-rich AGB stars, 2 OH/IR stars, and 2 B[e] supergiants with peculiar IR spectra. We find that the overwhelming majority of the sample AGB stars (with typical IR luminosities ~1.0E4 L_sun) have C-rich envelopes, while the O-rich objects are predominantly luminous RSGs with L_IR ~ 1.0E5 L_sun. We determine mean bolometric corrections to the stellar K-band flux densities and find that for carbon stars, the bolometric corrections depend on the infrared color, whereas for RSGs, the bolometric correction is independent of IR color. Our results reveal that objects previously classified as PNe on the basis of IR colors are in fact compact HII regions with very red IRS spectra that include strong atomic recombination lines and PAH emission features. We demonstrate that the IRS spectral classes in our sample separate clearly in infrared color-color diagrams that use combinations of 2MASS data and synthetic IRAC/MIPS fluxes derived from the IRS spectra. On this basis, we suggest diagnostics to identify and classify, with high confidence levels, IR-luminous evolved stars and HII regions in nearby galaxies using Spitzer and near-infrared photometry.

We describe Very Large Array observations of the extreme carbon star IRC+10216 at 8.4, 14.9, and 22.5 GHz made over a two year period. We find possible variability correlated with the infrared phase and a cm- to sub-millimeter wavelength spectral index very close to 2. The variability, observed flux densities, and upper limit on the size are consistent with the emission arising from the stellar photosphere or a slightly larger radio photosphere.

We have produced an atlas of Spitzer Infrared Spectrograph (IRS) spectra of mass-losing, evolved stars in the Large Magellanic Cloud. These stars were selected to have high mass-loss rates and so contribute significantly to the return of processed materials to the ISM. Our high-quality spectra enable the determination of the chemistry of the circumstellar envelope from the mid-IR spectral features and continuum. We have classified the spectral types of the stars and show that the spectral types separate clearly in infrared color-color diagrams constructed from 2MASS data and synthetic IRAC/MIPS fluxes derived from our IRS spectra. We present diagnostics to identify and classify evolved stars in nearby galaxies with high confidence levels using Spitzer and 2MASS photometry. Comparison of the spectral classes determined using IRS data with the IR types assigned based on NIR colors also revealed a significant number of misclassifications and enabled us to refine the NIR color criteria resulting in more accurate NIR color classifications of dust-enshrouded objects.

We present the spectral energy distribution (SED) of the pre-planetary nebula, IRAS 19475+3119 (I19475), from the optical to the far-infrared. We identify emission features due to crystalline silicates in the ISO SWS spectra of the star. We have fitted the SED of I19475 using a 1-D radiative transfer code, and find that a shell with inner and outer radii of 8.8X10^{16} and 4.4X10^{17}cm, and dust temperatures ranging from about 94K to 46K provide the best fit. The mass of this shell is greater than/equal to 1[34cm^{2}g^{-1}/kappa(100micron)][delta/200]M_Sun, where kappa(100micron) is the 100micron dust mass absorption coefficient (per unit dust mass), and delta is the gas-to-dust ratio. In agreement with results from optical imaging and millimeter-wave observations of CO emission of I19475, our model fits support an r^{-3} density law for its dust shell, with important implications for the interaction process between the fast collimated post-AGB winds and the dense AGB envelopes which results in the observed shapes of PPNs and PNs. We find that the observed JCMT flux at sub-millimeter wavelengths (850micron) is a factor ~ 2 larger than our model flux, suggesting the presence of large dust grains in the dust shell of I19475 which are not accounted for by our adopted standard MRN grain size distribution.

We report the first detection of X-ray emission in a pre-planetary nebula, Hen 3-1475. Pre-planetary nebulae are rare objects in the short transition stage between the Asymptotic Giant Branch and planetary nebula evolutionary phases, and Hen 3-1475, characterised by a remarkable S-shaped chain of optical knots, is one of the most noteworthy members of this class. Observations with the Advanced CCD Imaging Spectrometer (ACIS) onboard the Chandra X-Ray observatory show the presence of compact emission coincident with the brightest optical knot in this bipolar object, which is displaced from the central star by 2.7 arcsec along the polar axis. Model fits to the X-ray spectrum indicate an X-ray temperature and luminosity, respectively, of (4.3-5.7) 10^6 K and (4+/-1.4) 10^{31} (D/5 kpc)^2 erg s^{-1}, respectively. Our 3-sigma upper limit on the luminosity of compact X-ray emission from the central star in Hen 3-1475 is ~5 10^{31} (D/5 kpc)^2 erg s^{-1}. The detection of X-rays in Hen 3-1475 is consistent with models in which fast collimated post-AGB outflows are crucial to the shaping of planetary nebulae; we discuss such models in the context of our observations.

High-resolution near-infrared HST NICMOS (F160W, F222M) images and polarization (2 um) observations were made of four bipolar proto-planetary nebulae (PPNs): IRAS 17150-3224, IRAS 17441-2411, IRAS 17245-3951, and IRAS 16594-4656. The first three of these are viewed nearly edge-on, and for the first time the central stars in them are seen. Color maps reveal a reddened torus between the bipolar lobes in the edge-on cases, with bluer lobes. The polarization values are high, with maximum values ranging from 40 to 80%. The polarization patterns are basically centrosymmetric, with some deviations in the low polarization equatorial regions. For IRAS 17150-3224, circumstellar arcs are seen at 1.6 um, along with a newly-discovered loop in the equatorial region. Bright caps are seen at the end of the lobes, indicating that they are not open-ended. A distinct point-symmetric pattern is seen in the strengths of the polarization vectors, especially in IRAS 17150-3224. HST NICMOS observations provide a valuable complement to the WFPC2 visible images in deriving the basic structure of bipolar PPNs.

Using two-dimensional hydrodynamical simulations, we investigate the interaction of a collimated fast wind (CFW) interacting with a spherical asymptotic giant branch (AGB) wind as the mechanism for shaping proto-planetary nebulae and young planetary nebulae. In particular, we compare our simulations to the observations of an evolved PPN with multiple, highly collimated lobes, CRL 618.

In this paper (paper I) we present optical long-slit spectroscopy and imaging of the protoplanetary nebula CRL618. The optical lobes of CRL618 consist of shock-excited gas, which emits many recombination and forbidden lines, and dust, which scatters light from the innermost regions. From the analysis of the scattered Halpha emission, we derive a nebular inclination of i=24+-6 deg. The spectrum of the innermost part of the east lobe (visible as a bright, compact nebulosity close to the star in the Halpha HST image) is remarkably different from that of the shocked lobes but similar to that of the inner HII region, suggesting that this region represents the outermost parts of the latter. We find a non-linear radial variation of the gas velocity along the lobes. The largest projected LSR velocities (~80 km/s) are measured at the tips of the lobes, where the direct images show the presence of compact bow-shaped structures. The velocity of the shocks in CRL618 is in the range ~75-200 km/s, as derived from diagnostic line ratios and line profiles. We report a brightening (weakening) of [OIII]5007AA ([OI]6300AA) over the last ~10 years that may indicate a recent increase in the speed of the exciting shocks. From the analysis of the spatial variation of the nebular extinction, we find a large density contrast between the material inside the lobes and beyond them: the optical lobes seem to be `cavities' excavated in the AGB envelope by interaction with a more tenuous post-AGB wind. The electron density, with a mean value n_e~5E3-1E4 cm-3, shows significant fluctuations but no systematic decrease along the lobes, in agreement with most line emission arising in a thin shell of shocked material (the lobe walls) rather than in the post-AGB wind filling the interior of the lobes. (...)

We have used infrared polarimetric imaging with NICMOS to determine precisely the position of the star that illuminates (and presumably generated) the bipolar, pre-planetary reflection nebula RAFGL 2688 (the Egg Nebula). The polarimetric data pinpoint the illuminating star, which is not detected directly at wavelengths less than or equal to 2 microns, at a position well within the dark lane that bisects the nebula, 0.55" (about 550 AU) southwest of the infrared peak which was previously detected at the southern tip of the northern polar lobe. The inferred position of the central star corresponds to the geometric center of the tips of the four principle lobes of near-infrared H2 emission; identifying the central star at this position also reveals the strong point symmetric structure of the nebula, as seen both in the intensity and polarization structure of the polar lobes. The polarimetric and imaging data indicate that the infrared peak directly detected in the NICMOS images is a self-luminous source and, therefore, is most likely a distant binary companion to the illuminating star. Although present theory predicts that bipolar structure in pre-planetary and planetary nebulae is a consequence of binary star evolution, the separation between the components of the RAFGL 2688 binary system, as deduced from these observations, is much too large for the presence of the infrared companion to have influenced the structure of the RAFGL 2688 nebula.

We report NICMOS wide-band and polarimetric observations of four proto-planetary nebulae. Molecular hydrogen emission is detected near the ends of the bipolar lobes of IRAS 17150-3224, which is evidence for the interaction of a fast, collimated outflow with the remnant of the asymptotic giant branch wind.