A necklace of dense cores in the high-mass star forming region G35.20-0.74 N: ALMA observations

A necklace of dense cores in the high-mass star forming region G35.200.74 N: ALMA observations

Á. Sánchez-Monge I. Physikalisches Institut, Universität zu Köln, Zülpicher Str. 77, D-50937, Köln, Germany Osservatorio Astrofisico di Arcetri, INAF, Largo Enrico Fermi 5, I-50125, Firenze, Italy E-mail: sanchez@ph1.uni-koeln.de (Á. S.-M.)    M. T. Beltrán Osservatorio Astrofisico di Arcetri, INAF, Largo Enrico Fermi 5, I-50125, Firenze, Italy    R. Cesaroni Osservatorio Astrofisico di Arcetri, INAF, Largo Enrico Fermi 5, I-50125, Firenze, Italy    S. Etoka Hamburger Sternwarte, Gojenbergsweg 112, D-21029, Hamburg, Germany Jodrell Bank Centre for Astrophysics, School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, UK    D. Galli Osservatorio Astrofisico di Arcetri, INAF, Largo Enrico Fermi 5, I-50125, Firenze, Italy    M. S. N. Kumar Centro de Astrofísica da Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal    L. Moscadelli Osservatorio Astrofisico di Arcetri, INAF, Largo Enrico Fermi 5, I-50125, Firenze, Italy    T. Stanke ESO, Karl-Schwarzschilk-Strasse 2, D-85748 Garching bei München, Germany    F. F. S. van der Tak SRON Netherlands Institute for Space Research, PO Box 800, 9700, AV Groningen, The Netherlands Kapteyn Astronomical Institute, University of Groningen, 9700 AV Groningen, The Netherlands    S. Vig Dpt. of Earth and Space Science, Indian Institute of Space Science and Technology, Thiruvananthapuram, 695 547 Kerala, India    C. M. Walmsley Osservatorio Astrofisico di Arcetri, INAF, Largo Enrico Fermi 5, I-50125, Firenze, Italy Dublin Institute for Advanced Studies (DIAS), 31 Fitzwilliam Place, Dublin 2, Ireland    K.-S. Wang Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands    H. Zinnecker SOFIA Science Center, NASA Ames Research Center, Mailstop 232-12, Moffett Field, CA 94035, USA    D. Elia Istituto di Astrofisica e Planetologia Spaziali (IAPS-INAF), via Fosso del Cavaliere 100, I-00133 Roma, Italy    S. Molinari Istituto di Astrofisica e Planetologia Spaziali (IAPS-INAF), via Fosso del Cavaliere 100, I-00133 Roma, Italy    E. Schisano Infrared Processing and Analysis Center, Institute of Technology, Pasadena, CA 91125, USA
Received; accepted
Key Words.:
stars: formation – stars: massive – ISM: molecules – ISM: jets and outflows – ISM: individual objects: G35.200.74 N
Abstract

Context:The formation process of high-mass stars (with masses ) is still poorly understood, and represents a challenge from both the theoretical and observational points of view. The advent of the Atacama Large Millimeter Array (ALMA) is expected to provide observational evidence to better constrain the theoretical scenarios.

Aims:The present study aims at characterizing the high-mass star forming region G35.200.74 N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk.

Methods:We used the radio-interferometer ALMA to observe the G35.200.74 N region in the submillimeter continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g., HCO, CO), molecular outflows (e.g., SiO), and hot cores (e.g., CHCN, CHOH). These observations were complemented with infrared and centimeter data.

Results:The ALMA 870 m continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide; full width at half maximum) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1–10  and sizes 1600 AU (spatial resolution 960 AU). The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as HCO or CO is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2–2 for CHCN and 0.6–5 for CHOH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4–18 . Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A.

Conclusions:The elongated dust structure in G35.200.74 N is fragmented into a number of dense cores that may form high-mass stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming high-mass stars.

1 Introduction

