Submillimeter Imaging of G5.89-0.39

Subarcsecond Submillimeter Imaging of the Ultracompact HII Region G5.89-0.39


We present the first subarcsecond submillimeter images of the enigmatic ultracompact HII region (UCHII) G5.89-0.39. Observed with the SMA, the 875 m continuum emission exhibits a shell-like morphology similar to longer wavelengths. By using images with comparable angular resolution at five frequencies obtained from the VLA archive and CARMA, we have removed the free-free component from the 875 m image. We find five sources of dust emission: two compact warm objects (SMA1 and SMA2) along the periphery of the shell, and three additional regions further out. There is no dust emission inside the shell, supporting the picture of a dust-free cavity surrounded by high density gas. At subarcsecond resolution, most of the molecular gas tracers encircle the UCHII region and appear to constrain its expansion. We also find G5.89-0.39 to be almost completely lacking in organic molecular line emission. The dust cores SMA1 and SMA2 exhibit compact spatial peaks in optically-thin gas tracers (e.g. SO), while SMA1 also coincides with 11.9m emission. In CO(3-2), we find a high-velocity north/south bipolar outflow centered on SMA1, aligned with infrared H knots, and responsible for much of the maser activity. We conclude that SMA1 is an embedded intermediate mass protostar with an estimated luminosity of 3000 L and a circumstellar mass of  M. Finally, we have discovered an NH (3,3) maser 12″ northwest of the UCHII region, coincident with a 44 GHz CHOH maser, and possibly associated with the Br outflow source identified by Puga et al. (2006).

HII regions — stars: formation — infrared: stars — ISM: individual (W28 A2) — ISM: individual (G5.89-0.39) — ISM: individual (IRAS 17574-2403) — techniques: interferometric — submillimeter

1 Introduction

The formation of massive stars is a fundamental process in astrophysics that is not well understood. The identification of ultracompact HII regions (UCHIIs) by Wood & Churchwell (1989) provided an important step forward in the study of massive star formation. These objects are less than  cm in diameter, have a centimeter wavelength spectral index consistent with free-free emission, and are associated with far-infrared and submillimeter sources of high bolometric luminosity (Hunter et al., 2000; Thompson et al., 2006). Based on these characteristics, UCHIIs are believed to be ionized nebulae powered by young OB stars located at (or near) their centers and surrounded by high column densities of dust and molecular gas (Churchwell, 2002). However, because of the high extinction at optical and infrared wavelengths, the exciting star is often too deeply embedded for its photospheric emission to be observed. In fact, until recently the exciting stars had been identified for only a few percent of UCHII regions (Hanson et al., 2002), with the first (and perhaps best) example being G29.96-0.02 (Watson et al., 1997; Martín-Hernández et al., 2003). Recent infrared studies have identified candidate ionizing sources toward several more UCHIIs (Alvarez et al., 2004; Apai et al., 2005; Comerón et al., 2006; Bik et al., 2006).

One UCHII for which a candidate ionizing star has been proposed is W28 A2, also known as G5.89-0.39, which is a well-known shell-type UCHII. Feldt et al. (2003) have identified a near-infrared star located along the northeastern rim of the radio shell which they argue is of spectral type O5 or earlier (hereafter called Feldt’s star). G5.89-0.39 is associated with a strong bipolar molecular outflow originally identified in CO by Harvey & Forveille (1988) and having one of the largest mass outflow rates known (Churchwell, 1997). The outflow emission has been subsequently studied using single-dish CO and SiO (Klaassen et al., 2006; Choi et al., 1993; Acord et al., 1997) and interferometric CO, HCO (Watson et al., 2007) and SiO (Acord et al., 1997; Sollins et al., 2004) observations. The reported position angle of the outflow is notably different between the tracers, being nearly east-west in CO and HCO (position angle +84°), vs. northeast-southwest in SiO (position angle +30°). While masers of many kinds are seen toward and around G5.89-0.39 (the most recent images being: Hofner & Churchwell (1996) (HO), Stark et al. (2007) and Fish et al. (2005) (OH), and Kurtz et al. (2004) (Class I CHOH)), the spatial-velocity gradients seen do not lead to a clear consensus view on the origin of the outflow. Furthermore, despite the fact that the purported O5 star lies along the axis of the SiO velocity gradient, it does not reside midway between the red and blue peaks, nor does it reside at the center of the radio shell. Such an asymmetry of the powering star’s location with respect to the radio shell was predicted by Ball et al. (1992) in their attempt to reconcile the mid-infrared and radio morphologies. Their model predicts the presence of dense molecular gas close to the northeast edge of the UCHII. While single-dish mm/submm surveys have established that G5.89-0.39 is a molecular line-rich source in the 875 m window (Hatchell et al., 1998), the majority of the lines are from sulfur-bearing species rather than the heavy organic molecules more typically seen in hot cores (Thompson & MacDonald, 1999). Clearly, higher angular resolution studies of molecular gas and dust are required to test the various hypotheses on the nature of this complex massive star formation region and attempt to unify the observed phenomena.

Toward this end, we have obtained the first subarcsecond submillimeter (875 m) observations of G5.89-0.39. Our Submillimeter Array (SMA) images have a factor of seven higher angular resolution (in beam area) than the 1.4 mm Sollins et al. (2004) data, and are substantially more sensitive to dust emission and molecular lines. We interpret our unprecedented submillimeter images in the context of new images with comparable resolution constructed from the best available centimeter wavelength data in the VLA archive, 3 mm data from the new B-configuration of the Combined Array for Millimeter Astronomy (CARMA), and archival near- to mid-infrared data.

The distance to G5.89-0.39 used in the literature varies, although all papers agree that the source does not lie at the far kinematic distance. Acord et al. (1998) report a detection of the expansion of the nebula (though with a fairly low signal-to-noise of mas yr) in multiepoch VLA images, which implies a distance of 2.0 kpc. The kinematic distance based on the main-line OH maser velocity of +18 km s is  kpc (Fish et al., 2003). The kinematic distance based on the +9.3 km s LSR velocity of CS(2-1) (Bronfman et al., 1996) is 2.6 kpc, while the distance to the nearby W28 supernova remnant is 1.9 kpc (Velázquez et al., 2002). We will adopt a distance of 2.0 kpc in this paper.

2 Observations

2.1 Submillimeter Array (SMA)

The SMA6 observations were made with eight antennas in the extended configuration on 2006 June 05. Sky conditions were clear and the 225 GHz opacity measured by the tipping radiometer at the Caltech Submillimeter Observatory (CSO) varied from 0.065 to 0.081 during the observations. The SMA receivers have double sideband mixers with 2 GHz bandwidth centered at an intermediate frequency of 5 GHz. The front-end local oscillator was tuned to place the SiO(8-7) spectral line ( GHz) in the center of correlator chunk 23 in the upper sideband (USB), resulting in center frequencies of 336.46 GHz in lower sideband (LSB) and 346.46 GHz in USB. The data were taken with a channel width of 0.81 MHz (0.71 km s), but were subsequently Hanning smoothed for a final spectral resolution of 1.4 km s. The phase center was , (J2000), and the projected baseline lengths ranged from 20 to 230 k. Typical system temperatures were 200 to 300 K at transit and the total time on-source was 4 hours. The gain calibrators were NRAO530 (J1733-130) (13° distant) and J1911-210 (17° distant) and the bandpass calibrator was J2232+117. Flux calibration is based on observations of Ceres and a model of its brightness distribution using the Miriad task smaflux. Comparison of the derived flux of the quasars with SMA flux monitoring suggests the absolute flux calibration is good to within 15%. The estimated accuracy of the absolute coordinates is ″. The uncertainty in the relative positions is better than ″.

The data were calibrated in Miriad, then exported to AIPS where the line and continuum emission were separated with the task UVLSF (using only line-free channels to estimate the continuum). Self-calibration was performed on the continuum data, and solutions were transferred to the line data. The continuum and line data were imaged using natural weighting. All channels were cleaned to a flux density limit of 350 mJy beam, and individual data cubes were then extracted for each spectral line detected. After combining the calibrated LSB and USB continuum uv-data, the 1 rms noise level achieved in the continuum image is 6.4 mJy beam. The noise level in a single channel of the spectral line images is 110 mJy beam. The line data have been corrected for the half channel error in SMA velocity labeling discovered in November 2007. The synthesized beam is  at position angle -32°, corresponding to a linear scale of 1700 AU. The primary beam is ; these data are insensitive to smooth structures larger than about 10″.

We also recalibrated and imaged the 1.4 mm data from the SMA compact configuration originally published by Sollins et al. (2004) (the details of the observing setup can be found in that paper). Using natural weighting, we achieved a continuum rms of 24 mJy beam and a spectral line rms of 150 mJy beam (per 1.1 km s channel), with a beamsize of (P.A.=5.5°). The imaged spectral lines include SiO (5-4), HS (), and HCN (25-24).

2.2 James Clerk Maxwell Telescope (JCMT)

We have obtained newly-processed 850 m SCUBA imaging data of G5.89-0.39 observed at the JCMT7 on Mauna Kea, Hawaii. The data were obtained from the Canadian Astronomy Data Center (CADC) repository of SCUBA Legacy Fundamental Object Catalogue.8 The image of the source was constructed from a number of independent observations and has an effective beamsize of . Further information on the data processing is described by Di Francesco et al. (2008). The total flux density of the source was determined from the image to be 36 Jy. For comparison, the flux density reported in the table of the SCUBA Legacy Catalogue is also 36 Jy (see also Sandell, 1994). Due to the nature of the chopped observations, these values should be considered lower limits. In order to be conservative, we have added 36 Jy as the zero spacing flux density to the SMA uv data. With this addition, the total flux density recovered in the high resolution image is 9.6 Jy, or 27%.

