Free Access
Issue
A&A
Volume 648, April 2021
Article Number A100
Number of page(s) 15
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202040121
Published online 20 April 2021

© ESO 2021

1 Introduction

Despite the many theoretical and observational efforts aimed at understanding the formation process of early-type stars, only limited progress has been made. While there is growing evidence that B-type or late O-type stars, with masses up to ~20 M, may form through disk accretion (Beltrán & de Wit 2016), the number of early O-type (proto)stars (L≳105 L) with Keplerian rotation signatures is still very limited (see Beltrán (2020), and references therein). Therefore, it is not yet clear whether disk accretion is the only or main formation mechanism. These objects are rare, distant, and lie in complex regions inside rich stellar clusters. All these features make it difficult from an observational point of view to characterize the environment around the star and thus to compare theoretical predictions (e.g., Bonnell & Bate 2006; Krumholz et al. 2009; Tan et al. 2014; Motte et al. 2018; Kuiper & Hosokawa 2018; Kumar et al. 2020). It is also possible that the complexity of such an environment (Peters et al. 2010) makes any model proposed to date far too simple to be satisfactory.

In an attempt to shed light on this issue and, in particular, establish the role of circumstellar disks in the formation of O-type stars, Cesaroni et al. (2017) performed a mini-survey at 1 mm with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of ~0.′′2 of six luminous sources (≳5 × 104 L) associated with known hot molecular cores (HMCs). These observations have revealed the presence of velocity gradients consistent with rotation in almost all cores and different degrees of fragmentation, with two cores presenting no evidence of fragmentation.

One of the HMCs that showed no signs of fragmentation was G31.41+0.31 (hereafter G31). G31 is a well-known HMC located at 3.75 kpc (estimated from trigonometric parallax; Immer et al. 2019), with a luminosity of ~ 4.5 × 104 L (from Osorio et al. 2009, after scaling to the parallactic distance). The core is ~ 5″ away from an ultracompact (UC) HII region and harbors two point-like (< 0.′′07) free-free continuum sources close to the center separated by ~ 0.′′2 (Cesaroni et al. 2010, hereafter C10). Given the intensity of the molecular line and continuum emission, interferometric studies at millimeter wavelengths with subarcsecond resolution (~ 0.′′8; Cesaroni et al. 2011) were feasible even before the advent of ALMA. Such observations have confirmed that the core is chemically very rich, presenting prominent emission in a large number of complex organic molecules (Beltrán et al. 2005, 2009; Fontani et al. 2007; Calcutt et al. 2014; Rivilla et al. 2017; Mininni et al. 2020). They have been successfully used to trace the distribution and kinematics of the dense hot gas, and to estimate basic physical parameters (e.g., temperature and density; Beltrán et al. 2005). The presence of a clear velocity gradient along the NE–SW direction in the core has been established (Beltrán et al. 2004; Girart et al. 2009; Cesaroni et al. 2011). This gradient has been interpreted as rotational motion around embedded massive stars, whose existence is suggested by the presence of the two free-free continuum sources (C10). Polarization measurements by Girart et al. (2009) and Beltrán et al. (2019) have revealed a magnetic field shaped like an hourglass, with the symmetry axis oriented perpendicular to the velocity gradient. Moreover, inverse P Cygni profiles in a few molecular lines have been detected toward the core (Girart et al. 2009; Mayen-Gijon et al. 2014), indicating that infall motions are present. All these features are consistent with a scenario where the core is collapsing and rotating about the direction of the magnetic field axis.

In a recent study, Beltrán et al. (2018, hereafter B18) analyzed in detail the line and continuum ALMA observations at 1 mm with ~ 0.′′2 (~ 750 au) angular resolution initially presented in Cesaroni et al. (2017). This high angular resolution allowed us for the first time to clearly resolve the core, whose angular diameter is ~ 2′′. The most remarkable results of this study are the following.

The dust continuum emission has been resolved into two cores, a Main core that peaks close to the position of the two free-free continuum sources detected by C10 and a much weaker and smaller NE core.

The Main core imaged in the continuum is well resolved and looks monolithic and featureless, with no hint of fragmentation, despite the high mass (~ 70 M; Cesaroni 2019) and the diameter of ~8000 au.

The velocity gradient along the NE–SW direction observed in the Main core is due to a smooth change in velocity across the core, which appears to increase linearly with distance from the core center, consistent with solid-body rotation, as shown in the position–velocity (PV) plot in Fig. 8 of B18.

The rotation seems to speed up with increasing energy of the CH3CN transition, with the PV plot in the vibrationally excited CH3CN(12–11) v8 = 1 K, l = 6, 1 (Eup = 778 K) line being steeper than the same plot for the ground-state CH313$_{3}^{13}$CN(12–11) K = 2 (Eup = 97 K) line (see Fig. 9 of B18). This spin-up toward the center of the rotation is expected from conservation of angular momentum in a rotating and infalling structure.

The detection of red-shifted absorption toward the core center supports the presence of infall in the core. Interestingly, the velocity of the absorption feature increases with the line excitation energy, which suggests that infall is accelerating toward the core center (where the temperature is higher), consistent with gravitational collapse on a central source.

These results taken together pose serious questions about the nature of the core. On the one hand, the presence of a magnetic field uniformly threading the core (Girart et al. 2009; Beltrán et al. 2019), the linear increase in rotation velocity with radius, and the (apparent) homogeneity of the core, point to a monolithic core stabilized by the magnetic field and undergoing solid-body rotation, similar to the scenario proposed by the “monolithic collapse” theory (also known as turbulent core model) of McKee & Tan (2003). On the other hand, the presence of red-shifted absorption, the existence of two embedded massive stars (the two free-free sources), and the rotational speed-up toward the core center are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum, a scenario reminiscent of the “competitive accretion” model of Bonnell & Bate (2006). However, the striking monolithic appearance of the Main core may be the result of a combination of high opacity of the 220 GHz continuum emission and insufficient angular resolution. We can reasonably assume that the size of a typical fragment is of the same order as the separation between the two free-free sources detected inside the core, ~ 0.′′2 or ~ 750 au, comparable to the synthesized beam of our ALMA maps, which could have prevented B18 from properly imaging the fragmentation in the core, if any was present. At the same time we must take into account that the continuum emission observed at 220 GHz is optically thick, as proved by the high brightness temperature of the continuum peak (~ 130 K) compared to the dust temperature (≳150 K: Beltrán et al. 2005) and the presence of deep absorption against this bright continuum, seen in many lines (Rivilla et al. 2017; B18). The high opacity is likely to mask any inhomogeneity of the core.

To discriminate between the two scenarios depicted above for the G31 core, either a (marginally) stable monolithic core undergoing solid-body rotation or a dynamically collapsing core undergoing global fragmentation, we observed the G31 core with ALMA at higher angular resolution (~ 0.′′1, or 375 au) at 1.4 mm, which should allow us to resolve the single fragments and estimate their masses. In addition, we also observed the core at a slightly higher angular resolution (~ 0.′′075, or ~ 280 au) at 3.5 mm, where the dust continuum opacity may be up to about four times lower than at 1.4 mm. These observations should allow us to discard or confirm that the core is monolithic as it appears to be in the previous 1.4 mm observations. We also carried out Very Large Array (VLA) observations at 7 mm and 1.3 cm at an even higher angular resolution (0.′′ 05–0.′′07) to better study the nature of the free-free continuum sources detected by C10 and possibly detect other fainter sources embedded in the core.

This article is organized as follows. In Sect. 2, we describe the ALMA and VLA observations. In Sect. 3, we present the results and analyze the continuum emission at millimeter and centimeter wavelengths toward G31. In Sect. 4, we discuss the most likely scenario for G31, and the role of the magnetic field. Finally, in Sect. 5 we give our main conclusions.

2 Observations

2.1 ALMA

