Issue |
A&A
Volume 557, September 2013
|
|
---|---|---|
Article Number | A121 | |
Number of page(s) | 12 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/201321786 | |
Published online | 18 September 2013 |
Unveiling shocks in planetary nebulae⋆
1 Instituto de Astrofísica de Andalucía, IAA-CSIC, Glorieta de la Astronomía s/n, 18008 Granada, Spain
e-mail: mar@iaa.es
2 Instituto de Astrofísica de Canarias, IAC, Vía Láctea s/n, 38205 La Laguna, Spain
3 Departamento de Astrofísica, Universidad de La Laguna, 38200 La Laguna, Tenerife, Spain
4 Departamento de Física Aplicada, Facultade de Ciencias, Campus Lagoas-Marcosende s/n, Universidade de Vigo, 36310 Vigo, Spain
5 Consejo Superior de Investigaciones Científicas (CSIC), c/ Serrano 117, 28006 Madrid, Spain
6 Departament de Física i Enginyeria Nuclear, EUETIB, Universitat Politècnica de Catalunya, Comte d’Urgell 187, 08036 Barcelona, Spain
7 Departament d’Astronomia i Meteorologia, Universitat de Barcelona, Av. Diagonal 647, 08028 Barcelona, Spain
8 Instituto de Ciencias Nucleares, UNAM, Ciudad Universitaria, 04510 Mexico City, Mexico
Received: 26 April 2013
Accepted: 25 July 2013
The propagation of a shock wave into a medium is expected to heat the material beyond the shock, producing noticeable effects in intensity line ratios such as [O iii]/Hα. To investigate the occurrence of shocks in planetary nebulae (PNe), we have used all narrow-band [O iii] and Hα images of PNe available in the HST archive to build their [O iii]/Hα ratio maps and to search for regions where this ratio is enhanced. Regions with enhanced [O iii]/Hα emission ratio can be ascribed to two different types of morphological structures: bow-shock structures produced by fast collimated outflows and thin skins enveloping expanding nebular shells. Both collimated outflows and expanding shells are therefore confirmed to generate shocks in PNe. We also find regions with depressed values of the [O iii]/Hα ratio which are found mostly around density bounded PNe, where the local contribution of [N ii] emission into the F656N Hα filter cannot be neglected.
Key words: ISM: jets and outflows / ISM: general / planetary nebulae: general
Tables 1 and 2 are available in electronic form at http://www.aanda.org
© ESO, 2013
1. Introduction
The increase in the local pressure of the material in a region makes it to expand, producing a disturbance in the surroundings that may drive a shock front. Shocks are ubiquitous in the universe whenever supernova explosions, photoionized gas, stellar winds, and collimated outflows drive shock waves in the circumstellar or interstellar media.
Shocks of fast moving outflows with material in the surrounding medium should produce appreciable effects. The most obvious is probably the formation of a bow-shock structure, as those associated with jet-like outflows in Herbig-Haro objects such as the archetypical HH 34 and HH 111 (Reipurth et al. 1997, 2002; Raga et al. 2002). Planetary nebulae (PNe) also exhibit fast collimated outflows whose occurrence is found at very early stages of the nebular formation, including the late stages of the asymptotic giant branch (AGB) and the proto-PN phases. These collimated outflows can produce notable bow-shock features, as those detected in IC 4593 (Corradi et al. 1997a).
On microscopic scales, a shock wave is an irreversible process, where entropy is generated as ordered kinetic energy is dissipated into heat. Therefore, a fast outflow propagating into a low density medium drives a forward shock which produces an increase of the electron temperature, Te. The increase in temperature would enhance the emission in the [O iii] λ5007 Å line, which is very sensitive to Te, whereas the low density would reduce the emissivity of the Hα line, which is more sensitive to the electron density, Ne. Consequently, it is expected a region in front of the bow-shock at the tip of a collimated outflow where the [O iii]/Hα ratio is significantly enhanced. Such caps of enhanced [O iii]/Hα ratio have been revealed by HST images of IC 4634, a PN characterized by the presence of fast collimated outflows resulting in a point-symmetric morphology and clear bow-shock structures (Guerrero et al. 2008).
A dense, fast shell expanding into a tenuous, static medium can also propagate shocks into the surroundings. This effect is observed in bubbles blown by the wind of massive stars, such as S 308, the wind-blown bubble around the Wolf-Rayet star HD 50896, whose shell generates a forward shock that produces a notable offset between the [O iii] and Hα emissions (Gruendl et al. 2000). Another example is provided by NGC 7635, a.k.a. the Bubble Nebula, a bubble blown by the wind of the massive O6.5 IIIf star BD+60°2522, where the offset between these emissions is ≈3 × 1015 cm (Moore et al. 2002). The expansion of the different shells of PNe is also expected to generate shocks: one at the leading edge of the outer shell propagating into the unperturbed AGB wind and another at the leading edge of the bright rim expanding into the outer shell (e.g., Perinotto et al. 2004). Balick (2004) found a thin skin of enhanced [O iii]/Hα enveloping NGC 6543 in HST images of this nebula. The origin of this skin of enhanced [O iii]/Hα remained uncertain, but an observational artefact could be ruled out.
Fast collimated outflows are acknowledged to play an essential role in the nebular shaping in early stages of the PN formation (Sahai & Trauger 1998). In particular, the expansion of bow-shock features has drawn the attention as a mechanism for the formation of bipolar and multi-polar PNe (see the review by Balick & Frank 2002). Similarly, the expansion of the different nebular shells is critical for the evolution of PNe (Villaver et al. 2002; Perinotto et al. 2004). Both the shocks produced by fast collimated outflows and those produced by the shell expansion may have important dynamical effects in PNe, contributing decisively to the nebular shaping, evolution, and excitation.
Inspired by the results obtained for IC 4634 and NGC 6543, we have used the HST [O iii] and Hα images of PNe available in the Mikulski Archive for Space Telescopes (MAST) to investigate the occurrence of skins of bright [O iii]/Hα in PNe and its relationship with fast collimated outflows and expanding shells. The description of the search and analysis of HST images is presented in Sect. 2. Attending to their [O iii]/Hα ratio maps (Sect. 3), the PNe in this sample can be divided in four different types of which two of them show regions of enhanced [O iii]/Hα emission associated with bow-shocks of collimated outflows and with expanding nebular shells, respectively. The results are further discussed in Sect. 4 and a short summary is given in Sect. 5.
Fig. 1 Hα and [O iii]/Hα ratio maps of the type A PNe IC 4593, IC 4634, NGC 3918, NGC 6210, NGC 6543, NGC 6572, and NGC 7009. In the [O iii]/Hα ratio maps, regions of relatively bright Hα emission appear bright (white), whereas regions of strong [O iii] emission relative to Hα are shown dark (black). |
2. Hα and [O iii] images of PNe in the HST archive
We searched MAST for HST WFPC2 or WFC3 coeval Hα and [O iii] images of PNe available by March 2013. This search yielded Hα and [O iii] images for 103 PNe obtained through the F656N and F502N filters, respectively (Table 1). Since the adjacent [N ii] λλ6548,6584 Å emission lines can be included in the bandpass of the F656N Hα narrow-band filter, the potential contamination of the emission from these low-excitation lines on the Hα image needs to be discussed into further detail.
The bandpass of the WFC3 F656N Hα filter only includes the contribution from the [N ii] λ6548 Å emission line with a filter transmission at its wavelength which is ~5% that at the wavelength of the targeted Hα line (O’Dell et al. 2013). As for the WFPC2 F656N Hα filter, its bandpass is broader (Lim et al. 2010) and includes contributions of both [N ii] lines, with filter transmissions relative to that at the wavelength of the Hα line of ~33% for the [N ii] λ6548 Å component and ~5% for the [N ii] λ6584 Å line. Therefore, we can expect a 10% contamination by [N ii] emission in the WFC3 F656N filter when the [N ii] λ6584 Å line is 6 times brighter than the Hα line, but only ~2/3 as bright for the WFPC2 F656N filter.
Fig. 2 Same as Fig. 1 for the type B PNe NGC 1501, NGC 2371-2, NGC 3242, NGC 6153, NGC 6818, NGC 6826, NGC 6884, and NGC 7662. |
The [N ii] to Hα line ratio may be large for low-excitation PNe with bright [N ii] emission lines. To preserve our sample from sources whose Hα image may be significantly contaminated by the contribution of [N ii] emission, hereafter we will exclude from our analysis those objects observed by the WFPC2 instrument with a [N ii] λ6584 Å to Hα line ratio greater than unity. This includes H 1-9, H 2-15, Hb 5, IC 4406, K 4-55, M 2-14, Mz 3, NGC 2438, NGC 2440, NGC 2818, NGC 3132, and NGC 6537. The source NGC 6302, observed with the WFC3 instrument, will be further excluded because its extremely high [N ii] λ6584 Å to Hα line ratio. The [O iii]/Hα ratio maps of sources with a [N ii] λ6584 Å to Hα line ratio close to unity will be carefully scrutinized to check for features associated to [N ii] contamination. We note that the (Hα+[N ii])/[O iii] ratio maps of PNe notably enhance the [N ii] emission of low-ionization small-scale structures (Corradi et al. 1996). These features are easy to identify.
The images were downloaded from the HST archive and reduced using the standard pipeline procedure. For each PN, individual frames of the same epoch were combined to remove the cosmic rays using IRAF1 routines. The resulting Hα and [O iii] images were subsequently used to generate [O iii] to Hα ratio maps. Before producing these ratios, we checked the relative alignment of the images in the different filters. The WFC3 images were notably well aligned, but we noticed a consistent shift of 01–03 between the WFPC2 Hα and [O iii] images caused most likely by the different light path across the filters. Accordingly, the images in the different filters were aligned using stars in the field of view (including the central stars when available) and registered on the same grid. The accuracy of this procedure was usually better than 10% of a pixel, i.e., ≲0005.
The [O iii]/Hα ratio maps were examined to investigate the relative positions of the [O iii] and Hα emissions and to search for regions of enhanced values of the [O iii]/Hα ratio. The images of 18 PNe were found not suitable for this analysis: NGC 6853 (the Dumbbell Nebula) and NGC 7293 (the Helix Nebula) are much larger than the WFPC2 field of view, SaSt 1-1 and Th 4-1 are point sources, Hen 2-436, NGC 6833, and PC 11 are compact PNe whose emission is not properly resolved, and H 1-43, H 1-55, H 2-25, H 2-26, Hen 2-262, JaFu 2, K 3-76, KFL 7, Th 3-6, Th 3-14, and Wray 16-282 are weak in [O iii], resulting in poor quality [O iii]/Hα ratio maps. All these 18 PNe have been discarded from subsequent analysis.
3. PN classification based on the [O iii]/Hα ratio
The remaining 72 PNe in our sample present [O iii]/Hα ratio maps suitable for its analysis. The Hα and [O iii]/Hα ratio images of a representative number of PNe with adequate [O iii]/Hα ratio maps are presented in Figs. 1 through 4. Grey-scales in the [O iii]/Hα ratio have been chosen to highlight the locations where the [O iii] and Hα emissions show spatial differences. A close inspection of the Hα and [O iii]/Hα ratio images of these PNe has allowed us to classify them into four different types:
-
Type A PNe that present caps of enhanced[O iii]/Hα ratios at the tip of jet-like features. Theseregions of bright [O iii] emission can be interpretedas bow-shock features produced by fast collimated outflows. TheHα images and [O iii]/Hα ratio maps of the type A PNe IC 4593, IC 4634, NGC 3918, NGC 6210, NGC 6543, NGC 6572, and NGC 7009 are shown in Fig. 1.
-
Type B PNe surrounded by skins of enhanced [O iii]/Hα emission associated with nebular shells. In some double-shell PNe, both shells exhibit a skin of enhanced [O iii]/Hα. The list of type B PNe includes NGC 1501, NGC 2371-2, NGC 3242, NGC 6153, NGC 6818, NGC 6826, NGC 7009, and NGC 7662 whose Hα images and [O iii]/Hα ratio maps are shown in Fig. 2. Although some of these PNe have Fast Low-Ionization Emission Regions (the so-called FLIERs), the skins of enhanced [O iii]/Hα ratios described here are not related to these features.
-
Type C PNe where the [O iii]/Hα ratio decreases at the outer edge of the nebula, i.e., these PNe show the opposite behavior to type B PNe. Among these PNe are: IC 418, NGC 6720, BD+30°3639, MyCn 18, Hen 3-1357, NGC 6881, NGC 6886, M 1-20, M 2-14, M 2-43, and Hen 2-260. The Hα image and [O iii]/Hα ratio map of these PNe are shown in Fig. 3.
-
Type D PNe whose [O iii]/Hα ratio maps do not show the effects described above (Fig. 4). Some of them have flat [O iii]/Hα ratio maps (e.g., Wray 16-286), some have low signal-to-noise (S/N) ratio maps (e.g., M 3-14), and some of them show too complex structures (e.g., NGC 2392 and NGC 6369) to be clearly ascribed to type A, B, or C.
The four different types described above are not mutually exclusive. In particular, most type A PNe are also type B, i.e., they show regions of enhanced [O iii]/Hα ratios that are associated both with expanding nebular shells and bow-shock features produced by fast collimated outflows. For these PNe, we use the type AB. The only exception is that of NGC 6572, where a collimated outflow towards the north of the nebula produces a bow-shock feature with enhanced values of the [O iii]/Hα ratio, whereas this ratio is diminished around the main nebular shell. For this PN, we use the type AC.
The PNe with adequate [O iii]/Hα ratio maps are listed in Table 2. In this table we provide the classification of each PN and relevant information for a comparative study: the morphology or nebular shape (E = elliptical, B = bipolar, MP = quadrupolar or multi-polar, PS = point-symmetric, MS = multiple shell), the temperature of the central star, the nebular electron density, an averaged angular radius, the distance and linear radius, the [O iii] λ5007 and [N ii] λ6584 to Hβ line intensity ratios, and comments on the presence of collimated outflows, FLIERs, and ansae, or on the quality of the observations. Most of the data in this table have been obtained from Frew (2008). Additional references are given in the last column of Table 2.
4. Discussion
4.1. Statistical properties
A preliminary inspection of Table 2 provides hints of the varying nature of PNe belonging to the different types described in the previous section. These are further illustrated by the distributions of different parameters shown in Fig. 5.
Type D PNe is an heterogeneous group severely affected by several observational biases. Most type D PNe are small, with angular radii ≤4′′, and linear sizes ≤0.03 pc. The small size of these objects is clearly introducing an observational bias in the sample, as the fine morphological details of these objects cannot be properly resolved. Furthermore, the S/N of the images of a significant number of type D PNe is limited, thus hampering the detection of structures with contrasting values of the [O iii]/Hα ratio. On the other hand, some objects in this group are large and well resolved, but they have very complex morphologies (e.g., NGC 2392 and NGC 6369). Although these sources have low [N ii] λ6584 Å to Hα integrated line ratios, the presence of low-excitation features may result in local regions of significant [N ii] emission. The [N ii] contamination on the Hα would correspond to features of apparently diminished values of the [O iii]/Hα ratio. The complexity of the ratio maps may hide the possible effects of the expansion of nebular shells and collimated outflows in the surrounding material.
PNe of type C also tend to be small, with radii and linear sizes similar to, but slightly larger than, those of type D PNe. Notably, some of the type C PNe (e.g., NGC 2346 and NGC 6720) have the largest angular radii and linear sizes among the PNe in this sample. In contrast to type D PNe, the distribution of the effective temperature of the central stars of type C PNe clearly peaks at low values (Teff ≤ 50 000 K). PNe of type C are either young, with low Teff central stars (e.g., BD+30°3639 and IC 418), or have hot central stars which have already turned towards the white dwarf cooling track and reduced their stellar luminosities (e.g., the evolved PN NGC 6720 and the prominent bipolar PNe NGC 2346 and NGC 6886, Napiwotzki 1999; Méndez & Niemela 1981; Pottasch & Surendiranath 2005). In the first case, the stellar spectrum is not hard enough to maintain the oxygen twice ionized throughout the whole nebula, whereas in the second case the spectrum is very hard but the photon flux seen by material at large distances from the central star is small. In both cases, the ionization degree drops at the edge of the nebula and an exterior O+ region forms. Indeed, the reduced number of type C PNe with available HST WFPC2 images in the [O ii] λ3727 Å line consistently reveal bright emission in regions of low [O iii]/Hα. The corresponding high value of the [O ii]/[O iii] line ratio in these regions can be interpreted as a sudden drop in the ionization degree. The decrease in the ionization degree will also enhance the emission of the low-excitation [N ii] lines. Therefore, the reduced [O iii]/Hα values observed in the ratio maps of IC 418 (Ramos-Larios et al. 2012) and NGC 6720 (O’Dell et al. 2013) among others can be partially ascribed to the local contaminant contribution in these regions of [N ii] emission on the Hα image.
Fig. 5 Distributions of the angular radii, linear sizes, [N ii]/[O iii] ratios, and stellar effective temperatures of the different types of PNe defined in Sect. 3. |
Fig. 6 Radial profiles of the [O iii]/Hα ratio maps of the type A PNe IC 4593, NGC 6210, and NGC 7009 along selected bow-shock features associated with collimated outflows. |
PNe of type A and B are, on average, better resolved than type C PNe, with angular radii spanning between ~4′′ and ~20′′, and linear sizes greater than 0.02 pc, peaking at ~0.1 pc. The distribution of the effective temperature of their central stars (median value ~75 000 K) is somehow similar to that of type C PNe (median value ~60 000 K). More obvious differences are revealed by the distributions of their [N ii]/[O iii] ratios: those of type A and B peak at low values of the [N ii]/[O iii] ratio, with most data points below 0.1, whereas the distribution of type C PNe shows a notable tail towards high values of the [N ii]/[O iii] ratio. If we could exclude abundances effects, this result would imply that type C PNe have lower excitation than type A and B PNe. Finally, we remark that the occurrence of multiple-shell morphology and collimated outflows or FLIERs (suggestive of fast collimated outflows, although with expansion velocities comparable to these of the nebular shells) is the highest among PNe of type A and B.
4.2. [O iii]/Hα radial profiles
To investigate in more detail the physical properties of the regions where the variations of the [O iii]/Hα ratio take place, we have extracted radial profiles of representative PNe of types A, B, and C. These radial profiles, shown in Figs. 6–8, respectively, have been derived from flux calibrated HST images; however, we have not attempted to deredden the line ratios and, therefore, the values of the [O iii]/Hα ratio in these profiles do not correspond to absolute line intensity ratios. The relative enhancement (or decrement for type C PNe) of the [O iii]/Hα ratio and the width of the region with anomalous [O iii] emission derived from these profiles are listed in Table 3.
Fig. 7 Radial profiles of the [O iii]/Hα ratio maps of the type B PNe NGC 3242 a), b), c) and NGC 6826 d). The radial profiles of NGC 3242 correspond to the inner shell a), the outer shell b), and the north-west blister of the outer shell c). The radial profile of NGC 6826 d) corresponds to the outer shell. The [O iii]/Hα peaks of the outer shell of NGC 3242 are not shown in the profile of the blister of this shell c). |
The regions of enhanced [O iii]/Hα ratio associated with bow-shocks, outer shells, and blisters show the largest increase of this ratio, being typically larger by a factor of two. On the contrary, the enhancement of the [O iii]/Hα ratio at the inner shells of type B PNe is typically smaller. Those inner shells are embedded within the corresponding outer shells, and thus we interpret that the smaller enhancement of the [O iii]/Hα ratio is the result of the dilution of the emission from the skin of high [O iii]/Hα ratio by the emission from the outer shell.
Fig. 8 Radial profiles of the [O iii]/Hα ratio maps of the type C PNe IC 418, NGC 6720, and NGC 6886 along their minor axes at PA ~ 70°, 130°, and 45°, respectively. |
There are also notable differences among the fractional sizes of these regions. The peaks of enhanced [O iii] emission in bow-shocks, outer shells, and blisters show the narrowest fractional sizes, all being ≤4%. Meanwhile, the regions of enhanced [O iii]/Hα ratios associated with the inner shell of NGC 3242 and the absorptions in the [O iii]/Hα radial profiles of type C PNe are broader, with fractional widths above 10%. The physical sizes of all these regions are rather similar, ≈1016 cm, but for the outer shell of NGC 6826, which is especially thin (4 × 1015 cm). Even the outer regions of diminished [O iii]/Hα ratios of type C have similar physical sizes, except for NGC 6720, which is especially broad (≈1017 cm).
There are two interesting features in these radial profiles which are worth mentioning as they are indicative of the variations in excitation affecting the nebular material and of the geometry of these regions. First, the radial profiles of bow-shocks in Fig. 6 display a depression just inside the emission peak which is then followed by a bounce in this ratio further in. The regions of diminished values in the [O iii]/Hα ratio are associated with the prominent low-ionization features found inside these bow-shocks, whereas the increase in the value of the line ratio suggests a relative brightening of the [O iii] emission on larger spatial scales than the outer skin. Secondly, the [O iii]/Hα peak associated with outer shells is preceded by a ramp in this ratio, as can be seen in the radial profiles of the outer shells of NGC 3242 and NGC 6826 (Fig. 7). Projection effects can be responsible for this feature.
4.3. Origin of the high [O iii]/Hα skins
As mentioned in Sect. 4.1, the [O iii]/Hα ratio maps of type C PNe can be explained in terms of a nebular ionization structure in which the O++/H+ ratio decreases in the outermost nebular layers. The low effective temperature of the central stars of young PNe or the low luminosity of the central stars of bipolar or evolved PNe may result in these low ionization degrees.
Similarly, one might seek an explanation to the [O iii]/Hα ratio maps of type A and B PNe in terms of a nebular ionization structure: the outward increase in the [O iii] emission would be caused by a decrease of the O3+/H+ ratio and a subsequent increase in the O++/H+ ratio. If this were the case, one would expect low [O iii] emission throughout the nebula and its brightening at the outer edge, mimicking the structure of the [O ii] emission of type C PNe. However, it turns out that the brightness of the [O iii] λ5007 Å line is consistently high throughout the whole nebula in both type A and B PNe. Furthermore, the O3+/O++ fraction required to reproduce such a structure would require stellar temperature and luminosities much larger than those typical of PNe (the maximum O3+/O++ values observed in PNe are about 0.5). This indicates that the local enhancement of the [O iii]/Hα ratio in type A and B PNe is most probably caused by a sharp increase in temperature, which cannot be reproduced by photoionization models, unless an ad hoc additional heating mechanism working on small spatial scales produces a sharp raise in Te. Such an increase in Te would reduce the effective recombination rate of Hα, which is roughly proportional to , and increase the collisional excitation of O++, which is roughly proportional to T− 1/2 e− E/kT. The [O iii]/Hα line intensity ratio can be expressed as: (1)We plot in Fig. 9 the variation of the [O iii]/Hα line intensity ratio with the increase of electron temperature above a baseline temperature in the range typical of PNe, 8000 K ≤ Te ≤ 12 000 K. We note that these curves have not been computed using the above expression, but the exact line emissivities. For a baseline temperature of the nebular material of 10 000 K, the observed increase of the [O iii]/Hα line intensity ratio by factors of 2 to 3 implies a jump in Te of 2000–5000 K.
The extra heating can be provided by a forward shock propagating into a tenuous medium (Cox 1972) as it has been observed in collimated outflows in Herbig-Haro objects (e.g., HH 34, Rodríguez-González et al. 2012) and ring nebulae around WR stars (Gruendl et al. 2000). The extra heating in type A PNe is provided by the forward shock of a fast outflow propagating into the tenuous, outer regions of the PNe or the interstellar medium (ISM). The shocks produced by collimated outflows in Herbig-Haro objects enhance the emission of low-ionization lines (e.g., [O i], [S ii], or [N ii]), as predicted for steady, plane-parallel shock models (Hartigan et al. 1987), but the situation is notably different in PNe. Their outflows are photoionized by the central star, and this results in variations of the expected line ratios and spatial distribution of the emission from different lines that has been modeled for fast, compact knots traveling away from the central star by Raga et al. (2008). In particular, these models predict a relative brightening of the [O iii] emission with respect to that of Hα at larger distances from the ionizing source (their Fig. 9), in clear agreement with the radial profiles displayed in Fig. 6.
Fig. 9 Variation of the [O iii]/Hα line intensity ratio associated with the increase in electron temperature from to Te. Three different curves are plotted, depending on the value of the pre-shock temperature. |
Meanwhile, the expansion of a higher density shell into a lower density, outer shell or into the ISM produces a forward shock in type B PNe which heats the material in front of the shell. In this case, the pre-shock material is already highly ionized by the central star of the PN, and no enhancement of the emission from low-ionization lines is expected. Therefore, the analysis of [O iii]/Hα ratio maps provides the only strategy to investigate the occurrence of the shocks associated to expanding shells of PNe.
Hydrodynamical models of PN evolution indeed predict the development of two different shocks during their formation (e.g., Perinotto et al. 2004; Schönberner et al. 2005). The sudden ionization of the material ejected during the AGB produces a D-type ionization front at the leading edge of the outer nebular shell of multiple-shell PNe. This shock propagates outwards, into the ambient AGB wind, and a thin skin of shocked material is expected at the leading shock of the shell. Model simulations of the propagation of such shock kindly performed by Schönberner and Steffen suggest that the post-shock temperatures depend very sensitively on the density jump across the shock, although a high numerical resolution is required to resolve the shock. Similarly, another shock is set by the wind interaction at the discontinuity between the outer nebular shell and the bright rim pushed by the thermal pressure of the hot bubble. The available model simulations predict the formation of a thick rim early in the evolution of the PN, resulting in a shock too weak to produce the observed enhancement of the [O iii]/Hα ratios. With time, the rim sharpens and moves faster, and its leading shock becomes stronger (Perinotto et al. 2004).
These detailed hydrodynamical simulations determine that the velocity jump across the shock front at the leading edge of the outer shell is 5–40 km s-1. Across the shock front between the inner rim and the shell, the velocity jump is smaller, ≤10 km s-1. For typical electron densities of ~4000 cm-3, ~1000 cm-3, and ~50 cm-3 for the rim, outer shell, and surrounding unperturbed AGB wind, these velocity jumps imply temperature increases ≃4000 K which can easily produce the observed enhancement of the [O iii]/Hα ratios.
5. Summary
Using archival HST [O iii] and Hα images of PNe, we have found regions of enhanced values of the [O iii]/Hα ratio. These regions can be described as thin skins which envelop either bow-shock features associated with fast collimated outflows PNe or nebular shells of PNe with multiple-shell morphology. The enhancement of the [O iii]/Hα ratio is interpreted to be caused by a local increase of Te which brightens the emission in the [O iii] line and damps the emissivity of the Hα line. The amplitude of the variations in the values of the [O iii]/Hα ratio implies variations in Te of 2000–5000 K. The increase in Te is associated with forward shocks produced by the expansion of collimated outflows and nebular shells into surrounding media of lower density. These results illustrate the complexity of PNe and its suitability to studying the detailed physics of ionized plasmas.
Online material
Sample of PNe with HST WFPC2 or WFC3 Hα and [O iii] observations.
[O iii]/Hα type and basic parameters of planetary nebulae.
Acknowledgments
All of the data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant NNX09AF08G and by other grants and contracts. This work is funded by grants AYA 2005-01495 of the Spanish MEC (Ministerio de Educación y Ciencia), grant FGM-1747 of the Junta de Andalucía, and CSIC grant 2006-50-I135. M.A.G. and J.A.T. are supported by the Spanish MICINN (Ministerio de Ciencia e Innovación) grant AYA 2011-29754-C03-02, V.L. by grant AYA 2011-22614, L.F.M. by grant AYA 2011-30228-C03-01, and A.R. by grant AYA 2011-30228-C03, all of them co-funded with FEDER funds. J.A.T. acknowledges support by the CSIC JAE-Pre student grant 2011-00189, and PVF by CONACyT grant 167611 and DGAPA-PAPIIT (UNAM) grant IN105312. We kindly thanks Detlef Schönberner and Matthias Steffen for helpful discussion.
References
- Acker, A., Raytchev, B., Koeppen, J., & Stenholm, B. 1991, A&AS, 89, 237 [NASA ADS] [Google Scholar]
- Acker, A., Marcout, J., Ochsenbein, F., et al. 1992, The Strasbourg-ESO Catalogue of Galactic Planetary Nebulae. Parts I, II., eds. A. Acker, J. Marcout, F. Ochsenbein, et al., European Southern Observatory, Garching (Germany), 1047 [Google Scholar]
- Bässgen, M., Diesch, C., & Grewing, M. 1990, A&A, 237, 201 [NASA ADS] [Google Scholar]
- Balick, B. 2004, AJ, 127, 2262 [NASA ADS] [CrossRef] [Google Scholar]
- Balick, B., & Frank, A. 2002, ARA&A, 40, 439 [NASA ADS] [CrossRef] [Google Scholar]
- Bohigas, J. 2001, Rev. Mex. Astron. Astrophys., 37, 237 [Google Scholar]
- Cavichia, O., Costa, R. D. D., & Maciel, W. J. 2010, Rev. Mex. Astron. Astrophys., 46, 159 [NASA ADS] [Google Scholar]
- Corradi, R. L. M., Manso, R., Mampaso, A., & Schwarz, H. E. 1996, A&A, 313, 913 [NASA ADS] [Google Scholar]
- Corradi, R. L. M., Guerrero, M., Manchado, A., & Mampaso, A. 1997a, New A, 2, 461 [NASA ADS] [CrossRef] [Google Scholar]
- Corradi, R. L. M., Perinotto, M., Schwarz, H. E., & Claeskens, J.-F. 1997b, A&A, 322, 975 [NASA ADS] [Google Scholar]
- Corradi, R. L. M., Sabin, L., Miszalski, B., et al. 2011, MNRAS, 410, 1349 [NASA ADS] [CrossRef] [Google Scholar]
- Cox, D. P. 1972, ApJ, 178, 143 [NASA ADS] [CrossRef] [Google Scholar]
- Dennis, T. J., Cunningham, A. J., Frank, A., et al. 2008, ApJ, 679, 1327 [NASA ADS] [CrossRef] [Google Scholar]
- Ercolano, B., Morisset, C., Barlow, M. J., Storey, P. J., & Liu, X.-W. 2003, MNRAS, 340, 1153 [NASA ADS] [CrossRef] [Google Scholar]
- Ercolano, B., Wesson, R., Zhang, Y., et al. 2004, MNRAS, 354, 558 [NASA ADS] [CrossRef] [Google Scholar]
- Frew, D. 2008, Ph.D. Thesis, Macquarie University, Sydney, Australia [Google Scholar]
- Gruendl, R. A., Chu, Y.-H., Dunne, B. C., & Points, S. D. 2000, AJ, 120, 2670 [NASA ADS] [CrossRef] [Google Scholar]
- Guerrero, M. A., Miranda, Luis F., Riera, Angels, et al. 2008, ApJ, 683, 272 [NASA ADS] [CrossRef] [Google Scholar]
- Guiles, S., Bernard-Salas, J., Pottasch, S. R., & Roellig, T. L. 2007, ApJ, 660, 1282 [NASA ADS] [CrossRef] [Google Scholar]
- Hartigan, P., Raymond, J., & Hartmann, L. 1987, ApJ, 316, 323 [NASA ADS] [CrossRef] [Google Scholar]
- Hyung, S., & Aller, L. H. 1996, MNRAS, 278, 551 [NASA ADS] [CrossRef] [Google Scholar]
- Hyung, S., & Aller, L. H. 1997, ApJ, 491, 242 [NASA ADS] [CrossRef] [Google Scholar]
- Hyung, S., Aller, L. H., & Feibelman, W. A. 1994a, MNRAS, 269, 975 [NASA ADS] [Google Scholar]
- Hyung, S., Aller, L. H., & Feibelman, W. A. 1994b, PASP, 106, 745 [NASA ADS] [CrossRef] [Google Scholar]
- Hyung, S., Aller, L. H., & Feibelman, W. A. 1997, ApJS, 108, 503 [NASA ADS] [CrossRef] [Google Scholar]
- Hyung, S., Aller, L. H., & Feibelman, W. A. 1999, ApJ, 525, 294 [NASA ADS] [CrossRef] [Google Scholar]
- Hyung, S., Aller, L. H., Feibelman, W. A., & Lee, S.-J. 2001, ApJ, 563, 889 [NASA ADS] [CrossRef] [Google Scholar]
- Kwitter, K. B., & Henry, R. B. C. 1998, ApJ, 493, 247 [NASA ADS] [CrossRef] [Google Scholar]
- Lee, C. F., & Sahai, R. 2003, ApJ, 586, 319. [NASA ADS] [CrossRef] [Google Scholar]
- Lim, P. L., Quijada, M., Baggett, S., et al. 2010, Space Telescope Science Institute, Instrument Science Report WFPC2 2010-05 [Google Scholar]
- Méndez, R. H., & Niemela, V. S. 1981, ApJ, 250, 240 [NASA ADS] [CrossRef] [Google Scholar]
- Méndez, R. H., Kudritzki, R. P., Herrero, A., Husfeld, D., & Groth, H. G. 1988, A&A, 190, 113 [NASA ADS] [Google Scholar]
- Moore, B. D., Walter, D. K., Hester, J. J., et al. 2002, AJ, 124, 3313 [NASA ADS] [CrossRef] [Google Scholar]
- Napiwotzki, R. 1999, A&A, 350, 101 [Google Scholar]
- O’Dell, C. R., Ferland, G. J., Henney, W. J., & Peimbert, M. 2013, AJ, 145, 92 [NASA ADS] [CrossRef] [Google Scholar]
- Perinotto, M., Schönberner, D., Steffen, M., & Calonaci, C. 2004, A&A, 414, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Phillips, J. P. 2004, MNRAS, 353, 589 [NASA ADS] [CrossRef] [Google Scholar]
- Phillips, J. P., & Cuesta, L. 1998, A&AS, 133, 381 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pottasch, S. R., & Surendiranath, R. 2005, A&A, 432, 139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pottasch, S. R., Gathier, R., Gilra, D. P., & Wesselius, P. R. 1981, A&A, 102, 237 [NASA ADS] [Google Scholar]
- Pottasch, S. R., Bernard-Salas, J., Beintema, D. A., & Feibelman, W. A. 2004, A&A, 423, 593 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pottasch, S. R., Beintema, D. A., & Feibelman, W. A. 2005, A&A, 436, 953 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pottasch, S. R., Bernard-Salas, J., & Roellig, T. L. 2009, A&A, 499, 249 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Raga, A. C., Velázquez, P. F., Cantó, J., & Masciadri, E. 2002, A&A, 395, 647 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Raga, A.C., Riera, A., Mellema, G., Esquivel, A., & Velázquez, P. F. 2008, A&A, 489, 1141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ramos-Larios, G., Vázquez, R., Guerrero, M. A., et al. 2012, MNRAS, 423, 3753 [NASA ADS] [CrossRef] [Google Scholar]
- Reipurth, B., Hartigan, P., Heathcote, S., Morse, J. A., & Bally, J. 1997, AJ, 114, 757 [NASA ADS] [CrossRef] [Google Scholar]
- Reipurth, B., Heathcote, S., Morse, J., Hartigan, P., & Bally, J. 2002, AJ, 123, 362 [NASA ADS] [CrossRef] [Google Scholar]
- Rodríguez-González, A., Esquivel, A., Raga, A. C., et al. 2012, AJ, 143, 60 [NASA ADS] [CrossRef] [Google Scholar]
- Ruiz, M. T., Peimbert, A., Peimbert, M., & Esteban, C. 2003, ApJ, 595, 247 [NASA ADS] [CrossRef] [Google Scholar]
- Ruiz, N., Guerrero, M. A., Chu, Y.-H., & Gruendl, R. A. 2011, AJ, 142, 91 [NASA ADS] [CrossRef] [Google Scholar]
- Sabbadin, F. 1976, A&A, 52, 291 [NASA ADS] [Google Scholar]
- Sahai, R., & Trauger, J. T. 1998, AJ, 116, 1357. [NASA ADS] [CrossRef] [Google Scholar]
- Schönberner, D., Jacob, R., Steffen, M., et al. 2005, A&A, 431, 963 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Turatto, M., Cappellaro, E., Ragazzoni, R., Benetti, S., & Sabbadin, F. 2002, A&A, 384, 1062 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vázquez, R. 2012, ApJ, 751, 116 [NASA ADS] [CrossRef] [Google Scholar]
- Villaver, E., Manchado, A., & García-Segura, G. 2002, ApJ, 581, 1204 [NASA ADS] [CrossRef] [Google Scholar]
- Yuan, H.-B., Liu, X.-W., Péquignot, D., et al. 2011, MNRAS, 411, 1035 [NASA ADS] [CrossRef] [Google Scholar]
- Walton, N. A., Reay, N. K., Pottasch, S. R., & Atherton, P. D. 1986, New Insights in Astrophysics. Eight Years of UV Astronomy with IUE, 263, 497 [NASA ADS] [Google Scholar]
- Wesson, R., & Liu, X.-W. 2004, MNRAS, 351, 1026 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
All Figures
Fig. 1 Hα and [O iii]/Hα ratio maps of the type A PNe IC 4593, IC 4634, NGC 3918, NGC 6210, NGC 6543, NGC 6572, and NGC 7009. In the [O iii]/Hα ratio maps, regions of relatively bright Hα emission appear bright (white), whereas regions of strong [O iii] emission relative to Hα are shown dark (black). |
|
In the text |
Fig. 2 Same as Fig. 1 for the type B PNe NGC 1501, NGC 2371-2, NGC 3242, NGC 6153, NGC 6818, NGC 6826, NGC 6884, and NGC 7662. |
|
In the text |
Fig. 3 Same as Fig. 1 for the type C PNe. |
|
In the text |
Fig. 4 Same as Fig. 1 for a subsample of type D PNe. |
|
In the text |
Fig. 5 Distributions of the angular radii, linear sizes, [N ii]/[O iii] ratios, and stellar effective temperatures of the different types of PNe defined in Sect. 3. |
|
In the text |
Fig. 6 Radial profiles of the [O iii]/Hα ratio maps of the type A PNe IC 4593, NGC 6210, and NGC 7009 along selected bow-shock features associated with collimated outflows. |
|
In the text |
Fig. 7 Radial profiles of the [O iii]/Hα ratio maps of the type B PNe NGC 3242 a), b), c) and NGC 6826 d). The radial profiles of NGC 3242 correspond to the inner shell a), the outer shell b), and the north-west blister of the outer shell c). The radial profile of NGC 6826 d) corresponds to the outer shell. The [O iii]/Hα peaks of the outer shell of NGC 3242 are not shown in the profile of the blister of this shell c). |
|
In the text |
Fig. 8 Radial profiles of the [O iii]/Hα ratio maps of the type C PNe IC 418, NGC 6720, and NGC 6886 along their minor axes at PA ~ 70°, 130°, and 45°, respectively. |
|
In the text |
Fig. 9 Variation of the [O iii]/Hα line intensity ratio associated with the increase in electron temperature from to Te. Three different curves are plotted, depending on the value of the pre-shock temperature. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.