2.3 Very Large Array (VLA)

Several archival datasets from the NRAO9 Very Large Array (VLA) were calibrated and imaged in AIPS. The observing parameters of these data are summarized in Table 1. Some of the datasets were observed in continuum mode, while others were observed in line mode. In project AZ075, the NH (2,2) and (3,3) transitions were observed simultaneously, each with 6.3 MHz bandwidth (82 km s) with 0.195 MHz channel width (2.45 km s).

2.4 Combined Array for Millimeter Astronomy (CARMA)

Our CARMA10 observations were obtained on 12 December 2007 in the B-configuration with 13 antennas (five 10.4m antennas and eight 6m antennas, for a total of 78 baselines). The 225 GHz zenith opacity ranged from 0.19 to 0.23 during the observations. Typical (SSB) system temperatures were 200 to 400 K and the total time on-source was 1.8 hours. The receivers were tuned to an LO frequency of 92.85 GHz. Three continuum bands of 500 MHz width, each comprised of 15 channels, were correlated in each sideband, centered at 89.6, 90.7, and 91.2  GHz (LSB) and 94.5, 95.0, and 96.1  GHz (USB). Only the bands centered at 90.7, 91.2, 94.5, and 95.0  GHz were used due to the significantly reduced sensitivity on the 6m antennas at IFs above 2.4 GHz. The phase center was , (J2000). The gain calibrator was NRAO530 (J1733-130) (13° distant), the bandpass calibrator was 3C454.3, and the flux calibrator was Uranus. The data were calibrated in Miriad and then self-calibrated in AIPS. The combination of the LSB and USB data imaged with a robust weighting factor of -0.5 yielded a beamsize of at position angle 38°, and an image rms of 1.6 mJy beam.

2.5 Infrared images

Several archival images were obtained from the ESO archive, as listed in Table 2. In particular, the data include those published in Feldt et al. (2003) from which those authors identified the proposed exciting source which we call ”Feldt’s star”. The near-infrared (m) images were processed with standard techniques: each raw frame was reprojected to a common projection, and the raw frames combined with a median filter to remove outliers and cosmic rays. No sky subtraction was attempted, since we use these images to identify infrared sources and morphology, but not for photometric measurements. Although the relative pointing of raw frames was very accurate, the absolute astrometry of the data were incorrect by up to 10″. We registered the relatively large field-of-view ISAAC image to 2MASS; since the latter has astrometric uncertainty of less than 0.1″, this registration is limited by centroiding the 10 stars used, with a resulting uncertainty on the ISAAC frame position of 0.1″. We then registered the smaller field-of-view NACO frames to the ISAAC image using 5–10 stars in each frame, with a total resulting positional uncertainty of 0.2″. Our positions differ from those of Feldt et al. (2003) by 0.2″ (our position is nearly directly east of theirs). However, since the offset we find is equal to the uncertainty, for this paper we have chosen to align our near-IR images so that the position of Feldt’s star matches the value published by those authors (18:00:30.44, -24:04:00.9 (J2000)).

The mid-infrared TIMMI2 data required somewhat more detailed processing to remove slowly-changing flat field and detector effects over the course of the observation11. We constructed a ”rolling sky” sequence of frames by taking the median of each four frames in chronological order. As the raw frames are chopped and the telescope nodded, each frame consists of a positive sky image minus a shifted negative sky image, with large-scale thermal emission removed. The position of the source changes according to the telescope nod position, so the median of four consecutive frames contains no information from the sky, but information about time-dependent detector variations. The rolling sky is subtracted from each raw frame, those are masked, shifted, and added to form a mosaic image. Finally, that image had to be registered, which was done using the astrometry of the mid-infrared images in Feldt et al. (1999) since we could not find the Brooks image (Table 2) in any publication.

3 Results

3.1 Continuum images

The subarcsecond resolution SMA 875 m continuum image is shown in Figure 1a. For comparison, we overlay two contours from the lower resolution 1.4 mm image. The north-south extension first noted by Sollins et al. (2004) is now resolved into discrete sources at 875 m; these regions have been labeled for later reference. The submillimeter emission shows a shell-like morphology with an obvious hole in the center. The brightness of the shell varies around its circumference with three major peaks. The easternmost peak lies east of Feldt’s star. Some emission (SMA-N) extends off the north-northeast side of the shell, and a point source (SMA-S) is seen completely separate from the shell to the south-southwest. The point of lowest brightness in the UCHII shell is on the southeast side and may indicate an opening in the shell. Further to the south-southeast, we find faint filamentary emission that appears to trace a coherent structure. No 875 m continuum emission was found outside the field shown in Figure 1a.

Our 2 cm image is shown in Fig 1b. This is the highest resolution and sensitivity image to date of G5.89-0.39 at radio wavelengths shorter than 3.6 cm. As seen in the original 2 cm survey image of Wood & Churchwell (1989), the dominant feature is a shell-like structure surrounding a central cavity of diameter ″. The elongated low-level emission surrounding the shell evident along a position angle of ° has been known since the original 6 cm VLA A-configuration image of Zijlstra & Pottasch (1988). In our new image, we find a nearly-closed loop of emission extending along this direction on the northwest side of the shell (also see Acord et al., 1998). Like the 875 m image, the 2 cm image shows a break in the SE edge of the shell, suggesting a ”blow-out” of ionized gas in this direction. We note that the southern portion of the elongated ionized emission is not exactly coincident with the extension seen at 875 m, but instead lies to its SW.

For comparison, the SMA 875 m and VLA 2 cm images are shown alongside a wide range of other centimeter to near-infrared wavelength images in Figure 2. All of the radio through submillimeter images have been generated from data with a comparable range of uv spacings, have been restored with the same beam ( at P.A.=+35°), and are displayed with the same flux density color scale. The beam was chosen to match the beam of the poorest resolution dataset (3 mm). At the lowest frequency (6 cm), the central part of the UCHII region shows a roughly uniform appearance. As one moves to higher frequency, the source brightens, and the shell structure emerges and becomes increasingly distinct. This change in appearance is due to the decreasing optical depth of free-free emission as the observed frequency increases.

To further demonstrate the effect of free-free optical depth, we measured the flux density in a single beam at various positions in the UCHII region (using the images in Figure 2) as well as the total emission in the field within an 12″ box (see Table 3). The corresponding spectral energy distributions (SEDs) are plotted in Figure 3. At centimeter to millimeter wavelengths, the data are consistent with free-free emission with a turnover at about 2 cm. The excess emission at 875 m toward SMA1 and SMA2 suggests the presence of dust toward these positions. In order to quantify this component, we must estimate and remove the free-free contribution to the 875 m emission. The common analytic approximations for free-free emission given by Altenhoff et al. (1960) and Scheuer (1960) become less accurate above 100 GHz. Thus, we have used the numerical calculation of the Gaunt factor for free-free radiation along with Equations 2 and 15 given by Beckert et al. (2000). Using the non-linear fitting utility of Matlab (lsqcurvefit), we produced the SED fits shown in Figure 3. The resulting electron temperatures and emission measures range from 9000 K and  cm pc at the central position to 12000 K and  cm pc close to Feldt’s star. For the integrated emission, the corresponding values are 8500 K and  cm pc. Assuming the source has equal extent in the third dimension as it does in the plane of the sky (4″ or 0.039 pc) this average emission measure corresponds to an average electron density of  cm. However, this estimate does not take into account the shell-like structure of the UCHII, and in fact the width of the shell is only marginally resolved with resolution, implying that significantly higher are present in the shell.

Based on these results, we followed two methods for deriving a free-free emission map at 875 m and compared their results. The first and simplest method consisted of choosing a characteristic temperature (10000 K) and emission measure ( cm pc), computing the expected flux density scale factor between 1.3 cm and 875 m, and applying it uniformly to all pixels in the 1.3 cm image. The scale factor (0.715) corresponds to an effective spectral index of -0.154 from 1.3 cm to 875 m12. This spectral index is in good agreement with the values used in recent attempts to quantify the Galactic foreground contribution to CMB studies (see for example Dickinson et al., 2003). The extrapolated free-free image was then subtracted from the observed 875 m image in order to generate an image of dust emission, shown in Figure 4. Due to the subtraction, the image rms is increased by to 10 mJy beam, thus it has been masked at the level (50 mJy beam).

For the second method, we first computed the ratio of the 1.3 cm and 3 mm images, then solved for the emission measure as a function of pixel that would theoretically produce the observed ratio, assuming isothermal 10000 K gas. The result is an emission measure image (shown in Figure 4), which we then used to generate a synthetic free-free image at 875 m, which we finally removed from the observed 875 m image. The resulting dust image closely matched that generated by the simpler method, as seen in the “model difference” image in Figure 4. The largest differences are less than when compared with the rms in the observed 875 m image.

Five discrete areas of dust emission are evident, and their positions and flux densities are tabulated in Table 4. SMA1 and SMA2 are compact sources which lie on or near the UCHII shell, SMA-S is compact and lies outside of the shell, and SMA-N and SMA-E are more extended and also lie outside of the shell. As can be seen in the VLA images (Figure 2), there is no corresponding centimeter wavelength emission toward the positions denoted SMA-N, SMA-S, and SMA-E; thus we conclude that this emission is from dust. We note that the locations of SMA-N and SMA-S lie along an axis defined by the southern cluster of OH masers (Stark et al., 2007) and passing through the center of the UCHII shell. Further south, the filamentary structure of 875 m emission must also arise from dust as there is no evidence for it in the centimeter images, which have more than adequate sensitivity to have detected it (if it was free-free emission). We note that there is no evidence for dust emission near the center of the UCHII shell. The 875 m flux density from dust at the position of Feldt’s star is 10 mJy beam which is equal to the noise level. For comparison, if we use our astrometry for Feldt’s star, the flux density from dust would be 31 mJy beam.

