Magnetically-regulated fragmentation of a massive, dense and turbulent clump
F. Fontani, B. Commerçon, A. Giannetti, M.T. Beltrán, Á. Sánchez-Monge, L. Testi, J. Brand, P. Caselli, R. Cesaroni, R. Dodson, S. Longmore, M. Rioja, J.C. Tan, C.M. Walmsley
AAstronomy & Astrophysics manuscript no. fragment-16061-aanda-pdf c (cid:13)
ESO 2018October 16, 2018 L etter to the E ditor Magnetically-regulated fragmentation of a massive, dense andturbulent clump
F. Fontani , B. Commerçon , A. Giannetti , M.T. Beltrán , A. Sánchez-Monge , L. Testi , , , J. Brand , P. Caselli , R.Cesaroni , R. Dodson , S. Longmore , M. Rioja , , , J.C. Tan , and C.M. Walmsley INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, I-50125, Florence, Italy Ecole Normale Supérieure de Lyon, CRAL, UMR CNRS 5574, Université Lyon I, 46 Allée d’Italie, 69364, Lyon Cedex 07, France Max-Planck-Institut für Radioastronomie, auf dem Hügel 69, 53121 Bonn, Germany I. Physikalisches Institut, Universität zu Köln, Zülpicher Str. 77, 50937 Köln, Germany European Southern Observatory, Karl-Schwarzschild-Str 2, D-85748, Garching bei München, Germany Gothenburg Center for Advance Studies in Science and Technology, Chalmers University of Technology and University of Gothen-burg, SE-412 96 Gothenburg, Sweden INAF-Istituto di Radioastronomia and Italian ALMA Regional Centre, via P. Gobetti 101, I-40129, Bologna, Italy Max-Planck-Institüt für extraterrestrische Physik, Giessenbachstrasse 1, D-85748, Garching bei München, Germany International Center for Radio Astronomy Research, M468, University of Western Australia, 35, Stirling Hwy, Crawley, WesternAustralia, 6009, Australia Astrophysics Research Institute, Liverpool John Moores University, Liverpool, L3 5RF, UK CSIRO Astronomy and Space Science, 26 Dick Perry Avenue, Kensington WA 6151, Australia Observatorio Astronomico Nacional (IGN), Alfonso XII, 3 y 5, E-28014 Madrid, Spain Departments of Astronomy & Physics, University of Florida, Gainesville, FL 32611, USAReceived date; accepted date
ABSTRACT
Massive stars, multiple stellar systems and clusters are born from the gravitational collapse of massive dense gaseous clumps, and theway these systems form strongly depends on how the parent clump fragments into cores during collapse. Numerical simulations showthat magnetic fields may be the key ingredient in regulating fragmentation. Here we present ALMA observations at ∼ . (cid:48)(cid:48) resolutionof the thermal dust continuum emission at ∼
278 GHz towards a turbulent, dense, and massive clump, IRAS16061–5048c1, in a veryearly evolutionary stage. The ALMA image shows that the clump has fragmented into many cores along a filamentary structure.We find that the number, the total mass and the spatial distribution of the fragments are consistent with fragmentation dominatedby a strong magnetic field. Our observations support the theoretical prediction that the magnetic field plays a dominant role in thefragmentation process of massive turbulent clumps.
Key words.
Stars: formation – ISM: clouds – ISM: molecules – Radio lines: ISM
1. Introduction
High-mass stars, multiple systems, and clusters are born from thegravitational collapse of massive dense clumps (compact struc-tures with M ≥
100 M (cid:12) , and n (H ) ≥ cm − ) inside largemolecular clouds. Stars more massive than 8 M (cid:12) are expected toform either through direct accretion of material in massive coreswithin the clump that does not fragment further (e.g. McKee &Tan 2003, Tan et al. 2013), or as a result of a dynamical evolutionwhere several low-mass seeds competitively accrete matter in ahighly fragmented clump (Bonnell et al. 2004). In the latter sce-nario, each clump forms multiple massive stars and many lowermass stars: the unlucky losers in the competitive accretion com-petition. There is still vigorous debate on which of these scenar-ios is more likely to occur, and fragmentation appears to be par-ticularly important in this debate. Theoretical models and simu-lations show that the number, the mass, and the spatial distribu-tion of the fragments depend strongly on which of the main com-petitors of gravity is dominant. The main physical mechanisms Send o ff print requests to : F. Fontani, e-mail: [email protected] that oppose gravity during collapse are: intrinsic turbulence, ra-diation feedback, and magnetic pressure (e.g. Krumholz 2006,Hennebelle et al. 2011). Feedback from nascent protostellar ob-jects through outflows, winds and / or expansion of ionised re-gions (especially from newly born massive objects) can be im-portant in relatively evolved stages (Bate 2009), but even thenseems to be of only secondary importance (Palau et al. 2013).In a pure gravo-turbulent scenario, the collapsing clumpshould fragment into many cores, the number of which is com-parable to the total mass divided by the Jeans mass (Dobbs etal. 2005); on the other hand, fragmentation can be suppressedby temperature enhancement due to the gravitational energy ra-diated away from the densest portion of the clump that collapsesfirst (Krumholz 2006), or by magnetic support (Hennebelle etal. 2011). The work by Commerçon et al. (2011) has shown thatmodels with strong magnetic support predict fragments moremassive and less numerous than those predicted by the modelswith weak magnetic support. The crucial parameter in their 3-D simulations is µ = ( M / Φ ) / ( M / Φ ) crit , where ( M / Φ ) is the ra-tio between total mass and magnetic flux, and the critical value( M / Φ ) crit , i.e. the ratio for which gravity is balanced by the mag- Article number, page 1 of 7 a r X i v : . [ a s t r o - ph . GA ] A ug & A proofs: manuscript no. fragment-16061-aanda-pdf netic field (thus, for ( M / Φ ) crit > µ .Hiterto, studies of the fragmentation level in massive clumpsat the earliest stages of the gravitational collapse remain lim-ited. This investigation is challenging for several reasons: pris-tine massive clumps are rare and typically located at distanceslarger than 1 kpc, hence to reach the linear resolution requiredfor a consistent comparison with the simulations (about 1000A.U.) requires observations with sub-arcsecond angular resolu-tion. Furthermore, the small mass of the fragments expected inthe simulations (fractions of M (cid:12) ) requires extremely high sen-sitivities. In general, the few studies performed so far with sub-arcsecond angular resolution, or close to ∼ (cid:48)(cid:48) , reveal either lowfragmentation (e.g. Palau et al. 2013, Longmore et al. 2011), ormany fragments but too massive to be consistent with the gravo-turbulent scenario (Bontemps et al. 2010, Zhang et al. 2015).Furthermore, comparisons with models that assume the actualphysical conditions (temperature, turbulence) of the collapsingparent clump have not been published yet.In this letter, we report on the population of fragments de-rived in the image of the dust thermal continuum emission at ∼
278 GHz obtained with the Atacama Large Millimeter Ar-ray (ALMA) towards the source IRAS16061–5048c1, hereafterI16061c1, a massive ( M ∼ (cid:12) , Beltrán et al. 2006, Gian-netti et al. 2013) and dense (column density of H , N ( H ) ∼ . × cm − ) molecular clump located at 3.6 kpc (Fontaniet al. 2005). The clump was detected at 1.2 mm at low angu-lar resolution with the Swedish-ESO Submillimeter Telescope(SEST, panel (A) if Fig.1, Beltrán et al. 2006), and found to benot blended with nearby millimeter clumps, which allows us toavoid confusion in the identification of the fragments. Its largemass and column density make it a potential site for the for-mation of massive stars and rich clusters, according to obser-vational findings (Kau ff mann & Pillai 2010, Lopez-Sepulcre etal. 2010). The clump was classified as an infrared dark cloud be-cause undetected in the images of the Midcourse Space Experi-ment (MSX) infrared satellite, although more sensitive images ofthe Spitzer satellite revealed infrared emission at a wavelength of24 µ m (panel (A) in Fig. 1). Nevertheless, several observationalresults indicate that the possible embedded star formation ac-tivity is in a very early stage (Sanchez-Monge et al. 2013). Inparticular, the depletion factor of CO (ratio between expectedand observed abundance of CO) is 12. This provides strong evi-dence for the chemical youth of the clump, because what causesdepletion factors of CO larger than unity is the freeze-out of thismolecule onto dust grains, a mechanism e ffi cient only in coldand dense pre-stellar and young protostellar cores (Caselli etal. 1999, Emprechtinger et al. 2009, Fontani et al. 2012). Theobservations and the data reduction procedures are presented inSect. 2. Our results are shown in Sect. 3 and discussed in Sect. 4.
2. Observations and data reduction
Observations of I16061c1 with the ALMA array were per-formed during southern winter, 2015. The array was in config-uration C36-6, with maximum baseline of 1091 m. The phasecentre was at R.A. (J2000): 16 h m (cid:48)(cid:48)
61 and Dec (J2000): − ◦ (cid:48) (cid:48)(cid:48) . The total integration time on source was ∼ Fig. 1. (A) : Dust continuum emission map (dashed contours) obtainedwith the SEST telescope with an angular resolution of 24 (cid:48)(cid:48) at 1.2 mmtowards I16061c1 (Beltrán et al. 2006). The map is superimposed onthe Spitzer-MIPS image at 24 µ m (in units of MJy / sr). The circle indi-cates the ALMA primary beam at 278 GHz ( ∼ (cid:48)(cid:48) ). (B) : Enlargementof the rectangular region indicated in panel (A) , showing the contourmap of the thermal dust continuum emission at frequency 278 GHz de-tected with ALMA, in flux density units. The first contour level, andthe spacing between two adjacent contours, both correspond to the 3 σ rms of the image (0.54 mJy / beam). The cross marks the phase center.The ellipse in the bottom left corner shows the synthesized beam, andcorresponds to 0 . × . (cid:48)(cid:48) (Position Angle =
86 deg). The numbersindicate the twelve identified fragments (see Sect. 3). (C) : Simulationsof the thermal dust emission at 278 GHz predicted by the models ofCommerçon et al. (2011), which reproduce the gravitational collapse ofa 300 M (cid:12) clump, in case of strong magnetic support ( µ = (D) : Same as panel (C) for the case µ = (E) : Synthetic ALMA images of the modelspresented in panel (C) . The contours correspond to 0.54, 1.2, 2, 5, 10,30, and 50 mJy / beam. (F) : same as panel (E) for the case µ =
200 (weakmagnetic support). ∼ . (cid:48)(cid:48) ) is ∼ . ∼
30% of the total flux. The remaining ∼
70% is likelycontained in an extended envelope that is resolved out. Contin-uum was extracted by averaging in frequency the line-free chan-
Article number, page 2 of 7ontani et al.: Magnetically regulated fragmentation
Fig. 2.
Top panels : histograms showing the distribution of the peakfluxes ( F peak ) of the fragments identified in the ALMA image ofI16061c1 (black), and in the synthetic images for µ = µ =
200 (blue, right).
Bottom panels : Empirical Cumulative DistributionFunction (ECDF) of the quantities plotted in the top panels. The blackline corresponds to the ALMA data; the green and blue lines indicatethe strong and weak field cases projected on the 3 planes. In all pan-els, The di ff erent line style indicates the projection plane: solid = (x,y);dot-dashed = (x,z); dashed = (y,z). The vertical dotted line correspondsto 0.8 mJy, which is approximately 5 times the rms noise level in boththe real and synthetic maps. Note that the µ = µ =
200 model is strongly biased towards fragments withmasses lower than those observed. nels. The total bandwidth used is ∼ software (McMullinet al. 2007), and the final images were analysed following stan-dard procedures with the software MAPPING of the GILDAS package. The angular resolution of the final image is ∼ . (cid:48)(cid:48) (i.e. ∼
900 AU at the source distance). We were sensitive to point-like fragments of 0.06 M (cid:12) . Together with the continuum, we de-tected several lines among which N H + (3–2). These data will bepresented and discussed in a forthcoming paper. In this letter, weonly use the N H + (3–2) line to derive the virial masses, as wewill show in Sect. 3.
3. Results
B The ALMA map of the dust thermal continuum emission isshown in Fig. 1 (B) : we have detected several dense condensa-tions distributed in a filamentary-like structure extended east-west, surrounded by fainter extended emission. This structurehas been decomposed into twelve fragments (Fig. 1 (B) ). Thefragments have been identified following these criteria: (1) peakintensity greater than 5 times the noise level; (2) two partiallyoverlapping fragments are considered separately if they are sep-arate at their half peak intensity level. The minimum thresholdof 5 times the noise was adopted according to the fact that somepeaks at the edge of the primary beam are comparable to about4-5 times the noise level. We decided to use these criteria and The Common Astronomy Software Applications (CASA) softwarecan be downloaded at http: // casa.nrao.edu http: // / IRAMFR / GILDAS decompose the map into cores by eye instead of using decom-position algorithms (such as Clumpfind) because small changesin their input parameters could lead to big changes in the num-ber of identified clumps (Pineda et al. 2009). The main physicalproperties of the fragments derived from the continuum map, i.e.integrated and peak flux, size, and gas mass, and the methodsused to derive them, are described in Appendix A. The derivedparameters are listed in Table A-1. The mean mass of the frag-ments turns out to be 4.4 M (cid:12) , with a minimum value of 0.7 M (cid:12) and a maximum of ∼ (cid:12) . The diameters (undeconvolved forthe beam) range from 0.011 to 0.032 pc, with a mean value of0.025 pc.To investigate the stability of the fragments, we have cal-culated the virial masses M vir , i.e. the masses required for thecores to be in virial equilibrium, from the line widths observedin N H + (3–2). As stated in Sect. 2, in this work we use thismolecular transition only for the purpose to derive the level ofturbulence (the key ingredient of the models, together with themagnetic support) of the dense gas out of which the fragmentsare formed. The approach used to derive M vir is described in Ap-pendix A and the values obtained are reported in Table A-1. Theaverage ratio between M vir and M computed from the contin-uum emission is about 0.4, indicating that the gravity dominates,according to other ALMA studies of fragmentation (Zhang etal. 2015). However, the uncertainties due to the dust mass opac-ity coe ffi cient (see Eq. 1) can be of a factor 2-3, hence it is dif-ficult to conclude that the fragments are unstable. Moreover, theformula of the virial mass we use does not consider the magneticsupport. Because this latter is expected to be relevant, it is likelythat the fragments are closer to virial equilibrium and would nottend to fragment further. If one assumes, for example, the valueof 0.27 mG measured by Pillai et al. (2015) in another infrared-dark cloud, the ratio between virial mass and gas mass becomesabout 1, and the fragments would be marginally stable. A sim-ilar conclusion is given in Tan et al. (2013), were the dynamicsof four infrared-dark clouds similar to I16061c1 is performed.
4. Discussion and conclusion
We have simulated the gravitational collapse of I16061c1through 3D numerical simulations following Commerçon etal. (2011), adopting mass, temperature, average density, andlevel of turbulence of the parent clump very similar to those mea-sured (Beltrán et al. 2006, Giannetti et al. 2013). In particular,the Mach number setting the initial turbulence, has been derivedfrom the line width of C O (3–2) by Fontani et al. (2012). Be-cause these observations were obtained with angular resolutionof ∼ (cid:48)(cid:48) , and the critical density of the line ( ∼ × cm − ,Fontani et al. 2005) is comparable to the average density of theclump as a whole (Beltrán et al. 2006), the C O line width repre-sents a reasonable estimate of the intrinsic turbulence of the par-ent clump. The models were run for µ = µ =
200 (weakly magnetised case). Then, we have post-processed the simulations data and computed the dust emissionmaps at 278 GHz: the final maps obtained in flux density unitsat the distance of the source have been imaged with the CASAsoftware, in order to reproduce synthetic images with the sameobservational parameters as those of the observations. A detaileddescription of the parameters used for the numerical simulations,of the resulting maps, and how they have been post-processed, isgiven in Appendix B. Further descriptions of the models can befound also in Commerçon et al. (2011). To investigate possiblee ff ects of geometry, we have imaged the outcome of the simu-lations projected on three planes: (x,y), (x,z) and (y,z), where x Article number, page 3 of 7 & A proofs: manuscript no. fragment-16061-aanda-pdf is the direction of the initial magnetic field. As an example, inFigure 1 (C) and 1 (D) we show the results for the cases of strongand weak magnetic support, respectively, projected on the (y,z)plane, i.e. on a direction perpendicular to the magnetic field. Thesynthetic images obtained with CASA are shown in panels (E) and (F) . All the planes for both µ = µ =
200 are shown in Ap-pendix B, in Figs. B-1 and B-2, respectively. An important resultof the simulations (see Figure B-3 in Appendix B) is that thetotal flux seen by the interferometer in the µ = × yrs and then goes through a minimum andstarts gradually to increase. In the µ =
200 case by contrast, thedecrease is not reversed. We conclude from this that in the µ = µ = , while for the µ =
200 case we have analysed the syn-thetic images obtained at the end of the simulations. This corre-sponds to two di ff erent times: t = µ =
2; t = µ = µ =
200 show more fragments withlower peak fluxes, and the overall structure is more chaotic andnever filamentary, independently of the projection plane. Theidentification of the fragments and the derivation of their proper-ties in the synthetic images have been made following the samecriteria and procedures used for the ALMA data. Therefore, anysystematic error introduced by the assumptions made (e.g. theassumed dust temperature, gas-to-dust ratio, dust grain emissiv-ity) are the same both in the real and synthetic images, thus theydo not a ff ect their comparison. The statistical properties of thesynthetic core populations are summarised in Table B-1 of Ap-pendix B. We have also compared the cumulative distribution ofthe peak fluxes of the fragments in the observed and syntheticimages. The results are shown in Fig. 2. The case with µ = µ = µ = µ =
200 model is heav-ily biased below the data. Also, non-parametric statistical test(Anderson-Darling test) implies that all the µ =
200 cases can beexcluded as being drawn from the same parent distribution as theobserved values with a confidence level exceeding 99.8%. The µ = ∼
90% cannot be ex-cluded at the 2 σ level. The deficit of very strong and very weakpeaks in the real image may be due to a di ff erence between the µ values assumed in the simulations and the real one, or to someother unknown (or doubtful) initial assumption such as, e.g., thedensity profile or the homogeneous temperature of the collaps-ing clump.Based on the overall morphologies shown in Fig. 1 (and inFigs. B-1 and B-2), and on the statistical properties of the frag-ments reported in Fig. 2, undoubtedly the model that better re-produces the data is the one with µ =
2. Hence, with these new https: // almascience.eso.org / proposing / technical-handbook ALMA observations, compared with realistic 3D simulationsthat assume as initial conditions the properties of the parentclump, we demonstrate that the fragmentation due to self-gravityis dominated by the magnetic support, based on the evidencethat: (1) the overall morphology of the fragmenting region isfilamentary, and this is predicted only in case of a dominantmagnetic support; (2) the observed fragment mass distributionis most easily understood in simulations assuming substantialmagnetic support.
Acknowledgments.
This paper makes use of the followingALMA data: ADS / JAO.ALMA.2012.1.00366.S. ALMA is apartnership of ESO (representing its member states), NSF (USA)and NINS (Japan), together with NRC (Canada), NSC andASIAA (Taiwan), and KASI (Republic of Korea), in cooper-ation with the Republic of Chile. The Joint ALMA Observa-tory is operated by ESO, AUI / NRAO and NAOJ. We acknowl-edge the Italian-ARC node for their help in the reduction ofthe data. We acknowledge partial support from Italian Ministerodell’Istruzione, Universitá e Ricerca through the grant ProgettiPremiali 2012 − iALMA (CUP C52I13000140001) and fromGothenburg Centre of Advanced Studies in Science and Tech-nology through the program Origins of habitable planets . References
Bate, M.R. 2009, MNRAS, 392, 1363Beltrán, M.T., Brand, J., Cesaroni, R., Fontani, F., Pezzuto, S. et al. 2006, A&A,447, 221Beuther, H., Schilke, P., Sridharan, T.K., Menten, K.M., Walmsley, C.M.,Wyrowski, F. 2002, A&A, 566, 945Bleuler, A., Teyssier, R. 2014, MNRAS, 445, 4015Bonnell, I.A., Vine, S.G., Bate M.R. 2004, MNRAS, 349, 735Bontemps, S., Motte, F., Csengeri, T., Schneider, N. 2010, A&A, 524, 18Caselli, P., Walmsley, C.M., Tafalla, M., Dore, L., Myers, P. 1999, ApJ, 523,L165Commerçon, B., Hennebelle, P., Henning, T. 2011, ApJ, 742, L9Commerçon, B., Launhardt, R., Dullemond, C., Henning, Th. 2012, A&A, 545,98Commerçon, B., Debout, V., Teyssier, R. 2014, A&A, 563, 11Crutcher, R.M. 2013, ARA&A, 50, 29Dobbs, C.L., Bonnell, I.A., Clark, P.C. 2005, MNRAS, 360, 2Emprechtinger, M., Caselli, P., Volgenau, N.H., Stutzki, J., Wiedner, M.C. 2009,A&A, 493, 89Fontani, F., Beltrán, M.T., Brand, J., Cesaroni, R., Molinari, S. et al. 2005, A&A,432, 921Fontani, F., Giannetti, A., Beltrán, M.T., Dodson, R., Rioja, M. et al. 2012, MN-RAS, 423, 2342Fromang, S., Hennebelle, P., Teyssier, R. 2006, A&A, 457, 371Giannetti, A., Brand, J., Sanchez-Monge, Á., Fontani, F., Cesaroni, R. et al. 2013,A&A, 556, 16Hennebelle, P., Commerçon, B., Joos, M., Klessen, R.S., Krumholz, M., Tan,J.C., Teyssier, R. 2011, A&A, 528, 72Kau ff mann, J. & Pillai, T. 2010, ApJ, 723, L7Krumholz, M.R. 2006, ApJ, 641, L45Krumholz, M.R., Klein, R., McKee, C.F., O ff ner, S.S.R., Cunningham,A.J.A. 2009, Science, 323, 754Longmore, S.N., Pillai, T., Keto, E., Zhang, Q., Qiu, K. 2011, ApJ, 726, L97Lopez-Sepulcre, A., Cesaroni, R., Walmsley, C.M. 2010, A&A, 517, 66Maret, S., Hily-Blant, P., Pety, J., Bardeau, S., Reynier, E. 2011, A&A, 526, 47McKee, C. & Tan, J.C. 2003, ApJ, 585, 850McMullin, J.P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in Astro-nomical Data Analysis Software and Systems XVI, eds. R. A. Shaw, F. Hill,& D. J. Bell, ASP Conf. Ser., 376, 127Mouschovias, T.C. & Spitzer, L.J. ApJ, 210, 326Ossenkopf, V., & Henning, T. 1994, A&A, 291, 943Palau, A., Fuente, A., Girart, J.M., Estalella, R., Ho, P.T.P., Sánchez-Monge, A.,et al. 2013, ApJ, 762, 120Pillai, T., Kau ff mann, J., Tan, J.C., Goldsmith, P.F., Carey, S.J., Menten,K.M. 2015, ApJ, 799, 74Pineda, J.E., Rosolowsky, E.W., Goodman, A.A. 2009, ApJ, 699, L134Sanchez-Monge, Á, Beltrán, M.T., Cesaroni, R., Fontani, F., Brand, et al. 2013,A&A, 550, 21Tan, J.C., Kong, S., Butler, M.J., Caselli, P., Fontani, F. 2013, ApJ, 779, 96Teyssier, R. 2002, A&A, 385, 337Weiss, A., Hillebrandt, W., Thomas, H.-C., Ritter, H. 2005, Book review: Coxand Giuli’s Principles of Stellar Structure, Cambridge, UK: Princeton Pub-lishing Associates LtdZhang, Q., Wang, K., Xing, L., Jiménez-Serra, I. 2015, ApJ, 804, 141 Article number, page 4 of 7ontani et al.: Magnetically regulated fragmentation
Appendix A: Physical properties of the fragments
A-1. Derivation of the parameters – Integrated flux densities : The integrated flux densities ofthe fragments, S ν , have been obtained from the 3 σ rms con-tour in the continuum image. In the few cases for which the3 σ rms contours of two adjacent fragments are partly over-lapping (e.g. fragments 5, 6, and 7 in Fig. 1), the edges be-tween the two have been defined by eye at approximatelyhalf of the separation between the peaks. The results areshown in Table A-1. – gas masses : The gas mass of each fragment has been calcu-lated from the equation: M = g S ν d κ ν B ν ( T d ) , (1)where S ν is the integrated flux density, d is the distance to thesource, κ ν is the dust mass opacity coe ffi cient, g is the gas-to-dust ratio (assumed to be 100), and B ν ( T d ) is the Planckfunction for a black body of temperature T d . We adopted T d =
25 K, corresponding to the gas temperature derivedby Giannetti et al. (2013), assuming coupling between gasand dust (reasonable assumption at the high average densityof the clump). The dust mass opacity coe ffi cient was derivedfrom the equation κ ν = κ ν ( ν/ν ) β . We assumed β = κ ν = .
899 cm − gr − at ν =
230 GHz, according to Os-senkopf & Henning (1994). The largest mass derived is 9M (cid:12) , the smallest is ∼ . (cid:12) (see Table A-1). – Size : the size of each fragment has been estimated as thediameter of the circle with area equivalent to that encom-passed by the 3 σ rms contour level. The results are shown inTable A-1. The ALMA beam size is much smaller than thesize of the fragments, so that deconvolution for the beam isirrelevant to derive the source size. – Virial masses : the virial masses were derived in this way:first, we extracted the N H + (3–2) spectra from the 3 σ rmslevel of the 12 continuum cores. All the spectra are fitted inan automatic way using a procedure based on the integra-tion of the python module PyMC and the CLASS extensionWEEDS (Maret et al. maret). Then, the virial masses, M vir ,are computed from the formula M vir = r ∆ v M (cid:12) , (2)where r is the size of the fragment, and ∆ v is the line width athalf maximum of the average N H + (3–2) spectrum obtainedfrom the fitting procedure described above. The results areshown in Table A-1. Appendix B: Simulations and synthetic images
B-1. Methods and initial conditions for the numericalcalculations
We perform a set of two radiation-magneto-hydrodynamics cal-culations which includes the radiative feedback from the accret-ing protostars. We use the RAMSES code with the grey flux-limited-di ff usion approximation for radiative transfer and theideal MHD for magnetic fields (Commerçon et al. 2012, 2014,Teyssier 2002, Fromang et al. 2006). The initial conditions aresimilar to those used in Hennebelle et al. (2011) and Com-merçon et al. (2011) with slight modifications in order to matchroughly the observed properties of I16061c1. Note that the mod-els presented in this paper have been made with initial condi-tions very similar to those measured from previous observations in I16061c1, to perform an appropriate comparison with obser-vations for this specific source. Our aim is not to fine-tune theinitial conditions such that the models best reproduce the obser-vations. We consider an isolated spherical core of mass 300 M (cid:12) ,radius 0.25 pc and temperature 20 K. We assume a Plummer-likeinitial density profile ρ ( r ) = ρ c / (1 + ( r / r ) ), with ρ c = . × cm − and r = .
085 pc, and a factor 10 for the density con-trast between the center and the border of the core. Such a den-sity profile is suggested by observational findings (Beuther etal. 2002). The initial magnetic field is aligned with the x-axisand its intensity is proportional to the column density throughthe cloud (Hennebelle et al. 2011). In this paper, we investigatetwo degrees of magnetization, µ = µ = O observations (Fontaniet al. 2012). Following Hennebelle et al. (2011, Eq. 2 therein),the virial parameter is α vir ∼ .
72 for µ =
2, and α vir ∼ .
54 for µ =
200 (close to virial equilibrium in both cases). The two cal-culations, µ = µ = ff ect ofdi ff erent initial turbulent seeds is beyond the scope of this paper.The computational box has a 2563 resolution, and the grid is re-fined according to the local Jeans length (at least 8 cells / Jeanslength) down to 7 levels of refinement (minimum cell size of 13AU). Below 13 AU, the collapsing gas is described using sub-grids models attached to sink particles, similar to what is donein other studies (Krumholz et al. 2009). We use the sink parti-cle method presented in (Bleuler et al. 2014), though with slightmodifications on the checks performed for the sink creation. Thesink particles accrete the gas that sit in their accretion volumes(sphere of radius ∼
52 AU, 4 cells) and that is Jeans unstable. Weconsider that half of the mass accreted into the sink particles ac-tually goes into stellar material. The luminosity of the protostarsis then computed using mass-radius and mass-luminosity empir-ical relations of main sequence stars (e.g. Weiss et al. 2005) andis injected within the accretion volume in the computational do-main as a source term (e.g. Krumholz et al. 2009). We do notaccount for accretion luminosity.
B-2. Outcomes of the numerical calculations
We run the calculations until they reach a star formation e ffi -ciency (SFE) >
20% (where the star formation e ffi ciency corre-sponds to the ratio between the mass of the gas accreted into thesink particles and the total mass of the cloud). Again, the choiceof the times at which we stop the calculations is not aimed at bestreproducing the observed values. Model µ = = µ =
200 at time t = µ = (cid:12) ) while model µ =
200 has formed119 sink particles (for a total mass of 85 M (cid:12) ). Fig. B-3 shows thetime evolution after the first sink creation of: the SFE, the num-
Article number, page 5 of 7 & A proofs: manuscript no. fragment-16061-aanda-pdf
Table A-1.
Peak position (in R.A. and Dec. J2000), integrated flux (inside the 3 σ rms contour level), peak flux, diameter, mass, line width at halfmaximum, and virial mass of the 12 fragments identified in Fig. 1 (B) . The line widths are computed from the N H + (3–2) spectra extracted fromthe polygons defining the external profile of each fragment, as explained in Appendix A. Fragment R.A. (J2000) Dec. (J2000) S ν S peak ν D M ∆ v M vir h:m:s deg : (cid:48) : (cid:48)(cid:48) Jy Jy beam − pc M (cid:12) km s − M (cid:12) − − − − − − − − − − − − × µ =
2, which probablyreflects the temperature increase caused by the radiation of thesink particles.
B-3. Production of the synthetic images
We first postprocessed the RAMSES calculations results usingthe radiative transfer code RADMC-3D and the interface pre-sented in (Commerçon et al. 2012). We produced dust emissionmaps at 278 GHz (see Figs. B-1 and B-2). We do not accountfor the stellar luminosities in the synthetic images since the stel-lar irradiation is reprocessed in the envelope at millimetre wave-lengths. However, we attempted to create models accounting forprotostellar luminosities, and found results that do not changesignificantly at the wavelength considered. The ALMA syntheticimages of the numerical simulations have been then producedthrough the CASA software: first, synthetic visibilities have beencreated with the task simobserve , which have then been imagedwith the task simanalyze . To precisely reproduce the observa-tions, in the tasks we have used the same parameters of the ob-servations: integration time on source of 18 minutes, precipitablewater vapour of about 1.8 mm, array configuration C36-6, starthour angle of 2.4 hours (see Sect. 2). The population of frag-ments in the final synthetic images were derived following thesame procedure described in Sect. 4. http: // / dullemond / software / radmc-3d / Fig. B-1.
Top panels show the thermal dust continuum emission mapat frequency 278 GHz predicted by the models of Commerçon etal. (2011), which reproduce the gravitational collapse of a 300 M (cid:12) clump, in case of strong magnetic support ( µ =
2) at time t = / beam. Contour levels are0.6, 1, 5, 10, 30 and 50 mJy beam − in all bottom panels. Table B-1.
Statistical comparison between the fragment population de-rived from the ALMA image of I16061c1 shown in Fig. 1 and the sim-ulations presented in Figs. B1 and B2 of the Appendix. The derivationof the parameters obtained for both the observed and synthetic imagesis described in Sect. 3 and in Appendix A, respectively. S tot ν M tot N D mean S mean ν M mean Jy M (cid:12) pc Jy M (cid:12)
ALMA 0.52 53 12 0.025 0.042 4.42 µ = µ = µ = µ =
200 (x,y) 0.22 23 13 0.015 0.017 1.74 µ =
200 (x,z) 0.24 25 15 0.014 0.016 1.67 µ =
200 (y,z) 0.28 24 16 0.016 0.021 2.19Article number, page 6 of 7ontani et al.: Magnetically regulated fragmentation
Fig. B-2.
Same as Fig. B-1 for the case µ =
200 at time t = N u m b e r o f s i n k µ =2µ =200 S F E ( % ) Time (kyr) T o t a l f l u x ( J y ) (x,y)(x,z)(y,z) Fig. B-3.
From top to bottom: evolution with time of the number ofsink particles, of the SFE, and of the total flux emission at 278 GHz(within a total area of 80000 AU × ff erent linescorrespond to the di ffff