Disk and outflow signatures in Orion-KL: The power of high-resolution thermal infrared spectroscopyBased on observations of the ESO program 380.C-0380(A).

Disk and outflow signatures in Orion-KL: The power of high-resolution thermal infrared spectroscopythanks: Based on observations of the ESO program 380.C-0380(A).

H. Beuther Max-Planck-Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
name@mpia.de
   H. Linz Max-Planck-Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
name@mpia.de
   A. Bik Max-Planck-Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
name@mpia.de
   M. Goto Max-Planck-Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
name@mpia.de
   Th. Henning Max-Planck-Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
name@mpia.de
Key Words.:
Stars: formation – Stars: early-type – Accretion, accretion disks – Techniques: spectroscopic – ISM: jets and outflows – Stars: individual: Orion-BN, Orion source n
Abstract

Context:The Orion-KL region contains the closest examples of high-mass accretion disk candidates. Studying their properties is an essential step in studying high-mass star formation.

Aims:Resolving at high spatial and spectral resolution the molecular line emission in the immediate environment of the exciting sources to infer the physical properties of the associated gas.

Methods:We used the CRIRES high-resolution spectrograph mounted on the VLT to study the ro-vibrational CO/CO, the Pfund  and H emission between 4.59 and 4.72 m wavelengths toward the BN object, the disk candidate source n, and a proposed dust density enhancement IRC3.

Results:We detected CO absorption and emission features toward all three targets. Toward the BN object, the data partly confirm the results obtained more than 25 years ago by Scoville et al., however, we also identify several new features. While the blue-shifted absorption is likely due to outflowing gas, toward the BN object we detect CO in emission extending in diameter to  AU with a velocity structure close to the . Although at the observational spectral resolution limit, the CO line width of that feature increases with energy levels, consistent with a disk origin. If one attributes the extended CO emission also to a disk origin, its extent is consistent with other massive disk candidates in the literature. For source n, we also find the blue-shifted CO absorption likely from an outflow. However, it also exhibits a narrower range of redshifted CO absorption and adjacent weak CO emission, consistent with infalling motions. We do not spatially resolve the emission for source n. For both sources we conduct a Boltzmann analysis of the CO absorption features and find temperatures between 100 and 160 K, and H column densities of the order a few times  cm. The observational signatures from IRC3 are very different with only weak absorption against a much weaker continuum source. However, the CO emission is extended and shows wedge-like position velocity signatures consistent with jet-entrainment of molecular gas, potentially associated with the Orion-KL outflow system. We also present and discuss the Pfund  and H emission in the region.

Conclusions:This analysis toward the closest high-mass disk candidates outlines the power of high spectral and spatial resolution mid-infrared spectroscopy to study the gas properties close to young massive stars. We will extend qualitatively similar studies to larger samples of high-mass young stellar objects to constrain the physical properties of the dense innermost gas structures in more detail also in a statistical sense.

1 Introduction

Understanding the physical structure of massive accretion disks is one of the main unsolved problems in high-mass star formation. Although indirect, the main line of arguments for accretion disks stems from massive molecular outflow observations that identify collimated and energetic outflows from high-mass young stellar objects (YSOs, e.g., Henning et al. 2000; Beuther et al. 2002b; Zhang et al. 2005; Arce et al. 2007). Collimated jet-like outflow structures are usually attributed to massive accretion disks and magneto-centrifugal acceleration. Recent 2D and 3D magneto-hydrodynamical simulations of massive collapsing gas cores also result in the formation of massive accretion disks (Yorke & Sonnhalter, 2002; Krumholz et al., 2007, 2009). However, it is still unclear whether such massive disks are similar to their low-mass counterparts, hence dominated by the central YSO and in Keplerian rotation, or whether they are maybe self-gravitating non-Keplerian entities.

While studies at (sub)mm wavelengths are a powerful tool to mainly study the cold gas and dust components on spatial scales of the order 1000 AU (e.g., Cesaroni et al. 2007), such observations are not that well suited to investigate the inner and warmer components of massive rotating structures. In contrast to that, mid-infrared spectral lines, e.g., ro-vibrationally excited CO emission lines, can trace these warm gas components. However, the spectral and/or spatial resolution was mostly lacking because absorption features were dominating and hence prohibiting the detection of the accretion disks in emission. Several recent studies have further demonstrated the power of high-spectral and high-spatial resolution CO mid-infrared spectroscopy for disks around low-mass young stellar objects and Herbig Ae stars (e.g., Goto et al. 2006; Pontoppidan et al. 2008; van der Plas et al. 2009).

To achieve the highest angular and spectral resolution possible, we observed some of the closest massive disk candidates in Orion at a distance of 414 pc (Menten et al., 2007) – Orion-BN (the Becklin-Neugebauer Object), source n and IRC3 – in the CO transitions around 4.65 m with the CRyogenic high resolution InfraRed Echelle Spectrograph (CRIRES, Käufl et al. 2004) at the VLT. Both objects are well detected at mid-infrared wavelengths exhibiting various kinds of disk-signatures.

The Becklin-Neugebauer Object: Since its detection in the 1960s, the BN object is one of the archetypical high-mass YSOs (Becklin & Neugebauer, 1967; Henning et al., 1990). Scoville et al. (1983) observed the source in several frequency settings between 2 and 5 m and detected molecular emission from several CO isotopologues (fundamental and overtone emission), and they inferred that BN exhibits an outflow/wind as well as a highly confined region of molecular gas at high densities and temperatures of 3500 K. The estimated luminosity of the BN object is  L corresponding to a B0.5 main sequence star (Scoville et al., 1983). More recently, Jiang et al. (2005) observed BN in polarized near-infrared emission, and they also identified signatures caused by an embedded accretion disk. The BN object has a high velocity along the line of sight of 21 km s compared with the cloud velocity of around 5 km s (Scoville et al., 1993). This is consistent with the measured high proper motions of that object (e.g., Plambeck et al. 1995). Whether the BN object is expelled from the Trapezium system or during a disintegration of a bound system once containing source I, source n and the BN object itself is still a matter of debate (e.g., Tan 2004; Gómez et al. 2005; Zapata et al. 2009).

Source n: Based on a bipolar radio morphology and HO maser association, Menten & Reid (1995) suggested that this source may be one of the driving sources of the powerful molecular outflows within Orion-KL. Extended mid-infrared emission was observed perpendicular to the outflow axis (Greenhill et al., 2004; Shuping et al., 2004), and Luhman (2000) detected CO overtone emission. Both features are interpreted as likely being due to an accretion disk. Source n is believed to be in an evolutionary younger stage than the BN object, and the luminosity is estimated to be lower as well, of the order 2000 L (Greenhill et al., 2004).