3.2 Infrared images

Images at 1.7 m, m ( band) and m ( band) are shown in the bottom row of Figure 2; Feldt’s star is indicated by the open star symbol. The cometary shape of the mid-infrared emission toward G5.89-0.39, coinciding with the northern half of the UCHII shell, has been noted previously (Ball et al., 1992). However, the band image also shows evidence for two distinct peaks within the cometary emission. The position of the brighter band peak is coincident with that of SMA1 (i.e. within the combined astrometric uncertainties). Similarly, the fainter peak coincides with Feldt’s star. In agreement with Puga et al. (2006), we find that most of the band emission in the field shown in Figure 2 is continuum rather than H, with the exception of the sources A, D1, and D2 denoted by Puga et al. (2006). We also examined Spitzer IRAC data for this region but with its poorer angular resolution (″) the mid-IR emission appears as a featureless blob. Comparison of the narrow band 1.2 m image (not shown, see Table 2) with the 1.7 m image shows that all four of the bright stars in the upper half of the 1.7 m image shown in Figure 2 are also detected at 1.2 m, with the star nearest SMA1 being the faintest. It is likely that all four stars are foreground to the nebula.

3.3 Spectral line images

Over three dozen spectral lines were detected in our SMA observations (see Table 5). Eighteen of them were detected with a high signal-to-noise ratio and suffer minimal confusion from line blending. To inspect the spatial distribution of emission in these lines, we created total intensity moment zero images using the total observed velocity extent of the emission, determined for each line independently. Contour maps of these images are shown in Figures 5 and 6 superposed on the 875 m continuum image. In Figure 6, the images of two vibrational lines of SO have been averaged before displaying the contour map. In addition, we detect three faint vibrationally-excited lines of HCN which we have similarly averaged in order to create the contour map in Figure 5.

Figure 5 shows contour maps of a mixture of species, including high-density tracers, high column-density tracers and outflow tracers. As with the dust, none of these gas species show emission inside the UCHII shell, suggesting a cavity at this location. Only CO, CO, CHOH, CS, and HCO exhibit emission significantly (4″) outside of the UCHII ring. Particularly in the southern part of the images, these molecules appear to be tracing the filamentary structures seen in dust emission (§ 3.1). In some cases (CHOH and CS), the molecules follow the inner edge of this structure rather than the structure itself. Faint emission from HCO and CO is also found along this structure. In contrast, the near-IR H features A, D1, and D2 (Puga et al., 2006) appear along the outer boundary of this structure. The strongest emission from HCN, CO, and HCO traces a ridge along the northern and western edges of the UCHII. The latter molecule peaks a fraction of an arcsecond further outward and is the best tracer of the boundary of the mid-infrared emission. In addition, these three molecules show an extension toward the north-northeast which appears to be associated with the dust source SMA-N. Other molecules also appear in this vicinity, including CHOH (also seen at this location by Sollins et al. (2004)), SiO, and CS. This emission is fairly broad ( km s). An additional point source of gas emission very close to the location of the dust source SMA-E is seen in HCO, HCN, CHOH, and CO. Compact peaks of CO and HCN are also seen at the position of SMA2.

Figure 6 shows integrated intensity contour maps of emission in nine lines of SO, SO and their isotopologues in our SMA bandpass. On the east and west sides of the UCHII ring, the peak emission from the most abundant species appears to trace ridges that form a border outside of the free-free emission. The position angles of these ridges are -24° and -26° on the east and west sides, respectively. In contrast, on the north side, the line peaks coincide with the ridge of free-free emission. Toward the southeast side, there is no line emission, further suggesting an incomplete shell in this direction. The less abundant isotopologues show peaks at the positions of SMA1 and SMA2, as do the vibrationally-excited lines of SO and HCN (Fig. 5), indicating high column densities of warm gas. The vibrationally-excited lines also show an isolated peak 1″ east-northeast of Feldt’s star. The spectra toward SMA-N are significantly broader than toward SMA1 or SMA2 and the line wings in HCN and CS extend to redder velocities (see Figure 7). This is the same velocity shift detected in lower angular resolution data, including the OVRO observations of CHCN (Akeson & Carlstrom, 1996) and the SMA observations of SiO (5-4) (Sollins et al., 2004).

The molecular gas shows complicated velocity structure in many of the species. Figure 8 shows first moment maps of several species. Perhaps best seen in SO, the primary feature is a change from redshifted emission on the northeast side of the UCHII to the LSR velocity ( km s) emission on the northwest and southwest sides. The magnitude of the shift is approximately 4 km s. The gas located in the filamentary dust structure near the southern edge of the images lies at the LSR velocity. A blueshifted component is seen in SO at the southern edge of the UCHII. The CHOH and CS show redshifted emission toward SMA-N (see Fig 7). Redshifted emission is also seen in SiO (8-7) near SMA-N, while blueshifted emission is seen in the vicinity of SMA-S. At lower angular resolution these two spatial components were interpreted as a bipolar outflow by Sollins et al. (2004). Although these areas of SiO emission lie along an axis including Feldt’s star, the emission does not extend back to it. Having seen these results, we reanalyzed the lower frequency SMA data from Sollins et al. (2004) and present channel maps of three lines in Figure 9. Clearly, the SiO (5-4) emission is more complicated than a single spatial velocity gradient at position angle +28°. In fact, the emission in all three of these lines demonstrate the complexity of the molecular emission within  km s of the systemic velocity. One is thus driven to studying higher velocity gas in order to seek the origin of bipolar outflow emission in this region.

The high abundance of the CO molecule provides the best chance to study the high velocity molecular gas. Contour maps showing the kinematics of CO (3-2) are shown in Figure 10. In the midst of the complexity of this emission, we identify a linear bipolar structure at high velocities at position angle -4°. The axis defined by this structure passes through the dust core SMA1 and intercepts highly redshifted CO emission north of the UCHII along with a Class I methanol maser (component 2 of Kurtz et al. (2004)). Two redshifted clumps of CO emission correspond in position to the redshifted knots C1 and C2 of the near-IR H emission (Puga et al., 2006). A highly redshifted (+78 km s) water maser feature (Hofner & Churchwell, 1996) also lies along this axis near the far northern extent of the CO emission. Toward the south, the blueshifted knots of the near-IR H emission (component A) lie near the peak of the blueshifted CO. Water masers are also found along the southern portion of the bipolar axis, including a moderately blue-shifted component at  km s at the southernmost maser spot. A higher velocity blueshifted component was once seen at -61 km s in single-dish observations by Genzel & Downes (1977) but has not been detected in subsequent VLA observations, possibly due to variability. Most of the known OH masers, specifically the collections of spots denoted G5.89 Center and G5.89 South (Stark et al., 2007), are clustered in close proximity to the southern axis of this north-south bipolar structure and predominantly exhibit the appropriate highly blueshifted velocities (30-40 km s from the LSR). A notable exception is the southernmost group of masers in G5.89 South which occur the southern edge of the UCHII region and have a small redshifted velocity (from 11-15 km s). At the same position and velocity, we see compact emission in the lowest energy SO line, which appears as a peculiar red feature in the first moment map (Figure 8). The other collection of OH masers, G5.89 East, lies furthest from the bipolar axis and most of its members have velocities near the LSR, suggesting that it is unrelated to the north/south structure.

To further explore the north/south bipolar structure, we reimaged and analyzed the data from the Berkeley-Illinois-Maryland Array (BIMA) published by Watson et al. (2007). Of the four molecular transitions observed, the HCO (1-0) line had the widest velocity coverage (131 km s). With a resolution of , we find the same north-south velocity gradient in HCO emission as is seen in our CO (3-2) data and the maser data. Contour maps of the red and blueshifted BIMA HCO emission are also shown in Fig 10. The sense and direction of the velocity gradient in CO (3-2) and HCO match that seen in a CS (3-2) map obtained with a 17″ beam at the IRAM 30m telescope (Cesaroni et al., 1991). Finally, in Figure 10 we have also plotted the highest velocity emission from the SMA SiO (5-4) data, extending beyond the range shown in the channel maps in Figure 9. Again, the velocity gradient follows a position angle close to the other lines.

3.4 Nh (3,3) maser emission

The NH (3,3) transition shows a complex arrangement of weak emission outside the UCHII region and absorption toward the shell, consistent with the lower resolution data presented by Gomez et al. (1991) and Wood (1993). However, we have discovered an intense point source of emission located at 18:00:29.539, -24:03:52.65 (J2000), which is ″ northwest of the UCHII region (see Fig 11). The fitted line peak is at km s, but the emission is not spectrally resolved with the 2.45 km s velocity resolution of these data (the emission is spread over two channels). After correction for primary beam attenuation and using uniform weighting to construct the image, the peak flux density is  mJy beam, corresponding to a lower limit for the peak brightness temperature of  K at a resolution of . This temperature is significantly higher than expected for thermal excitation and indicates weak population inversion. Other examples of objects where this NH transition is found to be inverted include the high-mass star-forming regions DR21(OH) (Mangum & Wootten, 1994), NGC 6334 (Kraemer & Jackson, 1995; Beuther et al., 2007), W51 (Zhang & Ho, 1995), and IRAS 20126+4104 (Zhang et al., 1999). In all cases, these masers are located in or at the ends of high velocity protostellar outflows. In G5.89-0.39, the peak position and velocity of the NH maser are coincident (to within and 1.0 km s) with a Class I CHOH maser observed with the VLA (component 1 of Kurtz et al. (2004)), which is another outflow tracer. G5.89-0.39 is only the second object in which these two maser transitions have been found to coincide spatially, the other example being DR21(OH) (Mangum & Wootten, 1994). Interestingly, the NH maser is also located within 3″ (north-northeast) of the near-IR knot B of Puga et al. (2006). As shown in Figure 11, these objects lie close to the axis (position angle = -54°) of the Br outflow identified by Puga et al. (2006) as do the near-IR knots D1 and D2. Figure 11 summarizes these new results within the context of the CO (1-0) outflow map from Watson et al. (2007).

