Probing the Mass and Structure of the Ring Nebula in Lyra with SOFIA/GREAT Observations of the [CII] 158 micron Line
Key Words.:planetary nebulae: individual (NGC 6720) – stars: winds, outflows – circumstellar matter – galactic: ISM – atoms: [CII]
We have obtained new velocity-resolved spectra of the [CII] 158 m line towards the Ring Nebula in Lyra (NGC 6720), one of the best-studied planetary nebulae, in order to probe its controversial 3-dimensional structure and to estimate the mass of circumstellar material in this object. We used the Terahertz receiver GREAT aboard the SOFIA airborne telescope to obtain the [CII] spectra at eight locations within and outside the bright optical ring of NGC 6720. Emission was detected at all positions except for the most distant position along the nebula’s minor axis, and generally covers a broad velocity range, V km s (FWZI), except at a position along the major axis located just outside the optical ring, where it is significantly narrower (V km s). The one narrow spectrum appears to be probing circumstellar material lying outside the main nebular shell that has not been accelerated by past fast wind episodes from the central star, and therefore most likely comes from equatorial and/or low-latitude regions of this bipolar nebula. Along lines-of-sight passing within about 10″ of the nebular center, the CII column density is a factor 46 higher than the CO column density. The total mass of gas associated with the [CII] emission inside a circular region of diameter is at least 0.11. A significant amount of [CII] flux arises from a photodissociation region immediately outside the bright optical ring, where we find a CII to CO ratio of , lower than that seen towards the central region. Comparing our data with lower-quality CI spectra, which indicate similarly large CI/CO ratios in NGC 6720, we conclude that the bulk of elemental carbon in NGC 6720 is divided roughly equally between CII and CI, and that the emissions from these species are far more robust tracers of circumstellar material than CO in this object and other evolved planetary nebulae.
Planetary Nebulae (PNe) result when “ordinary” stars, with main-sequence masses 1-8, die extraordinary deaths. During their post-main-sequence evolutionary phase (i.e., as AGB stars), these stars eject half or more of their total mass in the form of nucleosynthetically-enriched material into the interstellar medium (ISM) – a process which dramatically alters the course of stellar evolution, and plays a key role in the chemical evolution of the Galaxy. PNe are formed when most of the stellar envelope has been ejected and the hot central star’s radiation field ionizes the ejecta.
PNe are almost a unique astrophysical class in terms of being bright at all wavelength bands from X-ray to radio, and are ideal laboratories for investigating important astrophysical phenomena such as nucleosynthesis, (magneto)hydrodynamics of colliding winds, astrochemistry, photodissociation, and the evolution of dust grains (via UV and shock-processing) because their geometries and kinematics are generally believed to be well understood (e.g., Gurzadyan 1997). The history of the mass-loss processes which result in a PN, beginning from its life as an AGB star, is written in the spatio-kinematic structure of the circumstellar ejecta surrounding the central star. Different geometries can require radically different formation models – for example, in the binary-driven model for aspherical PNe, ellipsoidal shapes result from sub-stellar mass companions, whereas bipolar shapes require interaction with a stellar-mass companion (Soker 1996).
However, inferring the 3D physical structure of a PN from its observed 2D projection on the sky can be a difficult task. This is because PNe appear in a dazzling variety of complex morphologies, with collimated, bipolar and multipolar shapes being dominant, and round ones being rather rare, as clearly seen in imaging surveys with HST (e.g., Sahai & Trauger 1998; Sahai, Morris & Villar 2011). Furthermore, although PNe shapes have traditionally been defined by optical imaging, the latter typically do not show where the bulk of the mass in these objects lies: in a photodissociation region (PDR) surrounding the bright optical nebula, of up to 1 , compared to 0.1 in ionized gas (c.f., Hollenbach & Tielens 1999). This PDR may represent most of the mass ejected during the AGB phase; some fraction of it is sculpted and compressed by fast post-AGB winds from the central star, resulting in the complex shapes of PNe seen in optical imagery. The bulk of the PDR material radiates most of its energy via thermal emission in the far-infrared.
At a distance of 0.7 kpc, NGC 6720 (M 57), the Ring Nebula in Lyra, is among the brightest and most extensively studied PNe (e.g., O’Dell et al. 2007: ODell07, and references therein, van Hoof et al. 2010). NGC6720 has an elliptical ring shape (size ) in optical emission lines. This bright ring lies at the center of a large “double” halo structure: a brighter inner one (size ) consisting of a system of limb-brightened loops, and an outer faint one which is irregular, knotty and round (″ diameter: Balick et al. 1992). Although it is considered a “classical” nebula because of the bright ring’s simple 2D shape, its detailed 3D structure remains controversial, in spite of many studies undertaken to elucidate it (see, e.g., ODell07).
We report here the results of a small program conducted as part of Basic Science observations on the Stratospheric Observatory For Infrared Astronomy (SOFIA), to use velocity-resolved observations of [CII] 158 m line emission for probing the mass and 3-D structure of NGC 6720. The use of the [CII] 158 micron line as a probe stemmed from our expectation that, because [CII] is present both in the ionized and atomic components of the circumstellar environment, and because this line has a low critical density (cm), it should be readily observable both within and outside the optically-bright PN ring.
The observations were performed with the (dual-channel) German Receiver for
Astronomy at Terahertz Frequencies (GREAT
We detected emission at all positions except #9, where the rms noise is 0.09 K (Table 1). The emission generally covers a broad velocity range, V50 km s(full-width at zero intensity: FWZI), except at position #4, where it is significantly narrower (V km s) (Fig. 2). An average spectrum from the central (#1) and four symmetric locations on the major and minor axes (#3,#5,#8,#8a) shows a symmetric profile with a central double-peaked core of width V km s, centered at km s, and a weaker high-velocity pedestal with a total width of about 100 km s at its base. The observed width of the high-velocity component is limited by the signal-to-noise, and is likely larger than 100 km s (e.g., by folding and averaging this profile around its center velocity, we reduce the noise and find clear evidence for emission extending over 120 km s). The line profile towards position #1, although noisy, clearly does not show even a hint of a minimum around the systemic velocity. Since an expanding hollow shell would produce a strong minimum, we conclude that C is not absent from the region interior to the bright ring.
We compute the [CII] column density at different locations in the nebula from the observed [CII] intensity, (erg s cm sr), using equation A.5 of Schneider et al. (2003), cm (Table 1). We have reasonably assumed that the density and temperature are high enough for thermalized emission – the upper level of the [CII] line is 91 K above ground, and the critical densities for collisional excitation by electrons (at K), H, and H are 46, and cm, respectively (Stacey et al. 2010, Schneider et al. 2003), to be compared with the electron density in the ionized nebula (500–700 cm, ODell07), and the density ( cm) and temperature (300 K) in the PDR of NGC 6720 (Liu et al. 2001). Hence, with erg s cm sr towards the nebular center, we find cm.
|Label||(″), (″)||K . km s||cm|
coordinates are J2000 RA=18:53:35.08, DEC=33:01:45.0
Since we have sampled the nebula sparsely with beams separated by , we estimate the total [CII] flux by approximating the [CII] emitting region as a central, circular region of diameter, surrounded by two concentric annuli. The mean radius of annulus 1 (2) is equal to the offsets of position #2 (#3; same for #’s 5,8 & 8a) from the center, and each has a width of . We derive the average intensity in each region from the average of the spectra sampling it, and multiply it by the region’s area to derive the flux. We thus find that the total [CII] flux, erg s cm. Assuming all of the carbon is in C, that the distance is 0.7 kpc, and taking for NGC 6720 (Liu et al. 2001, Table 5), we find that the total mass of gas associated with the [CII] emission inside a circular region of diameter is at least 0.114. Our mass estimate is a very conservative lower limit because (i) we have not taken into account the emission from locations beyond the main ring, such as position #4, and (ii) a comparable fraction of carbon is likely to be neutral (see below), We note that Liu et al. (2001) estimated a mass of 0.03 from their [CII] 158 m line flux of erg s cm for NGC 6720, measured using the ISO Long Wavelength Spectrometer (LWS) with a single 70″ aperture covering the main optical nebula. Clearly, a significant amount of [CII] flux lies outside the 70″ ISO/LWS aperture, i.e., in the PDR beyond the main optical nebula.
We compare our [CII] results with those from studies of CO and CI emission in NGC 6720. For CO, Bachiller at al. (1989: BBMG89) have obtained complete maps in the J=1–0 and 2–1 lines, whereas for CI, Bachiller at al. (1994: BHCF94) detect emission in the 609 m line with low S/N at a few (3 of 5) positions along the minor axis. The CO and CI column densities towards the nebular center are cm (BBMG89) and cm (BHCF94). Hence, both the CII/CO and the CI/CO abundance ratios are about 50 towards the nebular center, i.e., almost all of the nebular carbon is distributed roughly equally between C and C. Towards positions offset from the nebular center, the CII/CO abundance ratio, while still large, is not as extreme. We find (i) cm towards position #4, and (ii) cm near position #4 from BBMG89’s CO J=2–1 map, using their spectrum at offset position (). Since this spectrum (obtained with a 13″ beam) samples the strongest emission peak in the CO map, whereas our #4 spectrum ( beam), is offset by from that peak, our derived C to CO ratio here is at least 6.5. Our CII study strongly supports the hypothesis that the CO emission region in NGC 6720 is really a PDR (BHCF94). We note that the CO J=2–1 emission shows significant clumpy structure, and given the strong [CII] emission seen at position #4 where the CO emission is also strong, and the lack of [CII] emission at #9 (3 upper limit of cm assuming same velocity width as at position #4), where the CO emission is very weak, it is likely that the [CII] is similarly clumpy.
Noting that the CII column densities at all positions we have sampled in NGC 6720 (Table 1) exceed the maximum CO column density in NGC 6720 by large factors, and that the average CI/CO ratio is 10 (BHCF94), we conclude that CII and CI are far more robust tracers of circumstellar mass in this object than CO. If we assume that the abundances of CII and CI are roughly equal, then our circumstellar mass estimate from [CII] data increases to . Hence, given estimates of the current central star mass of 0.61–0.62 (Odell07), and that the mass of the outer halo has not been accounted for, it is likely that the mass of NGC 6720’s progenitor star exceeded 1 .
Models for the 3-D structure of NGC 6720 can be divided into two broad classes. The first (model 1) is a (roughly) prolate ellipsoid, having open ends along the long axis: thus a barrel shape, with an equatorial density enhancement and the axis of symmetry inclined at 30 to the line of sight (Masson 1990, Guerrero et al. 1997: GMC97). This model has been refined by Hiriart (2004) and ODell07 (using mapping of the 2.1 m H line, and multislit optical spectra, respectively) to one in which the barrel cross-section is elliptical and it is seen nearly pole-on. In the second class of models (model 2), NGC 6720 is seen as a nearly pole-on bipolar nebula based on multislit spectra along the major axis (Bryce et al. 1994). Kwok et al. (2008) compare deep emission-line images of NGC 6720 with another well-studied PN, NGC 6853, and suggest that both have the same biconical 3-D structure, except that NGC 6853 is oriented edge-on, whereas NGC 6720 is nearly pole-on. In the Kwok et al. (2008) model, the inner and outer halo together represent the projections of three pairs of bi-conical outflows with similar opening angles, viewed along their common symmetry axes. From their CO J=2-1 study, BBMG89 infer a structure which has elements of both class 1 and class 2: a hollow cylindrical shape with its axis inclined to the line-of-sight.
We propose that NGC 6720 is a PN with a barrel-shaped central region viewed (nearly) along its axis, which can reconcile the above model classes. Barrel-shaped central regions have recently been recognized as an important feature of many PNe in a comprehensive new morphological classification system based on HST images of 119 young PNe by Sahai, Morris & Villar (2011). NGC 6720 is also multipolar, with the “flower-petal” structures in the H and [NII] images seen in the “inner halo” of NGC 6720 being the limb-brightened peripheries of these lobes. The presence of such lobes is also indicated by the appearance of multiple, faint, red- and blue-shifted velocity components within the overall position-velocity ellipse seen in optical long-slit spectra reported by GMC97 and interpreted as resulting from “bubbles of material (that) protrude from the main nebula”. The NGC 6720 ring shows several extended filamentary structures which suggests that a significant fraction of the ring’s optical line emission comes from the bright basal regions of multiple lobes emanating from the central region; the continuation of these regions to slightly higher latitudes is seen in the fainter, more diffuse, emission in the H image immediately outside the bright ring (Fig. 1). This suggestion finds support in GMC97’s finding that “the velocity structure of the edge of the nebula shows a complex behaviour”, “reflecting the projected velocity of these clumps and bubbles”. Optical spectra of the flower-petal region (Fig. 8 of GMC97) show distinct narrow components separated by km s. In our model, these result from the inclined walls of lobes on the near- and far-side of the equatorial plane. These spectra also show the presence of point-symmetric pairs of knotty structures which are consistent with an origin in oppositely-directed multiple outflows, rather than bubbles created by a fast wind flowing through a fragmented shell as GMC97 propose, because it is very unlikely that a stochastic process (fragmentation) can create point-symmetrically located holes in a shell. The multipolar PN He 2-47 (Sahai 2000) provides us a rough idea of what NGC 6720’s lobes might look like if the symmetry axis of its central barrel was in (or nearly in) the sky-plane.
Although a detailed interpretation of all the published data on NGC 6720 with our proposed model is outside the scope of this paper, our [CII] data appear consistent with this model, provided the lobes are expanding at similar velocities. In this model, the [CII] emission at positions 3, 5, 8, and 8A, which likely includes emission from both the barrel-shaped central region and low-latitude regions of multipolar lobes, is expected to be roughly similar, as observed. The spectrum at position #4 provides insight into the circumstellar material that lies outside the main nebular ring and whose physical properties are not well known – the narrow profile (V km s) suggests that the [CII] emission is probing material that has not been accelerated by past fast wind episodes from the central star, and therefore most likely comes from equatorial and low-latitude regions of the slowly-expanding progenitor AGB wind. In contrast, the emission from positions #8 and #8a which, like that from #4, also samples material largely just outside the shell’s periphery but along the minor axis, is broad (V kms) and may be associated with the multipolar lobes whose peripheries show up in projection as the “flower-petals” in the H image (Fig. 1).
In summary, our results show that for NGC 6720 (and by inference other evolved PNe), C is a much more robust tracer of the circumstellar material than CO because it is far more abundant. We therefore expect future GREAT [CII] observations to play a major role in the study of evolved PNe.
Acknowledgements.This study is based on observations made with the NASA/DLR Stratospheric Observatory for Infrared Astronomy (SOFIA). SOFIA Science Mission Operations is operated by the Universities Space Research Association, Inc., under NASA contract NAS2-97001, and the Deutsches SOFIA Institut under DLR contract 50 OK 0901. RS’s contribution to the research described here was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA. Financial support was provided by NASA through a Long Term Space Astrophysics award to RS and MM, and a SOFIA award to RS.
- GREAT is a development by the MPI für Radioastronomie and the Universität zu Köln, in cooperation with the MPI für Sonnensystemforschung and the DLR Institut für Planetenforschung.
- Balick, B., Gonzalez, G., Frank, A., & Jacoby, G. 1992, ApJ, 392, 582
- Bachiller, R., Bujarrabal, V., Martin-Pintado, J., & Gomez-Gonzalez, J. 1989, A&A, 218, 252 (BBMG89)
- Bachiller, R., Huggins, P. J., Cox, P., & Forveille, T. 1994, A&A, 281, L93 (BHCF94)
- Bryce, M., Balick, B., & Meaburn, J. 1994, MNRAS, 266, 721
- Guerrero, M. A., Manchado, A., & Chu, Y.-H. 1997, ApJ, 487, 328 (GMC97)
- Gurzadyan, G. A. 1997, The Physics and Dynamics of Planetary Nebulae, XVI, 513 pp. 179, Springer-Verlag
- Heyminck, S., Graf, U., Güsten, R., et al. 2012, A&A, this volume
- Hiriart, D. 2004, PASP, 116, 1135
- Hollenbach, D. J., & Tielens, A. G. G. M. 1999, Reviews of Modern Physics, 71, 173
- Kwok, S., Chong, S.-N., Koning, N., Hua, T., & Yan, C.-H. 2008, ApJ, 689, 219
- Liu, X.-W., Barlow, M.J., Cohen, M. et al. 2001, MNRAS, 323, 343
- Masson, C. R. 1990, ApJ, 348, 580
- O’Dell, C. R., Sabbadin, F., & Henney, W. J. 2007, AJ, 134, 1679 (ODell07)
- Sahai, R. 2000, ApJ, 537, L43
- Sahai, R. & Trauger, J. T. 1998, AJ, 116, 1357
- Sahai, R., Morris, M. R., & Villar, G. G. 2011, AJ, 141, 134
- Schneider, N., Simon, R., Kramer, C., et al. 2003, A&A, 406, 915
- Soker, N. 1996, ApJ, 460, L53
- Stacey, G. J., Hailey-Dunsheath, S., Ferkinhoff, C., et al. 2010, ApJ, 724, 957
- van Hoof, P. A. M., van de Steene, G. C., Barlow, M. J., et al. 2010, A&A, 518, L137