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.     

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.     

Asymmetry and inhomogeneity in proto- and young planetary nebulae

I study some effects of aspherical mass loss during the last stages of the asymptotic giant branch (AGB) on the appearance of proto-planetary nebulae (proto-PNs) and young PNs. The aspherical mass loss can be small-scale inhomogeneities, and/or axially symmetric mass-loss geometry. I first examine the role of the dust opacity in the optical band on the appearance of proto-PNs. I conclude that large optical depths will be found in proto-PNs that are post-AGB stars having high equatorial mass-loss rates, which require a stellar binary companion for their existence. In these cases light from the central star will reach larger distances along and near the polar directions, leading to the appearance of an elongated reflection nebula. These proto-PNs will become bipolar PNs, i.e., PNs with two lobes and an equatorial waist between them, or extreme ellipticals, e.g., a ring but no lobes on the two sides of the equatorial plane. I then derive the conditions for the enhancement of non-radial density inhomogeneity by the propagation of the ionization front at the early PN stages. The ionization will proceed faster in the radial direction along low-density regions. The low-density regions will be heated earlier, and they will expand as a result of their higher pressures, reducing further their densities. The opposite occurs in high-density regions. The condition for this ionization instability to develop is that the ionization time difference between two directions at the same radius is longer than the sound crossing time between these two regions. This condition for the ionization front instability can be expressed as a condition on the mass-loss rate inhomogeneity, i.e., its dependence on direction.     

Far-ultraviolet spectroscopy of the hot DA white dwarf WD 2218+706 (DeHt5) with STIS

We report a study of the photospheric composition of the hot DA white dwarf WD 2218+706, which is also the central star of the old planetary nebula DeHt5. Helium is detected in the far-UV spectrum. In addition, the star clearly contains significant quantities of elements heavier than He at abundances generally a factor of 2 to 10 higher than those found in the archetypal heavy element-rich DA G191−B2B. This is the first detection of trace He using the He  ii λ 1640 line in an isolated DA white dwarf, but the low surface gravity is more indicative of a binary evolution route from the red giant branch rather than a path along the asymptotic giant branch (AGB) as a single star. However, the absence of any evidence for a companion star and the uncertainty in the measured mass for WD 2218+706 still allow the possibility of an origin along an AGB evolutionary track.
We reanalyse the existing optical spectra of WD 2218+706 using our latest pure H and heavy element-rich model atmospheres, obtaining a good match between the observed and synthetic spectra with either set of models. We find little evidence of any inconsistency in the temperature required to fit individual Balmer lines, as reported elsewhere for this star. Any discrepancies we see are confined to the H α line and the core of H β but they do not compromise our analysis.     

CCD images and long-slit spectroscopy of the ring nebula around θ Mus

We present the first digital CCD images and long-slit spectroscopy of the optical ring nebula around the Wolf–Rayet star θ Mus. The CCD images obtained through narrow-band filters centred at [O  iii ] and Hα show that the nebula has a filamentary structure, similar to supernova remnants, mainly seen in [O  iii ]. A spatial detachment between [O  iii ] and Hα images suggests excitation stratification, or multiple rings. An analysis of the physical conditions in the nebula was performed by means of long-slit CCD spectra. The spectral images show that the nebula is of low density and medium excitation. By means of quotients of recombination and collisional spectral line fluxes we determine that the principal excitation mechanism is photoionization. We have determined the electronic temperature and density, and chemical abundances for the oxygen at different sites within the nebula. Nebular chemical abundances are found to be similar to the Galactic ISM, indicating that the nebula is mainly composed of swept up material.     

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.     

A large radio nebula around P Cygni

We present a large set of radio observations of the luminous blue variable P Cygni. These include two 6-cm images obtained with MERLIN which spatially resolve the 6-cm photosphere, monitoring observations obtained at Jodrell Bank every few days over a period of two months, and VLA observations obtained every month for seven years. This combination of data shows that the circumstellar environment of P Cyg is highly inhomogeneous, that there is a radio nebula extending to almost an arcminute from the star at 2 and 6 cm, and that the radio emission is variable on a time-scale no longer than one month, and probably as short as a few days. This short-time-scale variability is difficult to explain. We present a model for the radio emission with which we demonstrate that the star has probably been losing mass at a significant rate for at least a few thousand years, and that it has undergone at least two major outbursts of increased mass loss during the past two millenia.     

Modeling the Dust Shell of V4334 Sgr

