Free Access
Issue
A&A
Volume 569, September 2014
Article Number A11
Number of page(s) 24
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201424032
Published online 09 September 2014

© ESO, 2014

1. Introduction

High-mass stars (O- and B-type stars with masses 8  M) are crucial for the understanding of many physical phenomena in the Galaxy; however, their first stages of formation are still poorly understood and represent a challenge from both a theoretical and observational point of view. A number of different theoretical scenarios have been proposed to explain the formation of OB-type stars (see Zinnecker & Yorke 2007 and Tan et al. 2014 for a review). Even with their differences, the most accepted theories (monolithic collapse in a turbulence-dominated core – Krumholz et al. 2009; competitive accretion driven by a stellar cluster – Bonnell & Bate 2006; Bondi-Hoyle accretion – Keto 2007) agree in the prediction of the formation of circumstellar disks around stars of all masses. However, while in the low-mass regime, circumstellar disks have been extensively identified and studied (e.g., Simon et al. 2000), there are only a few convincing disks (with sizes 1000 AU) found around B-type protostars (e.g., Schreyer et al. 2002; Cesaroni et al. 2005, 2007; Carrasco-González et al. 2012; Wang et al. 2012). Moreover, convincing evidence of centrifugally supported disks around O-type protostars still remains elusive. Infrared interferometric observations of the CO overtone bands suggest that disks might be associated with O-type stars (e.g., Kraus et al. 2010; de Wit et al. 2011; Boley et al. 2013), but (sub)millimeter images have only revealed the existence of huge (~0.1 pc) and massive (~100  M) cores undergoing solid-body rotation, the so-called “toroids” (e.g., Beltrán et al. 2011). The detection of circumstellar disks, with expected sizes 1000 AU (e.g., Cesaroni et al. 2007), requires angular resolutions < 0".1 (assuming a typical distance of 5 kpc). Thus, the lack of observed disks is expected to dramatically change with the advent of the Atacama Large Millimeter Array (ALMA), with which resolutions will be easily obtained.

thumbnail Fig. 1

a)Spitzer 4.5 μm image (in logarithmic scale) of the star forming region G35.200.74 N, overlaid with a contour map of the 870 μm (350 GHz) continuum emission obtained with ALMA. The IR image has been obtained by applying HiRes deconvolution (Velusamy et al. 2008) to the Spitzer/IRAC data. The submillimeter map has been corrected for primary beam attenuation. Contours start at 5σ, increasing in steps of 6σ (with σ = 1.8 mJy beam-1), up to a maximum level of 201.7 mJy beam-1. The dashed white circle denotes the primary beam of the ALMA 12 m antennas (~9′′), while the ALMA synthesized beam (, PA = 46°) is shown in the bottom-right corner of the panel. b) Close-up of the central region that shows the 870 μm continuum emission map obtained with ALMA. Six cores (from A to F) have been identified and marked in the image (see Table 2). Yellow crosses mark the position of OH masers (Hutawarakorn & Cohen 1999). The ALMA synthesized beam is shown in the bottom-right corner of the panel. The spatial scale is also indicated in both panels.

Open with DEXTER

In a progressive effort to confront the theoretical scenario of the formation of massive stars through disk-mediated accretion, and to set the stage for a follow-up quest for the more challenging disks around O-type stars, we carried out ALMA Cycle 0 observations towards two IR disk candidates around B-type protostars in G35.200.74 N and G35.03+0.35. Results on G35.03+0.35 are presented in a parallel paper by Beltrán et al. (2014), while first results on G35.200.74 N have been already published in Sánchez-Monge et al. (2013b). In this work, we present the ALMA results on G35.200.74 N in a more extensive and detailed way.

G35.200.74 is a star forming complex located at a distance of kpc (Zhang et al. 2009). The main site of high-mass star formation activity, known as G35.200.74 N (hereafter G35.20N), is associated with the IRAS source 18556+0136 and has a bolometric luminosity ~1–10 × 104 L (e.g., Gibb et al. 2003; Sánchez-Monge et al. 2013b; Zhang et al. 2013). In the infrared, the region is dominated by the presence of a butterfly-shaped nebula oriented NE–SW (see Fig. 1a). This direction is coincident with the orientation (PA58°) of a bipolar outflow observed in 12CO (e.g., Dent et al. 1985a; Gibb et al. 2003; Birks et al. 2006; López-Sepulcre et al. 2009; Qiu et al. 2013). The 12CO (1–0) line emission appears to trace also a N–S collimated flow (see Birks et al. 2006), coinciding with a thermal radio jet (e.g., Heaton & Little 1988; Gibb et al. 2003) seen also at IR wavelengths (e.g., Dent et al. 1985b; Walther et al. 1990; Fuller et al. 2001; De Buizer 2006; Zhang et al. 2013). The different orientations between the poorly collimated NE–SW outflow and the N–S jet have been interpreted as manifestations of the same flow undergoing precession (e.g., Little et al. 1998; Gibb et al. 2003), or as emission coming from a number of flows driven by the different sources forming in the G35.20N complex (e.g., Gibb et al. 2003; Birks et al. 2006).

Emission from dense gas tracers (e.g., NH3, C18O, CS, HCN, H13CO+), as well as continuum dust emission at millimeter wavelengths (e.g., Little et al. 1985; Brebner et al. 1987; Gibb et al. 2003; López-Sepulcre et al. 2009), have revealed the presence of a dense clump elongated perpendicular to the NE–SW outflow, with a velocity gradient from NW (~32 km s-1) to SE (~36 km s-1). This velocity gradient was first thought to originate from a large (~0.6 pc) flattened structure rotating about the axis of the NE–SW outflow. Subsequently, on the basis of their H13CO+ and H13CN observations, Gibb et al. (2003) propose that this is actually a fragmented rotating envelope containing multiple young stellar objects. The high angular resolution map at 870 μm obtained with ALMA (Sánchez-Monge et al. 2013b) reveals an elongated structure containing a chain of compact cores (see Fig. 1b), and confirms a significant level of clumpiness. The two strongest cores (named core A and core B in Sánchez-Monge et al. 2013b) are found close to the center of the elongated structure (core A to the NW and core B to the SE, separated by ~2′′). These cores, detected in several complex organic molecules, show velocity gradients roughly aligned with the elongated structure, i.e., perpendicular to the large-scale NE–SW outflow, but opposite to each other. The kinematics of core B, which is located at the center of symmetry of the outflow/jet, have been modeled as a Keplerian disk rotating about a central mass of ~18  M (Sánchez-Monge et al. 2013b).

Finally, emission from different maser species (H2O, OH, CH3OH) has been detected toward the G35.20N region (e.g., Brebner et al. 1987; Forster & Caswell 1989, 1999; Hutawarakorn & Cohen 1999; Vlemmings 2008; Sugiyama et al. 2008; Surcis et al. 2012). The highest angular-resolution observations show two main regions for the maser emission, separated by about (~5000 AU), and oriented SE-NW. Hutawarakorn & Cohen (1999) measured magnetic field strengths (through 1665 MHz OH Zeeman-splitting observations) between 2.5 mG to the NW, and +5.2 mG to the SE. Polarimetric observations of 6.7 GHz CH3OH masers (Surcis et al. 2012), indicate strong (5–12%) linear polarization, with the magnetic field oriented almost E–W, PA80–110°. Recent polarimetric dust observations carried out with the SMA (Qiu et al. 2013), revealed that the magnetic field is aligned with the long axis of the elongated dust structure to the NW roughly following the large-scale field (PA56°; Vallée & Bastien 2000), but changes its direction drastically by approximately 90° when moving to the south of the filament (at a position between cores B and C; see Fig. 2 in Qiu et al. 2013). The field strength is estimated to be ~1 mG, i.e., consistent with the strength measured via Zeeman-splitting observations (Hutawarakorn & Cohen 1999).

Given the complexity of this region and the (sometimes contradictory) findings of previous observations, the main goal of our ALMA observations was to address some fundamental questions such as the fragmentation level of the elongated structure, the kinematics of the dense gas at different scales and the number of outflows and powering sources. This paper is organized as follows. In Sect. 2, we describe our ALMA observations. In Sect. 3, we present the main results and reconstruct the spectral energy distribution of G35.20N from the millimeter to the near IR. In Sect. 4, we analyze and discuss our findings focusing on the kinematical properties of the structures detected. Finally, in Sect. 5 we summarize the main results and draw the conclusions.

Table 1

Molecular transitions studied in this paper.

2. ALMA observations

G35.20N was observed with ALMA in Cycle 0 between May and June 2012. The source was observed in Band 7 (~350 GHz) with the 16 antennas of the array in the extended configuration, i.e., baselines in the range 36–400 m, providing sensitivity to structures 2′′. The digital correlator was configured in four spectral windows (with dual polarization) of 1875 MHz and 3840 channels each, providing a resolution of ~0.4 km s-1. The four spectral windows covered the frequency ranges [336 849.57–338 723.83] MHz, [334 965.73–336 839.99] MHz, [348 843.78–350 718.05] MHz, and [346 891.29–348 765.56] MHz. The phase center of the observations is α(J2000) = 18h58m13, δ(J2000) = 01°40360. Flux, gain, and bandpass calibrations were obtained through observations of Neptune and J1751+096. The data were calibrated and imaged using CASA1 (McMullin et al. 2007). A continuum map was obtained from line-free channels and subtracted from the data, directly in the (u,v)-domain. Channel maps were created for some selected molecular transitions (see Table 1 and Sect. 3). Continuum and channel maps were created with the robust parameter of Brigg (1995) set equal to 0.5, as a compromise between resolution and sensitivity to extended sources. The resulting continuum image has a synthesized CLEANed beam of , PA = 46° and a rms noise of 1.8 mJy beam-1. The synthesized CLEANed beams of the spectral line maps vary from to with increasing frequency, with PA = 40–45°. The rms noise of each spectral channel (of 0.6 km s-1) varies between 5 and 20 mJy beam-1, with the largest values measured in channels with strong line emission. The spectral line images were analyzed using the CLASS and GREG programs of the GILDAS2 software package developed by the IRAM and Observatoire de Grenoble. All the resulting spectra have been smoothed to a resolution of 0.6 km s-1.

3. Results

3.1. Continuum emission