Interferometric observations of G31 at 1.4 mm in Band 6 and 3.5 mm in Band 3 were carried out with ALMA in Cycle 4 and Cycle 6 as part of projects 2016.1.00223.S (P.I.: M. Beltrán) and 2018.1.00252.S (P.I.: M. Beltrán). Observations at 1.4 mm were carried out in one of the more extended configurations, C40-7. The 1.4 mm observations were combined with those of the ALMA Cycle 2 project 2013.1.00489.S (P.I.: R. Cesaroni), which were observed in the extended C34-7/(6) configuration and with the same digital correlator configuration. We refer to Cesaroni et al. (2017) and B18 for detailed information on the Cycle 2 observations. Observations at 3.5 mm were performed in the very extended C40-9 and C43-9 configurations, and in the more compact C40-6 and C43-6 configurations (see Table 1 for more information). Observations at both wavelengths were initially carried out in Cycle 4 (project 2016.1.00223.S), between September 2016 and September 2017. After the project was completed, post-processing calibration and analysis of the 3.5 mm data revealed severe problems on the bandpass shape of all spectral windows for both extended and compact configurations. This raised serious doubts about the calibration and overall quality of the continuum map, in particular of the fluxes of the sources. Reprocessing of the data by ALMA resulted in a quality assessment failure (QA2_FAIL) and the scheduling blocks were put back in the observing queue of Cycle 5, but unfortunately were not scheduled so they could not be re-observed. After proposing again in Cycle 6, the new observations at 3.5 mm in both extended and compact configurations were finally carried out as part of the new project 2018.1.00252.S. In early 2020, after a thorough recalibration of the 3.5 mm Cycle 4 data1 by one of the co-authors of this work, L. Maud, the spectral baselines were finally properly corrected, and these data have been used in combination with those of Cycle 6 to obtain the continuum emission map at 3.5 mm.

The ALMA 1.4 mm C40-7 observations used baselines in the range 21–3637 m, providing sensitivity to structures < 1.′′4. The Cycle 2 C34-7/(6) observations used baselines in the range 40–1500 m. The digital correlator was configured in 13 spectral windows (SPWs), with SPW0 having a bandwidth of 1875 MHz and 3840 channels, while the rest used a 234.38 MHz bandwidth and 960 channels. This translates into a velocity resolution of ~ 1.3 km s−1 for SPW0 and ~ 0.66–0.77 km s−1 for the other SPWs. At 3.5 mm, the source was observed in compact and extended configurations. Because the data were observed twice in different cycles, the baselines and maximum recoverable scales vary even for the same configurations, so we report the different values for i) the extended configuration: baselines 41–12145 m (C40-9) and 83–16196 m (C43-9), and maximum recoverable scales < 1.′′75 (C40-9) and < 0.′′97 (C43-9); and ii) the compact configuration: baselines 15–3247 m (C40-6) and 41–3637 m (C43-6), and maximum recoverable scales < 5.′′7 (C40-6) and < 4.′′4 (C43-6). The digital correlator at 3.5 mm was configured in four SPWs of 1875 MHz and 3840 channels each, which translates into a velocity resolution of~3.4−4.0 km s−1.

The phase reference center of the observations was set to the position α(J 2000) = 18h 47m 34. s 315, δ(J 2000) = −01° 12′ 45.′′ 9. Flux and bandpass calibrations were achieved through observations of J1751+0939 and J1924−2914, while gain calibration through observations of J1851+0035. The data were calibrated and imaged using the CASA2 software package. Following ALMA Memo 5993, we conservatively assumed that the uncertainties on the absolute flux density calibration were ~ 6% at 3.5 mm and ~10% at 1.4 mm. Maps were created with the ROBUST parameter of Briggs (1995) set equal to 0. Further imaging and analysis were done with the GILDAS4 software package. The continuum was determined from SPW0, which is the broadest spectral window at 1.4 mm, centered at 218 GHz, and the first spectral window at 3.5 mm, centered at 85.3 GHz. The continuum was subtracted from the line emission using the STATCONT5 algorithm (Sánchez-Monge et al. 2018). The resulting synthesized CLEANed beam and rms noise of the continuum maps are given in Table 1.

Table 1

Summary of the continuum observations.

2.2 VLA

Interferometric continuum observations of G31 at 1.3 cm (K band) and 7 mm (Q band) were carried out with the VLA of the National Radio Astronomy Observatory (NRAO) as part of the project 18A−031 (PI: V.M. Rivilla). The observations were done in the extended A configuration between March and April 2018. The VLA K-band and Q-band observationsused baselines in the range 680–36 400 m. We used the WIDAR backend with 2 GHz baseband pairs using the 3-bit samplers forboth bands. The four basebands with 2.048 GHz of bandwidth were centered at the frequencies 19.2, 21.2, 23.2, 25.2 GHz (K band), and 41.0, 43.0, 45.0 and 47.0 GHz (Q band). Eachbaseband contains 16 sub-bands of 128 MHz, with 64 channels of 2 MHz, which translates into a velocity resolution of 23–33 km s−1 in the K band, and 12.5–15.0 kms−1 in the Q band. The phase reference center of the observations was set to the position α(J 2000) = 18h 47m 34. s 315, δ(J2000) = −01° 12′ 45.′′ 90. Flux, bandpass, and gain calibrations were achieved through observations of 1331+305, 1733-1304 (NRAO530), and 1851+0035, respectively. The gain-calibration cycles lasted 3 and 1.5 min for the K-band and Q-band observations, respectively. The data were calibrated and imaged using the VLA Pipeline CASA software package6, which is designed to perform automatic flagging and calibration for Stokes I continuum data. Conservatively, the absolute accuracy of the flux density scale is estimated to be ~5% at both wavelengths (Perley & Butler 2017). Further imaging and analysis was done with the GILDAS package. The resulting synthesized beam and rms noise of the continuum maps are given in Table 1.

By comparing the VLA and ALMA maps we realized that there was a systematic positional offset between them. Since both observations used the same phase calibrator 1851+0035 (J1851+0035), we cross-checked its coordinates in the two datasets and found an offset of Δα =0.′′021, Δδ = −0.′′0504. After contacting the VLA and ALMA staff, we reached the conclusion that the coordinates of the calibrator used by ALMA were more accurate than those used by the VLA. Therefore, the phase reference center of the VLA maps was shifted accordingly to the measured offset. After applying this correction, no significant offset was observed between the ALMA and VLA data.

The VLA A-Configuration at K band resolves out the more extended emission from the UC HII region located ~ 5″ away from the HMC (Cesaroni et al. 1994); owing to sidelobe effects and poor cleaning, spurious extended features contaminate the K-band image of the HMC. First, we filtered out this spurious extended emission by cutting away all baselines with length < 5 km. Then, we tried to improve the dynamical range of the K-band map by self-calibrating the emission of the two compact sources inside the HMC. Phase-only self-calibration improved the dynamical range of the image only marginally (a few percent), while the signal is too weak for amplitude self-calibration. Thus, our final K-band image has no self-calibration applied. We also attempted to self-calibrate the compact emission at Q band, but the signal-to-noise ratio is not sufficient even for phase-only self-calibration.

thumbnail Fig. 1

Left: ALMA map of the 1.3 mm continuum emission of the G31 HMC from Beltrán et al. (2019). The contours are − 5, 5, 10, 15, 20, 40, 80, 160, and 300 times 1σ, which is 1.2 mJy beam−1. The labels indicate the two cores, NE and Main, resolved by B18. The red cross indicates the position of the UC HII region imaged by Cesaroni et al. (1994). The synthesized beam of ~ 0.′′24 is shown in the lower left corner. Right: close-up of the central region around cores NE and Main that shows the ALMA 1.3 mm continuum emission image overlaid with the ALMA 1.4 mm continuum emission (white contours) observed with an angular resolution of ~ 0.′′1. White contours are 3, 6, 9, 15, 30, 60, 90, 120, and 165 times 1σ, which is 0.32 mJy beam−1. The ALMA 1.3 mm and 1.4 mm synthesized beam is shown in the lower left and right corner, respectively.

3 Results and analysis

3.1 Continuum emission of the Main and NE cores

Figure 1 (left panel) shows the dust continuum emission of the G31 HMC at 1.3 mm obtained by Beltrán et al. (2019) in their dust polarized emission observations carried out with ALMA at an angular resolution of ~ 0.′′24, which have clearly resolved the emission of G31 into two cores, Main and NE. In this figure, the Main core shows a monolithic and homogeneous appearance. The right panel shows the new high angular resolution (~ 0.′′1) observations at 1.4 mm overlaid on the ~ 0.′′24 observations. In the higher angular resolution map the emission from both cores is clearly detected and the emission of the Main core starts to show a less homogeneous morphology revealing the presence of clumpiness. In addition, the observationsalso reveal an elongation to the north of the Main core that was already visible but not previously resolved.