High-mass stars (O- and B-type stars with masses ) 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 K.-S. 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 ) 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 (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.

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.030.35. Results on G35.030.35 are presented in a parallel paper by Beltrán et al. (A&A, submitted), 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 B. 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 185560136 and has a bolometric luminosity 1–  (e.g., Gibb et al., 2003; Sánchez-Monge et al., 2013b; Zhang Y. 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 (P.A.58°) of a bipolar outflow observed in CO (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 CO (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 Y. 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., NH, CO, CS, HCN, HCO), 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) to SE (36 km s). 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 HCO and HCN 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  (Sánchez-Monge et al., 2013b).

Frequency
Observation (MHz) (K)
CHOCHO (27–26) A 335 015.90 443.5
CHOH (2–3) 335 133.57 044.7
CHOH (12–12) 335 560.21 192.7
CHOH (7–6) 335 582.02 079.0
CHOCHO (27–26) A 336 111.32 277.9
CHOH (14–15) 336 438.21 488.2
HCN (37–36) 336 520.08 306.9
CHOH (12–12) 336 865.15 197.1
CHOCHO (26–25) E 336 889.20 235.5
CHOCHO (26–25) A 336 918.18 235.5
CO (3–2) 337 061.05 032.4
CS (7–6) 337 396.69 050.2
CHOH (7–6 337 748.83 488.5
CHOH (7–6 337 969.44 390.1
HCS (10–9) 338 083.20 102.4
SO (18–18) 338 305.99 196.8
CHOH (7–6) 338 583.22 112.7
CHOH (7–6) 338 639.80 102.7
HCO (4–3) 346 998.34 041.6
CHOH (14–14) 347 188.28 254.3
SiO (8–7) 347 330.63 075.0
CHOCHO (28–27) E 347 628.34 306.8
CHOCHO (28–27) E 348 049.89 266.1
CHOCHO (28–27) A 348 065.97 266.1
HCS (10–9) 348 534.36 105.2
CHOH (14–14) 349 107.00 260.2
CHCN (19–18) 349 280.77
CHCN (19–18) 349 453.70
HNCO (16–15) 350 333.35 186.2
CHCN ( 350 552.44

The frequencies are obtained from the CDMS (Cologne Database for Molecular Spectroscopy; Müller et al. 2001) and JPL (Jet Propulsion Laboratory; Pickett et al. 1998) catalogs. is the upper level energy, and the Boltzmann constant.
The frequency is that of the =0 component.
The frequency is that of the =2, component.

Table 1: Molecular transitions studied in this paper
Figure 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), up to a maximum level of 201.7 mJy beam. 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.
Peak position Source Diameter
Core h m s (Jy beam) (Jy) (arcsec) (AU)
A 18 58 12.945 01 40 37.36 0.208 0.414 0900
B 18 58 13.025 01 40 35.92 0.076 0.838 1800
C 18 58 13.126 01 40 33.20 0.056 0.717 1600
D 18 58 12.801 01 40 40.00 0.083 0.571 1300
E 18 58 13.195 01 40 30.96 0.048 0.842 1900
F 18 58 12.699 01 40 41.12 0.049 1.181 2600

 Primary beam corrected peak intensity () and flux density () integrated within the 5 contour level.
 Deconvolved average diameter of the 50% contour level (following Sánchez-Monge et al., 2013a).

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

Finally, emission from different maser species (HO, OH, CHOH) 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 CHOH masers (Surcis et al., 2012), indicate strong (5–12%) linear polarization, with the magnetic field oriented almost E–W, P.A.80–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.

Figure 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 arcsec, yielding conversion factors of 8.9, 5.0, 11.3, 21.7, 15.7, 27.1 K Jy, 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.
Figure 3: Overlay of the 870 m continuum emission (red contour) on the integrated intensity maps (colors) of CO, CS, HCO, and SiO. The contour level of the continuum emission corresponds to the 9.0 mJy beam 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. 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°.

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. 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) = 18 58 13, (J2000) = 01 40 360. Flux, gain, and bandpass calibrations were obtained through observations of Neptune and J1751096. The data were calibrated and imaged using CASA111The Common Astronomy Software Applications (CASA) software can be downloaded at http://casa.nrao.edu (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 Section 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. 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) varies between 5 and 20 mJy beam, 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 GILDAS222http://www.iram.fr/IRAMFR/GILDAS 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.

Figure 4: Maps of the integrated intensity (contours) and velocity field (color maps in the bottom panels) of a) CHOH (2–3), b) CHOH (14–15), c) CHOH (12–12), and d) CHOH =1 (7–6) towards cores A and B. For the different panels, the starting (and increasing contour levels in terms of , in Jy beam km s) 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.
Figure 5: Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CHCN =2 at 349.426 GHz, b) CHOH (2–3), c) HCS (10–9), and d) HCN (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 km s, for the four panels. Other symbols as in Fig. 4.

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 probably333The 6.7 GHz CHOH 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 CHOH 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. 6.7-GHz CHOH 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 to 200 mJy beam, 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 versus 8.5 Jy beam) 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., CHCN, CHOH, CHOCHO, NHCHO, HCN, CHOH, CHOCHO). Out of the rest of the cores, only core C is likely associated with emission in complex molecules such as CHCN and CHOH. 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 CO, CS, HCO (see Fig. 3) and weak CHOH, 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 (CO, CS, HCO 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, and compare them with the dust continuum emission. For these molecules, the line emission coincides with the dust continuum emission, with CO and HCO being detected toward all the cores and tracing the elongated structure, and SiO and CS 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 HCO 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 HCO (1–0) and HCN (1–0) lines (e.g., Gibb et al., 2003). In Sect. 4.5, we will analyze the kinematics of this extended molecular structure.

Figure 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 CHCN, open circles to CHCN vibrationally excited, red pentagons to CHOH, and blue stars to transitions of other molecules. The typical error in size is shown in the top panel.
Figure 7: SiO (8–7) channel map of the G35.200.74 N region, averaged over 2 km s 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. Gray contour: 5 level of the 870 m ALMA continuum emission shown in Fig. 1. Black contour: 50% intensity level of the CHCN (19–18) integrated emission as shown in Figs. 21 and 5. The structures named ‘out-S’, ‘out-SW’ and ‘out-NE’ are discussed in Section 3.2.3, while the spectra of these three regions are shown in Fig. 8.

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: CHCN, CHCN, and vibrationally excited CHCN, as well as CHOH, CHOCHO, HCN, HNCO, HCS, and SO. 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 CHCN =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; R. Furuya, private communication). 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 23, 24 and 25 show a close-up view of the CHCN and CHOH transitions detected towards cores A, B and C, respectively. The observational setup (see Sect. 2) permits to cover up to the =9 transition of CHCN (19–18), and up to the =5 line of CHCN (19–18), with upper level energies of 745 K and 517 K, respectively. All of these CHCN lines have been detected in the three cores. Regarding the CHCN, only a few lines are not blended (e.g., =2, =4), while others clearly overlap with methyl cyanide or methanol lines (e.g., =3, =5). Several vibrationally excited CHCN 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 (CHOH) and methyl formate (CHOCHO) 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 HCS, HCN, SO or HNCO.