In Fig. 1a, we show the map of the ALMA 870 μm (350 GHz) continuum emission towards G35.20N overlaid on the Spitzer 4.5 μm image, whose resolution has been enhanced with HiRes (Velusamy et al. 2008). The submillimeter continuum emission appears to trace an elongated structure in the SE-NW direction, perpendicular to the 4.5 μm extended emission, and located across the waist of the butterfly-shaped nebula seen in the infrared. With an extension of ~0.15 pc and a mean width of ~0.013 pc (full width at half maximum), this is probably the densest part of the elongated structure observed at larger scales (0.21 pc long and 0.13 pc wide; Little et al. 1985; Brebner et al. 1987; Gibb et al. 2003). This structure is not smooth but fragmented in several smaller condensations. The 2D CLUMPFIND algorithm (Williams et al. 1994) identifies six of these, which we name cores A–F (see Fig. 1b), although additional substructure is seen in the map. The strongest of these, cores A and B, are located at the center of the elongated structure, and were already identified by Sánchez-Monge et al. (2013b). These two cores are associated with faint centimeter continuum emission likely tracing young H ii regions or radiojets (see Gibb et al. 2003; Sánchez-Monge et al. 2013b), with a number of 1665-MHz OH maser features (Brebner et al. 1987; Hutawarakorn & Cohen 1999), and with probably3 6.7 GHz CH3OH maser emission. Core C, located ~3′′ to the SE of core B, coincides with the millimeter source G35MM2 reported by Gibb et al. (2003). The ALMA 870 μm dust continuum emission shown in Fig. 1a is consistent with the Submillimeter Array (SMA) 880 μm image obtained by Qiu et al. (2013, see their Fig. 1c). There are some morphological differences to the southeast of the elongated structure, probably due to the vicinity of the primary beam edge in our ALMA observations.

In Table 2, we list for all the cores the coordinates of the peak position, the primary beam corrected fluxes, and the deconvolved sizes. The fluxes have been computed by integrating the emission within the 5σ contour level. The peak intensities range from 40 mJy beam-1 to 200 mJy beam-1, while the integrated fluxes are in the range 80–700 mJy. These fluxes are slightly lower (a factor of 2) than those measured by Qiu et al. (2013) in their SMA observations. However, we note that Qiu et al. determine the flux densities for larger dust entities containing at least two cores (e.g., their object MM1, with a total flux of 3 Jy, is the sum of cores A and B). We have convolved the ALMA image to the beam of the SCUBA 850 μm image (15′′; Qiu et al. 2013), and compared the peak fluxes of both images (~2.5 Jy beam-1 versus ~8.5 Jy beam-1) to conclude that we are recovering around 30% of the flux measured with SCUBA. The deconvolved diameters of the dust cores (see Table 2) have been computed from the 50% contour level applying Gaussian deconvolution (see procedure described in Table 3 of Sánchez-Monge et al. 2013a). The mean deconvolved diameter is ~1600 AU, with core A being the most compact (900 AU), and core F the most extended (2600 AU). The different cores appear regularly spaced, with a mean (projected) separation of ~ or ~0.023 pc.

3.2. Line emission

In Fig. 2, we show the spectra obtained by averaging the line emission over the cores in G35.20N (lower sideband: 334 965.73–338 723.83 MHz, and upper sideband: 346 891.29–350 718.05 MHz). Cores A and B show rich spectra, with emission of several complex organic molecules (e.g., CH3CN, CH3OH, CH3OCHO, NH2CHO, HC3N, C2H5OH, CH3OCHO). Out of the rest of the cores, only core C is likely associated with emission in complex molecules such as CH3CN and CH3OH. Core C was found to be associated with molecular emission at millimeter wavelengths in previous studies (e.g., Gibb et al. 2003; Birks et al. 2006). The rest of the cores are only associated with faint emission of more abundant species such as C17O, C34S, H13CO+ (see Fig. 3) and weak CH3OH, suggesting that these cores could either host only low-mass star formation, or be in an evolutionary phase prior to the hot core stage.

3.2.1. Large-scale molecular emission

Out of all the species detected in the region (mainly associated with cores A and B; see Sect. 3.2.2), there are four molecules (C17O, C34S, H13CO+ and SiO) that show a complex (extended) spatial distribution. Their emission appears partially resolved out by the interferometer, which is sensitive to angular scales 2′′ (see Sect. 2). In order to improve the detection of the extended emission, we generated new CLEANed maps for these four species after tapering the data at 1′′ resolution. In Fig. 3, we show the corresponding integrated emission (zeroth-order moment) maps, over the velocity range 23 to 42 km s-1, and compare them with the dust continuum emission. For these molecules, the line emission coincides with the dust continuum emission, with C17O and H13CO+ being detected toward all the cores and tracing the elongated structure, and SiO and C34S mostly concentrated towards cores A and B, consistent with the SMA images reported by Qiu et al. (2013). We estimate a missing flux of 50% from the comparison of our ALMA H13CO+ image and the map obtained by Qiu et al. using compact configurations of the SMA. Such an elongated structure has also been revealed at larger scales (~40′′) in the H13CO+ (1–0) and H13CN (1–0) lines (e.g., Gibb et al. 2003). In Sect. 4.5, we will analyze the kinematics of this extended molecular structure.

Table 2

Parameters of the cores detected in the 870 μm ALMA continuum map (see Fig. 1).

thumbnail Fig. 2

Spectra of the full frequency range surveyed by the ALMA observations, toward the six cores identified in the continuum image (see Fig. 1). The spectra are extracted by averaging the intensity over an area of 1.37, 2.35, 1.04, 0.55, 0.75, 0.44 arcsec2, yielding conversion factors of 8.9, 5.0, 11.3, 21.7, 15.7, 27.1 K Jy-1, for cores A to F respectively. The brightness temperature scale of cores A and B is different from that of the rest of the cores. The absorption seen in some spectra is likely produced by the filtering of extended emission.

Open with DEXTER

3.2.2. Compact molecular emission: hot cores

As shown in Fig. 2, cores A and B show a very rich chemistry with strong emission of many species typical of hot cores. Core C also shows emission in complex molecules, although with an intensity several times fainter than the two cores located at the center of the elongated structure. A detailed analysis of the molecular content towards these three cores will be the subject of a forthcoming paper. In the present study, we will only focus on a few hot core tracers, basically those that are clearly detected, and less blended with other species: CH3CN, CH313CN, and vibrationally excited CH3CN, as well as CH3OH, CH3OCHO, HC3N, HNCO, H2CS, and SO2. Similar to the large-scale molecular and continuum emission, we have estimated the missing flux in our maps of hot-core tracers by comparing the ALMA CH3CN K = 0,1 lines with the same lines observed with the SMA in compact configurations (Qiu et al. 2013) and with the Atacama Submillimeter Telescope Experiment (ASTE; Furuya, priv. comm.). The ALMA observations recover around 75% of the flux measured with ASTE. The fluxes between the SMA and ALMA images (after convolution to the same beam) are equal within a 8%.

Figures B.1B.3 show a close-up view of the CH3CN and CH3OH transitions detected towards cores A, B and C, respectively. The observational setup (see Sect. 2) permits to cover up to the K = 9 transition of CH3CN (19K–18K), and up to the K = 5 line of CH313CN (19K–18K), with upper level energies of 745 K and 517 K, respectively. All of these CH3CN lines have been detected in the three cores. Regarding the CH313CN, only a few lines are not blended (e.g., K = 2, K = 4), while others clearly overlap with methyl cyanide or methanol lines (e.g., K = 3, K = 5). Several vibrationally excited CH3CN transitions, with upper level energies ranging from ~600 K to ~1100 K, have been detected. The intensity of the methyl cyanide lines in core A is almost twice the intensity of the lines in core B, and ten times stronger than in core C. A large number of methanol (CH3OH) and methyl formate (CH3OCHO) lines have been detected in these three cores, but in the following, we have only considered those transitions that are not blended with other (strong) lines. The selected lines are listed in Table 1, and cover a range of upper level energies from 45 K to 488 K in the case of methanol, and from 236 K to 444 K for methyl formate. Finally, we have also selected some strong and isolated lines of different species such as H2CS, HC3N, SO2 or HNCO.

In Figs. 4, A.1 and A.2, we present the zeroth (integrated intensity) and first (velocity field) moment maps of different hot core species detected towards cores A and B: CH3CN K = 2 and 8, CH313CN K = 2, and CH3CN v8 = 1K,l = 3,1, as well as different transitions of CH3OH, CH3OCHO, H2CS, HNCO, SO2, and HC3N. The moment maps have been computed following the same procedure as in Sánchez-Monge et al. (2013b), using the velocity interval 25.9–38.0 km s-1 for core A, and 25.9–33.0 km s-1 for core B. As seen in the maps, the line emission peak in core A is perfectly coincident with the position of the dust emission peak. In both the molecular and continuum images, the emission of core A is barely resolved with a mean deconvolved diameter of ~1100 AU (~), obtained from a 2D Gaussian fit to the zero-order moment maps, and applying Gaussian deconvolution. In core B, the strongest peak of the molecular emission is typically located to the SE of the dust peak. For most of the transitions, the molecular emission of core B is clearly resolved (deconvolved size ~; ~2000 AU) in an elongated structure oriented in the SE-NW direction, with the dust peak located close to the center of the extended line emission. For both cores, the emission is more compact in the maps of high excitation energy transitions (e.g., CH3CN v8 = 1, or CH3OCHO at 335.015 GHz) or less abundant species (e.g., CH313CN K = 2), suggesting the existence of temperature and/or molecular abundance gradients within the cores. This is better seen in Fig. 6, where we plot the deconvolved size for different transitions against the upper level energy of the transition. For both cores, the trend is such that the higher the excitation energy of the transition is, the more compact the emission is. From the velocity field maps, we see a clear velocity gradient in both cores for all the species. The velocity gradient of core B is oriented in the (S)SE-(N)NW direction, while the velocity gradient in core A is closer to the N–S direction. We also note that the velocity gradient in both cores is reversed, with redshifted emission to the north for core B and to the south for core A. A more detailed analysis of the kinematics of these two cores is presented in Sect. 4.3.

Finally, in Fig. 5, we present the integrated intensity and velocity field maps of four transitions of hot-core species towards core C: CH3CN (192–182), CH3OH (22,1–31,2), H2CS (101,9–91,8), and HC3N (37–36). The moments were computed using the velocity interval 26.1–41.0 km s-1. Similar to core A, the molecular emission in core C is compact and coincident with the dust emission. Differently from cores A and B, we do not find a clear velocity gradient in any of the molecules.

thumbnail Fig. 3

Overlay of the 870 μm continuum emission (red contour) on the integrated intensity maps (colors) of C17O, C34S, H13CO+, and SiO. The contour level of the continuum emission corresponds to the 9.0 mJy beam-1 intensity (see Fig. 1). The intensity scale for the integrated emission for each molecule is shown in the color bar at the right of each panel, in K km s-1. The synthesized beam of the continuum emission (, PA = 46°) is shown in the bottom-left corner of each panel. For each molecule, we show the synthesized beam of the tapered maps in the bottom-right corner: a) , PA = 67°; b) , PA = 67°; c) , PA = 69°; and d) , PA = 69°.

Open with DEXTER

thumbnail Fig. 4

Maps of the integrated intensity (contours) and velocity field (color maps in the bottom panels) of a) CH3OH (22,1–31,2), b) CH3OH (147,8–156,9), c) 13CH3OH (121,11–120,12), and d) CH3OH vt = 1 (71,6–61,5) towards cores A and B. For the different panels, the starting (and increasing contour levels in terms of σ, in Jy beam-1 km s-1) are 5σ (5σ, with σ = 0.18), 5σ (10σ, with σ = 0.11), 5σ (5σ, with σ = 0.25), and 5σ (10σ, with σ = 0.19) for core A; and 5σ (5σ, with σ = 0.11), 5σ (5σ, with σ = 0.06), 3σ (3σ, with σ = 0.14), and 5σ (5σ, with σ = 0.11) for core B. Offsets are measured with respect to the phase center (see Sect. 2). The red crosses indicate the position of the dust emission peaks of cores A and B (see Table 2). The synthesized beam and spatial scale are shown in the bottom-right corner (upper panels). The energy of the upper level of each transition is indicated in the bottom-right corner (lower panels). The dashed line is intended to stress that the velocity intervals used to compute the moment maps are different for cores A and B.

