Within ten nearby d lt 450 pc Gould belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz and the gas column density structures quantified by the gradient of the column density NSUBHSUB. The selected regions covering several degrees in size are analysed at an effective angular resolution of 10 FWHM thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from NSUBHSUB 10SUP21SUP to10SUP23SUP cmSUP2SUP and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analysed in bins of column density using the novel statistical tool called histogram of relative orientations. Throughout this study we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region I.e. we assume a constant dust grain alignment efficiency independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing NSUBHSUB from mostly parallel or having no preferred orientation to mostly perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfvnic or subAlfvnic turbulence implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck. We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.
['ISM: general', 'ISM: magnetic fields', 'ISM: clouds', 'dust', 'extinction', 'submillimeter: ISM', 'infrared: ISM', 'Astrophysics - Astrophysics of Galaxies']
Planck intermediate results. XXXV. Probing the role of the magnetic field in the formation of structure in molecular clouds
{'Planck intermediate results. XXXV. Probing the role of the magnetic field in the formation of structure in molecular clouds': "Planck Collaboration: P. A. R. Ade 83 , N. Aghanim 56 , M. I. R. Alves 56 , M. Arnaud 70 , D. Arzoumanian 56 , M. Ashdown 66 ; 5 , J. Aumont 56 , C. Baccigalupi 81 , A. J. Banday 89 ; 8 , R. B. Barreiro 61 , N. Bartolo 27 ; 62 , E. Battaner 90 ; 91 , K. Benabed 57 ; 88 , A. Benoˆıt 54 , A. Benoit-L'evy 21 ; 57 ; 88 , J.-P. Bernard 89 ; 8 , M. Bersanelli 30 ; 46 , P. Bielewicz 78 ; 8 ; 81 , J. J. Bock 63 ; 9 , L. Bonavera 61 , J. R. Bond 7 , J. Borrill 11 ; 85 , F. R. Bouchet 57 ; 84 , F. Boulanger 56 , A. Bracco 56 , C. Burigana 45 ; 28 ; 47 , E. Calabrese 87 , J.-F. Cardoso 71 ; 1 ; 57 , A. Catalano 72 ; 69 , H. C. Chiang 24 ; 6 , P. R. Christensen 79 ; 33 , L. P. L. Colombo 20 ; 63 , C. Combet 72 , F. Couchot 68 , B. P. Crill 63 ; 9 , A. Curto 61 ; 5 ; 66 , F. Cuttaia 45 , L. Danese 81 , R. D. Davies 64 , R. J. Davis 64 , P. de Bernardis 29 , A. de Rosa 45 , G. de Zotti 42 ; 81 , J. Delabrouille 1 , C. Dickinson 64 , J. M. Diego 61 , H. Dole 56 ; 55 , S. Donzelli 46 , O. Dor'e 63 ; 9 , M. Douspis 56 , A. Ducout 57 ; 52 , X. Dupac 35 , G. Efstathiou 58 , F. Elsner 21 ; 57 ; 88 , T. A. Enßlin 76 , H. K. Eriksen 59 , D. Falceta-Gonc¸alves 82 ; 34 , E. Falgarone 69 , K. Ferri'ere 89 ; 8 , F. Finelli 45 ; 47 , O. Forni 89 ; 8 , M. Frailis 44 , A. A. Fraisse 24 , E. Franceschi 45 , A. Frejsel 79 , S. Galeotta 44 , S. Galli 65 , K. Ganga 1 , T. Ghosh 56 , M. Giard 89 ; 8 , E. Gjerløw 59 , J. Gonz'alez-Nuevo 16 ; 61 , K. M. G'orski 63 ; 92 , A. Gregorio 31 ; 44 ; 50 , A. Gruppuso 45 , J. E. Gudmundsson 24 , V. Guillet 56 , D. L. Harrison 58 ; 66 , G. Helou 9 , P. Hennebelle 70 , S. Henrot-Versill'e 68 , C. Hern'andez-Monteagudo 10 ; 76 , D. Herranz 61 , S. R. Hildebrandt 63 ; 9 , E. Hivon 57 ; 88 , W. A. Holmes 63 , A. Hornstrup 13 , K. M. Hu GLYPH<11> enberger 22 , G. Hurier 56 , A. H. Ja GLYPH<11> e 52 , T. R. Ja GLYPH<11> e 89 ; 8 , W. C. Jones 24 , M. Juvela 23 , E. Keihanen 23 , R. Keskitalo 11 , T. S. Kisner 74 , J. Knoche 76 , M. Kunz 14 ; 56 ; 2 , H. Kurki-Suonio 23 ; 40 , G. Lagache 4 ; 56 , J.-M. Lamarre 69 , A. Lasenby 5 ; 66 , M. Lattanzi 28 , C. R. Lawrence 63 , R. Leonardi 35 , F. Levrier 69 , M. Liguori 27 ; 62 , P. B. Lilje 59 , M. Linden-Vørnle 13 , M. L'opez-Caniego 35 ; 61 , P. M. Lubin 25 , J. F. Mac'ıas-P'erez 72 , D. Maino 30 ; 46 , N. Mandolesi 45 ; 28 , A. Mangilli 56 ; 68 , M. Maris 44 , P. G. Martin 7 , E. Mart'ınez-Gonz'alez 61 , S. Masi 29 , S. Matarrese 27 ; 62 ; 38 , A. Melchiorri 29 ; 48 , L. Mendes 35 , A. Mennella 30 ; 46 , M. Migliaccio 58 ; 66 , M.-A. Miville-Deschˆenes 56 ; 7 , A. Moneti 57 , L. Montier 89 ; 8 , G. Morgante 45 , D. Mortlock 52 , D. Munshi 83 , J. A. Murphy 77 , P. Naselsky 79 ; 33 , F. Nati 24 , C. B. Netterfield 17 , F. Noviello 64 , D. Novikov 75 , I. Novikov 79 ; 75 , N. Oppermann 7 , C. A. Oxborrow 13 , L. Pagano 29 ; 48 , F. Pajot 56 , R. Paladini 53 , D. Paoletti 45 ; 47 , F. Pasian 44 , L. Perotto 72 , V. Pettorino 39 , F. Piacentini 29 , M. Piat 1 , E. Pierpaoli 20 , D. Pietrobon 63 , S. Plaszczynski 68 , E. Pointecouteau 89 ; 8 , G. Polenta 3 ; 43 , N. Ponthieu 56 ; 51 , G. W. Pratt 70 , S. Prunet 57 ; 88 , J.-L. Puget 56 , J. P. Rachen 18 ; 76 , M. Reinecke 76 , M. Remazeilles 64 ; 56 ; 1 , C. Renault 72 , A. Renzi 32 ; 49 , I. Ristorcelli 89 ; 8 , G. Rocha 63 ; 9 , M. Rossetti 30 ; 46 , G. Roudier 1 ; 69 ; 63 , J. A. Rubi˜no-Mart'ın 60 ; 15 , B. Rusholme 53 , M. Sandri 45 , D. Santos 72 , M. Savelainen 23 ; 40 , G. Savini 80 , D. Scott 19 , J. D. Soler 56 ? , V. Stolyarov 5 ; 86 ; 67 , R. Sudiwala 83 , D. Sutton 58 ; 66 , A.-S. Suur-Uski 23 ; 40 , J.-F. Sygnet 57 , J. A. Tauber 36 , L. Terenzi 37 ; 45 , L. To GLYPH<11> olatti 16 ; 61 ; 45 , M. Tomasi 30 ; 46 , M. Tristram 68 , M. Tucci 14 , G. Umana 41 , L. Valenziano 45 , J. Valiviita 23 ; 40 , B. Van Tent 73 , P. Vielva 61 , F. Villa 45 , L. A. Wade 63 , B. D. Wandelt 57 ; 88 ; 26 , I. K. Wehus 63 , N. Ysard 23 , D. Yvon 12 , and A. Zonca 25  \n(A GLYPH<14> liations can be found after the references)  \nReceived 13 February 2015 / Accepted 25 May 2015", 'ABSTRACT': "Within ten nearby ( d < 450 pc) Gould Belt molecular clouds we evaluate statistically the relative orientation between the magnetic field projected on the plane of sky, inferred from the polarized thermal emission of Galactic dust observed by Planck at 353 GHz, and the gas column density structures, quantified by the gradient of the column density, N H. The selected regions, covering several degrees in size, are analysed at an e GLYPH<11> ective angular resolution of 10 0 FWHM,thus sampling physical scales from 0.4 to 40 pc in the nearest cloud. The column densities in the selected regions range from N H GLYPH<25> 10 21 to 10 23 cm GLYPH<0> 2 , and hence they correspond to the bulk of the molecular clouds. The relative orientation is evaluated pixel by pixel and analysed in bins of column density using the novel statistical tool called 'Histogram of Relative Orientations'. Throughout this study, we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region, i.e., we assume a constant dust grain alignment e GLYPH<14> ciency, independent of the local environment. Within most clouds we find that the relative orientation changes progressively with increasing N H, from mostly parallel or having no preferred orientation to mostly perpendicular. In simulations of magnetohydrodynamic turbulence in molecular clouds this trend in relative orientation is a signature of Alfv'enic or sub-Alfv'enic turbulence, implying that the magnetic field is significant for the gas dynamics at the scales probed by Planck . We compare the deduced magnetic field strength with estimates we obtain from other methods and discuss the implications of the Planck observations for the general picture of molecular cloud formation and evolution.  \nKey words. ISM: general, dust, magnetic fields, clouds - Infrared: ISM - Submillimetre: ISM", '1. Introduction': "The formation and evolution of molecular clouds (MCs) and their substructures, from filaments to cores and eventually to stars, is the product of the interaction between turbulence, magnetic fields, and gravity (Bergin & Tafalla 2007; McKee & Ostriker 2007). The study of the relative importance of these dynamical processes is limited by the observational techniques used to evaluate them. These limitations have been particularly critical when integrating magnetic fields into the general picture of MC dynamics (Elmegreen & Scalo 2004; Crutcher 2012; Heiles & Haverkorn 2012; Hennebelle & Falgarone 2012).  \nThere are two primary methods of measuring magnetic fields in the dense interstellar medium (ISM). First, observation of the Zeeman e GLYPH<11> ect in molecular lines provides the line-of-sight component of the field B k (Crutcher 2005). Second, polarization maps - in extinction from background stars and emission from dust - reveal the orientation of the field averaged along the line of sight and projected on the plane of the sky (Hiltner 1949; Davis & Greenstein 1951; Hildebrand 1988; Planck Collaboration Int. XXI 2014).  \nAnalysis of the Zeeman e GLYPH<11> ect observations presented by Crutcher et al. (2010) shows that in the di GLYPH<11> use ISM sampled by H i lines ( n H < 300 cm GLYPH<0> 3 ), the maximum magnetic field strength B max does not scale with density. This is interpreted as the e GLYPH<11> ect of di GLYPH<11> use clouds assembled by flows along magnetic field lines, which would increase the density but not the magnetic field strength. In the denser regions ( n H > 300 cm GLYPH<0> 3 ), probed by OH and CN spectral lines, the same study reports a scaling of the maximum magnetic field strength B max / n 0 : 65 H . The latter observation can be interpreted as the e GLYPH<11> ect of isotropic contraction of gas too weakly magnetized for the magnetic field to a GLYPH<11> ect the morphology of the collapse. However, given that the observations are restricted to pencil-like lines of sight and the molecular tracers are not homogeneously distributed, the Zeeman e GLYPH<11> ect measurements alone are not su GLYPH<14> cient to determine the relative importance of the magnetic field at the multiple scales within MCs.  \nThe observation of starlight polarization provides an estimate of the projected magnetic field orientation in particular lines of sight. Starlight polarization observations show coherent magnetic fields around density structures in MCs (Pereyra & Magalh˜aes 2004; Franco et al. 2010; Sugitani et al. 2011; Chapman et al. 2011; Santos et al. 2014). The coherent polarization morphology can be interpreted as the result of dynamically important magnetic fields. However, these observations alone are not su GLYPH<14> cient to map even the projected magnetic field morphology fully and in particular do not tightly constrain the role of magnetic fields in the formation of structure inside MCs.  \nThe study of magnetic field orientation within the MCs is possible through the observation of polarized thermal emission from dust. Far-infrared and submillimetre polarimetric observations have been limited to small regions up to hundreds of square arcminutes within clouds (Li et al. 2006; Matthews et al. 2014) or to large sections of the Galactic plane at a resolution of several degrees (Benoˆıt et al. 2004; Bierman et al. 2011). On the scale of prestellar cores and cloud segments, these observations reveal both significant levels of polarized emission and coherent field morphologies (Ward-Thompson et al. 2000; Dotson et al. 2000; Matthews et al. 2009).  \nThe strength of the magnetic field projected on the plane of the sky ( B ? ) can be estimated from polarization maps using the Davis-Chandrasekhar-Fermi (DCF) method (Davis 1951; Chandrasekhar & Fermi 1953). As discussed in Appendix D, it is  \nassumed that the dispersion in polarization angle & 1 is entirely due to incompressible and isotropic turbulence. Turbulence also a GLYPH<11> ects the motion of the gas and so broadens profiles of emission and absorption lines, as quantified by dispersion GLYPH<27> v k . In the DCF interpretation B ? is proportional to the ratio GLYPH<27> v k =& . Application of the DCF method to subregions of the Taurus MC gives estimates of B ? GLYPH<25> 10 GLYPH<22> Gin low-density regions and GLYPH<25> 25 to GLYPH<25> 42 GLYPH<22> G inside filamentary structures (Chapman et al. 2011). Values of B ? GLYPH<25> 760 GLYPH<22> G have been found in dense parts of the Orion MC region (Houde et al. 2009). Because of the experimental di GLYPH<14> culties involved in producing large polarization maps, a complete statistical study of the magnetic field variation across multiple scales is not yet available.  \nAdditional information on the e GLYPH<11> ects of the magnetic field on the cloud structure is found by studying the magnetic field orientation inferred from polarization observations relative to the orientation of the column density structures. Patterns of relative orientation have been described qualitatively in simulations of magnetohydrodynamic (MHD) turbulence with di GLYPH<11> erent degrees of magnetization. This is quantified as half the ratio of the gas pressure to the mean-field magnetic pressure (Ostriker et al. 2001; Heitsch et al. 2001), with the resulting turbulence ranging from sub-Alfv'enic to super-Alfv'enic. Quantitative analysis of simulation cubes, where the orientation of B is available directly, reveals a preferred orientation relative to density structures that depends on the initial magnetization of the cloud (Hennebelle 2013; Soler et al. 2013). Using simple models of dust grain alignment and polarization e GLYPH<14> ciency to produce synthetic observations of the simulations, Soler et al. (2013) showed that the preferred relative orientation and its systematic dependence on the degree of magnetization are preserved.  \nObservational studies of relative orientation have mostly relied on visual inspection of polarization maps ( e.g., Myers & Goodman 1991; Dotson 1996). This is adequate for evaluating general trends in the orientation of the field. However, it is limited ultimately by the need to represent the field orientation with pseudo-vectors, because when a large polarization map is to be overlaid on a scalar-field map, such as intensity or column density, only a selection of pseudo-vectors can be plotted. On the one hand, if the plotted pseudo-vectors are the result of averaging the Stokes parameters over a region, then the combined visualization illustrates di GLYPH<11> erent scales in the polarization and in the scalar field. On the other hand, if the plotted pseudo-vectors correspond to the polarization in a particular pixel, then the illustrated pattern is influenced by small-scale fluctuations that might not be significant in evaluating any trend in relative orientation.  \nTassis et al. (2009) present a statistical study of relative orientation between structures in the intensity and the inferred magnetic field from polarization measured at 350 GLYPH<22> m towards 32 Galactic clouds in maps of a few arcminutes in size. Comparing the mean direction of the field to the semi-major axis of each cloud, they find that the field is mostly perpendicular to that axis. Similarly, Li et al. (2013) compared the relative orientation in 13 clouds in the Gould Belt, calculating the main cloud orientation from the extinction map and the mean orientation of the intercloud magnetic field from starlight polarization. That study reported a bimodal distribution of relative cloud and field orientations; that is, some MCs are oriented perpendicular and some parallel to the mean orientation of the intercloud field. In both studies each cloud constitutes one independent observation  \nFig. 1. Magnetic field and column density measured by Planck towards the Taurus MC. The colours represent column density. The 'drapery' pattern, produced using the line integral convolution method (LIC, Cabral & Leedom 1993), indicates the orientation of magnetic field lines, orthogonal to the orientation of the submillimetre polarization.  \n<!-- image -->  \nof relative orientation, so that the statistical significance of each study depends on the total number of clouds observed. In a few regions of smaller scales, roughly a few tenths of a parsec, Koch et al. (2013) report a preferred orientation of the magnetic field, inferred from polarized dust emission, parallel to the gradient of the emission intensity.  \nBy measuring the intensity and polarization of thermal emission from Galactic dust over the whole sky and down to scales that probe the interiors of nearby MCs, Planck 2 provides an unprecedented data set from a single instrument and with a common calibration scheme, for studying the morphology of the magnetic field in MCs and the surrounding ISM, as illustrated for the Taurus region in Fig. 1. We present a quantitative analysis of the relative orientation in a set of nearby ( d < 450 pc) well-known MCs to quantify the role of the magnetic field in the formation of density structures on physical scales ranging from tens of parsecs to approximately one parsec in the nearest clouds.  \nThe present work is an extension of previous findings, as reported by the Planck collaboration, on their study of the polarized thermal emission from Galactic dust. Previous studies in-  \nclude an overview of this emission (Planck Collaboration Int. XIX 2014), which reported dust polarization percentages up to 20 % at low N H, decreasing systematically with increasing N H to a low plateau for regions with N H > 10 22 cm GLYPH<0> 2 . Planck Collaboration Int. XX (2014) presented a comparison of the polarized thermal emission from Galactic dust with results from simulations of MHD turbulence, focusing on the statistics of the polarization fractions and angles. Synthetic observations were made of the simulations under the simple assumption of homogeneous dust grain alignment e GLYPH<14> ciency. Both studies reported that the largest polarization fractions are reached in the most diffuse regions. Additionally, there is an anti-correlation between the polarization percentage and the dispersion of the polarization angle. This anti-correlation is reproduced well by the synthetic observations, indicating that it is essentially caused by the turbulent structure of the magnetic field.  \nOver most of the sky Planck Collaboration Int. XXXII (2014) analysed the relative orientation between density structures, which is characterized by the Hessian matrix, and polarization, revealing that most of the elongated structures (filaments or ridges) have counterparts in the Stokes Q and U maps. This implies that in these structures, the magnetic field has a well-defined mean direction on the scales probed by Planck . Furthermore, the ridges are predominantly aligned with the magnetic field measured on the structures. This statistical trend becomes more striking for decreasing column density and, as ex-  \npected from the potential e GLYPH<11> ects of projection, for increasing polarization fraction. There is no alignment for the highest column density ridges in the N H & 10 22 cm GLYPH<0> 2 sample. Planck Collaboration Int. XXXIII (2014) studied the polarization properties of three nearby filaments, showing by geometrical modelling that the magnetic field in those representative regions has a well-defined mean direction that is di GLYPH<11> erent from the field orientation in the surroundings.  \nIn the present work, we quantitatively evaluate the relative orientation of the magnetic field inferred from the Planck polarization observations with respect to the gas column density structures, using the histogram of relative orientations (HRO, Soler et al. 2013). The HRO is a novel statistical tool that quantifies the relative orientation of each polarization measurement with respect to the column density gradient, making use of the unprecedented statistics provided by the Planck polarization observations. The HRO can also be evaluated in both 3D simulation data cubes and synthetic observations, thereby providing a direct comparison between observations and the physical conditions included in MHD simulations. We compare the results of the HRO applied to the Planck observations with the results of the same analysis applied to synthetic observations of MHD simulations of super-Alfv'enic, Alfv'enic, and sub-Alfv'enic turbulence.  \nThus by comparison with numerical simulations of MHD turbulence, the HRO provides estimates of the magnetic field strength without any of the assumptions involved in the DCF method. For comparison, we estimate B ? using the DCF method and the related method described by Hildebrand et al. (2009) ( DCF + SF , for DCF plus structure function) and provide a critical assessment of their applicability.  \nThis paper is organized as follows. Section 2 introduces the Planck 353 GHz polarization observations, the gas column density maps, and the CO line observations used to derive the velocity information. The particular regions where we evaluate the relative orientation between the magnetic field and the column density structures are presented in Sect. 3. Section 4 describes the statistical tools used for the study of these relative orientations. In Sect. 5 we discuss our results and their implications in the general picture of cloud formation. Finally, Sect. 6 summarizes the main results. Additional information on the selection of the polarization data, the estimation of uncertainties a GLYPH<11> ecting the statistical method, and the statistical significance of the relative orientation studies can be found in Appendices A, B, and C, respectively. Appendix D presents alternative estimates of the magnetic field strength in each region.", '2.1. Thermal dust polarization': 'Over the whole sky Planck observed the linear polarization (Stokes Q and U ) in seven frequency bands from 30 to 353 GHz (Planck Collaboration I 2014). In this study, we used data from the High Frequency Instrument (HFI, Lamarre et al. 2010) at 353 GHz, the highest frequency band that is sensitive to polarization. Towards MCs the contribution of the cosmic microwave background (CMB) polarized emission is negligible at 353 GHz, making this the Planck map that is best suited to studying the spatial structure of the dust polarization (Planck Collaboration Int. XIX 2014; Planck Collaboration Int. XX 2014).  \nWe used the Stokes Q and U maps and the associated noise maps made from five independent consecutive sky surveys of the Planck cryogenic mission, which together correspond to the  \nDR3 (delta-DX11d) internal data release. We refer to previous Planck publications for the data processing, map making, photometric calibration, and photometric uncertainties (Planck Collaboration II 2014; Planck Collaboration V 2014; Planck Collaboration VI 2014; Planck Collaboration VIII 2014). As in the first Planck polarization papers, we used the International Astronomical Union (IAU) conventions for the polarization angle, measured from the local direction to the north Galactic pole with positive values increasing towards the east.  \nThe maps of Q , U , their respective variances GLYPH<27> 2 Q , GLYPH<27> 2 U , and their covariance GLYPH<27> QU are initially at 4. 0 8 resolution in HEALPix format 3 with a pixelization at N side = 2048, which corresponds to an e GLYPH<11> ective pixel size of 1. 0 7. To increase the signal-to-noise ratio (S / N) of extended emission, we smoothed all the maps to 10 0 resolution using a Gaussian approximation to the Planck beam and the covariance smoothing procedures described in Planck Collaboration Int. XIX (2014).  \nThe maps of the individual regions are projected and resampled onto a Cartesian grid by using the gnomonic projection procedure described in Paradis et al. (2012). The HRO analysis is performed on these projected maps.', '2.2. Column density': 'We used the dust optical depth at 353 GHz ( GLYPH<28> 353) as a proxy for the gas column density ( N H). The GLYPH<28> 353 map (Planck Collaboration XI 2014) was derived from the all-sky Planck intensity observations at 353, 545, and 857 GHz, and the IRAS observations at 100 GLYPH<22> m, which were fitted using a modified black body spectrum. Other parameters obtained from this fit are the temperature and the spectral index of the dust opacity. The GLYPH<28> 353 map, computed initially at 5 0 resolution, was smoothed to 10 0 to match the polarization maps. The errors resulting from smoothing the product GLYPH<28> 353 map, rather than the underlying data, are negligible compared to the uncertainties in the dust opacity and do not significantly a GLYPH<11> ect the results of this study.  \nTo scale from GLYPH<28> 353 to N H, following Planck Collaboration XI (2014), we adopted the dust opacity found using Galactic extinction measurements of quasars,  \nGLYPH<28> 353 = N H = 1 : 2 GLYPH<2> 10 GLYPH<0> 26 cm 2 : (1)  \nVariations in dust opacity are present even in the di GLYPH<11> use ISM and the opacity increases systematically by a factor of 2 from the di GLYPH<11> use to the denser ISM (Planck Collaboration XXIV 2011; Martin et al. 2012; Planck Collaboration XI 2014), but our results do not critically depend on this calibration .', '3. Analysed regions': "The selected regions, shown in Fig. 2, correspond to nearby ( d < 450 pc) MCs, whose characteristics are well studied and can be used for cloud-to-cloud comparison (Poppel 1997; Reipurth 2008). Their properties are summarized in Table 1, which includes: Galactic longitude l and latitude b at the centre of the field; field size GLYPH<1> l GLYPH<2> GLYPH<1> b ; estimate of distance; mean and maximum total column densities from dust, h N H i and max ( N H), respectively; and mean H2 column density from CO.  \nIn the table the regions are organized from the nearest to the farthest in three groups: (a) regions located at d GLYPH<25> 150 pc, namely Taurus, Ophiuchus, Lupus, Chamaeleon-Musca, and Corona Australis (CrA); (b) regions located at d GLYPH<25> 300 pc,  \nFig. 2. Locations and sizes of the regions selected for analysis. The background map is the gas column density, N H, derived from the dust optical depth at 353 GHz (Planck Collaboration XI 2014).  \n<!-- image -->  \nTable 1. Locations and properties of the selected regions  \n- c Using the line integral W CO over GLYPH<0> 10 < v k = (kms GLYPH<0> 1 ) < 10 from the Dame et al. (2001) survey and X CO = 1 : 8 GLYPH<2> 10 20 cm GLYPH<0> 2 K GLYPH<0> 1 km GLYPH<0> 1 s.  \nAquila Rift and Perseus; and (c) regions located at d GLYPH<25> 450 pc, IC 5146, Cepheus, and Orion.  \nAmong the clouds in the first group (all shown in Fig. 3, left column) are Ophiuchus and Lupus, which are two regions with di GLYPH<11> erent star-forming activities but are close neighbours within an environment disturbed by the Sco-Cen OB association (Wilking et al. 2008; Comer'on 2008). Chamaeleon-Musca is a region evolving in isolation, and it is relatively unperturbed (Luhman 2008). Taurus (see also Fig. 1), a cloud with low-mass star formation, appears to be formed by the material swept up by an ancient superbubble centred on the Cas-Tau group (Kenyon et al. 2008). Finally, CrA is one of the nearest regions with recent intermediate- and low-mass star formation, possibly formed by  \na high-velocity cloud impact on the Galactic plane (Neuhauser &Forbrich 2008).  \nIn the second group we consider Aquila Rift and Perseus, shown in Fig. 4 (left column). Aquila Rift is a large complex of dark clouds where star formation proceeds in isolated pockets (Eiroa et al. 2008; Prato et al. 2008). The Perseus MC is the most active site of on-going star formation within 300 pc of the Sun. It features a large velocity gradient and is located close to hot stars that might have impacted its structure (Bally et al. 2008).  \nIn the third group are IC 5146, Cepheus, and Orion, shown in Fig. 5 (left column). IC 5146 is an MC complex in Cygnus. It includes an open cluster surrounded by a bright optical nebulosity called the Cocoon nebula, and a region of embedded lower-  \nFig. 3. Left : Columm density map, log 10 ( N H = cm GLYPH<0> 2 ), overlaid with magnetic field pseudo-vectors whose orientations are inferred from the Planck 353 GHz polarization observations. The length of the pseudo-vectors is normalized so does not reflect the polarization fraction. In this first group, the regions analysed are, from top to bottom, Taurus, Ophiuchus, Lupus, Chamaeleon-Musca, and CrA. Right : HROs for the lowest, an intermediate, and the highest N H bin (black, blue, and red, respectively). For a given region, bins have equal numbers of selected pixels (see Sect. 4.1.1 and Appendix A) within the N H ranges labelled. The intermediate bin corresponds to selected pixels near the blue contours in the column density images. The horizontal dashed line corresponds to the average per angle bin of 15 GLYPH<14> . The widths of the shaded areas for each histogram correspond to the GLYPH<6> 1 GLYPH<27> uncertainties related to the histogram binning operation. Histograms peaking at 0 GLYPH<14> correspond to B ? predominantly aligned with isoN H contours. Histograms peaking at 90 GLYPH<14> and / or GLYPH<0> 90 GLYPH<14> correspond to B ? predominantly perpendicular to isoN H contours.  \n<!-- image -->  \nmass star formation known as the IC 5146 Northern 4 Streamer (Harvey et al. 2008). The Cepheus Flare, called simply Cepheus in this study, is a large complex of dark clouds that seems to belong to an even larger expanding shell from an old supernova remnant (Kun et al. 2008). Orion is a dark cloud complex with on-going high and low mass star formation, whose structure appears to be a GLYPH<11> ected by multiple nearby hot stars (Bally 2008).  \nWhen taking background / foreground emission and noise within these regions into account, pixels are selected for analysis according to criteria for the gradient of the column density (Appendix A.2) and the polarization (Appendix A.2).", '4.1.1. Histogram of relative orientations': "We quantify the relative orientation of the magnetic field with respect to the column density structures using the HRO (Soler et al. 2013). The column density structures are characterized by their gradients, which are by definition perpendicular to the isocolumn density curves (see calculation in Appendix B.1). The gradient constitutes a vector field that we compare pixel by pixel  \nFig. 3. (continued).  \n<!-- image -->  \nFig. 4. Same as Fig. 3 for the second group, Aquila Rift and Perseus.  \n<!-- image -->  \nto the magnetic field orientation inferred from the polarization maps.  \nIn practice we use GLYPH<28> 353 as a proxy for N H (Sect. 2.2). The angle GLYPH<30> between B ? and the tangent to the GLYPH<28> 353 contours is evaluated using 5  \nGLYPH<30> = arctan GLYPH<16> jr GLYPH<28> 353 GLYPH<2> ˆ E j ; r GLYPH<28> 353 GLYPH<1> ˆ E GLYPH<17> ; (2)  \nwhere, as illustrated in Fig. 6, r GLYPH<28> 353 is perpendicular to the tangent of the isoGLYPH<28> 353 curves, the orientation of the unit polarization pseudo-vector ˆ E , perpendicular to B ? , is characterized by the polarization angle  \n= 1 2 arctan( GLYPH<0> U ; Q ) ; (3)  \nand in Eq. (2), as implemented, the norm actually carries a sign when the range used for GLYPH<30> is between GLYPH<0> 90 GLYPH<14> and 90 GLYPH<14> .  \nThe uncertainties in GLYPH<30> due to the variance of the GLYPH<28> 353 map and the noise properties of Stokes Q and U at each pixel are characterized in Appendix B.  \nThe gradient technique is one of multiple methods for characterizing the orientation of structures in a scalar field. Other methods, which include the Hessian matrix analysis (Molinari et al. 2011; Planck Collaboration Int. XXXII 2014) and the inertia matrix (Hennebelle 2013), are appropriate for measuring the orientation of ridges, i.e., the central regions of filamentary structures. The gradient technique is sensitive to contours and in that sense it is better suited to characterizing changes in the relative orientation in extended regions, not just on the crests of structures (Soler et al. 2013; Planck Collaboration Int. XXXII 2014). Additionally, the gradient technique can sample multiple scales by increasing the size of the vicinity of pixels used for its calculation (derivative kernel; see Appendix B.1). Previous studies that assign an average orientation of the cloud (Tassis et al. 2009; Li et al. 2013) are equivalent to studying the relative orientation using a derivative kernel close to the cloud size.  \nThe selected pixels belong to the regions of each map where the magnitude of the gradient jr GLYPH<28> 353 j is greater than in a dif-  \nFig. 5. Same as Fig. 3 for the third group, IC 5146, Cepheus, and Orion.  \n<!-- image -->  \nfuse reference field (Appendix A). This selection criterion aims at separating the structure of the cloud from the structure of the background using the reference field as a proxy. For each region the selected reference field is the region with the same size and Galactic latitude that has the lowest average N H. See Appendix A.1.  \nIn addition to the selection on jr GLYPH<28> 353 j , we only consider pixels where the norms of the Stokes Q and U are larger than in the di GLYPH<11> use reference field, therefore minimizing the e GLYPH<11> ect of background / foreground polarization external to the cloud. The relative orientation angle, GLYPH<30> , is computed by using polarization measurements with a high S / Nin Stokes Q and U , i.e., only considering pixels with j Q j =GLYPH<27> Q or j U j =GLYPH<27> U > 3. This selection allows the unambiguous definition of ˆ E by constraining the uncertainty in the polarization angle (see Appendix A.2).  \nFig. 6. Schematic of the vectors involved in the calculation of the relative orientation angle GLYPH<30> .  \n<!-- image -->  \nOnce we have produced a map of relative orientations for selected pixels following Eq. (2), we divide the map into bins of N H containing an equal number of pixels and generate a histogram of GLYPH<30> for each bin. The shape of the histogram is used to evaluate the preferred relative orientation in each bin directly. A concave histogram, peaking at 0 GLYPH<14> , corresponds to the preferred alignment of B ? with the N H contours. A convex histogram, peaking at 90 GLYPH<14> and / or GLYPH<0> 90 GLYPH<14> , corresponds to the preferred orientation of B ? perpendicular to the N H contours.  \nThe HROs in each region are computed in 25 N H bins having equal numbers of selected pixels (10 bins in two regions with fewer pixels, CrA and IC 5146). The number of N H bins is determined by requiring enough bins to resolve the highest N H regions and at the same time maintaining enough pixels per N H bin to obtain significant statistics from each histogram. The typical number of pixels per bin of N H ranges from approximately 600 in CrA to around 4 000 in Chamaeleon-Musca. We use 12 angle bins of width 15 GLYPH<14> .  \nThe HROs of the first group of regions, the nearest at d GLYPH<25> 150 pc, are shown in the right-hand column of Fig. 3. For the sake of clarity, we only present the histograms that correspond to three bins, namely the lowest and highest N H and an intermediate N Hvalue. The intermediate bin is the 12th (sixth in two regions with fewer pixels, CrA and IC 5146), and it corresponds to pixels near the blue contour in the image in the left-hand column of Fig. 3. The widths of the shaded areas for each histogram correspond to the 1 GLYPH<27> uncertainties related to the histogram binning operation, which are greater than the uncertainties produced by the variances of Q , U , and GLYPH<28> 353 (Appendix B). The sharp and narrow features ('jitter') in the HROs are independent of these variances. They are the product of sampling the spatial correlations in the magnetic field over a finite region of the sky together with  \nTable 2. Fit of GLYPH<24> vs. log 10 ( N H = cm GLYPH<0> 2 ). a  \nthe histogram binning; these features average out when evaluating the relative orientation over larger portions of the sky (Planck Collaboration Int. XXXII 2014).  \nAlthough often asymmetric, most histograms reveal a change in the preferred relative orientation across N H bins. The most significant feature in the HROs of Taurus, Ophiuchus, and Chamaeleon-Musca is the drastic change in relative orientation from parallel in the lowest N H bin to perpendicular in the highest N H bin. In Lupus the behaviour at low N H is not clear, but at high N H it is clearly perpendicular. In contrast, CrA tends to show B ? as parallel in the intermediate N H bin, but no preferred orientation in the other N H bins.  \nThe HROs of the clouds located at d GLYPH<25> 300 pc, Aquila Rift and Perseus, are shown in the right-hand column of Fig. 4. They indicate that the relative orientation is usually parallel in the lowest N H bins and perpendicular in the highest N H bins.  \nThe HROs of the third group, located at d GLYPH<25> 400 GLYPH<0> 450 pc, IC 5146, Cepheus, and Orion, are presented in the right-hand column of Fig. 5. In both IC 5146 and Orion the HROs for the highest N H bins reveal a preferred orientation of the field perpendicular to the N H contours (Orion is quite asymmetric), whereas the HROs corresponding to the low and intermediate N H bins reveal a preferred alignment of the field with N H structures. This trend is also present, but less pronounced, in the Cepheus region.", '4.1.2. Histogram shape parameter GLYPH<24>': 'The changes in the HROs are quantified using the histogram shape parameter GLYPH<24> , defined as  \nGLYPH<24> = A c GLYPH<0> A e A c + A e ; (4)  \nwhere A c is the area in the centre of the histogram ( GLYPH<0> 22. GLYPH<14> 5 < GLYPH<30> < 22. GLYPH<14> 5) and A e the area in the extremes of the histogram ( GLYPH<0> 90. GLYPH<14> 0 < GLYPH<30> < GLYPH<0> 67. GLYPH<14> 5 and 67. GLYPH<14> 5 < GLYPH<30> < 90. GLYPH<14> 0). The value of GLYPH<24> , the result of the integration of the histogram over 45 GLYPH<14> ranges, is independent of the number of bins selected to represent the histogram if the bin widths are smaller than the integration range.  \nA concave histogram corresponding to B ? mostly aligned with N H contours would have GLYPH<24> > 0. A convex histogram corresponding to B ? mostly perpendicular to N H contours would have GLYPH<24> < 0. A flat histogram corresponding to no preferred relative orientation would have GLYPH<24> GLYPH<25> 0 .  \nThe uncertainty in GLYPH<24> , GLYPH<27>GLYPH<24> , is obtained from  \nGLYPH<27> 2 GLYPH<24> = 4 ( A 2 e GLYPH<27> 2 A c + A 2 c GLYPH<27> 2 A e ) ( A c + A e) 4 : (5)  \nFig. 7. Histogram shape parameter GLYPH<24> (Eqs. 4 and 5) calculated for the di GLYPH<11> erent N H bins in each region. The cases GLYPH<24> > 0 and GLYPH<24> < 0 correspond to the magnetic field oriented mostly parallel and perpendicular to the structure contours, respectively. For GLYPH<24> GLYPH<25> 0 there is no preferred orientation. The displayed values of C HRO and X HRO were calculated from Eq. (6) and correspond to the grey dashed line in each plot.  \n<!-- image -->  \nThe variances of the areas, GLYPH<27> 2 A e and GLYPH<27> 2 A c , characterize the jitter of the histograms. If the jitter is large, GLYPH<27>GLYPH<24> is large compared to j GLYPH<24> j and the relative orientation is indeterminate. The jitter depends on the number of bins in the histogram, but GLYPH<24> does not.  \nFigure 7 illustrates the change in GLYPH<24> as a function of log 10 ( N H = cm GLYPH<0> 2 ) of each bin. For most of the clouds, GLYPH<24> is positive in the lowest and intermediate N H bins and negative or close to zero in the highest bins. The most pronounced changes  \nin GLYPH<24> from positive to negative are seen in Taurus, ChamaeleonMusca, Aquila Rift, Perseus, and IC 5146.  \nThe trend in GLYPH<24> vs. log 10 ( N H = cm GLYPH<0> 2 ) can be fit roughly by a linear relation  \nGLYPH<24> = C HRO [log 10 ( N H = cm GLYPH<0> 2 ) GLYPH<0> X HRO] : (6)  \nThe values of C HRO and X HRO in the regions analysed are summarized in Table 2. For the clouds with a pronounced change in relative orientation the slope C HRO is steeper than about GLYPH<0> 0 : 5,  \nFig. 8. Same as Fig. 3 for the two test regions located directly east (ChamEast) and south of Chamaeleon-Musca (ChamSouth).  \n<!-- image -->  \nand the value X HRO for the log 10 ( N H = cm GLYPH<0> 2 ) at which the relative orientation changes from parallel to perpendicular is greater than about 21.7. Ophiuchus, Lupus, Cepheus, and Orion are intermediate cases, where GLYPH<24> definitely does not go negative in the data, but seems to do so by extrapolation; these tend to have a shallower C HRO and / or a higher X HRO.  \nThe least pronounced change in GLYPH<24> is seen in CrA, where GLYPH<24> is consistently positive in all bins, and the slope is very flat. We applied the HRO analysis to a pair of test regions (Fig. 8) with even lower N H values (log 10 ( N H = cm GLYPH<0> 2 ) < 21 : 6; see also Fig. 11 below) located directly south and directly east of the ChamaeleonMusca region. As in CrA, we find that B ? is mostly parallel to the N H contours, a fairly constant GLYPH<24> , and no indication of predominantly perpendicular relative orientation.', '4.2. Comparisons with previous studies': "The above trends in relative orientation between B ? and the N H contours in targeted MCs, where B ? tend to become perpendicular to the N H contours at high N H, agree with the results of the Hessian matrix analysis applied to Planck observations over the  \nwhole sky, as reported in Fig. 15 of Planck Collaboration Int. XXXII (2014).  \nEvidence for preferential orientations in sections of some regions included in this study has been reported previously. The Taurus region has been the target of many studies. Moneti et al. (1984) and Chapman et al. (2011) find evidence using infrared polarization of background stars for a homogeneous magnetic field perpendicular to the embedded dense filamentary structure. High-resolution submillimetre observations of intensity with Herschel find evidence of faint filamentary structures ('striations'), which are well correlated with the magnetic field orientation inferred from starlight polarization (Palmeirim et al. 2013) and perpendicular to the filament B211. Heyer et al. (2008) report a velocity anisotropy aligned with the magnetic field, which can be interpreted as evidence of the channeling e GLYPH<11> ect of the magnetic fields. But the magnetic field in B211 and the dense filamentary structures are not measured directly. As described above, using Planck polarization we find that B ? is mostly aligned with the lowest N H contours (20 : 8 < log 10 ( N H = cm GLYPH<0> 2 ) < 21 : 3) (see also Planck Collaboration Int. XXXIII 2014), although the aligned structures do not correspond to the striations,  \nwhich are not resolved at 10 0 resolution. However, we also find that at higher N H the relative orientation becomes perpendicular.  \nSimilar studies in other regions have found evidence of striations correlated with the starlight-inferred magnetic field orientation and perpendicular to the densest filamentary structures. The regions studied include Serpens South (Sugitani et al. 2011), which is part of Aquila Rift in this study, Musca (Pereyra & Magalh˜aes 2004), and the Northern Lupus cloud (Matthews et al. 2014). By studying Planck polarization in larger regions around the targets of these previous observations, we show that a systematic change in relative orientation is the prevailing statistical trend in clouds that reach log 10 ( N H = cm GLYPH<0> 2 ) & 21 : 5.", '5.1.1. Spatial distribution of the HRO signal': 'The maps obtained using Eq. (2) characterize the relative orientation in each region, without assuming an organization of the N H structures in ridges or filaments; HROs basically just sample the orientation of N H contours. However, the resulting maps of the relative orientation angle, shown for the Taurus region in Fig. 9, reveal that the regions that are mostly oriented parallel or perpendicular to the field form continuous patches, indicating that the HRO signal is not only coming from variations in the field or the N H contours at the smallest scales in the map.  \nHRO analysis on larger scales in a map can be achieved by considering a larger vicinity of pixels for calculating the gradient r GLYPH<28> 353. This operation is equivalent to calculating the nextneighbour gradient on a map first smoothed to the scale of interest. The results for relative orientation after smoothing to resolutions of 15 0 , 30 0 , and 60 0 are illustrated for the Taurus region in Fig. 9. Figure 10 shows that the corresponding HROs have a similar behaviour for the three representative N H values. These results confirm that the preferred relative orientation is not particular to the smallest scales in the map, but corresponds to coherent structures in N H.  \nA study of the preferred orientation for the whole cloud would be possible by smoothing the column density and polarization maps to a scale comparable to the cloud size. The statistical significance of such a study would be limited to the number of clouds in the sample and would not be directly comparable to previous studies of relative orientation of clouds, where elongated structures were selected to characterize the mean orientation of each cloud (Li et al. 2013).', '5.1.2. Statistical significance of the HRO signal': 'Our results reveal a systemic change of GLYPH<24> with N H, suggesting a systematic transition from magnetic field mostly parallel to N H contours in the lowest N H bins to mostly perpendicular in the highest N H bins of the clouds studied. The statistical significance of this change can in principle be evaluated by considering the geometrical e GLYPH<11> ects that influence this distribution. In Appendix C.1, using simulations of Q and U maps, we eliminate the possibility that this arises from random magnetic fields, random spatial correlations in the field, or the large scale structure of the field. In Appendix C.2 we simply displace the Q and U maps spatially and repeat the analysis, showing that the systemic trend of GLYPH<24> vs. N H disappears for displacements greater than 1 GLYPH<14> .  \nUsing a set of Gaussian models, Planck Collaboration Int. XXXII (2014) estimated the statistical significance of this transition in terms of the relative orientation between two vectors in 3D and their projection in 2D. As these authors emphasized, two vectors that are close to parallel in 3D would be projected as parallel in 2D for almost all viewing angles for which the projections of both vectors have a non-negligible length, but on the other hand, the situation is more ambiguous seen in projection for two vectors that are perpendicular in 3D, because they can be projected as parallel in 2D depending on the angle of viewing. The quantitative e GLYPH<11> ects are illustrated by the simulations in Appendix C.3, where we consider distributions of vectors in 3D that are mostly parallel, mostly perpendicular, or have no preferred orientation. The projection tends to make vector pairs look more parallel in 2D, but the distribution of relative orientations in 2D is quite similar though not identical to the distribution in 3D. In particular, the signal of perpendicular orientation is not erased and we can conclude that two projected vectors with nonnegligible lengths that are close to perpendicular in 2D must also be perpendicular in 3D. This would apply to the mostly perpendicular orientation for the highest bin of N H in the Taurus region.', '5.2. Comparison with simulations of MHD turbulence': "As a complement to observations, MHD simulations can be used to directly probe the actual 3D orientation of the magnetic field B with respect to the density structures. The change in the relative orientation with N H was previously studied in MHD simulations using the inertia matrix and the HRO analysis (Hennebelle 2013; Soler et al. 2013). Soler et al. (2013) showed that in 3D, the change in the relative orientation is related to the degree of magnetization. If the magnetic energy is above or comparable to the kinetic energy (turbulence that is sub-Alfv'enic or close to equipartition), the less dense structures tend to be aligned with the magnetic field and the orientation progressively changes from parallel to perpendicular with increasing density. In the super-Alfv'enic regime, where the magnetic energy is relatively low, there appears to be no change in relative orientation with increasing density, with B and density structures being mostly parallel.  \nSoler et al. (2013) describe 2D synthetic observations of the MHD simulations. The synthetic observations are produced by integrating the simulation cubes along a direction perpendicular to the mean magnetic field and assuming a homogeneous dust grain alignment e GLYPH<14> ciency GLYPH<15> = 1 : 0. The angular resolution of the simulation is obtained by assuming a distance d = 150 pc and convolving the projected map with a Gaussian beam of 10 0 full width at half maximum (FWHM). The trends in the relative orientation with N H seen in 3D are also seen using the these 2D synthetic observations. Given that sub-Alfv'enic or close to Alfv'enic turbulence does not significantly disturb the well-ordered mean magnetic field, the orientation of B perpendicular to the isodensity contours is projected well for lines of sight that are not close to the mean magnetic field orientation. In contrast, the projected relative orientation produced by super-Alfv'enic turbulence does not necessarily reflect the relative orientation in 3D as a result of the unorganized field structure.  \nThe direct comparison between the HROs of the regions in this study and of the synthetic observations is presented in Fig. 11. The trends in the relative orientation parameter, GLYPH<24> , show that the simulation with super-Alfv'enic turbulence does not undergo a transition in relative orientation from parallel to perpendicular for log 10 ( N H = cm GLYPH<0> 2 ) < 23. In contrast, most of the observed clouds show a decrease in GLYPH<24> with increasing N H, close  \nFig. 9. Maps of the absolute value of the relative orientation angle, j GLYPH<30> j , in the Taurus region. These maps are produced after smoothing the input maps to beam FWHMs of 10 0 , 15 0 , 30 0 , and 60 0 and then resampling the grid to sample each beam FWHM with the same number of pixels. The regions in red correspond to B ? close to perpendicular to N H structures. The regions in blue correspond to B ? close to parallel to N H structures. The black contour, corresponding to the N H value of the intermediate contour introduced in Fig. 3, provides a visual reference to the cloud structure.  \n<!-- image -->  \nto the trends seen from the simulations with Alfv'enic or subAlfv'enic turbulence for log 10 ( N H = cm GLYPH<0> 2 ) < 23. Furthermore, X HRO, the value of log 10 ( N H = cm GLYPH<0> 2 ) where GLYPH<24> goes through zero, is near 21.7, which is consistent with the behaviour seen in the simulations with super-Alfv'enic or Alfv'enic turbulence. Given that the physical conditions in the simulations ( GLYPH<27> v = 2 : 0 km s GLYPH<0> 1 and n = 500 cm GLYPH<0> 3 ) are typical of those in the selected regions ( GLYPH<27> v k is given in Table D.1), the similarities in the dependence of GLYPH<24> on N H suggest that the strength of the magnetic field in most of the regions analysed would be about the same as the mean magnetic fields in the Alfv'enic and sub-Alfv'enic turbulence simulations, which are 3.5 and 11 GLYPH<22> G, respectively. However, more precise estimates of the magnetic field strength coming directly from the HROs would require further sampling of the magnetization in the MHD simulations and detailed modelling of the e GLYPH<11> ects of the line-of-sight integration.  \nIndirectly, the presence of a N H threshold in the switch in preferred relative orientation between B ? and the N H structures hints that gravity plays a significant role. By contrast, for the regions with low average N H (i.e., CrA and the two test regions; Fig. 11), there is little change in GLYPH<24> and certainly no switch in the preferred relative orientation to perpendicular.", '5.3. Physics of the relative orientation': "The finding of dense structures mostly perpendicular to the magnetic field (and the small mass-to-flux ratios discussed in Appendix D.4) suggests that the magnetic field in most of the observed regions is significant for the structure and dynamics. However, discerning the underlying geometry is not obvious. As one guide, a frozen-in and strong interstellar magnetic field would naturally cause a self-gravitating, static cloud to become oblate, with its major axis perpendicular to the field lines, be-  \nFig. 10. HROs of the Taurus region after smoothing the input maps to beam FWHMs of 15 0 , 30 0 , and 60 0 , shown from left to right, respectively.  \n<!-- image -->  \ncause gravitational collapse would be restricted to occurring along field lines (Mouschovias 1976a,b). In the case of less dense structures that are not self-gravitating, the velocity shear can stretch matter and field lines in the same direction, thereby producing aligned structures, as discussed in Hennebelle (2013) and Planck Collaboration Int. XXXII (2014).  \nIf the MCs are isolated entities and the magnetic field is strong enough to set a preferred direction for the gravitational collapse, the condensations embedded in the cloud are not very likely to have higher column densities than their surroundings (Nakano 1998). This means that the formation of dense substructures, such as prestellar cores and stars, by gravitational collapse would be possible only if the matter decouples from the magnetic field. This is possible through the decoupling between neutral and ionized species (ambipolar di GLYPH<11> usion, Mouschovias 1991; Li & Houde 2008) or through removal of magnetic flux from clouds via turbulent reconnection (Lazarian & Vishniac 1999; Santos-Lima et al. 2012).  \nAlternatively, if we regard MCs not as isolated entities but as the result of an accumulation of gas by large-scale flows (Ballesteros-Paredes et al. 1999; Hartmann et al. 2001; Koyama & Inutsuka 2002; Audit & Hennebelle 2005; Heitsch et al. 2006), the material swept up by colliding flows may eventually form a self-gravitating cloud. If the magnetic field is strong the accumulation of material is favoured along the magnetic field lines, thus producing dense structures that are mostly perpendicular to the magnetic field. The inflow of material might eventually increase the gravitational energy in parts of the cloud, thereby producing supercritical structures such as prestellar cores.  \nFor supersonic turbulence in the ISM and MCs, density structures can be formed by gas compression in shocks. If the turbulence is strong with respect to the magnetic field (super-Alfv'enic), gas compression by shocks is approximately isotropic; because magnetic flux is frozen into matter, field lines are dragged along with the gas, forming structures that tend to be aligned with the field. If the turbulence is weak with respect to the magnetic field (sub-Alfv'enic), the fields produce a clear anisotropy in MHD turbulence (Sridhar & Goldreich 1994; Goldreich & Sridhar 1995; Matthaeus et al. 2008; Banerjee et al. 2009) and compression by shocks that is favored to occur along the magnetic field lines, creating structures perpendicular to the field. The cold phase gas that constitutes the cloud receives no information about the original flow direction because the mag-  \ntic field redistributes the kinetic energy of the inflows (Heitsch et al. 2009; Inoue & Inutsuka 2009; Burkhart et al. 2014). This seems to be the case in most of the observed regions, where the mostly perpendicular relative orientation between the magnetic field and the high column density structures is an indication of the anisotropy produced by the field.  \nThe threshold of log 10 ( N H = cm GLYPH<0> 2 ) GLYPH<25> 21 : 7 above which the preferential orientation of B ? switches to being perpendicular to the N H contours is intriguing. Is there a universal threshold column density that is independent of the particular MC environment and relevant in the context of star formation? In principle, this threshold might be related to the column density of filaments at which substructure forms, as reported in an analysis of Herschel observations (Arzoumanian et al. 2013), but the Planck polarization observations leading to B ? do not fully resolve such filamentary structures. In principle, this threshold might also be related to the column density at which the magnetic field starts scaling with density, according to the Zeeman e GLYPH<11> ect observations of B k (Figure 7 in Crutcher 2012). However, establishing such relationships requires further studies with MHD simulations to identify what densities and scales influence the change in relative orientation between B ? and N H structures and to model the potential imprint in B k observations and in B ? observations to be carried out at higher resolution.", '5.4. Effect of dust grain alignment': "Throughout this study we assume that the polarized emission observed by Planck at 353 GHz is representative of the projected morphology of the magnetic field in each region; i.e., we assume a constant dust grain alignment e GLYPH<14> ciency ( GLYPH<15> ) that is independent of the local environment. Indeed, observations and MHD simulations under this assumption (Planck Collaboration Int. XIX 2014; Planck Collaboration Int. XX 2014) indicate that depolarization e GLYPH<11> ects at large and intermediate scales in MCs might arise from the random component of the magnetic field along the line of sight. On the other hand, the sharp drop in the polarization fraction at N H > 10 22 cm GLYPH<0> 2 (reported in Planck Collaboration Int. XIX 2014), when seen at small scales, might be interpreted in terms of a decrease of GLYPH<15> with increasing column density (Matthews et al. 2001; Whittet et al. 2008).  \nA leading theory for the process of dust grain alignment involves radiative torques by the incident radiation (Lazarian & Hoang 2007; Hoang & Lazarian 2009; Andersson 2015). A crit-  \nFig. 11. Histogram shape parameter GLYPH<24> (Eqs. 4 and 5) calculated for the di GLYPH<11> erent N H bins in each region. Top: relative orientation in synthetic observations of simulations with super-Alfv'enic (blue), Alfv'enic (green), and sub-Alfv'enic (red) turbulence, as detailed in Soler et al. (2013). Middle: relative orientation in the regions selected from the Planck all-sky observations, from Fig. 7. The blue data points correspond to the lowest N H regions (CrA and the test regions in Fig. 8, ChamSouth and ChamEast) and the orange correspond to the rest of the clouds. Bottom: comparison between the trends in the synthetic observations (in colours) and the regions studied (grey). The observed smooth transition from mostly parallel ( GLYPH<24> > 0) to perpendicular ( GLYPH<24> < 0) is similar to the transition in the simulations for which the turbulence is Alfv'enic or sub-Alfv'eic.  \n<!-- image -->  \nical parameter for this mechanism is the ratio between the dust grain size and the radiation wavelength. As the dust column density increases, only the longer wavelength radiation penetrates the cloud and the alignment decreases. Grains within a cloud (without embedded sources) should have lower GLYPH<15> than those at the periphery of the same cloud. There is evidence for this from near-infrared interstellar polarization and submillimetre polarization along lines of sight through starless cores (Jones et al. 2015), albeit on smaller scales and higher column densities than considered here. If GLYPH<15> inside the cloud is very low, the observed  \npolarized intensity would arise from the dust in the outer layers, tracing the magnetic field in the 'skin' of the cloud. Then the observed orientation of B ? is not necessarily correlated with the column density structure, which is seen in total intensity, or with the magnetic field deep in the cloud.  \nSoler et al. (2013) presented the results of HRO analysis on a series of synthetic observations produced using models of how GLYPH<15> might decrease with increasing density. They showed that with a steep decrease there is no visible correlation between the inferred magnetic field orientation and the highN H structure, corresponding to nearly flat HROs.  \nThe HRO analysis of MCs carried out here reveals a correlation between the polarization orientation and the column density structure. This suggests that the dust polarized emission samples the magnetic field structure homogeneously on the scales being probed at the resolution of the Planck observations or, alternatively, that the field deep within highN H structures has the same orientation of the field in the skin.", '6. Conclusions': "We have presented a study of the relative orientation of the magnetic field projected on the plane of the sky ( B ? ), as inferred from the Planck dust polarized thermal emission, with respect to structures detected in gas column density ( N H). The relative orientation study was performed by using the histogram of relative orientations (HRO), a novel statistical tool for characterizing extended polarization maps. With the unprecedented statistics of polarization observations in extended maps obtained by Planck , we analyze the HRO in regions with di GLYPH<11> erent column densities within ten nearby molecular clouds (MCs) and two test fields.  \nIn most of the regions analysed we find that the relative orientation between B ? and N H structures changes systematically with N H from being parallel in the lowest column density areas to perpendicular in the highest column density areas. The switch occurs at log 10 ( N H = cm GLYPH<0> 2 ) GLYPH<25> 21 : 7. This change in relative orientation is particularly significant given that projection tends to produce more parallel pseudo-vectors in 2D (the domain of observations) than exist in 3D.  \nThe HROs in these MCs reveal that most of the high N H structures in each cloud are mostly oriented perpendicular to the magnetic field, suggesting that they may have formed by material accumulation and gravitational collapse along the magnetic field lines. According to a similar study where the same method was applied to MHD simulations, this trend is only possible if the turbulence is Alfv'enic or sub-Alfv'enic. This implies that the magnetic field is significant for the gas dynamics on the scales sampled by Planck . The estimated mean magnetic field strength is about 4 and 12 GLYPH<22> G for the case of Alfv'enic and sub-Alfv'enic turbulence, respectively.  \nWealso estimate the magnetic field strength in the MCs studied using the DCF and DCF + SF methods. The estimates found seem consistent with the above values from the HRO analysis, but given the assumptions and systematic e GLYPH<11> ects involved, we recommend that these rough estimates be treated with caution. According to these estimates the analysed regions appear to be magnetically sub-critical. This result is also consistent with the conclusions of the HRO analysis. Specific tools, such as the DCF and DCF + SF methods, are best suited to the scales and physical conditions in which their underlying assumptions are valid. The study of large polarization maps covering multiple scales calls for generic statistical tools, such as the HRO, for characterizing their properties and establishing a direct relation to the physical conditions included in MHD simulations.  \nThe study of the structure on smaller scales is beyond the scope of this work, however, the presence of gravitationally bound structures within the MCs, such as prestellar cores and stars, suggests that the role of magnetic fields is changing on di GLYPH<11> erent scales. Even if the magnetic field is important in the accumulation of matter that leads to the formation of the cloud, e GLYPH<11> ects such as matter decoupling from the magnetic field and the inflow of matter from the cloud environment lead to the formation of magnetically supercritical structures on smaller scales. Further studies will help to identify the dynamical processes that connect the MC structure with the process of star formation.  \nAcknowledgements. The development of Planck has been supported by: ESA; CNES and CNRS / INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MICINN, JA, and RES (Spain); Tekes, AoF, and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER / SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT / MCTES (Portugal); and PRACE (EU). A description of the Planck Collaboration and a list of its members, including the technical or scientific activities in which they have been involved, can be found at http://www.sciops.esa.int/index.php? project=planck&page=Planck\\_Collaboration . The research leading to these results has received funding from the European Research Council under the European Union's Seventh Framework Programme (FP7 / 2007-2013) / ERC grant agreement No. 267934.", 'References': "- Andersson, B.-G. 2015, in Astrophysics and Space Science Library, Vol. 407, Astrophysics and Space Science Library, ed. A. Lazarian, E. M. de Gouveia Dal Pino, & C. Melioli, 59 Arzoumanian, D., Andr'e, P., Peretto, N., & Konyves, V. 2013, A&A, 553, A119 Audit, E. & Hennebelle, P. 2005, A&A, 433, 1 Ballesteros-Paredes, J., Hartmann, L., & V'azquez-Semadeni, E. 1999, ApJ, 527, 285 Bally, J. 2008, Overview of the Orion Complex, ed. B. Reipurth, 459 Bally, J., Walawender, J., Johnstone, D., Kirk, H., & Goodman, A. 2008, The Perseus Cloud, ed. B. Reipurth, 308 Banerjee, R., V'azquez-Semadeni, E., Hennebelle, P., & Klessen, R. S. 2009, MNRAS, 398, 1082 Benoˆıt, A., Ade, P., Amblard, A., et al. 2004, A&A, 424, 571 Bergin, E. A. & Tafalla, M. 2007, ARA&A, 45, 339 Bierman, E. M., Matsumura, T., Dowell, C. D., et al. 2011, ApJ, 741, 81 Burkhart, B., Lazarian, A., Le˜ao, I. C., de Medeiros, J. R., & Esquivel, A. 2014, ApJ, 790, 130 Cabral, B. & Leedom, L. C. 1993, in Special Interest Group on GRAPHics and Interactive Techniques Proceedings., Special Interest Group on GRAPHics and Interactive Techniques Proceedings., 263-270 Chandrasekhar, S. & Fermi, E. 1953, ApJ, 118, 113 Chapman, N. L., Goldsmith, P. F., Pineda, J. L., et al. 2011, ApJ, 741, 21 Comer'on, F. 2008, The Lupus Clouds, ed. B. Reipurth, 295 Crutcher, R. M. 2005, in American Institute of Physics Conference Series, Vol. 784, Magnetic Fields in the Universe: From Laboratory and Stars to Primordial Structures., ed. E. M. de Gouveia dal Pino, G. Lugones, & A. Lazarian, 129-139 Crutcher, R. M. 2012, ARA&A, 50, 29 Crutcher, R. M., Nutter, D. J., Ward-Thompson, D., & Kirk, J. M. 2004, ApJ, 600, 279 Crutcher, R. M., Wandelt, B., Heiles, C., Falgarone, E., & Troland, T. H. 2010, ApJ, 725, 466 Dame, T. M., Hartmann, D., & Thaddeus, P. 2001, ApJ, 547, 792 Davis, L. 1951, Physical Review, 81, 890 Davis, Jr., L. & Greenstein, J. L. 1951, ApJ, 114, 206 Dotson, J. L. 1996, ApJ, 470, 566 Dotson, J. L., Davidson, J., Dowell, C. D., Schleuning, D. A., & Hildebrand, R. H. 2000, ApJS, 128, 335 Draine, B. T. 2011, Physics of the Interstellar and Intergalactic Medium Eiroa, C., Djupvik, A. A., & Casali, M. M. 2008, The Serpens Molecular Cloud, ed. B. Reipurth, 693 Elmegreen, B. G. & Scalo, J. 2004, ARA&A, 42, 211 Falceta-Gonc¸alves, D., Lazarian, A., & Kowal, G. 2008, ApJ, 679, 537 Franco, G. A. P., Alves, F. O., & Girart, J. M. 2010, ApJ, 723, 146 Goldreich, P. & Sridhar, S. 1995, ApJ, 438, 763\n- G'orski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759\n- Hartmann, L., Ballesteros-Paredes, J., & Bergin, E. A. 2001, ApJ, 562, 852\n- Harvey, P. M., Huard, T. L., Jørgensen, J. K., et al. 2008, ApJ, 680, 495 Heiles, C. & Haverkorn, M. 2012, Space Sci. Rev., 166, 293  \nHeitsch, F., Slyz, A. D., Devriendt, J. E. G., Hartmann, L. W., & Burkert, A.  \n2006, ApJ, 648, 1052  \nHeitsch, F., Stone, J. M., & Hartmann, L. W. 2009, ApJ, 695, 248  \nHeitsch, F., Zweibel, E. G., Mac Low, M.-M., Li, P., & Norman, M. L. 2001,  \nApJ, 561, 800  \n- Hennebelle, P. 2013, A&A, 556, A153\n- Hennebelle, P. & Falgarone, E. 2012, A&A Rev., 20, 55  \nHeyer, M., Gong, H., Ostriker, E., & Brunt, C. 2008, ApJ, 680, 420  \n- Hildebrand, R. H. 1988, QJRAS, 29, 327\n- Hildebrand, R. H., Kirby, L., Dotson, J. L., Houde, M., & Vaillancourt, J. E. 2009, ApJ, 696, 567\n- Hiltner, W. A. 1949, Science, 109, 165\n- Hoang, T. & Lazarian, A. 2009, ApJ, 697, 1316\n- Houde, M., Fletcher, A., Beck, R., et al. 2013, ApJ, 766, 49\n- Houde, M., Rao, R., Vaillancourt, J. E., & Hildebrand, R. H. 2011, ApJ, 733, 109\n- Houde, M., Vaillancourt, J. E., Hildebrand, R. H., Chitsazzadeh, S., & Kirby, L. 2009, ApJ, 706, 1504\n- Inoue, T. & Inutsuka, S.-i. 2009, ApJ, 704, 161\n- Jones, T. J., Bagley, M., Krejny, M., Andersson, B.-G., & Bastien, P. 2015, AJ, 149, 31\n- Kenyon, S. J., G'omez, M., & Whitney, B. A. 2008, Low Mass Star Formation in the Taurus-Auriga Clouds, ed. B. Reipurth, 405  \nKnude, J. & Hog, E. 1998, A&A, 338, 897  \n- Kobulnicky, H. A., Molnar, L. A., & Jones, T. J. 1994, AJ, 107, 1433\n- Koch, P. M., Tang, Y.-W., & Ho, P. T. P. 2013, ApJ, 775, 77  \nKoyama, H. & Inutsuka, S.-i. 2002, ApJ, 564, L97  \n- Kun, M., Kiss, Z. T., & Balog, Z. 2008, Star Forming Regions in Cepheus, ed. B. Reipurth, 136\n- Lamarre, J., Puget, J., Ade, P. A. R., et al. 2010, A&A, 520, A9  \nLazarian, A. & Hoang, T. 2007, MNRAS, 378, 910 Lazarian, A. & Vishniac, E. T. 1999, ApJ, 517, 700 Li, H., Gri GLYPH<14> n, G. S., Krejny, M., et al. 2006, ApJ, 648, 340 Li, H.-b., Fang, M., Henning, T., & Kainulainen, J. 2013, MNRAS, 436, 3707  \n- Li, H.-b. & Houde, M. 2008, ApJ, 677, 1151\n- Luhman, K. L. 2008, Chamaeleon, ed. B. Reipurth, 169\n- Martin, P. G., Roy, A., Bontemps, S., et al. 2012, ApJ, 751, 28\n- Matthaeus, W. H., Pouquet, A., Mininni, P. D., Dmitruk, P., & Breech, B. 2008, Physical Review Letters, 100, 085003\n- Matthews, B. C., McPhee, C. A., Fissel, L. M., & Curran, R. L. 2009, ApJS, 182, 143\n- Matthews, B. C., Wilson, C. D., & Fiege, J. D. 2001, ApJ, 562, 400\n- Matthews, T. G., Ade, P. A. R., Angil'e, F. E., et al. 2014, ApJ, 784, 116\n- McKee, C. F. & Ostriker, E. C. 2007, ARA&A, 45, 565\n- McKee, C. F., Zweibel, E. G., Goodman, A. A., & Heiles, C. 1993, in Protostars and Planets III, ed. E. H. Levy & J. I. Lunine, 327\n- Molinari, S., Schisano, E., Faustini, F., et al. 2011, A&A, 530, A133\n- Moneti, A., Pipher, J. L., Helfer, H. L., McMillan, R. S., & Perry, M. L. 1984, ApJ, 282, 508\n- Montier, L., Plaszczynski, S., Levrier, F., et al. 2015, A&A, 574, A135\n- Mouschovias, T. C. 1976a, ApJ, 206, 753\n- Mouschovias, T. C. 1976b, ApJ, 207, 141\n- Mouschovias, T. C. 1991, ApJ, 373, 169\n- Myers, P. C. & Goodman, A. A. 1991, ApJ, 373, 509\n- Nakano, T. 1998, ApJ, 494, 587\n- Nakano, T. & Nakamura, T. 1978, PASJ, 30, 671\n- Neuhauser, R. & Forbrich, J. 2008, The Corona Australis Star Forming Region,  \ned. B. Reipurth, 735  \n- Ostriker, E. C., Stone, J. M., & Gammie, C. F. 2001, ApJ, 546, 980\n- Palmeirim, P., Andr'e, P., Kirk, J., et al. 2013, A&A, 550, A38\n- Paradis, D., Dobashi, K., Shimoikura, T., et al. 2012, A&A, 543, A103\n- Pereyra, A. & Magalh˜aes, A. M. 2004, ApJ, 603, 584\n- Planck Collaboration XXIV. 2011, A&A, 536, A24  \nPlanck Collaboration I. 2014, A&A, 571, A1  \nPlanck Collaboration II. 2014, A&A, 571, A2  \nPlanck Collaboration V. 2014, A&A, 571, A5  \n- Planck Collaboration VI. 2014, A&A, 571, A6  \nPlanck Collaboration VIII. 2014, A&A, 571, A8  \nPlanck Collaboration XI. 2014, A&A, 571, A11  \nPlanck Collaboration Int. XIX. 2014, A&A, in press  \nPlanck Collaboration Int. XX. 2014, A&A, in press  \nPlanck Collaboration Int. XXI. 2014, A&A, in press  \nPlanck Collaboration Int. XXXII. 2014, A&A, submitted  \n- Planck Collaboration Int. XXXIII. 2014, A&A, submitted\n- Poppel, W. 1997, Fund. Cosmic Phys., 18, 1  \nPrato, L., Rice, E. L., & Dame, T. M. 2008, Where are all the Young Stars in  \nAquila?, ed. B. Reipurth, 18  \nReipurth, B. 2008, Handbook of Star Forming Regions, Volume I: The Northern Sky  \n- Santos, F. P., Franco, G. A. P., Roman-Lopes, A., Reis, W., & Rom'an-Z'u˜niga, C. G. 2014, ApJ, 783, 1\n- Santos-Lima, R., de Gouveia Dal Pino, E. M., & Lazarian, A. 2012, ApJ, 747, 21\n- Schlafly, E. F., Green, G., Finkbeiner, D. P., et al. 2014, ApJ, 786, 29 Serkowski, K. 1958, Acta Astron., 8, 135\n- Soler, J. D., Hennebelle, P., Martin, P. G., et al. 2013, ApJ, 774, 128\n- Spitzer, Jr., L. 1968, Interscience Tracts on Physics and Astronomy, 28 Sridhar, S. & Goldreich, P. 1994, ApJ, 432, 612\n- Straiˇzys, V., ˇ Cernis, K., & Bartaˇsi¯ut˙e, S. 2003, A&A, 405, 585  \nSugitani, K., Nakamura, F., Watanabe, M., et al. 2011, ApJ, 734, 63  \nTassis, K., Dowell, C. D., Hildebrand, R. H., Kirby, L., & Vaillancourt, J. E. 2009, MNRAS, 399, 1681  \nWard-Thompson, D., Kirk, J. M., Crutcher, R. M., et al. 2000, ApJ, 537, L135 Whittet, D. C. B., Hough, J. H., Lazarian, A., & Hoang, T. 2008, ApJ, 674, 304 Whittet, D. C. B., Prusti, T., Franco, G. A. P., et al. 1997, A&A, 327, 1194 Wilking, B. A., Gagn'e, M., & Allen, L. E. 2008, Star Formation in the GLYPH<26> Ophiuchi Molecular Cloud, ed. B. Reipurth, 351  \nZweibel, E. G. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 97, Polarimetry of the Interstellar Medium, ed. W. G. Roberge &  \nD. C. B. Whittet, 486", 'Appendix A: Selection of data': 'The HRO analysis is applied to each MC using common criteria for selecting the areas in which the relative orientation is to be assessed.', 'A.1. Gradient mask': 'The dust optical depth, GLYPH<28> 353, observed in each region, can be interpreted as  \nGLYPH<28> obs 353 = GLYPH<28> mc 353 + GLYPH<28> bg 353 + GLYPH<14>GLYPH<28> 353 ; (A.1)  \nwhere GLYPH<28> mc 353 is the optical depth of the MC, GLYPH<28> bg 353 is the optical depth of the di GLYPH<11> use regions behind and / or in front of the cloud (background / foreground), and GLYPH<14>GLYPH<28> 353 the noise in the optical depth map with variance GLYPH<27> 2 GLYPH<28> 353 = GLYPH<14> 2 GLYPH<28> 353 . The gradient of the optical depth can be then written as  \nr GLYPH<28> obs 353 = r GLYPH<28> mc 353 + r ( GLYPH<28> bg 353 + GLYPH<14>GLYPH<28> 353 ) : (A.2)  \nWe quantify the contribution of the background / foreground and the noise, r ( GLYPH<28> bg 353 + GLYPH<14>GLYPH<28> 353 ), by evaluating r GLYPH<28> 353 in a reference field with lower submillimetre emission. Given that the dominant contribution to the background / foreground gradient would come from the gradient in emission from the Galactic plane, for each of the regions analysed we chose a reference field of the same size at the same Galactic latitude and with the lowest average N H in the corresponding latitude band. We compute the average of the gradient norm in the reference field, D jr GLYPH<28> ref 353 j E , and use this value as a threshold for selecting the regions of the map where r GLYPH<28> 353 carries significant information about the structure of the cloud. We note that this threshold includes a contribution from the noise r GLYPH<14>GLYPH<28> 353 . The HROs presented in this study correspond to regions in each field where jr GLYPH<28> 353 j > D jr GLYPH<28> ref 353 j E .', 'A.2. Polarization mask': "The total Stokes parameters Q and U measured in each region can be interpreted as  \nQ obs = Q mc + Q bg + GLYPH<14> Q ; U obs = U mc + U bg + GLYPH<14> U ; (A.3)  \nwhere Q mc and U mc correspond to the polarized emission from the MC, Q bg and U bg correspond to the polarized emission from the di GLYPH<11> use background / foreground, and GLYPH<14> Q and GLYPH<14> U are the noise contributions to the observations, such that the variances GLYPH<27> 2 Q = GLYPH<14> 2 Q and GLYPH<27> 2 U = GLYPH<14> 2 U .  \nAs in the treatment of the gradient, we estimate the contributions of the background / foreground polarized emission and the noise using the rms of the Stokes parameters in the same reference field, Q ref rms and U ref rms . The HROs presented in this study correspond to pixels in each region where j Q j > 2 Q ref rms or j U j > 2 U ref rms . The 'or' conditional avoids biasing the selected values of polarization. This first selection criterion provides a similar sample to the alternative coordinate-independent criterion p Q 2 + U 2 > 2 p ( Q ref rms ) 2 + ( U ref rms ) 2 . This first criterion aims to distinguish between the polarized emission coming from the cloud and the polarized emission coming from the background / foreground estimated in the reference regions.  \nAdditionally, as a second criterion, our sample is restricted to polarization measurements where j Q j > 3 GLYPH<27> Q or j U j > 3 GLYPH<27> U . This aims to select pixels where the uncertainty in the polarization angle is smaller than the size of the angle bins used for the constructions of the HRO (Serkowski 1958; Montier et al. 2015). In terms of the total polarized intensity, P = p Q 2 + U 2 , and following equations. B.4 and B.5 in Planck Collaboration Int. XIX (2014), the second criterion corresponds to P =GLYPH<27> P > 3 and uncertainties in the polarized orientation angle GLYPH<27> < 10 GLYPH<14> .  \nThe fractions of pixels considered in each region, after applying the selection criteria described above, are summarized in Table A.1. The largest masked portions of the regions correspond to the gradient mask, which selects mostly those areas of column density above the mean column density of the background / foreground D N bg H E . The polarization mask provides an independent criterion that is less restrictive. The intersection of these two masks selects the fraction of pixels considered for the HRO analysis.", 'Appendix B: Construction of the histogram of relative orientations and related uncertainties': 'The HROs were calculated for 25 column density bins having equal numbers of selected pixels (10 bins in two regions with fewer pixels, CrA and IC 5146). For each of these HROs we use 12 angle bins of width 15 GLYPH<14> (see Sect. 4.1.1).', 'B.1. Calculation of r GLYPH<28> 353 and the uncertainty of its orientation': 'The optical depth gradient ( r GLYPH<28> 353) is calculated by convolving the GLYPH<28> 353 map with a Gaussian derivative kernel (Soler et al. 2013), such that r GLYPH<28> 353 corresponds to  \nr GLYPH<28> 353 = ( Gx GLYPH<10> GLYPH<28> 353) ˆ i + ( Gy GLYPH<10> GLYPH<28> 353) ˆ j = gx ˆ i + gy ˆ j ; (B.1)  \nwhere Gx and Gy are the kernels calculated using the x - and y -derivatives of a symmetric two-dimensional Gaussian function. The orientation of the isoGLYPH<28> 353 contour is calculated from the components of the gradient vector,  \nGLYPH<18> = arctan GLYPH<16> GLYPH<0> gx ; gy GLYPH<17> : (B.2)  \nBecause the calculation of the gradient through convolution is a linear operation, the associated uncertainties can be calculated using the same operation,so that  \nGLYPH<14> gx = Gx GLYPH<10> GLYPH<14>GLYPH<28> 353 ; GLYPH<14> gy = Gy GLYPH<10> GLYPH<14>GLYPH<28> 353 ; (B.3)  \nTable A.1. Selection of data a .  \nfrom which we obtain the GLYPH<27> 2 gx = GLYPH<14> 2 gx and GLYPH<27> 2 gy = GLYPH<14> 2 gy .  \nThe standard deviation of the angle GLYPH<18>  \ncan be written as  \nGLYPH<27>GLYPH<18> = s @GLYPH<18> @ gx ! 2 GLYPH<27> 2 gx + @GLYPH<18> @ gy ! 2 GLYPH<27> 2 gy ; (B.4)  \nwhich corresponds to  \nGLYPH<27>GLYPH<18> = 1 g 2 x + g 2 y q g 2 y GLYPH<27> 2 gx + g 2 x GLYPH<27> 2 gy : (B.5)  \nIn the application discussed here, the standard deviations in the GLYPH<28> 353 map within the selected areas are much less than a few percentage points, so their e GLYPH<11> ect on the estimate of the orientation of the gradient is negligible.', 'B.2.1. Uncertainties in the construction of the histogram': 'To estimate the uncertainty associated with the noise in Stokes Q and U , we produce 1000 noise realizations, Qr and Ur , using Monte Carlo sampling. We assume that the errors are normally distributed and are centred on the measured values Q and U with dispersions GLYPH<27> Q and GLYPH<27> U (Planck Collaboration II 2014; Planck Collaboration VI 2014; Planck Collaboration V 2014; Planck Collaboration VIII 2014). Given that GLYPH<27> QU is smaller than GLYPH<27> 2 Q and GLYPH<27> 2 U , it is justified to generate Qr and Ur independently of each other. We then introduce Qr and Ur in the analysis pipeline and compute the HRO using the corresponding GLYPH<28> 353 map in each region. The results, presented in Fig. B.1 for the Taurus region, show that the noise in Q and U does not critically a GLYPH<11> ect the shape of the HROs or the trend in GLYPH<24> . Together, the low noise in the maps of GLYPH<28> 353 and the selection criteria for the polarization measurements ensure that the HRO is well determined.  \nFig. B.1. HROs in the Taurus region that correspond to the indicated N H bins. The plotted values are obtained using the original GLYPH<28> 353 map at 10 0 resolution and maps of the Stokes parameters Qr and Ur , which correspond to 1000 random noise realizations. Each realization is generated using a Gaussian probability density function centred on the measured values Q and U with variances GLYPH<27> 2 Q and GLYPH<27> 2 U .  \n<!-- image -->  \nAnother source of uncertainty in the HRO resides in the histogram binning process. The variance in the k th histogram bin is given by  \nGLYPH<27> 2 k = hk 1 GLYPH<0> hk h tot ! ; (B.6)  \nwhere hk is the number of samples in the k th bin, and h tot is the total number of samples.  \nOf the two independent sources of uncertainty above, we find that the largest contribution comes from the binning process, so that these are the ones shown as the shaded uncertainty ranges in all figures of HROs, for example for Taurus in Fig. 3. Because of the large number of samples in each histogram bin, the uncertainties in the HRO do not significantly a GLYPH<11> ect the results of this study.', 'B.2.2. Uncertainties in the histogram shape parameter GLYPH<24>': 'As in the case of the HRO, the uncertainty in GLYPH<24> , as defined in Eq. (4), can be quantified using the random realizations introduced in the previous section. Figure B.2 shows the dependence of GLYPH<24> on log 10 N H obtained using Qr and Ur for the Taurus region. The small variations around the trend line indicate that the uncertainties in Q and U do not significantly a GLYPH<11> ect the trends discussed in this study, as expected from the behaviour of the histograms presented in Fig. B.1. The main source of uncertainty in the estimation of GLYPH<24> is related to the histogram binning, characterized by the error bars calculated using Eq. (5). As seen in Fig. 7 and reproduced in Fig. C.1, these are much larger than the dispersion of the values of GLYPH<24> in Fig. B.2.  \nFig. B.2. Histogram shape parameter, GLYPH<24> , as a function of log 10 ( N H = cm GLYPH<0> 2 ) in the Taurus region. The values are obtained using the GLYPH<28> 353 map at 10 0 resolution and maps of the Stokes parameters Qr and Ur that correspond to 1000 random-noise realizations. Each realization is generated using a Gaussian probability density function centred on the measured values Q and U with variances GLYPH<27> 2 Q and GLYPH<27> 2 U . By joining the values at each N H bin, we find a trend very close to the black line in Fig. 7, with little dispersion from the noise; much larger are the uncertainties in evaluating GLYPH<24> at each N H bin, as given in Fig. 7 but not shown here.  \n<!-- image -->', 'C.1. A product of chance?': 'To investigate various potential sources of the signal found in the HROs we use the GLYPH<28> 353 map at 10 0 resolution in each region in combination with Q and U maps that are produced with di GLYPH<11> erent recipes, each with 1000 realizations.  \nTo eliminate random fields, we use Q and U maps produced with a random realization of with a uniform distribution and with unit-length polarization pseudo-vectors. The results of this numerical experiment are shown in Fig. C.1 for the Taurus region. For each of the 1000 realizations, we join the values at each N H bin to show the trend lines in order to compare with the lines in Fig. 7.  \nFig. C.1. Histogram shape parameter, GLYPH<24> , as a function of log 10 ( N H = cm GLYPH<0> 2 ) in the Taurus region. Blue lines join the GLYPH<24> values obtained using the GLYPH<28> 353 map at 10 0 resolution and each of 1000 random realizations of Q and U maps corresponding to a uniform distribution of . Results in black and grey are from the analysis of the Planck data, as reported in Fig. 7.  \n<!-- image -->  \nTo eliminate the large-scale magnetic field as the source we use Q and U maps produced with random realizations of with  \na Gaussian distribution and unit-length polarization pseudovectors. The polarization angle distribution is centred at 0 = 0 GLYPH<14> with a standard deviation & = 45 GLYPH<14> . The results of this numerical experiment are shown in Fig. C.2 for the Taurus region.  \nFig. C.2. Like Fig. C.1 but for values of GLYPH<24> (red) obtained using 1000 random realizations of Q and U maps corresponding to a Gaussian distribution of centred on 0 = 0 GLYPH<14> and with standard deviation & = 45 GLYPH<14> .  \n<!-- image -->  \nTo eliminate random spatial correlations, we use Q and U maps produced from random realizations with a power spectrum P ( k ) / k GLYPH<11> m . For the spectral indices we adopted, GLYPH<11> m = GLYPH<0> 1 : 5, GLYPH<0> 2 : 5, and GLYPH<0> 3 : 5, correlations are introduced in the orientation of the polarization pseudo-vectors that are independent of the structure of matter. This test evaluates the statistical significance of the random sampling of spatial correlations in the magnetic field implied when we calculate the HROs in a finite region of the sky. The results of this numerical experiment are shown in Fig. C.3 for the Taurus region.  \nThe results of these numerical experiments show that the trends with N H found in the HROs and GLYPH<24> do not arise by chance from these potential sources.', 'C.2. A product of random correlations in the polarization maps?': 'The random realizations of the magnetic field presented above are useful for characterizing the behaviour of GLYPH<24> . However, they are bound to produce very di GLYPH<11> erent statistics compared to those in the real observations. In reality the orientations of B ? have some non-trivial correlation structure in the map. In principle, the di GLYPH<11> erence in the HROs might be due to the di GLYPH<11> erent correlations of the observed B ? and of the random realizations of B ? .  \nTo evaluate whether the signal present in the HROs arises from an actual physical relationship between B ? and GLYPH<28> 353, we introduce randomness by shifting the Stokes Q and U maps with respect to the GLYPH<28> 353 map and then calculate the corresponding HROs and GLYPH<24> as a function of log 10 ( N H = cm GLYPH<0> 2 ). The intrinsic statistical properties of the two maps are unchanged because the two maps are unchanged, only shifted.  \nIf the trend in GLYPH<24> as a function of log 10 ( N H = cm GLYPH<0> 2 ) were arising randomly, these trends would be unchanged even for significant shifts. Instead the results of this experiment, illustrated in Fig. C.4 for the Taurus region, show that the trends tend to disappear with increasing values of the size of the shift. For shifts of about 1 GLYPH<14> , the correlation between the magnetic field and the matter is still present, as expected from the results presented in Fig. 10, but for larger shifts the correlation is lost and the trend  \nFig. C.3. Like Fig. C.1 but for values of GLYPH<24> (green) obtained using 1000 random realizations of Q and U maps corresponding to a power spectrum P ( k ) / k GLYPH<11> m , with GLYPH<11> m = GLYPH<0> 1 : 5, GLYPH<0> 2 : 5, and GLYPH<0> 3 : 5.  \n<!-- image -->  \ndoes not survive. Over the many MCs studied, the nature of the trends at large shifts appears to be random.', 'C.3. Projection effects': 'We evaluate the statistical significance of the relative orientation between B ? and r GLYPH<28> 353 by considering the distribution of relative orientations between two vectors in 3D space compared to the distribution of relative orientation between their projections in 2D.  \nThe projection of a vector v onto a plane normal to the unit vector ˆ n can be written as  \nu = v GLYPH<0> ( v GLYPH<1> ˆ n ) ˆ n : (C.1)  \nConsider two unit vectors in 3D, v 1 and v 2, separated by an angle GLYPH<11> , such that  \ncos GLYPH<11> = v 1 GLYPH<1> v 2 : (C.2)  \nThe angle GLYPH<12> between the projections of these two vectors onto a plane normal to ˆ n , u 1, and u 2, can be written as  \ncos GLYPH<12> = u 1 GLYPH<1> u 2 j u 1 jj u 2 j = ( v 1 GLYPH<1> v 2 GLYPH<0> ( v 1 GLYPH<1> ˆ n )( v 2 GLYPH<1> ˆ n ) ) GLYPH<0> v 1 GLYPH<1> v 1 GLYPH<0> ( v 1 GLYPH<1> ˆ n ) 2 GLYPH<1> 1 = 2 GLYPH<0> v 2 GLYPH<1> v 2 GLYPH<0> ( v 2 GLYPH<1> ˆ n ) 2 GLYPH<1> 1 = 2 : (C.3)  \nFig. C.4. Histogram shape parameter, GLYPH<24> , as a function of log 10 ( N H = cm GLYPH<0> 2 ) calculated from the GLYPH<28> 353 map at 10 0 resolution in the Taurus region and the Stokes Q and U maps shifted in Galactic longitude and latitude by the values indicated, for GLYPH<1> l and GLYPH<1> b , respectively.  \n<!-- image -->  \nGiven a particular distribution of angles between the vectors v 1 and v 2, this expression, which is solved numerically, is useful for evaluating the resulting distribution of angles between the projected vectors u 1 and u 2. Without any loss of generality, we can assume that v 1 is oriented along the axis of a spherical coordinate system, such that v 1 = ˆ k , and that v 2 is oriented at polar angle GLYPH<18> and azimuth GLYPH<30> , such that v 2 = cos GLYPH<30> sin GLYPH<18> ˆ i + sin GLYPH<30> sin GLYPH<18> ˆ j + cos GLYPH<18> ˆ k . Then, we can simulate a distribution of cos GLYPH<11> by simulating a distribution of cos GLYPH<18> and generating GLYPH<30> with a uniform distribution.  \nWe thus generate a set of vectors v 2 that follow a particular distribution of cos GLYPH<11> . We choose three examples: a uniform distribution of relative orientation between v 2 and v 1, v 2 vectors that are mostly parallel to v 1, and v 2 vectors that are mostly perpendicular to v 1. The last two are Gaussian distributions centred at cos GLYPH<11> = 0 (for mostly parallel) or cos GLYPH<11> = GLYPH<6> 1 (for mostly perpendicular, given the periodicity of GLYPH<11> ) with a dispersion &GLYPH<11> .  \nThese distributions are shown in the top panel of Fig. C.5 for three values of &GLYPH<11> .  \nFig. C.5. Normalized distributions of relative orientations of random vectors in 3D (top) and of their projections in 2D (bottom). The curves correspond to mostly perpendicular vectors (red), uniform distribution of relative orientations (black), and mostly parallel vectors (green). The solid, dashed, and dotted lines correspond to dispersions of the relative orientation angle &GLYPH<11> = 18 GLYPH<14> , 36 GLYPH<14> , and 72 GLYPH<14> , respectively.  \n<!-- image -->  \nUsing Eq. (C.3) we calculate the distribution of the angles between the projected vectors and the results are shown in the bottom panel of Fig. C.5. Projection generally tends to make vector pairs look more parallel in 2D than they are in 3D, regardless of the orientation in 3D, in agreement with the results of the Gaussian models presented in Planck Collaboration Int. XXXII (2014). Nevertheless, the distribution of relative orientations in 2D is qualitatively similar (though not identical) to the distribution in 3D. In particular, the signature of perpendicular orientation is not erased. We conclude that the observation in 2D of B ? close to perpendicular to the column density structures is a direct indicator of the perpendicular configuration of B with respect to the density structures in 3D. In contrast, the parallel alignment of B ? with the structure, or no signs of preferential orientation, while suggestive of parallel orientation in 3D, does not unambiguously rule out the presence of some close-to-perpendicular orientations within the distribution in 3D.', 'Appendix D: Alternative estimates of magnetic field strength': 'Given the historic importance of the DCF method (Davis 1951; Chandrasekhar & Fermi 1953) and the related DCF + SF method  \n(Hildebrand et al. 2009), we have used them to estimate the magnetic field strength in each region. We also evaluated the mass-to-flux ratio, which is important in investigating the stability against gravitational collapse. We now provide a critical assessment of the applicability of these results.', 'D.1. Davis-Chandrasekhar-Fermi method': "The DCF method estimates the strength of B ? in a region using the dispersion of the polarization angles & . Assuming that the magnetic field is frozen into the gas and that the dispersion of the B ? orientation angles (or equivalently & ) is due to transverse incompressible Alfv'en waves, then  \nB DCF ? = p 4 GLYPH<25> GLYPH<26> GLYPH<27> v k & ; (D.1)  \nwhere GLYPH<27> v k is the dispersion of the radial velocity of the gas (Appendix D.3) and GLYPH<26> the gas mass density.  \nWe calculate & directly from Stokes Q and U using  \n& = q GLYPH<10> ( GLYPH<1> ) 2 GLYPH<11> (D.2)  \nand  \nGLYPH<1> = 1 2 arctan ( Q h U i GLYPH<0> h Q i U ; Q h Q i + h U i U ) ; (D.3)  \nwhere h : : : i denotes an average over the selected pixels in each map (Planck Collaboration Int. XIX 2014).", 'D.2. Davis-Chandrasekhar-Fermi plus structure function method': "As described by Hildebrand et al. (2009), the DCF + SF method characterizes the magnetic field dispersion about local structured fields by considering the di GLYPH<11> erence in angle, GLYPH<1> ( ' ) = ( x ) GLYPH<0> ( x + ' ), between pairs of B ? vectors separated by displacements ' in the plane of the sky. Assuming that the angle di GLYPH<11> erences are statistically isotropic (i.e., they depend only on ' = j ' j and not on the orientation of ' ), they can binned by distance, ' . From the N ( ' ) pairs of B ? vectors for that bin, the square of the second-order structure function is  \nS 2 2 ( ' ) = D [ GLYPH<1> ( x ; ' )] 2 E x = * 1 N ( ' ) N ( ' ) X i = 1 ( GLYPH<1> x ; i ) 2 + x ; (D.4)  \nas introduced by Kobulnicky et al. (1994) and FalcetaGonc¸alves et al. (2008). In terms of the Stokes parameters, each term in the sum can be written  \nGLYPH<1> x ; i = 1 2 arctan ( QiUx GLYPH<0> QxUi ; QiQx + UiUx ) ; (D.5)  \nwhere the subscripts x and i represent the central and displaced positions, respectively.  \nHildebrand et al. (2009) assume that B ( x ) is composed of a large-scale structured field, B 0( x ), and a random component, B r( x ), which are statistically independent of each other. Assuming that B 0( x ) is a smoothly varying quantity, its contribution to S 2 2 ( ' ) should increase in proportion to ' 2 for distances that are much smaller than the scales at which B 0 itself fluctuates. Additionally, assuming that turbulence occurs on scales that are small enough to completely decorrelate B r for the range  \nTable D.1. Magnetic properties of the selected regions. a  \nof scales probed by the displacements ' , Hildebrand et al. (2009) derived the approximation  \nS 2 2 ( ' ) = b 2 + m 2 ' 2 ; (D.6)  \nwhere the two terms on the right-hand side give the contributions from the random and large-scale magnetic fields, respectively. These assumptions are not necessarily valid for the range of scales in the Planck data.  \nHildebrand et al. (2009) used Eq. (D.6) to estimate b 2 as the intercept and then related b to the ratio of the random to the large-scale magnetic field strength, both projected onto the plane of the sky, through  \nD B 2 r ; ? E 1 = 2 B 0 ; ? = b p 2 GLYPH<0> b 2 ; (D.7)  \nwhere D B 2 r ; ? E 1 = 2 stands for the root mean square (rms) variations about the large-scale magnetic field, B 0 ; ? . The same assumptions that result in Eq. (D.1) - i.e., considering only incompressible and isotropic turbulence, magnetic fields frozen into the gas, and dispersion of the B ? orientation originating in transverse incompressible Alfv'en waves - lead to  \nD B 2 r ; ? E 1 = 2 B 0 ; ? = GLYPH<27> v k V A ; ? ; (D.8)  \nwhere  \nV A ; ? = B 0 ; ? p 4 GLYPH<25>GLYPH<26> (D.9)  \nis the speed of the transverse incompressible Alfv'en waves. Combining Eqs. (D.7) - (D.9) results in  \nB DCF + SF ? GLYPH<17> B DCF + SF 0 ; ? = p 4 GLYPH<25> GLYPH<26> GLYPH<27> v k p 2 GLYPH<0> b 2 b : (D.10)", 'D.3. Calculation': "The estimates of velocity dispersion GLYPH<27> v k and mass density GLYPH<26> are the same in the two methods (DCF and DCF + SF).  \nWeobtain GLYPH<27> v k from the most complete CO emission-line survey of the Milky Way currently available, that of Dame et al.  \n(2001). This data set consists of 488 000 spectra that beamsample (1 / 8 GLYPH<14> ) a set of MCs and the Galactic plane over a strip from 4 GLYPH<14> to 10 GLYPH<14> wide in latitude. We find that material in the MCs has centroid velocities in the range GLYPH<0> 10 < v k = (kms GLYPH<0> 1 ) < 10 in the first moment map. Furthermore, for the calculations below, we need to select pixels with su GLYPH<14> cient S / N in the polarization observations and so use the second criterion described in Appendix A.2. This combination of cuts selects the areas illustrated in Fig. D.1. Over these areas we calculate the average of the velocity dispersion from the second moment map and use this value as GLYPH<27> v k . This is tabulated along with other properties in Table D.1.  \nThe mass density, GLYPH<26> = GLYPH<22> n m p, is the product of the proton mass, m p, the mean molecular weight per hydrogen molecule, GLYPH<22> = 2 : 8, and the mean number density, n . We require a value of n , which for the discussions here we approximate to be 100 cm GLYPH<0> 3 . This is a typical value for MCs (Draine 2011) and in rough agreement with the column densities and cloud sizes presented in Table 1. In practice, n varies from one cloud to the next, and its estimation involves the assumption of a particular cloud geometry, resulting in additional uncertainties that are not considered in this study.  \nFor the DCF method, the dispersion of the orientation angle & corresponds to the centred second moment of the distribution evaluated in pixels within the selected region. Following the analysis presented in Planck Collaboration Int. XIX (2014) for S ( ' ), the variance on & can be expressed as  \nGLYPH<27> 2 & = 1 N 2 & 2 2 6 6 6 6 6 6 6 4 0 B B B B B @ N X i = 1 & i 1 C C C C C A 2 GLYPH<27> 2 + N X i = 1 & 2 i GLYPH<27> 2 i 3 7 7 7 7 7 7 7 5 ; (D.11)  \nwhere i and GLYPH<27> i are the orientation angle and uncertainty for each polarization pseudo-vector, N is the number of measurements, & i = i GLYPH<0> h i , and h i is the average orientation angle in each region. Like other quadratic functions, GLYPH<27> 2 & is biased positively when noise is present, leading to an overestimation of this quantity. However, given that we limit our analysis to polarization measurements with high S / N and that the uncertainties in GLYPH<27> v k are considerably larger, the bias correction does not have a significant e GLYPH<11> ect on our estimates. We apply Eq. (D.1) to estimate B DCF ? and propagate the errors to obtain the values listed in Table D.1.  \nFig. D.1. Line-of-sight centroid velocity v k maps inferred from CO emission (Dame et al. 2001). Areas shown have su GLYPH<14> cient S / N in the polarization observations (second polarization criterion in Appendix A.2) for the DCF and DCF + SF analyses and lie in the range GLYPH<0> 10 < v k = (kms GLYPH<0> 1 ) < 10.  \n<!-- image -->  \nFig. D.2. Structure function S 2( ' ) calculated from the Q and U maps of the selected regions following Eq. (D.4). The black lines indicate the fits given by Eq. (D.6). The vertical dashed line marks the common 10 0 resolution of the data used in the analysis; the data are correlated for low values of ' causing the drop in S 2( ' ).  \n<!-- image -->  \nFor the DCF + SF method, to estimate the parameter b , we first calculate S 2 2 ( ' ) for each of the regions. The resolution of the data used is 10 0 . We evaluate S 2 2 ( ' ) in steps of 3. 0 44 for lags 0 0 < ' < 34. 0 4 and in steps of 34. 0 4 for lags 40 0 < ' < 200 0 . For each considered lag ' , GLYPH<1> x ; i is evaluated pixel by pixel, considering all the pixels located in an annulus with radius ' . The term GLYPH<1> x ; i is only considered in the calculation of S 2 2 ( ' ) if there are at least three pixels in the annulus. As anticipated, the model from Eq. (D.6) does not agree with the data on all sampled scales, so the range of scales considered is limited to ' above the resolu-  \ntion of the data and between 50 0 and 200 0 , where the behaviour of S 2 2 ( ' ) is approximately linear in ' 2 . We then make a linear fit of S 2 2 ( ' ) using ' 2 as the independent variable and estimate the value of b and its uncertainty, GLYPH<27> b . The results are plotted in Fig. D.2.  \nUsing the calculated b values, 6 we apply Eq. (D.10) to estimate B DCF + SF ? , and the results are listed in Table D.1. As ex-  \npected, the DCF + SF method produces estimates of B ? that are about twice as large as those obtained with the DCF method. However, it is worth noting that because of the discrepancy between the model and the observations, the value of b depends on the selected ' -range and so propagates into a di GLYPH<11> erent value of B DCF + SF ? .", 'D.4. Mass-to-flux ratios': 'The critical value for the mass that can be supported against gravity by a magnetic flux GLYPH<8> can be estimated to first order for a uniform disk from ( M = GLYPH<8> )crit GLYPH<17> 1 = (2 GLYPH<25> G 1 = 2 ) (Nakano & Nakamura 1978). The precise value of the right-hand side changes for di GLYPH<11> erent geometries (e.g., Spitzer 1968; McKee et al. 1993). Stability can be assessed using  \nGLYPH<21> = ( M = GLYPH<8> ) ( M = GLYPH<8> )crit = 7 : 6 GLYPH<2> 10 GLYPH<0> 21 ( N mft H2 = cm GLYPH<0> 2 ) ( B mft =GLYPH<22> G) ; (D.12)  \nwhere N mft H2 and B mft are the H2 column density and magnetic field strength along a magnetic flux tube (Crutcher et al. 2004). A cloud is supercritical and prone to collapse under its own gravity, when GLYPH<21> > 1; otherwise, when GLYPH<21> < 1, the cloud is sub-critical, magnetically supported against gravitational collapse.  \nWhat is observable is GLYPH<21> obs, in which N mft H2 = B mft is replaced by N H2 = B ? . We evaluate GLYPH<21> obs by combining the value of GLYPH<10> N H2 GLYPH<11> computed from the integrated CO line emission and the conversion factor X CO = (1 : 8 GLYPH<6> 0 : 3) GLYPH<2> 10 20 cm GLYPH<0> 2 K GLYPH<0> 1 km GLYPH<0> 1 s (Dame et al. 2001) and B ? estimated with the DCF and DCF + SF methods. The X CO factor may show cloud-to-cloud or regional variations (Draine 2011), but we consider that these are not significant in comparison to the uncertainties involved in the estimation of B ? . The calculated values of GLYPH<21> obs are listed in Table D.1. They are consistent with being less than unity.  \nFrom this we might obtain GLYPH<21> by judicious deprojection, since  \nGLYPH<21> = ( N mft H2 = N H2) ( B ? = B mft ) GLYPH<2> GLYPH<21> obs GLYPH<17> f dp GLYPH<21> obs : (D.13)  \nHere, B ? is always less than B mft , pushing f dp below unity. The situation for column density depends on the geometry of the structure relative to the magnetic field. For a structure with the magnetic field along the short axis, N mft H2 < N H2, again lowering f dp.  \nStatistically, the mean mass-to-flux ratio can be related to the observed value by assuming a particular geometry of the cloud, an ellipsoid with equatorial radius a and centre-to-pole distance c , and a magnetic field oriented along the polar axis of the ellipsoid. For an oblate spheroid, flattened perpendicular to the orientation of the magnetic field, f dp = 1 = 3, yielding  \nM = GLYPH<8> = Z GLYPH<25>= 2 0 M cos GLYPH<12> GLYPH<8> = sin GLYPH<12> sin GLYPH<12> d GLYPH<12> = 1 3 ( M = GLYPH<8> )obs ; (D.14)  \nwhere GLYPH<12> is the inclination angle with respect to the line of sight, and the sin GLYPH<12> dependence in the flux comes from B = B ? sin GLYPH<12> , (Crutcher et al. 2004). For a sphere there is no cos GLYPH<12> dependence and f dp = 1 = 2. For a prolate spheroid elongated along the orientation of the field, the mass is multiplied by sin GLYPH<12> instead of cos GLYPH<12> , resulting in f dp = 3 = 4. Investigating which geometry, if any, is most relevant to the actual MCs and the magnetized structures detected within them is obviously important before firm conclusions can be drawn regarding gravitational instability.', 'D.5. Discussion': "Our estimates of the magnetic field strengths in the MCs analysed and the mass-to-flux ratios, presented in Table D.1, stem from the classic calculation presented by Chandrasekhar & Fermi (1953) and an updated interpretation presented by Hildebrand et al. (2009). Both of these methods assume very specific conditions, so (as we discuss below) this limits the conclusions that can be drawn from these estimates of the magnetic field strength. The deduced mass-to-flux ratios suggest that the clouds are potentially magnetically sub-critical, but we again need to be cautious about drawing conclusions from this application of the DCF and DCF + SF analyses alone.  \nIn the case of the DCF method, the values of B DCF ? are obtained by assuming that the structure of the magnetic field is the product only of incompressible Alfv'en waves, where the displacements are perpendicular to the direction of propagation. This is not the case for turbulence in MCs where the random component of the magnetic field can have any orientation. The dispersion measured about mean fields, assumed to be straight in DCF, may be much larger than should be attributed to MHD waves or turbulence, leading to an overestimation of & and to low values of B ? . Moreover true interstellar turbulence in MCs involves not only incompressible Alfv'en waves, but also compressible magneto-sonic waves, which do not satisfy Eq. D.1. Furthermore, depending on the scales examined, the magnetic field may have structures due to e GLYPH<11> ects such as di GLYPH<11> erential rotation, gravitational collapse, or expanding H ii regions.  \nIn the case of the DCF + SF method, the values of B DCF + SF ? are obtained by assuming a very specific model of the magnetic field , which is in principle just a first-order approximation. First, this model assumes that the e GLYPH<11> ect of the large-scale structured magnetic field, B 0, is to cause the square of the second-order structure function, S 2 2 ( ' ), to increase as ' 2 . This corresponds to a very specific correlation function for B 0. Second, this model assumes that the dispersion of the random component of the field, B r, is scale-independent, which is not realistic for the range of scales probed by Planck (Elmegreen & Scalo 2004; Hennebelle &Falgarone 2012).  \nIn addition, the magnetic field orientation deduced from the polarization angle in a particular direction is not generally that of the field at a single point along the line of sight. The observed polarization is the average of various field pseudo-vectors weighted by local dust emission along the line of sight. The net e GLYPH<11> ect of the integration of multiple uncorrelated components along the line of sight is an observed dispersion of the polarization angle that is smaller than the true 3D dispersion of the magnetic field orientation, thus leading to an overestimation of B ? . Myers & Goodman (1991) presented an analysis of this e GLYPH<11> ect in terms of the number of correlation lengths of the magnetic field along the line of sight through a cloud, which they calculated empirically.  \nHoude et al. (2009) presented an extension of the DCF + SF method that includes the e GLYPH<11> ect of signal integration along the line of sight and across the area subtended by the telescope beam. The extended method, also implemented in Houde et al. (2011, 2013), is based on the identification of the magnetized turbulence correlation length ( GLYPH<14> ) by means of the structure function of the polarization angles. In the case of the Planck 353 GHz observations, the angular resolution is not su GLYPH<14> cient to identify GLYPH<14> and the corrections would have to rely on rough estimates of this value and the depth of integration ( GLYPH<1> ). Following equation 29 in Houde et al. (2009), rough estimates GLYPH<14> GLYPH<25> 0 : 2 pc and GLYPH<1> GLYPH<25> 10 pc result in a few correlation lengths across the beam, corresponding to correction factors around 0.4. Coincidentally, such cor-  \nrection factors lead to values of B DCF + SF ? close to those of B DCF ? in Table D.1. But note that this correction relies on specific assumptions on the nature of the turbulence correlation function, it does not circumvent the necessity to observe the dust polarized emission with higher angular resolution to fully characterize the magnetic field.  \nMHDsimulations provide a potentially useful guide to what modifications might be introduced into the DCF formula to allow for inhomogeneity and line-of-sight averaging. Using synthetic observations of MHD simulations, Ostriker et al. (2001) showed that correcting Eq. (D.1) by a factor C GLYPH<25> 0 : 5 provides a good approximation to the actual magnetic field strength in cases where &  < 25 GLYPH<14> . However, the e GLYPH<11> ect of nonlinear amplitudes is uncertain (Zweibel 1996), and the method fails for values of &  > 25 GLYPH<14> , which is the case for all of the regions in this study (Table D.1). Falceta-Gonc¸alves et al. (2008) propose a method that is potentially valid for any value of & , based on a fit to the B DCF ? values obtained from maps at di GLYPH<11> erent resolutions, again concluding that the field should be lower than estimated with Eq. (D.1). However, this approach was not tested in MHD simulations that include gravity, which is the critical process that we aim to evaluate by using the DCF method.  \nAnother strong assumption is that the behaviour of the velocity and the magnetic field are represented well by the observed quantities GLYPH<27> v k and & for a particular set of scales, which might not necessarily be the case. Even if the power spectra of v and B are comparable in 3D, the integration along the line of sight is di GLYPH<11> erent for the two quantities. The dispersion GLYPH<27> v k is based on the emission-line profile v k tracing a gas species, and while the emission is directly integrated along the line of sight, the line profile is possibly a GLYPH<11> ected by radiative transfer and excitation e GLYPH<11> ects. The tracer of the magnetic field is the optically-thin polarized submillimetre emission of dust, and both the polarization and B ? are projected and integrated along the line of sight as pseudo-vectors.  \nFinally, in the estimates of B DCF ? and B DCF + SF ? , a common mean density GLYPH<26> has been adopted, while in practice GLYPH<26> is di GLYPH<11> erent from cloud to cloud. Direct estimation of the values of GLYPH<26> from the measured column densities h N H i and GLYPH<10> N H2 GLYPH<11> relies heavily on the geometrical modelling of the cloud and the filling factors of each species, introducing uncertainties that a GLYPH<11> ect the calculated values of B ? .  \nGiven the line-of-sight integration and the fact that &  > 25 GLYPH<14> uniformly, the calculated values could be considered as upper limits to the actual B ? . However, other shortcomings, such as the assumptions about the correlation structure and the uncertainties in the determination of the density, do not necessarily bias the estimate of the magnetic field strength towards high values. In conclusion, the values presented in Table D.1 should be viewed only as a reference and only applied with caution, given the many assumptions in both methods at the scales considered.  \n- 1 APC, AstroParticule et Cosmologie, Universit'e Paris Diderot, CNRS / IN2P3, CEA / lrfu, Observatoire de Paris, Sorbonne Paris Cit'e, 10, rue Alice Domon et L'eonie Duquet, 75205 Paris Cedex 13, France\n- 2 African Institute for Mathematical Sciences, 6-8 Melrose Road, Muizenberg, Cape Town, South Africa\n- 3 Agenzia Spaziale Italiana Science Data Center, Via del Politecnico snc, 00133, Roma, Italy\n- 4 Aix Marseille Universit'e, CNRS, LAM (Laboratoire d'Astrophysique de Marseille) UMR 7326, 13388, Marseille, France\n- 5 Astrophysics Group, Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, U.K.\n- 6 Astrophysics & Cosmology Research Unit, School of Mathematics, Statistics & Computer Science, University of KwaZulu-Natal, Westville Campus, Private Bag X54001, Durban 4000, South Africa\n- 7 CITA, University of Toronto, 60 St. George St., Toronto, ON M5S 3H8, Canada\n- 8 CNRS, IRAP, 9 Av. colonel Roche, BP 44346, F-31028 Toulouse cedex 4, France\n- 9 California Institute of Technology, Pasadena, California, U.S.A.\n- 10 Centro de Estudios de F'ısica del Cosmos de Arag'on (CEFCA), Plaza San Juan, 1, planta 2, E-44001, Teruel, Spain\n- 11 Computational Cosmology Center, Lawrence Berkeley National Laboratory, Berkeley, California, U.S.A.\n- 12 DSM / Irfu / SPP, CEA-Saclay, F-91191 Gif-sur-Yvette Cedex, France\n- 13 DTU Space, National Space Institute, Technical University of Denmark, Elektrovej 327, DK-2800 Kgs. Lyngby, Denmark\n- 14 D'epartement de Physique Th'eorique, Universit'e de Gen'eve, 24, Quai E. Ansermet,1211 Gen'eve 4, Switzerland\n- 15 Departamento de Astrof'ısica, Universidad de La Laguna (ULL), E38206 La Laguna, Tenerife, Spain\n- 16 Departamento de F'ısica, Universidad de Oviedo, Avda. Calvo Sotelo s / n, Oviedo, Spain\n- 17 Department of Astronomy and Astrophysics, University of Toronto, 50 Saint George Street, Toronto, Ontario, Canada\n- 18 Department of Astrophysics / IMAPP, Radboud University Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, The Netherlands\n- 19 Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, British Columbia, Canada\n- 20 Department of Physics and Astronomy, Dana and David Dornsife College of Letter, Arts and Sciences, University of Southern California, Los Angeles, CA 90089, U.S.A.\n- 21 Department of Physics and Astronomy, University College London, London WC1E 6BT, U.K.\n- 22 Department of Physics, Florida State University, Keen Physics Building, 77 Chieftan Way, Tallahassee, Florida, U.S.A.\n- 23 Department of Physics, Gustaf Hallstromin katu 2a, University of Helsinki, Helsinki, Finland\n- 24 Department of Physics, Princeton University, Princeton, New Jersey, U.S.A.\n- 25 Department of Physics, University of California, Santa Barbara, California, U.S.A.\n- 26 Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois, U.S.A.\n- 27 Dipartimento di Fisica e Astronomia G. Galilei, Universit'a degli Studi di Padova, via Marzolo 8, 35131 Padova, Italy\n- 28 Dipartimento di Fisica e Scienze della Terra, Universit'a di Ferrara, Via Saragat 1, 44122 Ferrara, Italy\n- 29 Dipartimento di Fisica, Universit'a La Sapienza, P. le A. Moro 2, Roma, Italy\n- 30 Dipartimento di Fisica, Universit'a degli Studi di Milano, Via Celoria, 16, Milano, Italy\n- 31 Dipartimento di Fisica, Universit'a degli Studi di Trieste, via A. Valerio 2, Trieste, Italy\n- 32 Dipartimento di Matematica, Universit'a di Roma Tor Vergata, Via della Ricerca Scientifica, 1, Roma, Italy\n- 33 Discovery Center, Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark\n- 34 Escola de Artes, Ciˆencias e Humanidades, Universidade de S˜ao Paulo, Rua Arlindo Bettio 1000, CEP 03828-000, S˜ao Paulo, Brazil\n- 35 European Space Agency, ESAC, Planck Science O GLYPH<14> ce, Camino bajo del Castillo, s / n, Urbanizaci'on Villafranca del Castillo, Villanueva de la Ca˜nada, Madrid, Spain\n- 36 European Space Agency, ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands\n- 37 Facolt'a di Ingegneria, Universit'a degli Studi e-Campus, Via Isimbardi 10, Novedrate (CO), 22060, Italy\n- 38 Gran Sasso Science Institute, INFN, viale F. Crispi 7, 67100 L'Aquila, Italy\n- 39 HGSFP and University of Heidelberg, Theoretical Physics Department, Philosophenweg 16, 69120, Heidelberg, Germany  \n- 40 Helsinki Institute of Physics, Gustaf Hallstromin katu 2, University of Helsinki, Helsinki, Finland\n- 41 INAF - Osservatorio Astrofisico di Catania, Via S. Sofia 78, Catania, Italy\n- 42 INAF -Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, Padova, Italy\n- 43 INAF - Osservatorio Astronomico di Roma, via di Frascati 33, Monte Porzio Catone, Italy\n- 44 INAF - Osservatorio Astronomico di Trieste, Via G.B. Tiepolo 11, Trieste, Italy\n- 45 INAF / IASF Bologna, Via Gobetti 101, Bologna, Italy\n- 46 INAF / IASF Milano, Via E. Bassini 15, Milano, Italy\n- 47 INFN, Sezione di Bologna, Via Irnerio 46, I-40126, Bologna, Italy\n- 48 INFN, Sezione di Roma 1, Universit'a di Roma Sapienza, Piazzale Aldo Moro 2, 00185, Roma, Italy\n- 49 INFN, Sezione di Roma 2, Universit'a di Roma Tor Vergata, Via della Ricerca Scientifica, 1, Roma, Italy\n- 50 INFN / National Institute for Nuclear Physics, Via Valerio 2, I-34127 Trieste, Italy\n- 51 IPAG: Institut de Plan'etologie et d'Astrophysique de Grenoble, Universit'e Grenoble Alpes, IPAG, F-38000 Grenoble, France, CNRS, IPAG, F-38000 Grenoble, France\n- 52 Imperial College London, Astrophysics group, Blackett Laboratory, Prince Consort Road, London, SW7 2AZ, U.K.\n- 53 Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, U.S.A.\n- 54 Institut N'eel, CNRS, Universit'e Joseph Fourier Grenoble I, 25 rue des Martyrs, Grenoble, France\n- 55 Institut Universitaire de France, 103, bd Saint-Michel, 75005, Paris, France\n- 56 Institut d'Astrophysique Spatiale, CNRS (UMR8617) Universit'e Paris-Sud 11, Bˆatiment 121, Orsay, France\n- 57 Institut d'Astrophysique de Paris, CNRS (UMR7095), 98 bis Boulevard Arago, F-75014, Paris, France\n- 58 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, U.K.\n- 59 Institute of Theoretical Astrophysics, University of Oslo, Blindern, Oslo, Norway\n- 60 Instituto de Astrof'ısica de Canarias, C / V'ıa L'actea s / n, La Laguna, Tenerife, Spain\n- 61 Instituto de F'ısica de Cantabria (CSIC-Universidad de Cantabria), Avda. de los Castros s / n, Santander, Spain\n- 62 Istituto Nazionale di Fisica Nucleare, Sezione di Padova, via Marzolo 8, I-35131 Padova, Italy\n- 63 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, California, U.S.A.\n- 64 Jodrell Bank Centre for Astrophysics, Alan Turing Building, School of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13 9PL, U.K.\n- 65 Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA\n- 66 Kavli Institute for Cosmology Cambridge, Madingley Road, Cambridge, CB3 0HA, U.K.\n- 67 Kazan Federal University, 18 Kremlyovskaya St., Kazan, 420008, Russia\n- 68 LAL, Universit'e Paris-Sud, CNRS / IN2P3, Orsay, France\n- 69 LERMA, CNRS, Observatoire de Paris, 61 Avenue de l'Observatoire, Paris, France\n- 70 Laboratoire AIM, IRFU / Service d'Astrophysique - CEA / DSM CNRS - Universit'e Paris Diderot, Bˆat. 709, CEA-Saclay, F-91191 Gif-sur-Yvette Cedex, France\n- 71 Laboratoire Traitement et Communication de l'Information, CNRS (UMR 5141) and T'el'ecom ParisTech, 46 rue Barrault F-75634 Paris Cedex 13, France\n- 72 Laboratoire de Physique Subatomique et Cosmologie, Universit'e Grenoble-Alpes, CNRS / IN2P3, 53, rue des Martyrs, 38026 Grenoble Cedex, France\n- 73 Laboratoire de Physique Th'eorique, Universit'e Paris-Sud 11 & CNRS, Bˆatiment 210, 91405 Orsay, France  \n- 74 Lawrence Berkeley National Laboratory, Berkeley, California, U.S.A.\n- 75 Lebedev Physical Institute of the Russian Academy of Sciences, Astro Space Centre, 84 / 32 Profsoyuznaya st., Moscow, GSP-7, 117997, Russia\n- 76 Max-Planck-Institut fur Astrophysik, Karl-Schwarzschild-Str. 1, 85741 Garching, Germany\n- 77 National University of Ireland, Department of Experimental Physics, Maynooth, Co. Kildare, Ireland\n- 78 Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland\n- 79 Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark\n- 80 Optical Science Laboratory, University College London, Gower Street, London, U.K.\n- 81 SISSA, Astrophysics Sector, via Bonomea 265, 34136, Trieste, Italy\n- 82 SUPA, Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, U.K.\n- 83 School of Physics and Astronomy, Cardi GLYPH<11> University, Queens Buildings, The Parade, Cardi GLYPH<11> , CF24 3AA, U.K.\n- 84 Sorbonne Universit'e-UPMC, UMR7095, Institut d'Astrophysique de Paris, 98 bis Boulevard Arago, F-75014, Paris, France\n- 85 Space Sciences Laboratory, University of California, Berkeley, California, U.S.A.\n- 86 Special Astrophysical Observatory, Russian Academy of Sciences, Nizhnij Arkhyz, Zelenchukskiy region, Karachai-Cherkessian Republic, 369167, Russia\n- 87 Sub-Department of Astrophysics, University of Oxford, Keble Road, Oxford OX1 3RH, U.K.\n- 88 UPMCUnivParis 06, UMR7095, 98 bis Boulevard Arago, F-75014, Paris, France\n- 89 Universit'e de Toulouse, UPS-OMP, IRAP, F-31028 Toulouse cedex 4, France\n- 90 University of Granada, Departamento de F'ısica Te'orica y del Cosmos, Facultad de Ciencias, Granada, Spain\n- 91 University of Granada, Instituto Carlos I de F'ısica Te'orica y Computacional, Granada, Spain\n- 92 Warsaw University Observatory, Aleje Ujazdowskie 4, 00-478 Warszawa, Poland"}