IRC3: The source n observations serendipitously covered the extended infrared source IRC3 (e.g., Dougados et al. 1993) which we present here as well. At 3.6 m wavelengths, IRC3 is elongated in the northeast-southwest direction (Dougados et al., 1993) and shows highly polarized near- to mid-infrared emission (Minchin et al., 1991). The observations are consistent with IRC3 being a dust density enhancement reprocessing light from another source, potentially IRC2 (Downes et al., 1981; Minchin et al., 1991; Dougados et al., 1993).

Figure 1 gives an overview of the region marking the sources discussed in the paper as well as the slit orientations (see also section 2). The nominal absolute positions for the three sources are listed in Table 1.

Figure 1: Overview image of the region. The grey-scale band image is taken from UKIDSS DR5 (Lawrence et al., 2007), the white stellar peak positions are due to saturation. The slits and important sources are marked. The cross denotes the position of IRc2.
Source R.A. Dec.
(J2000.0) (J2000.0)
BN object 05h35m14.12s -05d22m22.9s
Source n 05h35m14.35s -05d22m32.9s
IRC3 05h35m13.90s -05d22m30.0s
Table 1: Source positions (from Dougados et al. 1993)

2 Observations

We obtained high-resolution spectra between 4.6 and 4.7 m with CRIRES (Käufl et al., 2004)) mounted on UT1 at the VLT on Paranal, Chile. Two grating settings were selected (12/-1/n, =4662.1 and 12/-1/i, =4676.1) to observe the spectral interval covering the CO  [P(1)–P(5)/R(0)–R(8)] lines without gaps.

One slit covered the BN object with a position angle of 126 degrees east of north, whereas the second slit included source n and IRC3 with a position angle of 110 degrees east of north. For the strong source BN, only 40 secs (DIT=2 secs, NDIT=10) were required. For the second object with the weaker sources we had a total on-slit integration time of 20 minutes (DIT=10 secs, NDIT=2). For the BN observation a slit width of was used while for the observation of source n a slit with of was selected which correspond to spectral resolving powers of 100,000 and 50,000, respectively. The non-AO mode was applied since no natural guide star is available in the close environment. The infrared seeing measured from the spectra was 0.35″ during the BN observation and 0.45″ for the source n observation. The nod-throw of all the observations was set to 10″. To correct for the telluric absorption lines, attached to every science observation a telluric standard star (HR 1666 with spectral type A3III) was observed.

As extended emission in the CO lines was present in the observations, we corrected the frames for distortion in order to get a wavelength solution valid for the whole chip. Firstly, the chips are slightly rotated with respect to the slit (ranging from 0.05 degree for chip 3 to 0.45 degree in the case of chip 4). We measured the position of the brightest object as function of wavelength and calculate from the displacement in position the rotation angle. Secondly, after the rotation angle has been corrected, the curvature of the slit is corrected by measuring the position (central wavelength) of a sky-line as function of the spatial coordinate. This could be described with a 2nd degree polynomial. Using the IDL routines polywarp and poly_2d, the distortion was corrected.

After the distortion correction, the raw files are processed by the ESO CRIRES pipeline (version 1.10.1) in combination with the Gasgano software. The data are dark subtracted, flat field corrected as well as corrected for non-linearity. The wavelength calibration is done using the telluric emission lines in combination with a HITRAN model spectrum (Rothman et al., 2005). A cross-correlation of the spectra with the HITRAN spectra showed that the wavelength accuracy of the spectra is 0.5 km s.