Open with DEXTER

thumbnail Fig. 5

Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CH3CN K = 2 at 349.426 GHz, b) CH3OH (22,1–31,2), c) H2CS (101,9–91,8), and d) HC3N (37–36) towards core C. The starting (and increasing) contour levels for the four panels are 1σ (1σ), with σ equal to 0.18, 0.22, 0.20 and 0.14 Jy beam-1 km s-1, for the four panels. Other symbols as in Fig. 4.

Open with DEXTER

3.2.3. High-velocity SiO outflow emission

We have analyzed the SiO (8–7) line emission to search for outflow signatures, since this molecule is known to be one of the best jet/outflow tracers in star forming regions (e.g., Gusdorf et al. 2008a,b; López-Sepulcre et al. 2011; Sánchez-Monge et al. 2013d). Channel maps of the SiO (8–7) emission are displayed in Fig. 7, with the 870 μm continuum emission and hot core (CH3CN) line emission overlaid. The SiO emission appears to span a wide range of velocities, however, line blending in hot cores A and B hinders the identification of high-velocity emission at the position of the two cores.

To overcome this problem, we identified three different regions to search for high-velocity SiO emission. These regions (named “out-S”, “out-SW” and “out-NE”) are located far from the cores identified in Fig. 1, thus avoiding contamination from the chemically-rich hot cores (see Fig. 7). The spectra of the emission integrated over the three regions are shown in Fig. 8, together with the spectra of the C17O (3–2) emission as a reference of the dense gas emission. The SiO emission at the “out-NE” position, spanning a range from 32 up to 60 km s-1, is clearly red-shifted with respect to the velocities of cores A and B (31.9 and 29.8 km s-1, respectively), and to the C17O line emission (with the peak at ~32.3 km s-1). In the “out-S” position, the SiO emission spans a range of velocities from 0 up to 32 km s-1, i.e., blue-shifted with respect to the C17O line emission. Finally, at the “out-SW” position, the SiO line has a narrow profile similar to the C17O line, both in velocity and linewidth. This suggests that the SiO emission at this position is associated with to the large-scale dense gas emission, rather than with high-velocity outflows.

thumbnail Fig. 6

Deconvolved size (in arcsec) for different molecular transitions (listed in Table 1) versus the upper level energy of the transition for cores A (top panel) and B (bottom panel). In both panels, black circles correspond to CH3CN, open circles to CH3CN vibrationally excited, red pentagons to CH3OH, and blue stars to transitions of other molecules. The typical error in size is shown in the top panel.

Open with DEXTER

In Fig. 9, we show the zeroth-order moment maps of the SiO emission integrated over two different velocity ranges: (1) blue-shifted velocities in the range 25.2–28.6 km s-1; and (2) red-shifted velocities in the range 40.8–57.7 km s-1. These velocity ranges have been chosen to avoid contamination of the hot core species associated with cores A and B (see middle panel in Fig. 8). The blue-shifted emission is located in between cores A and B, and thus it is difficult to determine the powering source. On the other hand, SiO red-shifted emission seems to be clearly emanating from core A. While the configuration of the SiO high-velocity emission does not resemble a bipolar outflow, it is worth comparing the SiO map with that of the free-free 3.6 cm continuum obtained by Gibb et al. (2003) with the Very Large Array (see black contours in Fig. 9). One sees two radio peaks that are lined with core A and the red-shifted SiO emission, suggesting that both the radio continuum and SiO emission might be tracing a jet from the core, oriented in the SW-NE direction.

In the following we test if the SiO and centimeter continuum emission could be arising from the same outflow. Anglada (1996) found an observational relation between the outflow momentum rate, , and the radio continuum luminosity, Sνd2, of the jet associated with the outflow. This observational relation is in agreement with the shock induced ionization mechanism proposed by Curiel et al. (1989) to explain the radio emission of the jets. The relation reads as (1)where Ω/4π is an efficiency factor that indicates the fraction of the solid angle that is shocked (~0.1 for low-mass radiojets/outflows, Anglada 1996), Sν is the flux density at 5 GHz, and d the heliocentric distance. The flux density of the 3.6 cm continuum emission of the radiojet is ~0.25–0.40 mJy (depending on the sources/knots considered), which corresponds to a flux density of ~0.18–0.30 mJy at 5 GHz, assuming a typical spectral index of 0.6. This results in an outflow momentum rate of 2.7–4.5 × 10-3  M yr-1 km s-1. We compared this number with the outflow momentum rate that can be computed from the SiO high-velocity emission. We followed the same procedure described in López-Sepulcre et al. (2009), and derived an outflow momentum rate of 1.5 × 10-3 and 0.4 × 10-3  M yr-1 km s-1, for the red and blue-shifted emission lobes shown in Fig. 9, assuming a SiO abundance of 3 × 10-9 (Sánchez-Monge et al. 2013d). The derived outflow momentum rate for the red-shifted lobe agrees, within a factor of 2–3, with that obtained from the radio continuum emission. The small discrepancy may be related to the uncertainties involved in these calculations (e.g., SiO abundance, Ω/ 4π efficiency factor). Observations, sensitive to scales ~2′′, of different outflow tracers should be performed to confirm the scenario in which core A is powering an outflow-jet directed approximately NE-SW.

3.3. Spectral energy distribution

Figure 10 shows the spectral energy distribution (SED) of G35.20N, obtained from the Herschel infrared Galactic Plane Survey (Hi-GAL; Molinari et al. 2010a,b) images complemented with literature and archival data at wavelengths ranging from mid-infrared to millimeter. In particular, we used the following data: Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE; Benjamin et al. 2003); Midcourse Space Experiment (MSX; Price et al. 1999); MIPSGAL/Spitzer (Carey et al. 2009); Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010); APEX Telescope Large Area Survey of the Galaxy (ATLASGAL; Schuller et al. 2009); Mooney et al. (1995); De Buizer (2006); Zhang et al. (2013). We stress that inclusion of the Hi-GAL data at 70 and 160 μm, represents an important improvement with respect to previous estimates, as it allows us to sample the SED around the peak, i.e., the spectral region responsible for most of the bolometric luminosity.

The bolometric luminosity has been derived by integrating the SED over the full observed spectral distribution. We note that different spatial apertures have been used at each wavelength due to the different angular resolution of the data. We obtain ~3 × 104 L, thus confirming the value quoted in the rms database4, but less than the luminosity (~0.7–2.2 × 105 L) derived by Zhang et al. (2013) from a model taking into account the photons which may escape through the outflow cavity (the so-called flashlight effect), but considering a single high-mass protostar forming within a massive core of radius 0.1 pc (which seems to be inconsistent with the fragmentation seen in our high-angular resolution observations).

4. Analysis and discussion

4.1. Temperature and mass estimates

The large number of methyl cyanide (CH3CN) transitions detected towards cores A, B and C, together with the methanol (CH3OH) transitions found in almost all the cores, can be used to obtain an estimate of the temperature of the cores in G35.20N. We fitted the observed methyl cyanide and methanol spectra for all the cores using the CASA interface of the XCLASS software called myXCLASS5 described in Schilke et al. (2001), Comito et al. (2005) and Zernickel et al. (2012). The myXCLASS software includes entries from the CDMS (Cologne Database for Molecular Spectroscopy; Müller et al. 2001) and JPL (Jet Propulsion Laboratory; Pickett et al. 1998) catalogs. With this information, it is possible to simulate a spectrum of a molecule, by providing parameters such as the excitation temperature, abundance, size of the source, linewidth and velocity. Once an initial guess is set up, the parameters can be fitted using the program MAGIX (Modeling and Analysis Generic Interface for eXternal numerical codes; Bernst et al. 2011; Möller et al. 2013). Figures B.1 to B.4 show the observed and myXCLASS fitted CH3CN and CH3OH spectra towards the cores in G35.20N. We fitted each molecule with a single component, i.e., the five free parameters are the same for all the lines of a given species. A more complex model, taking into account temperature and density gradients, might be necessary to explain the variations in size shown in Fig. 6. Finally, we note that the myXCLASS software takes the line opacity into account. The fits do not include only ground state transitions of the two species, but also vibrationally and torsionally excited as well as transitions from isotopologues such as CHCN and 13CH3OH. We assumed the relative abundance [12C/13C] to be 50 for a galactocentric distance of 6.5 kpc (Wilson & Rood 1994). The parameters of the best fits are listed in Table 3. The temperature ranges typically between 100–300 K for cores A–D, in both methanol and methyl cyanide species. For cores A and C, the estimated CH3CN temperature is large (close to 300 K) but with a large uncertainty as shown in the χ2 plots (see Figs. B.1B.4). For core D, the CH3CN emission is too weak and cannot be fitted, while from the CH3OH line fitting we measure a temperature of ~120 K (cooler than the three strongest cores). The column densities range between 1–4 × 1016 cm-2 for CH3CN, and 1–6 × 1018 cm-2 for CH3OH.

thumbnail Fig. 7

SiO (8–7) channel map of the G35.200.74 N region, averaged over 2 km s-1 wide velocity intervals. The central velocity of each channel is indicated in the upper left corner, and the velocity to take as reference is +32 km s-1. Gray contour: 5σ level of the 870 μm ALMA continuum emission shown in Fig. 1. Black contour: 50% intensity level of the CH3CN (19–18) K = 2 integrated emission as shown in Figs. A.1 and 5. The structures named “out-S”, “out-SW” and ‘out-NE’ are discussed in Sect. 3.2.3, while the spectra of these three regions are shown in Fig. 8.

Open with DEXTER

The mass of each core has been computed from (2)with S870 μm the dust continuum flux density at 870 μm (see Table 2), and D the distance to the source. We have adopted a dust opacity of 1.75 cm2 g-1 at 870 μm (Ossenkopf & Henning 1994), with a gas-to-dust ratio of 100. We assume T = Tex, where Tex is obtained from the myXCLASS analysis. Since CH3CN and CH3OH may be biased towards warm gas (>100 K), but the dust may also be tracing cooler material, we have set a lower limit for the temperature of 30 K. In Table 4, we give the range of temperatures used and the derived masses. For cores E and F, without detection of hot core molecules, we assume a temperature range of 30–50 K. The masses of the cores in G35.20N are in the range 1–10  M.