The structure of the Main core is better highlighted at 3.5 mm as shown in Fig. 2. At this wavelength and with slightly higher angular resolution (~ 0.′′075), the emission is clearly resolved into four embedded cores, labeled A, B (following the nomenclature of C10), C, and D. To avoid confusion between the larger cores observed with an angular resolution of ~ 0.′′24 and the cores resolved with higher angular resolution, hereafter we use the term “source” for the latter. Sources A and B, which are located at the center of the core and are separated by ~ 0.′′20 (~ 750 au), coincide with the two free-free continuum emission sources first detected at 1.3 cm and 7 mm by C10. The morphologies of these two sources are quite similar to each other. Source C lies ~ 0.′′30 (~ 1100 au) east of source B, while source D is located ~ 0.′′45 (~ 1700 au) south of source A. Source D is located close to the intersection of two of the SiO outflows, the N–S and the E–W, detected by B18 (see Fig. 2). The source is located along the direction of the E–W SiO outflow and close to its center of symmetry, which suggests that source D is likely its driving source. The NE core and the northern elongation are marginally visible at 3.5 mm.

Figure 3 shows the continuum emission at 1.3 cm and 7 mm obtained with new VLA observations, which in addition to detecting sources A and B observed by C10, also detected the emission from sources C and D at 7 mm, while at 1.3 cm only source D was detected. At 1.3 cm, the emission of the sources is quite compact, with the exception of source A that shows a central compact component plus an elongation to the south. This elongation was marginally visible with the 1.3 cm observations of C10. Figure 4 shows the VLA 7 mm and 1.3 cm continuum emission overlaid on the ALMA 1.4 mm and 3.5 mm continuum emission. The agreement in the position of the sources at the different wavelengths is remarkable. We note that neither the NE core nor the elongation observed northward of the Main core at 1.4 mm and 3.5 mm have been detected at 7 mm or 1.3 cm.

thumbnail Fig. 2

Overlay of the ALMA 3.5 mm continuum emission (white contours) on the 1.4 mm continuum emission (colors) toward cores NE and Main in G31. The white contours are 3, 6, 10, 15, 20, 30, and 40 times 1σ, which is 0.16 mJy beam−1. The 1.4 mm and 3.5 mm synthesized beams are shown in the lower left and right corner, respectively. The names of the sources in which the emission of the Main core has been resolved (see Sect. 3.1 and Table 2) are indicated. Dashed black lines indicate the direction of the E–W and N–S outflows mapped in SiO by B18.

thumbnail Fig. 3

Overlay of the VLA 1.3 cm continuum emission (contours) on the 7 mm continuum emission (colors) toward cores NE and Main in G31. The dashed black contours are −3 and the white contours 3, 6, 10, 15, 20, 30, 50, and 80 times 1σ, which is 0.005 mJy beam−1. The 1.3 cm and 7 mm synthesized beams are shown in the lower left and right corner, respectively. The names of the cores in which the emission of the Main core was resolved are indicated.

thumbnail Fig. 4

Top panels: overlay of the VLA 7 mm (left) and 1.3 cm (right) continuum emission (contours) on the ALMA 1.4 mm continuum emission (colors) toward cores NE and Main core in G31. The black contours are −3 and the white contours 3, 6, 12, 30, 50, and 80 times 1σ, which is 0.025 mJy beam−1 at 7 mm and 0.005 mJy beam−1 at 1.3 cm. The 1.4 mm synthesized beam is shown in the lower left corner of the panels, and the 7 mm and 1.3 cm synthesized beams in the lower right corner. The names of the sources in which the emission of the Main core has been resolved are indicated. Bottom panels: same but overlaid on the ALMA 3.5 mm continuum emission (colors).

3.2 Parameters of the sources embedded in the Main core

Table 2 lists the position, flux density, and size of the four sources embedded in the Main core of G31 at the four different wavelengths. The positions of the emission peaks of each source at the different wavelengths coincide. Furthermore, sources A and B are the brightest at all wavelengths, and have similar peak intensities at millimeter wavelengths, while at 1.3 cm source A is twice as bright as source B. The integrated flux density of the sources at the different wavelengths was estimated by fitting elliptical 2D Gaussians in the image plane. The spectral energy distributions (SEDs) of the sources are shown in Fig. 5. At 1.4 mm, 7 mm, and 1.3 cm the integrated flux densities of sources A and B are clearly greater than those of sources C and D, but at 3.5 mm the emission of source C is stronger than that of sources A and B.

The size of the sources was estimated by fitting elliptical 2D Gaussians in the image plane. At the highest angular resolution achieved, that is ~ 0.′′05 (or ~180 au at the distance of G31) at 7 mm, all the sources were resolved and have deconvolved geometric mean diameters of ~ 200–350 au. At any wavelength but 1.3 cm the sizes of sources A and B are comparable. At 1.3 cm, the major axis of source A is greater than that of source B, probably due to the presence of the southern elongation visible in the maps (see Fig. 3), which could be associated with source A or could indicate the presence of an additional source. Therefore, the similarities between sources A and B are not only related to their peak intensities and integrated flux densities, but also to their morphologies and sizes.

3.2.1 Free-free contribution to the emission

As already mentioned, C10 observed the region at 3.6 cm, 1.3 cm, and 7 mm and resolved the emission into two sources named A and B. These authors discussed the nature of their continuum emission in terms of HII regions or thermal radio jets, and discarded the former hypothesis. An important outcome of the new data is the good match between the positions of the centimeter and the millimeter peaks revealed by the new VLA and ALMA observations (Fig. 4). Such a good match proves that the free-free emission from sources A and B is tracing two distinct (jets from) young stellar objects (YSOs) and not the lobes of the same radio jet. In the latter case one should see only one millimeter peak located right in the middle of the two free-freepeaks.

Cesaroni et al. (2010) fitted the SED of sources A and B assuming that the emission originates from thermal radio jets with a power law Sννα, where the spectral index α is 0.9 for source A and 1.3 for source B. These spectral indices are consistent with the values found for partially thick thermal radio jets (e.g., Anglada et al. 2018). The integrated flux emission, Sν, at 1.3 cm estimated by C10 for sources A and B is 0.76 and 0.46 mJy, respectively, and is consistent with the values measured by us (see Table 2). Therefore, using the same spectral index of 0.9 for source A and 1.3 for source B and extrapolating the 1.3 cm fluxes to 7 mm, 3.5 mm, and 1.4 mm, we estimated the expected free-free emission of the sources. We used the 1.3 cm flux instead of the 3.6 cm because the emission of the sources is better resolved and the angular resolution at 1.3 cm, with a synthesized beam of 0.′′095 × 0.′′084 (C10) is comparable to that of the ALMA observations (0.′′075–0.′′1). Table 3 shows the expected flux densities at the different wavelengths as well as the contribution of the free-free emission to the total emission. The contribution of the free-free emission for sources A and B is important at 7 mm and negligible at 3.5 mm and 1.4 mm. Figure 5 shows the SEDs of the sources, with a dashed line indicating the extrapolation of the 1.3 cm fluxes with a power law with an index of 0.9 for source A and of 1.3 for source B.

The sensitivity of our VLA observations is >6 times higherat 1.3 cm and >3 times at 7 mm than that of C10, and this allowed us to detect source D at both wavelengths, and source C at 7 mm. Withthe current data, it is not possible to estimate the spectral index of the free-free emission for source C. For source D, even though the source has been detected at both 1.3 cm and 7 mm, an accurate estimate of the spectral index is not possible because the emission at 7 mm could be contaminated by thermal dust emission. We therefore estimated the expected free-free flux by assuming two different spectral indices: α =2, for which the contribution of the free-free emission at millimeter wavelengths should be maximum, and α = 0.9, the lowest spectral index estimated by C10 for the sources embedded in the Main core, which should give the lowest free-free contribution at millimeter wavelengths. We note that if the spectral index were < 0.9, the free-free contributionat millimeter wavelengths would be even smaller. Table 3 gives the expected free-free fluxes and the contribution to the total flux emission, while in Fig. 5, we show with a dotted and a dashed lines the extrapolations of the 1.3 cm emission for α = 2 and 0.9, respectively. Taking into account that source D is probably the driving source of the SiO E–W outflow detected by B18, then its free-free emission is likely due to a thermal radio jet, which has α <1.3 (see Anglada et al. 2018 and references therein). In any case, independently of the spectral index used to extrapolate the centimeter emission to millimeter wavelengths, the expected free-free contribution at 1.4 mm and 3.5 mm is negligible.

