© ESO, 2014

1. Introduction

The SiO molecule is the most abundantly observed silicon-bearing molecular species in the interstellar medium (ISM), and it can be produced by the photo-evaporation of icy grain mantles (e.g. Turner 1991; Walmsley et al. 1999; Schilke et al. 2001; Godard & Cernicharo 2013) or by the destruction of grain cores in shocks (e.g. Martin-Pintado et al. 1992; Schilke et al. 1997; Le Picard et al. 2001; Gusdorf et al. 2008a,b). The high-velocity shocks from protostellar outflows can easily produce SiO, and therefore, observations of SiO emission are commonly used to trace such outflows (e.g. Beuther et al. 2002; Qiu et al. 2007; López-Sepulcre et al. 2011). Until recently, only fast molecular shocks were believed to be able to produce SiO by sputtering of the grain cores (typically C-shocks with ν_s> 25 km s^-1, but fast J-shocks can also form SiO, see e.g. Guillet et al. 2009). Nevertheless, we do observe significant SiO emission at low velocities and with a small velocity dispersion. For instance, single-dish SiO pointed observations towards the Cygnus-X massive dense clumps (MDCs) by Motte et al. (2007) had shown evidence of two different gas components contributing to the SiO emission: a broad component tracing the shocked gas from protostellar outflows and a narrower one (σ_ν ≲ 1.5 km s^-1, i.e. FWHM ≲ 3.5 km s^-1) interpreted at the time as possibly arising from the hot cores. More recently, some individual studies have detected narrow and widespread SiO emission in molecular clouds, even in regions where no feedback from the forming protostars has yet kicked in (e.g. Jiménez-Serra et al. 2010; Nguyên-Lu’o’ng et al. 2013; Sanhueza et al. 2013; Kauffmann et al. 2013; Louvet et al., in prep.). The origin of SiO is in these cases enigmatic.

Some models suggest that the formation of SiO from low-velocity shocks could be linked to the presence of the element Si in the pre-shock medium. If there is some Si already in the gas phase, shocks with velocities as low as ∼5 km s^-1 can produce SiO (e.g. Nguyên-Lu’o’ng et al. 2013). On the other hand, despite only considering silicates in grain cores, Caselli et al. (1997) demonstrated that grain-grain collisions dominate over sputtering for evaporation of the ices on the grain mantles for low-velocity shocks (10 − 15 km s^-1). In this context, if a fraction of Si is in the grain ice mantles, then shock velocities of ∼10 km s^-1 can desorb Si into the gas phase and effectively increase the abundances of SiO to match observations (e.g. Gusdorf et al. 2008b; Jiménez-Serra et al. 2008). If the Si is solely in the grain cores, then shock velocities above 25 km s^-1 are needed to form SiO (e.g. Schilke et al. 1997; Gusdorf et al. 2008a). However, none of these models include grain-grain interactions. To our knowledge, the most sophisticated account of these interactions has been given by the studies of Guillet et al. (2007, 2009, 2011) and Anderl et al. (2013). When included in shock models, these interactions (namely vaporisation, coagulation, and shattering) do significantly affect the structure of C-type shock waves propagating in high pre-shock density conditions (n_H> 10⁵ cm^-3). Because they result in the production of numerous small grains, the modelled shock layers are both thinner and hotter. Ultimately, even if collisional vaporisation becomes more efficient than sputtering at destroying grains, especially at slightly lower shock velocities (20−25 km s^-1), the overall result of these improved models is less bright and narrower SiO line profiles that do not match the observed intensities, unless a fraction of silicon-bearing material is placed in the grain mantles, similar to what Gusdorf et al. (2008b) had done.

It is vital to further investigate the origin of the measured SiO components, to be able to disentangle SiO outflows from other possible forms of shocks. However, most existing surveys observing SiO emission towards star-forming regions do not possess enough angular resolution or sensitivity to determine the exact origin of the SiO emission, and it becomes impossible to disentangle any emission that may not arise from outflows. On the one hand, ignoring other sources of shocks for the origin of SiO could introduce significant uncertainties on any outflow properties derived from SiO luminosities alone. On the other hand, pinpointing the exact origin of SiO emission can help to understand the physics and mechanisms for accretion at core-scales.

In this article, we present SiO (2 − 1) interferometric observations towards a sample of IR-quiet MDCs in Cygnus-X, located at a distance of 1.4 kpc (Rygl et al. 2012), which are known to have both broad and narrow SiO emission (Motte et al. 2007), and for which we have detailed information about the protostellar content and the existing outflows (Duarte-Cabral et al. 2013). We address the consequences of not resolving the SiO emission for inferring outflow properties, and we investigate the different possibilities for the origin of the observed SiO emission. Details on the source sample and observations are provided in Sect. 2. We describe our results in Sect. 3, which include a comparison of the spectral profiles of SiO with those obtained with single-dish observations, and a comparison of the SiO and CO outflow emission. We analyse the different trends observed with SiO and the emission morphologies in Sect. 4, and present our concluding remarks in Sect. 5.

2. Source sample

2.1. Cygnus-X massive dense clumps

We have studied the six IR-quiet MDCs in Cygnus-X, originally identified by Motte et al. (2007) and more recently observed with the PdBI in the 1 mm and 3 mm continuum emission (Bontemps et al. 2010) and ¹²CO (2 − 1) outflow emission (Duarte-Cabral et al. 2013). Table 1 lists the name and properties of each MDC, including the number of fragments and morphology.

These MDCs are found in a range of different environments, and each field, with the exception of CygX-N40, is subfragmented into several cores with envelope masses larger than ∼10 M_⊙ (as estimated from the 1mm emission by Bontemps et al. 2010; Duarte-Cabral et al. 2013). The CygX-N3 and N12 MDCs are situated to the west of the DR21 ridge and they both show elongated morphologies, consistent with being shaped by the winds from the nearby OB clusters (to the N and NW). Both are fragmented into several compact cores, with two high-mass cores per field (see Table 1). CygX-N40, N48, and N53 are located along the massive DR21 ridge, a highly dynamical star formation site (Schneider et al. 2010; Hennemann et al. 2012). The interferometric continuum data of CygX-N40 (from Bontemps et al. 2010) showed that this core is not centrally condensed, with its emission being extended and filtered out by the PdBI. Bontemps et al. (2010) had detected a single low-mass fragment in this region, with M_env< 2 M_⊙. CygX-N48 has a clustered star formation, not very centrally condensed, but with two massive compact sources, while CygX-N53 is centrally condensed and fragmented into two compact massive cores (with more than 20 M_⊙ each). Finally, CygX-N63 is a relatively isolated dense core to the south of DR21. It does not show signs of fragmentation and has nearly 50 M_⊙ within its inner ∼2500 AU, making it the most massive core of the sample.

In Duarte-Cabral et al. (2013), we confirmed the nature of these millimetre cores as high-mass Class 0-like protostars, based on their masses, luminosities, and outflow power/accretion rates. We identified the individual outflows powered by each massive protostar, and estimated their properties, concluding that the opening angles and morphologies are similar to those of low-mass objects, and they are 1−2 orders of magnitude more powerful, following a linear trend with the envelope mass.

2.2. Observations

The six MDCs were observed in 2004 with the PdBI in the 1.3 mm and 3.5 mm continuum emission and in four spectral units covering the ¹²CO (2−1), SiO (2−1), H¹³CO⁺ (1−0), and H¹³CN (1−0) lines. The observations in the continuum and in the H¹³CO⁺ (1 − 0) and H¹³CN (1 − 0) lines are reported in Bontemps et al. (2010) and Csengeri et al. (2011a). The ¹²CO (2 − 1) observations (shown in Duarte-Cabral et al. 2013) have been short-spaced using IRAM 30 m data, and will be used here to aid the identification of individual outflows. This article presents the SiO (2−1) observations, at 86.85 GHz, with an angular resolution of ∼3^′′ (i.e. ∼0.02 pc at 1.4 kpc distance), tracing interstellar shocks at core-scales (typical protostellar outflows being ∼0.1−0.4 pc in length).

