Abstract. The shaping of Pre-planetary nebulae (PPN) and Planetary Nebulae (PN) is believed to result from the interaction of a fast, collimated post-AGB wind (CFW) plowing into the slow, dense wind emitted during the AGB phase, followed by an isotropic tenuous wind during the PN phase (Sahai & Trauger 1998). The expanding shell which forms the PPN and later the PN is expected to produce extended X-ray emission due to the large speed of the fast wind and resulting high temperatures in the shocked gas. X-ray emission was only detected in 3 of 60 PNs observed with ROSAT, followed by a few more from CHANDRA and XMM (e.g. Guerrero et al. 2005). In the case of PPNs, there is so far only one confirmed X-ray detection (Sahai et al. 2003), although many have been observed with CHANDRA. The general problem of understanding the formation and shaping of PNs has been addressed analytically and numerically (see Balick & Frank 2002, and references therein). X-ray emission is one of the most direct probes of the fast wind and the interaction process which drives PN formation. However, this probe has not been fully exploited in previous modeling studies. Very recently, Akashi et al. (2006, hereafter ASB06) used analytical, self-similar, spherically symmetric models to address this problem. We are performing numerical simulations with the hydrodynamics code FLASH varying the basic parameters of the fast and slow wind over an extensive parameter grid and computing the X-ray emission as a function of these parameters and the time history of the fast wind.Keywords. planetary nebulae: general, X-rays: stars, ISM: jets and outflows, stars: mass loss, hydrodynamics
Simulations & ResultsUsing Flash, we solve the differential equations of ideal hydrodynamics, and include radiative losses using the cooling functions of Sutherland & Dopita (1993) and Dalgarno & McCray (1972). We let a fast post-AGB wind (mass loss rateṀ f , velocity v f ) expand into a spherically symmetric, slowly expanding AGB wind (Ṁ s , v s ). The values of these four parameters cover the range adopted from ASB06 for ease of comparison; we have expanded the range to fit models to observations. The expected X-ray flux is determined using the radial density and temperature profiles from the hydrodynamical simulations. The atomic database ATOMDB with IDL (Smith et al. 2001) is used to calculate the emissivity. We compute the total flux in X-rays, the spectrum and emission maps.Our 1-D spherical models are in good agreement with analytical solutions of the structure and kinematics of the hot bubble (Chevalier & Imamura 1983 resulting from the wind-wind interaction. These models also support many of the analytical results of ASB06, although interesting differences occur. For example, we find that the initial phase of low X-ray luminosity, relevant for the understanding of X-rays from PPNe, is significantly shorter than in ASB06. We believe that our results are more accurate because of the self-consistent treatment of radiative cooling and hydrodynamics in our models.
519h...