In conclusion, while the continuum emission at 7 mm still has an important contribution of free-free emission, especially for sources A and B, the emission at 3.5 mm and 1.4 mm is dominated by thermal dust emission for all sources.

Table 2

Position, flux densities, peak brightness temperature, and size of the sources embedded in the Main core of G31.41+0.31.

thumbnail Fig. 5

Spectral energy distribution of the four continuum sources detected in the Main core of G31. The dashed lines indicate the extrapolation of the fluxes at 3.6 cm or 1.3 cm with a power law of the type Sννα. The values of α are given in each box.

Table 3

Free-free continuum flux densities.

3.2.2 Mass estimates

The mass of the sources embedded in the Main core was estimated assuming that the emission is optically thin by using Mgas=gSνd2κνBν(Td),\begin{equation*} M_{\textrm{gas}} = \frac{g\,S_{\nu}\,d^2}{\kappa_{\nu}\,B_{\nu}(T_{\textrm{d}})} ,\end{equation*}(1)

where Sν is the flux density, d the distance to the source, κν the dust opacity coefficient, g the gas-to-dust ratio, and Bν (Td) the Planck function for a blackbody with dust temperature Td. We used thecontinuum emission at 3.5 mm because at this wavelength the emission is optically thinner than at 1.4 mm. Based on the findings of Sect. 3.2.1, we assumed that the continuum emission at 3.5 mm is pure dust thermal emission and estimated the masses Mgas, the mean H2 column densities NH2$N_{\textrm{H}_2}$, the mean H2 volume densities nH2$n_{\textrm{H}_2}$, and the mass surface density Σ of the four sources associated with the Main core (see Table 4). We adopted a dust opacity of 0.2 cm2 g−1 at 3.5 mm (by extrapolating, with a dust opacity index β = 1.8, the value at 1.3 mm of Ossenkopf & Henning 1994 for a Mathis et al. 1977 (MRN) distribution with thin ice mantles and a gas density of 108 cm−3) and a gas-to-dust mass ratio of 100. Recent works propose higher values of the gas-to-dust mass ratio, for example 150 (Draine 2011) or 162 (Peters et al. 2017). It is possible to estimate the masses and densities given in Table 4 with a different gas-to-dust mass ratio simply by multiplying the values by (new gas-to-dust mass ratio)/100. Assuming that gas and dust are coupled, we used the excitation temperature Tex estimated for methyl formate, CH3OCHO v = 0 and vt = 1 by B18 as Td. These authors fitted the CH3OCHO emission as a function of radius across the Main core, thus obtaining a temperature profile. Therefore, from the distance of the source to the center of the Main core, it is possible to estimate Tex using their Eq. (1). The temperature profile of B18 is valid for radii ≥0.′′22 (the angular resolution of their ALMA observations); therefore, we assumed that sources A and B, located at 0.′′ 22 and 0.′′ 04 from the central position of the temperature profile, have the same Tex of ~ 500 K, estimated at 0.′′22. Because the temperature profile was obtained with an angular resolution ~ 3 times lower than that of the 3.5 mm observations, part of the gas inside the beam could be slightly colder and Tex should be taken as a lower limit. This would decrease the mass and density estimates. In order to have an idea of the uncertainty introduced by the previous assumption on Tex, we recomputed the masses and densities increasing the temperatures by 100 K.

The masses of the sources range from ~15 to ~ 26 M (Table 4), suggesting that the cores are probably massive enough to form high-mass young stellar objects. The H2 column densities are higher than 1025 cm−2 for all sources. These values correspond to visual extinctions of > 5000 mag, confirming that the sources are deeply embedded in the core. The volume densities are also very high, several 109 cm−3. Using the NH2$N_{\textrm{H}_2}$ values we estimated the dust optical depth at 3.5 mm, and obtained values between 0.5 and 0.7 (see Table 4). These opacity values indicate that at 3.5 mm the dust emission is mostly optically thin and that the masses should not be largely underestimated.

Table 4

Parameters of the sources embedded in the Main core.

3.3 Other sources in the region

Figure 6 shows the 8.′′5 central region of the G31 HMC at the different wavelengths. The UC HII region located ~ 5″ northeast of the Main core is clearly detected at 1.3 cm, 7 mm, and 3.5 mm, albeit partially resolved out. The UC HII region appears spherical, with a diameter of ~ 1″ or 3750 au (0.018 pc). According to C10, the star ionizing the UC HII region is a O6 star.

In addition to the UC HII region and the sources associated with the Main core and the NE core, our observations at millimeter wavelengths revealed the presence of other sources in the region. The sources were identified by eye. We consider a source real when detected at least at a 3σ level and attwo different wavelengths. This last constraint has the purpose of avoiding the spurious detection of weak sources that are actually noise. The mass sensitivity estimated at 3.5 mm for a 3σ detection limit is 1.5–4.0 M, for a range of temperature of 20–50 K. At 1.4 mm the 3σ mass sensitivity is 0.1–0.3 M for the sametemperature range.

Following what was done for source A at 1.3 cm, we decided not to split the sources when elongated (e.g., source NE or source 4), unless their emission clearly shows well-separated peaks. The possible multiplicity of these sources will be investigated in greater detail with future observations at higher angular resolution. The positions of these other sources, the wavelengths at which they have been detected, and their flux densities at 1.4 mm and 3.5 mm can be found in Table 5. The NE core is included in the table. We note that none of these sources are detected at centimeter wavelengths. The total number of sources detected in G31 is 12 plus the UC HII region, but taking into account that some sources could actually be multiples, this number should be taken as a lower limit.

Table 5 shows the gas masses of these seven additionally detected sources. The masses were estimated from the continuum at 3.5 mm except for two sources not detected at this wavelength for which the mass has been estimated at 1.4 mm. We adopted a dust opacity of 1.0 cm2 g−1 at 1.4 mm (Ossenkopf & Henning 1994 for a MRN distribution with thin ice mantles and a gas density of 108 cm−3). For these sources, there are no temperature estimates. Therefore, to estimate the masses we used a range of temperatures from 50 to 100 K. The values of the masses (Table 5) suggest that all but source 4 and the NE core may form young stellar objects of low to intermediate mass.

4 Discussion

4.1 No monolithic core

The homogeneous, round, and seamless morphology of the Main core in G31 revealed by the ~ 0.′′2 observationsat 1.4 mm of B18 and 1.3 mm of Beltrán et al. (2019) was been confirmed by our higher angular resolution observations at 1.4 mm (0.′′1), 3.5 mm (0.′′075), and 7.0 mm (0.′′047). The Main core was fragmented into at least four sources within ~ 1″. This confirms what was already suggested by the detection of the free-free continuum sources close to the center (C10), that the apparent homogeneity of the core is due to a combination of the large opacity of the 1.3 mm and 1.4 mm dust continuum emission and insufficient angular resolution. Therefore, the new observations have allowed us to discriminatebetween the two scenarios proposed for the G31 Main core (see Sect. 1), and suggest that G31 is a dynamically collapsing core undergoing global fragmentation. This conclusion is consistent with the detection of red-shifted absorption toward the sources embedded in the core (e.g., Girart et al. 2009; B18, and Fig. 7), which suggests accretion of material from the common envelope. The analysis of the line emission, and in particular of the kinematics of the embedded sources, will be the subject of a forthcoming paper.

