# Warm gas in protostellar outflows. II. EHV jet and outflows from OMC-2/3

A. I. Gómez-Ruiz, A. Gusdorf, S. Leurini, K. M. Menten, S. Takahashi, F.

Wyrowski, R. Gusten

► To cite this version:

A. I. Gómez-Ruiz, A. Gusdorf, S. Leurini, K. M. Menten, S. Takahashi, et al. Warm gas in protostellar outflows. II. EHV jet and outflows from OMC-2/3. Astronomy and Astrophysics - A&A, 2019, 10.1051/0004-6361/201424156. [hal-02268368](https://hal.science/hal-02268368)

HAL Id: hal-02268368

https://hal.science/hal-02268368

Submitted on 10 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Warm gas in protostellar outflows

## II. Extremely high-velocity emission jet and outflows from OMC-2/3*

A. I. Gómez-Ruiz¹, A. Gusdorf², S. Leurini³, K. M. Menten⁴, S. Takahashi⁵,⁶,⁷, F. Wyrowski⁴, and R. Güsten⁴

¹ CONACYT-Instituto Nacional de Astrofísica, Óptica y Electrónica, Luis E. Erro 1, C.P. 72840, Tonantzintla, Puebla, México
e-mail: aigomez@inaoep.mx

² LERMA, UMR 8112 du CNRS, Observatoire de Paris, École Normale Supérieure, 61 Av. de l'Observatoire, 75014 Paris, France

³ INAF – Osservatorio Astronomico di Cagliari, Via della Scienza 5, 09047 Selargius (CA), Italy

⁴ Max-Planck-Institut für Radioastronomie (MPIfR), Auf dem Hügel 69, 53121 Bonn, Germany

⁵ Joint ALMA Observatory, Alonso de Córdova 3107, Vitacura, Santiago, Chile

⁶ NAOJ Chile Observatory, Alonso de Córdova 3788, Oficina 61B, Vitacura, Santiago, Chile

⁷ Department of Astronomical Science, School of Physical Sciences, SOKENDAI (The Graduate University for Advanced Studies), Mitaka, Tokyo 181-8588, Japan

Received 8 May 2014 / Accepted 7 June 2019

## ABSTRACT

**Context.** OMC-2/3 is one of the nearest embedded cluster-forming regions that includes intermediate-mass protostars at early stages of evolution. A previous CO (3–2) mapping survey towards this region revealed outflow activity related to sources at different evolutionary phases.

**Aims.** The present work presents a study of the warm gas in the high-velocity emission from several outflows found in CO (3–2) emission by previous observations, determines their physical conditions, and makes a comparison with previous results in low-mass star-forming regions.

**Methods.** We used the CHAMP+ heterodyne array on the APEX telescope to map the CO (6–5) and CO (7–6) emission in the OMC-2 FIR 6 and OMC-3 MMS 1-6 regions, and to observe $^{13}$CO (6–5) at selected positions. We analyzed these data together with previous CO (3–2) observations. In addition, we mapped the SiO (5–4) emission in OMC-2 FIR 6.

**Results.** The CO (6–5) emission was detected in most of the outflow lobes in the mapped regions, while the CO (7–6) was found mostly in the OMC-3 outflows. In the OMC-3 MMS 5 outflow, a previously undetected extremely high-velocity gas was found in CO (6–5). This extremely high-velocity emission arises from the regions close to the central object MMS 5. Radiative transfer models revealed that the high-velocity gas from MMS 5 outflow consists of gas with $n_{H_2} = 10^4–10^5$ cm$^{-3}$ and $T > 200$ K, similar to what is observed in young Class 0 low-mass protostars. For the other outflows, values of $n_{H_2} > 10^4$ cm$^{-3}$ were found.

**Conclusions.** The physical conditions and kinematic properties of the young intermediate-mass outflows presented here are similar to those found in outflows from Class 0 low-mass objects. Due to their excitation requirements, mid-$J$ CO lines are good tracers of extremely high-velocity gas in young outflows likely related to jets.

**Key words.** stars: formation – ISM: clouds – ISM: jets and outflows – shock waves – ISM: individual objects: OMC-2/3

## 1. Introduction

Intermediate-mass (IM) protostars, i.e. protostars whose mass is in the range from 2 to 8 $M_\odot$, provide crucial information on the star formation process because they link the low- and high-mass star formation (e.g. Crimier et al. 2010). IM protostars share some properties of massive protostars, such as their ability to ionize the surrounding gas, and that they exist in clusters. On the other hand, they present the advantage that they are more numerous and therefore are located at closer distances ($\le 1$ kpc) than most of the massive star-forming regions. Therefore, they could fill the gap between low- and high-mass star formation studies. Despite their relevance, very little is known about the formation and early evolutionary stages of IM protostars (e.g. class 0-type objects). Although some detailed studies have been

conducted in recent years (e.g., Beltrán 2015; Takahashi & Ho 2012; Takahashi et al. 2009, 2006; van Kempen et al. 2012; Crimier et al. 2010; Beltrán et al. 2008; Fuente et al. 2007, 2001; Zapata et al. 2004, 2005), little information is available, in particular about outflows from class 0-type sources.

Observations have shown that outflows from IM protostars (hereafter IM outflows) share some characteristics of outflows from low-mass protostars, such as the high collimation factor in early evolutionary phases (Beltrán et al. 2008). On the other hand, IM outflows are, in general, intrinsically more energetic than those driven by low-mass sources (Beltrán et al. 2008; Takahashi et al. 2008). The values of the outflow momentum rate, $F_{\text{out}}$, are higher in IM protostars than in low-mass objects, which may be an indication that IM young stellar objects (YSOs) accrete material faster (at a higher accretion rate) than low-mass ones (Beltrán et al. 2008).

Regarding their physical parameters, such as kinetic temperature ($T_{\text{kin}}$) and $H_2$ density ($n$), only for a few cases have they been determined via observations of more than one CO

<footnote>* Based on observations acquired with the Atacama Pathfinder Experiment (APEX). APEX is a collaboration between the Max-Planck-Institut für Radioastronomie, the European Southern Observatory, and the Onsala Space Observatory.</footnote>
transition (e.g., Cep-E and NGC 2071; Gómez-Ruiz et al. 2012; Chernin & Masson 1992; Hatchell et al. 1999). The best studied case is Cep-E, in which a multi-line CO study at high spectral resolution uncovered different excitation components related to the different kinematic structures (Lefloch et al. 2015). In the jet component, the emission from the low-$J$ CO lines (from $J=1-0$ to $J=5-4$) is dominated by a gas layer with $T_{\text{kin}}=80-100$ K, and $n=(0.5-1) \times 10^5$ cm$^{-3}$; while the high-$J$ CO lines (from $J=12-11$ to $J=16-15$) trace warmer and denser gas, with $T_{\text{kin}}=400-750$ K, and $n=(0.5-1) \times 10^6$ cm$^{-3}$. In the outflow cavity, the low-$J$ CO lines are dominated by a gas layer with $T_{\text{kin}}=55-85$ K, and density in the range $(1-8) \times 10^5$ cm$^{-3}$; while the high-$J$ lines are dominated by a hot, denser gas layer with $T_{\text{kin}}=500-1500$ K, and $n=(1-5) \times 10^6$ cm$^{-3}$. The terminal bowshock consists of gas with moderate excitation, with a temperature in the range $T_{\text{kin}}=400-500$ K, and density $n=(1-2) \times 10^6$ cm$^{-3}$. These observations show a complex excitation structure of the outflow, with the CO (5-4) line (close to what we later call mid-$J$ CO lines) contributing mostly to a relatively low-excitation component from the jet and cavity, and a moderate excitation component from the bowshock.

In a previous paper (Gómez-Ruiz et al. 2013, hereafter Paper I), we presented a study of the CO (6-5) and (7-6) emission (throughout the paper referred to as mid-$J$ CO transitions) in two low-mass class 0 outflows, namely L1448 and HH211. In both objects we found CO (6-5) and (7-6) emission tracing bipolar structures. In the case of L1448, extremely high-velocity emission (EHV), i.e. emission at relative velocities >50 km s$^{-1}$ with respect to the cloud velocity, was detected in both transitions. By performing a large velocity gradient (LVG) analysis, we inferred that the gas related with the outflow structures is dense (>10$^5$ cm$^{-3}$) and warm ($T > 200$ K). It was determined that the EHV emission shows higher densities and temperatures than the high-velocity emission. Based on the LVG results and the position-velocity distribution of the emitting gas, we concluded that these mid-$J$ CO transitions are good tracers of molecular material likely related to the primary jet and to the highly excited bow-shock region in protostellar outflows.

In the present paper we extend the study of mid-$J$ CO emission to outflows from IM protostars in the Orion molecular cloud-2/3 (OMC-2/3), including Class 0- and I-type objects. The main purpose of this paper is to present the study of the distribution of the warm gas traced by the CO (6-5) and (7-6) transitions in IM outflows and compare it with the low-mass case. By incorporating lower-$J$ CO transitions in our analysis, the excitation and physical conditions of IM outflows, as traced by these lines, were also studied.

## 2. Target regions

OMC-2/3 is located north of the Orion-A Giant Molecular Cloud (OMC 1). The distance to OMC 1 range from ~390 pc to 440 pc, given its 3-D structure determined by recent astrometry data (Kounkel et al. 2018). This region is one of the nearest cluster-forming sites with IM protostars (e.g. Takahashi et al. 2006, 2008, here after T08). In this region, filamentary molecular clouds were discovered in the millimeter/submillimeter continuum observations with the IRAM-30m telescope and the JCMT (Chini et al. 1997; Johnstone & Bally 1999). Within the filaments, millimeter/submillimeter continuum sources are found, presumably tracing prestellar and protostellar cores (e.g., Takahashi et al. 2009, 2012, 2013; Takahashi & Ho 2012; Chini et al. 1997). Follow-up CO and H$_2$ molecular line observations revealed the presence of a number of outflows originating from

<figure><img src="image_2.png" /><figcaption>**Fig. 1.** CO (6-5) outflow emission in the OMC-2 FIR 6 region. Blue contours are for the blue-shifted emission in the range from -12 to +8 km s$^{-1}$; while red contours are for the red-shifted emission from +14 to +30 km s$^{-1}$. The first contour is $3\sigma$, with contour spacing of $2\sigma$ ($\sigma=5.5$ and $2.4$ K km s$^{-1}$, for blue- and red-shifted emission, respectively). The black solid and dashed contours are the blue- and red-shifted CO (3-2) emission, respectively, as presented by T08. The black square shows the whole field covered by our CHAMP+ observations. The triangles show the positions of the continuum sources and the green dashed lines indicate the outflow orientation (both according to T08).</figcaption></figure>

some of the continuum sources (T08; Yu et al. 1997; Aso et al. 2000; Stanke et al. 2002; Williams et al. 2003). One of the most extensive molecular outflow surveys in the OMC-2/3 region was conducted by T08. By observing the CO (3-2) transition with the ASTE telescope, these authors detected 14 molecular outflows, indicative of a high ongoing star formation activity in this region.

From the CO (3-2) outflow survey of T08 and its follow-up observations, we selected the most promising sources in terms of outflow activity (bipolar morphology and high-velocity emission) and early evolutionary stage (class 0-type sources). Under these criteria, the regions around FIR (far-infrared) 6 in OMC-2 and MMS (millimeter source) 1-6 in OMC-3 (hereafter OMC-2 FIR 6 and OMC-3 MMS 1-6, respectively) were selected for follow-up observations in mid-$J$ CO lines. For these sources the derived luminosities are in the range of 6-251 and 9-150 $L_\odot$, respectively, which makes them very likely IM protostars (T08, Takahashi et al. 2009; Furlan et al. 2016). The envelope masses determined by T08 (>3 $M_\odot$) are consistent with values of IM protostars.

OMC-2 FIR 6 contains at least four millimeter continuum sources named FIR 6a, 6b, 6c, and 6d (see Fig. 1). T08 reported the FIR 6b source as being a likely Class I-type object (based on their SED), while the lack of sufficient information for the
<figure><img src="image_3.png" /></figure>

**Fig. 2.** CO (6–5) outflow emission in the OMC-3 MMS 1–6 region. Blue contours are for the blue-shifted emission in the range from –25 to +7 km s⁻¹; while red contours are for the red-shifted emission from +16 to +31 km s⁻¹. The first contour is 3σ and contour spacing is 2σ (σ = 1.9 and 0.8 K km s⁻¹, for blue- and red-shifted emission, respectively). The black solid and dashed contours are the blue- and red-shifted CO (3–2) emission, respectively, as presented by T08. The black square shows the whole field covered by our CHAMP+ observations. The triangles show the positions of the continuum sources and the green dashed lines indicate the outflow orientation (both according to T08).

rest of the sources prevents us from clearly assessing their evolutionary status. The CO (3–2) data towards OMC-2 FIR 6 region revealed that two molecular outflows are driven by FIR 6b and 6c (see Fig. 1). As a follow-up study of T08, *Shimajiri et al.* (2009) observed the OMC-2 FIR 6 region with the Nobeyama Millimeter Array in the 3 mm continuum, CO (1–0), and SiO (2–1) transitions. The two outflows found in CO (3–2) were confirmed in CO (1–0), in addition to a third outflow emanating from FIR 6d. The SiO (2–1) emission was found only at the blue-lobe of FIR 6c outflow, but no SiO emission was clearly detected towards the 6b outflow. From the line ratio between the CO (3–2) and (1–0) transitions, it was found that kinetic temperature in the outflow lobes is ~50 K (T08; *Shimajiri et al.* 2009).

The OMC-3 MMS 1-6 region consists of six millimeter continuum sources, named by Chini et al. (1997) as MMS (see Fig. 2). A clear identification of the evolutionary status through the SED was given for MMS 2 and MMS 5, which are class I- and class 0-type, respectively (T08; *Takahashi et al.* 2009). The MMS 2 source is actually a binary system of class I sources (T08; *Tsujimoto et al.* 2004). T08 reported clear bipolar outflow morphologies originating from MMS 2 and MMS 5 in the CO (3-2) emission. In the case of MMS 6, subarcsec resolution observations with the SMA in the CO (3–2) line by *Takahashi & Ho* (2012) revealed an extremely compact molecular outflow (lobe size of ≈800 AU) associated with the continuum peak named MMS 6-main. Furthermore, *Takahashi et al.* (2012) reported that MMS 6-main is most likely in the protostellar phase. No clear molecular outflow signatures are detected in the other MMS sources, so that they are likely in the prestellar phase (*Takahashi et al.* 2013).

### 3. APEX observations

In this section we describe the observations made with two instruments of the APEX observatory.

#### 3.1. CHAMP+ observations

Submillimeter observations towards OMC-3 MMS 1-6 and OMC-2 FIR 6 were performed with the MPIfR principal investigator (PI) instrument CHAMP⁺ (The Carbon Heterodyne Array of the MPIfR, *Kasemann et al.* 2006) on the APEX telescope (Atacama Pathfinder EXperiment, *Güsten et al.* 2006) during November 2008 and July 2010, respectively. CHAMP⁺ is a dual-color 2×7 pixel heterodyne array for operation in the 450 µm (low-frequency array, LFA) and 350 µm (high-frequency array, HFA) atmospheric windows. The optics allow simultaneous observations in both colors. Both subarrays present a hexagonal arrangement. The front-end was connected to a fast Fourier transform spectrometer (FFTS, *Klein et al.* 2006).

The CO (6–5) and (7–6) lines were mapped simultaneously in the on-the-fly (OTF) mode, with an ON time of 0.5 second per position, with steps of 3". For the OMC-2 FIR 6 region a map of 60"×120" was centered at the position of the FIR 6b source. For the OMC-3 MMS 1-6 region a map with of 280"×150" was centered at the position of MMS 5. Table 1 shows the coordinates (right ascension and declination) of the central position of the maps and their average rms. Additional single-pointing observations were done on two selected positions of OMC-2 FIR 6 (called FIR 6b and FIR 6c-B1; see Sect. 4) in a setup consisting of ¹³CO (6–5) and HCO⁺ (9–8). No emission was detected in the HCO⁺ (9–8) line, with upper limits of 186 mK ($T_A^*$), and therefore was not further discussed. For all observations, the backend was set to provide a total bandwidth of 2.8 GHz divided in 8192 channels for each pixel. The final spectra were re-sampled to 1.0 km s⁻¹ spectral resolution for all transitions in both sources.

During the observations, the precipitable water vapor (PWV) was in the range of 0.3–0.7 mm. The pointing was determined by CO (6–5) cross-scan observations on IK Tau. Pointing accuracy was always within 3". The focus was checked on Jupiter and Mars. The calibration was done by observing hot and cold
<table><caption>**Table 1.** Map central position and rms.</caption><thead><tr><th rowspan="2">Source</th><th rowspan="2">α(J2000)<br>(h m s)</th><th rowspan="2">δ(J2000)<br>(° ′′)</th><th colspan="2">rms. (K) <sup>(a)</sup></th></tr><tr><th><sup>12</sup>CO<br>(6–5)</th><th>SiO<br>(7–6)</th></tr></thead><tbody><tr><td>OMC-2 FIR 6b</td><td>05:35:23.4</td><td>-05:12:03</td><td>0.7</td><td>3.0</td><td>0.05</td></tr><tr><td>OMC-3 MMS 5</td><td>05:35:22.5</td><td>-05:01:15</td><td>0.5</td><td>1.9</td><td>—</td></tr></tbody></table>

**Notes.** <sup>(a)</sup>Average rms ($T_{MB}$) in a channel map of 1.0 km s⁻<fnref n="1" />.

<table><caption>**Table 2.** Characteristics of the observed transitions.</caption><thead><tr><th>Line</th><th>ν<sub>0</sub><br>GHz</th><th>A<br>10<sup>-5</sup> s<sup>-1</sup></th><th>E<sub>u</sub>/k<br>K</th><th>HPBW<br>"</th></tr></thead><tbody><tr><td>CO (6–5)</td><td>691.473</td><td>2.137</td><td>116.16</td><td>9.0</td></tr><tr><td>CO (7–6)</td><td>806.651</td><td>3.422</td><td>154.87</td><td>7.7</td></tr><tr><td><sup>13</sup>CO (6–5)</td><td>661.067</td><td>299.0</td><td>111.05</td><td>9.4</td></tr><tr><td>SiO (5–4)</td><td>217.107</td><td>51.96</td><td>31.26</td><td>30.5</td></tr><tr><td>HCO<sup>+</sup> (9–8)</td><td>802.458</td><td>4400</td><td>192.58</td><td>7.8</td></tr></tbody></table>

**Notes.** Molecular data taken from LAMDA (Schöier et al. 2005).

loads. In our observations the OFF position was not checked for emission, but this position was selected to avoid regions known to present emission in low–J CO transitions. However, due to an erroneous selection of the OFF position during the <sup>13</sup>CO(6–5) and HCO<sup>+</sup>(9–8) single-pointing observations of the FIR 6b position, contamination of the OFF position was found in these spectra. Standard data reduction procedures such as flagging bad spectra and base line subtraction were done with CLASS, while images were produced with the XY_MAP task in GReG, both programs being part of the GILDAS software<fnref n="1" />. The gridding procedure in XY_map was made to provide a final resolution equal to the half power beamwidth (HPBW) of the telescope at the observed frequency. Table 2 lists the observed transitions and rest frequencies ($v_0$), their corresponding HPBW, upper state energy ($E_u$), and Einstein coefficient (A).

The images and spectra are presented in main beam brightness temperature scale, $T_{MB} = T^*/\eta_{MB}$ ($\eta_{MB} = \eta_s/F_{eff}$), for which we have used a forward efficiency ($F_{eff}$) of 0.95 and a beam coupling efficiency ($\eta_s$) measured on Jupiter. During the OMC-2 observations the measured $\eta_s$ was 0.48 for both LFA and HFA, while during the OMC-3 observations $\eta_s$ was 0.48 and 0.45 for LFA and HFA, respectively. Based on the variation of the integrated mid–J CO flux observed in the line pointing sources, we estimated a 20% calibration uncertainty for both LFA and HFA (see also Yıldız et al. 2015).

## 3.2. APEX-1 observations

Complementary mapping observations in the SiO (5–4) line were carried out towards OMC-2 FIR 6b with the APEX-1 facility instrument in September 2010. The receiver was connected to a FFTS that provides a bandwidth of 1 GHz and 4096 channels, which at the rest frequency of SiO (5–4) resulted in a spectral resolution of 0.3 km s⁻<fnref n="1" />. The OTF map with a size of 70" × 70", in steps of 9", was centered at the position of the FIR 6b source. Therefore, the map covers the central parts of the FIR 6b

<table><caption>**Table 3.** Coordinates of the outflow peak positions.</caption><thead><tr><th>Position</th><th>α(J2000)<br>(h m s)</th><th>δ(J2000)<br>(° ′′)</th><th>v<sub>min</sub>,v<sub>max</sub><br>(km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>FIR 6b-B1</td><td>05:35:21.2</td><td>-05:12:15</td><td>5,8</td></tr><tr><td>FIR 6b-R1</td><td>05:35:23.7</td><td>-05:11:59</td><td>14,30</td></tr><tr><td>FIR 6c-B1</td><td>05:35:21.6</td><td>-05:13:00</td><td>-16,8</td></tr><tr><td>FIR 6c-B2</td><td>05:35:21.7</td><td>-05:13:00</td><td>-10,8</td></tr><tr><td>MMS 5-B1</td><td>05:35:21.3</td><td>-05:01:16</td><td>-95,7</td></tr><tr><td>MMS 5-R1</td><td>05:35:22.5</td><td>-05:01:14</td><td>16,75</td></tr><tr><td>MMS 2-B1</td><td>05:35:18.2</td><td>-05:00:32</td><td>-12,7</td></tr><tr><td>MMS 2-R1</td><td>05:35:17.7</td><td>-05:00:34</td><td>16,30</td></tr></tbody></table>

**Notes.** <sup>(*)</sup>Velocity boundaries of the outflow emission (see Sect. 4).

outflow lobes and the blue lobe of the FIR 6c outflow. The pointing accuracy was checked by CO (2–1) cross-scan observations on IK Tau. Pointing accuracy was always within 5". The focus was checked on Jupiter and Mars. The calibration was done by observing hot and cold loads. The data reduction and image processing are similar to the procedure described for the CHAMP+ observations, using CLASS and GReG. The final spectra were re-sampled to a 2.0 km s⁻<fnref n="1" /> spectral resolution. In order to convert to $T_{MB}$ scale, the nominal values of 0.95 and 0.75 were adopted for the forward and beam coupling efficiency, respectively. A calibration uncertainty of 20% is also expected for this receiver, based on the flux variation of the line pointing sources.

## 4. Results

In order to define the outflow velocity range we use the CO (6–5) channel maps and the spectra taken at different outflow positions. We defined the lower boundary limit ($v_{min}$) by looking for outflow morphologies clearly separated from the ambient diffuse emission. We then established $v_{max}$ from the spectra taken at the peak positions within the outflow’s lobes. In our case, $v_{max}$ was defined as the velocity in which the emission in the spectrum drops below the 3σ level. The velocity limits for all identified outflows and the positions of the outflow peaks are given in Table 3. Based on the measurements of the H<sup>13</sup>CO<sup>+</sup> (1–0) line by Aso et al. (2000), throughout the paper we assume the cloud LSR velocity ($V_{LSR}$) of +11 km s⁻<fnref n="1" />.

### 4.1. OMC-2 FIR 6 outflows

In the following we describe the results obtained for the CO (6–5) and CO (7–6), <sup>13</sup>CO<sup>+</sup> (6–5), and SiO (5–4), respectively.

#### 4.1.1. CO (6–5) and CO (7–6) emission

Figure 1 shows the CO (6–5) emission detected in the FIR 6 region, integrated along the blue- and red-shifted wings. Figure 1 also shows the CO (3–2) wing emission from T08. We found the two outflows related to FIR 6b and FIR 6c, reported previously in CO (3–2) by T08. Both outflow lobes from FIR 6b are detected in CO (6–5), and also the blue-shifted lobe from FIR 6c (the red-shifted lobe was not covered by our map). Four emission peaks were found at the outflow positions, which we list in Table 3. The CO (6–5) emission from FIR6b seems to peak closer to the central source than the CO (3–2) emission. Convolutions to the CO (3–2) beam confirms that this tendency is not an effect of the different angular resolutions. On the other hand, the peak position

<footnote n="1">¹ http://iram.fr/IRAMFR/PDB/gildas/gildas.html</footnote>
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** CO (6–5) sample spectra at selected positions (see Table 3) within FIR 6 outflows (1 km s⁻¹ spectral resolution). The dashed line indicates the cloud velocity.</figcaption></figure>

of the CO (6–5) emission in the blue lobe of FIR 6c outflow is consistent with the location of the CO (3–2) emission peak. No high-velocity emission tracing outflows from FIR 6a and FIR 6d was found. In the case of FIR 6d we are not able to unambiguously probe the outflow emission, since the outflow structures revealed by previous low-*J* CO observations are located right at the edge of our map. Around the northern edge of the map, blue- and red-shifted emission is also noticed; however, these structures are not classified as a clear outflow by T08, and therefore we exclude these structures in the analysis that follows.

Figure 3 shows the spectra taken at the position of the emission peak of each lobe. In all positions, the peak of the CO (6–5) profiles is found at the cloud systemic velocity ($V_{LSR} \approx 11.0$ km s⁻¹). Three positions show wings (FIR 6b-R1) and secondary peaks at high velocities (FIR 6c-B1 and FIR 6c-B2). At the FIR 6b-B1 position the line profile only extends down to $V_{LSR} = 5$ km s⁻¹.

The CO (7–6) emission was detected mostly at low velocities ($V_{LSR} = 8–14$ km s⁻¹), tracing the cloud and/or envelope emission. At the outflow peaks the CO (7–6) emission is barely noticed ($S/N < 5$) in the spectra.

### 4.1.2. ¹³CO (6–5) emission and opacity estimation

In Fig. 4 we show the ¹³CO (6–5) spectra taken at two positions within the FIR 6 region: the FIR 6b source and FIR 6c-B1. The ¹³CO (6–5) emission is only detected in about three channels, around the cloud velocity. All spectra in Fig. 4 show the

emission peaking at the systemic velocity ($V_{LSR} = 11$ km s⁻¹), with a FWHM of ~2 km s⁻¹. At the FIR 6b position an absorption feature is also present at $V_{LSR} = 9$ km s⁻¹, which is more evident in the external pixels. This absorption feature is likely due to contamination from the OFF position (see Sect. 2), and thus the ¹³CO (6–5) emission detected at the FIR 6b position should be taken as a lower limit.

Based on the ¹³CO (6–5) observations, we estimate the ¹²CO (6–5) line opacity at the FIR 6b and FIR 6c-B1 positions. We follow the formulation introduced in Paper I, and use the assumption that the ¹³CO emission is optically thin and that ¹²CO and ¹³CO have the same excitation temperature, which means that the ¹²CO to ¹³CO antenna temperature line ratio provide a measure of the ¹²CO opacity, $\tau_{12}$ (see Eq. (1) of Paper I). Our opacity calculations were made assuming a [¹²CO]/[¹³CO] abundance ratio of 50 (Kahane et al. 2018). We also note that the error of the ratios does not require additional calibration uncertainties, since the observations were made with the same instrument and within the same band. In Table 4 we show the opacities estimated for the FIR6b and FIR6c-B1 positions. These estimations are performed using the three channels around the cloud velocity. The $\tau_{12}$ values vary from ~2 to 10. Therefore we see that the ¹²CO is optically thick at the cloud velocity. On the other hand, in the outflow velocity range, where ¹³CO was not detected, the $3\sigma$ upper limits indicate that the emission can be optically thin: with an intensity ratio >50, $\tau_{12}$ is <0.4. Based on these results, in the calculations presented in the analysis, to determine the gas excitation and outflow properties, we assume that the ¹²CO (6–5) high-velocity emission is optically thin in the FIR 6 outflows, and in this way use the simple formulations for the optically thin case. Although we did not observe the ¹³CO (6–5) line in MMS 1–6 region, we also assume that the ¹²CO (6–5) high-velocity emission from the outflows in that region is optically thin.

### 4.1.3. SiO (5–4) emission

The SiO (5–4) emission is only detected to the SW of FIR 6b, in the velocity range from −12 to 15 km s⁻¹. The upper panel of Fig. 5 shows the SiO (5–4) integrated emission map. The SiO emission is located towards the region that corresponds to the blue-lobe of the FIR 6c outflow. The angular resolution (~30′′) of the SiO (5–4) observations is not sufficient to allow a clear identification of outflow structures. The lower panel of Fig. 5 shows the SiO (5–4) spectrum averaged over the area in which the emission is detected (indicated by the green polygon in the upper panel of Fig. 5). The SiO (5–4) averaged spectrum shows a velocity range −12 ≤ $V_{LSR}$ ≤ +15 km s⁻¹. The lower panel of Fig. 5 also shows a comparison with the CO (6–5) spectrum averaged over the same region, obtained from a map convolved to the same angular resolution. The velocity range is similar, suggesting that the SiO (5–4) and CO (6–5) are tracing the same outflow gas. However, while the CO profile shows a typical wing-like profile with a main peak at low velocities, the SiO is not wing-like and with a secondary peak at high velocities ($V_{LSR} \sim -8$ km s⁻¹). This behavior, also observed in other outflows, is indicative of the uniqueness of the SiO as a shock tracer (e.g., L1448: Paper I; Bachiller et al. 1990). In CO (6–5) an indication of secondary peaks at high velocities is only observed in the spectra at the original angular resolution (Fig. 3), in fact at about the same SiO (5–4) velocity range.

Our APEX SiO (5–4) data showed a SiO gas distribution similar to that found in previous interferometric SiO (2–1) observations by Shimajiri et al. (2009). In both cases the SiO emission was detected only towards the region that corresponds to the blue
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** $^{13}$CO (6–5) spectra taken at selected positions within FIR 6 region (1 km s$^{-1}$ spectral resolution). The figure shows the seven CHAMP+ pixels for each pointing observation. Left: observation centered on FIR 6b. Right: observation centered on FIR 6c-B1 position (offsets with respect to FIR 6c).</figcaption></figure>

<table><caption>**Table 4.** $^{12}$CO (6–5) opacity around cloud $V_{LSR}$ in FIR 6b and FIR 6c-B1.</caption><thead><tr><th>$V_{LSR}$ (km s<sup>-1</sup>)</th><th><sup>12</sup>CO/<sup>13</sup>CO (6–5) T<sub>MB</sub> ratio</th><th>Opacities</th></tr></thead><tbody><tr><td colspan="3" style="text-align: center;">FIR 6b</td></tr><tr><td>10</td><td>8.6(0.3)</td><td>5.8(0.3)</td></tr><tr><td>11</td><td>4.8(0.1)</td><td>10.4(0.2)</td></tr><tr><td>12</td><td>7.3(0.3)</td><td>6.8(0.3)</td></tr><tr><td colspan="3" style="text-align: center;">FIR 6c-B1</td></tr><tr><td>10</td><td>10.1(0.5)</td><td>4.9(0.2)</td></tr><tr><td>11</td><td>5.5(0.2)</td><td>9.1(0.3)</td></tr><tr><td>12</td><td>18.2(2.0)</td><td>2.5(0.3)</td></tr></tbody></table>

**Notes.** Parenthesis show statistical errors.

lobe of the FIR 6c outflow. Our larger map confirms that no SiO emission is clearly detected from the FIR 6b outflow. Similar cases where CO outflows do not present an SiO counterpart have been reported in the literature (see Zapata et al. 2006, for discussion).

## 4.2. OMC-3 MMS 1-6 outflows: CO (6–5) and CO (7–6) emission

Figure 2 shows the region covered by our CO (6–5) and (7–6) maps. We clearly identified the CO (6–5) outflow emission originating from the sources MMS 5 and MMS 2. We identified four outflow peaks, listed in Table 3. These two outflows were previously reported by T08. However the compact outflow detected in CO (3–2) by Takahashi & Ho (2012) emanating from MMS 6 was not detected in CO (6–5). The nondetection of this outflow is possibly due to beam dilution, since even with Nobeyama interferometric CO (1–0) observations by Takahashi et al. (2009), with an angular resolution of several arcseconds, the compact outflow was not detected (see also discussion by Takahashi & Ho 2012). No clear evidence of outflow activity towards the remaining MMS objects was found. The outflow emanating from MMS 5 seems to be orientated along the east-west direction,

while for the outflow from MMS 2 the observations cannot clearly disentangle a tendency in the outflow direction. We also notice two blobs of blue-shifted emission on the far side of the putative MMS 2 blue-shifted outflow lobe, which coincide with similar structures detected in CO (3–2). Since T08 and Williams et al. (2003) could not unambiguously relate these structures with MMS2 (or the other nearby millimeter sources), for consistency with the later physical parameters comparisons, we did not take these structures into account in the following calculations. The most intense CO (6–5) emission is related to the MMS 5 outflow, which shows an extended blue lobe, while the red lobe shows a more compact structure.

A remarkable result from our CO (6–5) observations is the detection of EHV emission from the MMS 5 outflow. The EHV range is defined as $V_{LSR} = -95$ to $-80$ km s$^{-1}$ for the blue-shifted lobe and $V_{LSR} = 60–70$ km s$^{-1}$ for the red-shifted lobe. In Fig. 6 we show the velocity integrated emission in the EHV range, overlaid on the near-infrared (NIR) emission (K-band, 2.12 µm, taken from T08), and the spectra taken at the position of the peak intensity of the EHV structures. The integrated intensity maps show that the EHV emission is located within a radius of ~10" from MMS 5. The blue-shifted EHV emission coincides with the inner collimated NIR structure. The CO (6–5) spectra clearly show secondary peaks related to these EHV features, with the red-shifted EHV component stronger than the blue-shifted one. The compact morphology and the characteristic line profiles suggest that these components are the equivalent of the so-called EHV bullets found in class 0 low-mass outflows (e.g. L1448-mm: Bachiller et al. 1990, Paper I). Although evidence for EHV emission was not found in the CO (3–2) spectra observed by T08, recent ALMA observations of CO (2–1) and SiO (5–4) have revealed the EHV emission from MMS 5 (Matsushita et al. 2019).

The CO (7–6) high-velocity emission was detected in the outflows associated with MMS 5 and MMS 2 (see Fig. 7). The outflow lobes from MMS 5 in the CO (7–6) emission show a similar structure to those traced in the CO (6–5) emission. An EHV component in the CO (7–6) line is barely noticed at $3\sigma$ only at the R1 position. In the case of the MMS 2 outflow, the CO (7–6) high-velocity emission is coming mainly from the red-shifted lobe, with the blue-shifted lobe only detected at the $4\sigma$ level.
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** *Upper panel:* SiO (5–4) map of the FIR 6 region, integrated from −12 to 15 km s⁻¹ (black contours), overlaid with the CO (6–5) outflow emission (blue and red contours; as in Fig. 1). The SiO contour step is 3σ, with the first contour at 3σ (σ = 0.38 K km s⁻¹). *Lower panel:* SiO (5–4) spectrum (black) averaged over the region indicated by the green polygon on the top panel. To increase the signal-to-noise ratio the spectral resolution was smoothed to 2 km s⁻¹. As comparison the CO (6–5) spectrum is shown in red.</figcaption></figure>

# 5. Analysis

I present the results of a radiative transfer analysis using two methods, as well as the determination of the outflow parameters.

## 5.1. Population diagram

Under the assumption that the gas is in LTE and that the ¹²CO emission is optically thin, we use the Boltzmann relation for the level populations to estimate the rotational temperature ($T_{\text{rot}}$) and the total column density ($N$) at selected positions within the outflow lobes, following the formulation presented in Paper I. The calculations were made for the emission from the outflow peaks reported in Table 3 (see also Figs. 1 and 2). The CO (3–2) data taken with ASTE by T08 were added to our APEX data. We convolved our CO (6–5) and (7–6) observations to the angular resolution of the CO (3–2) observations (i.e., 26′′). Table 5 shows the CO (3–2, 6–5, 7–6) integrated intensities of the selected positions, convolved to 26′′ resolution. The emission

is integrated between the outflow’s wing limits, $v_{\text{min}}$ and $v_{\text{max}}$, which we assume to be optically thin (see Sect. 4). In the calculation we included the calibration uncertainties of each line, whose value is 20% for all the cases (including the ASTE data; see T08). For future references, we report separately in Table 5 the intensities for the EHV range, although we do not provide calculations for this component, given that the CO (3–2) and CO (7–6) lines were mostly not detected.

Table 6 shows the results, $T_{\text{rot}}$ and $N$, obtained from the rotational diagram analysis. To within the errors, most of the positions present a rotational temperature in the range of 60–70 K, with the exception of the FIR 6b-B1 position, which shows a lower temperature of ~30 K. In terms of the column density, higher values are found in the OMC-2 FIR 6 positions.

## 5.2. Large velocity gradient

Radiative transfer calculations with RADEX (van der Tak et al. 2007) were made in the LVG approximation and plane-parallel geometry. The molecular data was retrieved from the LAMDA database<fnref n="2" /> (Schöier et al. 2005) and the collisional rate coefficients were adopted from Yang et al. (2010).

The RADEX calculations were used to reproduce simultaneously the observed CO (3–2)/CO (6–5), CO (3–2)/CO (7–6), and CO (6–5)/CO (7–6) integrated intensity ratios at different positions within the outflow lobes. The intensities are from Table 5, and we added the calibration uncertainties already reported; we note, however, that the ratios with the CHAMP+ lines themselves do not require additional calibration errors since they were observed simultaneously.

In Fig. 8, kinetic temperature ($T_{\text{kin}}$) vs. H₂ density ($n$) plots are shown for each outflow position. In most of the cases only lower limits to kinetic temperature and H₂ density can be provided, with values of $T_{\text{kin}}$ > 30–50 K and $n$ > 10⁴ cm⁻³. Only in the case of MMS 5-R1 the three ratios are able to set more stringent constraints to the density and a higher kinetic temperature lower limit: $n \approx 10^4–10^5$ cm⁻³, $T_{\text{kin}}$ > 200 K. These values indicate that the outflow driven by MMS 5 present physical conditions similar to those obtained in the low-mass class 0 outflows from L1448-mm and HH211 (see Paper I).

## 5.3. Outflow parameters

Using the velocity integrated CO (6–5) high-velocity maps (pixels with S/N above 3), we calculate the H₂ mass of the outflow gas, following the relations presented in Paper I between the total CO column density ($N$) and the velocity integrated emission, under the assumption of optically thin ¹²CO (shown to be the case, based on the opacity upper limit determined in Sect. 4 for the high-velocity gas at some outflow positions). The high-velocity maps are produced by integrating the emission between the minimum and maximum velocity of each outflow lobe, following the $v_{\text{min}}$,$v_{\text{max}}$ values reported in Table 3. In the calculations we used the temperatures obtained from the rotational diagram analysis (Table 6), the standard ISM value for the relative ¹²CO abundance, X(H₂/CO) ~ 10⁴ (e.g., Wilson & Matteucci 1992), and a distance to the source of 400 pc (which is close to the average of the distance range of Orion A 3D structure). As in Paper I, we point out that considering uncertainties such as data calibration, definition of the velocity boundary limits for the integrated emission, and opacity assumptions, the mass determinations are usually accurate within a factor of 2–3 (see, e.g., Cabrit & Bertout 1990).

<footnote n="2">² http://www.strw.leidenuniv.nl/~moldata/</footnote>
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Left: EHV CO (6–5) emission from MMS 5 outflow. The solid blue contours show the CO (6–5) emission integrated from $V_{LSR} = -91$ to $-78$ km s$^{-1}$ (blue-shifted EHV range), while the solid red contours show the CO (6–5) emission integrated from $V_{LSR} = 60–75$ km s$^{-1}$ (red-shifted EHV range). The first contour is $3\sigma$, with contour spacing in steps of $\sigma$ (2.2 K km s$^{-1}$). The dashed contours show the high-velocity CO (6–5) emission as presented in Fig. 2. The background greyscale image is the NIR (K-band, 2.12 µm) map taken from T08. Right: CO (6–5) spectrum taken at the B1 and R1 positions (see Table 3). The vertical dot-dashed lines show the boundaries of the EHV range, while the dashed line indicates the cloud velocity.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 7.** CO (7–6) high-velocity emission map centered on MMS 5. Blue and red contours show the blue- and red-shifted high-velocity emission (integrated in the same velocity ranges as in Fig. 2), respectively. Contour spacing is $1\sigma$, with $3\sigma$ as the first contour ($\sigma = 7.4$ K km s$^{-1}$ for blue-shifted emission and $\sigma = 6.7$ K km s$^{-1}$ for red-shifted emission). The background greyscale image is as in Fig. 4, the filled triangles as in Fig. 2. The inlay shows the CO (7–6) spectrum taken at the B1 and R1 positions, with the vertical lines as in Fig. 6.</figcaption></figure>

The mass detection limit of our CO (6–5) observations, assuming a $T_{ex}$ of 60 K (following Table 6) and the rms at 1 km s$^{-1}$, is $\sim 1 \times 10^{-5} M_{\odot}$ inside a beam of $9''$. For each outflow the $H_2$ mass is reported in Table 7. The total $H_2$ masses range from $\sim 1 \times 10^{-2} M_{\odot}$ to $1 \times 10^{-3} M_{\odot}$. We note that our mass estimates based on the CO (6–5) emission are lower by a factor of 4 with respect to those reported in Table 3 of T08 for MMS 5, MMS 2, and FIR 6c (blue-lobe), while it is lower by a factor of 9 for FIR 6b outflow. The higher discrepancy between the mass estimates from CO (3–2) and CO (6–5) in FIR 6b outflow may be mostly due to the fact that our CO (6–5) map misses the upper

<table><caption>**Table 5.** CO integrated intensities at selected positions.</caption><thead><tr><th rowspan="2">Position</th><th rowspan="2">Velocity (km s<sup>-1</sup>)<br>(v<sub>min</sub>, v<sub>max</sub>)</th><th colspan="3">∫ T<sub>MB</sub> dv(σ) (K km s<sup>-1</sup>) <sup>(a)</sup></th></tr><tr><th>3–2</th><th>6–5</th><th>7–6</th></tr></thead><tbody><tr><td>FIR 6b-B1</td><td>0, 8</td><td>32(2)</td><td>17(1)</td><td>&lt;4 <sup>(b)</sup></td></tr><tr><td>FIR 6b-R1</td><td>14, 30</td><td>22(1)</td><td>27(1)</td><td>20(4)</td></tr><tr><td>FIR 6c-B1</td><td>-16, 8</td><td>23(1)</td><td>27(2)</td><td>20(5)</td></tr><tr><td>FIR 6c-B2</td><td>-10, 8</td><td>33(1)</td><td>39(1)</td><td>20(5)</td></tr><tr><td>MMS 5-B1</td><td>-25, 7</td><td>36(2)</td><td>61(1)</td><td>42(4)</td></tr><tr><td>EHV</td><td>-95, -75</td><td>&lt;10 <sup>(b)</sup></td><td>11(1)</td><td>3(3)</td></tr><tr><td>MMS 5-R1</td><td>16, 31</td><td>6(2)</td><td>16(1)</td><td>8(2)</td></tr><tr><td>EHV</td><td>60, 75</td><td>&lt;7 <sup>(b)</sup></td><td>14(1)</td><td>8(2)</td></tr><tr><td>MMS 2-B1</td><td>-12, 7</td><td>16(2)</td><td>18(1)</td><td>12(3)</td></tr><tr><td>MMS 2-R1</td><td>16, 30</td><td>9(1)</td><td>11(1)</td><td>8(3)</td></tr><tr><td></td><td>31, 50</td><td>&lt;8 <sup>(b)</sup></td><td>7(1)</td><td>&lt;3 <sup>(b)</sup></td></tr></tbody></table>

**Notes.** Data convolved to $26''$ ASTE beam (i.e. CO (3–2) transition). <sup>(a)</sup>Integrated between $v_{min}, v_{max}$. Statistical errors in parenthesis. <sup>(b)</sup>$3\sigma$ upper limit.

half of the red-shifted lobe, a region that was included in the CO (3–2) mass estimation from T08. By recomputing the mass from the CO (3–2) emission in an area equivalent to our CO (6–5) map, the difference is reduced to a factor of 3. The lower masses obtained with respect to the CO (3–2) observations may be due the typical mass estimation uncertainties we already mention above; however, by looking at the maps it is clear in some cases (such as the MMS 2 outflow) that the emitting region in both lines is not exactly the same, with the CO (6–5) emission showing systematically more compact structures.

We also determine other parameters such as dynamical timescale ($t_d$), mass outflow rate ($\dot{M}$), mechanical force ($F_m$), kinetic energy ($E_k$), and mechanical luminosity ($L_{mech}$). These quantities are defined as follows: $t_d = R/V$, $\dot{M} = M/t_d$,
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Kinetic temperature ($T_{kin}$) versus H$_2$ density ($n$) plots from LVG calculations. The curves show the observed (3–2)/(6–5), (3–2)/(7–6), and (6–5)/(7–6) integrated intensity ratios (computed from Table 5). The dashed lines represent the errors in the line ratios. The positions are indicated by the labels in each panel.</figcaption></figure>

<table><caption>**Table 6.** Rotational diagram.</caption><thead><tr><th>Position</th><th><i>T</i><sub>rot</sub>(K)</th><th><i>N</i> (10<sup>16</sup> cm<sup>-2</sup>)</th></tr></thead><tbody><tr><td>FIR 6b-B1</td><td>37 ± 9</td><td>1.5 ± 0.8</td></tr><tr><td>FIR 6b-R1</td><td>69 ± 4</td><td>1.2 ± 0.1</td></tr><tr><td>FIR 6c-B1</td><td>67 ± 4</td><td>1.2 ± 0.1</td></tr><tr><td>FIR 6c-B2</td><td>60 ± 10</td><td>1.7 ± 0.5</td></tr><tr><td>MMS 5-B1</td><td>83 ± 12</td><td>2.2 ± 0.4</td></tr><tr><td>MMS 5-R1</td><td>95 ± 65</td><td>0.4 ± 0.3</td></tr><tr><td>MMS 2-R1</td><td>68 ± 5</td><td>0.48 ± 0.05</td></tr><tr><td>MMS 2-B1</td><td>63 ± 4</td><td>0.83 ± 0.09</td></tr></tbody></table>

**Notes.** We have used integrated intensities from Table 5. Statistical and calibration errors are taken into account.

$P = M \times \delta v_{max}$, $F_m = M \times V/t_d$, $E_k = M \times \delta v_{max}^2/2$, $L_{mech} = E_k/t_d$ (e.g. Beuther et al. 2002). In the above expressions, $R = R_{max}/\sin i$ and $V = \delta v_{max}/\cos i$, where $R_{max}$ and $\delta v_{max}$ are the maximum extension and maximum relative velocity of the lobes, and $i$ the angle of the lobe with respect to the plane of the sky. To compare with T08 we assumed $i = 45^\circ$ for all outflows. Table 7 shows our results. Here we also point out that considering the uncertainties in the mass estimation plus the geometry assumptions mentioned above, the kinematic parameters should be accurate within a factor of 10–30 (see, e.g., discussion in Cabrit & Bertout 1990). Within the uncertainties mentioned above our results are consistent with T08.

# 6. Discussion

In the following we first discuss our data in comparison with other similar observations and after that focus on the case of the EHV emission from MMS 5.

## 6.1. Comparison with other mid-J CO observations

The outflow properties determined for our set of IM outflows are not very different from the values found by other authors in low-mass outflows (e.g., Paper I; Curtis et al. 2010). In particular, $F_m$ and $E_k$ values for the IM outflow MMS 5 ($3 \times 10^{-3} M_\odot \text{ km s}^{-1} \text{ yr}^{-1}$ and $2 \times 10^{45} \text{ erg, respectively}$) are similar, within the typical errors, to those found in Paper I for the low-mass class 0 outflow L1448-mm ($\sim 8 \times 10^{-4} M_\odot \text{ km s}^{-1} \text{ yr}^{-1}$ and $2 \times 10^{45} \text{ erg, respectively}$). We therefore conclude that the IM outflow MMS 5 is as energetic as its low-mass counterpart in L1448-mm. However, we note that this may be only a particular case, since previous studies towards statistical significant samples have shown a trend of increasing energy and mechanical luminosity as the protostellar masses increase (e.g., Beuther et al. 2002).

van Kempen et al. (2016) presented CO (6–5) observations of a sample of six IM protostars. They determined outflow parameters employing methods similar to ours, allowing for a direct comparison. In general our results are in agreement with those of van Kempen et al. (2016). In particular, the outflow force is within the same range, except for the case of NGC 2071. More discrepant are the outflow masses which are higher for most of the sources in the van Kempen et al. (2016) sample. The latter is likely due to the larger size of the outflow lobes in that sample, which in turn may be an effect of the age. Given that we obtained a similar outflow force range as found for the van Kempen et al. (2016) sources that span a similar range of bolometric luminosities of the central objects, the correlation between the outflow force and the bolometric luminosity holds for both our sample and theirs.

## 6.2. The EHV collimated jet from MMS 5

Figures 7 and 6 show the CO (6–5) and CO (7–6) high-velocity and EHV emission, overlaid with the NIR emission from the
<table><caption>**Table 7.** Outflow properties derived from CO (6–5).</caption><thead><tr><th>Velocity range</th><th>δv<sub>max</sub><sup>(1)</sup><br>(km s<sup>-1</sup>)</th><th>R<sub>max</sub><sup>(2)</sup><br>(10<sup>-2</sup>)<br>(pc)</th><th>t<sub>d</sub><sup>(3)</sup><br>(10<sup>3</sup>)<br>(yr)</th><th>M(H<sub>2</sub>)<sup>(4)</sup><br>(10<sup>-3</sup>)<br>(M<sub>☉</sub>)</th><th>F<sub>m</sub><sup>(5)</sup><br>(10<sup>-5</sup>)<br>(M<sub>☉</sub> km s<sup>-1</sup> yr<sup>-1</sup>)</th><th>E<sub>k</sub><sup>(6)</sup><br>(10<sup>43</sup>)<br>(erg)</th><th>L<sub>mec</sub><sup>(7)</sup><br>(L<sub>☉</sub>)</th></tr></thead><tbody><tr><td colspan="8" style="text-align: center;">FIR 6b</td></tr><tr><td>Blue</td><td>8</td><td>3</td><td>4</td><td>6</td><td>2</td><td>0.7</td><td>0.02</td></tr><tr><td>Red</td><td>19</td><td>4</td><td>2</td><td>3</td><td>3</td><td>2</td><td>0.07</td></tr><tr><td>Total</td><td></td><td></td><td></td><td>9</td><td>5</td><td>2.7</td><td>0.09</td></tr><tr><td colspan="8" style="text-align: center;">FIR 6c</td></tr><tr><td>Blue</td><td>26</td><td>4</td><td>1</td><td>5</td><td>13</td><td>7</td><td>0.4</td></tr><tr><td colspan="8" style="text-align: center;">MMS 5</td></tr><tr><td>Blue</td><td>106</td><td>6</td><td>0.6</td><td>9</td><td>230</td><td>200</td><td>28</td></tr><tr><td>Red</td><td>59</td><td>4</td><td>0.6</td><td>5</td><td>70</td><td>30</td><td>5</td></tr><tr><td>Total</td><td></td><td></td><td></td><td>14</td><td>300</td><td>230</td><td>33</td></tr><tr><td colspan="8" style="text-align: center;">MMS 2</td></tr><tr><td>Blue</td><td>23</td><td>4</td><td>2</td><td>1</td><td>2</td><td>1</td><td>0.06</td></tr><tr><td>Red</td><td>19</td><td>4</td><td>2</td><td>0.4</td><td>0.5</td><td>0.3</td><td>0.01</td></tr><tr><td>Total</td><td></td><td></td><td></td><td>1.4</td><td>2.5</td><td>1.3</td><td>0.07</td></tr></tbody></table>

**Notes.** <sup>(1)</sup>Maximum relative velocity of the lobe. <sup>(2)</sup>Maximum extension of the lobe. <sup>(3)</sup>Dynamical time. <sup>(4)</sup>H<sub>2</sub> mass. <sup>(5)</sup>Mechanical force. <sup>(6)</sup>Kinetic energy. <sup>(7)</sup>Mechanical luminosity.

MMS 5 region. The main NIR features in this region are an extended monopolar structure and three knots to the west of it. These features are coincident with the CO (6–5) and CO (7–6) blue-shifted emission. On the other hand, no NIR structure is coincident with the red-shifted CO lobe, likely because of the known affect of extinction from the foreground gas in the red-shifted lobes (counter-jet) (see e.g., Noriega-Crespo & Raga 2012). The NIR monopolar structure is indeed well collimated (size ~8.7" × 1.7", hence a collimation factor ~5), which suggests the jet nature of the NIR emission. Since the knots to the west are not well aligned with this monopolar structure, it is possible that the jet from MMS 5 is precessing. As noted in Fig. 6, the EHV CO(6–5) blue-shifted emission seems to be coincident with the NIR collimated monopolar structure and the closest, brightest, NIR knot to the west of it. Therefore, the positional correlation with this NIR structure and knot may suggest the jet nature of the CO EHV emission. The higher angular resolution ALMA data reveal a collimated structure of the EHV gas in the CO (2–1) and SiO (5–4) transitions (Matsushita et al. 2019). The similar line profile and velocity range of the ALMA CO observations indicates that our mid–J CO observations are tracing the same outflow from MMS 5. However, the authors have found that the EHV gas and the NIR jet are not positionally coincident, suggesting that the H<sub>2</sub> knots may not be related to the CO and SiO jets revealed by their data.

To our knowledge Cep-E and NGC 2071 are the only other confirmed cases of EHV emission from IM outflows (Lefloch et al. 2015; Gómez-Ruiz et al. 2012; Chernin & Masson 1992; Hatchell et al. 1999), so the discovery of EHV gas from the MMS 5 outflow is an important addition to this kind of objects that would help to understand the jet phenomenon within the context of star formation theories, in particular for the IM case. Taking into account only the EHV emission and assuming that our CO (6–5) data resolved the major axis of the blue-shifted lobe, a dynamical timescale of ~150 yr is estimated for this

kinematic structure, which is within the range measured in other very young IM outflows (e.g., Zapata et al. 2006). This result highlights the relevance of the mid–J CO lines when searching for outflows in their early phase of evolution.

## 7. Summary and conclusions

The main results of this study on the mid–J CO lines from IM outflows can be summarized in this way:

- The CO (6–5) (and in some cases the CO (7–6)) emission in OMC-2 FIR 6 and OMC-3 MMS 1-6 was found tracing outflows related to class 0-type IM objects within those regions.

- The CO (6–5) line profiles particularly highlight the secondary peaks, likely related to shock structures.

- Extremely high-velocity CO (6–5) emission was detected from the MMS 5 outflow in the OMC-3 region. The CO (7–6) transitions was detected at the 3σ level from the red lobe of the same outflow.

- By comparison with previous observations, we found that at positions close to the central objects the CO (6–5) and (7–6) emission is stronger than the low–J CO (3–2) emission. This behavior is similar to what we found in low-mass outflows (L1448 and HH211).

- The kinematics and physical conditions of MMS 5 outflow suggest that it is the counterpart of the young class 0 low-mass outflows.

From these results we can conclude that the physical and kinematic properties of the young IM outflows presented here are similar to those found in outflows from Class 0 low-mass objects, and that given their excitation requirements, the mid–J CO lines are important tracers of young outflows likely related to jets.

*Acknowledgements.* We wish to thank all the APEX staff in Chile for their enthusiastic help during these observations. AIGR was supported through a stipend
from the International Max-Planck Research School for Astronomy and Astro-
physics at the universities of Bonn and Cologne, and is currently supported
by CONACYT under the program Cátedras CONACYT para jóvenes investi-
gadores. A.G. acknowledges support via the grant ANR-09-BLAN-0231-01 from
the French Agence Nationale de la Recherche as part of the SCHISM project.

References

Aso, Y., Tatematsu, K., Sekimoto, Y., et al. 2000, ApJS, 131, 465

Bachiller, R., Martin-Pintado, J., Tafalla, M., Cernicharo, J., & Lazareff, B. 1990, A&A, 231, 174

Beltrán, M. T. 2015, Ap&SS, 355, 283

Beltrán, M. T., Estalella, R., Girart, J. M., Ho, P. T. P., & Anglada, G. 2008, A&A, 481, 93

Beuther, H., Schilke, P., Sridharan, T. K., et al. 2002, A&A, 383, 892

Cabrit, S., & Bertout, C. 1990, ApJ, 348, 530

Chernin, L. M., & Masson, C. R. 1992, ApJ, 396, L35

Chini, R., Reipurth, B., Ward-Thompson, D., et al. 1997, ApJ, 474, L135

Crimier, N., Ceccarelli, C., Alonso-Albi, T., et al. 2010, A&A, 516, A102

Curtis, E. I., Richer, J. S., Swift, J. J., & Williams, J. P. 2010, MNRAS, 408, 1516

Fuente, A., Neri, R., Martín-Pintado, J., et al. 2001, A&A, 366, 873

Fuente, A., Ceccarelli, C., Neri, R., et al. 2007, A&A, 468, L37

Furlan, E., Fischer, W. J., Ali, B., et al. 2016, ApJS, 224, 5

Gómez-Ruiz, A. I., Gusdorf, A., Leurini, S., et al. 2012, A&A, 542, L9

Gómez-Ruiz, A. I., Wyrowski, F., Gusdorf, A., et al. 2013, A&A, 555, A8

Güsten, R., Nyman, L. Å., Schilke, P., et al. 2006, A&A, 454, L13

Hatchell, J., Fuller, G. A., & Ladd, E. F. 1999, A&A, 346, 278

Johnstone, D., & Bally, J. 1999, ApJ, 510, L49

Kahane, C., Jaber Al-Edhari, A., Ceccarelli, C., et al. 2018, ApJ, 852, 130

Kasemann, C., Güsten, R., Heyminck, S., et al. 2006, SPIE Conf. Ser., 6275, 19

Klein, B., Philipp, S. D., Krämer, I., et al. 2006, A&A, 454, L29

Kounkel, M., Covey, K., Suárez, G., et al. 2018, AJ, 156, 84

Lefloch, B., Gusdorf, A., Codella, C., et al. 2015, A&A, 581, A4

Matsushita, Y., Takahashi, S., Machida, M. N., & Tomisaka, K. 2019, ApJ, 871, 221

Noriega-Crespo, A., & Raga, A. C. 2012, ApJ, 750, 101

Schöier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 369

Shimajiri, Y., Takahashi, S., Takakuwa, S., Saito, M., & Kawabe, R. 2009, PASJ, 61, 1055

Stanke, T., McCaughrean, M. J., & Zinnecker, H. 2002, A&A, 392, 239

Takahashi, S., & Ho, P. T. P. 2012, ApJ, 745, L10

Takahashi, S., Saito, M., Takakuwa, S., & Kawabe, R. 2006, ApJ, 651, 933

Takahashi, S., Saito, M., Ohashi, N., et al. 2008, ApJ, 688, 344

Takahashi, S., Ho, P. T. P., Tang, Y.-W., Kawabe, R., & Saito, M. 2009, ApJ, 704, 1459

Takahashi, S., Saigo, K., Ho, P. T. P., & Tomida, K. 2012, ApJ, 752, 10

Takahashi, S., Ho, P. T. P., Teixeira, P. S., Zapata, L. A., & Su, Y.-N. 2013, ApJ, 763, 57

Tsujimoto, M., Koyama, K., Kobayashi, N., et al. 2004, PASJ, 56, 341

van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, A&A, 468, A627

van Kempen, T. A., Longmore, S. N., Johnstone, D., Pillai, T., & Fuente, A. 2012, ApJ, 751, 137

van Kempen, T. A., Hogerheijde, M. R., van Dishoeck, E. F., et al. 2016, A&A, 587, A17

Williams, J. P., Plambeck, R. L., & Heyer, M. H. 2003, ApJ, 591, 1025

Wilson, T. L., & Matteucci, F. 1992, A&ARv, 4, 1

Yang, B., Stancil, P. C., Balakrishnan, N., & Forrey, R. C. 2010, ApJ, 718, 1062

Yıldız, U. A., Kristensen, L. E., van Dishoeck, E. F., et al. 2015, A&A, 576, A109

Yu, K. C., Bally, J., & Devine, D. 1997, ApJ, 485, L45

Zapata, L. A., Rodríguez, L. F., Kurtz, S. E., & O'Dell, C. R. 2004, AJ, 127, 2252

Zapata, L. A., Rodríguez, L. F., Ho, P. T. P., et al. 2005, ApJ, 630, L85

Zapata, L. A., Ho, P. T. P., Rodríguez, L. F., et al. 2006, ApJ, 653, 398