skip to main content

Title: Common-envelope shaping of planetary nebulae – IV. From protoplanetary to planetary nebula

We present 2D hydrodynamical simulations of the transition of a protoplanetary nebula (PPN) to a planetary nebula for central stars in binary systems that have undergone a common-envelope event. After 1000 yr of magnetically driven dynamics (PPN phase), a line-driven stellar wind is introduced into the computational domain and the expansion of the nebula is simulated for another 10 000 yr, including the effects of stellar photoionization. In this study we consider central stars with main sequence (final) masses of 1 (0.569) and 2.5 (0.677) M⊙, together with a 0.6-M⊙ main-sequence companion. Extremely bipolar, narrow-waisted PPNe result in bipolar planetary nebulae, while the rest of the shapes mainly evolve into elliptical planetary nebulae. The initial magnetic field’s effects on the collimated structures, such as jets, tend to disappear in most of the cases, leaving behind the remnants of those features in only a few cases. Equatorial zones fragmented mainly by photoionization (1-M⊙progenitors), result in ‘necklace’ structures made of cometary clumps aligned with the radiation field. On the other hand, fragmentation by photoionization and shocked wind (2.5-M⊙progenitors) give rise to the formation of multiple clumps in the latitudinal direction, which remain within the lobes, close to the center, which are immersed and surrounded more » by hot shocked gas, not necessarily aligned with the radiation field. These results reveal that the fragmentation process has a dependence on the stellar-mass progenitor. This fragmentation is made possible by the distribution of gas in the previous post-common-envelope PPN as sculpted by the action of the jets.

« less
; ;
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 3822-3831
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT Current models predict that binary interactions are a major ingredient in the formation of bipolar planetary nebulae (PNe) and pre-planetary nebulae (PPNe). Despite years of radial velocity (RV) monitoring, the paucity of known binaries amongst the latter systems means data are insufficient to examine this relationship in detail. In this work, we report on the discovery of a long-period (P = 2654 ± 124 d) binary at the centre of the Galactic bipolar PPN IRAS 08005−2356 (V510 Pup), determined from long-term spectroscopic and near-infrared time-series data. The spectroscopic orbit is fitted with an eccentricity of 0.36 ± 0.05, which is similar to that of other long-period post-AGB binaries. Time-resolved Hα profiles reveal high-velocity outflows (jets) with deprojected velocities up to 231$_{-27}^{+31}$ km s−1 seen at phases when the luminous primary is behind the jet. The outflow traced by Hα is likely produced via accretion on to a main-sequence companion, for which we calculate a mass of 0.63 ± 0.13 M⊙. This discovery is one of the first cases of a confirmed binary PPN and demonstrates the importance of high-resolution spectroscopic monitoring surveys using large telescopes in revealing binarity among these systems.
  2. The 12C/13C ratio has been measured toward a sample of planetary nebulae (PNe) using millimeter observations of CO, HCN, HNC, CN, and other species, conducted with the 12 m antenna and the Submillimeter Telescope of the Arizona Radio Observatory. The observed nebulae spanned the entire lifetime of PNe, from ∼900 to 12,000 yr, and include well-known objects such as NGC 7293 (Helix), NGC 6720 (Ring), and NGC 2440, as well as relatively unexplored nebulae (M3–28, M2–48, and M3–55). In most cases, multiple molecules and transitions were used in the ratio determination, resulting in the most accurate values available to date, with 10%–40% uncertainties. The ratios found were unexpectedly low, lying in the range 12C/13C ∼1.0 ± 0.7–13.2 ± 4.9, with an average value of 3.7—drastically less than found in the envelopes of C-rich AGB stars, and, in some cases, lower than the minimum value achieved in equilibrium CNO burning. Such low values are expected for the two O-rich nebulae studied (M2–9 and M2–48), because of insufficient third dredge-up events. However, most of the PNe observed were clearly carbon-rich, as deduced from the large number of C-bearing molecules present in them. Because nucleosynthesis ceases in the PN stage, both the C/Omore »and the 12C/13C ratios must reflect abundances at the end of the AGB. These consistently low 12C/13C ratios, combined with the bipolar/multipolar morphologies of all planetary nebulae observed, suggest an explosive process involving proton-capture occurred at the AGB–PN transition.« less
  3. Using a code that employs a self-consistent method for computing the effects of photo-ionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive stars. We take into account changes in stellar properties and mass-loss over the star’s evolution. Our simulations show how various properties, such as the density and ionization fraction, change throughout the evolution of the star. The multi-dimensional simulations reveal the presence of strong ionization front instabilities in the main-sequence phase, similar to those seen in galactic ionization fronts. Hydrodynamic instabilities at the interfaces lead to the formation of filaments and clumps that are continually being stripped off and mixed with the low density interior. Even though the winds start out as completely radial, the spherical symmetry is quickly destroyed, and the shocked wind region is manifestly asymmetrical. The simulations demonstrate that it is important to include the effects of the photoionizing photons from the star, and simulations that do not include this may fail to reproduce the observed density profile and ionization structure of wind-blown bubbles around massive stars.
  4. ABSTRACT The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and the dense toroidal-shaped ejecta left over from common envelope (CE) evolution. Here we use the 3D hydrodynamic adaptive mesh refinement (AMR) code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the smoothed particle hydrodynamics (SPH) code phantom serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta pre-PN outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post-CE windsmore »may be aspherical, our models show that collimation via ‘inertial confinement’ will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock-focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC-driven bipolar lobes.« less
  5. ABSTRACT We discuss a theoretical model for the early evolution of massive star clusters and confront it with the ALMA, radio, and infrared observations of the young stellar cluster highly obscured by the molecular cloud D1 in the nearby dwarf spheroidal galaxy NGC 5253. We show that a large turbulent pressure in the central zones of D1 cluster may cause individual wind-blown bubbles to reach pressure confinement before encountering their neighbours. In this case, stellar winds energy is added to the hot shocked wind pockets of gas around individual massive stars that leads them to meet and produce a cluster wind in time-scales less than 105 yr. In order to inhibit the possibility of cloud dispersal, or the early negative star formation feedback, one should account for mass loading that may come, for example, from pre-main-sequence (PMS) low-mass stars through photoevaporation of their protostellar discs. Mass loading at a rate in excess of 8 × 10−9 M⊙ yr−1 per each PMS star is required to extend the hidden star cluster phase in this particular cluster. In this regime, the parental cloud remains relatively unperturbed, while pockets of molecular, photoionized and hot gas coexist within the star-forming region. Nevertheless, the most likely scenariomore »for cloud D1 and its embedded cluster is that the hot shocked winds around individual massive stars should merge at an age of a few million of years when the PMS star protostellar discs vanish and mass loading ceases that allows a cluster to form a global wind.« less