4 Discussion

In order to interpret our unprecedented subarcsecond submillimeter images, we have assembled a comprehensive dataset on this object. The multiwavelength nature of our work warrants a fresh examination of the morphology of the UCHII region, the nature of the ionizing source, the number and structure of the outflow(s), and the presence of young stellar objects in this important region.

4.1 Shape and structure of the ionized nebula

UCHII regions are divided into a number of morphological classes (Wood & Churchwell, 1989), with one of the rarer classes being “shell-type”. G5.89-0.39 is perhaps the most distinctive and well known member of this class with a well-defined central cavity with a radius approximately one-third that of the outer UCHII radius. Our new subarcsecond submillimeter observations have revealed that this cavity is devoid of dust, as would be expected when a strong ionizing source is present within the shell (Churchwell et al., 1990; Faison et al., 1998). A striking result from our molecular line images is that, at subarcsecond scales, no species are seen on the line of sight to Feldt’s star, indicating a low column density of obscuring material close to the star. This result is a vivid demonstration of why this star can be seen at near-IR wavelengths. In the discovery paper, it was noted that Feldt’s star is significantly offset from the center of the cavity, and the possibility of proper motion as an explanation was discussed. At present, there is no knowledge of the velocity of the star with respect to the UCHII region or the molecular gas. Another explanation that may apply to G5.89-0.39 comes from recent modeling of the dynamical expansion of UCHII regions in self-gravitating molecular clouds, which demonstrate that nearly spherical shells can be produced even if the ionizing star is off center (Mac Low et al., 2007). In any case, our data provide no counter evidence to the hypothesis that Feldt’s star is the ionizing source of the UCHII region.

Now we turn our attention to the cm emission just outside of the shell. Ever since the original VLA A-configuration 6 cm image (Zijlstra & Pottasch, 1988), the elongation of the low-level cm halo of G5.89-0.39 has been evident along a position angle of °. In our SMA spectral line images, we now see that the distribution of the molecular gas traced by SO follows this same position angle as it effectively encompasses and bounds the UCHII region on two sides (see the upper left panel of Fig. 6). The morphology of this line emission suggests that the molecular gas is being driven outward by the expansion of the UCHII region. The velocity shifts seen in the first moment maps (Fig 8) suggest that 2 km s is a lower limit to the expansion rate (because an additional component of motion could be present transverse to the line of sight). The fact that the molecular gas on the northwestern side of the UCHII peaks on the shell rather than outside of it, in contrast to the northeast and southwest sides, suggests that molecular gas is missing from the northwestern “cap”. This arrangement would explain the position angle of the northern elongation of the cm halo. The possibility that the ionization front is also “blowing out” in the southeastern direction was previously suggested by Ball et al. (1992). Our observations support this suggestion in two respects: we see a break in the 875 m continuum emission at the point where this southeastern axis intersects the UCHII shell, and we see extended source of dust emission (SMA-E) located radially outward from this break. The fact that nearly all of the dust emission lies within a few arcseconds of the outer edge of the UCHII region is a result that agrees with predictions from previous infrared observations (Harvey et al., 1994).

4.2 Location of Molecules

As seen in single dish spectra (Thompson & MacDonald, 1999), we find G5.89-0.39 to be almost completely lacking in organic molecular line emission while being rich in sulfur-bearing molecules. This result is in sharp contrast to observations of massive young stellar objects (with hot cores) with the same spectral setup and angular resolution with the SMA (see e.g. Brogan et al., 2007a). Based on JCMT line profiles of G5.89-0.39, Thompson & MacDonald (1999) suggest that HCN traces the envelope surrounding the UCHII while the sulfur and silicon-bearing species with broad lines trace an outflow. Our images confirm that HCN emission is located primarily along the northwestern half of the UCHII shell, and that much of the SiO and CS is found at larger distances from the UCHII. However, as described in sections 3.3 and 4.1, we find that SO and SO are also concentrated primarily around the shell, apparently constraining the expansion of the UCHII region. They do not show any obvious evidence of tracing the beginnings of the larger scale east/west outflow reported in CO (Klaassen et al., 2006; Watson et al., 2007). Instead, the optically-thin isotopologues of SO and SO appear to be tracing compact regions of emission at the LSR velocity, similar to CO and HCN. At least two of these compact regions of line emission are coincident with the dust cores SMA1 and SMA2, while a third may be associated with the Br outflow source described by Puga et al. (2006) (see Figures 5 and 6).

East of SMA-S, there is an arc-like structure seen in CO, CHOH, and CS just north of a similar filamentary structure seen in dust emission. Moreover, emission from the high column density tracers CO and HCO coincides with the dust emission. It is unclear whether this morphology is associated with the expansion of the UCHII region or is an independent structure. Remarkably, SMA-S itself is nearly free of line emission in our data. We find this to be an enigmatic source, similar to NGC6334I SMA4 (Hunter et al., 2006). It could be a very young protostellar object in which most common gaseous species are frozen onto grains. Following the method described in § 4.3, we estimate the total mass of SMA-S to be 12-20 M assuming the temperature is in the range of 15-20 K. This amount of mass is sufficient to be an intermediate mass pre-protostellar core. Future high-angular resolution observations in other species that remain in the cold gas phase longer in the evolutionary sequence, such as NH and HD (Flower et al., 2006), would be useful in studying this object.

4.3 The Nature of SMA1 and SMA2

There are a growing number of examples of chemical complexity on arcsecond scales in regions of massive star formation. One reason for this complexity is the superposition of multiple young stellar objects at slightly different LSR velocities, as in, for example Cep-A East, (Brogan et al., 2007b; Comito et al., 2007). The SMA images of G5.89-0.39 have a spatial resolution of AU, which is finer than the typical separation of members of protoclusters (Hunter et al., 2006). Thus we are able to resolve the emission from individual objects and compare their line strengths in various chemical species. The set of spectral line profiles in Figure 7 shows that the ratio of the strength of HCN to SO and SO is larger in SMA1 than in SMA2, as is the ratio of CO to SO and SO. The high critical density and excitation temperature of the submillimeter HCN transitions makes them an excellent tracer of warm dense gas, as is seen in many high-mass star forming regions observed at high angular resolution (Cyganowski et al., 2007; Mookerjea et al., 2007; Wright et al., 1996). Also, the strong detection of the CO (3-2) line requires a large column density of moderately warm gas. The relatively greater strength of these two lines in SMA1 combined with its detection at m, and association with the north/south bipolar outflow suggests that SMA1 likely contains a protostar. SMA2, being richer in sulfur-bearing molecules, likely has a greater fraction of its molecular excitation due to shocks, and it is less clear that it harbors a protostar. We have fit the line profiles toward SMA1 and SMA2 with single Gaussians (Table 6) and find that SMA2 is somewhat broader (by 1 km s) and bluer (by 0.5 km s), supporting our conclusions.

Unfortunately, it is difficult to estimate the luminosities of SMA1 and SMA2 due to the lack of arcsecond resolution images between 20m (Feldt et al., 1999) and 875 m. The m detection of SMA1 cannot accurately constrain the peak of the dust spectral energy distribution or the dust temperature. The best hope for further progress in the near future would be to measure the gas temperature. One can imagine constructing a rotation diagram from our dataset using the multiple transitions from SO. However, this molecule is clearly very optically thick, as its isotopologues rival the main line in brightness, and they have somewhat smaller linewidths. Thus, an optical depth correction is essential, but this can only be done using a pair of transitions from the same energy level. Although we do have one such pair (SO  and SO ), the comparison is complicated by the hyperfine structure of SO. Accurate temperature measurements of SMA1 and SMA2 will require a more extensive combination of spectral lines and isotopologues than are found in the present SMA data. In the meantime, a lower limit to the temperature can be obtained from the brightness temperature of the brightest line (in both cases it is an SO line). For SMA1, this value is 34 K, and for SMA2 it is 65 K. When compared to the brightness temperatures of the dust emission in Table 4, these values yield upper limits to the 875 m dust optical depth of 0.066 and 0.10. Using these values along with the dust continuum flux densities, one can derive an upper limit estimate for the mass. In the case of SMA1 and SMA2, we use the peak flux densities in order to assess the unresolved point source component at these positions. Following Equation 1 of Brogan et al. (2007b), in which  cm g and the gas to dust mass ratio is 100, we find upper limits of 6 M and 2 M for SMA1 and SMA2, respectively. For comparison, if the dust temperature is 100 K, the corresponding masses would be 1.6 M and 1.3 M.

