Issue |
A&A
Volume 604, August 2017
|
|
---|---|---|
Article Number | A60 | |
Number of page(s) | 25 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/201730648 | |
Published online | 07 August 2017 |
Exploring molecular complexity with ALMA (EMoCA): Detection of three new hot cores in Sagittarius B2(N)
1 Max-Planck Institut für Radioastronomie, Auf dem Hügel 69, 53121 Bonn, Germany
e-mail: bonfand@mpifr-bonn.mpg.de
2 Departments of Chemistry and Astronomy, University of Virginia, Charlottesville, VA 22904, USA
3 I. Physikalishes Institut, Universität zu Köln, Zülpicher Str. 77, 50937 Köln, Germany
Received: 17 February 2017
Accepted: 25 March 2017
Context. The Sagittarius B2 molecular cloud contains several sites forming high-mass stars. Sgr B2(N) is one of its main centers of activity. It hosts several compact and ultra-compact HII regions, as well as two known hot molecular cores (Sgr B2(N1) and Sgr B2(N2)) in the early stage of the high-mass star formation process, where complex organic molecules (COMs) are detected in the gas phase.
Aims. Our goal is to use the high sensitivity of the Atacama Large Millimeter/submillimeter Array (ALMA) to characterize the hot core population in Sgr B2(N) and thereby shed new light on the star formation process in this star-forming region.
Methods. We use a complete 3 mm spectral line survey conducted with ALMA to search for faint hot cores in the Sgr B2(N) region. The chemical composition of the detected sources and the column densities are derived by modeling the whole spectra under the assumption of local thermodynamic equilibrium. Population diagrams are constructed to fit rotational temperatures. Integrated intensity maps are produced to derive the peak position and fit the size of each molecule’s emission distribution. The kinematic structure of the hot cores is investigated by analyzing the line wing emission of typical outflow tracers. The H2 column densities are computed from ALMA and SMA continuum emission maps.
Results. We report the discovery of three new hot cores in Sgr B2(N) that we call Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5). The three sources are associated with class II methanol masers, well known tracers of high-mass star formation, and Sgr B2(N5), also with a UCHII region. Their H2 column densities are found to be between approximately 16 and 36 times lower than the one of the main hot core Sgr B2(N1). The spectra of these new hot cores have spectral line densities of 11 up to 31 emission lines per GHz above the 7σ level, assigned to 22–25 molecules plus 13–20 less abundant isotopologs. We derive rotational temperatures of approximately 140–180 K for the three new hot cores and mean source sizes of 0.4″ for Sgr B2(N3) and 1.0″ for Sgr B2(N4) and Sgr B2(N5). The chemical composition of Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) is very similar, but it differs from that of Sgr B2(N2). Finally, Sgr B2(N3) and Sgr B2(N5) show high-velocity wing emission in typical outflow tracers, with a bipolar morphology in their integrated intensity maps suggesting the presence of an outflow, like in Sgr B2(N1). No sign of an outflow is found around Sgr B2(N2) and Sgr B2(N4). We derive statistical lifetimes of 4 × 104 yr for the class II methanol maser phase and 6 × 104 yr for the hot core phase in Sgr B2(N).
Conclusions. The associations of the hot cores with class II methanol masers, outflows, and/or UCHII regions tentatively suggest the following age sequence: Sgr B2(N4), Sgr B2(N3), SgrB2(N5), Sgr B2(N1). The status of Sgr B2(N2) is unclear. It may contain two distinct sources, a UCHII region and a very young hot core.
Key words: stars: formation / ISM: individual objects: Sagittarius B2(N) / astrochemistry / ISM: molecules
© ESO, 2017
1. Introduction
The Sagittarius B2 molecular cloud (Sgr B2 hereafter) is one of the most prominent regions forming high-mass stars in our Galaxy, with a mass of ~10 in a diameter of ~40 pc (Lis & Goldsmith 1990). It is located in the central molecular zone, close to the Galactic center (~100 pc from the central supermassive black hole Sgr A∗ in projection) and at a distance of 8.34 ± 0.16 kpc from the Sun (Reid et al. 2014). The cloud contains several sites of on-going high-mass star formation. Its two major centers of activity are called Sgr B2(M)ain and Sgr B2(N)orth. They both host a cluster of (ultra)compact HII regions (Mehringer et al. 1993; Gaume et al. 1995; De Pree et al. 1998, 2015) as well as class II methanol masers (Caswell 1996); both phenomena provide strong evidence of on-going high-mass star formation. The high density of molecular lines and the continuum emission detected toward this region highlight the presence of a large amount of material to form new stars. Sgr B2(N) contains two dense and compact hot molecular cores, called Sgr B2(N1) and Sgr B2(N2) (Belloche et al. 2016), separated by ~5″ in the north-south direction (corresponding to 0.2 pc in projection, Belloche et al. 2008; Qin et al. 2011). They are both in the early stage of star formation when a high-mass protostar has already formed and has started to warm up its circumstellar envelope, and is producing ionizing radiation that creates an ultracompact HII region around it. The earliest stages of the formation process of high-mass stars are still poorly understood since, first, the high dust content of high-mass protostellar cores makes observations at infrared and shorter wavelengths impossible and, second, because of the small angular sizes of these faraway regions. However, it is in these stages that many molecules are formed, directly in the gas phase or at the surface of dust grains. These molecules can readily be studied by their rotational spectrum and now, with the Atacama Large Millimeter/submillimeter Array (ALMA), at exquisite angular resolution and sensitivity. While the newly ignited protostar increases the temperature of its surroundings, the molecules frozen in the ice mantles of the dust grains are released in the gas phase, explaining the large number of molecules detected toward Sgr B2.
Over the past five decades, nearly 200 molecules have been discovered in the interstellar medium (ISM) or in circumstellar envelopes of evolved stars1. Among them, approximately 60 are composed of six atoms or more and are called Complex Organic Molecules (COMs) in the field of astrochemistry (Herbst & van Dishoeck 2009). Most species have been discovered toward the warm and dense parts of star-forming regions and many of the first detections of interstellar molecules at radio and (sub)millimeter wavelengths were made toward Sgr B2, such as acetic acid CH3COOH (Mehringer et al. 1997), glycolaldehyde CH2(OH)CHO (Hollis et al. 2000), acetamide CH3CONH2 (Hollis et al. 2006), and aminoacetonitrile NH2CH2CN (Belloche et al. 2008); see also the overview by Menten (2004). Sgr B2(N) thus appears to be one of the best targets for studying COMs and searching for new molecules, and thus expanding our view of the chemical complexity of the ISM.
Through the investigation of the chemical composition of Sgr B2(N), we wish to characterize its hot core population to shed new light on the star formation process in this region. To this aim we analyze a spectral line survey recently conducted with ALMA in its cycles 0 and 1 at high angular resolution. We take advantage of the high sensitivity of the EMoCA survey (standing for “Exploring Molecular Complexity with ALMA”) to search for fainter hot cores in Sgr B2(N) in order to extend our view of the distribution of active star-forming regions in this outstanding cloud. The article is structured as follows. The observations and method of analysis are presented in Sect. 2. The results are given in Sect. 3 and discussed in Sect. 4. Finally the conclusions are presented in Sect. 5.
2. Observations and method of analysis
2.1. ALMA
The ALMA observations that provided the data used here target Sgr B2(N) with the phase center located half way between Sgr B2(N1) and Sgr B2(N2), at αJ2000 = 17h47m19.87s, δJ2000 = −28o22′16″. It is a complete spectral line survey between 84.1 GHz and 114.4 GHz, conducted at high angular resolution (~1.6″) and with a high sensitivity ~3 mJy/beam per 488 kHz (1.7 to 1.3 km s-1) wide channel. The spectral line survey is divided into five spectral setups, each one delivering four spectral windows. Details about the different setups and the data reduction are presented in Belloche et al. (2016). The size (HPBW) of the primary beam of the 12 m antennas varies between 69″ at 84 GHz and 51″ at 114 GHz (Remijan 2015).
2.2. Submillimeter Array (SMA)
Qin et al. (2011) observed the Sgr B2 region using the SMA in the compact and very extended configurations, reporting the first high-angular-resolution submillimeter continuum observations of this region. The continuum map used here was obtained at 342.883 GHz, with a synthesized beam of 0.4″× 0.24″ and a position angle of 14.4o. The map has been corrected for the primary beam attenuation. More details can be found in Qin et al. (2011).
2.3. Radiative transfer modeling of the line survey
Given the high densities observed in the Sgr B2(N) region (Belloche et al. 2008, 2014; Qin et al. 2011), it is appropriate to work under the local thermodynamic equilibrium (LTE) approximation. We use Weeds (Maret et al. 2011), which is part of the CLASS software2, to perform the line identification and modeling of the spectra, after correction for the primary beam attenuation. A synthetic spectrum is produced for each species solving the radiative transfer equation and taking into account the finite angular resolution of the interferometer, the line opacity, and line blending. Each molecule is modeled separately adjusting the following parameters: column density, rotational temperature, angular size of the emitting region (assumed to be Gaussian), velocity offset with respect to the systemic velocity of the source, and linewidth (full width at half maximum – FWHM). For each species, population diagrams are plotted to derive the rotational temperature and two-dimensional (2D) Gaussians are fitted to integrated intensity maps to measure the size of the emitting region. The linewidth and velocity offset are derived from one-dimensional (1D)-Gaussian fits to the lines detected in the spectra. All parameters are then adjusted manually until a good fit to the data is obtained. Finally, contributions from all species are added to obtain the complete synthetic spectrum. Through this modeling process, a line is assigned to a given molecule only if all lines from this molecule emitted in the frequency range of the survey are detected with the right intensity predicted by the model and if no line is missing in the observed spectrum. More information about the modeling procedure can be found in Belloche et al. (2016). The spectroscopic predictions used to model the spectra are the same as in Belloche et al. (2016, 2017), and Müller et al. (2016). They originate mainly from the CDMS and JPL catalogs (Endres et al. 2016; Pearson et al. 2010).
3. Results
3.1. Detection of three new sources
The high sensitivity of our ALMA data set allows us to search for fainter hot cores in the vicinity of Sgr B2(N1) and Sgr B2(N2). To this aim, we counted the channels with continuum-subtracted flux densities above the 7σ level (1σ ~ 3 mJy/beam) over the whole frequency range for each pixel in the field of view. In this count we excluded setup 3 which has the lowest angular resolution (HPBW > 2″, see Table 2 in Belloche et al. 2016) and we counted only once the channels located in frequency ranges overlapped by adjacent spectral windows. The analyzed ranges cover 47 296 channels, that is, 23.1 GHz. The result is presented in Fig. 1 as a contour map of the Sgr B2(N) region. On this map, the contours represent the number of channels with flux density above the 7σ threshold, that is, higher contours reflect the presence of more emission lines. The figure clearly shows the two main hot cores Sgr B2(N1) and Sgr B2(N2). It also reveals high spectral line densities toward three other positions, unveiling the presence of three new sources that we call Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5). The region inside the contour located south-east of Sgr B2(N5) shows emission lines only in setups 2 and 5 (4 spectral windows) with line intensities lower than toward Sgr B2(N5) but a line content that is very similar. We think this structure is a deconvolution artifact.
We fit 2D Gaussians to the channel count map with the GILDAS procedure GAUSS-2D in order to derive the peak position of the five sources. The results of the fits are listed in Table 1 and the position of the spectral line density peaks of each source is marked with a blue cross in Fig. 1. These positions are adopted as reference positions of the hot cores. Known HII regions (Gaume et al. 1995; De Pree et al. 2015) and class II methanol masers (Caswell 1996) are shown in Fig. 1 as green and red crosses, respectively. The distance of each hot core to the closest UCHII region and the distance to the closest class II methanol maser are given in Table 1. Sgr B2(N1), Sgr B2(N2), and Sgr B2(N5) are associated with UCHII regions. All the new hot cores, Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) happen to be associated with class II methanol masers.
![]() |
Fig. 1 Contour map of the number of channels with continuum-subtracted flux density above the 7σ level (1σ ~ 3 mJy/beam). The contour levels are: 500, 1000, 2000, 5000, 10 000, 20 000, 30 000, and 40 000. The blue crosses mark the peaks of spectral line density. Red crosses represent the 6.7 GHz class II methanol masers (Caswell 1996). Green crosses represent the compact and ultra-compact HII regions (Gaume et al. 1995; De Pree et al. 2015). The offsets are defined with respect to the phase center (see Sect. 2.1). |
Position and spectral line density of the hot cores embedded in Sgr B2(N), and distance to UCHII regions and class II methanol masers.
3.2. Continuum properties of the new sources
In order to calculate H2 column densities from the dust thermal emission arising from the hot cores embedded in Sgr B2(N) we need to know the dust mass opacity coefficient κν, which depends on the dust emissivity exponent β. To derive these parameters we need observations at at least two frequencies. Here we use the ALMA and SMA data after processing them in the following way: the flux densities measured in the ALMA continuum maps need to be corrected for the contribution of the free-free emission and the SMA continuum map has to be smoothed to the ALMA resolution.
3.2.1. ALMA data: correction for free-free emission
![]() |
Fig. 2 Continuum map of the Sgr B2(N) region obtained with ALMA at 108 GHz. Contour levels (positive in black solid line and negative in dashed line) start at five times the rms noise level, σ, of 3.0 mJy/beam and double in value up to 320σ. The filled ellipse shows the synthesized beam (1.65″× 1.21″, PA = −83.4o). The black cross represents the phase center. The blue crosses mark the positions of the five hot cores embedded in Sgr B2(N) derived from Fig. 1. The red crosses represent the average peak positions of the continuum emission derived from 2D-Gaussian fits to all ALMA continuum maps. The green crosses mark the compact and ultra-compact HII regions. The dotted red circle represents the size (HPBW) of the primary beam of the 12 m antennas at 108 GHz. The map is not corrected for the primary beam attenuation. |
We investigate 20 continuum maps obtained at different frequencies over the whole frequency range of the ALMA survey. As an example, Fig. 2 shows the continuum map obtained at 108 GHz with ALMA. Only four of the five hot cores embedded in Sgr B2(N) are detected above the 5σ level in this map. Over the whole frequency range, Sgr B2(N1) and Sgr B2(N2) show a strong continuum emission while Sgr B2(N4) and Sgr B2(N5) appear clearly weaker. Sgr B2(N3) is not detected at all; for this source we can only derive an upper limit to its H2 column density from the measurement of the noise level in the maps. In this case, noise histograms are plotted from all pixels using the command GO NOISE of the GILDAS software to derive the rms noise level in each map. For the other hot cores we measure the peak flux density (Jy/beam) by fitting a 2D-Gaussian to the continuum maps at different frequencies (using the GAUSS-2D task of the GILDAS software). The derived flux density is then corrected for the primary beam attenuation. We note that the compact continuum source detected to the north of Sgr B2(N2) corresponds to the UCHII region K4. It does not have a high spectral line density (see Fig. 1) and thus does not harbor a hot core.
The average peak position of the continuum emission associated to each hot core, derived from the 2D-Gaussian fit, is marked with a red cross in Fig. 2. These positions are also given in Table 2 as well as the distances to the closest hot core and UCHII region. The offset between the continuum peak position and the hot core reference position is at most one fourth of the beam.
Position and size of the ALMA continuum sources, and distances to hot cores and UCHII regions.
Usually (sub)millimeter continuum emission observed toward Sgr B2(N) is attributed to thermal emission from interstellar dust (Kuan & Snyder 1996), while at wavelength ≥1 cm, the continuum emission is dominated by free-free emission arising from HII regions. However in the ALMA frequency range considered in this paper, the free-free emission may significantly contribute to the 3 mm flux densities. Figure 2 shows that Sgr B2(N1), Sgr B2(N2), and Sgr B2(N5) are associated with UCHII regions. Furthermore the overall shape of the extended continuum emission detected with ALMA around Sgr B2(N1) and Sgr B2(N2) is similar to the shape of the 1.3 cm free-free emission reported by Gaume et al. (1995). The peak flux densities measured at 3 mm toward the hot cores need thus to be corrected for the contribution of the free-free emission. Details about this analysis are presented in Appendix A. The fraction of free-free emission estimated toward Sgr B2(N1) varies between ~17% at 85 GHz and ~9% at 114 GHz. The free-free contribution is higher toward Sgr B2(N2), from ~75% at 85 GHz, down to ~32% at 113.5 GHz, and it varies between 40% and 60% toward Sgr B2(N5). The flux densities measured in the 20 ALMA continuum maps and corrected for the free-free contribution are presented in Tables B.1–B.5. We do not make any correction for Sgr B2(N3) and Sgr B2(N4) because they are not associated with any known UCHII region.
3.2.2. SMA data: smoothing to the ALMA resolution
We use the SMA map obtained by Qin et al. (2011) to derive the dust emissivity index β in a joint ALMA/SMA analysis. To this aim we first need to smooth the SMA map to the ALMA resolution. We use the task GAUSS-SMOOTH of the GILDAS software. Figure 3 shows as an example the SMA map smoothed to the same resolution as the ALMA map shown in Fig. 2. We proceed in the same way to smooth the SMA map to the 20 different angular resolutions corresponding to the five ALMA setups (20 spectral windows, 4 per setup; see Belloche et al. 2016, for details about the angular resolution).
Figure 3 shows that only the two main hot cores Sgr B2(N1) and Sgr B2(N2) are detected in the SMA map. To derive peak flux densities toward these two sources we fit 2D Gaussians to the smoothed maps. For the other hot cores that are not detected, we measure in each map the noise level inside the polygons plotted in red in Fig. 3. Each square has a dimension of 6″, and an area of approximately twice the beam size is excluded in the middle. The flux densities or rms derived from the 20 maps are listed in Tables B.1–B.5. At 343 GHz, the contribution of the free-free emission to the flux density is lower than 0.1% and can safely be ignored (Qin et al. 2011).
3.2.3. SMA data: original resolution
To derive H2 column densities from the SMA map at its original resolution (~0.3″), we measure the peak flux densities toward Sgr B2(N1) and Sgr B2(N2) with 2D-Gaussian fits to the map. For the fainter hot cores Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) not detected in the SMA map, we measure the average noise level inside a polygon of 6′′2 excluding an area of approximately twice the beam size in the middle.
![]() |
Fig. 3 Continuum map of the Sgr B2(N) region obtained with the SMA at 343 GHz (Qin et al. 2011) and smoothed to the ALMA resolution. The new beam is shown in the bottom left corner (1.65″× 1.21″, PA = −83.4o). The contour levels (positive in solid line, negative in dashed line) start at 8σ and double up to 128σ with σ = 0.15 Jy/beam, the noise level measured inside the polygon defined around Sgr B2(N3). The black cross represents the SMA phase center. The red crosses mark the peak positions of the ALMA continuum emission derived from Fig. 2. The position of Sgr B2(N3) (blue cross) is derived from Fig. 1. The red dashed circle represents the size (HPBW) of the primary beam of the SMA 6 m antennas at 343 GHz. The map is corrected for the primary beam attenuation. |
3.2.4. H2 column densities
We assume a dust temperature Td ~ Tgas ~ 150 K (see Sect. 3.6). Neglecting the cosmic microwave background temperature, the radiative transfer equation can be written as follows: (1)with
the solid angle of the synthesized beam, Bν(Td) the Planck function at the dust temperature, and τ the dust opacity.
is the peak flux density measured in the continuum maps in Jy/beam.
is corrected for the primary beam attenuation and, in the case of the ALMA data, for the free-free contribution.
From this equation we can calculate the dust opacity as follows: (2)The results of the opacity calculations for both the ALMA and SMA data are summarized in Table 3. The emission is optically thin for all cores except Sgr B2(N1) at the SMA frequency. In the SMA map at its original resolution (0.3″), the dust emission of Sgr B2(N1) is optically thick and inconsistent with a temperature of 150 K. From Eq. (1) and assuming optically thick emission, we derive a lower limit to the dust temperature of 200 K at a scale of 0.3″.
The dust opacity, τ, is related to the H2 column density, NH2, via: (3)with μH2 = 2.8 the mean molecular weight per hydrogen molecule (Kauffmann et al. 2008), mH the mass of atomic hydrogen, and κν the dust mass opacity (in cm2 g-1) given by the power law:
(4)with χd = 100 the standard gas-to-dust ratio, β the dust emissivity exponent, and κ0 the mass absorption coefficient at frequency ν0.
We derive the dust emissivity index β from the ALMA and SMA data. From Eqs. (1), (3), and (4) we can write: (5)and
(6)which can also be written as:
(7)and
(8)Then β is given by:
(9)which only depends on the peak flux density measured on the continuum maps, the beam solid angle, the Planck function, and the frequency. We perform the calculations using the flux densities measured toward Sgr B2(N1) and Sgr B2(N2) in the ALMA maps and in the SMA maps smoothed to the ALMA resolution, excluding setup 3 (low angular resolution). We obtain β = 0.8 ± 0.1 and β = 1.2 ± 0.1 for Sgr B2(N1) and Sgr B2(N2), respectively.
As mentioned above, the continuum emission detected with SMA toward Sgr B2(N1) is very optically thick. Even if the emission in the SMA map smoothed to the ALMA resolution should be less optically thick and our analysis does not assume optically thin emission, it is likely that the value derived for β toward Sgr B2(N1) is underestimated due to these high opacities. It could also be that we have underestimated the contribution of free-free emission toward Sgr B2(N1) at the ALMA frequencies. Therefore, we consider the value obtained toward Sgr B2(N2) as more reliable, and we adopt β = 1.2 for all sources. Such a flat index suggests a dust opacity spectrum intermediate between the models of dust grains without ice mantles and those with thin ice mantles of Ossenkopf & Henning (1994, see their Fig. 5). We adopt a dust mass absorption coefficient κ0 = 1.6 cm2 g-1 (of dust) at λ0 = 1.3 mm, which is intermediate between these models. With this normalization, we obtain a methanol abundance of ~2 × 10-5 relative to H2 for Sgr B2(N2) (see Sect. 3.8), consistent with the peak gas-phase abundance of methanol predicted by our chemical models (Garrod et al. 2009; Müller et al. 2016). A higher value of κ0 would imply lower H2 column densities and, in turn, higher methanol abundances that would no longer be realistic.
Dust opacity of the five hot cores for a dust temperature of 150 K.
From Eqs. (1) and (3) we can now calculate the H2 column density for each hot core using the following equation: (10)Figure 4 shows the results of the H2 column density calculations at Td = 150 K, from both the SMA and ALMA data as a function of ALMA frequency. The figure shows that Sgr B2(N4) is not systematically detected over the whole frequency range of the ALMA survey. All upper limits correspond to five times the noise level. All the results are also listed in Tables B.1–B.5. One should keep in mind that the beam is slightly different as a function of frequency, especially in the frequency range covered by setup 3 which has the lowest angular resolution (HPBW > 2″). We thus calculate the average peak H2 column density of each hot core from both the SMA (smoothed to the ALMA resolution) and ALMA data excluding setup 3. Table 4 summarizes the results obtained from the ALMA data before and after correction for the free-free contribution and for the SMA data after and before smoothing to the ALMA resolution. This table shows that the H2 column densities of Sgr B2(N2), Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) are, respectively, 8, >36, 28, and 16 times lower than that of the main hot core Sgr B2(N1). Within the uncertainties, the column densities or upper limits obtained from the ALMA data for the faint hot cores Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) are consistent with the upper limits derived from the SMA smoothed maps.
![]() |
Fig. 4 H2 column densities as a function of frequency for the five hot cores embedded in Sgr B2(N). The values obtained based on the ALMA data before and after correction for the free-free contribution are shown as black and red crosses respectively. The green crosses represent the results obtained based on the SMA map at 343 GHz smoothed to the angular resolution of the ALMA maps. Error bars are calculated from the error on |
H2 column densities for a dust temperature of 150 K.
3.3. Line identification
![]() |
Fig. 5 Part of the continuum-subtracted ALMA spectra observed toward the hot cores embedded in Sgr B2(N). The spectra have been corrected for the primary beam attenuation, and shifted along the y axis for display purposes. The spectra of Sgr B2(N1) and Sgr B2(N2) have been divided by 10 and 6, respectively. The frequency axis corresponds to the systemic velocities derived in Sect. 3.7. |
Table 1 clearly shows that the spectral line density is much lower toward Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) than for Sgr B2(N1) and Sgr B2(N2), considerably reducing the occurrence of line blending. This is also illustrated in Fig. 5 where a portion of the ALMA spectrum of each hot core is shown.
Statistics of the lines detected toward the three new hot cores.
We use Weeds as described in Sect. 2.3 to perform the line identification and model the spectra observed toward the three new hot cores in order to derive their chemical composition. Table 5 shows the total number of emission lines detected above the 7σ level, counted manually toward the peak position of each hot core, excluding setup 3. It also summarizes the number of species identified so far and the fraction of remaining unidentified lines above 7σ (U-lines hereafter). The line density derived here is somewhat lower than reported in Table 1. In Table 1, we divided the number of channels derived from Fig. 1 by the typical FWHM of an emission line. This was a rough estimate because a faint line may have only a single channel emitting above the 7σ threshold while a strong line will have more channels above 7σ than the number of channels covered by its FWHM. This suggests that Sgr B2(N1) and Sgr B2(N2) probably contain less emission lines than indicated in Table 1.
Table 5 also shows that Sgr B2(N4) has a low spectral line density compared to the other hot cores although roughly the same number of molecules were identified. This difference can be explained by the low number of less abundant isotopologs and vibrationally excited states detected in Sgr B2(N4) compared to Sgr B2(N3) and Sgr B2(N5).
![]() |
Fig. 6 Spatial distribution of the molecules identified around Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5). In each panel, the red star represents the position of the hot core derived from Fig. 1. Blue crosses represent the mean peak position for each species. Error bars correspond to the standard deviation weighted by the uncertainties given by the GAUSS-2D routine. The radius of the dashed circle is given in the top right corner. The average peak position of HC3N, ν7 = 1 is shown in green. |
![]() |
Fig. 7 Integrated intensity maps of selected transitions toward Sgr B2(N3) (top row), Sgr B2(N4) (middle row), and Sgr B2(N5) (bottom row). For each map, the red cross shows the peak position of the emission and the black cross is the reference position of the hot core derived from Fig. 1. The blue ellipse represents the result of the Gaussian fit to the map while the red ellipse is the deconvolved emission size. The black filled ellipse represents the synthesized beam. The rms and contour levels are indicated in Table 6. |
3.4. Spatial distribution of the molecules
In order to get information about the spatial structure of the three new hot cores, we investigate for each identified molecule the integrated intensity maps produced from its vibrational ground state transitions that are well detected and free of contamination from other species. The position of the emission peak of each spectral line is derived from a 2D-Gaussian fit to its map using the GAUSS-2D task of the GILDAS software. Figure 6 shows the spatial distribution of the molecules identified so far toward Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5). Each cross represents the mean peak position of a given species. Figure 6a shows that almost all species identified toward Sgr B2(N3) peak within a distance of ~0.3″ from the reference position of the hot core. Only SO, CH3CCH, and HC3N peak beyond 0.5″ south-west of this position. In the case of Sgr B2(N4) (Fig. 6b) and Sgr B2(N5) (Fig. 6c), all molecules peak within distances of approximately 0.6″ and 0.4″ from the hot core position, respectively. For both Sgr B2(N4) and Sgr B2(N5), HC3N in its vibrational ground state shows broad lines contaminated by other species. For this reason we used transitions in its first vibrationally excited state, ν7 = 1. They peak close to the reference position of each hot core, in particular toward Sgr B2(N3) for which the emission peak of the vibrational ground state transitions is offset to the south-west by ~0.6″.
3.5. Size of the hot cores
For each hot core, we select the transitions that are well reproduced by the model, are not severely contaminated by other species, and have a high signal-to-noise ratio (typically ≥8 up to ~99) to fit 2D Gaussians to their integrated intensity maps and derive their emission size.
It is difficult to constrain the size of Sgr B2(N3) because the emission is unresolved for most species. Only four transitions of OCS and C2H5CN show spatially resolved emission. The first row of Fig. 7 shows their integrated intensity maps. The result of the Gaussian fit is displayed in blue and the red ellipse shows the deconvolved emission size. The middle row of Fig. 7 shows the maps of four transitions for which the emission of Sgr B2(N4) is resolved. In total for this source, 12 transitions from seven distinct species show resolved emission with a strong signal-to-noise ratio (≥8 up to ~16). Finally, the bottom row of Fig. 7 shows four maps produced for Sgr B2(N5). In total for this hot core, 22 transitions from six distinct species show resolved emission with a signal-to-noise ratio from ~10 to 99. We decided to focus only on species showing compact emission. For this reason we did not take into account the transitions of CH3CCH detected toward Sgr B2(N5), also spatially resolved but showing extended emission, with a deconvolved size of ~3.6″, around the position of the hot core.
Table C.1 gives the results of the Gaussian fits to the integrated intensity maps of the lines discussed above that show resolved emission. These results are also plotted in Figs. D.1–D.3 for Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5), respectively. For each transition, the deconvolved major and minor diameters of the emission (θmaj and θmin), given in Cols. 10 and 11 of Table C.1, allow us to calculate the average deconvolved size of the emitting region (). The results are given in Col. 13 of Table C.1. The mean deconvolved size of each hot core is derived from these values. We obtain sizes of 0.4 ± 0.1″ for Sgr B2(N3), 1.0 ± 0.3″ for Sgr B2(N4), and 1.0 ± 0.4″ for Sgr B2(N5).
3.6. Rotational temperature
Population diagrams are plotted based on transitions that are well detected and not severely contaminated by lines from other species to derive rotational temperatures. We use the following equation (Snyder et al. 2005) to compute ordinate values (upper level column density divided by statistical weight, Nu/gu), assuming optically thin emission in first approximation: (11)where kb is the Boltzmann constant, ν the frequency, W the integrated intensity in brightness temperature scale, h the Planck constant, c the speed of light, Aul the Einstein coefficient for spontaneous emission,
the beam filling factor, gu the statistical weight of the upper level, Eu the upper level energy, Trot the rotational temperature, Z the partition function, and Nmol the molecular column density.
We present here the results of this analysis for ethanol and methyl formate, two species that have many well detected lines spread over a large energy range (see Table 10). To derive the rotational temperature of CH3OCHO, we use both its ground and first vibrationally excited states, modeled with the same parameters (Nmol, Trot, source size, voff, Δv). Only the ground state transitions are used for ethanol. Figures 8a, c, 9a, c, and 10a, c show the population diagrams of both species for Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5), respectively. In Figs. 8b, d, 9b, d, and 10b, d we applied to both the observed and synthetic populations the opacity correction factor (see Goldsmith & Langer 1999; Snyder et al. 2005) using the opacities of our radiative transfer model. In all cases the opacity correction only slightly affects the population diagram (see values of τmax in Col. 5 of Table 7). Since a transition can be partially contaminated by other species, we also subtracted from the measured integrated intensities the contribution of contaminating molecules using our model that includes all species identified so far. The synthetic and observed points are closer to one another after this correction, however the synthetic data points are globally below the observed ones. This can be explained by our modeling procedure in which we try not to overestimate the observed peak flux density of each spectral line. The modeled spectrum can then well reproduce the observed spectrum in terms of peak flux density, but it will not necessarily exactly fit the whole line profile because we use a single linewidth to model all detected lines. In addition, one must keep in mind that even after removing the contamination from other species, the measured integrated intensities can still be affected by residual contamination from U-lines.
The synthetic data points are not affected by contamination from other species and should be strictly aligned. The residual dispersion of the synthetic data points seen in Figs. 8b, d, 9b, d, and 10b, d can be explained by the frequency boundaries set to integrate the intensity which are a compromise between covering the line as much as possible and limiting the contamination from other species emitting at nearby frequencies in the observed spectrum.
After correction for the opacity and contamination from other species, the observed data points are roughly aligned and can be fitted by a single straight line, meaning one single temperature component. The blue line in Figs. 8b, d, 9b, d, and 10b, d shows the weighted linear fit to the observed data points. According to Eq. (11) and considering the LTE approximation, the inverse of the slope gives us the kinetic temperature of the hot cores. The results are presented in Table 7 and show that all hot cores have temperatures around 140 K, also consistent with the kinetic temperature of Sgr B2(N2) (Müller et al. 2016; Belloche et al. 2016).
For all the species emitting numerous lines spread over a broad energy range, we use the same method to plot population diagrams and derive rotational temperatures. For the rest of the molecules, we set the temperature to 145 K in our LTE model, which allows in most cases to reproduce well the observed spectra of the three new hot cores.
![]() |
Fig. 8 Population diagrams of C2H5OH (panels a and b) and CH3OCHO (panels c and d) for Sgr B2(N3). The black points are computed using the integrated intensities of the observed spectrum while the red points are computed using the integrated intensities of our synthetic model. The error bars on the observed data points are 1σ uncertainties on Nu/gu. No correction is applied in panels a and c, while in panels b and d the optical depth correction has been applied to both the observed and synthetic populations and the contamination from all other species included in the full LTE model has been removed from the observed data points. The blue line is the weighted linear fit to the observed populations. |
3.7. Kinematic structure of the hot cores
We derived the systemic velocity of the faint hot cores embedded in Sgr B2(N) by fitting 1D Gaussians to numerous well detected transitions over the whole spectrum. We obtained a velocity of 74 km s-1 for Sgr B2(N3), 64 km s-1 for Sgr B2(N4), and 60 km s-1 for Sgr B2(N5). Sgr B2(N3) and Sgr B2(N4) have the same velocity as Sgr B2(N2) and Sgr B2(N1), respectively.
These fits also give information on the typical linewidth of each species. The three hot cores have a median linewidth of approximately 5.0 km s-1, ranging from 4.0 to 7.7 km s-1 for Sgr B2(N3), from 3.5 to 5.5 km s-1 for Sgr B2(N4), and from 3.5 to 6.8 km s-1 for Sgr B2(N5) (see Table 10).
Rotational temperatures derived from population diagrams.
The spectra observed toward the three new hot cores show some broader lines with wing emission at high velocities which could suggest the presence of outflows. Investigating these molecular outflows can provide information on the evolutionary stage during the star formation process. To assess the presence of outflows and characterize their properties, we investigate spectral lines from molecules known as typical tracers of outflows. Figure 11 presents the spectra of some of these lines detected toward Sgr B2(N3) and Sgr B2(N5). All of them show broad wing emission at blue- and red-shifted velocities compared to the systemic velocity of the source. The top row of Fig. 11 shows two transitions of SO and the OCS(8–7) transition observed toward Sgr B2(N3). The blue wing of the latter appears narrower because its boundaries have been defined in order to avoid contamination from another species (see green spectrum in Fig. 11). The bottom row of Fig. 11 presents the CS(2–1) transition and three transitions of SO observed toward Sgr B2(N5). The SO transition with the lowest upper energy level, SO(23–12), shows a deep absorption profile and strong broad wings in emission, like the CS(2–1) transition. For all lines, the wing boundaries are chosen to avoid contamination from the line core emission predicted by the LTE model (magenta spectrum in Fig. 11) and contamination from other species (green spectrum in Fig. 11). The velocity ranges used to integrate the blue- and red-shifted emission are summarized in Cols. 6 and 7 of Table 8.
The integrated intensity maps of blue- and red-shifted emission are presented in Fig. 12 for each line along with the map of the line core, or the continuum emission when the line core is in absorption. All maps show a bipolar morphology for both Sgr B2(N3) and Sgr B2(N5) with distinct blue and red lobes, shifted compared to the line core, which could suggest the presence of an outflow. The top row of Fig. 12 shows the integrated intensity maps of the SO lines observed toward Sgr B2(N3). Both maps are similar, showing the peak position of the blue wing clearly shifted north-east of the line core, while the red one is slightly shifted to the south. The map produced for the OCS(8–7) transition confirms the north-south velocity gradient observed in the SO maps. The bottom row of Fig. 12 presents the maps of the line wings observed toward Sgr B2(N5). The maps of the SO(23–12) and CS(2–1) transitions, with the line core seen in absorption, show a clear bipolar structure oriented NE-SW. The other two transitions of SO exhibit a different morphology, with the red lobe shifted east of the continuum emission. The reason for this behavior is unclear but we note that the velocities over which the red-shifted wing emission is integrated are smaller than the ones used for SO(23–12) and CS(2–1) (see Table 8).
Properties of molecular outflows toward Sgr B2(N3) and Sgr B2(N5).
![]() |
Fig. 11 Spectra of the lines investigated to search for outflows toward Sgr B2(N3) (top row) and Sgr B2(N5) (bottom row). The dashed vertical line marks the systemic velocity of the source and the high velocity wings are highlighted in blue and red. The magenta spectrum represents our best-fit model while the green spectrum shows the model that contains all the identified species. The rest frequency and upper level energy (in temperature units) of each transition are indicated in each panel. |
![]() |
Fig. 12 Integrated intensity maps of the lines shown in Fig. 11 toward Sgr B2(N3) (top row) and Sgr B2(N5) (bottom row). In each panel, maps of the blue- and red-shifted wings are presented in blue and red contours, respectively, overlaid on the integrated emission of the line core (black contours). For the line cores affected by absorption, the black contours represent the continuum emission. Rms noise levels and contour levels used for each map are listed in Table 9. Each cross corresponds to the peak position of the emission. The black triangle marks the position of the hot core (Sgr B2(N3) or Sgr B2(N5)) derived from Fig. 1. |
For each transtion the distance r between the peak positions of the blue- and red-shifted outflow lobes and the hot core reference positions are given in Table 8. In Cols. 8 and 9 we calculated the maximum outflow velocities, Vmax, for the blue and red lobes as the difference between the high end of the velocity range set to integrate the wing emission and the systemic velocity of the source. From these values we calculate dynamical timescale of each outflow lobe as assuming that the inclination of the outflow axis with respect to the line of sight is approximately 45o as the maps show two distinct blue and red lobes. The average dynamical times obtained from the blue- and red-shifted lobes are on the order of a few thousand years (Table 8).
Finally, we found no evidence for a bipolar structure around Sgr B2(N4) and Sgr B2(N2). In the case of Sgr B2(N4), the same lines as those investigated in Sgr B2(N3) and Sgr B2(N5) are too weak to show wing emission (see SO lines in Fig. D.4) or they show atypical shapes with blue-shifted wings but without red-shifted wings (see OCS lines in Fig. D.4). The same shape is also observed in the lines of other species not considered as typical tracers of outflows (see Fig. D.5). It is thus difficult to conclude whether this broad blue-shifted emission in Sgr B2(N4) is due to an outflow or to a fainter component emitting at lower velocity.
3.8. Chemical composition
We use Weeds as described in Sect. 2.3 to model the emission lines of ten (complex) organic molecules detected toward the three new hot cores. For each molecule, the source size, rotational temperature, and velocity and linewidth are derived as described in Sects. 3.5–3.7, respectively. All these parameters are listed in Table 10 along with the molecular column densities obtained from our best-fit LTE models toward the three hot cores. Column densities derived for Sgr B2(N2) are also shown for comparison (Müller et al. 2016; Belloche et al. 2016, 2017). We investigated the isotopologs 13CH3CN and CHCN instead of CH3CN because the vibrational ground state transitions of the latter are optically thick. We assume the isotopic ratio [CH3CN]/[13CH3CN] = [CH3CN]/[CH
CN] = 21 derived by Belloche et al. (2016) to obtain the column density of CH3CN. The resulting chemical composition of the four hot cores is displayed in Fig. 13a. Figure 13b shows the column densities normalized to the column density of C2H5CN.
We compute chemical abundances relative to H2 using the H2 column densities derived in Sect. 3.2. From the emission lines detected toward Sgr B2(N3), we derived a deconvolved source size of 0.4″ (see Sect. 3.5), therefore we use the upper limit on the H2 column density derived from the SMA map at its original resolution of 0.3″ (Table 4). This gives us lower limits to the chemical abundances of the molecules detected toward Sgr B2(N3). In the case of Sgr B2(N2) and Sgr B2(N5) we derived deconvolved source sizes >2.0″ from the ALMA continuum maps (see Table 2), larger than the ALMA angular resolution of 1.6″. For both hot cores, we thus extrapolate the H2 column density obtained at the ALMA angular resolution (Table 4) to the more compact region where the molecular emission comes from (~1.2″ and 1.0″ for Sgr B2(N2) and Sgr B2(N5) respectively, see Sect. 3.5). We assume spherical symmetry and a density profile proportional to r-1.5, which implies a column density scaling with r-0.5. Finally, the ALMA continuum maps yield an average deconvolved source size of 0.7″ for Sgr B2(N4), smaller than the ALMA resolution. Therefore, we correct the H2 column density obtained for Sgr B2(N4) in Sect. 3.2.4 for the beam dilution, that is, we multiply it by . The resulting H2 column densities are listed in Table 11 and used in Fig. 13c to derive chemical abundances.
Parameters of our best-fit LTE model.
Figure 13 shows that the three new hot cores have similar compositions but differ from Sgr B2(N2). Among the new hot cores, Sgr B2(N3) appears to be closer to Sgr B2(N5) than Sgr B2(N4) in terms of chemical content. Figure 13b shows indeed that Sgr B2(N3) and Sgr B2(N5) have, within a factor of 2, the same abundances relative to C2H5CN. It is also the case for Sgr B2(N4), except for CH3NCO and CH3SH, which are more abundant in this source. NH2CHO is not detected toward Sgr B2(N4). Relative to H2, Sgr B2(N4) has COM abundances roughly one order of magnitude below that of Sgr B2(N5), while Sgr B2(N2) and Sgr B2(N3) lie roughly one order of magnitude above (Fig. 13c). C2H5CN and NH2CHO are in particular much more abundant in Sgr B2(N2).
4. Discussion
4.1. Discovery of three new hot cores in Sgr B2(N)
In Sect. 3.1 we presented the detection of three new hot cores in Sgr B2(N), which we call Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5), based on the detection of high spectral line density regions in Fig. 1. While Sgr B2(N3) is not detected in the ALMA continuum maps, Sgr B2(N4), and Sgr B2(N5) show a faint continuum level with H2 column densities approximately 28 and 16 times lower than that calculated for Sgr B2(N1). The two main hot cores, Sgr B2(N1) and Sgr B2(N2), have high H2 column densities of 1.3 × 1025 cm-2 and 1.6 × 1024 cm-2 at a resolution of ~1.6″. In a previous analysis of Sgr B2(N) based on the same ALMA survey, Belloche et al. (2014) calculated a peak H2 column density at 98.8 GHz of 4.2 × 1024 cm-2 for Sgr B2(N2), for a beam size of 1.8″ × 1.3″, which is approximately three times higher than our result. This difference results from the dust mass opacity coefficient they assumed (κ98.8GHz = 4.3×10-3 cm2 g-1), ~1.3 times smaller than the value we derive for this frequency based on our combined analysis of the ALMA and SMA data sets. In addition, they assumed that all the continuum emission is due to dust while we took into account the contribution of the free-free emission. This results in an additional factor ~2 on the H2 column density. From the flux density measured toward Sgr B2(N2) in the SMA map at its original resolution we derived a H2 column density of 5.1 × 1024 cm-2. Based on the same map, Qin et al. (2011) calculated an H2 column density of 1.5 × 1025 cm-2 toward Sgr B2(N2). They used the same dust temperature, but they assumed optically thin emission and a dust mass opacity coefficient κ343 GHz = 6.8 × 10-3 cm2 g-1 that is 3.8 times smaller than the value we obtain at this frequency from our combined ALMA/SMA analysis. In addition, we derived a dust opacity of 0.61 for Sgr B2(N2) in the original SMA map (see Table 3), which implies that they underestimated the column density by a factor ~1.3. The two effects explain the factor ~3 difference between the two studies.
![]() |
Fig. 13 a) Column densities of various molecules detected toward four of the five hot cores embedded in Sgr B2(N) (see Table 10). b) Column densities normalized to the column density of C2H5CN. c) Chemical abundances with respect to H2. In each panel, lower and upper limits are indicated with arrows. |
Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) are associated with 6.7 GHz methanol masers (see Fig. 1). Indeed, Col. 8 of Table 1 shows that within the uncertainties, the peak positions of line density of the three hot cores are consistent with the positions of the class II methanol masers reported in Sgr B2(N) by Caswell (1996). The 6.7 GHz methanol transition is the strongest and most widespread of the class II methanol masers (Menten 1993). It is known to be one of the best tracers of star formation and thought to be associated exclusively with regions forming high-mass stars (Minier et al. 2003; Xu et al. 2008). Indeed, for a methanol transition to exhibit maser emission, both suitable physical conditions and a sufficient abundance of methanol are required. Sobolev et al. (1997) showed that pumping requires dust temperature >150 K, high methanol column densities (>2 × 1015 cm-2), and moderate densities (nH< 108 cm-3) to excite the 6.7 GHz methanol maser transition. For these reasons, low-mass stars are not expected to produce class II methanol masers and indeed such masers have not been detected toward regions forming low-mass stars so far (see, e.g., Pandian et al. 2008). The specific conditions for masers to exist make them powerful probes of high-mass star formation sites and confirm the nature of the new hot cores discovered in Sgr B2(N).
Methanol masers can also allow us to assess the evolutionary stages of these new sources as they are thought to trace evolutionary stages from the IR dark cloud (Pillai et al. 2006) to the UCHII phase. Urquhart et al. (2014) proposed an evolutionary sequence for high-mass star formation in which methanol masers and HII regions trace two different phases, with the masers probing an earlier stage in the high-mass star formation process. Walsh et al. (1997) showed that methanol maser emission is detectable before radio continuum emission, that is, before the formation of UCHII regions. The masers are thought to be associated with deeply embedded high-mass protostars not evolved enough to ionize the surrounding gas and produce a detectable HII region. The maser emission remains active after first, a hyper and then an ultra compact HII region has formed around the star for a significant portion of the UCHII phase, during which the methanol molecules are shielded from the central UV radiation by the warm dust in the UCHII region’s slowly expanding molecular envelope whose emission also provides its mid-IR pumping photons. Finally, the maser emission stops as the UCHII region expands. This picture is supported by recent radio observations with the greatly increased sensitivity of the newly expanded Karl G. Jansky Very Large Array (JVLA). Hu et al. (2016) find radio continuum emission in the vicinity of a third of their sample of 372 methanol masers. This is a significantly higher percentage than found by Walsh et al. (1998) with the Australia Compact Array (ATCA), which is due to the fact that the new JVLA images are ~20 times deeper; that is, they have a typical rms noise level 45 μJy beam-1 at 4–8 GHz compared to the ~1 mJy beam-1 (at 8.64 GHz) of the Walsh et al. ATCA data. Most likely, at the JVLA’s higher sensitivity, one is able to detect the emission from hypercompact HII regions surrounding the still accreting protostar exciting the masers (Keto 2007). Due to such regions’ compactness, their radio emission is very weak. However, an increasing number of such objects is now being detected with the JVLA (Rosero et al. 2016; Hu et al. 2016). van der Walt (2005) estimated the lifetime of class II methanol masers at between 2.5 × 104 and 4.5 × 104 yr (depending on the assumed IMF). After this period, the HII region then exists without the maser emission. These last considerations suggest that Sgr B2(N1) and Sgr B2(N2) are already more evolved than the new hot cores for which the maser emission still exists. Among them, Sgr B2(N5) is associated with both a UCHII region and a class II methanol maser, which might suggest that it is in a phase between Sgr B2(N3)/Sgr B2(N4) and Sgr B2(N1)/Sgr B2(N2), when a UCHII region has formed and coexists with maser emission.
4.2. Chemical composition of the new hot cores
Another way to assess the evolutionary stages of the hot cores is to compare their chemical compositions. Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) have spectral line densities of 31, 11, and 22 lines per GHz above 7σ, respectively (see Sect. 3.3), much lower than the two main hot cores, Sgr B2(N1) and Sgr B2(N2). Our LTE model allowed us to identify approximately 91%, 89%, and 93% of the lines detected above 7σ that have been assigned to 22, 23, and 25 main species, respectively (see Table 5). This is much less than the 52 species identified so far toward Sgr B2(N2) based on the EMoCA survey (Belloche, priv. comm.). Although Sgr B2(N3) has the highest spectral line density among the new hot cores, Sgr B2(N5) contains the largest number of identified molecules. However fewer isotopologs and vibrationally excited states of these molecules are detected toward Sgr B2(N5) than Sgr B2(N3). Sgr B2(N4)’s spectrum shows a low spectral line density and less species have been identified toward this source. Most of the remaining unidentified lines are thought to belong to vibrationally excited states of already known molecules for which the spectroscopic predictions are still missing (Belloche et al. 2013), but the presence of new molecules in the ALMA spectra is not excluded.
In order to assess whether the differences reported above are real or sensitivity limited, we need to compare the chemical content of the three new hot cores in more detail. In Sect. 3.8, we reported the chemical abundances of (complex) organic molecules detected toward the new hot cores plus Sgr B2(N2) (Fig. 13). Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) have a similar chemical composition but differ from Sgr B2(N2). According to our chemical model presented in Fig. 3 of Belloche et al. (2014), the chemical abundance of C2H5CN in the gas phase after sublimation is expected to decrease with time while the abundance of C2H3CN increases. The abundance of CH3CN is also expected to increase with time after sublimation. Table 12 and Fig. 13 show that Sgr B2(N2) has the lowest ratios of [C2H3CN]/[C2H5CN] and [CH3CN]/[C2H5CN] among the four hot cores, which could suggest that it is less evolved than the three new hot cores. According to our chemical model presented in Fig. 5 of Belloche et al. (2017), the chemical abundances of CH3NCO and HNCO in the gas phase are also expected to increase with time after sublimation. Sgr B2(N2) also shows low [CH3NCO]/[C2H5CN] and [HNCO]/[C2H5CN] ratios, consistent with the hypothesis made above. However, all these ratios involve C2H5CN in the denominator, which may dominate the general trend discussed above. A deeper analysis involving ratios of various molecules will be necessary to constrain the evolutionary stages of Sgr B2(N)’s hot cores from their chemistry in a more reliable way, especially because the conclusion drawn here that Sgr B2(N2) appears chemically younger than the three new hot cores is in contradiction with the conclusion drawn in Sect. 4.1 from the associations of the hot cores with methanol masers and/or UCHII regions (see discussion in Sect. 4.4, however).
One should also note that the chemical timescales of molecules within any single model may not relate directly to the dynamical age of the source, as pointed out by Garrod et al. (2008). This may be exacerbated by the physical differences between the sources, which will also affect the chemical timescales. Only a more explicit and individualized treatment of the time- and space-dependent chemistry in each source is likely to give an accurate explanation for the observed chemical differences.
H2 column densities used to derive chemical abundances.
Column density ratios of selected molecules.
4.3. Spatial distribution of the molecules toward the new hot cores
In Sect. 3.4 we investigated the spatial distribution of the molecules detected toward the three new hot cores. Most of the identified species show compact emission with emission peaks within distances of ~0.3″, 0.6″, and 0.4″ of the positions of Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5), respectively (see Fig. 6). This led us to define the reference position of each hot core as the position where its spectral line density peaks (Sect. 3.1). Within the uncertainties the emission peak of most species is consistent with this position. Only three molecules, SO, HC3N, and CH3CCH, show significant offsets between their emission peaks and the position of Sgr B2(N3). The transitions of SO detected toward Sgr B2(N3) have an emission peak located, on average, ~0.5″ south-west of the hot core (Fig. 6a and top row of Fig. 12). Figures 14a, b present the integrated intensity maps of two vibrational ground state transitions of CH3CCH and HC3N. They both show the same elongated shape as the SO maps. In both cases, the peak position of the emission given by the 2D-Gaussian fit to the map is shifted beyond 0.5″ from the hot core position due to this extended emission. On the contrary, transitions from within the first vibrationally excited state of HC3N have a compact emission that peaks at a position consistent with the position of Sgr B2(N3). The ν7 = 1 transitions of HC3N thus trace the hot core better than the ν= 0 transitions. The nature of the extended emission traced by the ν= 0 transitions of HC3N, SO and CH3CCH is unclear. It may be related to the outflow that is seen along the same direction (Fig. 12).
![]() |
Fig. 14 Integrated intensity maps of selected transitions detected in Sgr B2(N3). The contour levels start at 3σ (rms ~ 3 mJ/beam) and increase with a step of 3σ. In each map the blue ellipse shows the result of the 2D-Gaussian fit. The black cross represents the position of the hot core derived from Fig. 1 and the red cross marks the peak position derived from the fit. The black filled ellipse represents the synthesized beam. The upper level energy of each transition is indicated in temperature units in each panel. |
4.4. Physical properties of the new hot cores
In Sects. 3.5–3.7, we presented basic physical properties of the three new hot cores discovered in Sgr B2(N). They show similarities with one of the already known hot cores, Sgr B2(N2). All of them have kinetic temperatures of ~140–180 K. Their spectra show narrow lines with typical linewidth of ~5 km s-1 similar to the values derived by Belloche et al. (2016) for Sgr B2(N2), ranging from 4.7 to 6.5 km s-1. Sgr B2(N4) and Sgr B2(N5) have mean molecular emission sizes of ~1.0″, which is slightly smaller than the size of Sgr B2(N2) (~1.2″, with values ranging between 0.8″ and 1.5″; Belloche et al. 2016, 2017). Sgr B2(N3) is more compact, with a molecular emission size of 0.4″.
In Table 10 we presented the parameters of our best fit models for ten molecules. For each hot core, we have decided to adopt the mean angular size derived from the transitions showing resolved emission. As the hot cores appear to be resolved for only a few species (especially Sgr B2(N3) for which transitions from only two species could be used to derive the source size, see Table C.1), we decided to use a single source size for each hot core to model all molecules. One must keep in mind this last consideration while comparing the chemical composition of the hot cores in Sect. 3.8 because the column densities derived from the model strongly depend on the adopted source size. For instance, we used a source size of 1.0″ to model the spectrum toward Sgr B2 (N5), although the results of the 2D-Gaussian fit to the integrated intensity maps of C2H5CN transitions suggest an average emission size of ~0.5″ for this molecule (see Table C.1). The column density of C2H5CN in Sgr B2(N5) might thus be underestimated.
In Sect. 3.7 we highlighted bipolar structures in the integrated intensity maps of the wings of typical outflow tracers. The north-south velocity gradient observed toward Sgr B2(N3) in the maps of two SO lines and the OCS(8–7) transition most probably suggests the presence of an outflow (Fig. 12). The SO(23–12) and CS(2–1) transitions investigated toward Sgr B2(N5) show a clear bipolar structure also suggesting the presence of an outflow. However, what the other transitions of SO trace toward this source is less clear. Nevertheless our interpretation that Sgr B2(N3) and Sgr B2(N5) drive outflows is reinforced by the fact that H2O maser emission is found in the close vicinity of both sources (McGrath et al. 2004,see also Fig. D.6).
Higuchi et al. (2015) have recently reported the presence of a bipolar molecular outflow in the east–west direction in Sgr B2(N1), on the basis of the SiO(2–1) and SO2(124,8–133,11) transitions in the same data set used here. They derived an average dynamical time of ~5 × 103 years, similar to our results for Sgr B2(N3) and Sgr B2(N5)(~3 × 103 yr). The fact that no outflow has been detected toward Sgr B2(N4) does not necessarily mean a real lack of outflow motion but instead it might reflect the youth of the source compared to the others. The outflow structure in this case might be too small to be detected at the resolution of our ALMA survey.
Codella et al. (2004) explored the possibility of using molecular outflows to estimate the age of the associated source. They investigated a large survey of outflows toward UCHII regions. They proposed an evolutionary scheme for regions forming high-mass stars in which class II methanol masers appear before an outflow is detectable, then both coexist in the same phase before the maser switches off as the UCHII region expands. For these reasons, the fact that no outflow has been detected toward Sgr B2(N2) might suggest that it is more evolved than the other cores because a UCHII region is already detected in the source and no class II methanol maser has been reported. However, this hypothesis is in contradiction with the conclusion derived from the comparison of the chemical composition of Sgr B2(N2) and the three new hot cores. An alternative explanation could be that there are two sources in Sgr B2(N2). Indeed, there is an offset of 0.43″ (~3600 au) between the position of peak line density and the position of the UCHII region, K7 (see Fig. 1 and Table 1). A similar situation occurs in the star-forming region W51e2 at a distance of 5.4 kpc. This source contains two objects separated by 0.8″, a hypercompact HII region, W51e2-W, and a hot molecular core, W51e2-E (Shi et al. 2010; Ginsburg 2017). At the distance of Sgr B2(N), these two objects would have a separation of 0.5″, similar to the offset seen between K7 and Sgr B2(N2). Observations at higher angular resolution with ALMA will be necessary to test this scenario.
![]() |
Fig. 15 Proposed evolutionary sequence for all the hot cores embedded in Sgr B2(N) based on their associations with UCHII regions, class II methanol masers, and outflows. K7 is the UCHII region located close to Sgr B2(N2) and may be a distinct source. |
4.5. Star formation timescales and evolutionary sequence in Sgr B2(N)
In a previous analysis of Sgr B2(N), Belloche et al. (2013) used the UCHII regions reported by Gaume et al. (1995) to estimate a star formation rate of 0.028–0.039 M⊙ yr-1 averaged over 105 yr for the whole Sgr B2 complex. We assume a constant star formation rate to estimate the lifetime of Sgr B2(N)’s hot cores. Combining the results of Gaume et al. (1995) and De Pree et al. (2015), we count ten HII regions within the ALMA primary beam centered between Sgr B2(N1) and Sgr B2(N2) (see Fig. 2). Eight of these sources are UCHII regions. Only five hot cores are detected in the same area. Considering a lifetime of ~105 yr for the UCHII regions (Peters et al. 2010), we estimate a lifetime of 10 yr for Sgr B2(N)’s hot cores. In the same way we derive a statistical lifetime of approximately 4 × 104 yr for the class II 6.7 GHz methanol masers detected in this area, which is, given the low number statistics of the Sgr B2(N) region, surprisingly consistent with the lifetime derived by van der Walt (2005) for class II methanol masers (2.5 × 104 to 4.5 × 104 yr).
Based in Fig. 5 of Codella et al. (2004) and the statistical lifetimes calculated above, we propose in Fig. 15 an evolutionary sequence of the hot cores embedded in Sgr B2(N). Among the three new hot cores, Sgr B2(N5) would be the most evolved one because it is associated with both an outflow and a class II methanol maser, and has already entered the UCHII phase. Sgr B2(N4) appears to be the youngest core because it is only associated with class II methanol maser emission. With an associated methanol maser and a detected outflow but no associated UCHII region, Sgr B2(N3) would be in-between. As mentioned in Sect. 3.8, the status of Sgr B2(N2) is unclear. Its association with a UCHII region but no methanol maser and no outflow would suggest that it is more evolved than the new hot cores, but its chemical composition seems to suggest the opposite. As discussed in Sect. 4.4, one reason for this apparent contradiction could be that Sgr B2(N2) actually contains two distinct sources, one associated with the UCHII region K7, and therefore more evolved than all the other hot cores, and the other one, the hot core, with a line density peak shifted from K7, and still too young to show outflow emission on the scales probed with ALMA or to harbor a UCHII region.
As shown in Table 2, the five hot cores detected in Sgr B2(N) have peak column densities that differ by more than one order of magnitude. Therefore, they may form stars of different final masses and our attempt to classify them with a single evolutionary sequence should be considered as only tentative. A deeper analysis of their properties, in particular at higher angular resolution, will be necessary to improve our understanding of their respective evolutionary states.
There is one complication to the evolutionary sequence proposed above. As discussed in Sects. 1.1 and 4.1 of Belloche et al. (2008) and shown here in Fig. D.6, Sgr B2(N1) coincides with the centroid position of a powerful H2O maser compact (4″ × 2″ sized) outflow (Reid et al. 1988). For their (collisional) pumping, H2O masers require temperatures of ~400 K and densities of ~109 cm-3, much higher than the values derived in this paper for Sgr B2(N1). These conditions are met in the post-shock regions of fast (J) shocks (Elitzur et al. 1989; Hollenbach et al. 2013). The compact H2O maser outflow may originate from a different source than that which drives the UCHII region K2, a situation reminiscent of the archetypical UCHII region W3(OH), which has powerful OH and methanol masers in its expanding envelope (Menten et al. 1992) and is separated by 5″ from the multiple hot core W3(OH)-H2O (Wyrowski et al. 1999). Like Sgr B2(N1), the latter drives a powerful bipolar H2O maser outflow (Hachisuka et al. 2006), but shows no methanol maser emission.
We also mention that several H2O masers are associated with a faint peak of continuum emission in our ALMA data (see Fig. D.6), which also coincides with a region of moderately enhanced spectral line density located ~3″ west of Sgr B2(N5) (Fig. 1). This region probably harbors an additional (faint) hot core that is not associated with any UCHII region or class II methanol maser.
5. Conclusions
We used the 3 mm line survey EMoCA conducted with ALMA in its cylces 0 and 1 to search for new hot cores in the Sgr B2(N) region, taking advantage of the high sensitivity of these observations. We report the discovery of three new hot cores that we call Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5), located at (αJ2000 = 17h47m19.248s, δJ2000 = −28o22′14.91″), (αJ2000 = 17h47m19.528s, δJ2000 = −28o22′32.41″), and (αJ2000 = 17h47m20.047s, δJ2000 = −28o22′41.34″) respectively. We analyzed the line survey to characterize their chemical composition and physical structure, and we searched for outflows and associations with UCHII regions or methanol masers. Our main results are summarized as follows:
- 1.
Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5) have spectral line densities above 7σ of 31, 11, and 22 lines per GHz, respectively, which qualify them as hot cores. Approximately 91%, 89%, and 93% of these lines have been identified and assigned to 22, 23, and 25 main species, respectively.
- 2.
The typical linewidth of the three new hot cores is 5 km s-1. They have rotational temperatures of ~140–180 K and the maps of emission lines indicate mean source sizes of 0.4″ for Sgr B2(N3) and 1.0″ for Sgr B2(N4) and Sgr B2(N5).
- 3.
Assuming a dust temperature of 150 K similar to the rotational temperatures obtained from the emission lines, we derive a dust emissivity index β ~ 1.2 by combining ALMA and SMA data. We obtain H2 column densities of <5 × 1023, 3 × 1024, and 1 × 1024 cm-2 for Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5), respectively, for the sizes listed above.
- 4.
We report the detection of outflows in Sgr B2(N3) and Sgr B2(N5) based on the analysis of three SO lines, and the OCS(8–7) and CS(2–1) transitions. The outflow dynamical time is ~2.5×103 yr for both sources. No outflow is detected toward Sgr B2(N4) and Sgr B2(N2).
- 5.
Each new hot core is associated with a 6.7 GHz class II methanol maser. Sgr B2(N4) is also associated with a UCHII region.
- 6.
We derived the column densities and abundances of ten (complex) organic molecules toward Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5). The three sources share a similar chemical composition, with Sgr B2(N3) resembling Sgr B2(N5) slightly more than Sgr B2(N4). However, they differ from Sgr B2(N2) and several molecular ratios suggest that Sgr B2(N2) is chemically less evolved than the three new hot cores.
- 7.
Assuming a lifetime of 105 yr for UCHII regions, we derive statistical lifetimes of 4 × 104 yr for the class II methanol maser phase and 6 × 104 yr for the hot core phase in Sgr B2(N).
Given that their peak column densities differ by more than one order of magnitude, the five hot cores may form stars of different final masses. Still, their associations with class II methanol masers, outflows, and/or UCHII regions tentatively suggest the following age sequence from the youngest to the oldest source: Sgr B2(N4), Sgr B2(N3), Sgr B2(N5), Sgr B2(N1). The status of Sgr B2(N2) is puzzling. Its association with a UCHII region but no outflow and no methanol maser suggests that it should be the oldest source in this sequence. However, this contradicts its youth suggested by its chemical composition. An explanation may be that Sgr B2(N2) contains two distinct sources, as suggested by the small angular offset separating its embedded UCHII region from the line density peak of the hot core. On-going observations at higher angular resolution with ALMA will help understanding the status of this source.
See, e.g., http://www.astro.uni-koeln.de/cdms/molecules
Acknowledgments
We thank Sheng-Li Qin and Peter Schilke for providing us the SMA map of Sgr B2(N) in electronic form. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00017.S, ADS/JAO.ALMA#2012.1.00012.S. ALMA is a partnership of ESO (representing its member states), NSF (USA), and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO, and NAOJ. The interferometric data are available in the ALMA archive at https://almascience.eso.org/aq/. This work has been in part supported by the Deutsche Forschungsgemeinschaft (DFG) through the collaborative research grant SFB 956 “Conditions and Impact of Star Formation”, project area B3.
References
- Belloche, A., Menten, K. M., Comito, C., et al. 2008, A&A, 482, 179 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Belloche, A., Müller, H. S. P., Menten, K. M., Schilke, P., & Comito, C. 2013, A&A, 559, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Belloche, A., Garrod, R. T., Müller, H. S. P., & Menten, K. M. 2014, Science, 345, 1584 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Belloche, A., Müller, H. S. P., Garrod, R. T., & Menten, K. M. 2016, A&A, 587, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Belloche, A., Meshcheryakov, A. A., Garrod, R. T., et al. 2017, A&A, 601, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Caswell, J. L. 1996, MNRAS, 283, 606 [NASA ADS] [CrossRef] [Google Scholar]
- Codella, C., Lorenzani, A., Gallego, A. T., Cesaroni, R., & Moscadelli, L. 2004, A&A, 417, 615 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- De Pree, C. G., Goss, W. M., & Gaume, R. A. 1998, ApJ, 500, 847 [NASA ADS] [CrossRef] [Google Scholar]
- De Pree, C. G., Peters, T., Mac Low, M. M., et al. 2015, ApJ, 815, 123 [NASA ADS] [CrossRef] [Google Scholar]
- Elitzur, M., Hollenbach, D. J., & McKee, C. F. 1989, ApJ, 346, 983 [NASA ADS] [CrossRef] [Google Scholar]
- Endres, C. P., Schlemmer, S., Schilke, P., Stutzki, J., & Müller, H. S. P. 2016, J. Mol. Spectr., 327, 95 [NASA ADS] [CrossRef] [Google Scholar]
- Garrod, R. T., Widicus Weaver, S. L., & Herbst, E. 2008, ApJ, 682, 283 [NASA ADS] [CrossRef] [Google Scholar]
- Garrod, R. T., Vasyunin, A. I., Semenov, D. A., Wiebe, D. S., & Henning, T. 2009, ApJ, 700, L43 [Google Scholar]
- Gaume, R. A., Claussen, M. J., de Pree, C. G., Goss, W. M., & Mehringer, D. M. 1995, ApJ, 449, 663 [NASA ADS] [CrossRef] [Google Scholar]
- Ginsburg, A. 2017, ArXiv e-prints [arXiv:1702.06627] [Google Scholar]
- Goldsmith, P. F., & Langer, W. D. 1999, ApJ, 517, 209 [NASA ADS] [CrossRef] [Google Scholar]
- Hachisuka, K., Brunthaler, A., Menten, K. M., et al. 2006, ApJ, 645, 337 [NASA ADS] [CrossRef] [Google Scholar]
- Herbst, E., & van Dishoeck, E. F. 2009, ARA&A, 47, 427 [NASA ADS] [CrossRef] [Google Scholar]
- Higuchi, A. E., Hasegawa, T., Saigo, K., Sanhueza, P., & Chibueze, J. O. 2015, ApJ, 815, 106 [NASA ADS] [CrossRef] [Google Scholar]
- Hollenbach, D., Elitzur, M., & McKee, C. F. 2013, ApJ, 773, 70 [NASA ADS] [CrossRef] [Google Scholar]
- Hollis, J. M., Lovas, F. J., & Jewell, P. R. 2000, ApJ, 540, L107 [NASA ADS] [CrossRef] [Google Scholar]
- Hollis, J. M., Lovas, F. J., Remijan, A. J., et al. 2006, ApJ, 643, L25 [Google Scholar]
- Hu, B., Menten, K. M., Wu, Y., et al. 2016, ApJ, 833, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Kauffmann, J., Bertoldi, F., Bourke, T. L., Evans, I. N. J., & Lee, C. W. 2008, A&A, 487, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Keto, E. 2007, ApJ, 666, 976 [NASA ADS] [CrossRef] [Google Scholar]
- Kuan, Y.-J., & Snyder, L. E. 1996, ApJ, 470, 981 [NASA ADS] [CrossRef] [Google Scholar]
- Lis, D. C., & Goldsmith, P. F. 1990, ApJ, 356, 195 [NASA ADS] [CrossRef] [Google Scholar]
- Maret, S., Hily-Blant, P., Pety, J., Bardeau, S., & Reynier, E. 2011, A&A, 526, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- McGrath, E. J., Goss, W. M., & De Pree, C. G. 2004, ApJS, 155, 577 [NASA ADS] [CrossRef] [Google Scholar]
- Mehringer, D. M., Palmer, P., Goss, W. M., & Yusef-Zadeh, F. 1993, ApJ, 412, 684 [NASA ADS] [CrossRef] [Google Scholar]
- Mehringer, D. M., Snyder, L. E., Miao, Y., & Lovas, F. J. 1997, ApJ, 480, L71 [NASA ADS] [CrossRef] [Google Scholar]
- Menten, K. M. 1993, 412, 199 [Google Scholar]
- Menten, K. M. 2004, in The Dense Interstellar Medium in Galaxies, eds. S. Pfalzner, C. Kramer, C. Staubmeier, & A. Heithausen, Springer Proc. Phys., 91, 69 [Google Scholar]
- Menten, K. M., Reid, M. J., Pratap, P., Moran, J. M., & Wilson, T. L. 1992, ApJ, 401, L39 [NASA ADS] [CrossRef] [Google Scholar]
- Minier, V., Ellingsen, S. P., Norris, R. P., & Booth, R. S. 2003, A&A, 403, 1095 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Müller, H. S. P., Belloche, A., Xu, L.-H., et al. 2016, A&A, 587, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943 [NASA ADS] [Google Scholar]
- Pandian, J. D., Leurini, S., Menten, K. M., Belloche, A., & Goldsmith, P. F. 2008, A&A, 489, 1175 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pearson, J. C., Müller, H. S. P., Pickett, H. M., Cohen, E. A., & Drouin, B. J. 2010, J. Quant. Spec. Radiat. Transf., 111, 1614 [Google Scholar]
- Peters, T., Mac Low, M.-M., Banerjee, R., Klessen, R. S., & Dullemond, C. P. 2010, ApJ, 719, 831 [NASA ADS] [CrossRef] [Google Scholar]
- Pillai, T., Wyrowski, F., Menten, K. M., & Krügel, E. 2006, A&A, 447, 929 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Qin, S.-L., Schilke, P., Rolffs, R., et al. 2011, A&A, 530, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Reid, M. J., Schneps, M. H., Moran, J. M., et al. 1988, ApJ, 330, 809 [NASA ADS] [CrossRef] [Google Scholar]
- Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2014, ApJ, 783, 130 [Google Scholar]
- Remijan, A. J. 2015, IAU General Assembly, 22, 2257092 [NASA ADS] [Google Scholar]
- Rosero, V., Hofner, P., Claussen, M., et al. 2016, ApJS, 227, 25 [NASA ADS] [CrossRef] [Google Scholar]
- Shi, H., Zhao, J.-H., & Han, J. L. 2010, ApJ, 710, 843 [NASA ADS] [CrossRef] [Google Scholar]
- Snyder, L. E., Lovas, F. J., Hollis, J. M., et al. 2005, ApJ, 619, 914 [NASA ADS] [CrossRef] [Google Scholar]
- Sobolev, A. M., Cragg, D. M., & Godfrey, P. D. 1997, A&A, 324, 211 [NASA ADS] [Google Scholar]
- Urquhart, J. S., Moore, T. J. T., Csengeri, T., et al. 2014, MNRAS, 443, 1555 [NASA ADS] [CrossRef] [Google Scholar]
- van der Walt, J. 2005, MNRAS, 360, 153 [NASA ADS] [CrossRef] [Google Scholar]
- Walsh, A. J., Hyland, A. R., Robinson, G., & Burton, M. G. 1997, MNRAS, 291, 261 [NASA ADS] [CrossRef] [Google Scholar]
- Walsh, A. J., Burton, M. G., Hyland, A. R., & Robinson, G. 1998, MNRAS, 301, 640 [NASA ADS] [CrossRef] [Google Scholar]
- Wyrowski, F., Schilke, P., Walmsley, C. M., & Menten, K. M. 1999, ApJ, 514, L43 [NASA ADS] [CrossRef] [Google Scholar]
- Xu, Y., Li, J. J., Hachisuka, K., et al. 2008, A&A, 485, 729 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Estimating the level of free-free emission
Here we present the procedure used to estimate the contribution of free-free emission expected toward Sgr B2(N1), Sgr B2(N2), and Sgr B2(N5) in the ALMA continuum maps. As mentioned in Sect. 3.2.1, the overall shape of the extended continuum emission detected with ALMA is similar to the shape of the 1.3 cm free-free emission reported by Gaume et al. (1995) that suggests that the extended emission detected with ALMA is dominated by free-free emission. Figure A.1b shows the profile of the continuum emission along the direction going through Sgr B2(N2) indicated in red in Fig. A.1a. Because Sgr B2(N2) is located in the direction of extended 1.3 cm free-free emission, we decompose the profile shown in Fig. A.1b into a peak, which we attribute to dust emission associated with Sgr B2(N2), and a pedestal, which we attribute to the extended free-free emission. We proceed in the same way for the 20 ALMA continuum maps. The results are plotted in Fig. A.2. The weighted linear fit to the results excluding setup 3 gives us the correction to subtract from the measured peak density as a function of frequency. For the maps belonging to the frequency range of setup 3 (with the lowest angular resolution), we directly apply the corrections estimated from the plots of the continuum profile. The corrections applied here take into account the contamination from the extended continuum emission detected in the ALMA continuum maps but a correction for the free-free emission arising from the UCHII region K7 that is associated with Sgr B2(N2) still has
![]() |
Fig. A.1 a) Continuum map obtained at 86.5 GHz with ALMA. The contour levels are 20σ, 40σ, 60σ, 100σ, 140σ, and 160σ, with σ = 4.3 mJy. The black crosses mark the peak positions of the continuum emission observed toward Sgr B2(N1) and Sgr B2(N2) derived from Fig. 2. The white cross represents the phase center. The filled ellipse shows the clean beam (2.04″× 1.51″, PA = −83.3o). b) Profile of the continuum emission along the direction plotted in red in panel a. The blue line marks the estimated level of free-free emission. |
to be applied. We estimate the flux expected from K7 based on the flux reported by De Pree et al. (2015). They measured an integrated flux of 30 mJy at 44.2 GHz for a source size of 0.08″. We assume that the free-free emission of K7 is optically thin and thus that the measured flux is proportional to να, with α = −0.1. From this relation we extrapolate a flux of ~28 mJy for the frequency range of the ALMA survey. This value is added to the correction for Sgr B2(N2) derived from Fig. A.2.
We proceed in the same way for Sgr B2(N5), based on the flux reported by Gaume et al. (1995) toward Z10.24. They measured an integrated flux of 36 mJy at 22.4 GHz for a source size <0.25″. We extrapolate this value to the ALMA frequency range assuming optically thin emission, and we obtain ~31 mJy. This value is directly subtracted from the values of the peak flux density measured in the ALMA continuum maps toward Sgr B2(N5).
Finally in the case of Sgr B2(N1) we use the flux measured by De Pree et al. (2015) toward K2 and K3 at 44.2 GHz. They measured integrated flux densities of 80 mJy toward K2 for a source size of 0.12″ and 242 mJy toward K3 for a source size of 0.27″. K3 being located 1.3″ north-east of K2, we estimate that only~15% of the free-free emission arising from this UCHII region contributes to the total flux density measured in the ALMA synthesized beam of ~1.6″ centered on K2. We thus obtain a total free-free contribution of ~116 mJy that we subtract from the flux densities measured toward Sgr B2(N1) in the ALMA continuum maps.
![]() |
Fig. A.2 Fraction of free-free emission relative to the total flux density measured toward Sgr B2(N2). This does not include the contribution of K7. The blue line is the weighted linear fit to the results excluding setup 3 that is shown in red. |
Appendix B: Additional tables: continuum analysis
Peak flux densities and H2 column densities toward Sgr B2(N1).
Peak flux densities and H2 column densities toward Sgr B2(N2).
H2 column density upper limits toward Sgr B2(N3).
Peak flux densities and H2 column densities toward Sgr B2(N4).
Peak flux densities and H2 column densities toward Sgr B2(N5).
Appendix C: Additional table: source size
Results of elliptical 2D-Gaussian fits to the integrated intensity maps of the transitions with resolved emission.
Appendix D: Additional figures
![]() |
Fig. D.1 Results of the 2D-Gaussian fits to the integrated intensity maps showing resolved emission toward Sgr B2(N3). a), b) Fitted position. The error bars are the uncertainties given by the Gaussian fits. c), d) Fitted size. The dots represent the deconvolved major and minor diameters of the emission (c) and the resulting average source size (d). The squares represent the major and minor axes of the synthesized beam (c). The dashed lines represent the reference position of the hot core (a, b), and the mean deconvolved angular size (d). |
![]() |
Fig. D.4 Lines of typical outflow tracers detected toward Sgr B2(N4). The red spectrum represents our LTE model fit to the complete observed spectrum of the molecule. The green spectrum shows the model that contains all the molecules identified so far toward Sgr B2(N4). The dashed horizontal line shows the 3σ level. The systemic velocity of the source is marked with the dashed vertical line. The rest frequency and upper level energy (in temperature units) of each transition are indicated in each panel. |
![]() |
Fig. D.6 Continuum map of the Sgr B2(N) region obtained with ALMA at 108 GHz. Contour levels (positive in black solid line and negative in dashed line) start at five times the rms noise level, σ, of 3.0 mJy/beam and double in value up to 320σ. The filled ellipse shows the synthesized beam (1.65″× 1.21″, PA = −83.4o). The black cross represents the phase center. The red crosses mark the positions of the five hot cores embedded in Sgr B2(N). The blue crosses represent H2O masers (McGrath et al. 2004). The dotted red circle represents the size (HPBW) of the primary beam of the 12 m antennas at 108 GHz. |
All Tables
Position and spectral line density of the hot cores embedded in Sgr B2(N), and distance to UCHII regions and class II methanol masers.
Position and size of the ALMA continuum sources, and distances to hot cores and UCHII regions.
Results of elliptical 2D-Gaussian fits to the integrated intensity maps of the transitions with resolved emission.
All Figures
![]() |
Fig. 1 Contour map of the number of channels with continuum-subtracted flux density above the 7σ level (1σ ~ 3 mJy/beam). The contour levels are: 500, 1000, 2000, 5000, 10 000, 20 000, 30 000, and 40 000. The blue crosses mark the peaks of spectral line density. Red crosses represent the 6.7 GHz class II methanol masers (Caswell 1996). Green crosses represent the compact and ultra-compact HII regions (Gaume et al. 1995; De Pree et al. 2015). The offsets are defined with respect to the phase center (see Sect. 2.1). |
In the text |
![]() |
Fig. 2 Continuum map of the Sgr B2(N) region obtained with ALMA at 108 GHz. Contour levels (positive in black solid line and negative in dashed line) start at five times the rms noise level, σ, of 3.0 mJy/beam and double in value up to 320σ. The filled ellipse shows the synthesized beam (1.65″× 1.21″, PA = −83.4o). The black cross represents the phase center. The blue crosses mark the positions of the five hot cores embedded in Sgr B2(N) derived from Fig. 1. The red crosses represent the average peak positions of the continuum emission derived from 2D-Gaussian fits to all ALMA continuum maps. The green crosses mark the compact and ultra-compact HII regions. The dotted red circle represents the size (HPBW) of the primary beam of the 12 m antennas at 108 GHz. The map is not corrected for the primary beam attenuation. |
In the text |
![]() |
Fig. 3 Continuum map of the Sgr B2(N) region obtained with the SMA at 343 GHz (Qin et al. 2011) and smoothed to the ALMA resolution. The new beam is shown in the bottom left corner (1.65″× 1.21″, PA = −83.4o). The contour levels (positive in solid line, negative in dashed line) start at 8σ and double up to 128σ with σ = 0.15 Jy/beam, the noise level measured inside the polygon defined around Sgr B2(N3). The black cross represents the SMA phase center. The red crosses mark the peak positions of the ALMA continuum emission derived from Fig. 2. The position of Sgr B2(N3) (blue cross) is derived from Fig. 1. The red dashed circle represents the size (HPBW) of the primary beam of the SMA 6 m antennas at 343 GHz. The map is corrected for the primary beam attenuation. |
In the text |
![]() |
Fig. 4 H2 column densities as a function of frequency for the five hot cores embedded in Sgr B2(N). The values obtained based on the ALMA data before and after correction for the free-free contribution are shown as black and red crosses respectively. The green crosses represent the results obtained based on the SMA map at 343 GHz smoothed to the angular resolution of the ALMA maps. Error bars are calculated from the error on |
In the text |
![]() |
Fig. 5 Part of the continuum-subtracted ALMA spectra observed toward the hot cores embedded in Sgr B2(N). The spectra have been corrected for the primary beam attenuation, and shifted along the y axis for display purposes. The spectra of Sgr B2(N1) and Sgr B2(N2) have been divided by 10 and 6, respectively. The frequency axis corresponds to the systemic velocities derived in Sect. 3.7. |
In the text |
![]() |
Fig. 6 Spatial distribution of the molecules identified around Sgr B2(N3), Sgr B2(N4), and Sgr B2(N5). In each panel, the red star represents the position of the hot core derived from Fig. 1. Blue crosses represent the mean peak position for each species. Error bars correspond to the standard deviation weighted by the uncertainties given by the GAUSS-2D routine. The radius of the dashed circle is given in the top right corner. The average peak position of HC3N, ν7 = 1 is shown in green. |
In the text |
![]() |
Fig. 7 Integrated intensity maps of selected transitions toward Sgr B2(N3) (top row), Sgr B2(N4) (middle row), and Sgr B2(N5) (bottom row). For each map, the red cross shows the peak position of the emission and the black cross is the reference position of the hot core derived from Fig. 1. The blue ellipse represents the result of the Gaussian fit to the map while the red ellipse is the deconvolved emission size. The black filled ellipse represents the synthesized beam. The rms and contour levels are indicated in Table 6. |
In the text |
![]() |
Fig. 8 Population diagrams of C2H5OH (panels a and b) and CH3OCHO (panels c and d) for Sgr B2(N3). The black points are computed using the integrated intensities of the observed spectrum while the red points are computed using the integrated intensities of our synthetic model. The error bars on the observed data points are 1σ uncertainties on Nu/gu. No correction is applied in panels a and c, while in panels b and d the optical depth correction has been applied to both the observed and synthetic populations and the contamination from all other species included in the full LTE model has been removed from the observed data points. The blue line is the weighted linear fit to the observed populations. |
In the text |
![]() |
Fig. 9 As in Fig. 8 but for Sgr B2(N4). |
In the text |
![]() |
Fig. 10 As in Fig. 8 but for Sgr B2(N5). |
In the text |
![]() |
Fig. 11 Spectra of the lines investigated to search for outflows toward Sgr B2(N3) (top row) and Sgr B2(N5) (bottom row). The dashed vertical line marks the systemic velocity of the source and the high velocity wings are highlighted in blue and red. The magenta spectrum represents our best-fit model while the green spectrum shows the model that contains all the identified species. The rest frequency and upper level energy (in temperature units) of each transition are indicated in each panel. |
In the text |
![]() |
Fig. 12 Integrated intensity maps of the lines shown in Fig. 11 toward Sgr B2(N3) (top row) and Sgr B2(N5) (bottom row). In each panel, maps of the blue- and red-shifted wings are presented in blue and red contours, respectively, overlaid on the integrated emission of the line core (black contours). For the line cores affected by absorption, the black contours represent the continuum emission. Rms noise levels and contour levels used for each map are listed in Table 9. Each cross corresponds to the peak position of the emission. The black triangle marks the position of the hot core (Sgr B2(N3) or Sgr B2(N5)) derived from Fig. 1. |
In the text |
![]() |
Fig. 13 a) Column densities of various molecules detected toward four of the five hot cores embedded in Sgr B2(N) (see Table 10). b) Column densities normalized to the column density of C2H5CN. c) Chemical abundances with respect to H2. In each panel, lower and upper limits are indicated with arrows. |
In the text |
![]() |
Fig. 14 Integrated intensity maps of selected transitions detected in Sgr B2(N3). The contour levels start at 3σ (rms ~ 3 mJ/beam) and increase with a step of 3σ. In each map the blue ellipse shows the result of the 2D-Gaussian fit. The black cross represents the position of the hot core derived from Fig. 1 and the red cross marks the peak position derived from the fit. The black filled ellipse represents the synthesized beam. The upper level energy of each transition is indicated in temperature units in each panel. |
In the text |
![]() |
Fig. 15 Proposed evolutionary sequence for all the hot cores embedded in Sgr B2(N) based on their associations with UCHII regions, class II methanol masers, and outflows. K7 is the UCHII region located close to Sgr B2(N2) and may be a distinct source. |
In the text |
![]() |
Fig. A.1 a) Continuum map obtained at 86.5 GHz with ALMA. The contour levels are 20σ, 40σ, 60σ, 100σ, 140σ, and 160σ, with σ = 4.3 mJy. The black crosses mark the peak positions of the continuum emission observed toward Sgr B2(N1) and Sgr B2(N2) derived from Fig. 2. The white cross represents the phase center. The filled ellipse shows the clean beam (2.04″× 1.51″, PA = −83.3o). b) Profile of the continuum emission along the direction plotted in red in panel a. The blue line marks the estimated level of free-free emission. |
In the text |
![]() |
Fig. A.2 Fraction of free-free emission relative to the total flux density measured toward Sgr B2(N2). This does not include the contribution of K7. The blue line is the weighted linear fit to the results excluding setup 3 that is shown in red. |
In the text |
![]() |
Fig. D.1 Results of the 2D-Gaussian fits to the integrated intensity maps showing resolved emission toward Sgr B2(N3). a), b) Fitted position. The error bars are the uncertainties given by the Gaussian fits. c), d) Fitted size. The dots represent the deconvolved major and minor diameters of the emission (c) and the resulting average source size (d). The squares represent the major and minor axes of the synthesized beam (c). The dashed lines represent the reference position of the hot core (a, b), and the mean deconvolved angular size (d). |
In the text |
![]() |
Fig. D.2 As in Fig. D.1 but for Sgr B2(N4). |
In the text |
![]() |
Fig. D.3 As in Fig. D.1 but for Sgr B2(N5). |
In the text |
![]() |
Fig. D.4 Lines of typical outflow tracers detected toward Sgr B2(N4). The red spectrum represents our LTE model fit to the complete observed spectrum of the molecule. The green spectrum shows the model that contains all the molecules identified so far toward Sgr B2(N4). The dashed horizontal line shows the 3σ level. The systemic velocity of the source is marked with the dashed vertical line. The rest frequency and upper level energy (in temperature units) of each transition are indicated in each panel. |
In the text |
![]() |
Fig. D.5 As in Fig. D.4 but for species not considered as typical outflow tracers. |
In the text |
![]() |
Fig. D.6 Continuum map of the Sgr B2(N) region obtained with ALMA at 108 GHz. Contour levels (positive in black solid line and negative in dashed line) start at five times the rms noise level, σ, of 3.0 mJy/beam and double in value up to 320σ. The filled ellipse shows the synthesized beam (1.65″× 1.21″, PA = −83.4o). The black cross represents the phase center. The red crosses mark the positions of the five hot cores embedded in Sgr B2(N). The blue crosses represent H2O masers (McGrath et al. 2004). The dotted red circle represents the size (HPBW) of the primary beam of the 12 m antennas at 108 GHz. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.