Cometary Knots Research Articles

Planetary Nebula NGC 2818: Revealing its complex 3D morphology

ABSTRACT We carry out an advanced morpho-kinematic analysis of the Planetary Nebula (PN) NGC 2818, whose complex morphology is described by a basic bipolar component, filamentary structures and a knotty central region. We performed an upgrated 3D Morpho-kinematic (MK) model by employing the shape software, combining for the first time in PNe optical 2D spatially resolved echelle spectra and Fabry–Perot data cubes. The best-fitting 3D model of NGC 2818 successfully reconstructs the main morphology, considering one bipolar component, radial filamentary structures, and an equatorial component as the geometrical locus of the group of cometary knots. The model shows that the equatorial component has the lower expansion velocity of the system at 70 ± 20 km s−1. The velocity of the bipolar component is 120 ± 20 km s−1, while all the filamentary structures were found to expand at higher velocities of 180 ± 20 km s−1. Moreover, Fabry–Perot data revealed for the first time a north-eastern filament expanding at a mean velocity of 80 ± 20 km s−1, while its equivalent counterpart in the south-western region was confirmed by a new detected substructure in the echelle data. A new detected knotty structure at velocity −40 ± 20 km s−1 is also reported, as expelled material from the fragmented eastern lobe of the nebula. We interpret the overall structure of NGC 2818 as the result of the evolution of a binary system that underwent the common envelope phase, in conjunction with the ejections of a magnetized jet, misaligned with respect to the symmetry axis of the bipolar/elliptical shell.

Radiation driven-dust hydrodynamics in late-phase AGB stars

The interplay of stellar luminosity variations and dust hydrodynamics in Asymptotic Giant Branch (AGB) stars and the consequences for dust survival and mass loss remain elusive. In this work, we broadly investigate the role of dust and radiation hydrodynamics in forming dust and gas structures, heterogeneous clumps observable in AGB remnants and planetary nebulae (PNe). Of interest in this study are the spatial perturbations driven by instabilities in the space that the mass travels through. These spatial perturbations in the dust and gas field may be responsible for forming larger clumps, such as cometary knots, seen in the PNe phase. Previous studies have considered similar physics in dust-driven winds at shorter lengths and time scales, using either 1D simulations or 2D simulations with a single mixed particle–gas fluid. Here we present an Eulerian–Lagrangian method for studying this problem at larger length and time scales. Simulations are performed in 2D, solving the Euler equations with source terms resulting from the particle phase, represented by free Lagrangian points. Radiation coupling was implemented for the particle phase, modeling radiation heating and acceleration of the particles and subsequent coupling to the gas phase through non-continuum heat and momentum transfer models. Spatial perturbations of the dust and radiation fields were found to drive the formation of small dust mass clumps that survive to late times, though these remain below the size of those observed in many PNe.

Molecular Hydrogen Microstructures in Planetary Nebulae

Molecular hydrogen (H 2 ) emission is commonly detected in planetary nebulae (PNe), specially in objects with bipolar morphologies. New studies showed that H 2 gas is also packed in microstructures embedded in PNe of any morphological type. Despite the presence of H 2 in cometary knots being known for years, only in the last five years, much deeper imagery of PNe have revealed that H 2 also exists in other types of low-ionisation microstructures (LISs). Significant differences are found between the host PNe of cometary knots and other types of LISs, such as nebula age, central star temperature (evolutionary stage) and the absolute sizes of the microstructure itself.

OH+ emission from cometary knots in planetary nebulae

