Issue
A&A
Volume 542, June 2012
GREAT: early science results
Article Number L13
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201218929
Published online 10 May 2012

© ESO, 2012

1. Introduction

M17 SW is a giant molecular cloud at a distance of ~1.98 kpc (Xu et al. 2011), illuminated by a highly obscured (Av > 10 mag) cluster of several (among ≳ 100 stars) O/B stars (Beetz et al. 1976; Hanson et al. 1997), and it harbors several candidate young stellar objects (Povich et al. 2009). The large amount of observational data available in the literature and its nearly edge-on geometry make M17 SW one of the best-studied prototypes of a clumpy photon-dominated region (PDR) interface in the Galaxy. Studies of molecular and atomic emission indicate that the structure of the gas is highly clumped and not homogeneous (Stutzki et al. 1988; Stutzki & Güsten 1990; Meixner et al. 1992; Pérez-Beaupuits et al. 2010, and references therein), and the structure of its neutral and molecular gas seems to be dominated by magnetic rather than by thermal gas pressure, in contrast to many other PDR regions (Pellegrini et al. 2007). Temperatures of ~275 K were found toward the VLA 21cm continuum arc (Brogan & Troland 2001), and are associated with NH3 and highly excited 12CO emission (Güsten & Fiebig 1988; Harris et al. 1987).

Recent AKARI observations (Okada et al. 2010) with spatial resolutions between 39′′ and 57′′ showed that the [C II] 158 μm line emission is widespread in the M17 complex and peaks at the northern (M17 N) and southern (M17 SW) bars, in agreement with previous observations (Matsuhara et al. 1989; Stutzki et al. 1988). However, this and previous observations of the hot and ionized gas in M17 SW are limited in spatial resolution and extent (e.g. Harris et al. 1987; Stutzki et al. 1988; Genzel et al. 1988; Meixner et al. 1992; Howe et al. 2000). Therefore, in this work we present new observations (of an area 5.7×3.7\hbox{$5\farcm7 \times 3\farcm7$}) of hot molecular (12CO J = 13–12) and ionized atomic ([C II] 158 μm) gas, with spatial resolutions of ~19.8\hbox{$19\farcs8$} and ~15.6\hbox{$\farcs6$}, respectively, which advances existing work in M17 SW.

thumbnail Fig. 1

Left panels – intensity maps of [C II] 158 μm (top) and 12CO J = 13–12 (bottom) in M17 SW, integrated in the velocity ranges 0–40 km  s-1 and 14–28 km  s-1, respectively. The contour levels are the 10%, 25%, 50% (thick line), 75% and 90% of the peak emission. Dashed lines are the strip lines shown in Fig. 2. The central position (0′′, 0′′) is marked with a cross. The ultracompact H II region M17-UC1 and four H2O masers (Johnson et al. 1998) are marked by the black circle and plus symbols, respectively. Right panel – 21 cm continuum emission by Brogan & Troland (2001) with the overlaid contours of the velocity integrated (same as above) emission of 12CO J = 13–12 (green), [C II] (red), and the [C I] P233P1\hbox{$^3P_2-{^3P_1}$} 370 μm (blue, integrated in 9–30 km  s-1) from Pérez-Beaupuits et al. (2010). The contour levels (from thin to thick) are the 50%, 75% and 90% of the peak emission. The stars indicate the O and B ionizing stars (Beetz et al. 1976; Hanson et al. 1997). Dashed frames depict the beam center for the edges of the 6 OTF strips. Contour maps are smoothed to 20′′ resolution.

2. Observations

The observations were performed with the German REceiver for Astronomy at Terahertz frequencies (GREAT1, Heyminck et al. 2012) on board the Stratospheric Observatory For Infrared Astronomy (SOFIA). We used the dual-color spectrometer during its first Short Science flight on 05 April 2011 to simultaneously map the fine-structure transition of [C II] at 1900.536900 GHz (157.7 μm) and the 12CO J = 13–12 transition at 1496.922909 GHz (200.3 μm) toward M17 SW. The observations were performed in on-the-fly (OTF) total power mode. The area mapped consists of six strips, each covering 224′′ × 32′′ (Δα × Δδ) with a sampling of 8′′ (half the beamwidth at 1.9 THz). Hence, each strip consists of four OTF lines containing 28 points each. We integrated 1s per dump and 5 s for the off-source reference.

