The interaction of planetary nebulae and their AGB 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 ISM interaction strongly affects outer PN structures, with the dominant shaping occuring during the AGB phase. The simulations predict four stages of PN--ISM interaction whereby the PN is initially unaffected (1), then limb-brightened in the direction of motion (2), then distorted with the star moving away from the geometric centre (3) and finally so distorted that the object is no longer recognisable as a PN and may not be classed as such (4). Parsec-size shells around PN are predicted to be common. The structure and brightness of ancient PNe is 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 planetary nebulae in globular clusters.Comment: 22 pages, 16 figures, accepted by MNRAS (consists of 14 page journal paper and 8 page online-only appendix). Email C Wareing for high quality PDF versio

A discrete slug population model determined by egg production

Slugs are significant pests in agriculture (as well as a nuisance to gardeners), and it is therefore important to understand their population dynamics for the construction of efficient and effective control measures. Differential equation models of slug populations require the inclusion of large (variable) temporal delays, and strong seasonal forcing results in a non-autonomous system. This renders such models open to only a limited amount of rigorous analysis. In this paper, we derive a novel batch model based purely upon the quantity of eggs produced at different times of the year. This model is open to considerable reduction; from the resulting two variable discrete-time system it is possible to reconstruct the dynamics of the full population across the year and give conditions for extinction or global stability and persistence. Furthermore, the steady state temporal population distribution displays qualitatively different behavior with only small changes in the survival probability of slugs. The model demonstrates how small variations in the favorability of different years may result in widely different slug population fluctuations between consecutive years, and is in good agreement with field data

A model of Mira's cometary head/tail entering the Local Bubble

We model the cometary structure around Mira as the interaction of an AGB wind from Mira A, and a streaming environment. Our simulations introduce the following new element: we assume that after 200 kyr of evolution in a dense environment Mira entered the Local Bubble (low density coronal gas). As Mira enters the bubble, the head of the comet expands quite rapidly, while the tail remains well collimated for a 100 kyr timescale. The result is a broad-head/narrow-tail structure that resembles the observed morphology of Mira's comet. The simulations were carried out with our new adaptive grid code WALICXE, which is described in detail.Comment: 12 pages, 8 figures (4 in color). Accepted for publication in The Astrophysical Journa

A latitude-dependent wind model for Mira's cometary head

We present a 3D numerical simulation of the recently discovered cometary structure produced as Mira travels through the galactic ISM. In our simulation, we consider that Mira ejects a steady, latitude-dependent wind, which interacts with a homogeneous, streaming environment. The axisymmetry of the problem is broken by the lack of alignment between the direction of the relative motion of the environment and the polar axis of the latitude-dependent wind. With this model, we are able to produce a cometary head with a ``double bow shock'' which agrees well with the structure of the head of Mira's comet. We therefore conclude that a time-dependence in the ejected wind is not required for reproducing the observed double bow shock.Comment: 4 pages, 4 figures, accepted for publication in ApJ

Probing the Magnetized Interstellar Medium Surrounding the Planetary Nebula Sh 2-216

We present 1420 MHz polarization images of a 2.5 X 2.5 degree region around the planetary nebula (PN) Sh 2-216. The images are taken from the Canadian Galactic Plane Survey (CGPS). An arc of low polarized intensity appears prominently in the north-east portion of the visible disk of Sh 2-216, coincident with the optically identified interaction region between the PN and the interstellar medium (ISM). The arc contains structural variations down to the ~1 arcminute resolution limit in both polarized intensity and polarization angle. Several polarization-angle "knots" appear along the arc. By comparison of the polarization angles at the centers of the knots and the mean polarization angle outside Sh 2-216, we estimate the rotation measure (RM) through the knots to be -43 +/- 10 rad/m^2. Using this estimate for the RM and an estimate of the electron density in the shell of Sh 2-216, we derive a line-of-sight magnetic field in the interaction region of 5.0 +/- 2.0 microG. We believe it more likely the observed magnetic field is interstellar than stellar, though we cannot completely dismiss the latter possibility. We interpret our observations via a simple model which describes the ISM magnetic field around Sh 2-216, and comment on the potential use of old PNe as probes of the magnetized ISM.Comment: 25 pages, 4 figures. Accepted for publication in the Astrophysical Journa

A new mechanical stellar wind feedback model for the Rosette Nebula

The famous Rosette Nebula has an evacuated central cavity formed from the stellar winds ejected from the 2–6 Myr old codistant and comoving central star cluster NGC 2244. However, with upper age estimates of less than 110 000 yr, the central cavity is too young compared to NGC 2244 and existing models do not reproduce its properties. A new proper motion study herein using Gaia data reveals the ejection of the most massive star in the Rosette, HD 46223, from NGC 2244 occurred 1.73 (+0.34, −0.25) Myr (1σ uncertainty) in the past. Assuming this ejection was at the birth of the most massive stars in NGC 2244, including the dominant centrally positioned HD 46150, the age is set for the famous ionized region at more than 10 times that derived for the cavity. Here, we are able to reproduce the structure of the Rosette Nebula, through simulation of mechanical stellar feedback from a 40 Mₒ star in a thin sheet-like molecular cloud. We form the 135 000 Mₒ cloud from thermally unstable diffuse interstellar medium (ISM) under the influence of a realistic background magnetic field with thermal/magnetic pressure equilibrium. Properties derived from a snapshot of the simulation at 1.5 Myr, including cavity size, stellar age, magnetic field, and resulting inclination to the line of sight, match those derived from observations. An elegant explanation is thus provided for the stark contrast in age estimates based on realistic diffuse ISM properties, molecular cloud formation and stellar wind feedback

The Shapes of AGB Envelopes as Probes of Binary Companions

We describe how the large scale geometry of the circumstellar envelopes of asymptotic giant branch stars can be used to probe the presence of unseen stellar companions. A nearby companion modifies the mass loss by gravitationally focusing the wind towards the orbital plane, and thereby determines the shape of the envelope at large distances from the star. Using available simulations, we develop a prescription for the observed shapes of envelopes in terms of the binary parameters, envelope orientation, and type of observation. The prescription provides a tool for the analysis of envelope images at optical, infrared, and millimetre wavelengths, which can be used to constrain the presence of companions in well observed cases. We illustrate this approach by examining the possible role of binary companions in triggering the onset of axi-symmetry in planetary nebula formation. If interaction with the primary leads to axi-symmetry, the spherical halos widely seen around newly formed nebulae set limits on the companion mass. Only low mass objects may orbit close to the primary without observable shaping effects: they remain invisible until the interaction causes a sudden change in the mass loss geometry.Comment: 11 pages, 7 figures, to appear in MNRA