The brightest sources in the core are A and B, and are located close to its center. These two sources are also the closest to each other in the Main core, being separated by ~ 0.′′2 or ~750 au in projection. Their peak intensities are greater, at some wavelengths by a factor ≳5, than those of sources C and D. Sources A and B have very similar SEDs (see Fig. 5), with source A being slightly stronger at 3.5 mm, 7 mm, and 1.3 cm. Similarities between the two sources can also be found in their sizes as shown in Table 2. The diameter of the sources at 3.5 mm, where the emission is tracing mostly thermal dust emission, is ~ 800 au. We note that the size of the sources at the highest angular resolution achieved with our observations (0.′′ 047 at 7 mm) is only ~200 au. However, as discussed in Sect. 3.2.1, because the emission at 7 mm is highly contaminated by free-free emission, this should be considered the size of the thermal radio jets associated with the sources and not the size of the dusty cores. The gas masses and densities of sources A and B are also similar (Table 4). The masses (~ 15 M) are consistent with the sources being associated with high-mass (proto)stars. This is confirmed by their high mass surface densities (~ 270 g cm−2; see Table 4) and the detection of infalling gas toward them. The similar properties and the relative small separation of sources A and B might indicate a common origin. The sources could be members of a binary system with a gas mass ratio close to 1. Photometric and spectroscopic observations of more evolved OB stars indicate a high rate (> 50%) of close binaries (Chini et al. 2013), with similar masses of the binary members. Hydrodynamical simulations of star-formingclusters predict a high fraction of binaries with similar masses for systems with separations < 10 au (e.g., Lund & Bonnell 2018; Clarke 2007 and references therein), much smaller than that between source A and B. On the other hand, 3D radiation-hydrodynamic simulations of high-mass star formation by Krumholz et al. (2009) predict the formation of massive binary systems at the center of the collapsing cores with separationss between thebinary members similar to that of sources A and B, and with a final stellar mass ratio of ~ 0.7. Therefore, these observational and theoretical works suggest that massive binary systems have mass ratios close to 1, while our results suggest that this mass ratio might be set already at the very early evolutionary stages of high-mass star formation. We note that this conclusion holds if not only the gas masses, but also the (proto)stellar masses of the objects embedded in sources A and B are similar, as suggested by the simulations of Bate (2000). This study predicts that regardless of their initial mass ratio, massive binaries will rapidly evolve toward a mass ratio of 1. In order to test the hypothesis of binarity in G31, it is necessary to analyze the line emission toward the sources to search for dynamical features that link them or reveal infalling material from a circumbinary disk.

The other two sources embedded in the Main core, sources C and D, are also massive, as indicated by their gas masses of ~ 26 M. These two sources are more extended, especially at 3.5 mm, and colder than sources A and B (Table 4), which suggests that they could be in a slightly earlier evolutionary phase. The enhanced surface densities and top-heavy mass function of the fragments in the Main core are similar to the properties predicted by the hub-filament system model (Kumar et al. 2020) because hubs have amplified densities and low virial parameters.

thumbnail Fig. 6

Overlay of the ALMA 1.4 mm (a) and 3.5 mm (b), and VLA 7 mm (c) and 1.3 cm (d) continuum emission (contours) on the ALMA 1.3 mm continuum emission of the G31 HMC observed by Beltrán et al. (2019). White contour levels are the same as in Figs. 1, 2, and 4. Black contours are −3 times 1σ. The 1.3 mm synthesized beam is shown in the lower left corner and that at the other wavelengths in the lower right corner of each panel. The red cross gives the position of the UC HII region imaged by Cesaroni et al. (1994), while the white crosses indicate the positions of the four sources embedded in the Main core. The labels of the other sources detected in the region are indicated in red (see Sect. 3.3) and their positions given in Table 5.

Table 5

Positions, flux densities, and masses of other sources in the region.

thumbnail Fig. 7

Bottom panel: overlay of the ALMA 1.4 mm continuum emission (white contours) on the integrated intensity (moment 0) map (colors) of CH313$_3^{13}$CN K = 3 at an angular resolution of ~0.′′1. The contours are the same as in Fig. 1 (right panel). The synthesized beam isshown in the lower left corner. The white crosses indicate the positions of the NE core and the four sources embedded in the Main core. Top panel: spectra of the CH313$_3^{13}$CN K = 3 and CH3CNK = 6 transitions toward the 1.4 mm continuum emission peak of sources A (black) and B (red). The red-shifted absorption is clearly detected in both transitions. The vertical blue dashed line indicates the systemic LSR velocity of 96.5 km s.

thumbnail Fig. 8

Magnetic field segments (red lines) toward the G31 HMC, obtained by rotating by 90° the polarization segments at 1.3 mm from Beltrán et al. (2019), overlapped on the ALMA 1.4 mm continuum emission map (contours). The segments are sampled every ten pixels. The contours are the same as in Fig. 6a. Negative contours are dotted black. The synthesized beam is shown in the lower left corner. The blue and red crosses indicate the position of the millimeter and/or centimeter sources in the region.

4.2 Fragmentation and magnetic field

From the analysis of the 1.3 mm dust polarization maps (see Fig. 8), the G31 HMC represents one of the clearest examples to date of an hourglass magnetic field morphology in the high-mass regime (Girart et al. 2009), with the magnetic field oriented southeast to northwest that maintains its morphology down to scales of < 1000 au (Beltrán et al. 2019). The central region of the core was found to be moderately magnetically supercritical, with a mass-to-flux ratio λ about 1.4–2.2 above the critical value and an estimated strength of the magnetic field in the range 8–13 mG (Beltrán et al. 2019). The alignment of sources A, B, and C, D along a direction NE–SW roughly perpendicular to the axis of the magnetic field also supports a scenario of magnetically regulated fragmentation. The position of the four sources in the Main core, found close to the “waist” of the hourglass field, suggests the possibility that they formed from the fragmentation of a massive pseudo-disk created by the magnetic focusing of the collapsing gas (Galli & Shu 1993). Given the inclination of the magnetic axis of the hourglass with respect to the plane of the sky of ~44° (Beltrán et al. 2019), the relative positions of the four sources are consistent with formation in a large, inclined, flattened dense structure at the center of the G31 Main core.

Magnetohydrodynamical (MHD) simulations of the collapse of rotating, magnetized clouds (e.g., Hosking & Whitworth 2004; Machida et al. 2005a,b; Price & Bate 2008; Peters et al. 2011) show that fragmentation is generally inhibited bydynamically important levels of magnetization: for moderate to high values of the mass-to-flux ratio λ (~ 2–14 times the critical value), the number of fragments formed during collapse is roughly reduced by half with respect to the pure hydrodynamical case, while for values of mass-to-flux ratio closer to the critical value fragmentation is completely suppressed and only one single object is formed (Hennebelle et al. 2011; Commerçon et al. 2011). However, the exact outcome of the process of fragmentation also depends on other factors, such as the density profile of the collapsing core or cloud (e.g., Girichidis et al. 2011), radiative and mechanical (jet or outflow) feedback from the newly born stars, turbulence, and the amount of rotation of the initial configuration. A notable feature of the G31 Main core is its relatively high angular rotation, on the order of Ω ~5 × 10−12 s−1, if the velocity gradient of ~ 160 km s−1 pc−1 observed by B18 in the CH313$_{3}^{13}$CN(12–11) K = 2 transition is interpreted as solid-body rotation. With a mean density of ~ 4 × 107 cm−3, estimated from the mass of 70 M (Cesaroni 2019) and a radius of ~ 1.′′076 or ~4000 au (B18), the characteristic free-fall time of the Main core is tff ~ 104 yr, which implies Ω tff ~ 1. According to the MHD collapse simulations by Machida et al. (2005a,b), during the collapse of a model cloud with levels of magnetization and rotation similar to those of the G31 Main core (their model D: λ = 1.4, Ω tff = 3) any initial non-axisymmetric perturbation is strongly amplified, and the dense core of size ~ 100 au formed at the end of the isothermal phase is prone to bar instability and rotationally induced fragmentation. This is similar to the ionization-MHD simulations of Peters et al. (2011). To establish whether this scenario can be applied to the case of G31, the enigmatic interplay between rotation and magnetic fields in the evolution of the G31 Main core (Girart et al. 2009; B18) deserves consideration by further studies.