We have also calculated the source-averaged column density, Ncore, and volume density, ncore, under the assumption that the cores are spherically symmetric and homogeneous: (3)(4)where μ = 2.3 is the mean molecular weight, and θS is the deconvolved angular diameter (see Table 2). The results are given in Table 4. The mean column and volume densities are 1–10 × 1024 cm-2 and 2–9 × 108 cm-3, with core A being the densest. From these values and the CH3CN and CH3OH column densities (see Table 3), we derive a CH3CN fractional abundance of 0.10–1.1 × 10-8, 0.21–1.3 × 10-8 and 0.27–2.0 × 10-8, for cores A, B, and C, respectively. For CH3OH, we obtain 0.15–1.7 × 10-6, 0.45–2.6 × 10-6, 1.3–9.8 × 10-6 and 0.11–0.47 × 10-6, for cores A, B, C, and D. These abundances6 are similar to the values found toward other hot cores associated with high-mass and intermediate-mass YSOs (e.g., Charnley et al. 1992; Fuente et al. 2005, 2009; Sánchez-Monge et al. 2010), and only slightly lower than the abundances (10-8 for CH3CN, and 10-6 for CH3OH) measured in the strongest hot cores (e.g., Remijan et al. 2004; Bisschop et al. 2007; Mookerjea et al. 2007; Wang et al. 2010). We also computed the source-averaged surface density as , where the source diameter, DS, is taken from Table 2.

Table 3

Best myXCLASS fit parameters for the CH3CN and CH3OH species (see Figs. B.1B.4).

The G35.20N region is believed to be forming high-mass stars, based on the large luminosity (see Sect. 3.3), association with maser emission, and measurements of the mass of the embedded stars (Sánchez-Monge et al. 2013b). However, from Table 4, one sees that the values of the core masses are typically small, so one may wonder whether such a limited amount of gas is consistent with the formation of a massive star. We note that, as a matter of fact, there are quite a few B-type protostars associated with cores of a few  M, such as IRAS 20126+4104 (Cesaroni et al. 2005; Beltrán et al. 2011), W3 IRS5 (Wang et al. 2013) and G35.03+0.35 (Beltrán et al. 2014). It is possible that cores that “light” are only the remnants of the original envelopes enshrouding the protostars, while most of the material has already been accreted and/or dispersed. As suggested by Cesaroni (2005), the “light” cores might be forming only one or very few massive stars, unlike the “heavy” cores >100  M, which likely contain multiple high-mass stars. The surface densities measured toward the cores in G35.20N (see Table 4) are significantly higher than 0.15–1 g cm-2, which is the theoretical minimum value needed to form massive stars (Krumholz & McKee 2008; Butler & Tan 2012). We note, however, that the theoretical minimum surface density usually refers to clumps or clouds larger than the cores detected in G35.20N. We also compared the mass and size of the six cores identified in G35.20N with the empirical relation for a dense core to form a massive star as derived by Kauffmann & Pillai (2010, see also Elia et al. 2013. In Fig. 11, we see that the six cores as well as the whole elongated structure7 lay above the mass-size relation. However, we note that this relation is derived for clumps/clouds with sizes 0.1–10 pc, and it might be no longer a valid description for cores with sizes ~0.01 pc. Further observations are thus required to confirm if cores D to F may form high-mass stars, or only intermediate/low-mass counterparts.

Finally, we note that the mass of the whole elongated structure derived by us is several times less than that of ~130  M estimated by Gibb et al. (2003) and Birks et al. (2006) at lower angular resolution, assuming a dust temperature of 33 K and the opacities of Ossenkopf & Henning (1994). This indicates that our ALMA observations resolve out part of the large-scale, lower-density envelope, and are only sensitive to more compact and dense material, as suggested from the missing flux estimation done in Sect. 3.1. We also note that a fraction of the mass available in the elongated structure might be eventually accreted onto the different dense cores, and increase the final mass of the stars forming there.

thumbnail Fig. 8

Spectra of the SiO (8–7) emission (black line) toward different positions in the G35.20N star forming region: positions “out-NE” (top), “out-S” (middle) and “out-SW” (bottom) as indicated in Fig. 7. The emission has been obtained by integrating over areas of 1.34, 0.45 and 1.19 arcsec2, which yield to conversion factors of 8.5, 25.3 and 9.6 K Jy-1, respectively for the three panels. The red line corresponds to the C17O (3–2) line emission toward the three positions. In the middle panel, the blue dotted line corresponds to the SiO (8–7) line emission at the position of core A (over an area of 1.37 arcsec2, which yields to a conversion factor of 8.9 K Jy-1). The velocity component at ~16 km s-1 corresponds to two blended CH3CHO transitions. The two vertical dotted lines indicate the velocities of core A (31.8 km s-1) and core B (29.6 km s-1).

Open with DEXTER

thumbnail Fig. 9

Maps of the blue- (blue contours) and red-shifted (red contours) SiO (8–7) emission overlaid on the CH3CN (19–18) K = 2 velocity field (color scale) towards core A. The maps have been obtained by averaging the emission over the velocity intervals 25.2–28.6 km s-1 and 40.8–57.7 km s-1. Contour levels range from 0.05 to 0.27 in steps of 0.05 Jy beam-1 km s-1 and from 0.15 to 1.19 in steps of 0.15 Jy beam-1 km s-1, for blue and red-shifted emission, respectively. Black contours are a map of the 3.6 cm continuum emission (Gibb et al. 2003), with levels from 0.039 to 0.533 in steps of 0.026 mJy beam-1. The dashed line shows the direction of a possible jet associated with core A.

Open with DEXTER

thumbnail Fig. 10

Continuum spectrum of the G35.20N star forming region. Triangles denote lower limits. Different colors correspond to different surveys and instruments, as indicated in the figure.

Open with DEXTER

Table 4

Physical parameters of the cores detected at 870 μm (see Fig. 1).

thumbnail Fig. 11

Mass-radius relation for the cores (black crosses) and the elongated structure (blue cross) in G35.20N. The vertical length of the crosses denotes the uncertainty in the mass (see Table 4), while the horizontal length denotes a 10% error in the radius of the core (Table 2). The red solid line indicate the criteria for massive star formation by Kauffmann & Pillai (2010), with the gray area corresponding to the low-mass star formation region.

Open with DEXTER

4.2. Dynamical status of the cores

We estimated the virial masses of the cores identified in Fig. 1 as Mvir = 105 [ DS/ pc ] [ Δv/ km s-1 ] 2 M (e.g., MacLaren et al. 1988; Bertoldi & McKee 1992), where DS is the diameter of the core from Table 2, Δv is the line full width at half maximum (FWHM), and the numerical value accounts for a uniform density distribution across the core. Different numerical factors have to be considered for non-uniform density distributions, e.g., 95 for nr-1, and 63 for nr-2 (MacLaren et al. 1988). This results in a virial mass smaller by a factor 1.1 and 1.7, respectively. The value of Δv is obtained from an average of the CH3CN and CH3OH spectral line fits (see Table 3). The virial mass derived from these hot core tracers is listed in Col. 8 of Table 4. An estimate of the virial mass is also obtained by using the H13CO+ line, which is detected throughout the whole elongated structure (see Fig. 3) and hence towards all the cores. We note, however, that extended structures are heavily resolved out in our observations. This results in fake absorption features in the spectra of the cores, affecting our estimates of the line widths, which might be underestimated. In Cols. 9 and 10 of Table 4, we list the lower limits obtained for the H13CO+ line FWHM and the corresponding virial masses.

From Table 4 one sees that the virial masses are typically greater than (or similar to) the core masses. Since the former are lower limits (as explained above), we conclude that the virial parameter is Mvirial/Mcore1, with the cores possibly confined by the external pressure. It is worth noting, that other studies (e.g., Sánchez-Monge et al. 2013c; Kauffmann et al. 2013) have found similar results for the virial parameters when studying cores at scales ~0.5 pc. There is a global trend in which virial parameters close to unity are found at scales ~1–10 pc, whereas greater values of the virial parameter are obtained on smaller scales (~0.01–0.5 pc), suggesting that external pressure effects might be important mainly at small scales. These results provide some evidence in favor of the theoretical idea that on small scales, feedback from e.g., protostellar winds is sufficiently effective to keep regions out of virial equilibrium at scales below 0.5 pc.

4.3. Velocity field of cores A and B

As shown in Figs. 4, A.1 and A.2, the emission of almost all molecular species reveals a velocity gradient in cores A and B, oriented N–S in core A and SE–NW (i.e., along the major axis of the hot core) in core B. The velocity shift is approximately the same (~3–4 km s-1) for all species and transitions. In Figs. 12 and 13, we present position-velocity (PV) cuts along the direction where the velocity gradient is maximum: PA = 157° for core B and PA = 10° for core A. In the following we analyze these plots separately for the two cores.

thumbnail Fig. 12

Position-velocity plots along the direction with PA = 157° towards core B for different species and transitions. The offsets are measured from the dust peak position of core B (see Table 2), positive towards southeast. Contour levels start at 10%, increasing in steps of 10% of the peak: a) 0.53, b) 0.72, c) 0.08, d) 0.33, e) 0.10, f) 0.51, g) 0.56, h) 0.27, i) 0.09, j) 0.37, k) 0.02 and l) 0.52 Jy km s-1. The white solid line marks the border of the region where emission is expected for an edge-on Keplerian disk of radius (2600 AU) rotating about a 18  M central mass.

Open with DEXTER

thumbnail Fig. 13

Position-velocity plots along the direction with PA = 10° towards core A for different species and transitions. The offsets are measured from the dust peak position of core A (see Table 2), positive towards northeast. Contour levels start at 10%, increasing in steps of 10% of the peak: a) 0.87; b) 0.99; c) 0.13; d) 0.64; e) 0.23; f) 0.81; g) 0.76; h) 0.47; i) 0.27; j) 0.30; k) 0.06 and l) 0.34 Jy km s-1. The white solid line marks the border of the region where emission is expected for an edge-on Keplerian disk of radius (1500 AU) rotating about a 4  M central mass.

Open with DEXTER

4.3.1. Core B

Sánchez-Monge et al. (2013b) studied the velocity pattern of core B. The authors fitted the emission of several molecules in each velocity channel with a 2D Gaussian. The positions of the Gaussian peaks, which provide a picture of the velocity field (see Fig. 4 of Sánchez-Monge et al. 2013b), were fitted with a simple Keplerian disk model. From the best fit, the authors concluded that the velocity gradient in core B traces a Keplerian disk rotating about a central mass of 18  M, and inclined 19° with respect to the plane of the sky. The position angle of the disk is 157°. Comparing the molecular disk detected with ALMA with the centimeter continuum emission (Gibb et al. 2003) and the 4.5 μm Spitzer/IRAC emission (see Fig. 1a), the authors concluded that the disk in core B rotates around a central object which is the powering source of the NE–SW outflow/jet seen in the centimeter (Gibb et al. 2003), in the infrared (Dent et al. 1985b; Fuller et al. 2001; De Buizer 2006; Zhang et al. 2013) and in large-scale molecular outflow tracers (Dent et al. 1985a; Gibb et al. 2003; Birks et al. 2006).

