The interaction of planetary nebulae and their asymptotic giant branch progenitors with the interstellar medium

Interaction with the interstellar medium (ISM) cannot be ignored in understanding planetary nebula (PN) evolution and shaping. In an effort to understand the range of shapes observed in the outer envelopes of PNe, we have run a comprehensive set of three-dimensional hydrodynamic simulations, from the beginning of the asymptotic giant branch (AGB) superwind phase until the end of the post-AGB/PN phase. A 'triple-wind' model is used, including a slow AGB wind, fast post-AGB wind and third wind reflecting the linear movement through the ISM. A wide range of stellar velocities, mass-loss rates and ISM densities have been considered.
We find that ISM interaction strongly affects outer PN structures, with the dominant shaping occurring during the AGB phase. The simulations predict four stages of PN–ISM interaction whereby (i) the PN is initially unaffected, (ii) then limb-brightened in the direction of motion, (iii) then distorted with the star moving away from the geometric centre, and (iv) finally so distorted that the object is no longer recognizable as a PN and may not be classed as such. Parsec-size shells around PNe are predicted to be common. The structure and brightness of ancient PNe are largely determined by the ISM interaction, caused by rebrightening during the second stage; this effect may address the current discrepancies in Galactic PN abundance. The majority of PNe will have tail structures. Evidence for strong interaction is found for all known PNe in globular clusters.     

Why every bipolar planetary nebula is 'unique'

We present the many evolutionary routes that progenitors of bipolar planetary nebulae (BPNe) can take. Overall, there are about a hundred qualitatively different evolutionary routes, hence about a hundred qualitatively different types of BPNe. Within each type there are quantitative differences as well. Adding the dependence of the appearance on inclination, we find that the number of different apparent structures of BPNe is about equal to, or even larger than, the number of known BPNe and proto-BPNe. Accordingly we argue that every BPN is a 'unique' object in its appearance, but all can be explained within the binary model paradigm. Therefore, we request a stop to the attaching of adjectives such as 'unique', 'peculiar', and 'unusual' to BPNe and proto-BPNe, thereby removing the need to invoke a new model for almost every 'unusual' BPN. As a case study we try to build a binary model for the proto-BPN OH 231.8+4.2. In our preferred model the AGB Mira-type star has a main sequence companion of mass ∼1 M_⊙, orbital period of ∼5 yr, and eccentricity of ≳0.1.     

Morphological and kinematic signatures of a binary central star in the planetary nebula Hu 2-1

We present H α , [N  ii ] and [O  iii ] ground-based and HST archive images, VLA–A 3.6-cm continuum and H92 α emission-line data and high-resolution long-slit [N  ii ] spectra of the planetary nebula Hu 2-1. A large number of structural components are identified in the nebula: an outer bipolar and an inner shell, two pairs of collimated bipolar structures at different directions, monopolar bow-shock-like structures, and an extended equatorial structure within a halo. The formation of Hu 2-1 appears to be dominated by anisotropic mass ejection during the late-AGB stage of the progenitor and by variable, 'precessing' collimated bipolar outflows during the protoplanetary nebula and/or early planetary nebula phases. Different observational results strongly support the existence of a binary central star in Hu 2-1, among them (1) the observed point-symmetry of the bipolar lobes and inner shell, and the departures from axial symmetry of the bipolar lobes, (2) the off-centre position of the central star, (3) the detection of mass ejection towards the equatorial plane, and (4) the presence of 'precessing' collimated outflows. In addition, (5) an analysis of the kinematics shows that the systemic velocity of the bipolar outflows does not coincide with the systemic velocity of the bipolar shell. We propose that this velocity difference is a direct evidence of orbital motion of the ejection source in a binary system. From a deduced orbital velocity of ∼10 km s⁻¹, a semimajor axis of ∼ 9–27 au and period of ∼ 25–80 yr are obtained, assuming a reasonable range of masses. These parameters are used to analyse the formation of Hu 2-1 within current scenarios of planetary nebulae with binary central stars.     

Classification of planetary nebulae by their departure from axisymmetry