Our observations of CO (3-2) provide an independent estimate of the mass of SMA1 and SMA2. CO has a well-known abundance with respect to H of (Frerking et al., 1982) and has been used in single dish studies of other ultracompact HII regions (Hofner et al., 2000). At the moderately warm temperatures of SMA1 and SMA2, the possible confounding effect of depletion of CO onto grains should be minimal, as the desorption temperature is likely in the range of 15-40 K (Jørgensen et al., 2006; Doty et al., 2004). We have measured the integrated intensity of CO towards these objects, and used the equations of Mangum & Shirley (2006) to determine the column density in the optically-thin, Rayleigh-Jean limit. Using the peak brightness temperatures of the line, we estimate the optical depths to be small as long as the excitation temperature is above  K. In Table 7, we list the opacity, column density and mass (0.9-1.6 M) computed for excitation temperatures of 75 K and 150 K. The column density and mass have been corrected for opacity. The values at 75 K are only about 20% higher than the minimum values one would obtain if the excitation temperature was set equal to the upper state energy (32.3 K). A temperature of 100 K provides good agreement between the dust-derived and CO circumstellar mass, and also provides a greybody model consistent with the submillimeter and mid-IR flux densities at the position of SMA1. Indeed, based on mid-infrared images at 11.7 and 20 m, Feldt et al. (1999) derive a temperature of 120 K for the hot dust component at the position of SMA1. Using a temperature of 100 K, the total luminosity of SMA1 is  L. For reference, on the main sequence this luminosity corresponds to an early B star with a central stellar mass of 7.5 to 8.5 M depending on the mass-luminosity relation used (Demircan & Kahraman, 1991; Hilditch & Bell, 1987).

In any event, the circumstellar masses we obtain for SMA1 and SMA2 are comparable to the upper limits measured for the BN and IRc2 objects in Orion (Eisner & Carpenter, 2006), and to circumstellar masses surrounding intermediate mass protostars (Beltran et al., 2007; Neri et al., 2007). However, we emphasize that the nature of SMA2, including whether or not it contains a central heating source, remains unclear, particularly because it lacks a clear association with a bipolar outflow (unlike SMA1). The locations of SMA1 and SMA2 with respect to the UCHII region shell and their gas masses are broadly consistent with simulations of “secondary collapse” by Mac Low et al. (2007). In the Mac Low et al. (2007) model, gravitational instabilities in the material swept up in an expanding UCHII region shell can produce a second generation of collapsing cores in the shell with masses of a few M. Quantitative predictions for the dust emission from such cores are not made, but the prediction for the free-free emission that would be present at the boundary between the ionized gas and dense core of a few 100 mJy is consistent with the observed free-free emission in the vicinity of SMA1 and SMA2.

The lack of organic ”hot core” emission from SMA1 (and SMA2 if it also contains a protostar) has several possible explanations: (1) the protostar may be of sufficiently late type (low mass) to preclude the formation of a hot core. However, as already discussed the circumstellar mass is consistent with those of intermediate mass protostars, and such sources are capable of producing hot core line emission (see for example Fuente, 2008). (2) The protostar could be sufficiently young that it has not yet reached the hot core phase, but this seems less likely given the presence of the energetic north/south outflow from SMA1. (3) The progenitor of the UCHII region (presumably Feldt’s star) may have melted enough of the icy dust mantles in its vicinity during its formation that the reservoir of organic material was severely depleted for later generations of protostars. (4) Although the simulations of Mac Low et al. (2007) show that the column densities of the “secondary collapse” cores are not disrupted by the passage of the UCHII shock, the effect on the chemistry of such cores has not been assessed. Thus it is possible that the passage of the UCHII shock destroyed the fragile organic molecules. Higher angular resolution study of the dust emissivity and molecular line emission with ALMA in the future will help distinguish between these possibilities.

4.4 Origin of the outflows

What do our SMA data reveal about the massive outflow from G5.89-0.39? First of all, the 875 m line emission shows some evidence of a general expansion centered on the UCHII region. The location of the clump of molecular line emission associated with SMA-N combined with the ridge of dust, CO, and HCO emission in the south suggest that Feldt’s star is more likely than SMA1 to be the origin of this activity. Regarding the large-scale east/west CO outflow, there is no obvious bipolar structure tracing back to Feldt’s star or any of the dust cores along the position angles quoted by Klaassen et al. (2006) and Watson et al. (2007). This result is consistent with the “extinct jet” hypothesis of Klaassen et al. (2006), whereby the large scale outflow is a remnant flow from a previous generation of protostellar activity. However, our CO(3-2) and HCO(1-0) images do show direct evidence for a collimated bipolar outflow originating from SMA1 at position angle -4°, lending further weight to our identification of it as a protostar. This outflow matches the outflow direction postulated by Puga et al. (2006), as it naturally explains the arrangement of the near-IR H knots A and C. In addition, this outflow explains the geometry of most of the maser emission spots in this region, including HO, OH, and Class I CHOH masers. In particular, the locations and kinematics of the high-velocity HO maser components to the north and south of the UCHII region, along with the blueshifted OH masers to the south, follow the CO and HCO velocity field centered on SMA1 (see Fig. 10).

In addition to the outflow from SMA1, we also find tentative evidence for a second outflow in redshifted and blueshifted CO emission and blueshifted HCO emission, centered on the Br outflow origin, with a position angle of -54°, similar to that found by Puga et al. (2006). The near-IR knots to the northwest and southeast (B and D) are aligned roughly with this axis, as is the new NH (3,3) maser. It seems likely that all of these phenomena are outflow-related, and have a common driving source. We find a point source of emission from vibrationally-excited SO and HCN at this position (see Fig. 5 and 6), suggesting a central powering source at this location. While we do not find a compact dust source at the origin, it does lie at the  mJy level within a ridge of dust emission associated with the northwest edge of the UCHII region. We estimate an upper limit of  mJy for a point source at this position, which corresponds to a mass upper limit of  M (assuming a temperature of 65 K), comparable to gas masses found to be surrounding intermediate-mass protostars (Beltran et al., 2007; Neri et al., 2007). Another explanation for knots B and D and the (3,3) maser is that they are simply enhancements in the large-scale east/west outflow. The fact that the NH maser (and its associated CHOH maser) sit at the edge of the large blueshifted CO (1-0) outflow lobe (Figure 11) and emit near the LSR velocity is consistent with them originating from a velocity coherent column of shocked gas moving transverse to the line of sight, which is an ideal geometry for generating strong maser features (see e.g. Liljeström & Gwinn, 2000).

One question that arises from the work presented here is the relative ages of the compact north/south outflow from the protostar SMA1, the large scale east/west outflow, the UCHII region and Feldt’s star. We can estimate the timescale for the north/south outflow from SMA1. Using the velocity (w.r.t. the G5.89 LSR of 9 km s) and projected offset of the northern HO maser relative to SMA1, (69 km s and 0.073 pc, respectively), the outflow timescale is yr ( is the radial velocity of the outflow and is the angle between the line-of-sight and the outflow direction). Using the southern CO (3-2) emission, the timescale is about 1700 years, yielding an average timescale of  yr. For comparison, the timescale for the motion of H knot A is 400 yr (Puga et al., 2006). Although the timescale for the north/south CO outflow is only a lower limit on the outflow age, it is a factor of five discrepant from the 7700 yr age reported by Watson et al. (2007) based on BIMA CO (1-0) data. However, the velocity ranges used by Watson et al. (2007) in computing the age do not include the high velocity gas that is present in the north/south outflow and are more appropriate to the lower velocity east/west outflow. The timescale we derive for the north/south outflow is closer to the value of 2000 yr reported by Klaassen et al. (2006). This agreement makes sense, because the highest-velocity CO (3-2) emission is unresolved in their single-dish maps, since those spectral channels are completely dominated by the compact outflow from SMA1 rather than the larger, older east/west outflow.

The question remains as to the driving source of the east/west flow. The width of the east/west outflow combined with the confused velocity field within 15 km s of the LSR prevents any constraint on the outflow origin to better than a few arcseconds. As an additional source of uncertainty, if the driving source has a small relative motion of 1 km s, it could have moved during 7700 yr. Two possible candidates are Feldt’s star and SMA2. Lacking direct evidence for an internal heating source in SMA2, we find the most natural scenario for the east/west outflow is that it was driven by Feldt’s star prior to the creation of the UCHII region. The expansion of the UCHII region has disrupted the velocity field around it making it impossible to trace this outflow back to its origin. For the same reason, there could be additional bipolar outflows present (e.g. from other embedded YSOs) that are difficult to discern in this complex velocity field.

Finally, a possible clue regarding the relative age of the outflow from SMA1 and the UCHII region may be found in the fact that the southern lobe of the outflow only appears distinctly in the CO (3-2) contour maps in the regions outside the boundary of the UCHII region. One could argue that the outflow from SMA1 predates the UCHII region, which has recently expanded and disrupted the inner portions of the bipolar morphology, leading to the OH masers. Alternatively, if the UCHII region is significantly older than the proposed 600 yr age (Acord et al., 1998) and SMA1 is located behind it, then the southern lobe of the outflow from SMA1 may have recently drilled through the UCHII region, creating the OH masers on the front side, and continuing on south of the UCHII region. In this picture, the absence of high velocity redshifted OH features north of SMA1 would be explained by the high optical depth of the UCHII region at 1.6 GHz. We believe the former picture seems more physically plausible, given the lack of a disturbance in the ionized gas morphology along the outflow direction (i.e. the position angle of the low-level centimeter emission differs from the outflow by 24°). In either case, our observations provide further evidence that the OH masers in G5.89-0.39 are associated with a protostellar outflow, a conclusion originally reached by Zijlstra et al. (1990). Other examples of this correlation have been found, for example the TW object of W3(OH) (Argon et al., 2003). Proper motion measurements of the OH masers in G5.89-0.39 should provide valuable insight on this phenomenon.

5 Conclusions