Two absorption lines are present in the standard star: the Pfund  and the Humphreys  line. Due to the strong telluric absorption, these lines are not trivial to remove. Therefore, we first reduced the spectrum without correcting for the intrinsic absorption lines of the standard star. In the final reduced spectrum these lines become eminent as emission lines free from contamination by the atmosphere. The Pfund  line is also seen in our science object. However, the Humphreys  line is not present in the science spectrum of BN before division by the standard star. Therefore, it can be used to correct for the absorption lines of the standard star. We used a high resolution Kurucz model spectrum from an A0V star (http://kurucz.har vard.edu/stars.html) and scaled and shifted the model spectrum such that the Humphreys  profile would fit that of the observed Humphreys  profile in the standard star. Assuming that the Pfund  line scales in the same way as the Humphreys  line does, we divided with the model spectrum to remove the line contamination.

The spectra were corrected for the earth velocity to the local standard of rest using rvcorrect in IRAF. The velocity corrections applied for the two observing dates were -3.8 km s for the BN data (observed on 21st October 2007) and +43.2 km s for the source n/IRC3 observations (taken on 21st February 2008). The velocity relative to the local standard of rest of Orion varies between 2.5 and 9 km s (e.g., Comito et al. 2005), and we adopt the approximate value of +5 km s.

3 Results

Figure 2: CRIRES observations of the R and P CO line series around 4.65 m. The strong absorption features are the CO lines, and the emission line at 4.654 m is Pfund .

Toward all three sources we detected the whole suite of CO  lines present in the spectral window, the CO  lines from R(6) to R(13) that were not blended by the CO  lines, as well as the Pfund  line. Figure 2 presents the complete spectrum toward the BN object, and Table 2 gives an overview of the covered lines, their wavelengths and the lower-level energy state of the transitions (). In total, this setup covers a broad range of energy levels extending up to 504 K.

Line
m K
CO
R(9) 4.5876 249
R(8) 4.5950 199
R(7) 4.6024 155
R(6) 4.6090 116
R(5) 4.6176 83
R(4) 4.6254 55
R(3) 4.6333 33
R(2) 4.6412 17
R(1) 4.6493 5.5
R(0) 4.6575 0
P(1) 4.6742 5.5
P(2) 4.6826 17
P(3) 4.6912 33
P(4) 4.7000 55
P(5) 4.7088 83
CO
R(13) 4.6641 504
R(12) 4.6741 432
R(11) 4.6782 365
R(10) 4.6853 304
R(9) 4.6926 249
R(8) 4.7000 199
R(7) 4.7075 155
R(6) 4.7150 116
Pfund  4.6538
H S9 4.6947 7198

The temperatures listed for the Pfund  and H S9 lines do not indicate gas temperatures since these lines are not excited by collisions within thermal gas, but rather by ionizing photons and within shocked regions, respectively.

Table 2: Observed lines

3.1 The BN object

Figure 3: The CO R(8) to R(1) lines from top to bottom toward the BN object. Telluric and velocities are marked. The feature between 0 and 8 km s is a telluric artifact.

3.1.1 CO absorption and emission

Figure 3 presents a zoom compilation of CO data from the R(1) to R(8) line covering lower energy levels between 5.5 and 199 K. If one ignores the telluric line feature at 3 km s two broad absorption features can be identified at approximately and +8 km s. Since all CO lines are saturated, their peak absorption velocities are unreliable, and we refer to the CO data (Fig. 4). The peak velocity of the blue-shifted component is  km s extending from to  km s. The second absorption feature has its peak at  km s, close to the of the cloud, and extends from to  km s. Furthermore, red-shifted from the absorption we clearly identify a CO emission peaking in CO and CO at  km s and extending from 15 to 30 km s. The overall extent of the CO absorption and emission is from to  km s. It should be noted that while the of the different cloud components for the Orion-KL region vary between approximately 3 and 9 km s, Scoville et al. (1983) inferred that the corresponding velocity of the BN object is significantly larger around 21 km s (consistent with the different BN ejection scenarios, e.g., Tan 2004; Gómez et al. 2005). While the absorption features stem from the warmer protostellar envelopes and the surrounding cloud with usual temperatures of the order 100 K (see also Boltzmann analysis below), the ro-vibrational lines in emission can be caused by different processes. For example, fluorescence via UV photons or resonance scattering from strong infrared fields can excite these lines without significantly heating the gas (e.g., Blake & Boogert 2004; Ryde & Schöier 2001). Alternatively, the ro-vibrational lines could be caused by hotter gas components (see discussion in section 4). We note that the critical densities of these lines are of order  cm which practically implies that extended gas components can be hardly responsible for the emission.

Although the blue-shifted part of the spectrum with respect to the of the molecular cloud is slightly broader than the red-shifted part, nevertheless we clearly identify red-shifted absorption as well. To first order, the blue-shifted gas seen in absorption can be identified with outflowing gas from the region, whereas the red-shifted features belong to gas infalling in the direction of the central source. The outflowing gas with a maximum velocity relative to the of  km s is consistent with outflow wings often observed at mm wavelengths from young massive star-forming regions (e.g., Beuther et al. 2002b). It should be noted that the even broader outflow wings observed at mm wavelength toward Orion-KL exceeding  km s (e.g., Chernin & Wright 1996) are likely not related to the BN object but rather to one or more sources about south-east of BN (source I, source n and/or SMA1, e.g., Greenhill et al. 2004; Jiang et al. 2005; Bally 2008; Beuther & Nissen 2008).

Figure 4: CO optical depths toward the BN object.

How do these general features compare with the data published by Scoville et al. (1983) which were observed during several observing runs between December 1977 and February 1981, hence about 30 years prior to our observations. The general CO line structure with two strong absorption features plus one red-shifted emission peak are largely the same. Also the overall extent of the emission is quite similar. Compared to our measured values of , and  km s for the three components, Scoville et al. (1983) report for the corresponding features velocities of , and  km s. While two velocities agree well, the most blue-shifted absorption peak appears to have shifted a little bit between the two observations. However, given that their spectral resolution was more than a factor 2 lower than that of the new CRIRES data (7 versus 3 km s), we refrain from further interpretation of this. Furthermore, Scoville et al. (1983) identify two more absorption features, one at  km s and one at +30 km s. Regarding the  km s component we are not able to infer any changes because that features lies very close to the telluric emission which obscures any reliable signature there. However, Fig. 3 shows that we do not detect any additional absorption feature blue-shifted from the 26.6 km s emission. Therefore, the +30 km s absorption dip reported by Scoville et al. (1983) was either a transient feature or not significant with respect to the signal-to-noise ratio.

Since the CO data are so strongly saturated, for the following analysis we use the corresponding CO data covering the R(6) to R(13) transitions (Fig. 2). Since the absorption depth is related to the optical depth via we can directly estimate the optical depth of the CO absorption lines shown in Fig. 4 if the lines are spectrally resolved. Since the full width at zero intensity (FWZI) is of order 20 km s (see above), this criterion is fulfilled with our spectral resolution which corresponds to a velocity resolution of 3 km s. Except of the lowest CO R(6) line component at  km s all other observed CO absorption features have optical depths below 1. This allows us to estimate rotational temperatures via Boltzmann plots from the equivalent line width following the approach outlined in Scoville et al. (1983) also adopting their finite optical depth corrections. Using their equation (A8), the column density of the lower-level energy state is

with , and denoting the optical depth correction following the appendix in Scoville et al. (1983), the equivalent line width (in units of cm) and absorption oscillator strength, respectively. Because the absorption feature at 9 km s is affected by the telluric line, in the following we only work with the data from the  km s line. The parameters and the calculated column densities are listed in Table 3. In a Boltzmann distribution, the total column density is related to the temperature and the lower state column density via:

with the statistical weight of the lower-level transition and the partition function at the given temperature.

K km/s cm
BN@-15 km s
R(6) 116 3.70 6.0 1.18 13 6.0
R(7) 155 3.16 5.9 1.14 15 5.0
R(9) 249 1.90 5.9 1.08 19 2.9
R(10) 304 1.36 5.8 1.06 21 2.0
R(11) 365 0.69 5.8 1.04 23 1.0
R(12) 432 0.51 5.8 1.02 25 0.7
R(13) 504 0.27 5.8 1.02 27 0.4
source n @-7 km s comp.
R(7) 155 4.1 5.9 1.14 15 6.5
R(9) 249 2.7 5.9 1.08 19 4.1
R(10) 304 2.0 5.8 1.08 21 3.1
R(12) 432 1.3 5.8 1.06 25 2.0
R(13) 504 1.5 5.8 1.05 27 2.2
source n @+5 km s comp.
R(7) 155 4.6 5.9 1.23 15 7.8
R(9) 249 2.7 5.9 1.16 19 4.4
R(10) 304 1.8 5.8 1.12 21 2.8
R(12) 432 0.6 5.8 1.04 25 0.9
R(13) 504 0.1 5.8 1.00 27 0.1

The Table entries are explained in the main text.

Table 3: CO values for Boltzmann plots

Figure 5: Boltzmann plot for the CO v=1-0 lines toward the BN object. The x-axis shows the lower-level energies of the transitions and the y-axis presents the natural logarithm of the corresponding column densities divided by their statistical weights.

Figure 5 presents the corresponding Boltzmann plot where the lower-state column density divided by the statistical weight is plotted against the lower-state energy divided by . A linear fit to the data gives the rotational temperature and the column density of CO divided by the partition function . For the  km s CO component the fitted is  K. Scoville et al. (1983) calculated the rotation temperature for the more blue-shifted absorption component from the CO transitions, and they find 150 K there, which is approximately consistent with our new determination. The other fitted parameter is  cm. For 2-atomic molecules like CO the partition function can be approximated by (where is the rotation constant), and we get a CO column density of  cm. Using furthermore the CO to CO isotopologic ratio of 69 (e.g., Sheffer et al. 2007), we derive a total CO column density of the  km s component toward Orion-BN of  cm. These column density estimates are in excellent agreement with the results derived by Scoville et al. (1983) from the CO  lines.

How do these values compare to other observations? The BN object was detected at 1.3 mm wavelength by Blake et al. (1996) at a 0.15 Jy level with a spatial resolution of . Following Plambeck et al. (1995), about 50 mJy of that flux can be attributed to circumstellar dust emission. Assuming optically thin dust emission, an average dust temperature of 100 K, a dust opacity index of 2 and a standard gas-to-dust ratio of 100, we can calculate the total H column density (Hildebrand, 1983; Ossenkopf & Henning, 1994; Henning et al., 1995; Beuther et al., 2002a, 2005a; Draine et al., 2007). The derived H column density and mass within their syntesized beam are then  cm and 0.1 M, respectively. For comparison, using a standard CO-to-H ratio of , the H column density estimated here from the CO IR spectroscopy data is  cm. While the mm continuum data trace all velocity components along the line of sight, the near-infrared data trace selected velocity components. Judging from Figures 3 and 4, where we see more than 1 velocity component, the total column density traced by the CO data is more than a factor 2 higher. Furthermore, the CO absorption only traces the gas in the foreground of the near-infrared source reducing the traced gas by another factor 2. While both approaches are affected by systematics – e.g., the mm derived column densities can be wrong by a factor 5 depending on the assumptions on the dust properties and temperatures – other reasons are more important for some of the differences. In particular, the mm continuum emission is sensitive to the cold and warm dust emission whereas the near-infrared absorption of the CO transitions useable for our analysis traces mainly the warmer gas. Therefore, these data indicate that a large fraction of the gas is at relatively low temperatures.

3.1.2 Potential signatures from the BN disk?

While most of the observed emission is spatially unresolved, we find weak extended CO emission toward the BN object. Figure 6 presents a position-velocity diagram along the slit direction (PA of 126 east of north, see sec. 2) which corresponds to the proposed disk orientation (Jiang et al., 2005). We find CO emission around the rest velocity of the BN object of  km s extending approximately in both directions. Since we do not identify any strong velocity dispersion, this extended emission is unlikely to be due to an outflow, but it may potentially come from the proposed disk (Jiang et al., 2005). The measured extent of would then correspond at the given distance of Orion of 414 pc to an approximate disk diameter of  AU or a disk radius of 1650 AU. While such disk size would be comparably large with respect to typical low-mass disks (e.g., several reviews in Reipurth et al. 2007), it is consistent with measured sizes of rotating structures in other high-mass star-forming regions (e.g., Schreyer et al. 2002; Cesaroni et al. 2005; Beltrán et al. 2006; Beuther & Walsh 2008).


Figure 6: Position velocity diagram of the CO (R0) line along the slit direction (PA of 126 east of north, see sec. 2) plotted in logarithmic intensity scale. Around offset 0 the figure is dominated by the continuum emission, but outside the emission stems from the CO line. Positive offsets go in northwest direction.

While we do not identify extended emission is the rarer CO isotopologue, the CO & CO emission feature around 20 km s exhibits an additional interesting feature in the CO data (Fig. 4): While the emission feature is at the edge of the bandpass for the CO R(6) line, it remains undetected for the R(7) line and only comes up for the lines greater R(9). This emission feature is also visible in all CO lines (R(1) to R(8)) shown in Figure 3, but it does not exhibit any significant shape change there. Likely, the CO emission is optically thick and traces only an outer envelope that does not show big variations with excitation temperature. In contrast to that, the measured Gaussian line width of that component for the more optically thin CO R(9), R(10), R(11), R(12) and R(13) lines are 4.3, 5.3, 6.5, 6.4 and 6.6 km s, respectively. Although errors on the line width are difficult to quantify because the emission is at the flank of the strong absorption feature, and furthermore our spectral resolution is only 3 km s, the data are indicative of a line width increases with excitation temperature of the line that appears to saturate for the highest detectable transitions R(11) to R(13). While for a centrifugally supported disk with a central mass of 10 M the circular velocity at 3300 AU is relatively small of order  km s, the measured line width increase is consistent with a rotating disk structure because the inner region with larger rotation velocities should have higher temperatures caused by the exciting central star. One has to keep in mind that for purely thermal excitation, rotation would not cause such an effect because the upper levels – all in the 3000 K regime – do not exhibit big relative excitation differences. Hence, they do not trace significantly different regions of a disk then. In contrast to that, for lines that are excited by UV fluorescence such an effect is possible because the excitation mechanism via an electronic excited state tends to preserve the level population of the rotational states (e.g., Brittain et al. 2009). Therefore, in this case the R(9) to R(13) lines are sensitive to gas temperatures in a relatively broader range between 150 and 500 K (Table 2). An exact reproduction of the line width increase would require a detailed disk model, including its density and temperature structure. Since we cannot constrain these parameters well from our data, this is beyond the scope of this paper. Although the line width increase is below our nominal spectral resolution element, qualitatively the observations are consistent with a rotating disk structure.

3.1.3 The Pfund  line toward the BN object and H emission nearby BN

Figure 7 presents a zoom into the Pfund  line. The line shape consists of two components, one central Gaussian component and broad line wings. It is possible to fit the whole profile relatively well with a two-component Gaussian fit where the central component has a FWHM  km s and the broad component has a FWHM of  km s. The full width down to zero intensity is approximately 340 km s (between -170 and 170 km s). Close to the peak of the profile at around 25 km s, the spectrum exhibits a small dip which is likely an artifact from the telluric corrections (see section 2). Considering this, the peak of the central Gaussian fit at 16.5 km s is still consistent with the velocity derived for the BN object by Scoville et al. (1983) of 21 km s.

Figure 7: Hydrogen recombination Pfund  line toward the BN object. The histogram shows the data, and the three other lines present a two-component Gaussian fit to the data. Both Gaussians are shown separately as well as the combined fit. The FWHM of both components is given in the figure. The big dots reproduce the best fit to the Br line (their Figure 7 re-scaled to our normalized Pfund  spectrum) for an optically thin outflow with a velocity law by Scoville et al. (1983).

Figure 7 also presents as thick dots the best fit obtained by Scoville et al. (1983) for the Br line (their Figure 7 re-scaled to our normalized spectrum). Their model consists of a supersonic, optically thin decelerating outflow with a velocity law of . It is remarkable how well the shape of their Br line obtained 30 years ago corresponds to the shape of the newly observed Pfund  line. Hence, these Pfund  data are also consistent with their outflow model.

Figure 8: Position velocity digram of the H emission south-east of the BN-object associated with the infrared source IRC15 (Dougados et al., 1993).

Furthermore, approximately south-east of the BN object we detect H emission from the H S9 transition with  K. This H emission feature is spatially associated with the corresponding H knot in near-infrared H images (e.g., Nissen et al. 2007) as well as with the mid-infrared source IRc15 reported by Shuping et al. (2004). Figure 8 shows a position velocity diagram of this feature. While we do not resolve any spatial structure of the H knot, it shows a very broad velocity extent going to blue-shifted velocities of  km s, in excess of the CO absorption features measured toward the BN object. Since we are mainly interested in the BN object, source n and IRC3, we refrain from further analysis of this offset H emission.

3.2 Source n

3.2.1 CO absorption and emission

Due to the different observing dates, for source n (and IRC3 in the following section) the telluric lines are almost entirely shifted out of the CO spectrum. Figures 9 and 10 present the corresponding CO R-series lines and the detected CO lines with their associated optical depths. In comparison to the previous BN data, for source n we only identify one broad absorption feature in the CO data which is dominated by blue-shifted outflowing gas. It is interesting to note that the most blue-shifted absorption does not have a Gaussian shape but rather a more extended wing-like structure. In the framework of accelerated winds with distance from the driving star or disk, such a spectral behavior would be expected. These accelerated winds increase in velocity with distance from the star. Simultaneously, with increasing distance from the driving source the density of the surrounding gas and dust envelope also decreases, lowering the corresponding absorption depth at those velocities. Hence, in this picture, higher-velocity gas should show shallower absorption features in the spectra (e.g., Lamers & Cassinelli 1999). Furthermore, we also identify a red-shifted emission component at +15 km s. However, compared to the BN object, where the emission feature can be identified in all transitions, for source n it is more prominent in the higher excited lines. The total width of the CO absorption and emission is 90 km s ranging from approximately -65 to +25 km s (broader than for BN). The blue-shifted end is less well determined because of telluric line contamination. Nevertheless the blue-shifted outflow part of the spectrum extends  km s from the assumed of +5 km s. This value exceeds that measured toward BN by about 20 km s. Although at the edge of the telluric line contamination, we tentatively identify an additional discrete absorption feature at -35 km s.

Figure 9: The CO R(8) to R(1) lines from top to bottom toward source-n. Telluric and velocities are marked. The linear slopes and dips around -40 km s associated with the R(8), R(6) and R(4) lines (for the latter extending from -5 to -65 km s) are artifacts due to the telluric corrections.
Figure 10: CO optical depths toward source-n. The absorption feature at +25 km s of the R(7) line is likely an artifact due to telluric lines.

Since the CO line is again always saturated (Fig. 9), for the additional analysis we use the corresponding CO data (Fig. 10). Because of line-blending and telluric emission, we only measure five CO lines for source n. The broad, saturated CO absorption feature is clearly resolved into two components with approximate velocities at -7 km s and +5 km s. While the +5 km s shows the deeper absorption features for the lower energy R-lines, it is interesting to note that the single absorption peak observed for the CO R(13) line is associated with the -7 km s component, indicating hotter gas at more blue-shifted velocities. Following the approach outlined for BN in section 3.1.1, the derived optical depth of CO is 0.6 at its highest.

Since the telluric lines are far away in velocity space (see Fig. 9), for source n we can conduct the Boltzmann analysis for both absorption components separately. The measured equivalent widths for both components are listed in Table 3111The equivalent width is measured by simultaneous Gaussian fits to both absorption features.. Again calculating the CO v=1 column densities (Table 3) and producing Boltzmann plots (Fig. 11), we find that for both velocity components fits to all 5 data points are less good than in the case of the BN object. This is likely due to the fact that fitting 2 Gaussians to the broad R(13) line, which does not show two well separated absorption features anymore, may overestimate the contribution of the -7 km s component and hence underestimate the contribution from the +5 km s component. Therefore, we also fit only the 4 lower-energy transitions between R(7) and R(12) resulting in better fits (Fig. 11). As expected from the different behavior of the two absorption components with increasing energy levels, the derived rotation temperatures using the 4 lower energy lines of the -7 km s and +5 km s are 163 K and 103 K, respectively. As a comparison, rotational temperatures measured at submm wavelength from CHOH emission lines toward source n are around 200 K (Beuther et al., 2005b). Considering that these measurements are conducted with different molecules and very different observational techniques, the overall range of similar temperatures derived at mid-infrared and submm wavelengths is reassuring for the complementarity of such multi-wavelengths observations.

Fitting only the 4 lower-energy lines, the column density of CO divided by the partition function results in  cm and  cm for the -7 km s and +5 km s components, respectively. Approximating again the partition function by (section 3.1.1) we have and , and using the CO to CO isotopologic ratio of 69 (e.g., Sheffer et al. 2007), the derived total CO column densities toward source n are  cm and  cm, of the same order as for the BN object. With the CO-to-H conversion factor of , the corresponding H column densities are  cm and  cm. Adding these two column density values, the total H column density traced toward source n by this near-infrared observations is  cm, more than an order of magnitude below the H column densities derived with the Submillimeter Array at 865 m ( cm, Beuther et al. 2004). Similar to the BN-case discussed in section 3.1.1, this shows that the near-infrared data mainly trace the warmer gas whereas the (sub)mm continuum observations trace warm and cold gas components.

Figure 11: Boltzmann plots for the CO v=1-0 lines toward source-n for the -7 km s and the +5 km s components in the upper and lower panels, respectively. The x-axis shows the lower-level energies of the transitions and the y-axis presents the natural logarithm of the corresponding column densities divided by their statistical weights. In both panels, the dashed line shows a fit to all five data points, whereas the full line presents a fit to only the four lower transitions.

3.2.2 The Pfund  line

The Pfund  line is also detected toward source n (Fig. 12), and we can fit a Gaussian to the recombination line with FWHM of  km s, a width down to 0 intensity of  km s (from to  km s) and a central velocity of 0 km s. The line is covered by 2 CRIRES chips, and while the general Pfund  line shape is the same for both chips, the central dip at 0 km s cannot independently be reproduced. Therefore, the dip is likely only an artifact due to insufficient signal to noise.

The thermal line width of a  K hydrogen gas is 21.4 km s. Convolving that with the 6 km s spectral resolution, the observable thermal line width should be 22.2 km s. Therefore, the measured line width does not exceed much the thermal line width of an Hii region. Hence, the Pfund  emission from source n does not exhibit strong signatures from a wind but is rather consistent with a thermal Hii region. We also cannot exclude that the Pfund  emission toward source n is contaminated by more broadly distributed Orion nebula emission.

Figure 12: Hydrogen recombination Pfund  line of source n.

3.3 Irc3

As outlined in the Introduction, in contrast to the BN object and source n, the source IRC3 is unlikely to be a YSO, and it is not prominent in any typical hot core tracer (e.g., Blake et al. 1996; Beuther et al. 2005b). The source rather resembles a more extended dust density enhancement that reprocesses light from other sources, potentially IRC2. Figure 13 presents a part of the 2D slit spectrum clearly showing the extended nature of the CO emission toward IRC3 in contrast to the point-like continuum structure from source n. The continuum emission from IRC3 is very weak compared to that from source n.

Figure 13: Original data from the slit covering source n and IRC3. Source n is the bright continuum source at offset whereas IRC3 exhibits the extended CO emission around offset .
Figure 14: The CO R(8) to R(1) lines (except of the R(4) line from top to bottom toward IRC3. Since the emission is extended, this spectrum is an average over along the slit. Telluric and velocities are marked.

Figure 14 shows the CO lines extracted as an average spectrum of length toward IRC3. While we again see a broad absorption feature against the weaker background extending from  km s out to the telluric contamination at about  km s, IRC3 shows a broad and strong emission feature peaking at  km s. Since this emission is extended over several arcseconds, Figure 15 presents the position velocity diagrams of selected CO and all CO lines. In contrast to the BN object where the extended CO emission is just around the rest velocity of the star (Fig. 6), here we see a clear trend of increasing velocity with increasing distance from the continuum peak (the so-called Hubble law of outflows). Furthermore, in particular the CO position velocity diagrams show a twofold structure with one velocity increase to values  km s at offsets of and another increase to similar velocities at offset from the main continuum peak. To guide the eye, these wedge-like structures are sketched in the top-right panel of Figure 15). Such multiple wedge position velocity structures are what jet-bow-shock entrainment models for molecular outflows predict (e.g., Arce et al. 2007). Furthermore, similar to source n, the CO absorption shows also for IRC3 the extended wing-like features toward the most blue-shifted absorption. Similar to the pv-diagrams, where we see acceleration of the gas with distance from the source, these wing-like absorption can also be interpreted in the framework of accelerated winds (see also section 3.2.1 or Lamers & Cassinelli 1999). In summary, IRC3 is not only a dust density enhancement scattering the light from another source, it may also be part of the famous outflow emanating from the Orion-KL region. Based on these data, we cannot clearly define the driving source of the CO outflow feature. While the continuum peak of IRC3 could be driving an outflow, the CO emission may also be part of the large-scale outflow where different driving sources have been proposed for in the literature, e.g., source I, SMA1, or the disintegration of a bound system once containing sources I, n and the BN-object (e.g., Menten & Reid 1995; Gómez et al. 2005; Beuther & Nissen 2008; Zapata et al. 2009). While the position velocity structure of the CO lines largely resembles that of the CO lines (Fig. 15), it is interesting to note that the strongest CO feature is at higher velocities and relatively distant from the continuum peak (offset and velocity  km s. This difference is likely caused by the varying optical depth of the CO and CO transitions. Since the velocity of the emission feature is very close to the of the cloud (typically between 5 and 9 km s, e.g., Comito et al. 2005), self-absorption of the most abundant CO isotopologue may veil this feature, whereas the less abundant CO can be utilised to penetrate more deeply into the cloud and hence detect the emission better.

Figure 15: Position velocity cuts along the slit axis for IRC3 with a PA of 110 degrees east of north. The top row shows diagrams for selected CO lines, and the bottom row presents the CO data. For CO the scaling is adapted to highlight the weaker extended emission in contrast to the continuum at offset . The wedge-like outflow structures are sketched in the top-right panel.

4 Discussion and conclusion

High spectral resolution mid-infrared observations allow us to infer several important characteristics of some of the major sources within the Orion-KL region. The fundamental CO lines are largely dominated by absorption from the individual YSO envelopes and their surrounding gas cloud. These absorption features are blue- and red-shifted with respect to the of the molecular cloud indicating that outflowing and inflowing gas are simultaneously present toward our target sources. However, we also identify interesting emission features. Together they confirm the youth of the sources where infall and likely accretion are still ongoing.

BN object:

For the BN object, our data confirm, at double the spectral resolution, several of the assessments conducted already by Scoville et al. (1983). However, also discrepancies arise. For example, an absorption feature reported previously around +30 km s is not found in the new data, indicating either transient components or poor signal-to-noise in the older data. From a Boltzmann analysis, the rotational temperature is around 112 K, and we derive CO column densities of several times  cm, well in agreement with the older results based on CO v=2-0 by Scoville et al. (1983). Using standard CO-to-H conversion factors, these column densities are about an order of magnitude below H column density estimates based on mm continuum emission. As discussed in section 3.1.1, while systematics may account for some of this discrepancy, the main difference is that the mm continuum emission is sensitive to cold and warm dust, whereas the near-infrared absorption lines trace only the warm gas components.

We also identify extended CO emission, likely emanating from the close environment of BN. While the absorption stems from warm gas with temperatures of the order 100 K, getting the ro-vibrationally lines in emission implies already much higher temperatures for the corresponding gas components. The velocity range of the CO emission between 15 and 30 km s encompasses the previously inferred velocity of the BN object of 21 km s. Following Scoville et al. (1983), who covered more J-transitions (fundamental and overtone emission) than we do, the rotational temperature of the compact emitting gas is around 600 K, significantly higher than the above derived rotational temperature of the absorbing envelope gas. As outlined in section 3.1.1, high critical densities are required to produce this line emission by pure thermal collisional excitation, and other processes like UV fluorescence or resonance scattering from infrared emission may also contribute to the emission lines. Scoville et al. (1983) infer an approximate size for that emission of 20 AU, whereas we now resolve the emission extending to , corresponding to a diameter of  AU. This discrepancy can be explained by the large dynamical range between the compact emission close to the source itself and the comparably very weak more extended emission. The small size found by Scoville et al. (1983) can be attributed to the strong compact emission which also emits in the overtone bands, whereas our more sensitive new data also detect the weaker extended features due to the high dynamic range available with CRIRES. While the absolute ratio of the peak emission (CO plus continuum) to the extended CO emission is 210, the ratio of the continuum subtracted CO emission toward BN compared to the extended CO is still 70. Scoville et al. (1983) were not sensitive enough to detect this faint extended emission compared to the strong central source. Although at the spectral resolution limit of our observations, the measured line width increase of the CO emission feature with increasing excitation temperature is consistent with a disk origin of the CO emission. Therefore, these parameters make the CO emitting gas potentially to be associated with a disk around the BN object (see also Jiang et al. 2005). Further supporting the disk-outflow interpretation for the BN object are the Pfund  hydrogen recombination line data which show broad high-velocity line wings consistent with the decelerating outflow scenario proposed by Scoville et al. (1983).

Source n:

While bright sources like the BN object were already feasible to be observed a while ago (e.g., Scoville et al. 1983) weaker sources like source n or IRC3 allowed reasonable high-spectral-resolution spectroscopy only with the advent of recent instruments like CRIRES on the VLT. The general picture for source n is relatively similar. A single broad absorption feature extending approximately to  km s traces mainly the molecular outflow whereas red-shifted emission likely stems from an inner infalling and accreting envelope/disk. This picture is consistent with disk/outflow proposals for this source deduced from cm and mid-infrared wavelength imaging projects (e.g., Menten & Reid 1995; Greenhill et al. 2004; Shuping et al. 2004). The clearly resolved double-peaked CO structure allows to conduct the Boltzmann analysis for both components. We find that the colder component has higher H column densities ( K &  cm) compared with the second warmer component ( K &  cm). As discussed above for the BN object, differences between the infrared CO derived column densities and those estimated from mm wavelengths observations may arise because both tracers are sensitive to gas (and dust) at different temperatures. In contrast to the BN object, the Pfund  emission from source n is consistent with a thermal Hii region without a strong wind component to the line shape.

Irc3:

The observational signatures from the dust density enhancement IRC3 are very different compared to the two previously discussed sources. The continuum from IRC3 is much weaker, nevertheless we detect CO absorption between and  km s. However, more importantly, we clearly detect extended CO and CO emission over scales of . This extended CO emission shows a multiple wedge-like velocity structure consistent with jet-entrainment models of molecular outflows. Hence IRC3 may be not merely a dust density enhancement, but it may be part of the famous outflow from the Orion-KL region.

Limitations and future:

One shortcoming of the data is that for our primary targets, the BN object and source n, we could not spatially resolve the inner region of the emission as done in the lower-mass case presented in Goto et al. (2006). The reasons may be different for the two sources. Source n is probably still too young and too deeply embedded so that the envelope overwhelms any emission from the embedded disk itself. The BN object has the advantage that its velocity of rest is offset from that of the cloud by more than 10 km , and hence absorption could be less of a problem for such kind of source. However, BN is likely significantly more evolved, and the continuum-to-line ratio is so high that we cannot reasonably filter out the continuum emission. Hence, we cannot well study the inner region of the proposed disk. Nevertheless, observing higher J-transitions as well as the CO overtone emission with todays higher sensitivity compared to the Scoville et al. (1983) observations will likely constrain the proposed disk structure in more detail.

How to proceed now if one wants to do similar-type fundamental CO line studies of disks in high-mass star formation? On the one hand, it is important to not select too young sources because their envelopes will likely almost always “destroy” the emission signatures. There may exist exceptions where one views straight through the outflow cavity face-on toward the disk. On the other hand, for more evolved regions, the continuum emission can be very strong or maybe the remaining disk size can be reduced again making the spatial resolution a problem. Therefore, in addition to very careful target selections, adding AO to achieve the best spatial resolution will be a crucial element for such kind of studies in the coming years. Furthermore, in particular for the mostly saturated CO lines, it will be important to extend the spectral coverage to also observe higher excited CO lines, that will likely not saturate anymore, as well as CO overtone emission. This will allow us to better assess the hotter gas components and hence to conduct a more detailed analysis of the CO data themselves. On longer time-scales, the ELT with its proposed mid-infrared instrument METIS promises orders of magnitude progress in this field based on its superior sensitivity and spatial resolution. With this instrument, we will be truly capable to resolve the gas signatures of accretion disks around (high-mass) YSOs.

Acknowledgements.
We like to thank a lot the anonymous referee as well as the Editor Malcolm Walmsley for thorough reviews which helped improving the paper. H.B. acknowledges financial support by the Emmy-Noether-Program of the Deutsche Forschungsgemeinschaft (DFG, grant BE2578).

References

  • Arce et al. (2007) Arce, H. G., Shepherd, D., Gueth, F., et al. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 245–260
  • Bally (2008) Bally, J. 2008, in Astronomical Society of the Pacific Conference Series, Vol. 387, Massive Star Formation: Observations Confront Theory, ed. H. Beuther, H. Linz, & T. Henning, 158
  • Becklin & Neugebauer (1967) Becklin, E. E. & Neugebauer, G. 1967, ApJ, 147, 799
  • Beltrán et al. (2006) Beltrán, M. T., Cesaroni, R., Codella, C., et al. 2006, Nature, 443, 427
  • Beuther & Nissen (2008) Beuther, H. & Nissen, H. D. 2008, ApJ, 679, L121
  • Beuther et al. (2002a) Beuther, H., Schilke, P., Menten, K. M., et al. 2002a, ApJ, 566, 945
  • Beuther et al. (2005a) Beuther, H., Schilke, P., Menten, K. M., et al. 2005a, ApJ, 633, 535
  • Beuther et al. (2002b) Beuther, H., Schilke, P., Sridharan, T. K., et al. 2002b, A&A, 383, 892
  • Beuther & Walsh (2008) Beuther, H. & Walsh, A. J. 2008, ApJ, 673, L55
  • Beuther et al. (2004) Beuther, H., Zhang, Q., Greenhill, L. J., et al. 2004, ApJ, 616, L31
  • Beuther et al. (2005b) Beuther, H., Zhang, Q., Greenhill, L. J., et al. 2005b, ApJ, 632, 355
  • Blake & Boogert (2004) Blake, G. A. & Boogert, A. C. A. 2004, ApJ, 606, L73
  • Blake et al. (1996) Blake, G. A., Mundy, L. G., Carlstrom, J. E., et al. 1996, ApJ, 472, L49
  • Brittain et al. (2009) Brittain, S. D., Najita, J. R., & Carr, J. S. 2009, ApJ, 702, 85
  • Cesaroni et al. (2007) Cesaroni, R., Galli, D., Lodato, G., Walmsley, C. M., & Zhang, Q. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil, 197–212
  • Cesaroni et al. (2005) Cesaroni, R., Neri, R., Olmi, L., et al. 2005, A&A, 434, 1039
  • Chernin & Wright (1996) Chernin, L. M. & Wright, M. C. H. 1996, ApJ, 467, 676
  • Comito et al. (2005) Comito, C., Schilke, P., Phillips, T. G., et al. 2005, ApJS, 156, 127
  • Dougados et al. (1993) Dougados, C., Lena, P., Ridgway, S. T., Christou, J. C., & Probst, R. G. 1993, ApJ, 406, 112
  • Downes et al. (1981) Downes, D., Genzel, R., Becklin, E. E., & Wynn-Williams, C. G. 1981, ApJ, 244, 869
  • Draine et al. (2007) Draine, B. T., Dale, D. A., Bendo, G., et al. 2007, ApJ, 663, 866
  • Gómez et al. (2005) Gómez, L., Rodríguez, L. F., Loinard, L., et al. 2005, ApJ, 635, 1166
  • Goto et al. (2006) Goto, M., Usuda, T., Dullemond, C. P., et al. 2006, ApJ, 652, 758
  • Greenhill et al. (2004) Greenhill, L. J., Gezari, D. Y., Danchi, W. C., et al. 2004, ApJ, 605, L57
  • Henning et al. (1995) Henning, T., Michel, B., & Stognienko, R. 1995, Planet. Space Sci., 43, 1333
  • Henning et al. (1990) Henning, T., Pfau, W., & Altenhoff, W. J. 1990, A&A, 227, 542
  • Henning et al. (2000) Henning, T., Schreyer, K., Launhardt, R., & Burkert, A. 2000, A&A, 353, 211
  • Hildebrand (1983) Hildebrand, R. H. 1983, QJRAS, 24, 267
  • Jiang et al. (2005) Jiang, Z., Tamura, M., Fukagawa, M., et al. 2005, Nature, 437, 112
  • Käufl et al. (2004) Käufl, H.-U., Ballester, P., Biereichel, P., et al. 2004, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 5492, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. A. F. M. Moorwood & M. Iye, 1218–1227
  • Krumholz et al. (2007) Krumholz, M. R., Klein, R. I., & McKee, C. F. 2007, ApJ, 656, 959
  • Krumholz et al. (2009) Krumholz, M. R., Klein, R. I., McKee, C. F., Offner, S. S. R., & Cunningham, A. J. 2009, Science, 323, 754
  • Lamers & Cassinelli (1999) Lamers, H. J. G. L. M. & Cassinelli, J. P. 1999, Introduction to Stellar Winds (Cambridge University Press, June 1999.)
  • Lawrence et al. (2007) Lawrence, A., Warren, S. J., Almaini, O., et al. 2007, MNRAS, 379, 1599
  • Luhman (2000) Luhman, K. L. 2000, ApJ, 544, 1044
  • Menten & Reid (1995) Menten, K. M. & Reid, M. J. 1995, ApJ, 445, L157
  • Menten et al. (2007) Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515
  • Minchin et al. (1991) Minchin, N. R., Hough, J. H., McCall, A., et al. 1991, MNRAS, 248, 715
  • Nissen et al. (2007) Nissen, H. D., Gustafsson, M., Lemaire, J. L., et al. 2007, A&A, 466, 949
  • Ossenkopf & Henning (1994) Ossenkopf, V. & Henning, T. 1994, A&A, 291, 943
  • Plambeck et al. (1995) Plambeck, R. L., Wright, M. C. H., Mundy, L. G., & Looney, L. W. 1995, ApJ, 455, L189
  • Pontoppidan et al. (2008) Pontoppidan, K. M., Boogert, A. C. A., Fraser, H. J., et al. 2008, ApJ, 678, 1005
  • Reipurth et al. (2007) Reipurth, B., Jewitt, D., & Keil, K., eds. 2007, Protostars and Planets V
  • Rothman et al. (2005) Rothman, L., Jacquemart, D., Barbe, A., & et al. 2005, Journal of Quantitative Spectroscopy and radiative transfer, 96, 139
  • Ryde & Schöier (2001) Ryde, N. & Schöier, F. L. 2001, ApJ, 547, 384
  • Schreyer et al. (2002) Schreyer, K., Henning, T., van der Tak, F. F. S., Boonman, A. M. S., & van Dishoeck, E. F. 2002, A&A, 394, 561
  • Scoville et al. (1983) Scoville, N., Kleinmann, S. G., Hall, D. N. B., & Ridgway, S. T. 1983, ApJ, 275, 201
  • Scoville et al. (1993) Scoville, N. Z., Carlstrom, J. E., Chandler, C. J., et al. 1993, PASP, 105, 1482
  • Sheffer et al. (2007) Sheffer, Y., Rogers, M., Federman, S. R., Lambert, D. L., & Gredel, R. 2007, ApJ, 667, 1002
  • Shuping et al. (2004) Shuping, R. Y., Morris, M., & Bally, J. 2004, AJ, 128, 363
  • Tan (2004) Tan, J. C. 2004, ApJ, 607, L47
  • van der Plas et al. (2009) van der Plas, G., van den Ancker, M. E., Acke, B., et al. 2009, A&A, 500, 1137
  • Yorke & Sonnhalter (2002) Yorke, H. W. & Sonnhalter, C. 2002, ApJ, 569, 846
  • Zapata et al. (2009) Zapata, L. A., Schmid-Burgk, J., Ho, P. T. P., Rodríguez, L. F., & Menten, K. M. 2009, ApJ, 704, L45
  • Zhang et al. (2005) Zhang, Q., Hunter, T. R., Brand, J., et al. 2005, ApJ, 625, 864
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
Cancel
Loading ...
232665
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

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
Test description