We model the molecular emission from cometary knots in planetary nebulae (PNe) using a combination of photoionization and photodissociation region (PDR) codes, for a range of central star properties and gas densities. Without the inclusion of ionizing extreme ultraviolet (EUV) radiation, our models require central star temperatures $T_*$ to be near the upper limit of the range investigated in order to match observed H$_2$ and OH$^+$ surface brightnesses consistent with observations - with the addition of EUV flux, our models reproduce observed OH$^+$ surface brightnesses for $T_* \ge 100 \, {\rm kK}$. For $T_* < 80 \, {\rm kK}$, the predicted OH$^+$ surface brightness is much lower, consistent with the non-detection of this molecule in PNe with such central star temperatures. Our predicted level of H$_2$ emission is somewhat weaker than commonly observed in PNe, which may be resolved by the inclusion of shock heating or fluorescence due to UV photons. Some of our models also predict ArH$^+$ and HeH$^+$ rotational line emission above detection thresholds, despite neither molecule having been detected in PNe, although the inclusion of photodissociation by EUV photons, which is neglected by our models, would be expected to reduce their detectability.

Extreme Radio-wave Scattering Associated with Hot Stars

Abstract We use data on extreme radio scintillation to demonstrate that this phenomenon is associated with hot stars in the solar neighborhood. The ionized gas responsible for the scattering is found at distances up to from the host star, and on average must comprise ∼105 distinct structures per star. We detect azimuthal velocities of the plasma, relative to the host star, up to , consistent with warm gas expanding at the sound speed. The circumstellar plasma structures that we infer are similar in several respects to the cometary knots seen in the Helix and in other planetary nebulae. There the ionized gas appears as a skin around tiny molecular clumps. Our analysis suggests that molecular clumps are ubiquitous circumstellar features, unrelated to the evolutionary state of the star. The total mass in such clumps is comparable to the stellar mass.

BUBBLES AND KNOTS IN THE KINEMATICAL STRUCTURE OF THE BIPOLAR PLANETARY NEBULA NGC 2818

High-resolution Hubble Space Telescope (HST) archive imaging and high-dispersion spectroscopy are used to study the complex morphological and kinematical structure of the planetary nebula, NGC 2818. We analyze narrow band H{\alpha}, [O III], [N II], [S II] and He II images, addressing important morphological features. Ground-based longslit echelle spectra were obtained crossing NGC 2818 at five different positions to precisely determine kinematical features in the structure of the nebula. A distance of 2.5 kpc was used to determine physical scales. Constructing models to fit the data with modern computational tools, we find NGC 2818 is composed by: (1) a non-uniform bipolar structure with a semi-major axis of 0.92 pc (75 arcsec), possibly deformed by the stellar wind, (2) a 0.17 pc (14 arcsec) diameter central region, which is potentially the remnant of an equatorial enhancement, and (3) a great number of cometary knots. These knots are preferentially located inside a radius of 0.24 pc (20 arcsec) around the central star. The ma jor axis of the main structure is oriented at i \simeq 60\degree with respect to the line-of-sight and at PA = +89\degree on the plane of the sky. Expansion velocities of this nebula are Vpol = 105 km/s and Veq = 20 km/s, which lead to our estimate of the kinematical age of {\tau}k \simeq 8,400 \pm 3, 400 yr (assuming homologous expansion). Our observations do not support the idea that high velocity collimated ejections are responsible for the formation of microstructures inside the nebula. We determine the systemic velocity of NGC 2818 to be VHEL = +26 \pm 2 km/s.

NGC 6778: a disrupted planetary nebula around a binary central star

The planetary nebula (PN) NGC 6778 harbors a binary central star with a short orbital period and displays two systems of fast collimated outflows. In order to assess the influence of the evolution through a common-envelope phase of the binary system of NGC 6778 on its formation and shaping, we have used narrow-band images and high-dispersion long-slit spectra of the nebula to investigate its detailed morphology and kinematics. We find that the overall structure of NGC 6778 can be described as a bipolar PN. The equatorial ring is highly disrupted and many radial features (filamentary wisps and cometary knots) also evidence strong dynamical effects. There are clear connections between the bipolar lobes and the fast collimated outflows: the collimated outflows seem to arise from bright knots at the tips of the bipolar lobes, whereas the kinematics of the bipolar lobes is distorted. We suggest that the interaction of the fast collimated outflows of NGC 6778 with its nebular envelope has resulted in the disruption of the nebular shell and equatorial ring.

