Open Access
Issue
A&A
Volume 640, August 2020
Article Number A111
Number of page(s) 8
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202038024
Published online 21 August 2020

© C. Arce-Tord et al. 2020

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

It has been proposed that the magnetic field, which is pervasive in the interstellar medium (ISM), might play an important role in the dynamical evolution of star-forming clouds (e.g., Shu et al. 1987; Hennebelle et al. 2011; Commerçon et al. 2011; Crutcher 2012; Planck Collaboration Int. XXXV 2016; Planck Collaboration Int. XXXIII 2016; Matsushita et al. 2018; Beuther et al. 2018). Large-scale observations (~1 pc) of the magnetic field revealed a well-ordered structure in the low-density envelopes of molecular clouds, suggesting that parsec-scale envelopes are magnetically supported against gravitational collapse (e.g., Franco et al. 2010). Besides their possible role in supporting clouds against gravity, magnetic fields may also strongly affect the formation and evolution of substructures within the clouds. Indeed the magnetic fields sometimes have a consistent morphology throughout scales, as reported by the TADPOL survey (Hull et al. 2014), which compares the orientation of the magnetic field at 20″ and 2.5″ in 30 low-mass star-forming regions distant of 125 to 2650 pc. At core scales (~ 0.1 pc), the magnetically-regulated core-collapse models presume a dominant role of the magnetic field (Shu et al. 1987, 2004; Galli & Shu 1993a; Tomisaka 1998; Allen et al. 2003a,b). These models are consistent with a handful of observations (e.g., Chapman et al. 2013; Qiu et al. 2014), but due to a small fraction of consistent data, the magnetically dominant core-collapse does not seem to be the predominant mode of low- or high-mass star formation (Hull & Zhang 2019). At the circumstellar disk scales, ~200 au, it is predicted that the magnetic fields have a perpendicular component to the disks, thus having a leading role for the triggering of outflows (Blandford & Payne 1982; Camenzind 1990; Konigl & Pudritz 2000; Shu et al. 2000). No systematic relation, however, has been observed thus far between the magnetic field direction and the outflow orientation when these are probed in the low-mass regime at ~1000 au scales (Hull et al. 2013; Hull & Zhang 2019). Nevertheless, one low-mass study toward 12 Class 0 protostars by Galametz et al. (2018) shows that at the envelope scale (600–1500 au), the magnetic field is preferentially oriented either parallel or perpendicular to the direction of the outflows.

A few interferometric studies, such as Girart et al. (2009), Beuther et al. (2010), and Sridharan et al. (2014), have probed the high-mass star-forming regime. In particular, Zhang et al. (2014a) performed multi-scale observations of the magnetic field orientation in 14 high-mass star-forming regions. They show that the magnetic field at core scales is either perpendicular or parallel to the magnetic field orientation at clump scales. Similar results have also been reported by Koch et al. (2014) toward 50 star-forming regions, and by Ching et al. (2017) toward six massive dense cores located in DR21. One issue in determining the importance of magnetic fields for high-mass star formation is the lack of consistent analysis, in a large sample of sources, of the magnetic field morphology with respect to the density structure and gas dynamics. We collected all of this information in the high-mass star-forming region (HMSFR) W43-MM1. W43-MM1, which is at a distance of 5.5 kpc from the Sun (Zhang et al. 2014b), has been reported as a mini-starburst cluster due to its high star formation rate of ~6000 M Myr−1 (Louvet et al. 2014). W43-MM1 has been observed in polarized dust emission at 1.3 mm by Cortes & Crutcher (2006) with BIMA. They reported an ordered polarization pattern at ~4.5″ angular resolution and derived an on-the-plane of the sky magnetic field strength of 1.7 mG. Later, Sridharan et al. (2014) observed W43-MM1 in polarized dust emission at 345 GHz at an angular resolution of ~2.3″ with the SMA, and they reported a pinched morphology of the magnetic field with a field strength of ~6 mG. More recently, Cortes et al. (2016) conducted the first ALMA high-angular resolution map at 0.5″ of the magnetic fields in W43-MM1, which was performed using a single pointing, and they also report an organized magnetic field morphology with field strengths ranging from 0.2 to 9 mG.

In this article, we present observations of the polarized dust emission toward W43-MM1 in five pointings with the ALMA interferometer. We compare the orientation of the plane-of-the-sky magnetic field with the orientation of the cores and outflows identified in W43-MM1. Motte et al. (2018) report that about 130 dense cores were forming stars in this region. From this sample of cores, 27 are driving 46 outflow lobes, whose orientations have been traced by Nony et al. (2020) using CO(2-1) and SiO(5-4) emission lines with ALMA. In Sect. 2 we describe the observations, calibration, and imaging of our data. In Sect. 3, we calculate the angle differences among the magnetic field, the orientations of the outflows, and the position angle (PA) of the cores; additionally, we build their statistical distributions in the form of cumulative distribution functions (CDFs) that are further compared with synthetic populations of angle distributions. In Sect. 4 we discuss the physical implications of our results and conclude in Sect. 5.