All our maps are centered on RA(J2000) = 18:20:27.6 and Dec(J2000) = –16:12:00.9, which corresponds to the SAO star 161357. For better system stability we used a nearby reference position at offset (345′′, − 230′′). A pointed observation of this reference position against the known better reference (offset (1040′′, −535′′) Matsuhara et al. 1989) showed that the reference is clean of 12CO emission, but contains weak (<20% of the peak emission) and relatively narrow (FWHM ~ 15 km  s-1) [C II] emission. All [C II] spectra presented here were corrected for this missing flux.

Pointing was established with the SOFIA optical guide cameras, and was accurate to better than 10′′. Because the acousto-optical and fast Fourier transform spectrometers (Klein et al. 2012) operated in parallel give redundant information, we base the following analysis on the data from the latter, which provided 1.5 GHz bandwidth with about 212 kHz (~0.03 km  s-1) of spectral resolution. Because during this commissioning flight the instrument showed some random gain variations, we discarded 15–20% of the spectra by filtering the spectra for outliers in the total-power IF level and for outliers in the noise rms (obtained after subtracting a third-order polynomial), retaining only reliable data. The calibration of these data to antenna temperature was performed with the kalibrate task from the kosma_software package (Guan et al. 2012). We then reduced and imaged the data further with the CLASS90 packages, which is part of the GILDAS2 software. Using the beam efficiencies (ηc) 0.51 for [C II] and 0.54 for 12CO J = 13–12, and the forward efficiency (ηf) of 0.95 (Heyminck et al. 2012), we converted all data to main beam brightness temperature scale, TB=ηf×TA/ηc\hbox{$T_{\rm B}=\eta_{\rm f}\times T_{\rm A}^{*}/\eta_{\rm c}$}.

thumbnail Fig. 2

Left – strip lines of the molecular (12CO) and atomic ([C I], [C II]) gas at Δδ = 30′′ (PA 90°) across the ionization front of M17 SW. Right – strip lines at PA 63° (Δδ = Δα/2). The lower-J lines of 12CO, and [C I] are from Pérez-Beaupuits et al. (2010, their Fig. 3). All these profiles are in units of K  km  s-1 . The VLA 21 cm continuum and H I optical depth (integrated between 0 and 30 km  s-1) by Brogan & Troland (2001) are in units of 500 × Jy beam-1 and 8 × τ(HI) km  s-1 , respectively. The offset, Δα = 0′′ in RA, is the same as in Fig. 1.

3. Results

The velocity-resolved spectra of the new GREAT/SOFIA observations reveal a much more complex structure than expected based on our understanding of M17 SW as a highly clumped PDR, with a significant fraction of the [C II] emission not at all associated with the molecular PDR material.

Figures 1–4 display the GREAT data compared with supplementary data of typical PDR tracers from the literature. Figure 1 shows the velocity-integrated emission of [C II] between 0 and 40 km  s-1 (peak ~ 1070 K  km  s-1 ) and 12CO J = 13–12 [14–28 km  s-1, peak ~ 180 K  km  s-1] (left panels). While the high-excitation CO emission follows the hot dense gas structure seen in previous submm line studies (but see the discussion below), the [C II] distribution reveals surprises: its velocity-integrated distribution is quite shallow and far more extended than predicted for the stratification in classical PDR models. This is demonstrated in Fig. 1a and the intensity cuts across the PDR in Fig. 2, showing strong [C II] emission over the whole extent of the mapped area.