The central star V4334 Sgr (Sakurai's Object) of the planetary nebula PN G010.4+04.4 underwent in 1995–1996 the rare event of a very late helium flash.It is only one of two such events during the era of modern astronomy (the other event was V605 Aql = Nova Aql 1919). All other prominentobjects of that type originate from events several thousands of years ago (e.g. A30, A78). Hence, only snapshots can be modeled for those objects.V4334 Sgr allows for the first time a dynamic consideration of the formation of the dust shell from the beginning.We present here a model which is able to describe the complete photometric behaviorof the object, including the fine structure dips of the optical light curve during the first two years of the mass loss and the dust formation.     

Backflow in post-asymptotic giant branch stars

We derive the conditions for a backflow toward the central star(s) of circumstellar material to occur during the post-asymptotic giant branch (post-AGB) phase. The backflowing material may be accreted by the post-AGB star and/or its companion, if such exists. Such a backflow may play a significant role in shaping the descendant planetary nebula, by, among other things, slowing down the post-AGB evolution, and by forming an accretion disc which may blow two jets. We consider three forces acting on a slowly moving mass element: the gravity of the central system, radiation pressure, and fast wind ram pressure. We find that for a significant backflow to occur, a slow dense flow should exist, such that the relation between the total mass in the slow flow, M _i , and the solid angle it covers Ω, is given by     , where     . The requirement for both a high mass-loss rate per unit solid angle and a very slow wind, such that it can be decelerated and flow back, probably requires close binary interaction, hence this process is rare.     

The blue supergiant Sher 25 and its intriguing hourglass nebula

The blue supergiant Sher 25 is surrounded by an asymmetric, hourglass-shaped circumstellar nebula. Its structure and dynamics have been studied previously through high-resolution imaging and spectroscopy, and it appears dynamically similar to the ring structure around SN 1987A. Here, we present long-slit spectroscopy of the circumstellar nebula around Sher 25, and of the background nebula of the host cluster NGC 3603. We perform a detailed nebular abundance analysis to measure the gas-phase abundances of oxygen, nitrogen, sulphur, neon and argon. The oxygen abundance in the circumstellar nebula  (12 + log O/H = 8.61 ± 0.13 dex)  is similar to that in the background nebula (8.56 ± 0.07), suggesting that the composition of the host cluster is around solar. However, we confirm that the circumstellar nebula is very rich in nitrogen, with an abundance of 8.91 ± 0.15, compared to the background value of 7.47 ± 0.18. A new analysis of the stellar spectrum with the fastwind model atmosphere code suggests that the photospheric nitrogen and oxygen abundances in Sher 25 are consistent with the nebular results. While the nitrogen abundances are high, when compared to stellar evolutionary models, they do not unambiguously confirm that the star has undergone convective dredge-up during a previous red supergiant phase. We suggest that the more likely scenario is that the nebula was ejected from the star while it was in the blue supergiant phase. The star's initial mass was around  50 M_⊙  , which is rather too high for it to have had a convective envelope stage as a red supergiant. Rotating stellar models that lead to mixing of core-processed material to the stellar surface during core H-burning can quantitatively match the stellar results with the nebula abundances.     

The evolution of a rapidly accreting helium white dwarf to become a low-luminosity helium star

We have examined the evolution of merged low-mass double white dwarfs which become low-luminosity (or high-gravity) extreme helium stars. We have approximated the merging process by the rapid accretion of matter, consisting mostly of helium, on to a helium white dwarf. After a certain mass is accumulated, a helium shell flash occurs, the radius and luminosity increase and the star becomes a yellow giant. Mass accretion is stopped artificially when the total mass reaches a pre-determined value. As the helium-burning shell moves inwards with repeating shell flashes, the effective temperature gradually increases as the star evolves towards the helium main sequence. When the mass interior to the helium‐burning shell is approximately 0.25 M_⊙, the star enters a regime where it is pulsationally unstable. We have obtained radial pulsation periods for these models.
These models have properties very similar to those of the pulsating helium star V652 Her. We have compared the rate of period change of the theoretical models with that observed in V652 Her, as well as with its position on the Hertzsprung–Russell diagram. We conclude that the merger between two helium white dwarfs can produce a star with properties remarkably similar to those observed in at least one extreme helium star, and is a viable model for their evolutionary origin. Such helium stars will evolve to become hot subdwarfs close to the helium main sequence. We also discuss the number of low-luminosity helium stars in the Galaxy expected for our evolution scenario.     

Light element abundances in carbon-enhanced metal-poor stars

We model the evolution of the abundances of light elements in carbon-enhanced metal-poor (CEMP) stars, under the assumption that such stars are formed by mass transfer in a binary system. We have modelled the accretion of material ejected by an asymptotic giant branch star on to the surface of a companion star. We then examine three different scenarios: one in which the material is mixed only by convective processes, one in which thermohaline mixing is present and a third in which both thermohaline mixing and gravitational settling are taken in to account. The results of these runs are compared to light element abundance measurements in CEMP stars (primarily CEMP- s stars, which are rich in s -processes elements and likely to have formed by mass transfer from an AGB star), focusing on the elements Li, F, Na and Mg. None of the elements is able to provide a conclusive picture of the extent of mixing of accreted material. We confirm that lithium can only be preserved if little mixing takes place. The bulk of the sodium observations suggest that accreted material is effectively mixed but there are also several highly Na and Mg-rich objects that can only be explained if the accreted material is unmixed. We suggest that the available sodium data may hint that extra mixing is taking place on the giant branch, though we caution that the data are sparse.     

X-ray binaries

Summary The various types and classes of X-ray binary are reviewed high-lighting recent results. The high mass X-ray binaries (HMXRBs) can be used to probe the nature of the mass loss from the OB star in these systems. Absorption measurements through one orbital cycle of the supergiant system X1700-37 are well modelled by a radiation driven wind and also require a gas stream trailing behind the X-ray source. In Cen X-3 the gas stream is accreted by the X-ray source via an accretion disk. Changes in the gas stream can cause the disk to thicken and the disk to obscure the X-ray source. How close the supergiant is to corotation seems to be as much a critical factor in these systems as how close it is to filling its Roche lobe. In the Be star X-ray binaries a strong correlation between the neutron stars rotation period and its orbital period has been explained as due to the neutron star being immersed in a dense, slow moving equatorial wind from the Be star. For the X-ray pulsars in the transient Be X-ray binaries a centrifugal barrier to accretion is important in determining the X-ray lightcurve and the spin evolution. The X-ray orbital modulations from the low mass X-ray binaries, LMXRBs, include eclipses by the companion and/or periodic dipping behaviour from structure at the edge of the disk. The corresponding optical modulations show a smooth sinusoidal like component and in some cases a sharp eclipse by the companion. The orbital period of the LMXRB XB1916-05 is 1% longer in the optical compared to that given by the X-ray dip period. The optical period has been interpreted as the orbital period, but this seems inconsistent with the well established view of the origin of the X-ray modulations in LMXRB. A new model is presented that assumes the X-ray dip period is the true orbital period. The 5.2 h eclipsing LMXRB XB2129+47 recently entered a low state and optical observations unexpectedly reveal an F star which is too big to fit into the binary. This is probably the first direct evidence that an X-ray binary is part of a hierarchical triple. Finally the class of X-ray binaries containing black hole candidates is reviewed focusing on the value of using X-ray signatures to identify new candidates.     

Composite spectra XIX: HD 208253, a long‐period binary whose hot secondary has enhanced strontium and barium

We separate and analyse the component spectra of the composite‐spectrum binary HD 208253. We find that the cool primary is an evolving star of spectral type G7 III, while its hot secondary is an early‐A dwarf. The giant is currently near the lowest point of the red‐giant branch and is slightly less luminous than its dwarf companion. We provide a set of precise radial‐velocity measurements for both stars. The double‐lined orbit which we derive from them shows that the component mass ratio is close to unity (q = 1.05 ± 0.01). We deduce the physical properties of both stars, determine their respective masses to be 2.75 ± 0.07 Me (giant) and 2.62 ± 0.07 Me (dwarf), and show that the orbit's inclination is within a degree or two of 68°. The spectrum of the A‐type component has quite component has quite narrow lines (we infer a rotational velocity of 18 km s^–1), though since the period of the orbit is well over 1 year that component cannot be in synchronous rotation. An intriguing property of the dwarf is its enhanced Sr and Ba, though it does not exhibit the other spectral peculiarities that would signal a classical Am star. While by no means unique amongst the multitude of oddities exhibited by A and early‐F stars, this dwarf which we have uncovered in a long‐period binary offers valuable constraints and challenges to stellar‐evolution theory. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The eccentricities of the barium stars

We investigate the eccentricities of barium (Ba  ii ) stars formed via a stellar wind accretion model. We carry out a series of Monte Carlo simulations using a rapid binary evolution algorithm, which incorporates full tidal evolution, mass loss and accretion, and nucleosynthesis and dredge-up on the thermally pulsing asymptotic giant branch. We follow the enhancement of barium in the envelope of the accreting main-sequence companion and dilution into its convective envelope once the star ascends the giant branch.
The observed eccentricities of Ba  ii stars are significantly smaller than those of an equivalent set of normal red giants but are nevertheless non-zero. We show that such a distribution of eccentricities is consistent with a wind accretion model for Ba  ii star production with weak viscous tidal dissipation in the convective envelopes of giant stars. We successfully model the distribution of orbital periods and the number of observed Ba  ii stars. The actual distribution of eccentricities is quite sensitive to the strength of the tides, so that we are able to confirm that this strength is close to, but less than, what is expected theoretically and found with alternative observational tests. Two systems – one very short-period but eccentric, and one long-period and highly eccentric – still lie outside the envelope of our models, and so require a more exotic formation mechanism. All our models, even those which were a good fit to the observed distributions, overproduced the number of high-period barium stars, a problem that could not be solved by some combination of the three parameters: tidal strength, tidal enhancement and wind accretion efficiency.     

IRAS 01005＋7910：a High Galactic Latitude Post—AGB Star

IRAS 01005＋7910 is a cold IRAS source. We present its optical identification, photometric and spectroscopic observation results. Its optical counterpart is classified as a B2Ie star with V magnitude 10.85. Its Hα line shows the P Cygni profile.According to its location in the Galaxy (b=16.6), we consider it to be a post-AGB star or a proto-planetary nebula.     

A unique Galactic planetary nebula with a [WN] central star

We report the discovery of the first probable Galactic [WN] central star of a planetary nebula (CSPN). The planetary nebula candidate was found during our systematic scans of the AAO/UKST Hα Survey of the Milky Way. Subsequent confirmatory spectroscopy of the nebula and central star reveals the remarkable nature of this object. The nebular spectrum shows emission lines with large expansion velocities exceeding 150 km s⁻¹, suggesting that perhaps the object is not a conventional planetary nebula. The central star itself is very red and is identified as being of the [WN] class, which makes it unique in the Galaxy. A large body of supplementary observational data supports the hypothesis that this object is indeed a planetary nebula and not a Population I Wolf–Rayet star with a ring nebula.     

IPHAS discoveries of young stars towards Cyg OB2 and its southern periphery

We report on the discovery of over 50 strong Hα emitting objects towards the large OB association Cyg OB2 and the H  ii region DR 15 on its southern periphery. This was achieved using the INT Photometric Hα Survey of the Northern Galactic Plane (IPHAS), combined with follow-up spectroscopy using the MMT multi-object spectrometer HectoSpec. We present optical spectra, supplemented with optical r ',  i ' and H α photometry from IPHAS, and near-infrared J ,  H and K photometry from Two Micron All Sky Survey. The position of the objects in the ( J − H ) versus ( H − K ) diagram strongly suggests most of them are young. Many show Ca  ii infrared triplet emission indicating that they are in a pre-main-sequence phase of evolution of T Tauri and Herbig Ae nature. Among these, we have uncovered pronounced clustering of T Tauri stars roughly a degree south of the centre of Cyg OB2, in an arc close to the H  ii region DR 15, and the radio ring nebula G79.29+0.46, for which we discuss its candidacy as a luminous blue variable. The emission-line objects towards Cyg OB2 itself could be the brightest most prominent component of a population of lower mass pre-main-sequence stars that has yet to be uncovered. Finally, we discuss the nature of the ongoing star formation in Cyg OB2 and the possibility that the central OB stars have triggered star formation in the periphery.     

The magnetosphere of oscillating neutron stars in general relativity

Just as a rotating magnetized neutron star has material pulled away from its surface to populate a magnetosphere, a similar process can occur as a result of neutron-star pulsations rather than rotation. This is of interest in connection with the overall study of neutron star oscillation modes but with a particular focus on the situation for magnetars. Following a previous Newtonian analysis of the production of a force-free magnetosphere in this way Timokhin et al., we present here a corresponding general-relativistic analysis. We give a derivation of the general relativistic Maxwell equations for small-amplitude arbitrary oscillations of a non-rotating neutron star with a generic magnetic field and show that these can be solved analytically under the assumption of low current density in the magnetosphere. We apply our formalism to toroidal oscillations of a neutron star with a dipole magnetic field and find that the low current density approximation is valid for at least half of the oscillation modes, similarly to the Newtonian case. Using an improved formula for the determination of the last closed field line, we calculate the energy losses resulting from toroidal stellar oscillations for all of the modes for which the size of the polar cap is small. We find that general relativistic effects lead to shrinking of the size of the polar cap and an increase in the energy density of the outflowing plasma. These effects act in opposite directions but the net result is that the energy loss from the neutron star is significantly smaller than suggested by the Newtonian treatment.