A New Radio Molecular Line Survey of Planetary Nebulae:
Key Words.:astrochemistry – ISM: molecules – planetary nebulae: individual (BD+303639, NGC 7027, NGC 6445, NGC 7008, NGC 6720, NGC 6853, NGC 6772, NGC 7293, NGC 6781) – radio lines: ISM
Certain planetary nebulae contain shells, filaments, or globules of cold gas and dust whose heating and chemistry are likely driven by UV and X-ray emission from their central stars and from wind-collision-generated shocks. We present the results of a survey of molecular line emission in the 88-236 GHz range from nine nearby (1.5 kpc) planetary nebulae spanning a range of UV and X-ray luminosities, using the 30 m telescope of the Institut de Radioastronomie Millimétrique. Rotational transitions of thirteen molecules, including CO isotopologues and chemically important trace species, were observed and the results compared with and augmented by previous studies of molecular gas in PNe. Lines of the molecules HCO, HNC, HCN, and CN, which were detected in most objects, represent new detections for five planetary nebulae in our study. Specifically, we present the first detections of CO (1-0, 2-1), HCO, CN, HCN, and HNC in NGC 6445; HCO in BD+303639; CO (2-1), CN, HCN, and HNC in NGC 6853; and CO (2-1) and CN in NGC 6772. Flux ratios were analyzed to identify correlations between the central star and/or nebular UV and X-ray luminosities and the molecular chemistries of the nebulae. This analysis reveals a surprisingly robust dependence of the HNC/HCN line ratio on PN central star UV luminosity. There exists no such clear correlation between PN X-rays and various diagnostics of PN molecular chemistry. The correlation between HNC/HCN ratio and central star UV luminosity demonstrates the potential of molecular emission line studies of PNe for improving our understanding of the role that high-energy radiation plays in the heating and chemistry of photodissociation regions.
Planetary nebulae (PNe) arise from outflowing stellar mass during the late evolutionary stages of intermediate-mass stars (0.8-8.0 M). These stars represent a significant stellar population in the Galaxy (Blöcker, 2001; Edwards et al., 2014). They progress through the main sequence and red giant branch (RGB) evolutionary phases, where core H and He fusion occur respectively, into the shell-burning asymptotic giant branch (AGB) phase. Once the star arrives there, slow AGB winds, originating in shocks and pulsations and driven by radiation pressure on dust, remove the bulk of the stellar envelope at mass loss rates of 10-10 M/yr (Zack & Ziurys, 2013; Blöcker, 2001). The presence of a close companion star can accelerate and otherwise profoundly affect this mass-loss process and the evolution of the resulting envelope (e.g., Marco & Izzard, 2017, and references therein).
When the AGB envelope depletes to a mass of 10 M, mass loss ceases and the star progresses to the post-AGB, and soon after, it becomes the central star of a PN (CSPN). Here, the cold, dusty AGB envelope is suddenly exposed to the hot (30-200 kK) post-fusion core’s ionizing UV and (in some cases) X-ray emission, which photodissociates and ionizes the envelope gas. This newly ionized gas constitutes the PN. Fast winds from the pre-white dwarf star (or its companion) may also slam into the slower moving ejected envelope, generating shocks that heat plasma to temperatures exceeding 10 K (Kastner et al., 2012, and references therein). Hot bubbles that form from this wind interaction have been found to be X-ray-luminous, thereby producing another ionization mechanism within the PN. Interactions from the winds that form the bubbles are also thought to continue shaping the nebula throughout its lifetime (Balick & Frank, 2002; Huarte-Espinosa et al., 2012). Atomic gas and dust within the nebula limit the penetration depth of incident UV photons from the CSPN, leaving the outer layers of the nebula insulated from them. Only higher energy (X-ray) photons, 0.5 keV, can penetrate the dense, neutral gas of the PN to ionize the cold, molecule-rich outer shells (Tielens & Hollenbach, 1985).
Millimeter CO and infrared H lines gave the first view of PN molecular gas 3 decades ago and have remained commonplace probes into the shells of ejected mass since (Zuckerman & Gatley, 1988; Huggins & Healy, 1989; Bachiller et al., 1991; Huggins et al., 1996, and references therein). Due to its high abundance, low excitation requirements and the low critical densities of its rotational transitions, CO is the most commonly observed and more widely reliable molecular species found in PNe (Huggins et al., 1996). Infrared lines of H are also observed from PNe, with detections of H mainly confined to bipolar (axisymmetric) nebulae (Kastner et al., 1996; Zuckerman & Gatley, 1988, and references therein). The complex interplay of PN central star radiation and composition of the proto-PN produce a rich environment for molecular chemistry that has served as motivation for various radio molecular line surveys (e.g., Bachiller et al., 1997; Zhang et al., 2008; Edwards et al., 2014; Schmidt et al., 2018). Whereas the foregoing surveys were largely restricted to well-studied objects, the recent molecular line surveys of Schmidt & Ziurys (2016, 2017a, 2017b) extended the PN sample coverage to younger objects.
In this paper we present molecular line surveys of nine PNe, obtained with the Institut de Radioastronomie Millimétrique 30 m telescope within the frequency range 88-236 GHz. These observations were specifically intended to explore the utility of potential tracers of high-energy irradiation of molecular gas. We report new detections of molecules and/or molecular transitions in five of these PNe, and place these results in the context of previously reported PN molecular line detections and measurements. We then evaluate the integrated flux ratios of observed molecular lines to study correlations between the molecular chemistry and high-energy radiation properties of the CSPNe. Specifically, we strive to find tracers of non-LTE chemistry due to X-irradiation of molecular gas and the effects of CSPN UV emission on the photodissociation regions (PDR) within PNe.
2 Observations and Data Reduction
|log(ergs s)||Source||log(ergs s)||(keV)|
PN data obtained from Montez Jr. et al. (2015) unless otherwise specified. X-Ray sources and median X-ray photon energies (E) from Kastner et al. (2012), except NGC 7027 and BD+30 from Montez Jr. et al. (2015). Upper limit X-ray luminosities for non-detections by Montez Jr. (unpublished).
P = point-like, D = diffuse, N = non-detection
The sample of nine PNe studied here were selected from the larger Chandra Planetary Nebula Survey PN sample (Kastner et al., 2012; Freeman et al., 2014). The ChanPlaNs survey targeted 35 objects within 1.5 kpc of Earth with the aim to identify and categorize PN X-ray emission, so as to understand the mechanisms giving rise to X-rays within PNe (Kastner et al. 2012). The ChanPlaNS survey mostly included high-excitation PNe and yielded X-ray detections from central star (point-like) and nebular (diffuse) X-ray sources.
We selected those ChanPlaNS PNe accessible to the IRAM 30 m telescope that had previously been detected333The lone exception, NGC 7008, was included in this survey due to its large X-ray luminosity. in CO and/or H. Physical properties of these nine PNe and their central stars are listed in Table 1. These properties include the morphology, angular radius in arcsec, distance, average radius, age, and ionized mass of each PN, along with the spectral type, temperature, and mass of the CSPN as compiled from the literature. The ultraviolet luminosity, X-ray luminosity, median energy of emission, and X-ray source properties for each of the planetaries are listed in Table 2.
Data were collected for the nine listed PNe during an observing run on the IRAM 30 m telescope on Pico Veleta from 5th to 10th June, 2012. The positions for each targeted PN are listed in Table 3 and illustrated in Figure 1. The observations were either centered on the PN central stars or the nebular rims; we obtained observations at two positions within NGC 6720, one toward the central star and one toward the rim. Frequencies in the regime of 88 GHz to 236 GHz were chosen to target key molecules for tracing chemistry. The molecular transitions detected during the course of the survey are listed in Table 4.
|Integration Time (hrs)|
|NGC 6720 Rim||185331.90||330125.00||2.98|
Integration times sampled from HCN, HCO, and HNC spectral region observations for each PN.
We used the Eight MIxer Receiver (EMIR), which is able to observe simultaneously in the 1 mm and 3 mm bands. Dual polarization mode was used whenever each specific frequency setting made it possible, and the lines from the two polarizations were averaged after checking that their intensities were compatible. We recorded the data with the FTS200 backend, providing a resolution of 200 kHz, which translates into velocity resolutions of 0.5 and 0.25 km s in the 88 and 236 GHz regimes, respectively, and a sufficient coverage of the lines listed above.
The beam sizes (Half Power Beam Widths) and beam efficiencies for the telescope ranged from 10.4” and 0.60, respectively, for the CO lines at 236 GHz to 27.8” and 0.80, respectively, for the HCN line at 88.6 GHz. Saturn and Mars were used to focus the telescope, and pointing errors were corrected to an accuracy of 3” by performing frequent measurements of nearby pointing calibrators (namely quasars such as QSO B1730-130 and QSO J0854+2006, but also Neptune, the ionized region K3-50A, and the planetary nebula NGC 7027). The subreflector was wobbled with a throw of 120 every 2 seconds to provide stable and flat baselines. The data were calibrated in units of main-beam temperature, T by frequent (15-20 min) calibration scans using the chopper wheel method. The T values so obtained were then re-scaled using sources with stable, well-known fluxes (CW Leo and NGC 7027). The estimated flux calibration accuracy is about 20%.
Integration times varied from an hour for CO lines to 6 hours for HCN, HCO, and HNC. Representative integration times are listed in Table 3. Weak lines were observed for significantly longer durations, up to 33 hours. System temperatures ranged between 102.4 K and 445.9 K, depending on the band and weather conditions.
Values of frequencies of hyperfine transitions obtained from the CDMS database. Line frequencies were measured in the laboratory by Müller et al. (2001), with updated catalog by Müller et al. (2005).
Asterisks indicate frequencies representative of hyperfine transition complexes.
Scan averaging and baseline subtraction of measurements were performed with the line analysis software available in CLASS and through targeted inspection. Baseline oscillations were observed in some higher frequency spectra, but did not significantly affect measurements of individual molecular transitions.
3 Results and Analysis
3.1 Molecular Line Measurements
Observations of the thirteen molecular lines detectable in the 88-236 GHz range provide new data for the nine PNe in our survey, including the well observed NGC 7027. We identified and compiled a list of these molecular transitions, with lines of primary interest and their measured intensities summarized in Tables LABEL:Tpeaks1 and LABEL:Tpeaks2, respectively, for objects that have been well observed in the past and those for which this survey yielded a significant number of new detections. The tables list peak flux (T) as defined by the main beam temperature, integrated line intensity (), reported parameters for lines detected in past studies, and the (1 ) formal errors for each measurement. The values of temperature, intensity, and their errors were obtained from the Gaussian fitting procedure in CLASS.
3.2 Survey Overview
We present the first detections of CO (1-0, 2-1), HCO, CN, HCN, and HNC in NGC 6445 (Section 3.3.6) as well as HCO in BD+303639 (Section 3.3.5). Transitions of CO (2-1), CN, HCN, and HNC were also detected for the first time in NGC 6853 (Section 3.3.8) as well as CO (2-1), CN, and HCO in NGC 6772 (Section 3.3.9). The individual PNe discussed in Sections 3.3.1-3.3.4 have been the subject of previous radio molecular line surveys (Bachiller et al., 1997; Edwards et al., 2014; Zhang et al., 2008), and the present survey provided only a few new detections for these objects. Among the molecular ions surveyed, HCO was detected in most of the PNe, including the aforementioned new detections in three objects. However, the only PN unambiguously detected in CO and NH was NGC 7027 (Section 3.3.1), which was previously known to exhibit emission from both molecules (Hasegawa & Kwok, 2001; Zhang et al., 2008).
Figures 2 through 14 illustrate the emission lines detected in each objects and compare their relative strengths. Across PNe with detected CN, intensities of the strongest line scale more or less consistently with lesser CN hyperfine lines. With this in mind, we only list parameters for the brightest CN (1-0) hyperfine components in Tables LABEL:Tpeaks1 and LABEL:Tpeaks2. Hydrogen recombination lines were also detected in two PNe, NGC 7027 and BD+303639, and these lines are displayed in Figure 5. Both objects are very young and compact, with high ionized gas densities capable of producing strong atomic hydrogen emission. All CO isotopologue lines for NGC 7008 (Section 3.3.7) were also previously unobserved, however the detections obtained here are likely due to ISM gas along the line of sight to the PN.
We find line strengths for the PNe in our survey compare favorably with previous observations (e.g., Bachiller et al., 1997). Because many of the objects were observed at different positions, we limit direct comparisons to a pair of representative PNe. In NGC 7027, we find line ratios of CO/CO and HCO/CO that are within 8% of those reported by Bachiller et al. (1997), while CN/CO is within 23%. Similarly, we noted a difference of only 15% in the reported CO/CO line ratio in NGC 6781, despite an offset of 50 between the positions observed within this PN in the two surveys.
3.3 Individual Planetary Nebulae
3.3.1 Ngc 7027
The data obtained for NGC 7027 yield by far the richest molecular line spectrum among the objects included in this survey. This young, bright nebula has a notably hot CSPN (=175,000 K, Kastner et al., 2012; Frew, 2008) and displays an abundance of atomic and molecular emission lines (Hasegawa & Kwok, 2001). With a dynamical age of only 700 years (Latter et al., 2000) and a high CSPN mass (0.7 M), NGC 7027 is a rapidly evolving nebula (Zhang et al., 2008). The morphology of the Hii region is that of a prolate ellipsoidal PN, with a distinct waist and symmetrical outflow lobes visible in the infrared and X-ray (Graham et al., 1993; Kastner et al., 2001). The clover-like symmetrical shell of molecular hydrogen forms a wispy and filamentary structure about the elongated shell of ionized gas (Latter et al., 2000; Cox et al., 2002). The H emission further traces out the thin PDR, separating the inner ionized gas from the expanding molecular envelope, which has been modeled spatio-kinematically by Santander-García et al. (2012) and determined to have a total molecular mass of roughly 1.3 . The presence of point-symmetric holes in the H region also indicates the presence of collimated outflows, a common phenomenon in young PNe (Cox et al., 2002).
All thirteen molecular transitions surveyed have been detected in NGC 7027 and their individual spectra are displayed in Figures 2 and 3, which illustrate detections of neutral and ionized molecules, respectively. From the shapes of the spectral emission lines, several characteristics of this PN can be discerned. The distinct double-peaked structure found in CO isotopologues and 1 mm lines traces the molecular gas with high line-of-sight velocities. Because the molecular shell of NGC 7027 extends to ¿70 (Hasegawa & Kwok, 2001), the 10 beam of the 30 m at 1 mm misses much of the outer envelope gas emitting near the systemic velocity. In contrast, the larger beam size of our 3 mm observations more completely samples the entire shell of the nebula, resulting in a larger contribution from material expanding along the plane of the sky and, hence, flatter-topped line profiles. The same explanation applies to the difference between the line profiles measured here and those measured by Zhang et al. (2008) with the larger beam of the ARO 12 m telescope. Indeed, our CO isotopologue spectral line shapes are consistent with 30 m observations by Herpin et al. (2002).
The central velocity of NGC 7027 was estimated by averaging the widths of each profile and assuming symmetric expansion. The CO spectral line profiles show blueshifted and redshifted edges at 10.7 and 36.9 km s, respectively. We thereby obtain a systemic velocity of 23.8 km s and an expansion velocity of 13.1 km s. These results are consistent with those of Graham et al. (1993), Cox et al. (2002), and Herpin et al. (2002), all of whom determined the systemic velocity of NGC 7027 to be approximately 25 km s and found molecular outflow velocities of 15 km s. The nitrogen-bearing molecules, however, present a slightly faster expanding shell, with blueshifted and redshifted outflow velocities of 5 and 40 km s, respectively. The velocity difference may be attributed to the N-bearing molecules being more abundant in the collimated, higher-velocity outflows present in NGC 7027.
We also observe evidence of self-absorption in the HCO spectrum and possibly in CO (1-0); on the blueshifted side of the HCO emission line, there appear absorption features at 0.8, 3.2, and 11.0 km s that dip below the averaged continuum, which has been baseline-subtracted (Figure 4). We believe these are circumstellar features and not interstellar as there is no observed extended emission in wide CO maps or self-absorption in higher CO J lines (Santander-García et al., 2012). Beyond the absorption features, an extended blueshifted wing is apparent within the HCO and CO (1-0) spectra (Figure 4). Extended CO wings were first detected in NGC 7027 by Bujarrabal et al. (2001). While the CO wing is less extensive, the HCO wing extends to -40 km s, indicating the presence of a molecular outflow with velocity 65 km s, i.e., nearly 5 times that of the overall molecular expansion velocity. This velocity is strictly a lower limit, dependent on the angle of observation. This suggests that low excitation, approaching regions of the outer molecular envelope are absorbing the line and continuum emission from molecules of high abundance in the inner regions of the PN.
We confirm the detection of the ion NH made by Zhang et al. (2008). Formation of NH could be a result of the strong presence of X-rays from the CSPN, or due to shocks in the pPN phase.
While the double-peaked structure that is characteristic of an expanding envelope of molecular gas can be seen in most of the emission lines observed, it is notably absent in HCN, due to its hyperfine splitting at closely spaced frequencies (88.630, 88.631, and 88.634 GHz). Line broadening from expansion blends these hyperfine features. Close hyperfine lines may also be blended and thus unresolved in CO, CO, and NH. CN displays hyperfine structure as well (Table 4), but at generally separable frequencies for both the N=2-1 and N=1-0 transitions. Two of the three hyperfine components of CO are detected (Figure 3).
Further analysis of the spectra revealed detections of hydrogen recombination lines at 92.0 (H41) and 222 (H38) GHz (Figure 5). H41 has a measured peak intensity of 0.278 K and integrated flux of 4.12 K km s, and was previously identified by Zhang et al. (2008). H38 has a peak intensity of 0.096 K and an integrated flux of 2.63 K km s.
3.3.2 Ngc 6720
Commonly known as the Ring Nebula, NGC 6720 (M57) has been well studied in molecular gas (Bachiller et al., 1997; Edwards et al., 2014). Its dynamical age is roughly 7,000 years (O’Dell et al., 2007) and it does not have detectable X-ray emission (Kastner et al., 2012). We observed two positions within NGC 6720, one toward the CSPN, and one just outside the optically bright nebula (Figure 1). Emission lines for these regions are displayed in Figures 6 and 7, respectively. Given the extensive CO survey of NGC 6720 performed by Bachiller et al. (1989a), we chose not to cover higher frequency wavebands (including three of the CO isotopologues) in our survey observations, and as such only six molecular transitions were detected.
The CSPN pointing (Figure 6) encompasses molecular gas associated with both the approaching and receding layers of the molecular shell, and as a result this pointing yields broad, double-peaked line profiles. We measure blueshifted and redshifted velocity components of -16.2 km s and 4.6 km s, respectively, at the central star of NGC 6720. These values agree well with the CO measurements of Edwards et al. (2014), though they also detect a redshifted component at 15.72.6 km s.
The spectra for the position offset from the CSPN yield line profiles that are dominated by emission near the systemic velocity of the PN, consistent with the expectation that the velocity of the bulk of the gas along the nebular rim is perpendicular to the line of sight. The extended tail toward the blueshifted side of the nebula indicates that the emission is dominated by the forward-directed side of the molecular region. At this rim position, emission peaks appear at -6.7 km s and -0.8 km s. Another peak was identified at 15.3 km s, but only in H-bearing molecules. This spectral feature in the Rim position corresponds to the receding emission component detected in molecular spectra obtained by Edwards et al. (2014). Although we did not observe CO (2-1) in NGC 6720, the observed profiles match the structure seen by Bachiller et al. (1989a) in their (=-40”, =-20”) position. Further, the average systemic velocity of NGC 6720 across both positions (-0.8 km s) agrees well with the measurements of Bachiller et al. (1997).
3.3.3 Ngc 7293
At a distance of only 0.20 kpc (Gaia Collaboration et al., 2018) and a radius of 0.46 pc, the Helix Nebula (NGC 7293) represents the largest and nearest PN in this study. It is a relatively evolved object with a dynamical age of 12,000 years (Schmidt & Ziurys, 2017a), and contains a hard, point-like X-ray source (Guerrero et al., 2001). Early radio molecular line observations were carried out by Huggins et al. (1992) and Bachiller et al. (1997), with many additional studies performed across the expansive molecular emission region of the PN (Young et al., 1999; Zack & Ziurys, 2013; Zeigler et al., 2013; Schmidt & Ziurys, 2017a; Schmidt et al., 2018).
Our observations of NGC 7293 were performed toward the edge of the molecular envelope (Figure 1). While all the molecular transitions studied in this paper have previously been detected in this PN, this work adds integrated flux values for these lines in a previously unobserved region of the object. As in previous single-dish mapping studies, there is significant structure present in the spectra (Figure 8). In particular, four distinct peaks appear in the CO (1-0, 2-1) transitions. This complex line profile can be attributed to multiple, distinct molecular knots along the line of sight sampled by the 30 m beam. Due to the large angular diameter of NGC 7293, the IRAM telescope beam measured only a small portion of the extensive molecular region. As such, it is only possible to make a comparison with a nearby region (=-435”, =75”), as observed by Zack & Ziurys (2013), Schmidt & Ziurys (2017a), and Schmidt et al. (2018). Their spectra revealed three CO velocity components at -11, -21, and -29 km s and -13, -21, -29 km s in HNC and HCN, respectively. This lines up well with our spectra, in which strong velocity components are observed at roughly -13, -23, and -32 km s, with two other potential features at -8 and -19 km s. Observations by Schmidt et al. (2018) of CN at the nearby position (-435, 75), as well as seven other positions across the Helix, detected hyperfine structure indicative of optically thin emission; our data are too noisy to confirm these results.
3.3.4 Ngc 6781
NGC 6781 is a highly evolved PN with an estimated dynamical age of at least 20,000 years (Ueta et al., 2014), making it the oldest object included in our survey. As one of our three surveyed PNe not detected in X-rays, NGC 6781 acts as a control for any potential X-ray-induced chemistry effects.
Radio molecular line observations of NGC 6781 were previously carried out by Bachiller et al. (1997). New observations of the PN presented here probe a different region of the expanding shell. NGC 6781 was detected in all molecular transitions except CO (Figure 9). The molecular line profiles of NGC 6781 display complex structures. At least four distinct velocity components were detected across the lines. Blueshifted and redshifted components at -0.3 and 29.4 km s respectively were detected for all molecular transitions. This suggests a systemic velocity of 14.6 km s. A central component at 14.7 km s was also ubiquitously detected, with a fourth line at 21.4 km s identified in most of the transitions. These velocity profiles generally match up with those presented in Bachiller et al. (1997), who found velocity extremes of CO from -5 to 39 km s.
The strong, narrow absorption component located at 6.5 km s in the CO line profiles is present across all detected transitions and is likely due to the presence of molecular cloud emission in the off-source (background sky) reference position.
Also known as Campbell’s Star, BD+303639 (hereafter BD+30), one of the youngest and most compact PN in our survey, has only just begun ionizing its envelope of ejected material within the past 1,000 years (Li et al., 2002). The Wolf-Rayet-type central star (WC9) has the coolest core of the PNe discussed in this paper. Nevertheless, BD+30 is a luminous object; at 10 ergs s, its CSPN UV luminosity is rivaled only by NGC 7027. It is also among the brightest diffuse PN X-ray sources (Kastner et al., 2012), and has been studied extensively across the electromagnetic spectrum (Freeman & Kastner, 2016). Still, it remains poorly characterized in molecular emission.
The structure of BD+30 is that of an elliptical nebula, with asymmetric CO bullets along possibly precessing jets (Bachiller et al., 1991; Akras & Steffen, 2011; Gussie & Taylor, 1995, and references therein). These misaligned jets may also be responsible for a hot bubble of X-ray-emitting material bounded by the ionized inner shell and outer dusty region (Akras & Steffen, 2011; Freeman & Kastner, 2016). An extensive multi-wavelength 3D structure modeling study was performed by Freeman & Kastner (2016).
In our survey data for BD+30, only HCO was detected (Figure 10). BD+30 subtends a diameter of 4” (Frew, 2008) and its CO emission is compact (Bachiller et al., 1991); hence, beam dilution may have a deleterious effect on CO and other molecules surveyed. The HCO line profile displays a fairly broad structure, with two distinct peaks. We find the expansion velocity of these regions within BD+30 to be -46.0 km s and 38.4 km s respectively. Compared with the observations by Bachiller et al. (1991), who determined CO bullet velocities of -63 km s and 41 km s, a puzzling discrepancy emerges. The redshifted HCO component matches the CO velocity, but the approaching velocity component differs by 17 km s, a significant offset. The velocity extension of the HCO emission suggests that the molecular ion is not contained solely within the CO bullets but spread across the PN. Future followup with interferometer observations will be necessary to confirm this interference.
As in the case of NGC 7027, examination of the spectra revealed the presence of strong H41 and H38 emission lines (Figure 5), with peak intensities 0.076 K, 0.033 K and integrated flux values of 1.02 K km s, 1.02 K km s, respectively. These H recombination lines will not be further discussed in this paper, but they are consistent with the presence of a large mass of ionized hydrogen in BD+30.
3.3.6 Ngc 6445
NGC 6445, or the Little Gem, is the most distant PN in this survey, at 1.38 kpc (Frew et al., 2016). It appears as a bipolar nebula of dynamical age roughly 5,000 years (Phillips, 1984) and its CSPN has a temperature of =170,000 K. NGC 6445 is at a comparable effective temperature and CSPN mass to NGC 7027, though it has far lower CSPN luminosity and larger nebula as a consequence of its more evolved state. The ChanPlaNs survey established that NGC 6445 harbors a compact X-ray source of uncertain physical origin, with coronal emission from a binary companion being a likely source (Montez Jr. et al., 2015). The CO isotopologues along with OH and OH were the only molecules previously detected in NGC 6445 (Huggins & Healy, 1989; Sun & Sun, 2000; Aleman et al., 2014).
Here, we present the detection of five new molecular species: CO, HCO, CN, HCN, and HNC. The shapes of these emission lines vary, with the nitrogen-bearing molecules displaying the broadest profiles (Figure 11). For those lines that are well-detected, double peaked structure is apparent and provides an estimated systemic velocity for the PN of 18.6 km s and outflow velocity of 22.0 km s. The lack of a blueshifted edge in the CO lines suggests weaker CO emission from the approaching side of NGC 6445. Nevertheless, these velocity measurements are comparable with those of Sun & Sun (2000) who measured a CO (1-0) systemic velocity of 20 km s and an expansion velocity of 33 km s. A feature near the systemic velocity at 20.4 km s can also be seen in the HCO and CN lines. The CO (1-0) transition also shows a dip in its emission line spectrum at 115.276 GHz that is likely the result of over-subtracted ISM molecular emission along the line of sight. There is some hint of NH emission at 0 km s, within the velocity range of the other well-detected emission lines, but given the absence of emission at 40 km s – where all the other lines peak – we do not regard this as a detection of NH.
3.3.7 Ngc 7008
NGC 7008 is a relatively young nebula, with a dynamical age of 5,700 years (Schmidt & Ziurys, 2016). Despite not being previously detected in molecular line emission, it was included in this survey for its strong X-ray emission; as in the case of NGC 6445, the CSPN is detected as a point-like X-ray source with an unconstrained emission origin (Montez Jr. et al., 2015). Previous molecular line measurements have been restricted to upper limits placed on HCO and HCN (Schmidt & Ziurys, 2016).
Observations were made along the northern edge of the molecular shell of NGC 7008 (Figure 1), yielding detections of narrow CO lines. Given the narrow widths and varying systemic velocities of the detected lines (Figure 12), we believe these emission lines of CO isotopologues originate from the ISM rather than from the nebula. Followup observations of NGC 7008 are necessary to establish the origin of CO emission towards the PN.
3.3.8 Ngc 6853
NGC 6853, also known as M27 or the Dumbbell Nebula, is a large PN at 0.74 pc in diameter that subtends 7,500 arcsec on the sky. In near-IR H emission, the Dumbbell Nebula consists of dense knots and ‘streamers’ that appear to radiate out from the CSPN (Kastner et al., 1996; Salama et al., 2000). Numerous atomic lines in the optical such as [O I] and [N II] further trace out the eponymous ‘dumbbell’ shape that reflect the PN’s bipolar structure (Edwards et al., 2014). In the ChanPlaNs survey, NGC 6853 was detected as a point-like source of X-ray emission, with intensity and low median energy as expected from photospheric emission due to the hot (=135,000 K) CSPN, given its small distance of 0.43 kpc (Frew et al., 2016).
CO was first detected in NGC 6853 by Huggins et al. (1996), who identified a complex filamentary structure much like that seen in H. Further observations were carried out by Salama et al. (2000) and more recently by Edwards et al. (2014). In the latter, the study focused on a position offset from the CSPN (=-68”, =-63”), where transitions of CO, HCO, and CS were identified. The region of the nebula targeted here is centered on the CSPN (Figure 1). H imaging revealed this line-of-sight intersects a large column of clumpy, molecular gas (Kastner et al., 1996).
In our 30 m spectra, a well defined double-peaked structure is observed, with blueshifted and redshifted components of -50.3 km s and 4.9 km s, respectively (Figure 13). When compared with the velocities of emission line components detected by Edwards et al. (2014) towards the edge of the PN (-35 km s and -5 km s, respectively) the velocity components we detect hint at the 3D structure of the PN. The measured velocity components indicate a systemic velocity of -22.7 km s and expansion velocity of 27.6 km s for the expanding molecular shell of NGC 6853.
3.3.9 Ngc 6772
One of the older objects in this sample, with a dynamical age of 11,000 years (Ali et al., 2012), the morphology of NGC 6772 is that of an elliptical shell with apparent distortion along the NE to SW axis. This structure may be a result of interactions between the expanding shell and the ISM (Schmidt & Ziurys, 2016). NGC 6772 is one of three PNe in this paper not detected in X-rays in the ChanPlaNS survey (Kastner et al., 2012). First detected in CO (1-0) by Zuckerman et al. (1990), Schmidt & Ziurys (2016) have more recently identified emission lines of CO (2-1) and HCN.
Our molecular line data were collected along the southern edge of the envelope (Figure 1). The 30 m spectra yield detections of both previously detected CO transitions, as well as CO (2-1), CN, and HCO (Figure 14). The broad peaks of the spectral line profiles terminate at a redshifted velocity of approximately 44.2 km s. This is consistent with the value determined by Schmidt & Ziurys (2016). Only CO (1-0) displays a narrow blueshifted peak at roughly 19.3 km s, as well as an apparent narrow feature near 0 km s that was also identified by Schmidt & Ziurys (2016). A possible explanation for this narrow feature is interstellar gas along the line of sight to NGC 6772.
|Object||Integrated Line Intensity||Uncertainty|
|(K km s)|
|NGC 6720 Rim||0.290||0.041|
Integrated line intensity values of the HNC/HCN ratio in observed PNe.
4.1 HNC/HCN Ratio
The HNC/HCN abundance ratio has been studied for the insight it provides into the physical conditions within the photodissociation regions of cold molecular clouds. Activation energy barriers and the existence of a critical temperature for the conversion between HNC and HCN, both formed predominantly by dissociative recombination of HCNH, have long been discussed (Schilke et al., 1992; Bachiller et al., 1997; Schmidt & Ziurys, 2017a). For semi-ionized gas near the PDR boundary, heightened temperatures should favor production and survival of HCN, decreasing the HNC/HCN ratio. The lack of convergence on the fundamental mechanisms and parameters that establish the ratio is apparent in the current literature, however (Graninger et al., 2014, and references therein).
Schilke et al. (1992) studied HNC/HCN as a diagnostic of local gas heating in the Orion Molecular Cloud (OMC-1), and in particular, explored the temperature dependence of the HNC+HHCN+H reaction. They found the HNC/HCN ratio to be 1 at temperatures of 10 K and noted that the ratio decreases in OMC-1 as the gas temperature increases from 10 to 150 K. Schilke et al. (1992) then assumed activation energies for both the HNC+H and HNC+O reactions of 200 K in order to reproduce the observed ratios. Using gas-grain and gas-phase chemical modeling, Graninger et al. (2014) furthered the study of HNC/HCN in OMC-1, where they identified the critical temperature range of the gas cloud at which the HNC/HCN ratio would begin to decrease from unity, i.e. 20-40 K, and placed upper limits on the activation barrier energy of well below 1200 K.
The conclusions from these foregoing observational studies are broadly consistent with those of Jin et al. (2015), who analyzed the HNC/HCN ratio across dark clouds, protostellar objects, and Hii regions. Jin et al. (2015) found the mean ratio decreased gradually across these evolutionary stages. They concluded that the HNC/HCN ratio can trace the evolutionary stages of massive star formation from cold molecular clouds (20 K) to the warmer molecular gas (100 K) associated with Hii regions generated by young, massive (OB) stars. Similarly, the ring-like vs. centrally peaked morphologies of HNC and HCN, respectively, within the disk orbiting TW Hya appear to place constraints on the characteristic temperature necessary to efficiently convert HNC into HCN at 25 K (Graninger et al., 2015).
The abundances of HNC and HCN in PNe can be expected to change due to thermal and kinetic processes as the envelope of AGB ejecta evolves to post-AGB and, with the CSPN reaching ionization temperatures, into PN stages (Schmidt & Ziurys, 2017a). As it is more stable, HCN is favored in the AGB stage where LTE chemistry dominates the inner envelope of the star (Schmidt & Ziurys, 2017a). Herpin et al. (2002) also observe lower abundances of HNC in post-AGB stages and claim this is due to increased ion-molecule interactions. By the proto-PN phase, abundances of HNC and HCN are expected to restabilize to unity (Schmidt & Ziurys, 2017a; Kimura et al., 2012). In more evolved PNe, models indicate the increased ionization fractions should result in enhanced abundances of HNC relative to HCN (Bachiller et al., 1997).
4.1.1 HNC/HCN Ratio with CSPN L
Eight PNe were observed in both the HNC and HCN during our molecular line survey. We find that the observed HNC/HCN line intensity ratio ranged from 0.03 to 0.72 within the objects sampled (Table 5).
Notably, this analysis has revealed a previously unknown correlation between the HNC/HCN line intensity ratio and the UV luminosity of the central star. In Figure 15, we plot the HNC/HCN ratios obtained from our 30 m data, as well as data from other nearby PNe taken from the literature, versus CSPN UV luminosity. The HNC/HCN line ratios for the PNe sample decrease from roughly unity to 0.027, as the CSPN UV luminosity increases. Linear regression for the data for our survey objects (omitting NGC 7027) supports the presence of a power-law relationship between the HNC/HCN ratio and CSPN UV luminosity, with a best-fit slope of and a correlation coefficient of . There exists no dependence of the HNC/HCN ratio on the type or presence of an X-ray source.
For the HNC/HCN line ratio to decrease steadily with UV luminosity, high energy photons must be either heating the gas within the PDR or selectively photodissociating HNC (Aguado et al., 2017). Under the former scenario, as CSPN UV luminosity rises, the resulting CSPN-induced ionization evidently raises the local temperature to exceed the activation energy for conversion of HNC into HCN+H. Graninger et al. (2015) suggest that when the local PDR is heated to the range of 25-40 K, the HNC+H reaction will lead to a decreased HNC/HCN ratio. Under this scenario, PDR heating would therefore be a direct result of photoelectrons generated by UV CSPN emission. Recent observations targeting proto-brown dwarfs observed a similar decrease in the HNC/HCN abundance ratio with bolometric luminosity for the case of protostars (Riaz et al., 2018). Riaz et al. (2018) suggest that this trend is due to the increased efficiency of heating of molecular gas, resulting in declining HNC abundance with increasing protostellar luminosity.
Another potential mechanism that might drive the HNC/HCN ratio is the difference in photodissociation rates of HNC and HCN. Calculations by Aguado et al. (2017), taking into account the photodissociation cross sections of each molecule as well as variation in the incident UV radiation field, find that selective photodissociation can increase the rate of destruction of HNC by a factor of 2-10 relative to that of HCN. While high effective temperature sources, such as CSPNe, have fewer low-energy photons and thus a lower effect on photodissociation rates than cool sources, the interaction between the strength of the evolving UV source and the photochemical reaction rates of HNC and HCN could be sufficient to drive the observed trend.
Modeling of HNC and HCN generation in cold, UV-irradiated molecular regions is necessary to identify the relevant processes that are most important in determining HNC/HCN variation in molecule-rich PNe. In addition, having established the correlation of HNC/HCN ratio with UV luminosity for an ensemble of PNe, it is necessary to investigate whether and how the HNC/HCN ratio depends on the local UV flux within individual objects.
4.1.2 HNC/HCN Ratio with Nebular Age
In the past, Herpin et al. (2002) suggested that the emergence of a strong UV source early in a PN’s history generates a significant drop in HCN with respect to HNC through PDR processing. Recent work on the HNC/HCN ratio in young PNe found that there is insignificant change with dynamical age (Schmidt & Ziurys, 2017a). These results were taken to indicate that the ratio does not vary significantly with nebular evolutionary state, and after the proto-PN stage, remains fixed. Because our survey samples PNe spanning a somewhat larger range of dynamical age, we have revisited the question of whether HNC/HCN correlates with PN age (Figure 16). Linear regression to our 30 m survey data (omitting NGC 7027) yields a correlation coefficient , suggesting a marginal correlation. When including data from the literature (pink and purple dots in Figure 16) (Schmidt & Ziurys, 2017a; Edwards & Ziurys, 2014, 2013), however, we find no correlation between HNC/HCN ratio and PN age. Nevertheless, the trend obtained from our data appears to match reasonably well with these previous PN survey data. These results support the notion that a continued rise in HNC with respect to HCN occurs throughout PN evolution. Due to the nature of PN age estimates, however, uncertainties in the ages are difficult to define and thus have not been included in the plot. Additional HNC and HCN observations of evolved PNe are necessary to confirm the tentative trend indicated in Figure 16.
4.2 Diagnostic of X-irradiated Gas
We investigated several potential tracers of X-ray-induced chemistry among the molecular emission lines. Line ratios were compared against X-ray and UV luminosity of each central star, as well as the ratio of the X-ray and UV luminosities. No significant trends emerged from this analysis. We focus here on HCO, which has long been proposed as a diagnostic of X-irradiation of molecular gas (Deguchi et al., 1990). Specifically, HCO is generated via X-ray ionization of molecular gas through the reaction CO+HHCO+H or by the molecular ion H through the reaction H+COHCO+H (Zhang et al., 2008). An enhanced abundance of HCO may then act as a tracer for X-ray emission in the PNe, and we would expect line ratios correlated with CSPN or nebular X-ray luminosity.
Figure 17 plots the ratio of HCO to CO against CSPN UV and X-ray emission. CO was chosen as it acts as a proxy of the molecular mass and is not coupled to X-irradiation. The absence of a clear trend in the left-hand panel of Figure 17 may suggest that X-rays are not directly driving HCO production in PNe. On the other hand, Figure 17 (center panel) also does not reveal a clear trend of enhanced HCO in the presence of strong UV from CSPNe. Hence, it remains to establish the irradiation source that is primarily responsible for the production of this molecular ion in PNe.
4.3 Co and NH
In PNe with hot central stars or high energy emission, X-rays have the potential to drive the chemical pathways available to the molecular gas. CO is produced in higher abundances in the presence of X-ray dominated regions through the reaction C+OHCO+H (Spaans & Meijerink, 2007). An enhanced abundance of CO then might suggest formation through this X-ray driven reaction. The main chemical pathway for the production of NH, i.e., H+NNH+H, also requires X-ray or cosmic ray irradiation as a means to generate the requisite H ion. (Zhang et al., 2008; Bell et al., 2007).
Of the nine observed PNe, transitions of CO and NH were detected only in the spectra of NGC 7027. None of the other objects surveyed here have had confirmed detections of CO, except NGC 6781, despite its potential importance in the generation of HCO (Bell et al., 2007). Considering the CO line was quite weak in NGC 7027 compared to CO, the molecule may be present in other X-ray emitting PNe, but at a level below our survey sensitivity. Given that NGC 7027 harbors an unusually luminous diffuse X-ray source and an unusually large mass of molecular gas, these detections hint at the viability of CO and NH diagnostics of X-irradiation of molecular gas. Shocks may also play a role in NH formation, however, and future observations are required to distinguish these two mechanisms.
This study presents observations of key molecular transitions across the 88-236 GHz range toward nine planetary nebulae with the IRAM 30 m radio telescope. We report new detections for five of the PNe in the survey across thirteen molecular lines included in our spectral coverage. Our survey yielded new detections of molecular transitions of CO, CN, HCN, HNC, and/or HCO in NGC 6445, NGC 6853, and NGC 6772. The molecular ion HCO was also detected in BD+303639 for the first time. In addition, confirmation of previously detected molecules were made in NGC 7027, NGC 6720, NGC 7293, and NGC 6781.
Our analysis of line ratios that are potential diagnostics of high-energy irradiation has revealed a strong but previously unrecognized anticorrelation between the HNC/HCN line intensity ratio and CSPN UV luminosity. This anticorrelation, which persists over nearly 3 orders of magnitude in central star UV luminosity, provides strong evidence that HNC/HCN acts as a tracer for heating of the nebular molecular gas by UV photons. Our survey additionally found a marginal correlation between the HNC/HCN line ratio and PN age.
Other irradiation tracers explored in our PN molecular line survey yield more ambiguous results. Though predicted to increase in abundance with X-irradiation, we find the intensity of HCO line emission in PNe shows no clear correlation with PN X-ray luminosity. Detections of CO and NH in our survey were limited to the particularly molecule-rich and luminous NGC 7027, which had been previously detected in lines of both species. This suggests that CO and NH trace especially intense high-energy irradiation of molecular gas.
Acknowledgements.This work was supported in part by a subcontract issued to RIT under NASA ADAP grant #80NSSC17K0057 to STScI (PI: B. Sargent), and by the Spanish MINECO within the program AYA2016-78994-P. It is also based on observations carried out with the IRAM telescopes. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). JB wishes to acknowledge useful discussions with Pierre Hily-Blant and Thierry Forveille during his Fall 2018 residency at IPAG under the support of the Chateaubriand Fellowship of the Office for Science & Technology of the Embassy of France in the United States.
- Aguado et al. (2017) Aguado, A., Roncero, O., Zanchet, A., Agúndez, M., & Cernicharo, J. 2017, The Astrophysical Journal, 838, 33
- Akras & Steffen (2011) Akras, S., & Steffen, W. 2011, Monthly Notices of the Royal Astronomical Society, 423, 925
- Aleman et al. (2014) Aleman, I., Ueta, T., Ladjal, D., et al. 2014, Astronomy & Astrophysics, 566, A79
- Ali et al. (2012) Ali, A., Sabin, L., Snaid, S., & Basurah, H. M. 2012, Astronomy & Astrophysics, 541, A98
- Bachiller et al. (1989a) Bachiller, R., Bujarrabal, V., Martin-Pintado, J., & Gómez-González, J. 1989a, Astronomy & Astrophysics, 218, 252
- Bachiller et al. (1997) Bachiller, R., Forveille, T., Huggins, P., & Cox, P. 1997, Astronomy & Astrophysics, 324, 1123
- Bachiller et al. (1991) Bachiller, R., Huggins, P., Cox, P., & Forveille, T. 1991, Astronomy & Astrophysics, 247, 525
- Bachiller et al. (1989b) Bachiller, R., Planesas, P., Martin-Pintado, J., Bujarrabal, V., & Tafalla, M. 1989b, Astronomy & Astrophysics, 210, 366
- Balick & Frank (2002) Balick, B., & Frank, A. 2002, Annual Review of Astronomy and Astrophysics, 40, 439
- Bell et al. (2007) Bell, T. A., Whyatt, W., Viti, S., & Redman, M. P. 2007, Monthly Notices of the Royal Astronomical Society, 382, 1139
- Blöcker (2001) Blöcker, T. 2001, Astrophysics and Space Science, 275, 1
- Bujarrabal et al. (2001) Bujarrabal, V., Castro-Carrizo, A., Alcolea, J., & Contreras, C. S. 2001, Astronomy & Astrophysics, 377, 868
- Cox et al. (2002) Cox, P., Huggins, P., Maillard, J.-P., et al. 2002, Astronomy & Astrophysics, 384, 603
- Deguchi et al. (1990) Deguchi, S., Izumiura, H., Kaifu, N., et al. 1990, The Astrophysical Journal, 351, 522
- Edwards et al. (2014) Edwards, J. L., Cox, E., & Ziurys, L. 2014, The Astronomical Journal, 791, 79
- Edwards & Ziurys (2013) Edwards, J. L., & Ziurys, L. M. 2013, The Astrophysical Journal Letters, 770, L5
- Edwards & Ziurys (2014) —. 2014, The Astrophysical Journal Letters, 794, L27
- Freeman et al. (2014) Freeman, M., Jr, R. M., Kastner, J. H., et al. 2014, The Astrophysical Journal, 794, 99
- Freeman & Kastner (2016) Freeman, M. J., & Kastner, J. H. 2016, The Astrophysical Journal Supplement Series, 226, 15
- Frew (2008) Frew, D. J. 2008, PhD thesis, Macquarie University, NSW, Australia
- Frew et al. (2016) Frew, D. J., Parker, Q. A., & Bojičić, I. S. 2016, Monthly Notices of the Royal Astronomical Society, 455, 1459
- Gaia Collaboration et al. (2018) Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1
- Górny et al. (1997) Górny, S. K., Stasińnska, G., & Tylenda, R. 1997, Astronomy & Astrophysics, 318, 256
- Graham et al. (1993) Graham, J. R., Serabyn, E., Herbst, T. M., et al. 1993, The Astronomical Journal, 105, 250
- Graninger et al. (2015) Graninger, D., Öberg, K. I., Qi, C., & Kastner, J. 2015, The Astrophysical Journal Letters, 807, L15
- Graninger et al. (2014) Graninger, D. M., Herbst, E., Öberg, K. I., & Vasyunin, A. I. 2014, The Astrophysical Journal, 787, 74
- Guerrero et al. (2001) Guerrero, M. A., Chu, Y.-H., Gruendl, R. A., Williams, R. M., & Kaler, J. B. 2001, The Astrophysical Journal, 553, L55
- Gussie & Taylor (1995) Gussie, G. T., & Taylor, A. 1995, Monthly Notices of the Royal Astronomical Society, 273, 801
- Hasegawa & Kwok (2001) Hasegawa, T. I., & Kwok, S. 2001, The Astrophysical Journal, 562, 824
- Herpin et al. (2002) Herpin, F., Goicoechea, J. R., Pardo, J. R., & Cernicharo, J. 2002, The Astrophysical Journal, 577, 961
- Huarte-Espinosa et al. (2012) Huarte-Espinosa, M., Frank, A., Balick, B., et al. 2012, Monthly Notices of the Royal Astronomical Society, 424, 2055
- Huggins et al. (1992) Huggins, P., Bachiller, R., Cox, P., & Forveille, T. 1992, The Astrophysical Journal, 401, L43
- Huggins et al. (1996) —. 1996, Astronomy & Astrophysics, 315, 284
- Huggins & Healy (1989) Huggins, P., & Healy, A. P. 1989, The Astrophysical Journal, 346
- Jin et al. (2015) Jin, M., Lee, J.-E., & Kim, K.-T. 2015, The Astrophysical Journal, 219, 1
- Kahane et al. (1992) Kahane, C., Cernicharo, J., Gómez-González, J., & Guélin, M. 1992, Astronomy & Astrophysics, 256, 235
- Kastner et al. (2012) Kastner, J. H., Montez Jr., R., Balick, B., et al. 2012, The Astronomical Journal, 144, 58
- Kastner et al. (2001) Kastner, J. H., Vrtilek, S. D., & Soker, N. 2001, The Astrophysical Journal, 550, L189
- Kastner et al. (1996) Kastner, J. H., Weintraub, D. A., Gatley, I., Merrill, K. M., & Probst, R. G. 1996, The Astrophysical Journal, 462, 777
- Kimura et al. (2012) Kimura, R. K., Gruenwald, R., & Aleman, I. 2012, Astronomy & Astrophysics, 541, A112
- Latter et al. (2000) Latter, W. B., Dayal, A., Bieging, J. H., et al. 2000, The Astrophysical Journal, 539, 783
- Li et al. (2002) Li, J., Harrington, J. P., & Borkowski, K. J. 2002, The Astronomical Journal, 123, 2676
- Marco & Izzard (2017) Marco, O. D., & Izzard, R. G. 2017, Publications of teh Astronomical Society of Australia, 34, 35
- Montez Jr. et al. (2015) Montez Jr., R., Kastner, J. H., Balick, B., et al. 2015, The Astrophysical Journal, 800, 8
- Müller et al. (2005) Müller, H. S. P., Schlöder, F., Stutzki, J., & Winnewisser, G. 2005, Journal of Molecular Structure, 742, 215
- Müller et al. (2001) Müller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, Astronomy & Astrophysics, 370, L49
- O’Dell et al. (2002) O’Dell, C. R., Balick, B., Hajian, A. R., Henney, W. J., & Burkert, A. 2002, The Astronomical Journal, 123, 3329
- O’Dell et al. (2007) O’Dell, C. R., Sabbadin, F., & Henney, W. J. 2007, The Astronomical Journal, 134, 1679
- Phillips (1984) Phillips, J. P. 1984, Astronomy & Astrophysics, 137, 92
- Riaz et al. (2018) Riaz, B., Thi, W.-F., & Caselli, P. 2018, Monthly Notices of the Royal Astronomical Society, 481, 4662
- Sabbadin et al. (1984) Sabbadin, F., Falomo, R., & Ortolani, S. 1984, Astronomy & Astrophysics, 137, 177
- Sahai et al. (2011) Sahai, R., Morris, M. R., & Villar, G. G. 2011, The Astronomical Journal, 141, 134
- Salama et al. (2000) Salama, A., Kessler, M. F., Leech, K., & Schulz, B., eds. 2000, An Infrared/Millimeter Study of the Dumbbell Nebula, ed. A. Salama, M. F. Kessler, K. Leech, & B. Schulz, Vol. 456, European Space Agency
- Santander-García et al. (2012) Santander-García, M., Bujarrabal, V., & Alcolea, J. 2012, Astronomy & Astrophysics, 545, A114
- Schilke et al. (1992) Schilke, P., Walmsley, C. M., des Forets, G. P., et al. 1992, Astronomy & Astrophysics, 256, 595
- Schmidt et al. (2018) Schmidt, D. R., Zack, L. N., & Ziurys, L. M. 2018, The Astrophysical Journal, 864, 31
- Schmidt & Ziurys (2016) Schmidt, D. R., & Ziurys, L. M. 2016, The Astrophysical Journal, 817, 175
- Schmidt & Ziurys (2017a) —. 2017a, The Astrophysical Journal, 835, 79
- Schmidt & Ziurys (2017b) —. 2017b, The Astrophysical Journal, 850, 123
- Spaans & Meijerink (2007) Spaans, M., & Meijerink, R. 2007, The Astrophysical Journal, 664, L23
- Stanghellini & Shaw (2008) Stanghellini, L., & Shaw, R. A. 2008, The Astrophysical Journal, 689, 194
- Sun & Sun (2000) Sun, J., & Sun, Y. 2000, Science in China, 43, 217
- Tielens & Hollenbach (1985) Tielens, A. G. G. M., & Hollenbach, D. J. 1985, The Astrophysical Journal, 291, 722
- Ueta et al. (2014) Ueta, T., Ladjal, D., Exter, K. M., et al. 2014, Astronomy & Astrophysics, 565, A36
- Young et al. (1999) Young, K., Cox, P., Huggins, P., Forveille, T., & Bachiller, R. 1999, The Astrophysical Journal, 522, 387
- Zack & Ziurys (2013) Zack, L. N., & Ziurys, L. M. 2013, The Astrophysical Journal, 765, 112
- Zeigler et al. (2013) Zeigler, N. R., Zack, L. N., Woolf, N. J., & Ziurys, L. M. 2013, The Astrophysical Journal, 778, 16
- Zhang et al. (2008) Zhang, Y., Kwok, S., & Dinh-V-Trung. 2008, The Astrophysical Journal, 678, 328
- Zuckerman & Gatley (1988) Zuckerman, B., & Gatley, I. 1988, The Astrophysical Journal, 324, 501
- Zuckerman et al. (1990) Zuckerman, B., Kastner, J. H., & Gatley, I. 1990, The Astrophysical Journal, 356, L59