Further indications of a magnetically controlled fragmentation come from considering the Jeans mass of the Main core. Assuming that the temperature of the Main core before the formation of the embedded sources was ~ 20 K, as suggested by observations of high-mass starless cores (e.g., Fontani et al. 2015), and assuming an initial uniform density similar to the mean density of the Main core of ~ 4 × 107 cm−3, the initial Jeans mass in the Main core would be ~ 0.1 M. The assumption of a density at the onset of fragmentation similar to the current value is based on the fact that the Jeans mass has a much stronger dependence on temperature than on density, and that the initial density of the Main core could not have been much lower than the current density. Moreover, the current temperature is the result of the heating of the gas by the newly formed stars, and the gas heating occurred on a timescale much shorter than the evolutionary timescale of the Main core that characterizes the density increase. Therefore, the density variation during the heating process should be negligible. Clearly, a cloud with these properties would have produced a large number of fragments of small mass, unlike what was observed in G31. This region appears to have turned into a power-law density profile by global gravitational contraction, dragging in the magnetic field initially present in the cloud, and producing in this process a relatively small number of massive fragments, mostly concentrated in the central part. Although highly speculative, these considerations support the view that the G31 Main core has undergone a phase of magnetically regulated fragmentation characterized by a reduced fragmentation efficiency, eventually leading to the formation of a small number of relatively massive dense cores, consistent with MHD simulations (e.g., Machida et al. 2005a,b; Commerçon et al. 2011; Peters et al. 2011; Myers et al. 2013).

Whether or not the process of fragmentation is now completed in the G31 Main core is difficult to assess. According to Krumholz &McKee (2008), radiative feedback effects suppress further fragmentation in cores with a surface density higher than ~ 1 g cm−2 as the radiation from the embedded protostars heats the gas providing efficient contrast to further local collapse. For the G31 Main core, the surface density estimated for a mass of 70 M and a radius of ~ 4000 au is quite high, ~ 12 g cm−2. This suggests that thermal (and possibly also mechanical) feedback from the embedded sources may have halted the fragmentation of the Main core.

thumbnail Fig. 9

Overlay of the ALMA 1.4 mm continuum emission (grayscale), the VLA combined C & D Array continuum emission at 1.3 cm (white contours) from Cesaroni et al. (1994), and the Herschel Hi-GAL 70 μm emission (black contours) on the Spitzer GLIMPSE 8.0 μm emission (colors). White contours range from 4 mJy beam−1 to 60 mJy beam−1 in steps of 7 mJy beam−1. Black contours range from 2000 MJy sterad−1 to 152 000 MJy sterad−1 in steps of 30000 MJy sterad−1.

4.3 The global picture of the G31.41+0.31 high-mass star-forming region

Figure 9 shows a large-scale image of the G31.41+0.31 high-mass star-forming region at different wavelengths: infrared, millimeter, and centimeter. The white contours show the continuum emission at 1.3 cm observed with the C & D configurations of the VLA by Cesaroni et al. (1994). The centimeter emission, which is dominated by the free-free emission of the UC HII region, is very extended and has a cometary shape with two peaks located close to the “head”. The strongest and most compact peak is located to the east and is likely tracing the position of the O6 star ionizing the UC HII region (C10), while the secondary and more elongated peak is located to the south close to the position of the G31 HMC. This secondary peak could indicate a zone where the UV photons escaping from the central part of the UC HII region ionize the dense material associated with the HMC, creating a sort of bright-rimmed ridge. The “tail” of the cometary free-free emission points northward, where the region is probably less dense. This is further supported by the fact that the tail of the free-free emission coincides with the peak of the infrared emission at 8.0 μm from Spitzer, which suggests that the photons from the stars embedded in the G31 HMC escape through the lower-density region. The peak of the infrared emission at 70 μm, as seen by Herschel, is also located north of the G31 HMC, but closer to it than the peak of the 8 μm emission, consistent with the lower dust opacity at 70 μm. This seems to confirm that the heating sources are indeed embedded inside the HMC.

Zooming in toward the position of the HMC, we find a small (proto)cluster of millimeter sources, probably in an earlier evolutionary stage with respect to the UC HII region. This small (proto)cluster consists of 12 young stellar objects, mostly concentrated in a small region of ~ 5″ around the central position of the Main core, as seen in Fig. 6. In addition to the four sources embedded in the Main core (A, B, C, and D), there are six additional sources (1, 2, 3, 4, 5, and the NE core) very close to it and located in what appears as streams or filaments of material pointing to the Main core. This would suggest a scenario similar to that described by the competitive accretion model of Bonnell & Bate (2006), that predicts that a molecular cloud initially fragments in low-mass cores of Jeans mass that form stars that compete to accrete mass from the common gas reservoir. Protostars located near the center of the gravitational potential accrete at a higher rate because of a stronger gravitational pull, and thus experience a faster mass growth. It would also be consistent with the “fragmentation-induced starvation” scenario of Peters et al. (2010). According to this model, the dense accretion flow onto the central massive object becomes gravitationally unstable leading to fragmentation and formation of secondary objects in the flow that starve the central star of accreting material. In G31, at the center of the cloud, we find a small cluster of four massive young stellar objects, sources A, B, C, and D, still actively accreting, as demonstrated by the detection of red-shifted absorption and molecular outflows associated with them.

From an observational point of view, the large-scale scenario in the G31 star-forming region resembles that observed in some of the densest and most massive HMCs forming protoclusters: e.g., Sgr B2(N) (Sánchez-Monge et al. 2017, 2019) and W51 (Goddi et al. 2020), where complex dusty streams and filaments converging onto the central compact sources have been detected. These filaments are dense and relatively small, with spatial scales of ~ 0.1 pc for Sgr B2(N) (Sánchez-Monge et al. 2017) and <2000 au for W51 (Goddi et al. 2020). Sánchez-Monge et al. (2019) have analyzed the velocity field along the filaments in Sgr B2(N) and found velocity gradients consistent with accretion to the center of the cores, similarly to what Peretto et al. (2013) have found in much younger (infrared-dark) high-mass star-forming cores. For W51, Goddi et al. (2020) have not been able to measure thegas kinematics along the dusty streams because of the large opacity of the 1.3 mm dust continuum emission. In G31 we searched for velocity gradients along the dust streams in which sources 1, 2, 3, 4, 5, and the NE core seem to be located to identify kinematic signatures of accretion of material onto the central region from the larger-scale cloud. Unfortunately, the current high angular resolution line data do not have enough sensitivity to allow us to study the velocity field, and we cannot completely confirm that this scenario is in place in G31.

5 Conclusions

We carried out ALMA observations at 1.4 mm and 3.5 mm of the high-mass star-forming region G31 at angular resolutions of ~ 0.′′1 (~ 375 au) and ~ 0.′′075 (~ 280 au), respectively.The goal of these observations was to establish whether the Main core in G31 is as monolithic as it appears at ~ 1 mm and ~ 0.′′2 angular resolution (B18, Beltrán et al. 2019) or whether its homogeneous appearance is due to a combination of both large dust opacity and insufficient angular resolution. We also carried out VLA observations at 7 mm and 1.3 cm at ~ 0.′′05 (~ 190 au) and ~ 0.′′07 (~ 260 au) angular resolutions, respectively, to better study the nature of the free-free continuum sources previously detected by C10.

The millimeter continuum emission of the Main core has been clearly resolved into at least four sources, A, B, C, D, within ~ 1″. Sources A, B, and D have also been detected at centimeter wavelengths. In all three cases, the free-free emission is likely associated with thermal radio jets. In addition to these four sources, eight more have been detected in the region within 12″, including the UC HII region located ~ 5″ northeast from the Main core. The UC HII has been detected at 1.3 cm, 7 mm, and 3.5 mm and shows a spherical and compact morphology, with a diameter of ~ 1″ or 3750 au (0.018 pc) at the distance of the region.

The masses of the four sources embedded in the Main core, estimated at 1.4 mm, range from ~ 15 to ~ 26 M and their number densities are several 109 cm−3. The deconvolved radii of the dust emission of the sources, estimated at 3.5 mm, range from ~ 400 to 500 au, while the radii of the free-free emission, estimated at 7 mm or 1.3 cm, range from ~ 100 to 175 au. The mass surface densities of all the sources are very high, with values of 250–350 g cm−2.