We propose a scheme to classify planetary nebulae (PNe) according to their departure from axisymmetric structure. We consider only departure along and near the equatorial plane, i.e. between the two sides perpendicular to the symmetry axis of the nebula. We consider six types of departure from axisymmetry: (1) PNe where the central star is not at the centre of the nebula; (2) PNe having one side brighter than the other; (3) PNe having unequal size or shape of the two sides; (4) PNe where the symmetry axis is bent, e.g. the two lobes in a bipolar PN are bent toward the same side; (5) PNe where the main departure from axisymmetry is in the outer regions, e.g. an outer arc; and (6) PNe that show no departure from axisymmetry, i.e. any departure, if it exists, is on scales smaller than the scale of blobs, filaments and other irregularities in the nebula. PNe that possess more than one type of departure are classified by the most prominent type. We discuss the connection between departure types and the physical mechanisms that may cause them, mainly resulting from the influence of a stellar binary companion. We find that ∼50 per cent of all PNe in the analysed sample possess large-scale departure from axisymmetry. This number is larger than that expected from the influence of binary companions, namely ∼25–30 per cent. We argue that this discrepancy comes from many PNe where the departure from axisymmetry, mainly unequal size, shape or intensity, results from the presence of long-lived and large (hot or cool) spots on the surface of their asymptotic giant branch progenitors. Such spots locally enhance the mass-loss rate, leading to a departure from axisymmetry, mainly near the equator, in the descendent PN.     

High-speed knots in the hourglass-shaped planetary nebula Hubble 12

We present a detailed kinematical analysis of the young compact hourglass-shaped planetary nebula Hb 12. We performed optical imaging and long-slit spectroscopy of Hb 12 using the Manchester echelle spectrometer with the 2.1-m San Pedro Mártir telescope. We reveal, for the first time, the presence of end caps (or knots) aligned with the bipolar lobes of the planetary nebula shell in a deep [N  ii ]λ6584 image of Hb 12. We measured from our spectroscopy radial velocities of  ∼120 km s⁻¹  for these knots.
We have derived the inclination angle of the hourglass-shaped nebular shell to be ∼65° to the line of sight. It has been suggested that Hb 12's central star system is an eclipsing binary which would imply a binary inclination of at least 80°. However, if the central binary has been the major shaping influence on the nebula, then both nebula and binary would be expected to share a common inclination angle.
Finally, we report the discovery of high-velocity knots with Hubble-type velocities, close to the core of Hb 12, observed in Hα and oriented in the same direction as the end caps. Very different velocities and kinematical ages were calculated for the outer and inner knots showing that they may originate from different outburst events.     

The birth of a planetary nebula around the carbon star IRC+10216

Our current understanding of the evolution of solar-type stars suggests that after a period as a red giant star, during which mass loss occurs continuously in the form of a stellar wind, a period of intense mass loss known as a superwind occurs, during which a significant fraction of the envelope of the star is ejected into space, forming the material from which a planetary nebula (PN) will be constructed. It has been suggested that this superwind ejects material from the star in a toroidal or disc-like fashion, rather than isotropically. Here we present Hubble Space Telescope optical images of a toroidal superwind caught in the act: our images of the carbon star IRC+10216, which is believed to be in the final stages of red giant evolution, show that most of its optical emission is a bipolar reflection nebula. We show that the full spectral energy distribution and these images can be modelled as an equatorially enhanced dusty superwind, providing the first direct observational support for the toroidal superwind model, and supporting the 'interacting winds' model of PN formation.     

Radially aligned clumps and tails in planetary nebulae

We study the formation of radially aligned condensations and tails through the compression of material inside ionization shadows at early ionization phases of planetary nebulae. A dense clump, formed before ionization starts, forms an ionization shadow behind it. The surroundings, which are ionized before the shadow, have a higher temperature, and as a result compress the material in the shadow, forming a compressed tail. If the compressed tail crosses a dense shell, a dense condensation (clump) is formed there. At later stages this condensation is ionized and observed as a bright knot, radially aligned with the inner clump. We find that for the shadow to be effective, the clump should be already present as the ionization by the central star starts, and its density enhancement should be by a factor of ≳ 5. We propose this mechanism as an explanation for the radially aligned condensations recently found in the planetary nebula IC 4593.     

Isolated versus common envelope dynamos in planetary nebula progenitors

The origin, evolution and role of magnetic fields in the production and shaping of proto-planetary nebulae (PPNe) and planetary nebulae (PNe) are a subject of active research. Most PNe and PPNe are axisymmetric with many exhibiting highly collimated outflows; however, it is important to understand whether such structures can be generated by isolated stars or require the presence of a binary companion. Towards this end, we study a dynamical, large-scale α−Ω interface dynamo operating in a 3.0 M_⊙ Asymptotic Giant Branch (AGB) star in both an isolated setting and a setting in which a low-mass companion is embedded inside the envelope. The back reaction of the fields on the shear is included and differential rotation and rotation deplete via turbulent dissipation and Poynting flux. For the isolated star, the shear must be resupplied in order to sufficiently sustain the dynamo. Furthermore, we investigate the energy requirements that convection must satisfy to accomplish this by analogy to the Sun. For the common envelope case, a robust dynamo results, unbinding the envelope under a range of conditions. Two qualitatively different types of explosion may arise: (i) magnetically induced, possibly resulting in collimated bipolar outflows and (ii) thermally induced from turbulent dissipation, possibly resulting in quasi-spherical outflows. A range of models is presented for a variety of companion masses.     