Our subarcsecond submillimeter images of the ultracompact radio source G5.89-0.39 (W28 A2) have shed new light on this enigmatic source. By using a comprehensive set of lower-frequency images, we have modeled and removed the free-free emission and find five residual sources of dust emission. With no dust emission located inside the shell, our observations support the previously-proposed picture of a dust-free cavity located inside a shell-like UCHII region with warm, high-density gas and dust tracing its periphery. Two of the compact dust objects, SMA1 and SMA2, exhibit compact spatial peaks in one or more of the optically thin tracers SO, SO, and CO. In CO (3-2) emission, we have identified a well-collimated, high-velocity outflow from SMA1 at position angle -4° and conclude that it is an embedded intermediate-mass protostar, surrounded by  M of circumstellar material at a temperature of approximately 100 K, and corresponding to the brightest source in the m image. The outflow from SMA1 explains much of the near-IR H and centimeter wavelength maser emission in the region. We also find tentative evidence for a second CO outflow associated with the Br outflow identified by Puga et al. (2006). The position angle of this outflow points toward the location of a new NH (3,3) maser that we have discovered 12″ northwest of the UCHII region. The origin of this outflow is marked by compact emission in vibrationally-excited transitions of SO and HCN. Regarding the dust source SMA2, it is unclear whether it harbors a central heating source or is the result of a strong shock. Nonetheless, the masses and locations of SMA1 and SMA2 are broadly consistent with the Mac Low et al. (2007) model of secondary collapse in material swept up in the expanding shells of UCHII regions. Located outside the UCHII shell, we detect a cold dust object (SMA-S) that is remarkably free of line emission. Assuming a low temperature of 15-20K, we estimate a mass of 12-20 M, which suggests it is likely to be an intermediate mass pre-protostellar core. Beyond SMA-S, we detect the beginning of a filamentary structure of dust and gas emission that extends for several arcseconds eastward before turning northeast. Whether this structure is related in any way to the UCHII region is an open question.