The broader distribution is consistent with the lower angular resolution data from the KAO (Stutzki et al. 1988), at the time interpreted to result from the deeper UV-penetration into a clumpy medium. However, the new GREAT/SOFIA data, resolving the spectra in velocity, show a more complex scenario. We find that (in M17 SW) the [C II] emission (Eu = 91.21 K, and ncrit ~ 5 × 103 cm-3) traces – in addition to the dense PDR gas west of the ionization front – also a widespread (toward the east), probably more diffuse atomic gas component that is invisible in CO. A sample of this is shown in Fig. 3 (left panel), displaying in gray scale the intensity distribution of the lower velocity components of [C II] and its (non)correlation with CO and H I. Comparing individual spectra (Fig. 4), we see that the span of observed velocities in [C II] is much wider than that of any other PDR tracing species. In the spectra at offset (0, 0), a position well in front of the PDR interface (Fig. 1b), for example, [C II] still shows a bright and amazingly wide line with ΔV ~ 30 km  s-1, including velocities not detected in molecular (CO) or (denser) atomic (neutral carbon [C I]) gas column density tracers. Most interestingly, those velocities (e.g., 4–10 km  s-1) exhibit appreciable opacities in atomic hydrogen (Brogan & Troland 2001 – note that H I is observed in absorption against the 21 cm background continuum, hence its detectability across the map depends on the continuum distribution). The broad and intricate structure of the [C II] spectral line (Fig. 4) can also be the signature of several clumps or layers of ionized carbon gas at “redder” and “bluer” velocities than the other tracers.

Figure 1b provides evidence for a C +  →  C stratification across the PDR: the bulk of the [C I] emission (defined as >75% of the peak emission) avoids the bulk emission of [C II], with their peaks being separated by ~70′′ in the NE-SW direction, i.e., 0.67 pc at the distance of ~1.98 kpc (Xu et al. 2011). In a previous analysis the extended [C I] emission (west from the ionization front), and its particular peak emission in front of the bulk CO emission (if stratified, the [C I] peak should be observed before the CO peak, from the direction of the PDR interface), was argued to emerge from the interclump regions of a very clumpy medium, but could also result from a partial face-on illumination of the molecular clouds (e.g., Stutzki et al. 1988; Stutzki & Güsten 1990; Meixner et al. 1992; Pérez-Beaupuits et al. 2010). In the context of the clumpiness of the PDR and the newly discovered association of much of the [C II] emission with atomic gas (see below) at velocities not matched by any molecular material, the [C II] peak cannot be interpreted as the classical stratification expected in a homogeneous PDR3, but instead has to be caused by excitation gradients with increasing optical depth into the atomic and molecular cloud.

The excitation study of the warm PDR layers, including the CO excitation, will be the subject of a forthcoming publication. Here we briefly elaborate on the new GREAT data: the 12CO J = 13–12 line follows a similar spatial distribution as the lower-J lines reported in Pérez-Beaupuits et al. (2010), although the whole 12CO J = 13–12 emission is shifted toward the ionization front. Its peak emission is ~16′′ (0.15 pc) southeast of the 12CO J = 6–5 peak. The J = 13–12 line is already fainter than the lower-J transitions, indicating subthermal excitation.

thumbnail Fig. 3

Selected 1 km  s-1 wide channel maps (gray scale) of the [C II] velocity-integrated (K  km  s-1) in the velocity ranges 9.5–10.5 km  s-1 and 18.5–19.5 km  s-1. The contours show the integrated emission in the corresponding velocity channels of [C I], 12CO J = 2–1, and τ(HI) in 20% steps (from thin to thick lines) of the peak channel integrated values (top right in the maps). All maps but τ(H I) were smoothed to a resolution of 20′′.

4. Discussion