The observations were performed in track-sharing mode with two targets per track for the following pairs: CygX-N48/N53, CygX-N3/N40, and CygX-N12/N63. The D configuration track observations were performed between June and October 2004 (five antennas with baselines between 24 m and 82 m). The C configuration tracks were obtained in November and December 2004 (six antennas in 6Cp with baselines from 48 m to 229 m). The bright nearby quasar 2013+370 was used as a phase calibrator, and the evolved star MWC349 as a flux calibrator.

The maps were cleaned using the natural weighting, and the resulting synthesised beam and root mean square (rms) in the continuum are summarised in Table 1, together with the field names and centres of phase. The pixel size is of 0.62^′′. The cleaning components were searched across the whole area of the primary beam. No support for cleaning was used to avoid introducing any bias into the resulting emission maps.

3. Results

3.1. Profile of the SiO emission

Motte et al. (2007) presented single-dish pointed observations of SiO (2−1) emission with the IRAM 30 m telescope (beam size of ∼29^′′) towards the sample of IR-quiet cores we study here. Figure 1 shows the SiO average profiles from our PdBI observations as dark grey shadowed spectra (averaged within the central 30^′′ × 30^′′) overlaid with the dotted blue spectra from the single-dish observations.

There is a good correspondence between the two profiles for 3 out of 6 sources (CygX-N3, N48 and N63). For these three MDCs, the PdBI recovers most of the flux detected with the single-dish (95%, 80%, and >100% for CygX-N3, N48, and N63, respectively; see Table 2). The fact that we recover slightly more emission with the PdBI with respect to the IRAM 30 m telescope in N63 can be explained by the uncertainties on the calibration of the IRAM 30 m data (∼10%), by some loss of power towards the edges of the single-dish primary beam, and by the fact that the area we used to extract the average PdBI profiles is a rectangular area of 30^′′ width, instead of a circular beam shape. In CygX-N53, even though the shape of the profile is similar, we miss some of the systemic velocity emission with the PdBI. Nevertheless, we still recover ∼70% of the single-dish flux with the PdBI. This indicates that most of the SiO emission is compact on these four sources, and not filtered out with the PdBI.

The two most striking exceptions are CygX-N12 and N40, where a significant fraction of the single-dish flux (∼65−70%) is missing with the PdBI, when only the central 30^′′ emission is taken, either because some flux is filtered out in the interferometric observations or because the single-dish observations suffer from contamination from strong emission sitting in the immediacies of the single-dish primary beam. In order to discriminate between these possibilities, we have integrated the PdBI spectra for these two regions considering a larger area of 50^′′ × 50^′′ (light grey shadowed spectra in the panels of these two sources). Doing so for CygX-N12 recovers not only the flux of the single-dish observations, but also the line shape. We conclude that this is due to the contamination from the region east of the single-dish primary beam, with strong outflow emission (see central panels of Fig. 2). Doing this for CygX-N40, despite a considerable improvement, still does not recover the peak flux detected with the single-dish at systemic velocities, suggesting that there is a considerable amount of extended SiO emission at systemic velocities which is filtered by the PdBI (with a need for short-spacings). On the other hand, it recovers some line-wings that were not seen with the single-dish.

Fig. 1 Average spectra for the six MDCs of our sample (labelled on the top-left corners). The left side of each subfigure shows the comparison of the single-dish (SD) IRAM 30 m pointed observations (Motte et al. 2007) in blue dashed lines, and the spectra of our PdBI SiO observations integrated over the central 30′′ × 30′′ (equivalent to the IRAM 30 m beam) are shadowed in grey. The vertical black dashed lines show the systemic velocity of each core, as measured with N2H+ by Bontemps et al. (2010). The right side shows the same PdBI SiO emission as the left, overplotted with a two-component Gaussian fit to illustrate the existence of narrow (in blue) and broad (in green) components of emission (whose sum is shown in red). For CygX-N12 and N40, we show a shadowed spectrum (in light grey) which corresponds to the PdBI emission integrated over an area of 50′′ × 50′′ (to assess a possible contamination from strong SiO emission outside the SD primary beam).

Fig. 2 Comparison of the high-velocity CO and SiO emission in CygX-N3 (left), N12 (centre), and N40 (right). The top panels show the high-velocity CO emission in contours (with the velocity ranges as in Duarte-Cabral et al. 2013) over the 1 mm continuum emission in greyscale. The primary beam of the CO and 1mm observations is shown as a dashed circle, and the synthetic beam is in the lower-left corner. Contours are stepped by 0.25 Jy beam-1 km s-1, starting at 2.0 (2.75) Jy beam-1 km s-1 for N3 blue (red), 1.5 Jy beam-1 km s-1 for N12, and 1.75 (1.25) Jy beam-1 km s-1 for N40 blue (red). The lower panels show the high-velocity SiO emission in contours (with velocity ranges as in Table 2) over the 3 mm continuum emission in greyscale for the same sources. The primary beams are shown as dashed circles and the synthetic beams are plotted in the lower-left corners. Contours are stepped by 0.1 Jy beam-1 km s-1, starting at 0.15 Jy beam-1 km s-1 for N3, 0.2 (0.05) Jy beam-1 km s-1 for N12 blue (red), and 0.2 Jy beam-1 km s-1 for N40. The most massive protostars are marked with yellow stars (and labelled in the lower panels), and the least massive fragments from Bontemps et al. (2010) are shown as green crosses. The arrows are the directions of the outflows inferred from Duarte-Cabral et al. (2013), and the dashed lines show the respective outflow cones.

Fig. 3 Same as Fig. 2 for CygX-N48 (left), N53 (centre), and N63 (right). In N48 we mark with circled crosses two Spitzer sources: ERO2 (Marston et al. 2004) and IS-1, which is a resolved source detected in the IRAC bands at 3 and 4 μm, and unresolved at 8 μm and onwards. IS-1 is likely responsible for the outflow directed east-west. Contour key for CO emission, top panels: N48 (left) contour steps of 0.5 Jy beam-1 km s-1 starting at 4 (3) Jy beam-1 km s-1 for blue (red); N53 (centre) contour steps of 1.5 Jy beam-1 km s-1 starting at 8.5 Jy beam-1 km s-1; N63 (right) contour steps of 1.5 Jy beam-1 km s-1 starting at 5.0 Jy beam-1 km s-1. Contours of SiO emission, lower panels: N48 (left) contour steps of 0.1 Jy beam-1 km s-1, starting at 0.2 (0.15) Jy beam-1 km s-1 for blue (red); N53 (centre) contour steps of 0.15 (0.25) Jy beam-1 km s-1, starting at 0.2 (0.15) Jy beam-1 km s-1 for blue (red); N63 (right) contour steps of 0.5 Jy beam-1 km s-1, starting at 0.2 Jy beam-1 km s-1, for both blue- and redshifted emission.

Based on the single-dish SiO emission, Motte et al. (2007) discussed the existence of two components of SiO emission. Our observations confirm that the average spectra of most sources can be accurately represented by the emission of two Gaussian components (Fig. 1), with the exception of CygX-N53 which shows a clear asymmetric profile, with broad redshifted emission (reaching a terminal velocity of ∼30 km s^-1). In Table 2 we show the systemic velocities, ν₀, of the different regions as measured by Bontemps et al. (2010) using N₂H⁺ pointed observations, also performed with the IRAM 30 m telescope (and therefore within an equivalent beam size). The uncertainties on these velocities are of 0.1 km s^-1. These are shown as dashed vertical lines in Fig. 1. For the MDCs where two Gaussian components fit well the observed average spectra, we present, in Table 2, the peak intensity, I_peak; the peak velocity, ν_peak; and the velocity dispersion, σ_ν, for both the narrow and the broad components. The peak of the narrow component of the SiO emission is only slightly different from the systemic velocities (less than 1 km s^-1 difference), which is not significant considering the linewidths of the emission and the uncertainties on the peak velocity. This is with the exception of N12, where there is an important shift between the two (of 2 km s^-1). We will discuss possible explanations for this shift later in the article (Sect. 4.2.1).