In Figs. 4, 21 and 22, we present the zeroth (integrated intensity) and first (velocity field) moment maps of different hot core species detected towards cores A and B: CHCN =2 and 8, CHCN =2, and CHCN =1 =3,1, as well as different transitions of CHOH, CHOCHO, HCS, HNCO, SO, and HCN. 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 for core A, and 25.9–33.0 km s 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., CHCN =1, or CHOCHO at 335.015 GHz) or less abundant species (e.g., CHCN =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: CHCN (19–18), CHOH (2–3), HCS (10–9), and HCN (37–36). The moments were computed using the velocity interval 26.1–41.0 km s. 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.

Figure 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 arcsec, which yield to conversion factors of 8.5, 25.3 and 9.6 K Jy, respectively for the three panels. The red line corresponds to the CO (3–2) line emission toward the three positions. In 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 arcsec, which yields to a conversion factor of 8.9 K Jy). The velocity component at 16 km s corresponds to two blended CHCHO transitions. The two vertical dotted lines indicate the velocities of core A (31.8 km s) and core B (29.6 km s).

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 (CHCN) 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 CO (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, is clearly red-shifted with respect to the velocities of cores A and B (31.9 and 29.8 km s, respectively), and to the CO line emission (with the peak at 32.3 km s). In the ‘out-S’ position, the SiO emission spans a range of velocities from 0 up to 32 km s, i.e., blue-shifted with respect to the CO line emission. Finally, at the ‘out-SW’ position, the SiO line has a narrow profile similar to the CO 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.

Figure 9: Maps of the blue- (blue contours) and red-shifted (red contours) SiO (8–7) emission overlaid on the CHCN (19–18) =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 and 40.8–57.7 km s. Contour levels range from 0.05 to 0.27 in steps of 0.05 Jy beam km s and from 0.15 to 1.19 in steps of 0.15 Jy beam km s, 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. The dashed line shows the direction of a possible jet associated with core A.

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, and (2) red-shifted velocities in the range 40.8–57.7 km s. 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, , 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 is an efficiency factor that indicates the fraction of the solid angle that is shocked (0.1 for low-mass radiojets/outflows, Anglada 1996), is the flux density at 5 GHz, and 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  yr km s. 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 and 0.4  yr km s, for the red and blue-shifted emission lobes shown in Fig. 9, assuming a SiO abundance of (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, 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.

Figure 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.
CHCN CHOH
v v size v v size
Core (km s) (km s) (arcsec) (K) (cm) (km s) (km s) (arcsec) (K) (cm)
A 31.8 5.2 0.24 309 31.8 4.0 0.34 192
B 29.6 2.8 0.31 148 29.9 2.9 0.36 177
C 32.4 8.1 0.12 280 32.0 5.4 0.09 176
D 31.2 1.5 0.07 120
Table 3: Best myXCLASS fit parameters for the CHCN and CHOH species (see Figs. 2326).
Core (K) () ( cm) ( cm) (g cm) (km s) () (km s) ()
A 030–250 11.2–1.0 40.2–3.6 44.3–4.0 156.0–14.0 4.6 09.7 3.2 4.7
B 030–150 05.4–0.9 04.7–0.8 02.6–0.4 18.9–3.1 2.9 07.7 2.9 7.7
C 030–200 03.1–0.4 03.7–0.5 02.4–0.3 13.7–1.8 6.8 37.7 2.6 5.5
D 030–100 03.7–0.9 07.0–1.7 05.6–1.4 24.8–6.0 1.5 01.5 2.9 5.6
E 30–50 02.6–1.4 02.3–1.2 01.2–0.7 8.2–4.4 1.8 3.1
F 30–50 04.5–2.4 02.0–1.1 00.8–0.4 7.5–4.0 1.0 1.3

 Temperature derived from the rotational diagram method when possible (see Sect. 4.1).
 Mass, , column density, , volume density, , and surface density, , estimated using a dust temperature equal to , and a dust opacity of 1.75 cm g at 870 m (see Sect. 4.1).

Table 4: Physical parameters of the cores detected at 870 m (see Fig. 1).

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 Y. 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  , thus confirming the value quoted in the RMS database444The RMS database is available at http://www.ast.leeds.ac.uk/cgi-bin/RMS/RMS_DATABASE.cgi, but less than the luminosity (0.7– ) derived by Zhang Y. 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 (CHCN) transitions detected towards cores A, B and C, together with the methanol (CHOH) 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 myXCLASS555myXCLASS can be downloaded from https://www.astro.uni-koeln.de/projects/schilke/myXCLASSInterface 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 23 to 26 show the observed and myXCLASS fitted CHCN and CHOH 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 CHOH. We assumed the relative abundance [C/C] 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 CHCN temperature is large (close to 300 K) but with a large uncertainty as shown in the plots (see Figs. 2326). For core D, the CHCN emission is too weak and cannot be fitted, while from the CHOH line fitting we measure a temperature of 120 K (cooler than the three strongest cores). The column densities range between 1– cm for CHCN, and 1– cm for CHOH.

The mass of each core has been computed from

(2)

with the dust continuum flux density at 870 m (see Table 2), and the distance to the source. We have adopted a dust opacity of 1.75 cm g at 870 m (Ossenkopf & Henning, 1994), with a gas-to-dust ratio of 100. We assume , where is obtained from the myXCLASS analysis. Since CHCN and CHOH 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 .

We have also calculated the source-averaged column density, , and volume density, , under the assumption that the cores are spherically symmetric and homogeneous:

(3)
(4)

where is the mean molecular weight, and is the deconvolved angular diameter (see Table 2). The results are given in Table 4. The mean column and volume densities are 1– cm and 2–9 cm, with core A being the densest. From these values and the CHCN and CHOH column densities (see Table 3), we derive a CHCN fractional abundance of 0.10–1.1, 0.21–1.3 and 0.27–2.0, for cores A, B, and C, respectively. For CHOH, we obtain 0.15–1.7, 0.45–2.6, 1.3–9.8 and 0.11–0.47, for cores A, B, C, and D. These abundances666We note that missing flux may affect differently the spatial distributions of gas and dust, introducing some uncertainties in the reported abundances. 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 ( for CHCN, and for CHOH) measured in the strongest hot cores (e.g., Remijan et al., 2004; Bisschop et al., 2007; Mookerjea et al., 2007; Wang K.-S. et al., 2010). We also computed the source-averaged surface density as , where the source diameter, , is taken from Table 2.

The G35.20N region is believed to be forming high-mass stars, based on the large luminosity (see Section 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 , such as IRAS 201264104 (Cesaroni et al., 2005; Beltrán et al., 2011), W3 IRS5 (Wang K.-S. et al., 2013) and G35.030.35 (Beltrán et al. A&A, submitted). 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 , 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, 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 ). In Fig. 11, we see that the six cores as well as the whole elongated structure777The 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 . 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  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 Section 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.

Figure 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.

4.2 Dynamical status of the cores

We estimated the virial masses of the cores identified in Fig. 1 as   (e.g., MacLaren et al., 1988; Bertoldi & McKee, 1992), where is the diameter of the core from Table 2, 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 , and 63 for (MacLaren et al., 1988). This results in a virial mass smaller by a factor 1.1 and 1.7, respectively. The value of is obtained from an average of the CHCN and CHOH spectral line fits (see Table 3). The virial mass derived from these hot core tracers is listed in column 8 of Table 4. An estimate of the virial mass is also obtained by using the HCO 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 columns 9 and 10 of Table 4, we list the lower limits obtained for the HCO 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 1, 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, 21 and 22, 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) 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: P.A.=157° for core B and P.A.=10° for core A. In the following we analyze these plots separately for the two cores.

Figure 12: Position-velocity plots along the direction with P.A.=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. 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  central mass.
Figure 13: Position-velocity plots along the direction with P.A.=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. 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  central mass.

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 , 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 Y. 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 with respect to the systemic velocity of 30 km s) at small offsets from the dust continuum peak (expected to trace the star position), and small velocities (1 km s) 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 12 (or 2600 AU), rotating about an 18  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 CHCN, HCN or SO have the emission peak to the southeast at blue-shifted velocities (28 km s), whereas species such as CHOH, CHOCHO, and HCS peak to the northwest at red-shifted velocities (31 km s). 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 CHCN =8, CHCN =2 or CHOCHO (27–26) (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 201264104 (e.g., Cesaroni et al. A&A, submitted).

Figure 14: Free-fall timescale to rotational period ratio () versus gas mass () 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.030.35 (Beltrán et al. A&A, submitted), 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 , containing a star of mass (e.g., 0  and 10 ; see Section 4.3.3).

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 ″ (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 Section 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 (with respect to the systemic velocity of 32 km s) and the low-velocity emission (1.5 km s) 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 at a radius of 07 (or 1500 AU), we derive a mass of  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

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 () versus the gas mass (), for a sample of 20 high-mass YSOs. In this expression, is the free-fall timescale of the gas mass calculated from (sub)mm continuum emission, given by the expression ( yr). This quantity is an estimate of the timescale of collapse of the gas under its own gravity. The rotational period, ( 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 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 of 2 km s and 1.2 km s, 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 for a disk of mass rotating about a star of mass . This can be obtained from the previous expressions of and , by substituting the rotation velocity with :

(5)

The dashed line shown in the figure correspond to =0  and =10 . 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 t 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 (CHCN =2, CHOH (7–6) HCS (10–9) and HCN (37–36)). The shifts are 4 km s over an extent of 0.02 pc for core B, and 2.5 km s over 0.01 pc for core A (see Figs. 422 and Figs. 13 and 12), which correspond to a velocity gradient at scales 0.01 pc of about 250 km s pc and 200 km s pc 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  (4 ), for core B (core A): 2.5 km s (1.8 km s).

Figure 15: Image of the H 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.
Figure 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.

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 Y. 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 H 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 UKIDSS image of the H 2.12 m line, 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 H 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 H emission seen in Fig. 15 is difficult to explain with a single precessing outflow. In particular, the branch of H emission extending to the east might be associated with the red lobe of the putative jet arising from core A (see Fig. 9).

Figure 17: Top: Overlay of the 870 m continuum emission (black contour) on the velocity field (colors) of the CO (3–2) line, after tapering the line emission at 1″. The contour level of the continuum emission is the 9.0 mJy beam level (see Fig. 1). The velocity scale is shown in the color bar at the top of the panel, in km s. Other symbols as in Fig. 3 Bottom: Position-velocity plot along the major axis of the elongated structure (P.A.=60°), centered at the position of core B, for CO (contours) and CHCN (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 and 32 km s, associated with the hot core emission of cores B and A, respectively.
Figure 18: As Figure 17, but for HCO (4–3) line.

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 NH 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 HCO (1–0) and HCN (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 CO (3–2) and HCO (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), while the central and northwestern regions have emission at velocities 31–32 km s. 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 (P.A.=60°). In these plots, the velocity appears to shift gradually from 37 km s close to cores C and E (southeast), to 31 km s 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 () perpendicular to the disk, whereas recent polarimetric observations (Qiu et al., 2013) suggest that the field is mainly aligned with the elongated structure. We will come back to the role of the 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 CHOH masers (Hutawarakorn & Cohen, 1999; Surcis et al., 2012), which indicates that the field is oriented in the (N)E-(S)W direction (see Fig. 5 of Surcis et al. 2012) over 03 (660 AU). The change of the direction of 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, . 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 line FWHM.  km s pc 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 ) and the disk surface (, where  pc is the length of the elongated structure). In order to obtain , the condition for instability, the line width must be  km s, 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.

Figure 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 HCO and CO ALMA line data presented in this work, respectively. Black five-point stars: NH (1,1) and CS (2–1) line data from Little et al. (1985). Black filled circles: HCN (1–0) and HCO (1–0) line data from Gibb et al. (2003). Black eight-point star: CO (2–1) line data from López-Sepulcre et al. (2009).
Figure 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, increasing in steps of 0.39 Jy beam until 1.885 Jy beam, then levels increase in steps of 1.98 Jy beam until the peak: 12 Jy beam. (b) Same as Fig. 1b.

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 K. et al., 2014), over which it starts to fragment. The theoretical separation of the fragments is given by the expression

(7)

where 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  km s 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 , a size 0.15 pc0.013 pc, and a gas density of  g cm. In addition, this line width of 0.7 km s results in a critical mass per unit length of 40  pc, while the measured mass per unit length in our filament is 200  pc. Such an initial value of 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.340.06 (Wang K. 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., HCO, CO), outflow (e.g., SiO) and hot core (e.g., CHCN, CHOH) 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 1-10 , 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 (CHCN) and methanol (CHOH) spectral lines and derived temperatures of around 150–250 K for the three cores and relative abundances of 0.2–2 for CHCN and 0.6–5 for CHOH.

  • The two densest and most chemically-rich cores (core 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  for core A and 18  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 UKIDSS H 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 HCO and CO 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.

Acknowledgements.
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 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

  • André et al. (2010) André, P., Men’shchikov, A., Bontemps, S., et al. 2010, A&A, 518, L102
  • Anglada (1996) Anglada, G. 1996, in Radio Emission from the Stars and the Sun, eds. A. R. Taylor, & J. M. Paredes, ASP Conf. Ser., 93, 3
  • Beltrán et al. (2011) Beltrán, M. T., Cesaroni, R., Neri, R., & Codella, C. 2011, A&A, 525, A151
  • Benjamin et al. (2003) Benjamin, R. A., Churchwell, E., Babler, B. L., et al. 2003, PASP, 115, 953
  • Bernst et al. (2011) Bernst, I., Schilke, P., Moeller, T., et al. 2011, in Astronomical Data Analysis Software and Systems XX, eds. I. N. Evans, A. Accomazzi, D. J. Mink, & A. H. Rots, ASP Conf. Ser., 442, 505
  • Bertoldi & McKee (1992) Bertoldi, F., & McKee, C. F. 1992, ApJ, 395, 140
  • Birks et al. (2006) Birks, J. R., Fuller, G. A., & Gibb, A. G. 2006, A&A, 458, 181
  • Bisschop et al. (2007) Bisschop, S. E., Jørgensen, J. K., van Dishoeck, E. F., & de Wachter, E. B. M. 2007, A&A, 465, 913
  • Boley et al. (2013) Boley, P. A., Linz, H., van Boekel, R., et al. 2013, A&A, 558, A24
  • Bonnell & Bate (2006) Bonnell, I. A., & Bate, M. R. 2006, MNRAS, 370, 488
  • Brebner et al. (1987) Brebner, G. C., Heaton, B., Cohen, R. J., & Davies, S. R. 1987, MNRAS, 229, 679
  • Brigg (1995) Briggs, D. 2005, Ph.D. Thesis, New Mexico Inst. of Mining and Technology
  • Busquet et al. (2013) Busquet, G., Zhang, Q., Palau, A., et al. 2013, ApJ, 764, L26
  • Butler & Tan (2012) Butler, M. J., & Tan, J. C. 2012, ApJ, 754, 5
  • Carey et al. (2009) Carey, S. J., Noriega-Crespo, A., Mizuno, D. R., et al. 2009, PASP, 121, 76
  • Carrasco-González et al. (2012) Carrasco-González, C., Galván-Madrid, R., Anglada, G., et al. 2012, ApJ, 752, L29
  • Cesaroni (2005) Cesaroni, R. 2005, in Massive Star Birth: a Crossroads of Astrophysics, IAU Symposium 227, eds. R. Cesaroni, M. Felli, E. Churchwell, & M. Walmsley (Cambridge University Press), 59
  • Cesaroni et al. (2005) Cesaroni, R., Neri, R., Olmi, L., et al. 2005, A&A, 434, 1039
  • Cesaroni et al. (2007) Cesaroni, R., Galli, D., Lodato, G., Walmsley, C. M., & Zhang, Q. 2007, in Protostars and Planets V, eds. B. Reipurth, D. Jewitt, & K. Keil (University of Arizona Press), 197
  • Charnley et al. (1992) Charnley, S. B., Tielens, A. G. G. M., & Millar, T. J. 1992, ApJ, 399, L71
  • Comito et al. (2005) Comito, C., Schilke, P., Phillips, T. G., et al. 2005, ApJS, 156, 127
  • Curiel et al. (1989) Curiel, S., Rodriguez, L. F., Canto, J., & Torrelles, J. M. 1989, Rev. Mexicana Astron. Astrofis., 17, 137
  • De Buizer (2006) De Buizer, J. M. 2006, ApJ, 642, L57
  • de Wit et al. (2011) de Wit, W. J., Hoare, M. G., Oudmaijer, R. D., et al. 2011, A&A, 526, L5
  • Dent et al. (1985a) Dent, W. R. F., Little, L. T., Kaifu, N., Ohishi, M., & Suzuki, S. 1985a, A&A, 146, 375
  • Dent et al. (1985b) Dent, W. R. F., Little, L. T., Sato, S., Ohishi, M., & Yamashita, T. 1985b, MNRAS, 217, 217
  • Elia et al. (2013) Elia, D., Molinari, S., Fukui, Y., et al. 2013, ApJ, 772, 45
  • Forster & Caswell (1989) Forster, J. R., & Caswell, J. L. 1989, A&A, 213, 339
  • Forster & Caswell (1999) Forster, J. R., & Caswell, J. L. 1999, A&AS, 137, 43
  • Fuente et al. (2005) Fuente, A., Neri, R., & Caselli, P. 2005, A&A, 444, 481
  • Fuente et al. (2009) Fuente, A., Castro-Carrizo, A., Alonso-Albi, T., et al. 2009, A&A, 507, 1475
  • Fuller et al. (2001) Fuller, G. A., Zijlstra, A. A., & Williams, S. J. 2001, ApJ, 555, L125
  • Gibb et al. (2003) Gibb, A. G., Hoare, M. G., Little, L. T., & Wright, M. C. H. 2003, MNRAS, 339, 1011
  • Gusdorf et al. (2008a) Gusdorf, A., Cabrit, S., Flower, D. R., & Pineau Des Forêts, G. 2008a, A&A, 482, 809
  • Gusdorf et al. (2008b) Gusdorf, A., Pineau Des Forêts, G., Cabrit, S., & Flower, D. R. 2008b, A&A, 490, 695
  • Heaton & Little (1988) Heaton, B. D., & Little, L. T. 1988, A&A, 195, 193
  • Hutawarakorn & Cohen (1999) Hutawarakorn, B., & Cohen, R. J. 1999, MNRAS, 303, 845
  • Kauffmann et al. (2013) Kauffmann, J., Pillai, T., & Goldsmith, P. F. 2013, ApJ, 779, 185
  • Kauffmann & Pillai (2010) Kauffmann, J., & Pillai, T. 2010, ApJ, 723, L7
  • Keto (2007) Keto, E. 2007, ApJ, 666, 976
  • Kraus et al. (2010) Kraus, S., Hofmann, K.-H., Menten, K. M., et al. 2010, Nature, 466, 339
  • Krumholz & McKee (2008) Krumholz, M. R., & McKee, C. F. 2008, Nature, 451, 1082
  • Krumholz et al. (2009) Krumholz, M. R., Klein, R. I., McKee, C. F., Offner, S. S. R., & Cunningham, A. J. 2009, Science, 323, 754
  • Little et al. (1985) Little, L. T., Dent, W. R. F., Heaton, B., Davies, S. R., & White, G. J. 1985, MNRAS, 217, 227
  • Little et al. (1998) Little, L. T., Kelly, M. L., & Murphy, B. T. 1998, MNRAS, 294, 105
  • López-Sepulcre et al. (2009) López-Sepulcre, A., Codella, C., Cesaroni, R., Marcelino, N., & Walmsley, C. M. 2009, A&A, 499, 811
  • López-Sepulcre et al. (2011) López-Sepulcre, A., Walmsley, C. M., Cesaroni, R., et al. 2011, A&A, 526, L2
  • MacLaren et al. (1988) MacLaren, I., Richardson, K. M., & Wolfendale, A. W. 1988, ApJ, 333, 821
  • McMullin et al. (2007) McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in Astronomical Data Analysis Software and Systems XVI, eds. R. A. Shaw, F. Hill, & D. J. Bell, ASP Conf. Ser., 376, 127
  • Molinari et al. (2010a) Molinari, S., Swinyard, B., Bally, J., et al. 2010a, PASP, 122, 314
  • Molinari et al. (2010b) Molinari, S., Swinyard, B., Bally, J., et al. 2010b, A&A, 518, L100
  • Möller et al. (2013) Möller, T., Bernst, I., Panoglou, D., et al. 2013, A&A, 549, A21
  • Mookerjea et al. (2007) Mookerjea, B., Casper, E., Mundy, L. G., & Looney, L. W. 2007, ApJ, 659, 447
  • Mooney et al. (1995) Mooney, T., Sievers, A., Mezger, P. G., et al. 1995, A&A, 299, 869
  • Müller et al. (2001) Müller, H. S. P., Thorwirth, S., Roth, D. A., & Winnewisser, G. 2001, A&A, 370, L49
  • Nagasawa (1987) Nagasawa, M. 1987, Progress of Theoretical Physics, 77, 635
  • Ossenkopf & Henning (1994) Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943
  • Peretto et al. (2013) Peretto, N., Fuller, G. A., Duarte-Cabral, A., et al. 2013, A&A, 555, A112
  • Pickett et al. (1998) Pickett, H. M., Poynter, R. L., Cohen, E. A., et al. 1998, J. Quant. Spec. Radiat. Transf., 60, 883
  • Price et al. (1999) Price, S. D., Egan, M. P., & Shipman, R. F. 1999, Astrophysics with Infrared Surveys: A Prelude to SIRTF, 177, 394
  • Qiu et al. (2013) Qiu, K., Zhang, Q., Menten, K. M., Liu, H. B., & Tang, Y.-W. 2013, ApJ, 779, 182
  • Remijan et al. (2004) Remijan, A., Sutton, E. C., Snyder, L. E., et al. 2004, ApJ, 606, 917
  • Sánchez-Monge et al. (2010) Sánchez-Monge, Á., Palau, A., Estalella, R., et al. 2010, ApJ, 721, L107
  • Sánchez-Monge et al. (2013a) Sánchez-Monge, Á., Beltrán, M. T., Cesaroni, R., et al. 2013a, A&A, 550, A21
  • Sánchez-Monge et al. (2013b) Sánchez-Monge, Á., Cesaroni, R., Beltrán, M. T., et al. 2013b, A&A, 552, L10
  • Sánchez-Monge et al. (2013c) Sánchez-Monge, Á., Palau, A., Fontani, F., et al. 2013c, MNRAS, 432, 3288
  • Sánchez-Monge et al. (2013d) Sánchez-Monge, Á., López-Sepulcre, A., Cesaroni, R., et al. 2013d, A&A, 557, A94
  • Schilke et al. (2001) Schilke, P., Benford, D. J., Hunter, T. R., Lis, D. C., & Phillips, T. G. 2001, ApJS, 132, 281
  • Schreyer et al. (2002) Schreyer, K., Henning, T., van der Tak, F. F. S., Boonman, A. M. S., & van Dishoeck, E. F. 2002, A&A, 394, 561
  • Schuller et al. (2009) Schuller, F., Menten, K. M., Contreras, Y., et al. 2009, A&A, 504, 415
  • Simon et al. (2000) Simon, M., Dutrey, A., & Guilloteau, S. 2000, ApJ, 545, 1034
  • Sugiyama et al. (2008) Sugiyama, K., Fujisawa, K., Doi, A., et al. 2008, PASJ, 60, 23
  • Surcis et al. (2012) Surcis, G., Vlemmings, W. H. T., van Langevelde, H. J., & Hutawarakorn Kramer, B. 2012, A&A, 541, A47
  • Tan et al. (2014) Tan, J. C., Beltran, M. T., Caselli, P., et al. 2014, arXiv:1402.0919
  • Vallée & Bastien (2000) Vallée, J. P., & Bastien, P. 2000, ApJ, 530, 806
  • Velusamy et al. (2008) Velusamy, T., Marsh, K. A., Beichman, C. A., Backus, C. R., & Thompson, T. J. 2008, AJ, 136, 197
  • Vlemmings (2008) Vlemmings, W. H. T. 2008, A&A, 484, 773
  • Walther et al. (1990) Walther, D. M., Aspin, C., & McLean, I. S. 1990, ApJ, 356, 544
  • Wang K.-S. et al. (2010) Wang, K.-S., Kuan, Y.-J., Liu, S.-Y., & Charnley, S. B. 2010, ApJ, 713, 1192
  • Wang K. et al. (2011) Wang, K., Zhang, Q., Wu, Y., & Zhang, H. 2011, ApJ, 735, 64
  • Wang K.-S. et al. (2012) Wang, K.-S., van der Tak, F. F. S., & Hogerheijde, M. R. 2012, A&A, 543, A22
  • Wang K.-S. et al. (2013) Wang, K.-S., Bourke, T. L., Hogerheijde, M. R., et al. 2013, A&A, 558, A69
  • Wang K. et al. (2014) Wang, K., Zhang, Q., Testi, L., et al. 2014, MNRAS, 439, 3275
  • Williams et al. (1994) Williams, J. P., de Geus, E. J., & Blitz, L. 1994, ApJ, 428, 693
  • Wilson & Rood (1994) Wilson, T. L., & Rood, R. 1994, ARA&A, 32, 191
  • Wright et al. (2010) Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868
  • Zernickel et al. (2012) Zernickel, A., Schilke, P., Schmiedeke, A., et al. 2012, A&A, 546, A87
  • Zhang B. et al. (2009) Zhang, B., Zheng, X. W., Reid, M. J., et al. 2009, ApJ, 693, 419
  • Zhang Y. et al. (2013) Zhang, Y., Tan, J. C., De Buizer, J. M., et al. 2013, ApJ, 767, 58
  • Zinnecker & Yorke (2007) Zinnecker, H., & Yorke, H. W. 2007, ARA&A, 45, 481

Appendix A Hot core intensity and velocity maps

In Figs. 21 and 22, we present the zeroth (integrated intensity) and first (velocity field) moment maps of various hot-core tracers detected towards cores A and B. In Section 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.

Figure 21: Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CHCN =2 at 349.426 GHz, b) CHCN =8 at 349.024 GHz, c) CHCN =2 at 349.254 GHz, and d) CHCN =1 =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 km s) 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.
Figure 22: Maps of the integrated intensity (upper panels) and velocity field (bottom panels) of a) CHOCHO (26–25), b) CHOCHO (27–26), c) HCS (10–9), d) HNCO (16–15), e) SO (18–18), and f) HCN (37–36) towards cores A and B. For the different panels, the starting (and increasing contour levels in terms of , in Jy beam km s) 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.

Appendix B ChOH and CHCN myXCLASS fits

In Figs. 2326 we show the spectra of the bulge of methanol (CHOH, black line in the top panels) and methyl cyanide (CHCN, 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 Section 4.1). The procedure searches, by minimizing the , for the best fit of five parameters: size (in arcsec), temperature (in K), column density (in cm), linewidth (in km s) and LSR velocity (in km s). The blue lines in the small panels show the 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 CHCN and CHOH), thus avoiding possible degeneracies (see Zernickel et al., 2012) between the parameters of the fit (e.g., column density and size).