The PV-plots towards core B (see Fig. 12) are consistent with the interpretation of Sánchez-Monge et al. (2013b), as they present a “butterfly” pattern typical of Keplerian rotation (see e.g., Cesaroni et al. 2005), with the emission extending to high velocities (±5 km s-1 with respect to the systemic velocity of ~30 km s-1) at small offsets from the dust continuum peak (expected to trace the star position), and small velocities (±1 km s-1) at large offsets. The white pattern shown in the figure, encompasses the region where emission is expected for an edge-on Keplerian disk of radius 1.̋2 (or 2600 AU), rotating about an 18  M central mass. This pattern is roughly consistent with the data, taking into account that it is obtained for zero line width and infinite angular resolution.

The PV-plots appear to change depending on the molecule. Species such as CH3CN, HC3N or SO2 have the emission peak to the southeast at blue-shifted velocities (~28 km s-1), whereas species such as CH3OH, CH3OCHO, and H2CS peak to the northwest at red-shifted velocities (~31 km s-1). We find that N-bearing molecules are typically stronger to the southeast of the continuum peak, while O-bearing molecules are stronger to the northwest. These differences could be a consequence of chemical variations across the disk, perhaps due to the interaction with nearby sources. An asymmetry was also observable in the Keplerian disk model fit presented in Sánchez-Monge et al. (2013b).

The Keplerian pattern is less evident in transitions like CH3CN K = 8, CH313CN K = 2 or CH3OCHO (277,21–267,20) (see Figs. 12c, e, and k), because these are detected only on one side of the disk. This suggests that not only the molecular abundance, but also the density distribution is not axially symmetric, similar to the case of the disk around the B-type protostar IRAS 20126+4104 (e.g., Cesaroni et al. 2014).

4.3.2. Core A

Core A is less resolved than core B, which hinders an analysis like that carried out by Sánchez-Monge et al. (2013b). The velocity gradient in this case is almost N–S. If this is tracing a rotating disk, one should see a jet/outflow perpendicular to it, i.e., in the E–W direction. The only clearly confirmed outflow with a similar orientation, is the one in the NE–SW direction observed with an angular resolution 5′′ (e.g., Gibb et al. 2003; Birks et al. 2006). As discussed later in Sect. 4.4, this outflow is likely powered by a precessing jet arising from core B (see also Sánchez-Monge et al. 2013b). In Sect. 3.2.3, we analyzed the SiO (8–7) transition covered in our ALMA observations. The combination of the SiO emission together with the radio continuum emission (see Fig. 9) suggests the presence of a jet/outflow centered in core A and oriented approximately ENE-WSW, i.e., roughly perpendicular to the velocity gradient in the core. This lends support to the hypothesis that core A might also be a rotating disk.

The PV-plots of core A shown in Fig. 13, reveal a compact source, with a small velocity gradient. Although the evidence of Keplerian rotation for core A is less clear than that in core B, in a few species the pattern looks slightly elongated (see panels a, b, f, and g), with the emission close to the center extending up to velocities as high as ±5 km s-1 (with respect to the systemic velocity of ~32 km s-1) and the low-velocity emission (±1.5 km s-1) detected at relatively large offsets (± with respect to the continuum peak). Assuming that one is indeed observing a Keplerian disk with a rotation velocity of 1.5 km s-1 at a radius of 0.̋7 (or 1500 AU), we derive a mass of ~4  M for the putative central star. Higher angular resolution data are needed to confirm whether core A is undergoing Keplerian rotation and derive a reliable estimate of the central mass.

4.3.3. Comparison with other disk candidates

thumbnail Fig. 14

Free-fall timescale to rotational period ratio (tff/trot) versus gas mass (Mgas) of known rotating disks or toroids. Open circles correspond to the objects studied by Beltrán et al. (2011), blue dot corresponds to source G35.03+0.35 (Beltrán et al. 2014), and the red dots correspond to cores A and B of G35.20N (this work). The connecting line corresponds to the different gas mass estimates for different temperatures (see Table 4). The gray areas show the regions where circumstellar disks and toroids are found. The dashed lines correspond to spheres of mass Mgas, containing a star of mass M (e.g., 0  M and 10  M; see Sect. 4.3.3).

Open with DEXTER

In Fig. 14, we follow the approach of Beltrán et al. (2011) to compare the rotating structures found in G35.20N with those found around B-type and O-type stars. These authors plot the ratio of the free-fall timescale to the rotational period (tff/trot) versus the gas mass (Mgas), for a sample of 20 high-mass YSOs. In this expression, tff is the free-fall timescale of the gas mass calculated from (sub)mm continuum emission, given by the expression (=1.83 × 107 [ Mgas/M ] -0.5 [ RS/ pc ] 1.5 yr). This quantity is an estimate of the timescale of collapse of the gas under its own gravity. The rotational period, trot = 2πRS/Vrot (=6.15 × 106 [ RS/ pc ] [ Vrot/ km s-1 ] -1 yr), is instead obtained from the observed velocity gradient and is basically the timescale needed for the disk to re-adjust its internal structure to the newly accreted material. Beltrán et al. (2011) find that rotating structures around B-type stars likely tracing circumstellar disks have higher tff/trot values than the massive toroids. In Fig. 14, we have added the results of cores A and B in G35.20N to the plot of Beltrán et al. (2011). In our calculations we consider a rotation velocity Vrot of 2 km s-1 and 1.2 km s-1, and a radius of 0.01 pc and 0.005 pc for cores B and A, respectively. The two cores of G35.20N are located in the region of Keplerian disks, which lends additional support to our hypotheses. In Fig. 14, we also show the theoretical curves corresponding to the value of tff/trot for a disk of mass Mgas rotating about a star of mass M. This can be obtained from the previous expressions of tff and trot, by substituting the rotation velocity with : (5)The dashed line shown in the figure correspond to M = 0 M and M = 10 M. From these it is clear that the objects lying in the “circumstellar disks” region are characterized by dynamically important levels of rotation, and have rotation timescales of the order (or shorter) than the dynamical timescale of collapse of the gas under its own self-gravity. Conversely, objects in the region labeled “massive toroids” are gas masses possibly supported by a combination of thermal pressure, turbulence, and magnetic fields with a little amount of rotation that however plays no role in their support. Thus, for these objects, the ratio tff/trott can take arbitrarily small values below the horizontal dashed line that corresponds to a homogeneous gas sphere in rotational equilibrium.

Finally, we estimated the velocity gradients of cores A and B, from the velocity shifts measured for the best resolved transitions (CH3CN K = 2, CH3OH (71,6–61,5) H2CS (101,9–91,8) and HC3N (37–36)). The shifts are ~4 km s-1 over an extent of ~0.02 pc for core B, and ~2.5 km s-1 over ~0.01 pc for core A (see Figs. 4A.2 and Figs. 13 and 12), which correspond to a velocity gradient at scales ~0.01 pc of about 250 km s-1 pc-1 and 200 km s-1 pc-1 for cores A and B, respectively. These values are consistent with the expected Keplerian velocity at a radius of 0.01 pc (0.005 pc) and around a central mass of 16  M (4  M), for core B (core A): ~2.5 km s-1 (~1.8 km s-1).

thumbnail Fig. 15

Image of the H2 2.12 μm line emission towards G35.20. The white and green contours are, respectively, maps of the 4.5 μm emission from the Spitzer/GLIMPSE survey (Benjamin et al. 2003) and 350 GHz continuum emission from our observations. The dotted curves outline the bipolar pattern.

Open with DEXTER

thumbnail Fig. 16

Same as Fig. 15, with overlaid also the 3.6 cm continuum map of Gibb et al. (2003) (cyan contours). The resolution of the IRAC image (white contours) has been enhanced by HiRes deconvolution (Velusamy et al. 2008). The dashed curve outlines the shape of the jet, which is bending by almost ~90° at an offset of ~11′′ to the north of core B.

Open with DEXTER

4.4. Precessing jet

As previously discussed, the velocity fields of cores A and B are suggestive of rotation about an axis oriented E–W and NE–SW, respectively. The latter direction is consistent with the axis of the bipolar nebula in Fig. 1a and corresponding outflow (see Sect. 1). This seems to confirm that the velocity gradient in core B is indeed due to rotation about a NE–SW axis. However, we know that several starforming cores are present in G35.20N, and consequently the butterfly-shaped nebula in Fig. 1a might be the result of the superposition of multiple outflows arising from these cores. Indeed, the morphology of both the mid-IR emission (see Fig. 2 of Sánchez-Monge et al. 2013b; Zhang et al. 2013) and 3.6 cm free-free continuum emission (Gibb et al. 2003) appears to indicate the existence of a jet from core B directed N–S, and hence not perpendicular to the velocity gradient in core B. Here, we argue that the non orthogonality between this jet and the core velocity gradient can be explained if the jet is precessing about an axis directed NE–SW. Little et al. (1998) were the first to suggest that such a bipolar structure could be created by a precessing jet, and Sánchez-Monge et al. (2013b) have proposed that a binary system at the center of the Keplerian disk in core B could be causing the precession.

To confirm this scenario, we have inspected the morphology of the H2 emission (a typical jet tracer) at different scales: from the scale of the core to that of the bipolar nebula, and beyond. For this purpose, we retrieved the UWISH2 image of the H2 2.12 μm line (Froebrich et al. 2011), shown in Fig. 15. One clearly sees that the line emission outlines a bipolar structure (dotted curves) that is the ideal extrapolation to a larger scale, of the butterfly-shaped nebula seen in the IRAC 4.5 μm image (white contours). Even more instructive is the small-scale structure of the H2 jet illustrated in Fig. 16. Putting together the information provided by the free-free 3.6 cm continuum emission (cyan contours), and the 2.12 μm and 4.5 μm images, one sees that the jet from core B is directed N–S, up to a distance of 8′′ (0.085 pc) from the origin, but at that point it sharply bends to the east (see dashed curve). This behavior is strongly suggestive of precession about an axis significantly inclined with respect to the N–S direction.

From the previous evidence we conclude that the morphology of the jet/outflow from core B is compatible with precession of the jet itself about an axis roughly perpendicular to the plane of the circumbinary disk in core B, as hypothesized by Sánchez-Monge et al. (2013b). We note that this interpretation does not rule out the possibility that other outflows are also present in this star forming region, as suggested by the complex structure of the SiO emission (see Sect. 3.2.3). Moreover, the morphology of the large-scale H2 emission seen in Fig. 15 is difficult to explain with a single precessing outflow. In particular, the branch of H2 emission extending to the east might be associated with the red lobe of the putative jet arising from core A (see Fig. 9).

thumbnail Fig. 17

Top: overlay of the 870 μm continuum emission (black contour) on the velocity field (colors) of the C17O (3–2) line, after tapering the line emission at 1′′. The contour level of the continuum emission is the 9.0 mJy beam-1 level (see Fig. 1). The velocity scale is shown in the color bar at the top of the panel, in km s-1. Other symbols as in Fig. 3. Bottom: position-velocity plot along the major axis of the elongated structure (PA = − 60°), centered at the position of core B, for C17O (contours) and CH3CN (gray scale). Vertical dotted lines mark the position of the cores identified in the continuum. The two horizontal dashed lines mark the velocities 30 km s-1 and 32 km s-1, associated with the hot core emission of cores B and A, respectively.