We can distinguish two different spatial regimes in the PDR with embedded star cores of active star formation. The strip line along Δδ = 30′′ (PA 90°, Fig. 2A) traces the embedded star-formation, going through the ultracompact H II region UC1 (Δα = −30′′), the embedded H2O maser positions farther west, and the weaker additional H II region seen as a separate clump at Δα = −110′′ in the 21 cm continuum map and matching the secondary peak of [C II], [C I], and the 12CO J = 2–1 lines (Figs. 1b and 2A). A dominant internal heating is likely suppressing the layering here. The strip with the position angle 63° (Fig. 2B) instead traces the classical PDR. For a homogeneous PDR one expects a stratified layering of [C II], [C I], high-J CO, and low-J CO. However, Stutzki et al. (1988) showed the structure to be very clumpy so that no stratification should be observed. Our observations now show that high-J CO and low-J CO peak at similar depths in the cloud while the integrated [C II] peak is shifted by ≈ 0.19 pc relative to the CO peaks. The [C II] shift cannot even be explained by homogeneous PDR models, much less by a clumpy medium.

thumbnail Fig. 4

Line profiles at selected positions along our cut across the PDR at PA 63° (Δδ = Δα/2). All data but τ(H I) were smoothed to a spatial and a spectral resolution of 20′′ and ~0.6 km  s-1, respectively.