2 Observations

We present five pointings with ALMA in full polarization targeting the most massive dense cores in W43-MM1, which were previously identified by Louvet et al. (2014) based on IRAM/Plateau de Bure data. The ALMA observations were obtained between April and May 2016 using an array with 35 antennas. We used the standard frequency setup for continuum polarization in band 6 (i.e., four spectral windows centered at 224.984, 226.984, 239.015, and 241.015 GHz). Each spectral window has 64 channels of 31.250 MHz each, corresponding to a total bandwidth of 1.875 GHz per spectral window. The maximum recoverable scale is ~10.6″. The resulting continuum images have an equivalent frequency of 233 GHz (or 1.3 mm) with an angular resolution of 0.55″ × 0.49″ (or ~2700 au, considering the distance of the target). The calibration, the imaging, and the analysis were done using version 5.1.2. of the Common Astronomical Software Applications (CASA, McMullin et al. 2007).

To improve the signal-to-noise ratio (S/N) of our images, we performed a phase self-calibration on each of the five calibrated datasets. The images were improved iteratively through four rounds of phase-only self-calibration using the Stokes I image as a model. The final imaging was performed with the tclean task from CASA using a Hogbom deconvolver and a Briggs weighting scheme witha robust parameter of 0.5. The Stokes U and Q maps were imaged individually for each of the five fields with the tclean task using the same parameters as above. After the cleaning, a primary beam correction was applied to each map. Table 1 shows the coordinates for each pointing and the root mean square (rms) noise level for each Stokes parameter. To provide an overview of the region, we present in Fig. 1 a mosaicked Stokes I image1.

The magnetic field maps were computed using the Stokes Q and U intensity maps for each pointing. For this, we created the linear polarized intensity map, P = Q2+U2$\!\sqrt{Q^2 + U^2}$. However, even when the Stokes Q and U maps could have negative values, P is always positive. In order to correct for this, we applied the debiasing method proposed by Vaillancourt (2006), which employs a Bayesian approach instead of the classical frequentist method, and gives a better estimation of the polarized emission for low signal-to-noise measurements, ≲3σP, where σP is the rms noise level in the polarization maps. We wielded the 3σP debiased intensity value as a threshold to select the pixels used to derive the polarization angle maps, which were computed as χ = 0.5arctan(UQ). Finally, we assumed alignment of the grains with respect to the magnetic field due to the “radiative alignment torques” mechanism (RATs), through which an external radiation field causes the dust grain to spin-up, contributing to an efficient alignment between its angular momentum and the magnetic field line (Lazarian 2007; Hoang & Lazarian 2009; Andersson et al. 2015). As a consequence, the aligned grains emit thermal radiation, which is polarized perpendicular to the magnetic field, letting us infer the magnetic field morphology onto the plane of the sky by rotating the polarization angle by 90°.

3 Results and analysis