Open with DEXTER

thumbnail Fig. 18

As Fig. 17, but for H13CO+ (4–3) line.

Open with DEXTER

4.5. Nature of the elongated structure

The shape and orientation of the elongated structure in Fig. 1 suggest that this might be a 2D flattened structure seen edge-on, lying along the waist of the bipolar nebula observed at IR wavelengths. The existence of a large-scale, rotating disk-like structure has been indeed proposed by several authors. In particular, Little et al. (1985) and Brebner et al. (1987) revealed a large (~1 or 0.6 pc) elongated structure in NH3 and CS, oriented SE–NW, i.e., perpendicular to the large-scale outflow, and with a velocity gradient, which they interpret as a large rotating disk/toroid. More recently, Gibb et al. (2003) studied the G35.20N star forming complex with the BIMA interferometer in the H13CO+ (1–0) and H13CN (1–0) transitions, with an angular resolution of ~10′′. These authors find evidence of fragmentation across the same structure that they interpret as a fragmented, rotating envelope containing multiple young stellar objects, rather than a smooth rotating disk.

An alternative to the previous hypothesis, one should also consider the trivial possibility that what looks like a filament on the plane of the sky, is indeed a filament also in space. Recent continuum Herschel/Hi-GAL and spectral line observations of the interstellar medium have revealed that molecular clouds are highly filamentary (e.g., André et al. 2010; Molinari et al. 2010b; Busquet et al. 2013), which makes the filament hypothesis very plausible.

In the following we will discuss the nature of the elongated structure, confronting the two scenarios depicted above: 2D disk-like “pancake” seen edge-on, versus 1D filament.

First of all, it is important to determine the velocity field of the elongated structure, and verify whether this could be rotating about the outflow axis. For this purpose we need a suitable tracer, detected all over the structure. As already shown in Fig. 3 and discussed in Sect. 3.2.1, a few molecules show extended emission which is spatially coincident with the dust 870 μm continuum emission in Fig. 1. Despite the limited sensitivity of our observations to angular scales 2′′, the C17O (3–2) and H13CO+ (4–3) transitions can be imaged toward most of the elongated structure, and may hence be used to characterize its kinematics. In the top panels of Figs. 17 and 18, we present the 870 μm continuum emission map overlaid on the first order moment (velocity field) maps of the two lines. The southeastern part of the filament is associated with red-shifted emission (~36–37 km s-1), while the central and northwestern regions have emission at velocities ~31–32 km s-1. This velocity shift is more visible in the bottom panels of the two figures, which show the PV-plots of the two species along the major axis of the elongated structure (PA = −60°). In these plots, the velocity appears to shift gradually from ~37 km s-1 close to cores C and E (southeast), to ~31 km s-1 close to cores D and F (northwest).

The existence of a smooth velocity gradient seems to suggest that we are observing an edge-on, large-scale, rotating “pancake”. Under this hypothesis one can plot the specific angular momentum as a function of distance from the center (assumed to be the position of core B), complementing our data with data from the literature (Little et al. 1985; Gibb et al. 2003; López-Sepulcre et al. 2009) that allow us to sample the putative rotating “pancake” on scales as large as ~2 pc. In Fig. 19 we show this plot, which clearly indicates that angular momentum is dissipated from the large to the small scale. Magnetic braking might be a plausible explanation of this effect, but this requires a magnetic field (B) perpendicular to the disk, whereas recent polarimetric observations (Qiu et al. 2013) suggest that the B field is mainly aligned with the elongated structure. We will come back to the role of the B field later on.

Another problem with the “pancake” scenario is that core B would be counter-rotating with respect to the large-scale velocity field. Even more striking is that cores A and B appear to be kinematically inconsistent with respect to the putative “pancake”, as one can see by comparing the gray-scale map (hot-core tracers) to the contour map (large-scale molecular emission) in Figs. 17 and 18. A discrepancy between the velocities of the cores and that of the structure enshrouding them is difficult to justify in a rotating “pancake”, but can be explained if the two cores are created by the dynamical interaction of different filaments. Is there any evidence of such a dynamical environment around G35.20?

Figure 20 presents a large-scale (~5 or 3 pc) view of the G35.20N star forming region. The Spitzer/IRAC 8.0 μm image (Benjamin et al. 2003) reveals a series of elongated dark features converging towards the position of G35.20N, e.g., the darkest feature, likely associated with the densest filament, towards the northwest, a curved arc-like structure towards the northeast, and two elongated structures to the south perpendicular to the main dark feature. The dust emission at 250 μm (from Hi-GAL; Molinari et al. 2010a) matches the infrared dark features quite well, with the strongest emission associated with G35.20N, where the 8.0 μm emission also peaks. The overall structure (panel a) resembles that of the SDC335.5790.272 infrared dark cloud (Peretto et al. 2013), with multiple filaments converging and accreting matter towards the center, where a number of massive stars are forming. Our ALMA observations might be imaging the most prominent filament, where the accreting material is converging.

As previously noted, to discriminate an edge-on “pancake” from a filament, one could use the direction of the magnetic field, which is expected to be perpendicular to the disk plane and parallel to the filament. Recent SMA polarization measurements of the submillimeter continuum emission from G35.20N (Qiu et al. 2013) have found that the magnetic field on scales of ~15′′ (0.16 pc) is almost parallel to the elongated structure, a result that supports the filament hypothesis. However, in this case the rotation axes of cores A and B should also be parallel to the filament, whereas they are significantly inclined with respect to it. The direction of the rotation axes is instead consistent with the small-scale direction of the magnetic field, obtained from observations of OH and CH3OH masers (Hutawarakorn & Cohen 1999; Surcis et al. 2012), which indicates that the B field is oriented in the (N)E-(S)W direction (see Fig. 5 of Surcis et al. 2012) over ~0.̋3 (660 AU). The change of the direction of B from the large to the small scale suggests that turbulence could play an important role in regulating the fragmentation of G35.20N at scales ~0.1 pc. Turbulence might also explain the drastic change in the magnetic field orientation (by approximately 90°) toward the southeast of the elongated structure (see Fig. 2 by Qiu et al. 2013).

One should also establish whether a big “pancake” could be unstable against perturbations and hence produce the observed fragments. This can be investigated by estimating the Toomre stability parameter for disks, Q. Assuming solid-body rotation (consistent with the PV-plots in Figs. 17 and 18), this can be expressed as , with Ω angular velocity, Σ surface density, and ΔV line FWHM. Ω ≃ 40 km s-1 pc-1 is estimated from the slope of the PV-plots in Figs. 17 and 18, while Σ is given by the ratio between the mass of the elongated structure (30  M) and the disk surface (πR2, where 2 R = 0.15 pc is the length of the elongated structure). In order to obtain Q ≤ 1, the condition for instability, the line width must be ΔV ≤ 0.7 km s-1, a reasonable value for the initial conditions of the unperturbed gas, prior to star formation. We hence conclude that the “pancake” should break up into fragments.

thumbnail Fig. 19

Specific angular momentum (top panel) and velocity gradient (bottom panel) at different radii traced by the available molecular line data. Small red and blue stars correspond to the H13CO+ and C17O ALMA line data presented in this work, respectively. Black five-point stars: NH3 (1,1) and CS (2–1) line data from Little et al. (1985). Black filled circles: H13CN (1–0) and H13CO+ (1–0) line data from Gibb et al. (2003). Black eight-point star: C18O (2–1) line data from López-Sepulcre et al. (2009).

Open with DEXTER

thumbnail Fig. 20

a) Overlay of the Hi-GAL 250 μm continuum emission (white contours) on the Spitzer/IRAC 8.0 μm image (colors) of the G35.200.74 region. Contour levels start at 65 mJy beam-1, increasing in steps of 0.39 Jy beam-1 until 1.885 Jy beam-1, then levels increase in steps of 1.98 Jy beam-1 until the peak: 12 Jy beam-1. b) Same as Fig. 1b.

Open with DEXTER

The filament hypothesis may also explain the observed distribution of cores A–F. A filament (or isothermal gas cylinder) supported by “turbulence pressure” against its self-gravity has a maximum mass per unit length equal to (6)(Nagasawa 1987; Wang et al. 2014), over which it starts to fragment. The theoretical separation of the fragments is given by the expression (7)where ΔV is the line FWHM due to both turbulent at the thermal broadening, and relates to the velocity dispersion by . The line width to be used in these equations should be the initial value, before fragmentation, and not the current value, measured by us, which is likely to be broadened by additional turbulence injected by the newly formed stars. In practice, one sees that a value of ΔV = 0.7 km s-1 is required for the theoretical expression to match the observed spacing between the cores (~0.023 pc), after considering a filament with a mass ~30  M, a size 0.15 pc×0.013 pc, and a gas density of ρ ≃ 10-16 g cm-3. In addition, this line width of 0.7 km s-1 results in a critical mass per unit length of ~40  M pc-1, while the measured mass per unit length in our filament is ~200  M pc-1. Such an initial value of ΔV is plausible and we conclude that the filament should undergo fragmentation. Similar structures undergoing fragmentation have been found in other star forming regions such as G11.110.12 and G28.34+0.06 (Wang et al. 2011, 2014).

Based on these results, both a rotating “pancake’ and a filament appear to be unstable and hence consistent with the observed clumpiness. However, one should take into account also an intriguing feature of the elongated structure, namely the regular spacing between the cores. The distance between two adjacent cores lies between 0.020 pc and 0.036 pc, quite a narrow range of values. It seems unlikely that a random 2D distribution of fragments over a “pancake” can produce such a regular distribution, when seen in an edge-on projection on the plane of the sky. In conclusion, for this reason and the orientation of the magnetic field, we tend to favor the filament interpretation of the elongated structure.

5. Summary

We have observed the G35.20N star forming region with ALMA at 870 μm (350 GHz), achieving an angular resolution ~ (~900 AU at 2.19 kpc), and covering a broad frequency range including dense gas (e.g., H13CO+, C17O), outflow (e.g., SiO) and hot core (e.g., CH3CN, CH3OH) tracers. Our conclusions can be summarized as follows:

  • The continuum emission shows an elongated dust structure (length ~0.15 pc and width ~0.013 pc), likely tracing the densest part of the elongated structure observed at larger scales, and perpendicular to the infrared and molecular large-scale outflow (e.g., Gibb et al. 2003). This structure is fragmented in six identified dense cores, although additional substructure is seen. The masses of fragments are between 110  M, and have sizes of around 1600 AU. The cores appear regularly spaced with a mean separation of ~0.023 pc.

  • Three out of the six dense cores show strong emission in complex organic molecules typical of hot cores. We fitted a large number of methyl cyanide (CH3CN) and methanol (CH3OH) spectral lines and derived temperatures of around 150–250 K for the three cores and relative abundances of 0.2–2 × 10-8 for CH3CN and 0.6–5 × 10-6 for CH3OH.

  • The two densest and most chemically-rich cores (cores A and B) show a coherent velocity field with a velocity gradient almost aligned with the dust elongated structure. Interestingly, the orientation of both velocity gradients is opposite, with red-shifted emission toward the north for core B and toward the south for core A. These velocity gradients are consistent with Keplerian disks rotating about central masses of 4  M for core A and 18  M for core B, i.e., intermediate/high-mass stars.

  • Core B seems to power a precessing jet/outflow oriented in the NE-SW direction. The combination of the IRAC 4.5 μm, the UWISH2 H2 2.12 μm and the VLA 3.6 cm images together with our ALMA continuum and line maps suggests that the outflow/jet associated with core B is undergoing precession, in agreement with the binary system hypothesized by Sánchez-Monge et al. (2013b). Our ALMA high-velocity SiO (7–6) map together with the VLA 3.6 cm continuum map suggest that core A might be powering an outflow in the east-west direction. Further observations of outflow/shock tracers are required to confirm this scenario.

  • The emission of dense gas tracers such as H13CO+ and C17O is also extended and coincident with the dust emission. This large molecular structure is consistent, in morphology and velocity, with the SMA maps reported by Qiu et al. (2013). Based on the velocity field, the orientation of the magnetic field (as reported by Qiu et al. 2013), and the regularly spaced fragmentation (~0.023 pc), we interpret the elongated dust structure as the densest part of a 1D filament fragmenting and forming high-mass stars.