In [C II] we detect all velocity components that were discussed in the H I channel maps by Brogan & Troland (2001, their Figs. 4, including a weak foreground cloud at 7 km  s-1, shocked gas in front of the cluster that is blown toward the observer at 11–17 km  s-1, and the main PDR velocity of ~21 km  s-1. In the higher density gas (west of the H II region), [C II] seems to be optically thick, showing self-absorption dips at the peak velocities of H I and the molecular lines. The 12CO J = 2–1 emission also traces the optically thin shocked gas down to 10 km  s-1 (Fig. 4), the 12CO J = 13–12 line is hardly excited there. A good match between [C II] and 12CO J = 2–1 can be seen at intermediate velocities (Fig. 3, right panel) where the shifted [C I] emission with respect to [C II] is also evident.

We can compare the column density contributions visible in the different velocity components using the column density for the H I gas by Brogan & Troland (2001) and the C+ column density (e.g., Schneider et al. 2003, their Eq. (A.4)) responsible for the observed [C II]emission. For the gas temperature and the density we assumed 250 K and 104 cm-3 at the offset (0, 0), and 150 K and 106 cm-3 for the offsets deeper in the cloud according to Pérez-Beaupuits et al. (2010).

To compare our results with the H I data we estimated [C II] and H I column densities at selected positions along the southern PDR (Fig. 4) from the [C II] emission and H I optical depths, integrating the spectra in the velocity ranges 11–17 km  s-1 and 17–24.5 km  s-1. Selection of these velocity ranges is described by Brogan et al. (1999). In the range 11–17 km  s-1 we found N([C II]) (with uncertainties of ~20%) between 6.5 × 1017 cm-2 at offset position (0,0) and 1.6 × 1018 cm-2 at (−60, −30), while the corresponding H I column densities vary (with uncertainties 7–16%) between 5.5 × 1021 cm-2 and 4.7 × 1021 cm-2, respectively. In the velocity range 17–24.5 km  s-1 the [C II] column densities are between 1.2 × 1018 cm-2 and 2.9 × 1018 cm-2 at the same offset positions, and we found N(H I) between 3.7 × 1021 cm-2 and 4.4 × 1021 cm-2. The variation in H I column densities is less than 20% in both velocity ranges, indicating a relatively homogeneous H I interface seen also by [C II]. The latter, instead, presents a much stronger column variation (40%–60%) depending on the ambient conditions.

Assuming the normal cosmic carbon abundance of X(C/H) = 1.2 × 10-4 (Wakelam & Herbst 2008) and complete ionization of the carbon, we obtained gas column densities N(H) = (5.4−13) × 1021 cm-2 for the low-velocity component and N(H) = (10−24) × 1021 cm-2 for the high-velocity component. The [C II] emission of the low-velocity component in front of the PDR is therefore consistent with pure H I gas, while for all other components (above 17 km  s-1) at least half of the [C II] emission stems from molecular gas, i.e., classic PDR material. The gas in front of the southern PDR must be atomic to a high degree, indicating non-equilibrium chemistry (e.g. Stoerzer & Hollenbach 1998).

5. Final remarks

The integrated intensity maps do not tell us about the origin of the emission lines. Only narrow channel maps of the velocity-resolved [C II] spectra show a correlation with the diffuse gas components. From the comparison of different tracers in different velocity components we found that for the regions of embedded star formation the internal heating is always the dominant process, providing a very good match of all tracers without significant layering. The good spatial correlation of CO, [C II], and other tracers at particular velocity components/channels confirms the clumpy PDR picture from Stutzki & Güsten (1990). However, a significant fraction of the [C II] velocity components is not associated (spatially) with the dense (high-J CO) or the diffuse (H Iand CO J = 2–1) PDR tracers. The structure neither matches a stratified nor a clumpy PDR model. The [C II] emitting gas is only partially explained by atomic gas. Optically thin, hot molecular gas must contribute more than half of the emission.


1

GREAT is a development by the MPI für Radioastronomie and the KOSMA/ Universität zu Köln, in cooperation with the MPI für Sonnensystemforschung and the DLR Institut für Planetenforschung.

3

In a clumpy PDR, with many internal surfaces, the [C II] peak should be even closer to the [C I] peak, becoming indistinguishable from it depending on clumpiness of the source and spatial resolution of the observations.

Acknowledgments

We are grateful to C. Brogan for providing the 21 cm and H I data. We thank the referee for the very careful reading of the manuscript and constructive comments. Results are partially based on observations made with the NASA/DLR Stratospheric Observatory for Infrared Astronomy. SOFIA Science Mission Operations are conducted jointly by the Universities Space Research Association, Inc., under NASA contract NAS2-97001, and the Deutsches SOFIA Institut under DLR contract 50 OK 0901. We also gratefully acknowledge the outstanding support by the observatory staff. This work was supported by the German Deutsche Forschungsgemeinschaft, DFG project SFB 956C.

References

  1. Beetz, M., Elsaesser, H., Weinberger, R., & Poulakos, C. 1976, A&A, 50, 41 [NASA ADS] [Google Scholar]
  2. Brogan, C. L., & Troland, T. H. 2001, ApJ, 560, 821 [NASA ADS] [CrossRef] [Google Scholar]
  3. Brogan, C., Troland, T., Roberts, D., & Crutcher, R. 1999, ApJ, 515, 304 [NASA ADS] [CrossRef] [Google Scholar]
  4. Genzel, R., Harris, A. I., Stutzki, J., & Jaffe, D. T. 1988, ApJ, 332, 1049 [NASA ADS] [CrossRef] [Google Scholar]
  5. Guan, X., Stutzki, J., Graf, U., et al. 2012, A&A, 542, L4 [Google Scholar]
  6. Güsten, R., & Fiebig, D. 1988, A&A, 204, 253 [NASA ADS] [Google Scholar]
  7. Hanson, M. M., Howarth, I. D., & Conti, P. S. 1997, ApJ, 489, 698 [NASA ADS] [CrossRef] [Google Scholar]
  8. Harris, A. I., Stutzki, J., Genzel, R., et al. 1987, ApJ, 322, L49 [NASA ADS] [CrossRef] [Google Scholar]
  9. Heyminck, S., Graf, U., Güsten, R., et al. 2012, A&A, 542, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Howe, J. E., Ashby, M. L. N., Bergin, E. A., et al. 2000, ApJ, 539, L137 [NASA ADS] [CrossRef] [Google Scholar]
  11. Johnson, C. O., Depree, C. G., & Goss, W. M. 1998, ApJ, 500, 302 [NASA ADS] [CrossRef] [Google Scholar]
  12. Klein, B., Hochgürtel, S., Krämer, I., et al. 2012, A&A, 542, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Matsuhara, H., Nakagawa, T., Shibai, H., et al. 1989, ApJ, 339, L67 [NASA ADS] [CrossRef] [Google Scholar]
  14. Meixner, M., Haas, M., Tielens, A., Erickson, E., & Werner, M. 1992, ApJ, 390, 499 [NASA ADS] [CrossRef] [Google Scholar]
  15. Okada, Y., Kawada, M., Murakami, N., et al. 2010, A&A, 514, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Pellegrini, E. W., Baldwin, J. A., Brogan, C. L., et al. 2007, ApJ, 658, 1119 [NASA ADS] [CrossRef] [Google Scholar]
  17. Pérez-Beaupuits, J. P., Spaans, M., Hogerheijde, M. R., et al. 2010, A&A, 510, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Povich, M. S., Churchwell, E., Bieging, J. H., et al. 2009, ApJ, 696, 1278 [NASA ADS] [CrossRef] [Google Scholar]
  19. Schneider, N., Simon, R., Kramer, C., et al. 2003, A&A, 406, 915 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Stoerzer, H., & Hollenbach, D. 1998, ApJ, 495, 853 [NASA ADS] [CrossRef] [Google Scholar]
  21. Stutzki, J., & Güsten, R. 1990, ApJ, 356, 513 [NASA ADS] [CrossRef] [Google Scholar]
  22. Stutzki, J., Stacey, G. J., Genzel, R., et al. 1988, ApJ, 332, 379 [NASA ADS] [CrossRef] [Google Scholar]
  23. Wakelam, V., & Herbst, E. 2008, ApJ, 680, 371 [NASA ADS] [CrossRef] [Google Scholar]
  24. Xu, Y., Moscadelli, L., Reid, M. J., et al. 2011, ApJ, 733, 25 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1

Left panels – intensity maps of [C II] 158 μm (top) and 12CO J = 13–12 (bottom) in M17 SW, integrated in the velocity ranges 0–40 km  s-1 and 14–28 km  s-1, respectively. The contour levels are the 10%, 25%, 50% (thick line), 75% and 90% of the peak emission. Dashed lines are the strip lines shown in Fig. 2. The central position (0′′, 0′′) is marked with a cross. The ultracompact H II region M17-UC1 and four H2O masers (Johnson et al. 1998) are marked by the black circle and plus symbols, respectively. Right panel – 21 cm continuum emission by Brogan & Troland (2001) with the overlaid contours of the velocity integrated (same as above) emission of 12CO J = 13–12 (green), [C II] (red), and the [C I] P233P1\hbox{$^3P_2-{^3P_1}$} 370 μm (blue, integrated in 9–30 km  s-1) from Pérez-Beaupuits et al. (2010). The contour levels (from thin to thick) are the 50%, 75% and 90% of the peak emission. The stars indicate the O and B ionizing stars (Beetz et al. 1976; Hanson et al. 1997). Dashed frames depict the beam center for the edges of the 6 OTF strips. Contour maps are smoothed to 20′′ resolution.

In the text
thumbnail Fig. 2

Left – strip lines of the molecular (12CO) and atomic ([C I], [C II]) gas at Δδ = 30′′ (PA 90°) across the ionization front of M17 SW. Right – strip lines at PA 63° (Δδ = Δα/2). The lower-J lines of 12CO, and [C I] are from Pérez-Beaupuits et al. (2010, their Fig. 3). All these profiles are in units of K  km  s-1 . The VLA 21 cm continuum and H I optical depth (integrated between 0 and 30 km  s-1) by Brogan & Troland (2001) are in units of 500 × Jy beam-1 and 8 × τ(HI) km  s-1 , respectively. The offset, Δα = 0′′ in RA, is the same as in Fig. 1.

In the text
thumbnail Fig. 3

Selected 1 km  s-1 wide channel maps (gray scale) of the [C II] velocity-integrated (K  km  s-1) in the velocity ranges 9.5–10.5 km  s-1 and 18.5–19.5 km  s-1. The contours show the integrated emission in the corresponding velocity channels of [C I], 12CO J = 2–1, and τ(HI) in 20% steps (from thin to thick lines) of the peak channel integrated values (top right in the maps). All maps but τ(H I) were smoothed to a resolution of 20′′.

In the text
thumbnail Fig. 4

Line profiles at selected positions along our cut across the PDR at PA 63° (Δδ = Δα/2). All data but τ(H I) were smoothed to a spatial and a spectral resolution of 20′′ and ~0.6 km  s-1, respectively.

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.