While the broad emission whose average width at the base reaches ≳30 km s^-1 suggests unequivocally outflow shocks, Motte et al. (2007) considered the possibility of the narrower component originating from shocks inside a hot core and/or lower velocity outflows.The only way to distinguish between a high-velocity or low-velocity shock origin is by carefully inspecting the spectra along the maps, and investigate the spatial distribution of the systemic-velocity and broad emissions (see Sects. 3.2 and 3.3). The spectra throughout the maps vary significantly from the average spectra, and examples of the different line profiles found in each MDC can be found in Appendix A, Fig. A.2.

3.2. Distribution of high-velocity SiO versus CO

To investigate how effectively SiO is tracing protostellar outflows, we compare the spatial distribution of the high-velocity SiO emission with that of CO (Duarte-Cabral et al. 2013). Table 2 shows the velocity ranges used to estimate the SiO integrated intensities and luminosities at high-velocities (high-ν range), which is in essence the emission offset from the systemic velocity of the cloud by more than ±3 km s^-1 (corresponding to roughly 3σ_v from the fits to the N₂H⁺ (1−0) spectra by Bontemps et al. 2010). Figures 2 and 3 show the high-velocity CO and SiO emission in the top and bottom panels, respectively. In these figures we also show the outflow cones (dotted lines) and directions (arrows). Although with a different resolution (∼1^′′) and field of view (see dashed circles in Figs. 2 and 3), the CO is used as a reference to identify the different outflow lobes and respective driving sources.

From a fast inspection of these figures we can see a number of similarities in the distribution of the high-velocity SiO and high-velocity CO outflow emission. In particular, for CygX-N3 (left column of Fig. 2) and CygX-N63 (right column of Fig. 3), the correspondence is quite clear both spatially and in terms of relative intensity of the blue- and redshifted outflow lobes.

In other regions, the correspondence in terms of spatial extent and intensities is poorer. This is the case, for instance, of CygX-N53 (central column of Fig. 3) where SiO is almost only redshifted while CO is roughly symmetric. This is puzzling. Perhaps this is due to the powerful red lobe encountering more material which is efficiently shocked to form and entrain SiO. On the other hand, the blueshifted SiO traces only vaguely the blueshifted emission of CO. In addition, the redshifted emission seen towards the west could be due to an additional outflow, since there is a blueshifted counterpart in the same direction towards the east, at the edge of the primary beam. However, although it may be that the corresponding CO emission is outside our field of view, the absence of such outflow in CO makes it hard to understand its origin. It could be that it arises from CygX-N53 MM2, but this is not the outflow we consider for this source. Duarte-Cabral et al. (2013) suggest a tentative detection of a very compact outflow using CO emission for which we estimate the respective SiO luminosities here (and hence these are only to be taken as an upper limit). For CygX-N48 (left column of Fig. 3), there is also some discrepancy between the intensities of the two molecules, namely on the outflow powered by MM2. In this region, only 25% of the high-velocity SiO luminosity arises from the outflows of MM1 and MM2, and most high-velocity outflow emission in this field is due to the outflow from IS-1 (IS-1 stands for Infrared Source 1, and corresponds to a resolved source detected in the IRAC bands at 3 and 4 μm, and unresolved at 8 μm and onwards, at RA 20^h39^m02.927^s and Dec 42^°22^′07.32^′′; using Herschel, Hennemann et al. (in prep.) define this source as a more evolved object). This demonstrates the need to identify individual outflows to perform good evaluations of individual outflow power (unresolved observations of this region would have assigned the emission from IS-1 outflow to the stronger millimetre peak, MM1).

Fig. 4 Integrated intensity of the SiO emission at ambient velocities with the contribution from the broad outflow emission subtracted, for CygX-N3 (left), N12 (centre), and N40 (right). The top row panels show the narrow emission as greyscale, with contours from the high-velocity (broad) SiO emission from Fig. 2 overlaid. The two lower rows show the integrated narrow SiO emission in contours (with steps of 0.1 Jy beam-1 km s-1, starting at 0.1 Jy beam-1 km s-1), with the colour scale showing the respective moment maps: the velocity field (first moment map) in the middle row, and the velocity dispersion (second moment map) in the last row. Annotations are as in Fig. 2.

Finally, other regions have high-velocity SiO emission at slightly different positions to where we had previously detected the CO outflows. For instance, in CygX-N12 (central column of Fig. 2), because of the larger field of view of the SiO observations, any outflow at the edge of the SiO maps was not covered by CO. Close to the sources, the high-velocity emission from CO and SiO is similar (redshifted), and the larger extent of the SiO emission to the east could indicate that the red CO outflow becomes blue farther away, which suggests that this is an outflow cone close to the plane of the sky. Although it could be that the SiO blue lobe is simply part of another outflow (and unrelated with the CO red outflow lobe), its small velocity range (∼10 km s^-1) and large spatial extent (∼20^′′, i.e. ∼0.13 pc) support the hypothesis of an outflow close to the plane of the sky. There is also some blue high-velocity SiO emission to the north which could be due to a second outflow not detected in CO (and which would follow the opposite direction of the south-directed CO lobe, marked with a dashed line). Using CO only, in Duarte-Cabral et al. (2013), we had interpreted this region as having two overlapping outflows, both E-W directed. With SiO, this picture could perhaps be revised with the existence of a N-S outflow, but because the emission does not trace back clearly to any of the two sources, we have refrained from doing so. Instead, and for the remainder of the paper, we will take into account only the outflow directed E-W. In CygX-N40 (right column of Fig. 2), the outflow emission we had detected with CO in Duarte-Cabral et al. (2013) close to MM1 has no SiO counterpart (perhaps it is too weak to be detected within out noise levels), while the two weak high-velocity lobes to the north-east and north-west of MM1 are detected with both molecules. This, along with the existing weak 3 mm peak, lead us to think that there is a second low-mass source (M_env ≲ 1 M_⊙), to the north of MM1, and which could be responsible for the additional E-W CO and SiO outflow. We will refer to this source as CygX-N40 MM2, whose position and estimate of the F_CO (as done for the other sources in Duarte-Cabral et al. 2013) is shown in Table. 3. In this region, only 10% of the total PdBI SiO luminosity is close to the central sources (MM2 in particular), 90% of which is high-velocity emission.

3.3. Distribution of the narrow SiO emission

We have investigated the distribution of the SiO emission close to the systemic velocities of the cloud, and whose main contribution is adequately fitted by a narrow Gaussian (see spectra in Fig. 1). Since sometimes the narrow SiO emission is along the same line of sight as some high-velocity gas, separating the contribution of the outflows from the emission at systemic velocities becomes critical. In an exercise to try and remove this contribution, we have created datacubes of the narrow emission of each field, by subtracting the contribution from a broad outflow profile whenever there was significant high-velocity emission. To do so, we selected all the pixels where the integrated high-velocity SiO emission was above the rms noise level, and fitted a double Gaussian profile to the spectrum at each pixel. The fitted broad Gaussian component was then subtracted from the original spectra. By doing so, we do not make an a priori assumption on the profile of the narrow component. Instead, we remove the broad contribution and analyse what is left from this process (see Fig. A.1 in Appendix A for an example spectrum before and after the subtraction of the broad component). Figures 4 and 5 show the spatial distribution of the moment maps (integrated intensities, velocity field, and velocity dispersion) of the SiO emission at ambient cloud velocities as a result of this exercise. The three moments were calculated using the central 6 km s^-1 velocity range, centred on the ν₀ value given in Table 2. Since the narrow component (as measured from the average spectra) has a  km s^-1, this velocity range is enough to contain more than 95% of the narrow emission. Naturally, this method is only correct if/when the broad outflow emission is Gaussian, and when both narrow and broad components are well distinguished, but this is not always the case. Therefore, the resulting narrow emission maps may still show some residual broad emission (see the velocity dispersion maps of Figs. 4 and 5). For the purpose of these figures, we degraded slightly the spatial resolution of the integrated intensity maps (i.e. the zeroth-moment maps), by convolving with a Gaussian kernel of two arcseconds, merely so as to smooth the appearance of the contours, namely around individual pixels where the broad component was not successfully removed¹. We did not do this spatial convolution for the velocity field and velocity dispersion maps, and the latter can be used to spot the positions where this method was not successful (pixels that still have large velocity dispersions, with  km s^-1).