1

The Common Astronomy Software Applications (CASA) software can be downloaded at http://casa.nrao.edu

3

The 6.7 GHz CH3OH maser observations carried out with the Japanese VLBI Network (JVN; Sugiyama et al. 2008) and the European VLBI Neatwork (EVN; Surcis et al. 2012) were performed without phase-referencing and it is thus impossible to establish the absolute position of the masers. However, the CH3OH maser distribution seems to match the OH maser distribution (Brebner et al. 1987; Hutawarakorn & Cohen 1999); since the latter are associated with cores A and B, it seems likely that the methanol masers are also associated with these two cores.

4

The rms database is available at http://www.ast.leeds.ac.uk/cgi-bin/RMS/RMS_DATABASE.cgi

6

We note that missing flux may affect differently the spatial distributions of gas and dust, introducing some uncertainties in the reported abundances.

7

The mass of the whole elongated structure imaged with ALMA has been computed by adding the masses of the cores to the mass of the extended emission. The flux density of the whole structure is ~3.3 Jy, with ~2 Jy coming from the cores (see Table 2). We assume that the remaining flux (~1.3 Jy) comes from gas with a temperature of 30–50 K to derive a total mass of the structure (extended emission plus dense cores) of 21–37  M.

Acknowledgments

We are grateful to the Italian ARC node for the usage of their computer facilities during the cleaning and imaging process. We are also grateful to Ray Furuya for providing us the ASTE spectrum, and to Keping Qiu for providing the SMA spectra. Á.S.-M. is grateful to Peter Schilke, Thomas Möller and Alexander Zernickel for helping with the analysis of the myXCLASS software. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00275.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.

References

Online material

Appendix A: Hot core intensity and velocity maps

In Figs. A.1 and A.2, we present the zeroth (integrated intensity) and first (velocity field) moment maps of various hot-core tracers detected towards cores A and B. In Sect. 3.2.2, we describe how the moments have been obtained, while the symbols shown in these two figures are described in the caption of Fig. 4.

thumbnail Fig. A.1

Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CH3CN K = 2 at 349.426 GHz; b) CH3CN K = 8 at 349.024 GHz; c) CHCN K = 2 at 349.254 GHz; and d) CH3CN v8 = 1K,l = 3,1 at 350.552 GHz towards cores A and B. For the different panels, the starting (and increasing contour levels in terms of σ, in Jy beam-1 km s-1) are 5σ (10σ, with σ = 0.15), 5σ (10σ, with σ = 0.15), 3σ (3σ, with σ = 0.10), and 5σ (10σ, with σ = 0.10) for core A; and 5σ (5σ, with σ = 0.09), 4σ (4σ, with σ = 0.09), 2σ (2σ, with σ = 0.06), and 3σ (3σ, with σ = 0.06) for core B. See details and description of symbols in Fig. 4.

Open with DEXTER

thumbnail Fig. A.2

Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CH3OCHO (266,20–256,19); b) CH3OCHO (277,21–267,20); c) H2CS (101,9–91,8), d) HNCO (161,16–151,15); e) SO2 (184,14–183,15); and f) HC3N (37–36) towards cores A and B. For the different panels, the starting (and increasing contour levels in terms of σ, in Jy beam-1 km s-1) are 5σ (10σ, with σ = 0.06), 3σ (3σ, with σ = 0.04), 5σ (10σ, with σ = 0.16), 5σ (10σ, with σ = 0.10), 5σ (5σ, with σ = 0.20), and 5σ (10σ, with σ = 0.11) for core A; and 3σ (3σ, with σ = 0.03), 1σ (1σ, with σ = 0.02), 5σ (5σ, with σ = 0.09), 5σ (5σ, with σ = 0.06), 5σ (5σ, with σ = 0.10), and 5σ (10σ, with σ = 0.07) for core B. See details and description of symbols in Fig. 4.

Open with DEXTER

Appendix B: CH3OH and CH3CN myXCLASS fits

In Figs. B.1B.4 we show the spectra of the bulge of methanol (CH3OH, black line in the top panels) and methyl cyanide (CH3CN, black line in the bottom panels) lines observed toward cores A, B, C and D. For this last core, the methyl cyanide lines are too weak and cannot be fitted. The observed spectra have been fitted (red lines) using myXCLASS (see Sect. 4.1). The procedure searches, by minimizing the χ2, for the best fit of five parameters: size (in arcsec), temperature (in K), column density (in cm-2), linewidth (in km s-1) and LSR velocity (in km s-1). The blue lines in the small panels show the χ2 of the fits for different values of the parameters. In Table 3, we list the values of the parameters obtained in the best fits. We note that we have simultaneously fitted optically thick and thin transitions (e.g., ground state, vibrationally excited and isotopologues of CH3CN and CH3OH), thus avoiding possible degeneracies (see Zernickel et al. 2012) between the parameters of the fit (e.g., column density and size).

thumbnail Fig. B.1

Large panels: observed CH3OH (top panel) and CH3CN (bottom panel) lines toward core A in G35.200.74 N. The red line shows the best fitted obtained by myXCLASS (see Sect. 4.1). The fit in the top panel includes ground state CH3OH and 13CH3OH and torsionally excited CH3OH lines, while the fit in the bottom panel includes ground state CH3CN and CHCN and vibrationally excited CH3CN lines in the frequency ranges shown in the panels. Small panels: χ2 values obtained from the fits minimization. From left to right the panels show the size, the temperature, the column density, the linewidth and the LSR velocity.

Open with DEXTER

thumbnail Fig. B.2

Same as Fig. B.1 for core B.

Open with DEXTER

thumbnail Fig. B.3

Same as Fig. B.1 for core C.

Open with DEXTER

thumbnail Fig. B.4

Same as Fig. B.1 for core D, and only for CH3OH lines.

Open with DEXTER

All Tables

Table 1

Molecular transitions studied in this paper.

Table 2

Parameters of the cores detected in the 870 μm ALMA continuum map (see Fig. 1).

Table 3

Best myXCLASS fit parameters for the CH3CN and CH3OH species (see Figs. B.1B.4).

Table 4

Physical parameters of the cores detected at 870 μm (see Fig. 1).

All Figures

thumbnail Fig. 1

a)Spitzer 4.5 μm image (in logarithmic scale) of the star forming region G35.200.74 N, overlaid with a contour map of the 870 μm (350 GHz) continuum emission obtained with ALMA. The IR image has been obtained by applying HiRes deconvolution (Velusamy et al. 2008) to the Spitzer/IRAC data. The submillimeter map has been corrected for primary beam attenuation. Contours start at 5σ, increasing in steps of 6σ (with σ = 1.8 mJy beam-1), up to a maximum level of 201.7 mJy beam-1. The dashed white circle denotes the primary beam of the ALMA 12 m antennas (~9′′), while the ALMA synthesized beam (, PA = 46°) is shown in the bottom-right corner of the panel. b) Close-up of the central region that shows the 870 μm continuum emission map obtained with ALMA. Six cores (from A to F) have been identified and marked in the image (see Table 2). Yellow crosses mark the position of OH masers (Hutawarakorn & Cohen 1999). The ALMA synthesized beam is shown in the bottom-right corner of the panel. The spatial scale is also indicated in both panels.

Open with DEXTER
In the text
thumbnail Fig. 2

Spectra of the full frequency range surveyed by the ALMA observations, toward the six cores identified in the continuum image (see Fig. 1). The spectra are extracted by averaging the intensity over an area of 1.37, 2.35, 1.04, 0.55, 0.75, 0.44 arcsec2, yielding conversion factors of 8.9, 5.0, 11.3, 21.7, 15.7, 27.1 K Jy-1, for cores A to F respectively. The brightness temperature scale of cores A and B is different from that of the rest of the cores. The absorption seen in some spectra is likely produced by the filtering of extended emission.

Open with DEXTER
In the text
thumbnail Fig. 3

Overlay of the 870 μm continuum emission (red contour) on the integrated intensity maps (colors) of C17O, C34S, H13CO+, and SiO. The contour level of the continuum emission corresponds to the 9.0 mJy beam-1 intensity (see Fig. 1). The intensity scale for the integrated emission for each molecule is shown in the color bar at the right of each panel, in K km s-1. The synthesized beam of the continuum emission (, PA = 46°) is shown in the bottom-left corner of each panel. For each molecule, we show the synthesized beam of the tapered maps in the bottom-right corner: a) , PA = 67°; b) , PA = 67°; c) , PA = 69°; and d) , PA = 69°.

Open with DEXTER
In the text
thumbnail Fig. 4

Maps of the integrated intensity (contours) and velocity field (color maps in the bottom panels) of a) CH3OH (22,1–31,2), b) CH3OH (147,8–156,9), c) 13CH3OH (121,11–120,12), and d) CH3OH vt = 1 (71,6–61,5) towards cores A and B. For the different panels, the starting (and increasing contour levels in terms of σ, in Jy beam-1 km s-1) are 5σ (5σ, with σ = 0.18), 5σ (10σ, with σ = 0.11), 5σ (5σ, with σ = 0.25), and 5σ (10σ, with σ = 0.19) for core A; and 5σ (5σ, with σ = 0.11), 5σ (5σ, with σ = 0.06), 3σ (3σ, with σ = 0.14), and 5σ (5σ, with σ = 0.11) for core B. Offsets are measured with respect to the phase center (see Sect. 2). The red crosses indicate the position of the dust emission peaks of cores A and B (see Table 2). The synthesized beam and spatial scale are shown in the bottom-right corner (upper panels). The energy of the upper level of each transition is indicated in the bottom-right corner (lower panels). The dashed line is intended to stress that the velocity intervals used to compute the moment maps are different for cores A and B.