Sources A and B, separated by ~750 au, dominate the flux density of the Main core and are those located close to its center. The millimeter emission peak coincides with the centimeter emission peak in both sources, proving that the free-free emission from sources A and B is tracing two distinct jets from young stellar objects and not the lobes of the same radio jet. Sources A and B have similar SEDs, masses, and sizes, and both are undergoing collapse, as suggested by the detection of red-shifted absorption toward them (see Fig. 7). Altogether, this indicates that the sources are associated with high-mass (proto)stars. This is confirmed by the high mass surface densities. Sources C and D, which are massive as well, have a more extended morphology and lower temperatures and could be high-mass young stellar objects in an earlier evolutionary phase. Source D is associated with centimeter emission likely arising from a thermal radio jet.

The observations confirm that the collapsing Main core in G31 has undergone fragmentation and therefore that the previously assessed monolithic appearance of G31 is a consequence of both insufficient angular resolution and dust opacity of the core. The low level of fragmentation suggests that the G31 Main core has undergone a phase of magnetically regulated evolution characterized by a reduced fragmentation efficiency, eventually leading to the formation of a small number of relatively massive dense cores, as found by theoretical models and simulations.

Acknowledgements

We thank the anonymous referee for their useful comments that have improved the manuscript. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2013.1.00489.S, ADS/JAO.ALMA#2016.0.00223.S, and ADS/JAO.ALMA#2018.0.00252.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. V.M.R. acknowledges support from the Comunidad de Madrid through the Atracción de Talento Investigador Modalidad 1 (Doctores con experiencia) Grant (COOL: Cosmic Origins Of Life; 2019-T1/TIC-15379), and from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 664931. H.B. acknowledges support from the European Research Council under the Horizon 2020 Framework Program via the ERC Consolidator Grant CSF-648405. H.B. also acknowledges support from the Deutsche Forschungsgemeinschaft in the Collaborative Research Center (SFB881) “The Milky Way System” (subproject B1). C.G. acknowledges support by the ERC Synergy Grant “BlackHoleCam: Imaging the Event Horizon of Black Holes” (Grant 610058). R.K. acknowledges financial support via the Emmy Noether Research Group on Accretion Flows and Feedback in Realistic Models of Massive Star Formation funded by the German Research Foundation (DFG) under grant no. KU 2849/3-1 and KU 2849/3-2. M.S.N.K. acknowledges the support from FCT - Fundacão para a Ciência e a Tecnologia through Investigador contracts and exploratory project (IF/00956/2015/CP1273/CT0002)