THE CIRCUMSTELLAR ENVIRONMENT OF R CORONAE BOREALIS: WHITE DWARF MERGER OR FINAL-HELIUM-SHELL FLASH?

In 2007, R Coronae Borealis (R CrB) went into an historically deep and long decline. In this state, the dust acts like a natural coronagraph at visible wavelengths, allowing faint nebulosity around the star to be seen. Imaging has been obtained from 0.5 to 500 micron with Gemini/GMOS, HST/WFPC2, Spitzer/MIPS, and Herschel/SPIRE. Several of the structures around R CrB are cometary globules caused by wind from the star streaming past dense blobs. The estimated dust mass of the knots is consistent with their being responsible for the R CrB declines if they form along the line of sight to the star. In addition, there is a large diffuse shell extending up to 4 pc away from the star containing cool 25 K dust that is detected all the way out to 500 micron. The SED of R CrB can be well fit by a 150 AU disk surrounded by a very large diffuse envelope which corresponds to the size of the observed nebulosity. The total masses of the disk and envelope are 10^-4 and 2 M(Sun), respectively, assuming a gas-to-dust ratio of 100. The evidence pointing toward a white-dwarf merger or a final-helium-shell flash origin for R CrB is contradictory. The shell and the cometary knots are consistent with a fossil planetary nebula. Along with the fact that R CrB shows significant Lithium in its atmosphere, this supports the final-helium-shell flash. However, the relatively high inferred mass of R CrB and its high fluorine abundance support a white-dwarf merger.

Modelling the warm H2 infrared emission of the Helix nebula cometary knots

Molecular hydrogen emission is commonly observed in planetary nebulae. Images taken in infrared H2 emission lines show that at least part of the molecular emission is produced inside the ionised region. In the best-studied case, the Helix nebula, the H2 emission is produced inside cometary knots (CKs), comet-shaped structures believed to be clumps of dense neutral gas embedded within the ionised gas. Most of the H2 emission of the CKs seems to be produced in a thin layer between the ionised diffuse gas and the neutral material of the knot, in a mini photodissociation region (PDR). However, PDR models published so far cannot fully explain all the characteristics of the H2 emission of the CKs. In this work, we use the photoionisation code \textsc{Aangaba} to study the H2 emission of the CKs, particularly that produced in the interface H^+/H^0 of the knot, where a significant fraction of the H2 1-0S(1) emission seems to be produced. Our results show that the production of molecular hydrogen in such a region may explain several characteristics of the observed emission, particularly the high excitation temperature of the H2 infrared lines. We find that the temperature derived from H2 observations even of a single knot, will depend very strongly on the observed transitions, with much higher temperatures derived from excited levels. We also proposed that the separation between the H_alpha and NII peak emission observed in the images of CKs may be an effect of the distance of the knot from the star, since for knots farther from the central star the NII line is produced closer to the border of the CK than H_alpha.

A “FIREWORK” OF H2KNOTS IN THE PLANETARY NEBULA NGC 7293 (THE HELIX NEBULA)