We thank C. Watson for providing access to BIMA data on this region. CJC is supported by a National Science Foundation Graduate Research Fellowship. RI was supported in part during this research by a NASA Spitzer fellowship to UVa. This research has made use of NASA’s Astrophysics Data System Bibliographic Services and the SIMBAD database operated at CDS, Strasbourg, France, and is based in part on observations collected at the European Southern Observatory, Paranal, Chile. This research used the facilities of the Canadian Astronomy Data Centre operated by the National Research Council of Canada with the support of the Canadian Space Agency. TRH thanks K. T. Constantikes for introducing him to the power of Matlab.
Figure 1: (a) Greyscale and thin contours: SMA continuum image at 875 m with a resolution of  at position angle -32°. Contour levels are -20, 20, 39, 65, 130, 195, 325, 455, 585 mJy beam. Thick contours: SMA continuum image at 1.4 mm with a resolution of at position angle +5°. Contour levels are -50, 50, 250 mJy beam. (b) Greyscale and thin contours: VLA 2 cm image with a resolution of  at position angle -1°. Contour levels are -0.9, 0.9, 3, 7.5, 15, 30, 60, 90, 150, 210 mJy beam. Thick contours: SMA continuum image at 1.4 mm. The respective beams are shown in the lower left corners. The star symbol marks the location of Feldt’s star (Feldt et al., 2003). The large crosses denote regions of excess 875 m emission above that expected from free-free emission extrapolated from cm wavelengths (see Figure 4 and Table 4).
Figure 2: Continuum images of G5.89-0.39 from low frequency in the upper left to high frequency in the lower right. The top six images have been restored with the beamsize of the 3 mm image ( at P.A.=+35°), and are shown on the same colorscale (0.0 to 1.1 Jy beam). The white crosses mark the positions of the submillimeter dust sources identified in Figure 4 and Table 4. The star symbol marks the position of Feldt’s star (Feldt et al., 2003).
Figure 3: Spectral energy distribution at four positions. Lower black line and points: SMA1; red line and points: SMA2; purple line and points: center of shell; blue line and points: east of Feldt’s star (see Table 3). The upper black line and points are the flux densities integrated over a 12″ square box. The solid lines are models constructed from free-free emission (Beckert et al., 2000), plus (in the case of SMA1, SMA2 and integrated) dust emission with a spectrum, shown as a dash-dot line. The upper open squares are the SCUBA peak flux density and the integrated flux densities from the 1.4 mm SMA data and the 3 mm CARMA data. The error bars reflect the uncertainty values in Table 3.
Figure 4: Top left panel: Observed SMA continuum image at 875 m restored with a resolution of  at position angle +35°. The beam is shown in the lower left corner. Contour levels are -0.03, 0.06, 0.12, 0.24, 0.48, 0.96 Jy beam. Top right panel: Estimated dust contribution to the observed flux density, computed by scaling the 1.3 cm image by 0.715, and removing it from the 875 m image. The 875 m image was masked at the 50 mJy level prior to subtraction. Bottom left panel: Model of the emission measure vs. position computed by comparing the flux ratio between the 1.3 cm and 3 mm images on a pixel-by-pixel basis. At 1.3 cm, the model has a mean optical depth of 0.115, and a maximum optical depth of 0.268. Bottom right panel: The difference between the dust image shown in the top right panel and the dust image produced by using the free-free model shown in the bottom left panel. In all panels, the star symbol marks the location of Feldt’s star, and the black crosses denote the positions of sources of dust emission (see Table 4). The triangles mark near-IR H emission components (A, D1, D2), and the yellow asterisk is the Br outflow origin, all from Puga et al. (2006).
Figure 5: Integrated intensity contour maps of submillimeter spectral lines overlaid on the SMA 875 m continuum image. The black crosses mark the positions of the dust sources SMA1, SMA2, SMA-N, SMA-E, and SMA-S (see Table 4). The yellow triangles mark the near-IR H knots (A, D1, D2), and the yellow asterisk marks the origin of the Br outflow (Puga et al., 2006). Contour levels are 3, 4, 5, 6, 7, 8, 10, 12, 14, 16, 18, 20 the base level listed in Table 5. The number below the species name is the energy above ground for the lower level in units of cm. The image for HCN is the average image of the two lines in the table. The image for HC is the average image of the three lines in the table.
Figure 6: Integrated intensity contour maps of submillimeter SO and SO spectral lines overlaid on the SMA 875 m continuum image. The black crosses mark the positions of the dust sources SMA1, SMA2, SMA-N, SMA-E, and SMA-S (see Table 4). The yellow triangles mark the near-IR H knots (A, D1, D2), and the yellow asterisk marks the origin of the Br outflow (Puga et al., 2006). Contour levels are 3, 4, 5, 6, 7, 8, 10, 12, 14, 16, 18, 20 the base level listed in Table 5. The number below the species name is the energy above ground for the lower level in units of cm. The image for SO  is the average image of the two lines in the table.
Figure 7: Spectral line profiles extracted from a single data cube pixel at the positions of the dust continuum sources SMA1, SMA2, and SMA-N. The dotted vertical line marks the LSR velocity of +9 km s. The SO line is , the SO line is , the HCN line is 38-37, and the CHOH line is .
Figure 8: Plotted in color scale are the first moment maps of various transitions from the 875 m SMA dataset. The black contour is from the 3 mm CARMA continuum image (0.1 Jy beam). The crosses mark the dust sources from Table 4 and the asterisk marks the origin of the Br outflow (Puga et al., 2006). The star symbol marks the position of Feldt’s star (Feldt et al., 2003). The SMA beam ( at P.A.=-32°) is drawn at the lower left.
Figure 9: Contour channel maps of three spectral lines from the 1.4 mm SMA dataset: HCN 25-24, SiO 5-4, and H, observed with a beamsize of at P.A.=+5°. Each panel is comprised of eight channels of width 1.1 km s centered at the velocity listed; there are no gaps in velocity coverage between panels. Contour levels are -7, -5, -3, 3, 5, 7, 10, 15, 20, 30, 40, 50, 60 Jy beam km s for HCN, Jy beam km s for SiO, and Jy beam km s for HS. In each panel, the greyscale image is the 217 GHz continuum, and the dotted lines mark the position angles reported (counter-clockwise from north) in SiO (5-4) (Sollins et al., 2004) and CO (1-0) (Watson et al., 2007). The position angle of the single-dish CO (3-2) outflow (Klaassen et al., 2006) is similar to that shown for CO (1-0).
Figure 10: (a) Contour map of CO (3-2) emission observed with a beamsize of with P.A.=-37°. The velocity ranges of integration are +80 to +40 km s (red contours), +40 to +20 km s (orange), -4 to -41 km s (green), -41 to -81 km s (blue). The contour levels are 15, 20, 25, 30, 37.5, 45, 52.5, 60, 67.5 Jy beam km s. (b) Contour map of the high velocity HCO (1-0) emission from the dataset of Watson et al. (2007) with a beamsize of (P.A. = 1.4°); red contours: +30 to +55 km s with levels 4.5, 6.3, 8.1, 9.9 Jy beam km s; blue contours: -20 to -45 km s with levels 2, 2.75, 3.5, 4.25 Jy beam km s. (c) Contour map of the high velocity SiO (5-4) emission from the Sollins et al. (2004) data with a beamsize of at P.A.=+5°, and identical velocity ranges as for HCO. Levels are 4, 7, 10, 14, 18 Jy beam km s. In all panels, the black contours are from the 3 mm continuum image (0.1 and 0.5 Jy beam). Filled squares are OH masers from Stark et al. (2007), open squares are HO masers from Hofner & Churchwell (1996), triangles are the near-IR H knots from Puga et al. (2006), and the pentagon is a Class I CHOH maser position (component 2 from Kurtz et al. (2004)). The dashed line indicates the outflow axis associated with the dust source SMA1 while the dotted line represents the orientation of the Br outflow and the asterisk marks it origin (Puga et al., 2006).
Figure 11: Contour maps of CO (1-0) emission from the dataset of Watson et al. (2007) observed with a beamsize of with P.A.=+1°. The velocity ranges of integration are +20 to +58 km s (dashed contours), and 0 to -37 km s (solid contours). These ranges begin at slightly higher velocities and extend to higher velocities than those shown in Watson et al. (2007). Levels are 5, 10, 20, 40, 80, 160 Jy beam km s. Triangles are the near-IR H knots from Puga et al. (2006), pentagons are Class I CHOH masers from Kurtz et al. (2004), crosses are the 875 m dust sources, and the “X” is the NH (3,3) maser from this work. The dashed line indicates the outflow axis associated with the dust source SMA1 while the dotted line represents the orientation of the Br outflow source (asterisk) (Puga et al., 2006).
Project Date Freq. Wave. sensitivity Beam
Code YYYY-MM-DD (GHz) (cm) Observation Type Config. (mJy beam) ″  ″ (P.A.) Calibrators
AW15813 1986-04-27 5 6.0 continuum A 1.3 (-17°) 3C286, J1748-253
AA227 1998-06-24 8.4 3.6 continuum BnA 0.17  (+73°) 3C286, J1730-130
AW550 2001-05-12 15 2.0 recomb. line B 0.30  (-1°) J2251+158
AZ075 1995-11-17 22 1.3 NH (2,2), (3,3) B 1.814 (+22°) 3C286, J1730-130
Table 1: Summary of VLA archival data
telescope instrument filter/band observation date program PI
VLT-Yepun NACO H 2.122m 2004/06/15 073.C-0178(A) Feldt
VLT-Yepun NACO Br 2.166m 2004/06/15 073.C-0178(A) Feldt
VLT-Yepun NACO L’ 3.8m 2004/06/15 073.C-0178(A) Feldt
VLT-Yepun NACO K 2.1m 2004/06/15 073.C-0178(A) Feldt
NTT SOFI Br 1999/08/01 63.H-0292(A) Kaper
NTT SOFI Br-continuum 1999/08/01 63.H-0292(A) Kaper
NTT SOFI 1.215m 1999/08/01 63.H-0292(A) Kaper
VLT-Antu ISAAC 1.71m 2002/07/09 69.C-0448(A) Bik
3.6m TIMMI2 11.9m 2003/07/20 71.C-0342(A) Brooks
Table 2: ESO Archival Infrared Images
Flux density (Jy beam)15
Position 6 cm 3.6 cm 2 cm 1.3 cm 3 mm 875 m
SMA1 0.205 0.449 0.706 0.863 0.771 1.11
SMA2 0.151 0.375 0.652 0.812 0.740 0.987
center of shell16 0.202 0.386 0.535 0.614 0.528 0.349
east of Feldt’s star17 0.199 0.429 0.673 0.829 0.765 0.632
total (within a 12″ box) 2.36 4.93 7.37 8.81 7.51 6.44
calibration uncertainty 5% 5% 5% 10% 10% 15%
Table 3: Comparitive continuum intensities at various positions in the UCHII region
Peak flux Peak brightness Integrated
Position density18 temperature flux density19
Source () () Jy beam K Jy
SMA1 18 00 30.37 -24 04 00.27 4.2 20
SMA2 18 00 30.30 -24 04 02.37 3.3 21
SMA-N 18 00 30.48 -24 03 58.87 2.6
SMA-S 18 00 30.29 -24 04 05.03 3.2
SMA-E 18 00 30.62 -24 04 03.14 1.6
Table 4: 875 m positions and flux densities of dust sources
Species Transition Frequency Sideband E Base contour
(GHz) (cm) Jy beam km s
CHOH 335.58200 lower 43.7 4.02
SO 336.08923 lower 180.6 3.76
HCN 37-36 336.52008 lower 202.1 2.40
SO 336.55381 lower 88.1 2.87
SO 336.66958 lower 159.1 2.76
CHOH 336.86511 lower 125.7 3.55
CO 3-2 337.0611322 lower 11.2 1.65
SO 337.19862 lower 44.7 3.73
HC 337.34469 lower 425.6 0.90
CS 7-6 337.39646 lower 33.8 4.46
SO 345.51966 upper 32.7 3.30
SO 345.55309 upper 28.2 ..
SO 345.5850023 upper 105.1 2.81
HCN 38-37 345.60901 upper 213.3 ..
SO 345.65129 upper 24.4 ..24
SO 345.67879 upper 21.1 ..25
CO 3-2 345.79599 upper 11.5 41.4
CHOH 345.90397 upper 219.7 ..26
CHOH 345.91919 upper 307.8 ..27
SO 345.92935 upper 112.5 3.26
CHOH 346.20278 upper 68.5 ..28
NS 346.22014 upper 37.8 ..29
HC 346.45573 upper 436.8 ..
SO 346.52388 upper 102.7 ..30
SO 346.52848 upper 43.2 ..31
SO  346.59179 upper 643.6 1.62
SO 346.65217 upper 105.3 5.16
HC 346.94912 upper 437.1 -
HCO 4-3 346.99834 upper 17.4 3.10
SiO 8-7 347.33056 upper 40.6 5.77
Table 5: Spectral lines detected in the 875 m SMA data
Species E Amplitude32 Velocity FWHM Amplitude33 Velocity FWHM
(cm) (Jy beam) (km s) (km s) (Jy beam) (km s) (km s)
SO 32.7 1.38 0.06 9.06 0.11 5.27 0.28 1.69 0.08 8.64 0.13 6.08 0.34
SO 105.1 1.03 0.06 9.32 0.14 4.92 0.35 1.06 0.06 8.82 0.11 4.18 0.27
HCN 213 1.17 0.10 8.95 0.17 3.96 0.42 0.83 0.06 8.08 0.17 4.77 0.45
SO 112.5 1.32 0.06 9.30 0.11 4.83 0.26 1.51 0.06 9.12 0.12 5.80 0.30
SO 105.3 2.06 0.06 7.76 0.09 5.96 0.22 4.33 0.06 8.08 0.06 7.94 0.14
SO 180.6 2.28 0.06 8.44 0.07 5.69 0.17 2.36 0.05 8.13 0.07 7.22 0.17
HCN 202.1 1.06 0.09 9.03 0.17 4.00 0.42 0.90 0.05 8.51 0.10 3.47 0.25
SO 88.1 1.81 0.05 9.43 0.06 4.42 0.15 1.91 0.06 8.78 0.08 5.37 0.20
SO 159.1 1.67 0.05 8.59 0.08 5.38 0.20 2.08 0.07 8.17 0.10 5.90 0.25
CO 11.2 1.05 0.07 9.18 0.13 3.99 0.32 0.59 0.08 8.82 0.29 4.54 0.71
SO 44.7 1.48 0.07 9.76 0.09 4.10 0.21 1.74 0.06 8.66 0.08 4.80 0.19
weighted average 8.94 0.03 4.94 0.07 8.42 0.03 6.07 0.07
Table 6: Results of Gaussian fits of lines toward SMA1 and SMA2
Integrated K K
Source Intensity34 N(CO) H mass35 N(CO) H mass36
Name (Jy km s) (neper) log(cm) (M) (neper) log(cm) (M)
SMA1 0.24 12.79 0.11 12.97
SMA2 0.13 12.70 0.06 12.89
Table 7: Gas masses from CO (3-2)


  1. affiliation: NRAO, 520 Edgemont Rd, Charlottesville, VA, 22903
  2. affiliation: NRAO, 520 Edgemont Rd, Charlottesville, VA, 22903
  3. affiliation: NRAO, 520 Edgemont Rd, Charlottesville, VA, 22903
  4. affiliation: University of Virginia, Astronomy Dept., P.O. Box 3818, Charlottesville, VA, 22903-0818
  5. affiliation: University of Wisconsin, Madison, WI 53706
  6. The Submillimeter Array (SMA) is a collaborative project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute of Astronomy & Astrophysics of Taiwan and is funded by the Smithsonian Institution and the Academia Sinica.
  7. The James Clerk Maxwell Telescope is operated by the Joint Astronomy Centre on behalf of the Science and Technology Facilities Council of the United Kingdom, the Netherlands Organisation for Scientific Research, and the National Research Council of Canada.
  9. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under agreement by the Associated Universities, Inc.
  10. Support for CARMA construction was derived from the states of California, Illinois, and Maryland, the Gordon and Betty Moore Foundation, the Kenneth T. and Eileen L. Norris Foundation, the Associates of the California Institute of Technology, and the National Science Foundation. Ongoing CARMA development and operations are supported by the National Science Foundation under a cooperative agreement, and by the CARMA partner universities.
  11. See the TIMMI2 documentation for more information.
  12. While it should be noted that the scale factor determined in this manner is slightly dependent on the electron temperature, it changes by % for temperatures in the range of 7500-15000 K.
  13.  Data originally published by Wood & Churchwell (1989).
  14.  Line free channels from both IFs combined.
  15. All images have been restored to the beamsize of the 3 mm image:  at P.A. = +35°.
  16. J2000 position: 18:00:30.402, -24:04:01.39
  17. J2000 position: 18:00:30.463, -24:04:00.97
  18. values taken from the free-free subtracted image with resolution of
  19. footnotemark:
  20. source is unresolved
  21. source is unresolved
  22. uncertainties include the 15% calibration uncertainty
  23. comprised of 4 hyperfine components
  24. blended with adjacent lines
  25. blended with adjacent lines
  26. blended with adjacent lines
  27. blended with adjacent lines
  28. blended with adjacent lines
  29. blended with adjacent lines
  30. blended with adjacent lines
  31. blended with adjacent lines
  32. The conversion from flux density to brightness temperature is approximately 15.0 K/Jy in both sidebands.
  33. The conversion from flux density to brightness temperature is approximately 15.0 K/Jy in both sidebands.
  34. Measured at the positions in Table 4; includes absolute calibration uncertainty
  35. comprised of 14 components
  36. Assuming the abundance ratio of CO to H to be


  1. Acord, J. M., Churchwell, E., & Wood, D. O. S. 1998, ApJ, 495, L107
  2. Acord, J. M., Walmsley, C. M., & Churchwell, E. 1997, ApJ, 475, 693
  3. Akeson, R. L., & Carlstrom, J. E. 1996, ApJ, 470, 528
  4. Altenhoff W., Mezger P.G., Wendker H., Westerhout G., 1960, Veröff. Univ.-Sternwarte Bonn 59, 48
  5. Alvarez, C., Feldt, M., Henning, T., Puga, E., Brandner, W., & Stecklum, B. 2004, ApJS, 155, 123
  6. Apai, D., Linz, H., Henning, T., & Stecklum, B. 2005, A&A, 434, 987
  7. Argon, A. L., Reid, M. J., & Menten, K. M. 2003, ApJ, 593, 925
  8. Ball, R., Arens, J. F., Jernigan, J. G., Keto, E., & Meixner, M. M. 1992, ApJ, 389, 616
  9. Beckert, T., Duschl, W. J., & Mezger, P. G. 2000, A&A, 356, 1149
  10. Beltran, M. T., Estalella, R., Girart, J. M., Ho, P. T. P., & Anglada, G. 2007, ArXiv e-prints, 712, arXiv:0712.1757
  11. Beuther, H., Walsh, A. J., Thorwirth, S., Zhang, Q., Hunter, T. R., Megeath, S. T., & Menten, K. M. 2007, A&A, 466, 989
  12. Bik, A., Kaper, L., & Waters, L. B. F. M. 2006, A&A, 455, 561
  13. Brogan, C. L., Hunter, T. R., Indebetouw, R., Chandler, C. J., Shirley, Y. L., Rao, R., & Sarma, A. P. 2007a, Ap&SS, 380
  14. Brogan, C. L., Chandler, C. J., Hunter, T. R., Shirley, Y. L., & Sarma, A. P. 2007b, ApJ, 660, L133
  15. Bronfman, L., Nyman, L.-A., & May, J. 1996, A&AS, 115, 81
  16. Cesaroni, R., Walmsley, C. M., Koempe, C., & Churchwell, E. 1991, A&A, 252, 278
  17. Choi, M., Evans, N. J., II, & Jaffe, D. T. 1993, ApJ, 417, 624
  18. Churchwell, E., Wolfire, M. G., & Wood, D. O. S. 1990, ApJ, 354, 247
  19. Churchwell, E. 2002, ARA&A, 40, 27
  20. Churchwell, E. 1997, ApJ, 479, L59
  21. Comerón, F., Pasquali, A., & Torra, J. 2006, A&A, 457, 553
  22. Comito, C., Schilke, P., Endesfelder, U., Jiménez-Serra, I., & Martín-Pintado, J. 2007, A&A, 469, 207
  23. Cyganowski, C. J., Brogan, C. L., & Hunter, T. R. 2007, AJ, 134, 346
  24. Demircan, O., & Kahraman, G. 1991, Ap&SS, 181, 313
  25. Dickinson, C., Davies, R. D., & Davis, R. J. 2003, MNRAS, 341, 369
  26. Di Francesco, J., Johnstone, D., Kirk, H., MacKenzie, T., & Ledwosinska, E. 2008, ApJS, 175, 277
  27. Doty, S. D., Schöier, F. L., & van Dishoeck, E. F. 2004, A&A, 418, 1021
  28. Eisner, J. A., & Carpenter, J. M. 2006, ApJ, 641, 1162
  29. Faison, M., Churchwell, E., Hofner, P., Hackwell, J., Lynch, D. K., & Russell, R. W. 1998, ApJ, 500, 280
  30. Feldt, M., et al. 2003, ApJ, 599, L91
  31. Feldt, M., Stecklum, B., Henning, T., Launhardt, R., & Hayward, T. L. 1999, A&A, 346, 243
  32. Fish, V. L., Reid, M. J., Argon, A. L., & Zheng, X.-W. 2005, ApJS, 160, 220
  33. Fish, V. L., Reid, M. J., Wilner, D. J., & Churchwell, E. 2003, ApJ, 587, 701
  34. Flower, D. R., Pineau Des Forêts, G., & Walmsley, C. M. 2006, A&A, 456, 215
  35. Frerking, M. A., Langer, W. D., & Wilson, R. W. 1982, ApJ, 262, 590
  36. Genzel, R., & Downes, D. 1977, A&AS, 30, 145
  37. Gomez, Y., Rodriguez, L. F., Garay, G., & Moran, J. M. 1991, ApJ, 377, 519
  38. Hanson, M. M., Luhman, K. L., & Rieke, G. H. 2002, ApJS, 138, 35
  39. Harvey, P. M., & Forveille, T. 1988, A&A, 197, L19
  40. Harvey, P. M., Lester, D. F., Colome, C., Smith, B., Monin, J.-L., & Vauglin, I. 1994, ApJ, 433, 187
  41. Hatchell, J., Thompson, M. A., Millar, T. J., & MacDonald, G. H. 1998, A&AS, 133, 29
  42. Hilditch, R. W., & Bell, S. A. 1987, MNRAS, 229, 529
  43. Hofner, P., & Churchwell, E. 1996, A&AS, 120, 283
  44. Hofner, P., Wyrowski, F., Walmsley, C. M., & Churchwell, E. 2000, ApJ, 536, 393
  45. Hunter, T. R., Churchwell, E., Watson, C., Cox, P., Benford, D. J., & Roelfsema, P. R. 2000, AJ, 119, 2711
  46. Hunter, T. R., Brogan, C. L., Megeath, S. T., Menten, K. M., Beuther, H., & Thorwirth, S. 2006, ApJ, 649, 888
  47. Jørgensen, J. K., Johnstone, D., van Dishoeck, E. F., & Doty, S. D. 2006, A&A, 449, 609
  48. Fuente, A. 2008, Ap&SS, 313, 135
  49. Klaassen, P. D., Plume, R., Ouyed, R., von Benda-Beckmann, A. M., & Di Francesco, J. 2006, ApJ, 648, 1079
  50. Kraemer, K. E., & Jackson, J. M. 1995, ApJ, 439, L9
  51. Kurtz, S., Hofner, P., & Álvarez, C. V. 2004, ApJS, 155, 149
  52. Liljeström, T., & Gwinn, C. R. 2000, ApJ, 534, 781
  53. Mac Low, M.-M., Toraskar, J., Oishi, J. S., & Abel, T. 2007, ApJ, 668, 980
  54. Mangum, J. ⁢ Shirley, Y. 2006,
  55. Mangum, J. G., & Wootten, A. 1994, ApJ, 428, L33
  56. Martín-Hernández, N. L., Bik, A., Kaper, L., Tielens, A. G. G. M., & Hanson, M. M. 2003, A&A, 405, 175
  57. Mookerjea, B., Casper, E., Mundy, L. G., & Looney, L. W. 2007, ApJ, 659, 447
  58. Neri, R., et al. 2007, A&A, 468, L33
  59. Puga, E., Feldt, M., Alvarez, C., Henning, T., Apai, D., Le Coarer, E., Chalabaev, A., & Stecklum, B. 2006, ApJ, 641, 373
  60. Sandell, G. 1994, MNRAS, 271, 75
  61. Scheuer, P. A. G. 1960, MNRAS, 120, 231
  62. Sollins, P. K., et al. 2004, ApJ, 616, L35
  63. Stark, D. P., Goss, W. M., Churchwell, E., Fish, V. L., & Hoffman, I. M. 2007, ApJ, 656, 943
  64. Thompson, M. A., Hatchell, J., Walsh, A. J., MacDonald, G. H., & Millar, T. J. 2006, A&A, 453, 1003
  65. Thompson, M. A., & MacDonald, G. H. 1999, A&AS, 135, 531
  66. Velázquez, P. F., Dubner, G. M., Goss, W. M., & Green, A. J. 2002, AJ, 124, 2145
  67. Watson, C., Churchwell, E., Zweibel, E. G., & Crutcher, R. M. 2007, ApJ, 657, 318
  68. Watson, A. M., Coil, A. L., Shepherd, D. S., Hofner, P., & Churchwell, E. 1997, ApJ, 487, 818
  69. Wood, D. O. S. 1993, Massive Stars: Their Lives in the Interstellar Medium, 35, 108
  70. Wood, D. O. S., & Churchwell, E. 1989, ApJS, 69, 831
  71. Wright, M. C. H., Plambeck, R. L., & Wilner, D. J. 1996, ApJ, 469, 216
  72. Zhang, Q., & Ho, P. T. P. 1995, ApJ, 450, L63
  73. Zhang, Q., Hunter, T. R., Sridharan, T. K., & Cesaroni, R. 1999, ApJ, 527, L117
  74. Zijlstra, A. A., & Pottasch, S. R. 1988, A&A, 196, L9
  75. Zijlstra, A. A., Pottasch, S. R., Engels, D., Roelfsema, P. R., Hekkert, P. T.-L., & Umana, G. 1990, MNRAS, 246, 217
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description