References

  1. Anglada, G., Rodríguez, L. F., & Carrasco-González, C. 2018, A&ARv, 26, 3 [NASA ADS] [CrossRef] [Google Scholar]
  2. Bate, M. R. 2000, MNRAS, 314, 33 [Google Scholar]
  3. Beltrán, M. T. 2020, ArXiv e-prints [arXiv:2005.06912] [Google Scholar]
  4. Beltrán, M. T., & de Wit, W. J. 2016, A&ARv, 24, 6 [Google Scholar]
  5. Beltrán, M. T., Cesaroni, R., Neri, R., et al. 2004, ApJ, 601, L190 [Google Scholar]
  6. Beltrán, M. T., Cesaroni, R., Neri, R., et al. 2005, A&A, 435, 901 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Beltrán, M. T., Codella, C., Viti, S., Neri, R., & Cesaroni, R. 2009, ApJ, 690, L93 [Google Scholar]
  8. Beltrán, M. T., Cesaroni, R., Rivilla, V. M., et al. 2018, A&A, 615, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Beltrán, M. T., Padovani, M., Girart, J. M., et al. 2019, A&A, 630, A54 [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bonnell, I. A., & Bate, M. R. 2006, MNRAS, 370, 488 [NASA ADS] [CrossRef] [Google Scholar]
  11. Calcutt, H., Viti, S., Codella, C., et al. 2014, MNRAS, 443, 3157 [Google Scholar]
  12. Cesaroni, R. 2019, A&A, 631, A65 [CrossRef] [EDP Sciences] [Google Scholar]
  13. Cesaroni, R., Churchwell, E., Hofner, P., Walmsley, C. M., & Kurtz, S. 1994, A&A, 288, 903 [NASA ADS] [Google Scholar]
  14. Cesaroni, R., Hofner, P., Araya, E., & Kurtz, S. 2010, A&A, 590, A50 [Google Scholar]
  15. Cesaroni, R., Beltrán, M. T., Zhang, Q., Beuther, H., & Fallscheer, C. 2011, A&A, 533, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Cesaroni, R., Sánchez-Monge, Á., Beltrán, M. T., et al. 2017, A&A, 602, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Chini, R., Barr, A., Buda, L. S., et al. 2013, Central Eur. Astrophys. Bull., 37, 295 [Google Scholar]
  18. Clarke, C. J. 2007, IAUS 240, 337, eds. W. I. Hartkopf, E. F. Guinan, & P. Harmanec (Cambridge University Press) [Google Scholar]
  19. Commerçcon, B., Hennebelle, P., & Henning, Th. 2011, ApJ, 742, L9 [Google Scholar]
  20. Draine, B. T. 2011, in Physics of the Interstellar and Intergalactic Medium, Princeton Series in Astrophysics (Princeton: Princeton University Press) [Google Scholar]
  21. Fontani, F., Pascucci, I., Caselli, P., et al. 2007, A&A, 470, 639 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Fontani, F., Busquet, G., Palau, A., et al. 2015, A&A, 575, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Galli, D., & Shu, F. H. 1993, ApJ, 417, 220 [Google Scholar]
  24. Girart, J. M., Zhang, Q., Rao, R., & Estalella, R. 2009, Science, 324, 1408 [Google Scholar]
  25. Girichidis, P., Federrath, C., Banerjee, R., & Klessen, R. S. 2011, MNRAS, 413, 2741 [NASA ADS] [CrossRef] [Google Scholar]
  26. Goddi, C., Ginsburg, A., Maud, L. T., Zhang, Q., & Zapata, L. A. 2010, ApJ, 905, 25 [Google Scholar]
  27. Hennebelle, P., & Joos, M., et al. 2011, A&A, 528, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Hosking, J. G., & Whitworth, A. P. 2004, MNRAS, 347, 1001 [Google Scholar]
  29. Immer, K., Li,J., et al. 2019, A&A, 632, A123 [CrossRef] [EDP Sciences] [Google Scholar]
  30. Krumholz, M. R., & McKee, C. F. 2008, Nature, 451, 1082 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  31. Krumholz, M. R., Klein, R. I., McKee, C. F., et al. 2009, Science 323, 754 [Google Scholar]
  32. Kuiper, R., & Hosokawa, T. 2018, A&A, 616, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Kumar, M. S. N., Palmeirim, P., Arzoumanian, D., & Inutsuka, S. I. 2020, A&A, 642, A87 [EDP Sciences] [Google Scholar]
  34. Lund, K., & Bonnell, I. A. 2018, MNRAS, 479, 2235 [Google Scholar]
  35. Machida, M. N., Matsumoto, T., Tomisaka, K., & Hanawa, T. 2005a, MNRAS, 362, 369 [Google Scholar]
  36. Machida, M. N., Matsumoto, T., Tomisaka, K., & Hanawa, T. 2005b, MNRAS, 362, 382 [Google Scholar]
  37. Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 [NASA ADS] [CrossRef] [Google Scholar]
  38. Mayen-Gijon, J. M., Anglada, G., Osorio, M., et al. 2014, MNRAS, 437, 3766 [Google Scholar]
  39. McKee, C. F., & Tan, J. C. 2003, ApJ, 585, 850 [Google Scholar]
  40. Mininni, C., Rivilla, V. M., et al. 2020, A&A, 644, A84 [EDP Sciences] [Google Scholar]
  41. Motte, F., Bontemps, S., & Louvet, F. 2018, ARA&A, 56, 41 [Google Scholar]
  42. Myers, A. T., McKee, C. F., Cunningham, A. J., Klein, R. I., & Krumholz, M. R. 2013, ApJ, 766, 97 [Google Scholar]
  43. Osorio, M., Anglada, G., Lizano, S., & D’Alessio, P. 2009, ApJ, 694, 29 [Google Scholar]
  44. Ossenkopf, V.,& Henning, Th. 1994, A&A, 291, 943 [Google Scholar]
  45. Peretto, N., Fuller, G. A., Duarte-Cabral, A., et al. 2013, A&A, 555, A112 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Perley, R. A., & Butler, B. J. 2017, ApJS, 230, 7 [NASA ADS] [CrossRef] [Google Scholar]
  47. Peters, T., Klessen, R. S., Mac Low, M.-M., & Banerjee, R. 2010, ApJ, 725, 134 [Google Scholar]
  48. Peters, T., Banerjee, R., Klessen, R. S., & Mac Low, M.-M. 2011, ApJ, 729, 72 [Google Scholar]
  49. Peters, T., Zhukovska, S., Naab, Th., et al. 2017, MNRAS, 467, 4322 [Google Scholar]
  50. Price, D., J., & Bate, M. R. 2008, MNRAS, 385, 1820 [Google Scholar]
  51. Rivilla, V. M., Cesaroni, R., et al. 2017, A&A, 598, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Schilke, P., Schmiedeke, A., et al. 2017, A&A, 604, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Schilke, P., Ginsburg, A., Cesaroni, R., & Schmiedeke, A. 2018, A&A, 609, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Schilke, P., et al. 2019, A&A, 628, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Tan, J. C., Caselli, P., et al. 2014, in Protostars and Planets VI, 149, eds. H. Beuther, R. Klessen, C. Dullemond, & Th. Henning, (University of Arizona Press) [Google Scholar]

1

The Cycle 4 data were re-calibrated using an updated version of CASA (5.4.0-70) with the respective update of the ALMA pipeline that had improved features compared to the initial assessments made. Further manual changes were made in the reduction, especially the line-free selections for continuum reduction after the already improved findcont process of the pipeline. The combined reprocessing resulted in corrected baselines of the spectral windows and a confirmed flux scale with respect to the Cycle 6 data.

2

The CASA package is available at http://casa.nrao.edu/

4

The GILDAS package is available at http://www.iram.fr/IRAMFR/GILDAS

All Tables

Table 1

Summary of the continuum observations.

Table 2

Position, flux densities, peak brightness temperature, and size of the sources embedded in the Main core of G31.41+0.31.

Table 3

Free-free continuum flux densities.

Table 4

Parameters of the sources embedded in the Main core.

Table 5

Positions, flux densities, and masses of other sources in the region.

All Figures

thumbnail Fig. 1

Left: ALMA map of the 1.3 mm continuum emission of the G31 HMC from Beltrán et al. (2019). The contours are − 5, 5, 10, 15, 20, 40, 80, 160, and 300 times 1σ, which is 1.2 mJy beam−1. The labels indicate the two cores, NE and Main, resolved by B18. The red cross indicates the position of the UC HII region imaged by Cesaroni et al. (1994). The synthesized beam of ~ 0.′′24 is shown in the lower left corner. Right: close-up of the central region around cores NE and Main that shows the ALMA 1.3 mm continuum emission image overlaid with the ALMA 1.4 mm continuum emission (white contours) observed with an angular resolution of ~ 0.′′1. White contours are 3, 6, 9, 15, 30, 60, 90, 120, and 165 times 1σ, which is 0.32 mJy beam−1. The ALMA 1.3 mm and 1.4 mm synthesized beam is shown in the lower left and right corner, respectively.

In the text
thumbnail Fig. 2

Overlay of the ALMA 3.5 mm continuum emission (white contours) on the 1.4 mm continuum emission (colors) toward cores NE and Main in G31. The white contours are 3, 6, 10, 15, 20, 30, and 40 times 1σ, which is 0.16 mJy beam−1. The 1.4 mm and 3.5 mm synthesized beams are shown in the lower left and right corner, respectively. The names of the sources in which the emission of the Main core has been resolved (see Sect. 3.1 and Table 2) are indicated. Dashed black lines indicate the direction of the E–W and N–S outflows mapped in SiO by B18.

In the text
thumbnail Fig. 3

Overlay of the VLA 1.3 cm continuum emission (contours) on the 7 mm continuum emission (colors) toward cores NE and Main in G31. The dashed black contours are −3 and the white contours 3, 6, 10, 15, 20, 30, 50, and 80 times 1σ, which is 0.005 mJy beam−1. The 1.3 cm and 7 mm synthesized beams are shown in the lower left and right corner, respectively. The names of the cores in which the emission of the Main core was resolved are indicated.

In the text
thumbnail Fig. 4

Top panels: overlay of the VLA 7 mm (left) and 1.3 cm (right) continuum emission (contours) on the ALMA 1.4 mm continuum emission (colors) toward cores NE and Main core in G31. The black contours are −3 and the white contours 3, 6, 12, 30, 50, and 80 times 1σ, which is 0.025 mJy beam−1 at 7 mm and 0.005 mJy beam−1 at 1.3 cm. The 1.4 mm synthesized beam is shown in the lower left corner of the panels, and the 7 mm and 1.3 cm synthesized beams in the lower right corner. The names of the sources in which the emission of the Main core has been resolved are indicated. Bottom panels: same but overlaid on the ALMA 3.5 mm continuum emission (colors).

In the text
thumbnail Fig. 5

Spectral energy distribution of the four continuum sources detected in the Main core of G31. The dashed lines indicate the extrapolation of the fluxes at 3.6 cm or 1.3 cm with a power law of the type Sννα. The values of α are given in each box.

In the text
thumbnail Fig. 6

Overlay of the ALMA 1.4 mm (a) and 3.5 mm (b), and VLA 7 mm (c) and 1.3 cm (d) continuum emission (contours) on the ALMA 1.3 mm continuum emission of the G31 HMC observed by Beltrán et al. (2019). White contour levels are the same as in Figs. 1, 2, and 4. Black contours are −3 times 1σ. The 1.3 mm synthesized beam is shown in the lower left corner and that at the other wavelengths in the lower right corner of each panel. The red cross gives the position of the UC HII region imaged by Cesaroni et al. (1994), while the white crosses indicate the positions of the four sources embedded in the Main core. The labels of the other sources detected in the region are indicated in red (see Sect. 3.3) and their positions given in Table 5.

In the text
thumbnail Fig. 7

Bottom panel: overlay of the ALMA 1.4 mm continuum emission (white contours) on the integrated intensity (moment 0) map (colors) of CH313$_3^{13}$CN K = 3 at an angular resolution of ~0.′′1. The contours are the same as in Fig. 1 (right panel). The synthesized beam isshown in the lower left corner. The white crosses indicate the positions of the NE core and the four sources embedded in the Main core. Top panel: spectra of the CH313$_3^{13}$CN K = 3 and CH3CNK = 6 transitions toward the 1.4 mm continuum emission peak of sources A (black) and B (red). The red-shifted absorption is clearly detected in both transitions. The vertical blue dashed line indicates the systemic LSR velocity of 96.5 km s.

In the text
thumbnail Fig. 8

Magnetic field segments (red lines) toward the G31 HMC, obtained by rotating by 90° the polarization segments at 1.3 mm from Beltrán et al. (2019), overlapped on the ALMA 1.4 mm continuum emission map (contours). The segments are sampled every ten pixels. The contours are the same as in Fig. 6a. Negative contours are dotted black. The synthesized beam is shown in the lower left corner. The blue and red crosses indicate the position of the millimeter and/or centimeter sources in the region.

In the text
thumbnail Fig. 9

Overlay of the ALMA 1.4 mm continuum emission (grayscale), the VLA combined C & D Array continuum emission at 1.3 cm (white contours) from Cesaroni et al. (1994), and the Herschel Hi-GAL 70 μm emission (black contours) on the Spitzer GLIMPSE 8.0 μm emission (colors). White contours range from 4 mJy beam−1 to 60 mJy beam−1 in steps of 7 mJy beam−1. Black contours range from 2000 MJy sterad−1 to 152 000 MJy sterad−1 in steps of 30000 MJy sterad−1.

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.