Figure 1 shows the 1.3 mm continuum emission of W43-MM1, which hosts two clusters of protostars separated by ~0.9 pc: the main region (pointings 1, 2, and 3) and the southwest (SW) region (pointings 4 and 5). The continuum emission shown here is very similar to that previously presented by Cortes et al. (2016) and Motte et al. (2018). The latter identified 131 cores in the two clusters. Those falling within the one-third area of the primary beams of our five pointings are superposed over Fig. 1. Figures 2 and 3 display the morphology of the magnetic field in each of the five pointings, which were imaged independently. In the main region, the large-scale magnetic field (~1 pc) shows a smooth and ordered morphology. In the SW region, the information about the magnetic field orientation is coarser: there are only eight independent locations where the magnetic field orientation could be derived. Linking the large-scale morphology of the magnetic field to its small-scale features goes beyond the scope of this work, and is the subject of a forthcoming publication (Louvet et al., in prep.). We measured the position angle (PA) of the magnetic field at the location of the cores identified by Motte et al. (2018). To ensure reliable linear polarization measurements, and following ALMA recommendations, we restrained our analysis to within one-third of the primary beam for each pointing. In total, 52 cores out of the 131 cores of Motte et al. (2018) fall into these trusted areas. Of these 52 cores, 29 display polarized thermal dust emission. We computed the magnetic field orientation for each core as BF = 0.5 arctan(UQ) + π∕2 (see Sect. 2), where Q and U are the mean Stokes Q and U intensities averaged within each core, whose area is defined by a 2D Gaussian fit of its Stokes I emission. The polarization angle uncertainties were estimated through error propagation as σχ = 0.5 × (σPP), where P is the mean polarized intensity within the corresponding core and σP = (<Q>×σQ)2+(<U>×σU)2/P$\sqrt{(<Q>\,{\times}\,\sigma_{\mathrm{Q}})^2 &#x002B; (<U>\,{\times}\,\sigma_{\mathrm{U}})^2 }/ P $ is the rms of the linearly polarized emission.

We compare the magnetic field orientation at the scale of the cores with archival data that provide the orientation of the outflows driven by the cores (Nony et al. 2020) and the PA of the major axis of the cores (Motte et al. 2018). The absolute values of these three sets of angle orientation, together with the relative orientations among the magnetic field lines at the location of the cores, the outflow orientations, and the PA of the cores are presented in Table 2. These differences are represented in the form of CDFs in Fig. 4 and compared with Monte Carlo simulations to see whether there is a favored orientation. These simulations select 100 000 pairs of random 3D vectors withan angle difference within a given range. These vectors are then projected onto the plane of the sky to measure their apparent angle difference. We considered the following five ranges of 3D angle differences: 0°–20° (parallel); 20°–50°; 50°–70°; 70°–90° (perpendicular); and 0°–90°(random). To further investigate the distributions’ tendencies, we performed Kolmogorov-Smirnov (KS) tests between the observed and the simulated populations. The KS test is a nonparametric test that quantitatively evaluates the difference between the cumulative distributions of two data sets. We chose the KS test over the chi-squared test as the former is preferable to compare nonequally sampled data. Thus, we computed the p-value (p) that evaluates if our data and a given synthetic population are drawn from the same parent distribution, with low values of p corresponding to different populations.Following a conservative rule of thumb, we rejected the hypothesis that the two populations are drawn from the same distribution when p < 0.05.

thumbnail Fig. 1

Continuum emission at 1.3 mm of the W43-MM1 HMSFR. Contours start at 3σ with steps of3σ, where σ = 2.4 mJy beam−1. Black circles illustrate one third of the primary beam for each of the five pointings, which are further presented in Figs. 2 and 3. The magenta ellipses indicate the cores presented by Motte et al. (2018).

Table 1

Central coordinates in J2000 and noise values for each pointing.

thumbnail Fig. 2

Dust polarization semi-vectors (rotated by 90° to showthe inferred magnetic field orientation) over the continuum emission at 1.3 mm for pointings 1, 2, and 3 toward W43-MM1. Continuum contours levels are at 3, 5, 10, 50, and 100 σI (σI levels in Table 1). The semi vectors in red and blue show the magnetic field orientation where the polarized intensity exceeds a noise level of 3σP and 5σP, respectively(see Table 1). The semi-vectors are plotted every three pixels, which correspond to a Nyquist spatial frequency of four vectors per synthesized beam (two in each dimension). The red and blue cones represent the red-shifted and blue-shifted outflow lobes, respectively. Dashes gray circles represent one third of the ~24′′ primary beam, within which we performed the analysis. The solid green ellipse shows the synthesized beam of 0.55′′ × 0.49′′, PA = −79.4°.

thumbnail Fig. 3

Following Fig. 2: dust polarization semi-vectors (rotated by 90° to show the inferred magnetic field orientation) over the continuum emission at 1.3 mm for pointings 4 and 5. Continuum contours levels are at 3, 5, 10, and 40 σI (σI levels in Table 1).

3.1 Comparison of the magnetic field orientation with that of the core major axis

We compared the PA of the 29 cores falling within the one-third region of the primary beams of our five pointings with the mean PA of the magnetic field inside the cores. The CDF of these differences in angle is presented in Fig. 4a, together with the CDFs arising from the five different synthetic angle simulations. The CDF of the data clearly shows that the orientation of the cores with respect to the magnetic field orientation is only reproduced by the 20–50° population. The results of the KS tests, which compare the observations with the synthetic populations, are reported in Table 3. Here, we confirm that the magnetic field orientation with respect to the cores PA in our data is not random: All of the synthetic populations (including the random distribution), except for the 20–50° distribution,are rejected by the KS test.

3.2 Comparison of the magnetic-field orientation with that of the outflows axis

Nony et al. (2020) studied the outflows associated with the cores of W43-MM1 by observing the 12CO(2-1) and the SiO(5-4) molecular emission lines with an ALMA plus ACA mosaic at an angular resolution of ~0.45′′, which is similar to the present study. In total, they report 46 individual outflow lobes that are associated with 27 cores. On the one hand, about one-fourth of these outflows are monopolar. They could be truly monopolar, as observed in some nearby objects (e.g., HH30, Louvet et al. 2018), or the emission arising from the companion outflow could be absorbed by the surrounding medium. On the other hand, most of the bipolar outflows are not straight, meaning that there is a shift between the PA of the red-shifted lobe with respect to the PA of the blue-shifted lobe (see Figs. 2 and 3). The authors suggest that this could be due to subfragmentation within the cores and/or deflection of the outflows. For these reasons, we consider each outflow lobe independently2. Eighteen of the 27 cores that have one or more outflows are detected in polarization. In Table 2, we report the orientation of the outflows associated with these 18 cores. The orientation of one outflow is defined as the angle of the line linking the core location, which is defined as the continuum peak, and the furthest knot of the outflow (see Nony et al. 2020 for more details about the outflow characteristics).

The CDF of the orientation differences between the magnetic field and the outflows is presented in Fig. 4b. In Table 3 we show the results of the KS tests, which compare the observations with the synthetic populations. Based on these tests, we can reject the hypothesis that the outflows are oriented parallel or perpendicular to the magnetic field. Instead, the KS tests favor either a random orientation of the outflows with respect to the magnetic field lines or the outflows are oriented 50–70° with respectto the magnetic field lines.

Table 2

Position angles of the cores, outflows, and magnetic field in W43.

thumbnail Fig. 4

Cumulative distribution function (CDF, black points) of the angle difference between the magnetic field orientation and the PA of the cores (a), the magnetic field and the outflows (b), and the outflows and the PA of the cores (c) in W43-MM1. The black histogram shows the CDF with 5° bin-widths. The brown bars show the errors on the CDF for each bin. Colored curves are the CDFs resulting from Monte Carlo simulations of different 3D angles projected onto the plane of the sky. The different intervals of angles considered for the synthetic populations are indicated on the bottom right of the plots.

Table 3

Kolmogorov-Smirnov tests for the observed and simulated angle distributions.

3.3 Comparison of the orientation of the outflows with that of the core major axis

The comparison between the PA of the outflows and the PA of the cores is presented in the CDF in Fig. 4c. Based on visual inspection, the observational CDF seems to be consistent with the random synthetic population, suggesting that there is no specific orientation of the outflows with respect to the elongation of the cores. However, when comparing the observations with the synthetic populations through the KS tests, while the random population obtains the highest statistical weight with p ~0.5, the 50–70° population cannot be discarded(p ~0.2).

4 Discussion

Our observations recover an ordered magnetic field pattern in the main protocluster of the W43-MM1 HMSFR (see Fig. 2), plus a few detections in its SW protocluster (see Fig. 3). This is in perfect agreement with previous observations of polarized dust emission toward the W43-MM1 HMSFR (Cortes & Crutcher 2006; Sridharan et al. 2014; Cortes et al. 2016). In the following subsections, we discuss the direction of the magnetic field compared to the orientation of the cores and outflows.

4.1 The alignment of the cores with respect to the magnetic field

Magnetic fields may contribute to support clouds against gravity, thus indirectly affecting the evolution of individual cores. Furthermore, magnetic fields could also directly affect the evolution of cores. In the so-called magnetically regulated core-collapse models (e.g., Galli & Shu 1993a,b; Allen et al. 2003a,b), magnetic fields are dynamically important and dominate the dynamics of the core-collapse by deflecting the infalling gas toward the equatorial plane to form a flattened structure, known as a pseudodisk, around the central protostar. The pseudodisks are not supported by rotation but by magnetic pressure and they develop orthogonally to the magnetic field direction. They have sizes of a few thousands of astronomical units, on the order of the angular scales probed in the present study. Therefore, if magnetically regulated core-collapse were at work, one would expect to observe flattened central regions – the pseudodisks – which are generically referred to as cores in the present study and oriented perpendicular to the magnetic field lines. The very few existing observational studies, which have reported this type of behavior, had very limited samples: Davidson et al. (2011) studied 350 μm polarization observations taken at the CSO toward three low-mass Class 0 protostars, and they report a magnetic field orientation, in loose agreement with the magnetically-controlled core-collapse predictions; Chapman et al. (2013) found the magnetic field orientation to be perpendicular to the pseudodisks in four low-mass cores, using 350 μm polarization observations taken at the CSO; Qiu et al. (2014) report a similar result toward the G240.31+0.07 high-mass core upon using the SMA at 0.88 mm. In contrast, our results, which were obtained with a sample of 29 cores, are not consistent with the prediction of the magnetically-regulated core-collapse models. Instead, as is shown in Fig. 4a, the major axes of the cores are not preferentially oriented orthogonal to the magnetic field, but rather they tend to be oriented 20–50° with respect to the field. Having a nonrandom orientation of the cores with respect to the magnetic field orientation demonstrates that the magnetic field and the dense material are well coupled.

It is interesting to point out that ideal MHD simulations of gravitational collapse of cores have shown that such configurations, where the rotation axis of the core is not parallel to the magnetic field direction, reduce the magnetic braking and favor the formation of circumstellar disks (e.g., Joos et al. 2012). More specifically, recent nonideal MHD simulations have reported that such configurations would engender warped disks, implying that the PA of the pseudo-disks would differ from the PA of the inner disks (see, e.g., the Fig. 3-middle row of Tsukamoto et al. 2018).

4.2 The orientation of the outflows with respect to the magnetic field

One common feature of the two families of models explaining the presence of outflows – the disk-wind models (see, e.g., Blandford & Payne 1982) and the entrainment models (see e.g., Masson & Chernin 1993) – is that they both need the circumstellar disk to be orthogonal to the magnetic field lines, causing the ejection of matter along the magnetic field lines. Hence, if the models are valid, one expects to see alignment between the outflows and the magnetic field near the disks. Carrasco-González et al. (2010) first found an alignment between the outflow-lobe associated to the massive YSO HH 80/81 and the magnetic field in its jet, inferred from polarized synchrotron emission measured with the VLA, up to ~0.5 pc from the driving source. Beuther et al. (2010) studied the high-mass disk-outflow system IRAS 18089-1732 and found that the magnetic field, which was measured from CO(3–2) polarized emission, is aligned with the outflow orientation from small core scales (~7000 au) to larger outflow scales of ~36 000 au. Sridharan et al. (2014) also reported alignment between the magnetic field, which was measured from dust polarized emission, and one outflow in W43-MM1 by using SMA observations at an angular resolution of ~3″. However, this interpretation has since been disproved by a more recent study using ALMA at ~0.5″, which revealed that the outflow seen with the SMA toward W34-MM1 was the sum of 12 outflow lobes with different orientations (Nony et al. 2020). Also, Zhang et al. (2014a) studied a sample of 14 clumps located in different regions where they could compare theorientation of the magnetic field lines with respect to the axes of the outflows. While they refrained from interpreting the detailed structure of the obtained distribution, they report no strong correlation between the outflow and the magnetic field orientations.The lack of alignment between outflows and the magnetic field is also observed in low-mass star-forming regions at the ~1000 au scale (see the recent review by Hull & Zhang 2019). At the envelope scale (~600–1500 au), Galametz et al. (2018) report that the magnetic field is either preferentially observed aligned or perpendicular to the outflow direction when studying a sample of 12 Class 0 envelopes in nearby clouds. They interpret this bimodality by considering that the cases that showed an alignment between the magnetic field and the rotation axis might be the result of a strong magnetic field. They also attribute this bimodality to the fact that most of the cores with a magnetic field perpendicular to the outflows are binaries.

In this work, which is based on a sample of 28 outflow lobes, we find two configurations to be compatible with the observational CDF of the angle difference between the orientation of the outflows and the magnetic field orientation (see Sect. 4.2).

The outflows could be randomly oriented to the magnetic field lines, which is in line with the observations toward low-mass star-forming regions (Hull & Zhang 2019). Assuming that the models for outflow launching are accurate, where gas from the accretion disk should be ejected along magnetic field lines, a random orientation would imply that most of our (high-mass) protostars form in an environment where the magnetic field is too weak to maintain a consistent orientation from the ~2700 au scales that we are probing, down to the 0.1–10 au scales where outflows are launched (Louvet et al. 2018). It would also imply that the orientation of the disks is not controlled by the magnetic field, but by another mechanism provoking angular momentum redistribution, such as interactions in multiple systems or the randomization of the disk orientation by accretion during the early phases of protostellar evolution (Lee et al. 2017).

Alternatively, the outflows could be oriented 50–70° to the magnetic field lines. It is striking to note that if the elongation of the cores corresponds to the major axes of the underlying disks, and considering that outflows propagate orthogonally to their hosting disk, a population of cores that is oriented 20–50° as reported in Sect. 4.1 would translate into a population of outflows oriented 40–70° with respect to the magnetic field lines, which is in excellent agreement with the outflow being preferentially oriented 50–70° to the magnetic field. It is this coherent view, arising from independent correlations, that we might be witnessing here. Nevertheless, we acknowledge that there is a better match for the CDF presenting the 20–50° angle differences between the PA of the cores and the magnetic field orientations than the CFD presenting the 50–70° angle differences between the outflow orientations and the magnetic field. Such a loss of statistical significance is expected if, in some cores, the magnetic field is not strong enough to govern the angular momentum from the core scale down to the outflow ejection location at 0.1–10 au.

4.3 The orientation of the outflows with respect to their hosting cores

The analysispresented in Sect. 3.3 asserts that two configurations are consistent with our data, that is, either the outflows are randomly oriented or oriented 50–70° to their hosting cores. We discuss in turn these two possibilities.

We first comment on the possibility of a random orientation of the outflows with respect to their hosting cores. Recent studies of protostellar objects forming low-mass stars have shown that outflows are usually observed perpendicular to their underlying disks at small scales (<200 au, see e.g., Louvet et al. 2016). Building on the hypothesis that this also happens for high-mass protostars, there are two ways of interpreting this result. If the core elongations are representative of the major axes of underlying disks, the random orientation of the outflowsto the cores would indicate that the outflows get randomly deflected at scales below our angular resolution (~2700 au). Such deflections have been observed in a few outflows in W43-MM1 at a larger scale by Nony et al. (2020) and, therefore, could also occur below our resolution. Nevertheless, it is reasonable to argue that a deflection is more likely to occur with a small angle and that the statistical impact of deflections should be minimal. Therefore, deflections are unlikely to result in a random orientation of the outflows with respect to their hosting cores. If outflows are indeed randomly oriented with respect to the core PA, a more likely explanation is that the elongation of the cores is not representative of the PA of underlying disks. In such a case, the angular momentum of the disks is set by processes acting at scales smaller than the resolution of our observations.

Commenting on the possibility that outflows are oriented 50–70° to their hosting cores, it is important to note that this type of configuration does not match any expectations from the models of launching and/or propagation of outflows. In order to come up with a more physical explanation, we tested bimodal synthetic populations. The fiducial synthetic distributions of angles presented in Sect. 3 only contain single populations, meaning, for example, that in our “50°–70°” distribution all of the angle differences are between 50° and 70° in the 3D space. Instead, the bimodal synthetic populations have a percentage n% of 3D vector pairs with angle differences within a certain angle range and the complementary fraction of vector pairs, 1 − n%, within another angle range. We found a good correlation of the observations with a bimodal population where 85% of the outflows are randomly oriented and 15% are orthogonally oriented to their respective core, which is in line with the assumption from Sect. 4.2. However, given the size of our sample, it is not statistically possible to discriminate between the 50–70° orientation, a random orientation, or an 85% random population plus a 15% orthogonal population.

5 Conclusions

We report thermal dust polarized emission toward a sample of dense cores in the high-mass protostellar cluster W43-MM1, at an angular resolution of ~0.50′′ (~ 0.01 pc, or 2700 au), using ALMA observations in Band 6. We compare the orientation of the magnetic field with archival data of the orientation of 29 dense cores and 28 outflow lobes in this region.

The major axes of the cores are not randomly oriented with respect to the magnetic field, showing that the magnetic field is well coupled to the dense material composing the cores. Instead, the cores are compatible with an orientation of 20–50° with respect to the magnetic field. If confirmed, this result rules out the magnetically-controlled core-collapse models in which a flattened envelope, or pseudodisk, is expected to develop orthogonally to the magnetic field lines. The outflows are oriented 50–70° with respect to the magnetic field (i.e., the orientation of the cores plus 90°) or randomly oriented. Given our statistics, we could not discriminate between these two possibilities.

We propose that, in some cases, the magnetic field at the scale of cores is strong enough to set the orientation of the disk, resulting in an outflow versus magnetic field orientation which is coherent with that of the core versus magnetic field. In other cases, the magnetic field at the scale of cores is not strong enough to set the orientation of the disk. The latter is then controlled by other mechanisms, such as angular momentum redistribution and/or gravitational interaction in multiple systems, which results in a random orientation of the outflow with respect to the core scale magnetic field orientation.

Acknowledgements

We thank the anonymous referee for her/his helpful comments that helped us clarify the article. F.K.L. thanks Maud Galametz and Sylvie Cabrit for fruitful discussions about the collapse of cores and the propagation of outflows, respectively. F.K.L. acknowledges the support of the Fondecyt program N° 3170360. F.M. acknowledges the support of the Joint ALMA Observatory Visitor Program. C.L.H.H. acknowledges the support of both the NAOJ Fellowship as well as JSPS KAKENHI grants 18K13586 and 20K14527. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.01020.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST 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. This work was supported by the Programme National de Physique Stellaire and Physique et Chimie du Milieu Interstellaire (PNPS and PCMI) of CNRS/INSU (with INC/INP/IN2P3) co-funded by CEA and CNES. G.G. and L.B. acknowledge support from CONICYT project AFB 170002.

References

  1. Allen, A., Li, Z.-Y., & Shu, F. H. 2003a, ApJ, 599, 363 [NASA ADS] [CrossRef] [Google Scholar]
  2. Allen, A., Shu, F. H., & Li, Z.-Y. 2003b, ApJ, 599, 351 [NASA ADS] [CrossRef] [Google Scholar]
  3. Andersson, B. G., Lazarian, A., & Vaillancourt, J. E. 2015, ARA&A, 53, 501 [NASA ADS] [CrossRef] [Google Scholar]
  4. Beuther, H., Vlemmings, W. H. T., Rao, R., & vander Tak, F. F. S. 2010, ApJ, 724, L113 [Google Scholar]
  5. Beuther, H., Mottram, J. C., Ahmadi, A., et al. 2018, A&A, 617, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 [NASA ADS] [CrossRef] [Google Scholar]
  7. Camenzind, M. 1990, Rev. Mod. Astron., 3, 234 [Google Scholar]
  8. Carrasco-González, C., Rodríguez, L. F., Anglada, G., et al. 2010, Science, 330, 1209 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  9. Chapman, N. L., Davidson, J. A., Goldsmith, P. F., et al. 2013, ApJ, 770, 151 [NASA ADS] [CrossRef] [Google Scholar]
  10. Ching, T.-C., Lai, S.-P., Zhang, Q., et al. 2017, ApJ, 838, 121 [NASA ADS] [CrossRef] [Google Scholar]
  11. Commerçon, B., Hennebelle, P., & Henning, T. 2011, ApJ, 742, L9 [NASA ADS] [CrossRef] [Google Scholar]
  12. Cortes, P., & Crutcher, R. M. 2006, ApJ, 639, 965 [Google Scholar]
  13. Cortes, P. C., Girart, J. M., Hull, C. L. H., et al. 2016, ApJ, 825, L15 [NASA ADS] [CrossRef] [Google Scholar]
  14. Crutcher, R. M. 2012, ARA&A, 50, 29 [NASA ADS] [CrossRef] [Google Scholar]
  15. Davidson, J. A., Novak, G., Matthews, T. G., et al. 2011, ApJ, 732, 97 [NASA ADS] [CrossRef] [Google Scholar]
  16. Ferreira, J. 1997, A&A, 319, 340 [Google Scholar]
  17. Franco, G. A. P., Alves, F. O., & Girart, J. M. 2010, ApJ, 723, 146 [Google Scholar]
  18. Galametz, M., Maury, A., Girart, J. M., et al. 2018, A&A, 616, A139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Galli, D., & Shu, F. H. 1993a, ApJ, 417, 220 [NASA ADS] [CrossRef] [Google Scholar]
  20. Galli, D., & Shu, F. H. 1993b, ApJ, 417, 243 [NASA ADS] [CrossRef] [Google Scholar]
  21. Girart, J. M., Beltrán, M. T., Zhang, Q., Rao, R., & Estalella, R. 2009, Science, 324, 1408 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  22. Hennebelle, P., Commerçon, B., Joos, M., et al. 2011, A&A, 528, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Hoang, T., & Lazarian, A. 2009, ApJ, 697, 1316 [Google Scholar]
  24. Hull, C. L. H., & Zhang, Q. 2019, Front. Astron. Space Sci., 6, 3 [NASA ADS] [CrossRef] [Google Scholar]
  25. Hull, C. L. H., Plambeck, R. L., Bolatto, A. D., et al. 2013, ApJ, 768, 159 [NASA ADS] [CrossRef] [Google Scholar]
  26. Hull, C. L. H., Plambeck, R. L., Kwon, W., et al. 2014, ApJS, 213, 13 [NASA ADS] [CrossRef] [Google Scholar]
  27. Joos, M., Hennebelle, P., & Ciardi, A. 2012, A&A, 543, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Koch, P. M., Tang, Y.-W., Ho, P. T. P., et al. 2014, ApJ, 797, 99 [NASA ADS] [CrossRef] [Google Scholar]
  29. Konigl, A., & Pudritz, R. E. 2000, Protostars and Planets IV (Tucson, AZ: University of Arizona Press), 759 [Google Scholar]
  30. Lazarian, A. 2007, J. Quant. Spectr. Rad. Transf., 106, 225 [NASA ADS] [CrossRef] [Google Scholar]
  31. Lee, J. W. Y., Hull, C. L. H., & Offner, S. S. R. 2017, ApJ, 834, 201 [NASA ADS] [CrossRef] [Google Scholar]
  32. Louvet, F., Motte, F., Hennebelle, P., et al. 2014, A&A, 570, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Louvet, F., Dougados, C., Cabrit, S., et al. 2016, A&A, 596, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Louvet, F., Dougados, C., Cabrit, S., et al. 2018, A&A, 618, A120 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Masson, C. R., & Chernin, L. M. 1993, ApJ, 414, 230 [NASA ADS] [CrossRef] [Google Scholar]
  36. Matsushita, Y., Sakurai, Y., Hosokawa, T., & Machida, M. N. 2018, MNRAS, 475, 391 [NASA ADS] [CrossRef] [Google Scholar]
  37. McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, ASP Conf. Ser., 376, 127 [Google Scholar]
  38. Motte, F., Nony, T., Louvet, F., et al. 2018, Nat. Astron., 2, 478 [NASA ADS] [CrossRef] [Google Scholar]
  39. Nony, T., Motte, F., Louvet, F., et al. 2020, A&A, 636, A38 [CrossRef] [EDP Sciences] [Google Scholar]
  40. Planck Collaboration Int. XXXIII. (2016), A&A, 586, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Planck Collaboration Int. XXXV. (2016), A&A, 586, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Qiu, K., Zhang, Q., Menten, K. M., et al. 2014, ApJ, 794, L18 [Google Scholar]
  43. Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, 23 [Google Scholar]
  44. Shu, F. H., Najita, J. R., Shang, H., & Li, Z.-Y. 2000, Protostars and Planets IV (Tucson, AZ: University of Arizona Press), 789 [Google Scholar]
  45. Shu, F. H., Li, Z.-Y., & Allen, A. 2004, ApJ, 601, 930 [NASA ADS] [CrossRef] [Google Scholar]
  46. Sridharan, T. K., Rao, R., Qiu, K., et al. 2014, ApJ, 783, L31 [Google Scholar]
  47. Tomisaka, K. 1998, MNRAS, 298, 797 [NASA ADS] [CrossRef] [Google Scholar]
  48. Tsukamoto, Y., Okuzumi, S., Iwasaki, K., Machida, M. N., & Inutsuka, S. 2018, ApJ, 868, 22 [NASA ADS] [CrossRef] [Google Scholar]
  49. Vaillancourt, J. E. 2006, PASP, 118, 1340 [Google Scholar]
  50. Zhang, Q., Qiu, K., Girart, J. M., et al. 2014a, ApJ, 792, 116 [NASA ADS] [CrossRef] [Google Scholar]
  51. Zhang, B., Moscadelli, L., Sato, M., et al. 2014b, ApJ, 781, 89 [NASA ADS] [CrossRef] [Google Scholar]

1

The mosaic was produced with data that were not self-calibrated.

2

We also performed the analysis treating bipolar outflows as single data points and obtained similar results as the ones presented in Sects. 3.2 and 3.3.

All Tables

Table 1

Central coordinates in J2000 and noise values for each pointing.

Table 2

Position angles of the cores, outflows, and magnetic field in W43.

Table 3

Kolmogorov-Smirnov tests for the observed and simulated angle distributions.

All Figures

thumbnail Fig. 1

Continuum emission at 1.3 mm of the W43-MM1 HMSFR. Contours start at 3σ with steps of3σ, where σ = 2.4 mJy beam−1. Black circles illustrate one third of the primary beam for each of the five pointings, which are further presented in Figs. 2 and 3. The magenta ellipses indicate the cores presented by Motte et al. (2018).

In the text
thumbnail Fig. 2

Dust polarization semi-vectors (rotated by 90° to showthe inferred magnetic field orientation) over the continuum emission at 1.3 mm for pointings 1, 2, and 3 toward W43-MM1. Continuum contours levels are at 3, 5, 10, 50, and 100 σI (σI levels in Table 1). The semi vectors in red and blue show the magnetic field orientation where the polarized intensity exceeds a noise level of 3σP and 5σP, respectively(see Table 1). The semi-vectors are plotted every three pixels, which correspond to a Nyquist spatial frequency of four vectors per synthesized beam (two in each dimension). The red and blue cones represent the red-shifted and blue-shifted outflow lobes, respectively. Dashes gray circles represent one third of the ~24′′ primary beam, within which we performed the analysis. The solid green ellipse shows the synthesized beam of 0.55′′ × 0.49′′, PA = −79.4°.

In the text
thumbnail Fig. 3

Following Fig. 2: dust polarization semi-vectors (rotated by 90° to show the inferred magnetic field orientation) over the continuum emission at 1.3 mm for pointings 4 and 5. Continuum contours levels are at 3, 5, 10, and 40 σI (σI levels in Table 1).

In the text
thumbnail Fig. 4

Cumulative distribution function (CDF, black points) of the angle difference between the magnetic field orientation and the PA of the cores (a), the magnetic field and the outflows (b), and the outflows and the PA of the cores (c) in W43-MM1. The black histogram shows the CDF with 5° bin-widths. The brown bars show the errors on the CDF for each bin. Colored curves are the CDFs resulting from Monte Carlo simulations of different 3D angles projected onto the plane of the sky. The different intervals of angles considered for the synthetic populations are indicated on the bottom right of the plots.

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.