We present a deep and wide field-of-view (4'x 7') image of the planetary nebula (PN) NGC 7293 (the Helix Nebula) in the 2.12 micron H2 v=1-0 S(1) line. The excellent seeing (0.4'') at the Subaru Telescope, allows the details of cometary knots to be examined. The knots are found at distances of 2.2'-6.4' from the central star (CS). At the inner edge and in the inner ring (up to 4.5' fromthe CS), the knot often show a `tadpole' shape, an elliptical head with a bright crescent inside and a long tail opposite to the CS. In detail, there are variations in the tadpole shapes, such as narrowing tails, widening tails, meandering tails, or multi-peaks within a tail. In the outer ring (4.5'-6.4' from the CS), the shapes are more fractured, and the tails do not collimate into a single direction. The transition in knot morphology from the inner edge to the outer ring is clearly seen. The number density of knots governs the H2 surface brightness in the inner ring: H2 exists only within the knots. Possible mechanisms which contribute to the shaping of the knots are discussed, including photo-ionization and streaming motions. A plausible interpretation of our images is that inner knots are being overrun by a faster wind, but that this has not (yet) reached the outer knots. Based on H2 formation and destruction rates, H2 gas can survive in knots from formation during the late asymptotic giant branch (AGB) phase throughout the PN phase. These observations provide new constraints on the formation and evolution of knots, and on the physics of molecular gas embedded within ionized gas.

THE OUTFLOWS AND THREE-DIMENSIONAL STRUCTURE OF NGC 6337: A PLANETARY NEBULA WITH A CLOSE BINARY NUCLEUS

NGC 6337 is a member of the rare group of planetary nebulae where a close binary nucleus has been identified. The nebula's morphology and emission line profiles are both unusual, particularly the latter. We present a thorough mapping of spatially resolved, long-slit echelle spectra obtained over the nebula that allows a detailed characterization of its complex kinematics. This information, together with narrow band imagery is used to produce a three dimensional model of the nebula using the code SHAPE. The 3-D model yields a slowly expanding toroid with large density fluctuations in its periphery that are observed as cometary knots. A system of bipolar expanding caps of low ionization are located outside the toroid. In addition, an extended high velocity and tenuous bipolar collimated outflow is found emerging from the core and sharply bending in opposite directions, a behavior that cannot be accounted for by pure magnetic launching and collimation unless the source of the outflow is precessing or rotating, as could be expected from a close binary nucleus.

Merged Ionization/Dissociation Fronts in Planetary Nebulae

The hydrogen ionization and dissociation front around an ultraviolet radiation source should merge when the ratio of ionizing photon flux to gas density is sufficiently low and the spectrum is sufficiently hard. This regime is particularly relevant to the molecular knots that are commonly found in evolved planetary nebulae, such as the Helix Nebula, where traditional models of photodissociation regions have proved unable to explain the high observed luminosity in H_2 lines. In this paper we present results for the structure and steady-state dynamics of such advection-dominated merged fronts, calculated using the Cloudy plasma/molecular physics code. We find that the principal destruction processes for H_2 are photoionization by extreme ultraviolet radiation and charge exchange reactions with protons, both of which form H_2^+, which rapidly combines with free electrons to undergo dissociative recombination. Advection moves the dissociation front to lower column densities than in the static case, which vastly increases the heating in the partially molecular gas due to photoionization of He^0, H_2, and H^0. This causes a significant fraction of the incident bolometric flux to be re-radiated as thermally excited infrared H_2 lines, with the lower excitation pure rotational lines arising in 1000 K gas and higher excitation H_2 lines arising in 2000 K gas, as is required to explain the H_2 spectrum of the Helix cometary knots.

VLT/near-infrared integral field spectrometer observations of molecular hydrogen lines in the knots of the planetary nebula NGC 7293 (the Helix Nebula)

Knots are commonly found in nearby planetary nebulae (PNe) and star forming regions. Within PNe, knots are often found to be associated with the brightest parts of the nebulae and understanding the physics involved in knots may reveal the processes dominating in PNe. As one of the closest PNe, the Helix Nebula (NGC 7293) is an ideal target to study such small-scale (~300 AU) structures. We have obtained infrared integral spectroscopy of a comet-shaped knot in the Helix Nebula using SINFONI on the Very Large Telescope at high spatial resolution (50-125 mas). With spatially resolved 2 micron spectra, we find that the H2 rotational temperature within the cometary knots is uniform. The rotational-vibrational temperature of the cometary knot (situated in the innermost region of the nebula, 2.5 arcmin away from the central star), is 1800 K, higher than the temperature seen in the outer regions (5-6 arcmin from the central star) of the nebula (900 K), showing that the excitation temperature varies across the nebula. The obtained intensities are reasonably well fitted with 27 km s-1 C-type shock model. This ambient gas velocity is slightly higher than the observed [HeII] wind velocity of 13 km s-1. The gas excitation can also be reproduced with a PDR (photo dominant region) model, but this requires an order of magnitude higher UV radiation. Both models have limitations, highlighting the need for models that treats both hydrodynamical physics and the PDR.

Infrared Observations of the Helix Planetary Nebula

We have mapped the Helix (NGC 7293) planetary nebula (PN) with the IRAC instrument on the Spitzer Space Telescope. The Helix is one of the closest bright PNs and therefore provides an opportunity to resolve the small-scale structure in the nebula. The emission from this PN in the 5.8 and 8 μm IRAC bands is dominated by the pure rotational lines of molecular hydrogen, with a smaller contribution from forbidden line emission such as [Ar III] in the ionized region. The IRAC images resolve the cometary knots, which have been previously studied in this PN. The tails of the knots and the radial rays extending into the outer regions of the PN are seen in emission in the IRAC bands. IRS spectra on the main ring and the emission in the IRAC bands are consistent with shock-excited H_2 models, with a small (~10%) component from photodissociation regions. In the northeast arc, the H_2 emission is located in a shell outside the Hα emission.

The Origin and Physical Properties of the Cometary Knots in NGC 7293

On the basis that the cometary knots observed in the Helix Nebula form as a result of larger parent clouds breaking up due to Rayleigh-Taylor instability induced by radiative acceleration of the clouds, we compute characteristics of the cometary knots and of the parent clouds as well. Present observations place constraints on the positions, velocities, and sizes of the parent clouds. Requiring the clouds to produce cometary knots that are stable places further constraints on the properties of parent clouds. We formulate those constraints and show how they further limit the predicted properties of the cometary knots and bring them into agreement with their observed properties.

3D radiation hydrodynamics: Interacting photo-evaporating clumps

We present the results of a new radiation hydrodynamics code called Maartje. This code describes the evolution of a flow in three spatial dimensions using an adaptive mesh, and contains a combination of a ray tracer and an atomic physics module to describe the effects of ionizing radiation. The code is parallelized using a custom threadpool library. We present an application in which we follow the ionization of two dense spherical clumps which are exposed to an ionizing radiation field from a 50 000 K black body. We study various configurations in which one of the clumps shields the other from the ionizing photons. We find that relatively long-lived filamentary structures with narrow tails are formed. This raises the possibility that cometary knots (such as are found in the Helix Nebula) may be the result of the interaction of an ionizing radiation field with an ensemble of clumps, as opposed to the identification of a single knot with a single clump.

Molecular Hydrogen in the Ring Nebula: Clumpy Photodissociation Regions

ABSTRACTWe present a 0.″65 resolution H2 1–0 S(1) 2.122 μm image of the Ring Nebula (NGC 6720), which was taken with the Near Infrared Imager at the WIYN 3.5 m telescope on Kitt Peak. The high resolution of the H2 observation is sufficient to reveal the finer structure of the molecular material in this nebula. The morphology of the molecular emission is compared to that of the ionized emission from the Ring Nebula as seen by the Hubble Space Telescope (HST; Heii, [Oiii], and [Nii]), and it is clear that the dark clumps seen by HST match the locations of clumpy H2 emission, suggesting that these clumps are similar to the cometary knots seen in the Helix Nebula. As with the Helix, the clumpy H2 emission from the main ring of the Ring Nebula is contained within the optically bright ionized nebula, implying that the molecular gas is shielded inside dense condensations. Comparison of the observed H2 average surface brightnesses for the Ring Nebula [(1.5 ± 0.5) × 10-4 ergs cm−2 s−1 sr−1] with time‐dependent models of the expected H2 emission from planetary nebulae (PNe) shows that it is consistent with H2 excitation in photodissociation regions (PDRs), confirming previous suggestions. Comparison of the Ring Nebula H2 emission with a younger PN, NGC 2346, and an older PN, the Helix Nebula, suggests an evolution in H2 surface brightness consistent with the time‐dependent PDR models. Moreover, the knots of molecular gas appear to become more isolated as the PN evolves, consistent with optical studies of knots in PNe.

High resolution molecular hydrogen imaging of the Ring Nebula

The relative morphologies and structures of molecular and ionized gas emission from planetary nebulae (PNe) allow a better understanding of the nature and evolution of these objects. The classical paradigm for the structure of PNe is that of an ionized gas bubble bounded by neutral gas and molecules. However, it has been shown that molecular gas exists within ionized regions, leading to a re-evaluation of the classic structure. In the Helix Nebula (NGC 7293) dense condensations known as cometary knots are known to exist in the main ionized nebula. The molecules in these knots are shielded from the ionizing radiation and thus survive within the ionized zone. Another PN in which H2 emission is seen to originate from within the ionized nebula in NGC 6720 (the Ring Nebula).

Large-Scale Extended Emission around the Helix Nebula: Dust, Molecules, Atoms, and Ions

We present new observations of the ionized gas, molecular gas, and cool dust in the Helix Nebula (NGC 7293). The ionized gas is observed in the form of an Himage, which is constructed using images from the Southern HSky Survey Atlas. The molecular emission was mapped using the H2 v =1 ! 0 S(1) line at 2.122 lm. The far-infrared (FIR) observations were obtained using ISOPHOT on the Infrared Space Observ- atory. The Hobservations are more sensitive than previous measurements and show the huge extent of the Helix, confirming it as a density-bounded nebula and showing previously unseen point-symmetric structures. The H2 observations show that the molecular gas follows the distribution of molecular material shown in pre- vious work. The molecular emission is confined to that part of the nebula seen in the classic optical image. Furthermore, comparison of the H2 emission strength with time-dependent models for photodissociation regions (PDRs) shows that the emission arises from thermal excitation of the hydrogen molecules in PDRs and not from shocks. The FIR observations, at 90 and 160 lm, represent mostly contributions from thermal dust emission from cool dust grains but include a small contribution from ionized atomic lines. Comparison of the FIR emission with the Hobservation shows that the dust and ionized gas are coincident and extend to � 1100 00 radius. This equates to a spatial radial extent of more than 1 pc (assuming a distance to the Helix of � 200 pc). Assuming that the outer layers of the circumstellar shell have spherical symmetry, radiative transfer modeling of the emission in Hgives a shell mass of � 1.5 M� . However, the modeling does not cover the outermost part of the shell (beyond � 600 00 radius), and therefore this is a lower limit for the shell mass. Moreover, the models suggest the need for very large dust grains, with � 80% of the dust mass in grains larger than 3.5 lm. Comparison of these new observations with previous observations shows the large-scale stratifi- cation of the Helix in terms of ionized gas and dust, as well as the coexistence of molecular species inside the ionized zones, where molecules survive in dense condensations and cometary knots.

Photoevaporating Flows from the Cometary Knots in the Helix Nebula (NGC 7293)

We explain the Ha emission of the cometary knots in the Helix Nebula (NGC 7293) with an analytical model that describes the emission of the head of the globules as a photoevaporated flow produced by the incident ionizing radiation of the central star.We compare these models with the Ha emission obtained from the HST archival images of the Helix Nebula. From a comparison of the Ha emission with the predictions of the analytical model we obtain a rate of ionizing photons from the central star of about 5e45 s^-1, which is consistent with estimates based on the total Hb flux of the nebula. We also model the tails of the cometary knots as a photoevaporated wind from a neutral shadow region produced by the diffuse ionizing photon field of the nebula. A comparison with the HST images allows us to obtain a direct determination of the value of the diffuse ionizing flux. We compare the ratio of diffuse to direct stellar flux as a function of radius inside an HII region with those obtained from the observational data through the analytical tail and head wind model. The agreement of this model with the values determined from the observations of the knots is excellent.