The enigma of the oldest 'nova': the central star and nebula of CK Vul

CK Vul is classified as, amongst others, the slowest known nova, a hibernating nova or a very late thermal pulse object. Following its eruption in ad 1670, the star remained visible for 2 yr. A 15-arcsec nebula was discovered in the 1980s, but the star itself has not been detected since the eruption. We here present radio images which reveal a 0.1-arcsec radio source with a flux of 1.5 mJy at 5 GHz. Deep Hα images show a bipolar nebula with a longest extension of 70 arcsec, with the previously known compact nebula at its waist. The emission-line ratios show that the gas is shock-ionized, at velocities  >100 km s⁻¹  . Dust emission yields an envelope mass of  ∼5 × 10⁻² M_⊙  . Echelle spectra indicate outflow velocities up to 360 km s⁻¹. From a comparison of images obtained in 1991 and 2004 we find evidence for expansion of the nebula, consistent with an origin in the 1670 explosion; the measured expansion is centred on the radio source. No optical or infrared counterpart is found at the position of the radio source. The radio emission is interpreted as thermal free–free emission from gas with   T _e∼ 10⁴ K  . The radio source may be due to a remnant circumbinary disc, similar to those seen in some binary post-AGB stars. We discuss possible classifications of this unique outburst, including that of a sub-Chandrasekhar mass supernova, a nova eruption on a cool, low-mass white dwarf or a thermal pulse induced by accretion from a circumbinary disc.     

The optical spectrum of the planetary nebula NGC 6543

With the Hamilton echelle spectrograph at the Lick Observatory, emission-rich spectral lines of the planetary nebula NGC 6543 were secured in the wavelength range from 3550 to 10 100 Å. We chose two bright regions, ∼8 arcsec east and ∼13 arcsec north of the central star, the physical conditions and chemical abundances of which may differ as a result of the different physical characteristics involving the mass ejection of different epochs. By combining Hamilton echelle observations with archive UV data secured with the International Ultraviolet Explorer ( IUE ), we obtain improved diagnostics and chemical compositions for the two observed regions. The diagnostic diagram gives the average value of T _e=8000∼8300 K, and the electron number density near N _e∼5000 cm⁻³ for most ions, while some low-excitation lines indicate much higher temperatures, i.e. T _e∼10 000 K. With the construction of a photoionization model, we try to fit the observed spectra in a self-consistent way: thus, for most elements, we employ the same chemical abundances in the nebular shell; and we adopt an improved Sobolev approximation model atmosphere for the hydrogen-deficient Wolf–Rayet type central star. Within the observational errors, the chemical abundances do not seem to show any positional variation except for helium. The chemical abundances of NGC 6543 appear to be the same as in average planetary nebulae. The progenitor star may have been an object of one solar mass, most of the heavier elements of which were less plentiful than in the Sun.     

The determination of progenitor masses for planetary nebulae

We have evaluated the likely progenitor masses M _PG of nebulae having elliptical, circular and bipolar morphologies, using observed ratios between the populations of these sources, and deduced central star mass functions. We find that most bipolar nebulae (BPNe) are likely to arise from progenitors having mass M _PG>2.3 M_⊙ and spectral types earlier than A3.2, whilst circular sources are associated with progenitors of mass 1.0 M_⊙< M _PG<1.2 M_⊙ and spectral range G1.9–F7.8 . Elliptical sources arise from intermediate-mass progenitors. The procedures employed to determine these values are relatively insensitive to uncertainties in scaleheights and population ratios, and completely insensitive to uncertainties in the distance scale. They are, however, dependent upon the precise forms adopted for the initial–final and central star mass functions, and we discuss the sensitivity of M _PG to uncertainties in these functions.     

Radio observations of the planetary nebula around the OH/IR star OH 354.88-0.54 (V1018 Sco)

We present radio observations of the unique, recently formed, planetary nebula (PN) associated with a very long-period OH/IR variable star V1018 Sco that is unequivocally still in its asymptotic giant branch phase. Two regions within the optical nebula are clearly detected in non-thermal radio continuum emission, with radio spectral indices comparable to those seen in colliding-wind Wolf–Rayet binaries. We suggest that these represent shocked interactions between the hot, fast stellar wind and the cold nebular shell that represents the PN's slow wind moving away from the central star. This same interface produces both synchrotron radio continuum and the optical PN emission. The fast wind is neither spherical in geometry nor aligned with any obvious optical or radio axis. We also report the detection of transient H₂O maser emission in this nebula.