Fig. 5 Same as Fig. 4, for CygX-N48 (left), N53 (centre), and N63 (right). Top panels have the narrow SiO emission in greyscale and the high-velocity (broad) SiO emission in contours as in Fig. 3. The two lower rows show the integrated narrow SiO emission in contours (with steps of 0.1 Jy beam-1 km s-1, starting at 0.1 Jy beam-1 km s-1), with the respective moment maps in colour scale: the velocity field (first moment map) in the middle row, and the velocity dispersion (second moment map) in the last row. Annotations are as in Fig. 3.

In CygX-N3 (left column of Fig. 4), based on the location of the CO outflows, we have determined that about 60% of the total PdBI SiO luminosity is not constrained to the location of known outflows (as defined by the existence of SiO high-velocity emission). There is, therefore, a significant fraction of PdBI SiO emission that is not clearly associated with any outflow, and which has narrow line profiles (σ_ν ∼ 1.2 km s^-1, i.e. FWHM ∼ 2.8 km s^-1). The velocity field of this narrow component shows a large velocity gradient across the sources (of 2 km s^-1 in 0.05 pc), with a shear-like structure aligned with the continuum filamentary structure. Interestingly, the velocity gradient seen in SiO is opposite to that of the H¹³CO⁺ flows seen by Csengeri et al. (2011a, also with a velocity gradient of ∼2 km s^-1.

In CygX-N12 (central column of Fig. 4), 45% of the PdBI SiO luminosity arises from narrow low-velocity gas (σ_ν ∼ 0.8 km s^-1, i.e. FWHM ∼ 1.8 km s^-1). This includes the narrow blueshifted emission towards the southern edge of the blueshifted outflow lobe (see blue emission in the velocity field panel of Fig. 4, or narrow spectra from Fig. A.2), which accounts for the missed emission from the inner 30^′′ but recovered in the 50^′′ area. In this particular case, the narrow SiO emission is not spatially overlapping with an outflow and yet it is clearly associated with one, perhaps indicating that we are seeing the post-shock gas at the edges of an outflow cavity. The remaining ambient gas has a spatial distribution slightly extended along the N-S direction around the central sources, with a rather complex velocity field which makes it hard to compare with the (also complex) velocity field from H¹³CO⁺ (Csengeri et al. 2011a).

For CygX-N40 (right column of Fig. 4), a significant amount of the total PdBI emission arises from systemic velocities (54%), in the form of extended narrow emission, with a velocity dispersion reaching as low as σ_ν ∼ 0.3 km s^-1 (i.e. FWHM ∼ 0.7 km s^-1). Some of the SiO emission is associated with DR21(OH)N1 (or CygX-N41 from Motte et al. 2007), seen as a 3mm continuum peak to the south of MM1. The existence of negative side-lobes in SiO between MM1 and the arc-shaped emission around MM1 and MM2 means that despite recovering a good fraction of the emission observed with single-dish, the emission is quite extended and filtered out by the PdBI. Therefore, the fidelity of the imaging may not be very good (and the arc structure may not be real). Nevertheless, the existence of such emission is significant. It is not associated with the central cores, nor with outflows, and its origin is to be elucidated.

For CygX-N48 (left column of Fig. 5), even though most of the PdBI SiO luminosity is associated with low-velocity gas (∼70%), there is a significant spatial overlap of the narrow SiO emission with high-velocity SiO, suggesting that a fraction of the low-velocity SiO emission could be associated with the outflows of the region (e.g. the emission associated with the IS-1, MM1, and/or MM2 outflows). Nevertheless, there is some narrow SiO emission which lies outside outflow areas (∼20%). The origin of this emission could be associated with the low-velocity shocks detected in this region by Csengeri et al. (2011b) in N₂H⁺ (see Sect. 4.2.2 below).

In CygX-N53 (central column of Fig. 5) there is only a relatively small amount of PdBI SiO emission at systemic velocities (∼35%). This emission shows a SE-NW elongation that has a velocity gradient of 2 km s^-1 in 0.1 pc, with a direction reminiscent of the rotation of the MDC and the envelope (as seen in H¹³CO⁺; see Fig. 5 of Csengeri et al. 2011a). The narrow SiO emission to the south, however, could be tracing the lower-velocities of the redshifted outflowing gas.

Finally, in CygX-N63 (right column of Fig. 5), the narrow emission is broken up in two areas. One to the north of (and including) MM1, which shows a velocity gradient across MM1 (E-W), with a direction consistent with that seen in H¹³CO⁺ (Csengeri et al. 2011a). This emission is coincident with the blueshifted outflow lobe, and the sharp velocity gradient of ∼2 km s^-1 in 0.05pc is coincident with the eastern edge of the lobe, perpendicular to the direction of the outflow. This could indicate that the SiO is tracing the cavity walls of the outflow. Farther south, there is another E-W extended emission that has narrow linewidths (σ_ν ∼ 0.8 km s^-1, i.e. FWHM ∼ 1.8 km s^-1). While the central peak could be interpreted as being part of the redshifted outflow walls, the origin of the extensions to the east and west are less clear.

3.4. Narrow SiO column densities

To estimate the order of magnitude of the column densities of the SiO narrow component, we have used the IRAM 30 m spectra of SiO (from Motte et al. 2007). Since these spectra include a contribution from broad outflow emission, we performed double Gaussian fittings to extract the narrow component. The H₂ volume densities in the MDCs are of the order of ∼10⁶ − 10⁷ cm^-3 (∼10⁶ cm^-3 being the average MDC densities, while 10⁷ cm^-3 is more representative of the average density of the protostellar envelopes; see e.g. Louvet et al. 2014), which means that we are above the SiO critical densities, and hence we can assume LTE to estimate the column densities (as e.g. in Nguyên-Lu’o’ng et al. 2013). Assuming a temperature of 20−40 K, we estimate that the observed beam-averaged SiO integrated intensities for the narrow component correspond to SiO column densities of the order of 0.1−4 × 10¹² cm^-2, where the optical depths are largely below one. We note, however, that this corresponds to the beam-averaged SiO column densities, within the IRAM 30 m beam. From the PdBI, we can see that the emission is more compact than the 30^′′ beam, and therefore these column densities are merely a lower limit for the effective column densities. The PdBI emission having linear dimensions typically of the order of ∼15^′′ (i.e. ∼0.1 pc) implies that the beam filling factors are easily low enough to bring the column densities to values closer to ∼10¹³ cm^-2. Since the H₂ column densities in the MDCs range between 6 × 10²³ cm^-2 and 2 × 10²⁴ cm^-2, we retrieve SiO abundances of the order of 0.5−1.0 × 10^-11 for the narrow SiO emission (which are similar to the abundance values of narrow SiO emission found in IRDCs by, e.g. Jiménez-Serra et al. 2010; Sanhueza et al. 2013).

Fig. 6 Observed LSiO against the corresponding FCO (from Duarte-Cabral et al. 2013). The light blue filled circles correspond to the total LSiO and total FCO for each of the six fields, and the blue dashed line is the respective linear fit in a log-scale. The open circles are for the resolved individual outflows (integration spatially constrained to location of the individual outflows), with the FCO estimated in Duarte-Cabral et al. (2013; Table 3, Col. 6). The red open circles have LSiO estimated using the entire velocity range (Table 3, Col. 7), and the respective linear fit to the log-values is shown as a dashed red line, which translates into a linear correlation of LSiO ∝ FCO. The green open circles show the LSiO only from the outflow high-velocity wing emission (Table 3, Col. 8).

4. Analysis

4.1. SiO outflow luminosity

To investigate the relation between the SiO luminosity, L_SiO, and the respective CO outflow momentum flux, F_CO, we have estimated the L_SiO associated with a given outflow, whenever possible. In practice, this involved the usage of polygons to limit the areas where to estimate the SiO flux. These areas were determined based on the extent and morphology of the SiO emission, but with the help of the CO emission presented in Duarte-Cabral et al. (2013), as the CO observations allow a more accurate identification of the driving sources of the different outflow lobes. We estimated this L_SiO for each source, first using the total velocity range, and then by restraining the calculations to the high-velocity SiO emission.

Table 3 summarises the results from this exercise, for the sample of sources for which we have been able to individually estimate the SiO luminosity. Columns 1 to 6 indicate the source names and positions, the envelope mass and bolometric luminosities retrieved from SED fittings, and the CO outflow momentum flux (the last three are from Duarte-Cabral et al. 2013). Columns 7 and 8 show the SiO luminosities associated with individual outflows, and the second-to-last column shows the total SiO luminosity in the entire PdBI field (equivalent to the integrated SiO emission that would be observed over the same field by a single-dish telescope). The last column shows the SiO luminosity in the entire PdBI field that arises from narrow emission (i.e. calculated from the integrated intensity maps shown in Figs. 4 and 5).

The first result from this estimate demonstrates that there is no 1:1 relation between L_SiO and F_CO. Nevertheless, the two stronger outflows in the sample (from the two most massive sources, CygX-N53 MM1 and CygX-N63 MM1) are also those showing the most luminous SiO emission. This is illustrated well in Fig. 6, where we can see that the two stronger outflows (with F_CO> 2 × 10^-3M_⊙ km s^-1 yr^-1) also have strong SiO luminosities (both when taking the entire fields, and the individual sources). However, for the remaining fields, a correlation is less evident. When measuring the SiO luminosities on entire fields, which in essence mimics observations of SiO that do not resolve individual outflows, we see no correlation of SiO luminosity with the corresponding outflow power (the dashed blue line in Fig. 6 shows the tentative fit to the light blue circles, with a weak dependency of ${\mathit{L}}_{\mathrm{SiO}}\mathrm{\propto }{\mathit{F}}_{\mathrm{CO}}^{\hspace{0.17em}\hspace{0.17em}\mathit{\alpha }\lesssim \mathrm{0.3}}$). When we are able to distinguish individual outflows, then the SiO luminosities do appear to correlate linearly with the outflow momentum fluxes measured from CO lines (see dashed red line in Fig. 6, which corresponds to L_SiO ∝ F_CO). This result suggests that SiO luminosities associated with individual outflows could statistically be used as a proxy for the outflow power. However, as shown in Table 3, there is a significant scatter of more than one order of magnitude for the SiO luminosities, for similar values of F_CO. The larger scatter in L_SiO is also visible in Fig. 7, where we show the relation between both L_SiO (in blue) and F_CO (in red) with the M_env, for the protostars of our sample. This figure shows that both L_SiO and F_CO scale linearly with envelope mass, when individual outflows can be resolved, in spite of the larger apparent scatter from the L_SiO. We caution that these results are based on our small sample, and need to be confirmed with larger samples.

Fig. 7 Observed LSiO (same as the red open circles in Fig. 6) and FCO, in blue and red circles, respectively, with respect to the envelope masses (Menv) for individual sources. The second red point from the left corresponds to CygX-N40 MM1, for which we have an FCO measurement, but no LSiO. The dashed red line shows the linear correlation of FCO = 7.5 × 10-5Menv, originally found by Bontemps et al. (1996) for low-mass protostars, and confirmed in Duarte-Cabral et al. (2013) to extend to the high-mass regime. The dashed blue line shows the linear correlation of LSiO ≃ 5 × 106Menv that we find for the blue points. This figure accentuates the fact that LSiO of individual outflows follows the same trend as FCO (both being ∝Menv), even though it introduces more scatter with respect to the linear relation found for FCO.

4.2. Nature of the narrow SiO emission

The fraction of SiO emission at low velocities is highly variable (Table 3), and while in regions with strong outflows (e.g. CygX-N53 and CygX-N63) the SiO high-velocity emission accounts for 80% of the full SiO luminosity, other sources have the inverse situation, where most of the SiO emission arises from narrow emission at systemic velocities. Such emission could be due to shocks other than outflows. To help clarify the nature of the narrow SiO emission, we summarise the results from Sect. 3.3 in Table 4, by describing the emission in the regions according to four distinct cases: (1) the existence of spatially overlapping broad and narrow emission; (2) the existence of broad emission with little/no narrow emission associated; (3) the existence of narrow emission adjacent to (but not coincident with) broad emission; and (4) the existence of narrow emission with absolutely no broad emission in the vicinity. We will look in more detail into each case in the following sections.

Fig. 8 Fraction of the total SiO luminosity that arises from the narrow component, fnarrow, as a function of the mass of the most massive protostellar envelope which is currently forming in each MDC. This plot accentuates the trend of a higher fraction of the emission in the form of narrow emission for the MDCs which have not yet gathered most of their mass in the envelopes of the forming protostars.

4.2.1. Narrow SiO emission from outflows

Strong high-velocity outflows (as those powered by the two most massive sources CygX-N53 and CygX-N63), have naturally a strong impact on the ambient gas, and can easily form and entrain SiO to high velocities. There are, however, several factors that come into play when SiO is formed in such a high-velocity shock, and these factors can change the observed line profile.

The typical profile of SiO produced in high-velocity shocks has a highly asymmetric profile, including a broad and a narrow component (from 1D models, e.g. Gusdorf et al. 2008a; Anderl et al. 2013). This typical profile explains case (1) of Table 4, where broad and narrow emission are spatially overlapping. While such models predict the narrow profiles to be shifted by tens of km s^-1 from the systemic velocities, projection effects can contribute to significantly decrease the observed velocity shifts of the shocked gas with respect to the ambient velocities. On the other hand, projection effects can also have another consequence, which is the fact that the profile observed is a result of the superposition of several ideal 1D profiles. One of the results of this can be the symmetrisation of the asymmetric profiles, i.e. the narrow component will become less obvious and the emission we see is mostly a broad (near Gaussian) profile. This explains the profiles of case (2) of Table 4.

If, however, the composition of the pre-shock gas is such that the silicon is solely in the grain cores, the high-velocity shocks will produce SiO with narrower profiles (Gusdorf et al. 2008a; Anderl et al. 2013). Therefore, the existence of narrow emission offset from the systemic velocities, with no broad emission associated, might be the signature of a fraction of SiO as emanating from the grain cores. Another way to have narrow SiO emission associated with high-velocity shocks, is as a result of the mixing of the shocked gas with the colder ambient gas, which can still form SiO but with lower velocity dispersion (e.g. Lefloch et al. 1998; Codella et al. 1999). Even though the timescales for the deceleration of the shocked gas or alternatively the oxidisation of SiO or its depletion back into the grain ice mantles are short (both of the order of ∼10⁴ yr, e.g. Martin-Pintado et al. 1992; Codella et al. 1999), the high-mass protostars in these MDCs are Class 0 equivalents, and therefore also relatively young. Furthermore, the outflows being a continuous process, they have the potential to continuously replenish the decaying decelerated SiO in outskirts of outflows. In this case, we expect to see narrow SiO emission at systemic velocities, in the vicinity of high-velocity shocks, for example, along outflow walls. Alternatively, and as suggested by Jiménez-Serra et al. (2004), the presence of narrow SiO at systemic velocities could also be explained by recent sputtering of the grain mantles containing a small fraction of Si/SiO and therefore could be indicative of young shocks (or simply tracing of the most recently shocked material). All three options could explain configuration (3) of Table 4.

Finally, another way that we could achieve narrow profiles from a high-velocity shock would be with a strong UV field from the protostar, illuminating either the pre- or post-shock region, effectively decreasing the thickness of the SiO emitting-region, and subsequently narrowing its line profiles (Gusdorf et al., in prep.). However, our sample of protostars (some of which host hot cores, e.g. CygX-N53 MM1 and CygX-N63 MM1) are Class 0 equivalents (with L_bol ∼ 100−500 L_⊙, see Duarte-Cabral et al. 2013, and Table 3), and therefore the radiation from the forming protostars is not yet strong enough to ionise the gas outside the inner protostellar envelopes/hot-core regions (which are unresolved by our observations), ruling out this hypothesis.

In our study, we do find narrow SiO line profiles associated with high-velocity outflow shocks. On one hand, we have narrow SiO emission at systemic velocities close to high-velocity outflow lobes, which could be due to the mixing of shocked gas and ambient gas along outflow walls. The most obvious candidates for this specific case in our sample are CygX-N40, next to the red outflow lobe; CygX-N53, along the red outflow cone, and to the west of the blue outflow lobe; and CygX-N63, around the blue outflow lobe, and at the tip of the red outflow. In CygX-N48 this process can also be in play, for instance along the edges of the red outflow from MM2 and next to the outflow from MM1. On the other hand, we also detected narrow SiO emission offset from systemic velocities, which could be an indication of a specific environment where the Si is in grain cores. This is the case of CygX-N12, where there are narrow line profiles (FWHM ∼ 1 km s^-1) associated with the east-bound outflow (see Fig. A.2) and slightly blueshifted from the systemic velocities.

4.2.2. Other narrow emission

Despite the existence of SiO narrow emission associated with outflows, we also detect some narrow SiO emission, at systemic velocities, that does not have any broad outflow emission in the immediacies. This is our configuration (4).

This is most striking for CygX-N40 where most of the SiO emission is narrow, extended, and offset (by ≳0.1 pc) from the forming protostars, hence definitely not associated with outflows (≲0.05 pc) nor any hot-core (<0.005 pc). For this specific region, we have also investigated whether this extended SiO emission could be due to external radiation, but this seems unlikely as we have found no evidence of higher dust temperatures surrounding the N40 core (see e.g. Hennemann et al. 2012). It is perhaps worth noting that this region sits close to the convergence point between the DR21 ridge and another filament thought to be feeding material into the main ridge (Schneider et al. 2010). The shocks from this infall of material into the MDC could explain the existence of extended SiO emission (in line with what had been interpreted for other regions with extended narrow SiO emission with FWHM of 0.8−2 km s^-1, e.g. by Jiménez-Serra et al. 2010; Nguyên-Lu’o’ng et al. 2013; Sanhueza et al. 2013).

In CygX-N3 and N48, even though the narrow low-velocity emission has a significant overlap/adjacency to outflow regions, the SiO narrow emission has a similar spatial distribution to the interface regions of the small-scale converging flows observed here (Csengeri et al. 2011a,b), suggesting that the SiO could also have been formed from the low-velocity shocks of such flows. Interestingly, in both cases, the SiO narrow emission has a velocity gradient which does not mimic the velocity shears seen in high-density tracers (H¹³CO⁺ and N₂H⁺, Csengeri et al. 2011a,b). The mismatched velocity field could be explained by the fact that the SiO is not formed on the flows themselves, but as part of the post-shock material. If this post-shock material is originally at systemic velocities, and is dense enough to retain its original velocity field despite the external flow, the SiO would retain the original velocity of the (denser) local gas.

In Cyg-X N48, Csengeri et al. (2011b) detected CH₃CN as well, and suggest that it is tracing the gas shocked by such converging low-velocity flows. In Fig. 9 we show how the narrow SiO emission (in contours) correlates with this CH₃CN emission. Some of the CH₃CN emission could be associated with high-velocity outflows as they peak in regions adjacent to outflow lobes (e.g. the peak next to the outflow from IS-1). Nevertheless, it is interesting to note that the interface between the two converging flows (marked with a blue line) coincides with the most massive protostar of this MDC, and could be responsible for the high star formation activity in this region. This interface has both CH₃CN and narrow SiO emission, supporting the idea that they could both potentially arise not only from the several high-velocity outflows, but also from the shocks from small-scale converging flows.

In CygX-N12, CygX-N53, and CygX-N63, apart from the narrow emission that is associated with outflows (either blueshifted or in outflow cavity walls), there are small extensions with narrow emission at systemic velocities where no outflows co-exist. This SiO narrow emission extensions in CygX-N53, N63, and N12 could, similarly to N3 and N48, be pinpointing the convergence zones of the global collapse of the MDC material (or perhaps even the collapse of larger scale cloud structures) onto the protostellar envelopes.

Fig. 9 Map of the CygX-N48 integrated intensity of CH3CN from (Csengeri et al. 2011b) in greyscale, overlaid with the SiO systemic velocity emission from this paper in contours (for both molecules the velocity integration range was from −6.5 to −0.5 km s-1). The red ellipses show the regions with narrow SiO emission adjacent to the clear velocity shear found in N2H+, shown here as a blue line. Even though there are strong and compact SiO peaks coincident with high-velocity emission as seen in Fig. 5, the weaker and more extended SiO systemic velocity emission could be associated with the lower-velocity shocks from this velocity shear.

5. Discussion and conclusions

5.1. SiO as a probe of outflow shocks

Our study confirms that SiO is a good indicator of outflows in massive star-forming regions. However, we advise caution when using the SiO luminosity as a proxy of outflow power when individual outflows cannot be resolved. Doing so in our sample of sources would have induced an overestimation of the outflow power which could reach more than an order of magnitude (e.g. in CygX-N3, N40, and N48). This is mostly due either to the existence of a narrow component of SiO emission, often extended, which may not be due to outflow shocks. Nevertheless, the stronger the outflows, the more the outflow SiO luminosity will dominate over any other contributing SiO emission, resulting in a tighter correlation with outflow power.

These problems are attenuated when we are able to separate the individual outflows, where we tentatively recover a linear relation between F_CO and the total L_SiO, but because SiO needs to be formed and entrained by the outflows (unlike CO that needs only to be entrained), the SiO outflow emission depends not only on the efficiency of the outflows in forming and entraining SiO (which is likely linked to the outflow power and, as such, the evolutionary stage of the protostar), but it also depends on the amount of material encountered and shocked by the outflows along their propagating axis as they pierce the cloud. The scatter is therefore more significant for the relation of SiO luminosity with envelope mass, than it is when using the CO outflow momentum fluxes.

5.2. Origin of the narrow SiO emission

From our sample of sources, we have realised that MDCs forming the least massive stars have their SiO luminosity dominated by a narrow component of SiO at systemic velocities (see Fig. 8). This narrow component, however, does seem to exist in all cases. While for the most massive cores it seems mostly related to outflow walls or compact narrow emission close to the protostellar envelopes (at systemic velocities), the least massive sources have a tendency to show more extended SiO narrow emission. Such narrow SiO emission could have various origins, and in the following we shall explore some of the possibilities.

One hypothesis is that it arises from the mixing of the gas due to the interaction of different outflows from a population of low-mass stars undetected in the PdBI continuum images. We note, however, that the sensitivity in the 1mm continuum observations is equivalent to ∼1 M_⊙, and that all the potential low-mass protostars with envelope masses above ∼1 M_⊙ are listed by Bontemps et al. (2010) and shown as crosses in our figures. The positions of these candidate low-mass protostars are not associated with any CO high-velocity outflow (which we would expect to detect as we did in CygX-N40), and they also do not lie near most of the narrow emission detected. Judging by the number of low-mass fragments detected we do not expect a significant underlying population of very low-mass protostars (<1 M_⊙). Even though we cannot rule out this hypothesis completely, we consider that this scenario is thus unlikely.

Another hypothesis is that the shocks that may create the narrow SiO emission where no nearby high-velocity shocks exist could be related to converging flows shocking with the cloud’s material. In our sample, these shocks could be due to the small-scale converging flows reported by Csengeri et al. (2011a,b) or due to the large scale collapse of the DR21 ridge as suggested by Schneider et al. (2010). This is in line with what had been suggested for other high-mass star-forming regions by e.g. Jiménez-Serra et al. (2010); Nguyên-Lu’o’ng et al. (2013). Perhaps seeing a dominant SiO narrow emission is actually a sign that there is a significant and ongoing infall of material from beyond the MDC scale (≳0.1−0.3 pc) to the envelope scale (≲0.02 pc), which could lead to a further increase of the masses of existing individual cores or to the formation of new ones. The most striking case of our sample is undoubtedly CygX-N40 where the inner envelope is not yet massive at all, but there are SiO shocks everywhere around it. We note that in spite of comprising nearly ∼100 M_⊙ inside less than 0.15 pc (Motte et al. 2007), Cygx-N40 has a rather shallow density profile, with no high-mass protostars yet formed in the centre. Perhaps N40 is at an earlier stage than the other sources, where the material of the MDC is yet to become centrally condensed and form one or more high-mass protostellar cores, and the SiO emission is the evidence of the shocks of the infall of large-scale flows and the collapsing MDC itself.

One of the problems with having SiO at low velocities not associated with outflows is to understand if and how low-velocity shocks could form SiO. As mentioned by Nguyên-Lu’o’ng et al. (2013), SiO can be efficiently formed by low-velocity shocks (<10 km s^-1) provided that there is already some Si in the gas phase or in ice mantles. Having Si surviving in the gas phase in the cold and dense environment of an MDC is, however, not very likely since Si is expected to either deplete back onto the grain ice mantles, or to be oxidised in SiO₂, itself being subsequently depleted back onto the grain ice mantles, both processes leading to its disappearance from the gas phase in about ∼10⁴ yr (e.g. Martin-Pintado et al. 1992; Codella et al. 1999; Gusdorf et al. 2008a). Nevertheless, it is possible to admit that the converging flows are feeding the MDCs with pre-shock material originating from the diffuse medium, which contains Si in the gas phase due to the exposition to the radiation field of nearby OB stars, or even that in some cases there is some SiO released into the gas phase due to shocks from earlier high-velocity outflows. This SiO could then easily be maintained in the gas phase by lower velocity shocks. If this is the case, shocks with velocities below 5 km s^-1 could start producing enough SiO to be detected. We retrieve SiO column densities of the order of 0.5−1.0 × 10¹³ cm^-2 (see Sect. 3.4) for the narrow SiO emission (similar to what was found by Nguyên-Lu’o’ng et al. 2013), which is consistent with the possibility that SiO has been formed with low-velocity shocks. The low-velocity shears detected in the Cygnus-X MDCs have a line of sight velocity difference around 2−3 km s^-1 (Csengeri et al. 2011a,b). Although this lies in the very lower limits to produce any SiO at all, these are lower limits for the actual velocity shears because projection effects are most certainly in play. Furthermore, and as mentioned in Sect. 4.2.2, since SiO molecules form in the post-shock gas, the actual velocities of SiO and the respective velocity dispersions can be different from those of the shocks, supporting the idea that the low-velocity converging flows in the region (Csengeri et al. 2011a,b) could be responsible for the widespread SiO narrow emission, despite the different observed velocities.

We therefore propose an evolutionary picture for the existence of SiO at systemic velocities and with a narrow line profile, not associated with outflows. In this scenario, the least centrally condensed MDCs are interpreted as being in an earlier stage of the collapse and the SiO would be tracing shocks from the global collapse of cloud onto the MDC. This would be the case for CygX-N40, which is a MDC with only a couple of low-mass protostars (<2 M_⊙) at its centre, a very weak outflow, and where 90% of the SiO emission arises from the outskirts. This could also be the case of, for instance, IRDC G028.23-00.19 (Sanhueza et al. 2013), but higher resolution observations would be necessary to confirm this. As the MDC collapses, the SiO emission becomes more confined to the close surroundings of denser parts, tracing the post-shock material from the shocks of the infalling MDC directly against the dense central cores. This is in essence the equivalent of the shocks from the small-scale low-velocity converging flows seen in CygX-N3, N12, and N48. At later stages, the MDC infall is less prominent and a more significant fraction of the MDC mass is already gathered in the central massive protostars. At this stage, the SiO luminosity is largely dominated by the powerful outflows, and there is only a weak narrow component which would point to the narrow SiO emission arising from the shocks from the very last remnants of the collapse of the MDC against the protostellar envelopes – this could be the case of CygX-N53 and N63. In this scenario, we would not expect to see a great amount of extended narrow SiO emission at systemic velocities in low-mass star-forming regions (apart from the narrow SiO emission intrinsically originating from high-velocity shocks), because the simple gravitational collapse of the clumps is unlikely to create shocks with enough speed to produce SiO (infall velocities being typically ∼1 km s^-1), and these regions are generally less dynamic and less dense environments, less prone to shocks in the ISM capable of forming SiO. This is consistent with the fact that the SiO emission detected so far in low-mass star-forming regions is associated with high-velocity protostellar outflows, either as broad SiO emission, or as narrow emission with a velocity peak offset from systemic velocities (e.g. Gueth et al. 1998; Nisini et al. 2007; Gusdorf et al. 2011; Gómez-Ruiz et al. 2013).

Even though this dynamical view would predict a smaller amount of widespread SiO narrow emission at systemic velocities in low-mass star-forming regions, it is sill puzzling to understand why we do not seem to observe any. If low-velocity shocks are able to form SiO (whether from outflows or other turbulent shocks in the ISM), then they should also be present in the low-mass case. One could think of two ways to explain the lack of detection of this component: either the composition of the gas is different in massive star-forming regions (e.g. with a greater amount of Si on grain mantles or in the gas phase), making it easier to form SiO with low-velocity shocks, or there is a density threshold where the lower velocity shocks become efficient at forming SiO only if the volume densities are sufficiently high. Hence it would be interesting to test the minimum densities needed to be able to form SiO through low-velocity shocks. This will be addressed in upcoming papers from Gusdorf et al. and Louvet et al. (in prep.). Observationally, progress can be made by observing higher energy transitions of SiO to help constrain all the characteristics of the shocks (velocities, magnetic field, pre-shock densities, chemistry, UV field, etc.).

Acknowledgments

A.D.C. acknowledges funding from the European Research Council for the FP7 ERC starting grant project LOCALSTAR. Part of the work by A.D.C. was accomplished under funding by the project PROBeS, and N.S. is supported by the project STARFICH, both funded by the French National Research Agency (ANR). T.C. acknowledges financial support for the ERC Advanced Grant GLOSTAR under contract No. 247078. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain). The data reduction and analysis made use of the GILDAS software (http://www.iram.fr/IRAMFR/GILDAS).

References

Appendix A: Examples of the SiO spectra

Fig. A.1 Procedure used to remove the broad high-velocity SiO emission from the datacubes. Top: example of the line fitting of a spectrum where both broad and narrow emission exist (corresponding to the blue spectrum of CygX-N3 shown in Fig. A.2). The two Gaussian fits are shown in green (for the broad component) and blue (for the narrow component). The red line shows the sum of the two Gaussian components. Bottom: same spectrum as the top panel, here with the respective broad Gaussian fit subtracted.

Fig. A.2 SiO emission (in units of Tmb) of the six MDCs studied. For each MDC we show, on the left, a greyscale showing the total integrated SiO emission and, on the right, three overlaid spectra extracted from the datacubes. The grey, blue, and red circles in the left-hand panels show the positions of the grey, blue, and red spectra displayed on the right. The grey spectra exemplify the narrow profile of the SiO emission, and the blue and red spectra are taken at positions where there is significant wing emission. Sources and outflow directions are as labelled in Figs. 2 and 3.

All Tables

All Figures

Fig. 1 Average spectra for the six MDCs of our sample (labelled on the top-left corners). The left side of each subfigure shows the comparison of the single-dish (SD) IRAM 30 m pointed observations (Motte et al. 2007) in blue dashed lines, and the spectra of our PdBI SiO observations integrated over the central 30′′ × 30′′ (equivalent to the IRAM 30 m beam) are shadowed in grey. The vertical black dashed lines show the systemic velocity of each core, as measured with N2H+ by Bontemps et al. (2010). The right side shows the same PdBI SiO emission as the left, overplotted with a two-component Gaussian fit to illustrate the existence of narrow (in blue) and broad (in green) components of emission (whose sum is shown in red). For CygX-N12 and N40, we show a shadowed spectrum (in light grey) which corresponds to the PdBI emission integrated over an area of 50′′ × 50′′ (to assess a possible contamination from strong SiO emission outside the SD primary beam).

In the text

Fig. 2 Comparison of the high-velocity CO and SiO emission in CygX-N3 (left), N12 (centre), and N40 (right). The top panels show the high-velocity CO emission in contours (with the velocity ranges as in Duarte-Cabral et al. 2013) over the 1 mm continuum emission in greyscale. The primary beam of the CO and 1mm observations is shown as a dashed circle, and the synthetic beam is in the lower-left corner. Contours are stepped by 0.25 Jy beam-1 km s-1, starting at 2.0 (2.75) Jy beam-1 km s-1 for N3 blue (red), 1.5 Jy beam-1 km s-1 for N12, and 1.75 (1.25) Jy beam-1 km s-1 for N40 blue (red). The lower panels show the high-velocity SiO emission in contours (with velocity ranges as in Table 2) over the 3 mm continuum emission in greyscale for the same sources. The primary beams are shown as dashed circles and the synthetic beams are plotted in the lower-left corners. Contours are stepped by 0.1 Jy beam-1 km s-1, starting at 0.15 Jy beam-1 km s-1 for N3, 0.2 (0.05) Jy beam-1 km s-1 for N12 blue (red), and 0.2 Jy beam-1 km s-1 for N40. The most massive protostars are marked with yellow stars (and labelled in the lower panels), and the least massive fragments from Bontemps et al. (2010) are shown as green crosses. The arrows are the directions of the outflows inferred from Duarte-Cabral et al. (2013), and the dashed lines show the respective outflow cones.

In the text

Fig. 3 Same as Fig. 2 for CygX-N48 (left), N53 (centre), and N63 (right). In N48 we mark with circled crosses two Spitzer sources: ERO2 (Marston et al. 2004) and IS-1, which is a resolved source detected in the IRAC bands at 3 and 4 μm, and unresolved at 8 μm and onwards. IS-1 is likely responsible for the outflow directed east-west. Contour key for CO emission, top panels: N48 (left) contour steps of 0.5 Jy beam-1 km s-1 starting at 4 (3) Jy beam-1 km s-1 for blue (red); N53 (centre) contour steps of 1.5 Jy beam-1 km s-1 starting at 8.5 Jy beam-1 km s-1; N63 (right) contour steps of 1.5 Jy beam-1 km s-1 starting at 5.0 Jy beam-1 km s-1. Contours of SiO emission, lower panels: N48 (left) contour steps of 0.1 Jy beam-1 km s-1, starting at 0.2 (0.15) Jy beam-1 km s-1 for blue (red); N53 (centre) contour steps of 0.15 (0.25) Jy beam-1 km s-1, starting at 0.2 (0.15) Jy beam-1 km s-1 for blue (red); N63 (right) contour steps of 0.5 Jy beam-1 km s-1, starting at 0.2 Jy beam-1 km s-1, for both blue- and redshifted emission.

In the text

Fig. 4 Integrated intensity of the SiO emission at ambient velocities with the contribution from the broad outflow emission subtracted, for CygX-N3 (left), N12 (centre), and N40 (right). The top row panels show the narrow emission as greyscale, with contours from the high-velocity (broad) SiO emission from Fig. 2 overlaid. The two lower rows show the integrated narrow SiO emission in contours (with steps of 0.1 Jy beam-1 km s-1, starting at 0.1 Jy beam-1 km s-1), with the colour scale showing the respective moment maps: the velocity field (first moment map) in the middle row, and the velocity dispersion (second moment map) in the last row. Annotations are as in Fig. 2.

In the text

Fig. 5 Same as Fig. 4, for CygX-N48 (left), N53 (centre), and N63 (right). Top panels have the narrow SiO emission in greyscale and the high-velocity (broad) SiO emission in contours as in Fig. 3. The two lower rows show the integrated narrow SiO emission in contours (with steps of 0.1 Jy beam-1 km s-1, starting at 0.1 Jy beam-1 km s-1), with the respective moment maps in colour scale: the velocity field (first moment map) in the middle row, and the velocity dispersion (second moment map) in the last row. Annotations are as in Fig. 3.

In the text

Fig. 6 Observed LSiO against the corresponding FCO (from Duarte-Cabral et al. 2013). The light blue filled circles correspond to the total LSiO and total FCO for each of the six fields, and the blue dashed line is the respective linear fit in a log-scale. The open circles are for the resolved individual outflows (integration spatially constrained to location of the individual outflows), with the FCO estimated in Duarte-Cabral et al. (2013; Table 3, Col. 6). The red open circles have LSiO estimated using the entire velocity range (Table 3, Col. 7), and the respective linear fit to the log-values is shown as a dashed red line, which translates into a linear correlation of LSiO ∝ FCO. The green open circles show the LSiO only from the outflow high-velocity wing emission (Table 3, Col. 8).

In the text

Fig. 7 Observed LSiO (same as the red open circles in Fig. 6) and FCO, in blue and red circles, respectively, with respect to the envelope masses (Menv) for individual sources. The second red point from the left corresponds to CygX-N40 MM1, for which we have an FCO measurement, but no LSiO. The dashed red line shows the linear correlation of FCO = 7.5 × 10-5Menv, originally found by Bontemps et al. (1996) for low-mass protostars, and confirmed in Duarte-Cabral et al. (2013) to extend to the high-mass regime. The dashed blue line shows the linear correlation of LSiO ≃ 5 × 106Menv that we find for the blue points. This figure accentuates the fact that LSiO of individual outflows follows the same trend as FCO (both being ∝Menv), even though it introduces more scatter with respect to the linear relation found for FCO.

In the text

Fig. 8 Fraction of the total SiO luminosity that arises from the narrow component, fnarrow, as a function of the mass of the most massive protostellar envelope which is currently forming in each MDC. This plot accentuates the trend of a higher fraction of the emission in the form of narrow emission for the MDCs which have not yet gathered most of their mass in the envelopes of the forming protostars.

In the text

Fig. 9 Map of the CygX-N48 integrated intensity of CH3CN from (Csengeri et al. 2011b) in greyscale, overlaid with the SiO systemic velocity emission from this paper in contours (for both molecules the velocity integration range was from −6.5 to −0.5 km s-1). The red ellipses show the regions with narrow SiO emission adjacent to the clear velocity shear found in N2H+, shown here as a blue line. Even though there are strong and compact SiO peaks coincident with high-velocity emission as seen in Fig. 5, the weaker and more extended SiO systemic velocity emission could be associated with the lower-velocity shocks from this velocity shear.

In the text

Fig. A.1 Procedure used to remove the broad high-velocity SiO emission from the datacubes. Top: example of the line fitting of a spectrum where both broad and narrow emission exist (corresponding to the blue spectrum of CygX-N3 shown in Fig. A.2). The two Gaussian fits are shown in green (for the broad component) and blue (for the narrow component). The red line shows the sum of the two Gaussian components. Bottom: same spectrum as the top panel, here with the respective broad Gaussian fit subtracted.

In the text

Fig. A.2 SiO emission (in units of Tmb) of the six MDCs studied. For each MDC we show, on the left, a greyscale showing the total integrated SiO emission and, on the right, three overlaid spectra extracted from the datacubes. The grey, blue, and red circles in the left-hand panels show the positions of the grey, blue, and red spectra displayed on the right. The grey spectra exemplify the narrow profile of the SiO emission, and the blue and red spectra are taken at positions where there is significant wing emission. Sources and outflow directions are as labelled in Figs. 2 and 3.

In the text