2018MNRAS.480...75N

To understand how the circumstellar environments of post-Active Giant Branch (AGB) stars develop into planetary nebulae, we initiate a systematic study of 2D axisymmetric hydrodynamic simulations of protoplanetary nebula (pPN) with a modified ZEUS code. The aim of this first work is to compare the structure of prolate ellipsoidal winds into a stationary ambient medium where both media can be either atomic or molecular. We specifically model the early twin-shock phase which generates a decelerating shell. A thick deformed and turbulent shell grows when an atomic wind expands into an atomic medium. In all other cases, the interaction shell region fragments into radial protrusions due to molecular cooling and chemistry. The resulting fingers eliminate any global slip parallel to the shell surface. This rough surface implies that weak shocks are prominent in the excitation of the gas despite the fast speed of advance. This may explain why low-excitation molecular hydrogen is found towards the front of elliptical pPN. We constrain molecular dissociative fractions and time-scales of fast H₂ winds and the pPN lifetime with wind densities ∼10⁵ cm^–3 and shock speeds of 80∼200 km s^–1. We identify a variety of stages associated with thermal excitation of H₂ near-infrared emission. Generated line emission maps and position-velocity diagrams enable a comparison and distinction with post-AGB survey results. The 1→0 S(1)&2→1 S(1) lines are lobe-dominated bows rather than bipolar shells.