Open with DEXTER
In the text
thumbnail Fig. 5

Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CH3CN K = 2 at 349.426 GHz, b) CH3OH (22,1–31,2), c) H2CS (101,9–91,8), and d) HC3N (37–36) towards core C. The starting (and increasing) contour levels for the four panels are 1σ (1σ), with σ equal to 0.18, 0.22, 0.20 and 0.14 Jy beam-1 km s-1, for the four panels. Other symbols as in Fig. 4.

Open with DEXTER
In the text
thumbnail Fig. 6

Deconvolved size (in arcsec) for different molecular transitions (listed in Table 1) versus the upper level energy of the transition for cores A (top panel) and B (bottom panel). In both panels, black circles correspond to CH3CN, open circles to CH3CN vibrationally excited, red pentagons to CH3OH, and blue stars to transitions of other molecules. The typical error in size is shown in the top panel.

Open with DEXTER
In the text
thumbnail Fig. 7

SiO (8–7) channel map of the G35.200.74 N region, averaged over 2 km s-1 wide velocity intervals. The central velocity of each channel is indicated in the upper left corner, and the velocity to take as reference is +32 km s-1. Gray contour: 5σ level of the 870 μm ALMA continuum emission shown in Fig. 1. Black contour: 50% intensity level of the CH3CN (19–18) K = 2 integrated emission as shown in Figs. A.1 and 5. The structures named “out-S”, “out-SW” and ‘out-NE’ are discussed in Sect. 3.2.3, while the spectra of these three regions are shown in Fig. 8.

Open with DEXTER
In the text
thumbnail Fig. 8

Spectra of the SiO (8–7) emission (black line) toward different positions in the G35.20N star forming region: positions “out-NE” (top), “out-S” (middle) and “out-SW” (bottom) as indicated in Fig. 7. The emission has been obtained by integrating over areas of 1.34, 0.45 and 1.19 arcsec2, which yield to conversion factors of 8.5, 25.3 and 9.6 K Jy-1, respectively for the three panels. The red line corresponds to the C17O (3–2) line emission toward the three positions. In the middle panel, the blue dotted line corresponds to the SiO (8–7) line emission at the position of core A (over an area of 1.37 arcsec2, which yields to a conversion factor of 8.9 K Jy-1). The velocity component at ~16 km s-1 corresponds to two blended CH3CHO transitions. The two vertical dotted lines indicate the velocities of core A (31.8 km s-1) and core B (29.6 km s-1).

Open with DEXTER
In the text
thumbnail Fig. 9

Maps of the blue- (blue contours) and red-shifted (red contours) SiO (8–7) emission overlaid on the CH3CN (19–18) K = 2 velocity field (color scale) towards core A. The maps have been obtained by averaging the emission over the velocity intervals 25.2–28.6 km s-1 and 40.8–57.7 km s-1. Contour levels range from 0.05 to 0.27 in steps of 0.05 Jy beam-1 km s-1 and from 0.15 to 1.19 in steps of 0.15 Jy beam-1 km s-1, for blue and red-shifted emission, respectively. Black contours are a map of the 3.6 cm continuum emission (Gibb et al. 2003), with levels from 0.039 to 0.533 in steps of 0.026 mJy beam-1. The dashed line shows the direction of a possible jet associated with core A.

Open with DEXTER
In the text
thumbnail Fig. 10

Continuum spectrum of the G35.20N star forming region. Triangles denote lower limits. Different colors correspond to different surveys and instruments, as indicated in the figure.

Open with DEXTER
In the text
thumbnail Fig. 11

Mass-radius relation for the cores (black crosses) and the elongated structure (blue cross) in G35.20N. The vertical length of the crosses denotes the uncertainty in the mass (see Table 4), while the horizontal length denotes a 10% error in the radius of the core (Table 2). The red solid line indicate the criteria for massive star formation by Kauffmann & Pillai (2010), with the gray area corresponding to the low-mass star formation region.

Open with DEXTER
In the text
thumbnail Fig. 12

Position-velocity plots along the direction with PA = 157° towards core B for different species and transitions. The offsets are measured from the dust peak position of core B (see Table 2), positive towards southeast. Contour levels start at 10%, increasing in steps of 10% of the peak: a) 0.53, b) 0.72, c) 0.08, d) 0.33, e) 0.10, f) 0.51, g) 0.56, h) 0.27, i) 0.09, j) 0.37, k) 0.02 and l) 0.52 Jy km s-1. The white solid line marks the border of the region where emission is expected for an edge-on Keplerian disk of radius (2600 AU) rotating about a 18  M central mass.

Open with DEXTER
In the text
thumbnail Fig. 13

Position-velocity plots along the direction with PA = 10° towards core A for different species and transitions. The offsets are measured from the dust peak position of core A (see Table 2), positive towards northeast. Contour levels start at 10%, increasing in steps of 10% of the peak: a) 0.87; b) 0.99; c) 0.13; d) 0.64; e) 0.23; f) 0.81; g) 0.76; h) 0.47; i) 0.27; j) 0.30; k) 0.06 and l) 0.34 Jy km s-1. The white solid line marks the border of the region where emission is expected for an edge-on Keplerian disk of radius (1500 AU) rotating about a 4  M central mass.

Open with DEXTER
In the text
thumbnail Fig. 14

Free-fall timescale to rotational period ratio (tff/trot) versus gas mass (Mgas) of known rotating disks or toroids. Open circles correspond to the objects studied by Beltrán et al. (2011), blue dot corresponds to source G35.03+0.35 (Beltrán et al. 2014), and the red dots correspond to cores A and B of G35.20N (this work). The connecting line corresponds to the different gas mass estimates for different temperatures (see Table 4). The gray areas show the regions where circumstellar disks and toroids are found. The dashed lines correspond to spheres of mass Mgas, containing a star of mass M (e.g., 0  M and 10  M; see Sect. 4.3.3).

Open with DEXTER
In the text
thumbnail Fig. 15

Image of the H2 2.12 μm line emission towards G35.20. The white and green contours are, respectively, maps of the 4.5 μm emission from the Spitzer/GLIMPSE survey (Benjamin et al. 2003) and 350 GHz continuum emission from our observations. The dotted curves outline the bipolar pattern.

Open with DEXTER
In the text
thumbnail Fig. 16

Same as Fig. 15, with overlaid also the 3.6 cm continuum map of Gibb et al. (2003) (cyan contours). The resolution of the IRAC image (white contours) has been enhanced by HiRes deconvolution (Velusamy et al. 2008). The dashed curve outlines the shape of the jet, which is bending by almost ~90° at an offset of ~11′′ to the north of core B.

Open with DEXTER
In the text
thumbnail Fig. 17

Top: overlay of the 870 μm continuum emission (black contour) on the velocity field (colors) of the C17O (3–2) line, after tapering the line emission at 1′′. The contour level of the continuum emission is the 9.0 mJy beam-1 level (see Fig. 1). The velocity scale is shown in the color bar at the top of the panel, in km s-1. Other symbols as in Fig. 3. Bottom: position-velocity plot along the major axis of the elongated structure (PA = − 60°), centered at the position of core B, for C17O (contours) and CH3CN (gray scale). Vertical dotted lines mark the position of the cores identified in the continuum. The two horizontal dashed lines mark the velocities 30 km s-1 and 32 km s-1, associated with the hot core emission of cores B and A, respectively.

Open with DEXTER
In the text
thumbnail Fig. 18

As Fig. 17, but for H13CO+ (4–3) line.

Open with DEXTER
In the text
thumbnail Fig. 19

Specific angular momentum (top panel) and velocity gradient (bottom panel) at different radii traced by the available molecular line data. Small red and blue stars correspond to the H13CO+ and C17O ALMA line data presented in this work, respectively. Black five-point stars: NH3 (1,1) and CS (2–1) line data from Little et al. (1985). Black filled circles: H13CN (1–0) and H13CO+ (1–0) line data from Gibb et al. (2003). Black eight-point star: C18O (2–1) line data from López-Sepulcre et al. (2009).

Open with DEXTER
In the text
thumbnail Fig. 20

a) Overlay of the Hi-GAL 250 μm continuum emission (white contours) on the Spitzer/IRAC 8.0 μm image (colors) of the G35.200.74 region. Contour levels start at 65 mJy beam-1, increasing in steps of 0.39 Jy beam-1 until 1.885 Jy beam-1, then levels increase in steps of 1.98 Jy beam-1 until the peak: 12 Jy beam-1. b) Same as Fig. 1b.

Open with DEXTER
In the text
thumbnail Fig. A.1

Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CH3CN K = 2 at 349.426 GHz; b) CH3CN K = 8 at 349.024 GHz; c) CHCN K = 2 at 349.254 GHz; and d) CH3CN v8 = 1K,l = 3,1 at 350.552 GHz towards cores A and B. For the different panels, the starting (and increasing contour levels in terms of σ, in Jy beam-1 km s-1) are 5σ (10σ, with σ = 0.15), 5σ (10σ, with σ = 0.15), 3σ (3σ, with σ = 0.10), and 5σ (10σ, with σ = 0.10) for core A; and 5σ (5σ, with σ = 0.09), 4σ (4σ, with σ = 0.09), 2σ (2σ, with σ = 0.06), and 3σ (3σ, with σ = 0.06) for core B. See details and description of symbols in Fig. 4.

Open with DEXTER
In the text
thumbnail Fig. A.2

Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CH3OCHO (266,20–256,19); b) CH3OCHO (277,21–267,20); c) H2CS (101,9–91,8), d) HNCO (161,16–151,15); e) SO2 (184,14–183,15); and f) HC3N (37–36) towards cores A and B. For the different panels, the starting (and increasing contour levels in terms of σ, in Jy beam-1 km s-1) are 5σ (10σ, with σ = 0.06), 3σ (3σ, with σ = 0.04), 5σ (10σ, with σ = 0.16), 5σ (10σ, with σ = 0.10), 5σ (5σ, with σ = 0.20), and 5σ (10σ, with σ = 0.11) for core A; and 3σ (3σ, with σ = 0.03), 1σ (1σ, with σ = 0.02), 5σ (5σ, with σ = 0.09), 5σ (5σ, with σ = 0.06), 5σ (5σ, with σ = 0.10), and 5σ (10σ, with σ = 0.07) for core B. See details and description of symbols in Fig. 4.

Open with DEXTER
In the text
thumbnail Fig. B.1

Large panels: observed CH3OH (top panel) and CH3CN (bottom panel) lines toward core A in G35.200.74 N. The red line shows the best fitted obtained by myXCLASS (see Sect. 4.1). The fit in the top panel includes ground state CH3OH and 13CH3OH and torsionally excited CH3OH lines, while the fit in the bottom panel includes ground state CH3CN and CHCN and vibrationally excited CH3CN lines in the frequency ranges shown in the panels. Small panels: χ2 values obtained from the fits minimization. From left to right the panels show the size, the temperature, the column density, the linewidth and the LSR velocity.

Open with DEXTER
In the text
thumbnail Fig. B.2

Same as Fig. B.1 for core B.

Open with DEXTER
In the text
thumbnail Fig. B.3

Same as Fig. B.1 for core C.

Open with DEXTER
In the text
thumbnail Fig. B.4

Same as Fig. B.1 for core D, and only for CH3OH lines.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.