We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

We present the properties of magnetic fields (B-fields) in two clumps (clump 1 and clump 2), located at the waist of the bipolar H II region Sh2-201, based on JCMT SCUBA-2/POL-2 observations of 850 $\mu$m polarized dust emission. We find that B-fields in the direction of the clumps are bent and compressed, showing bow-like morphologies, which we attribute to the feedback effect of the H II region on the surface of the clumps. Using the modified Davis-Chandrasekhar-Fermi method we estimate B-fields strengths of 266 $\mu$G and 65 $\mu$G for clump 1 and clump 2, respectively. From virial analyses and critical mass ratio estimates, we argue that clump 1 is gravitationally bound and could be undergoing collapse, whereas clump 2 is unbound and stable. We hypothesize that the interplay between thermal pressure imparted by the H II region, B-field morphologies, and the various internal pressures of the clumps (such as magnetic, turbulent, and gas thermal pressure), has the following consequences: (a) formation of clumps at the waist of the H II region; (b) progressive compression and enhancement of the B-fields in the clumps; (c) stronger B-fields will shield the clumps from erosion by the H II region and cause pressure equilibrium between the clumps and the H II region, thereby allowing expanding I-fronts to blow away from the filament ridge, forming bipolar H II regions; and (d) stronger B-fields and turbulence will be able to stabilize the clumps. A study of a larger sample of bipolar H II regions would help to determine whether our hypotheses are widely applicable.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) passed national acceptance and is taking pilot cycle of 'Shared-Risk' observations. The 19-beam receiver covering 1.05-1.45 GHz was used for most of these observations. The electronics gain fluctuation of the system is better than 1\% over 3.5 hours, enabling enough stability for observations. Pointing accuracy, aperture efficiency and system temperature are three key parameters of FAST. The measured standard deviation of pointing accuracy is 7.9$''$, which satisfies the initial design of FAST. When zenith angle is less than 26.4$^\circ$, the aperture efficiency and system temperature around 1.4 GHz are $\sim$ 0.63 and less than 24 K for central beam, respectively. The measured value of these two parameters are better than designed value of 0.6 and 25 K, respectively. The sensitivity and stability of the 19-beam backend are confirmed to satisfy expectation by spectral HI observations toward N672 and polarization observations toward 3C286. The performance allows FAST to take sensitive observations in various scientific goals, from studies of pulsar to galaxy evolution.

Magnetic fields are ubiquitous in our Universe, but remain poorly understood in many branches of astrophysics. A key tool for inferring astrophysical magnetic field properties is dust emission polarimetry. The James Clerk Maxwell Telescope (JCMT) is planning a new 850$\mu$m camera consisting of an array of 7272 paired Microwave Kinetic Inductance Detectors (MKIDs), which will inherently acquire linear polarization information. The camera will allow wide-area polarization mapping of dust emission at 14$^{\prime\prime}$-resolution, allowing magnetic field properties to be studied in a wide range of environments, including all stages of the star formation process, Asymptotic Giant Branch stellar envelopes and planetary nebula, external galaxies including starburst galaxies and analogues for the Milky Way, and the environments of active galactic nuclei (AGN). Time domain studies of AGN and protostellar polarization variability will also become practicable. Studies of the polarization properties of the interstellar medium will also allow detailed investigation of dust grain properties and physics. These investigations would benefit from a potential future upgrade adding 450$\mu$m capability to the camera, which would allow inference of spectral indices for polarized dust emission in a range of environments. The enhanced mapping speed and polarization capabilities of the new camera will transform the JCMT into a true submillimetre polarization survey instrument, offering the potential to revolutionize our understanding of magnetic fields in the cold Universe.

We report high angular resolution observations of linearly polarized dust emission towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A (hereafter, IRAS4A) using the Karl G. Jansky Very Large Array (JVLA) at K (11.5-16.7 mm), Ka (8.1-10.3 mm), and Q bands (6.3-7.9 mm), and using the Atacama Large Millimeter Array (ALMA) at Band 6 (1.2 mm) and Band 7 (0.85-0.89 mm). On 100-1000 AU scales, all of these observations consistently trace the hourglass shaped magnetic field topology as shown in the previous studies. In the innermost 100 AU region of IRAS4A1, the polarization position angles (E-field) detected at 6.3-16.7 mm are consistent, however, are nearly 90 degrees offset from those detected at 1.2 mm and 0.85-0.89 mm. Such a 90 degree offset may be explained by that the inner $\sim$100 AU area is optically thick at wavelengths shorter than $\sim$1.5 mm, whereby the observations probe the absorption of aligned dust against the weakly or unpolarized warm dust emission from the innermost region. This can also consistently explain why the highest angular resolution ALMA images at Band 7 show that the polarization percentage increases with dust brightness temperature in the inner $\sim$100 AU region of IRAS4A1. Following this interpretation and assuming that the dust grains are aligned with the magnetic fields, the inferred magnetic field position angle based on the 90$^{\circ}$ rotated at 6.3-7.9 mm in the central peak of IRAS4A1 is $\sim-22^{\circ}$, which is approximately consistent with the outflow direction $\sim-9^{\circ}$.

The HH 211 protostellar system is currently the youngest Class 0 system found with a rotating disk. We have mapped it at ~ 50 au (0.16") resolution, studying its magnetic field morphology with dust polarization in continuum at 232 and 358 GHz and its kinematics in C18O J=2-1 line. A flattened envelope extending out to ~ 400 au from the disk is detected in the continuum and C18O, slightly misaligned with the disk by 8 degree. It is spiraling inwards and expected to transform into a rotating disk at ~ 20 au, consistent with the disk radius estimated before. It appears to have a constant specific angular momentum and itself can result from an inside-out collapse of an extended envelope detected before in NH$_3$. In the flattened envelope, the polarization is mainly due to the magnetically aligned dust grains, inferring a highly pinched poloidal field morphology there. Thus, both the kinematics and field morphology support that the flattened envelope is a pseudodisk formed as the infalling gas is guided by the field lines to the equatorial plane. Interestingly, a point symmetric polarization distribution is also seen in the flattened envelope, implying that the pinched field lines also have a significant toroidal component generated by the rotation. No significant loss of angular momentum and thus no clear magnetic braking are detected in the flattened envelope around the disk probably because of the large misalignment between the axis of the rotation and the axis of the magnetic field in the cloud core.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We report Submillimeter Array molecular line observations in the 345 GHz band of five massive dense cores, Cyg-N38, Cyg-N43, Cyg-N48, Cyg-N51, and Cyg-N53 in the DR21 filament. The molecular line data reveal several dynamical features of the cores: (1) prominent outflows in all cores seen in the CO and SiO lines, (2) significant velocity gradients in Cyg-N43 and Cyg-N48 seen in the H13CN and H13CO+ lines suggesting 0.1-pc-scale rotational motions, and (3) possible infalls in Cyg-N48 found in the SiO and SO lines. Comparing the molecular line data and our dust polarization data in Ching et al. (2017), we find that the gradients of line-of-sight velocities appear to be randomly oriented relative to the plane-of-sky magnetic fields. Our simulations suggest that this random alignment implies parallel or random alignment between the velocity gradients and magnetic fields in the three dimensional space. The linewidths of H13CN emission are consistently wider than those of H13CO+ emission in the 3''-10'' detectable scales, which can be explained by the existence of ambipolar diffusion with maximum plane-of-the-sky magnetic field strengths of 1.9 mG and 5.1 mG in Cyg-N38 and Cyg-N48, respectively. Our results suggest that the gas dynamics may distort the magnetic fields of the cores of into complex structures and ambipolar diffusion could be important in dissipating the magnetic energies of the cores.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We performed the new JVLA full polarization observations at 40-48 GHz (6.3-7.5 mm) towards the nearby ($d$ $=$147$\pm$3.4 pc) Class 0 YSO IRAS 16293-2422, and compare with the previous SMA observations reported by Rao et al. (2009; 2014). We observed the quasar J1407+2827 which is weakly polarized and can be used as a leakage term calibrator for $<$9 GHz observations, to gauge the potential residual polarization leakage after calibration. We did not detect Stokes Q, U, and V intensities from the observations of J1407+2827, and constrain (3-$\sigma$) the residual polarization leakage after calibration to be $\lesssim$0.3\%. We detect linear polarization from one of the two binary components of our target source, IRAS\,16293-2422\,B. The derived polarization position angles from our observations are in excellent agreement with those detected from the previous observations of the SMA, implying that on the spatial scale we are probing ($\sim$50-1000 au), the physical mechanisms for polarizing the continuum emission do not vary significantly over the wavelength range of $\sim$0.88-7.5 mm. We hypothesize that the observed polarization position angles trace the magnetic field which converges from large scale to an approximately face-on rotating accretion flow. In this scenario, magnetic field is predominantly poloidal on $>$100 au scales, and becomes toroidal on smaller scales. However, this interpretation remains uncertain due to the high dust optical depths at the central region of IRAS\,16293-2422\,B and the uncertain temperature profile. We suggest that dust polarization at wavelengths comparable or longer than 7\,mm may still trace interstellar magnetic field. Future sensitive observations of dust polarization in the fully optically thin regime will have paramount importance for unambiguously resolving the magnetic field configuration.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

Polarized emission is detected in two young nearly edge-on protostellar disks in 343 GHz continuum at ~ 50 au (~ 0.12") resolution with Atacama Large Millimeter/submillimeter Array. One disk is in HH 212 (Class 0) and the other in HH 111 (early Class I) protostellar system. Polarization fraction is ~ 1%. The disk in HH 212 has a radius of ~ 60 au. The emission is mainly detected from the nearside of the disk. The polarization orientations are almost perpendicular to the disk major axis, consistent with either self-scattering or emission by grains aligned with a poloidal field around the outer edge of the disk because of optical depth effect and temperature gradient; the presence of a poloidal field would facilitate the launching of a disk wind, for which there is already tentative evidence in the same source. The disk of HH 111 VLA 1 has a larger radius of ~ 220 au and is thus more resolved. The polarization orientations are almost perpendicular to the disk major axis in the nearside, but more along the major axis in the farside, forming roughly half of an elliptical pattern there. It appears that toroidal and poloidal magnetic field may explain the polarization on the near and far side of the disk, respectively. However, it is also possible that the polarization is due to self-scattering. In addition, alignment of dust grains by radiation flux may play a role in the farside. Our observations reveal a diversity of disk polarization patterns that should be taken into account in future modeling efforts.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We present Submillimeter Array 880 $\mu$m dust polarization observations of six massive dense cores in the DR21 filament. The dust polarization shows complex magnetic field structures in the massive dense cores with sizes of 0.1 pc, in contrast to the ordered magnetic fields of the parsec-scale filament. The major axes of the massive dense cores appear to be aligned either parallel or perpendicular to the magnetic fields of the filament, indicating that the parsec-scale magnetic fields play an important role in the formation of the massive dense cores. However, the correlation between the major axes of the cores and the magnetic fields of the cores is less significant, suggesting that during the core formation, the magnetic fields below 0.1 pc scales become less important than the magnetic fields above 0.1 pc scales in supporting a core against gravity. Our analysis of the angular dispersion functions of the observed polarization segments yields the plane-of-sky magnetic field strengths of 0.4--1.7 mG of the massive dense cores. We estimate the kinematic, magnetic, and gravitational virial parameters of the filament and the cores. The virial parameters show that in the filament, the gravitational energy is dominant over magnetic and kinematic energies, while in the cores, the kinematic energy is dominant. Our work suggests that although magnetic fields may play an important role in a collapsing filament, the kinematics arising from gravitational collapse must become more important than magnetic fields during the evolution from filaments to massive dense cores.

We present Submillimeter Array (SMA) polarization observations of the CO $J$ = 3--2 line toward the NGC1333 IRAS 4A. The CO Stokes $I$ maps at an angular resolution of $\sim$1$\arcsec$ reveal two bipolar outflows from the binary sources of the NGC 1333 IRAS 4A. The kinematic features of the CO emission can be modeled by wind-driven outflows at $\sim$ 20$\arcdeg$ inclined from the plane of the sky. Close to the protostars the CO polarization, at an angular resolution of $\sim$$2\farcs3$, has a position angle approximately parallel to the magnetic field direction inferred from the dust polarizations. The CO polarization direction appears to vary smoothly from an hourglass field around the core to an arc-like morphology wrapping around the outflow, suggesting a helical structure of magnetic fields that inherits the poloidal fields at the launching point and consists of toroidal fields at a farther distance of outflow. The helical magnetic field is consistent with the theoretical expectations for launching and collimating outflows from a magnetized rotating disk. Considering that the CO polarized emission is mainly contributed from the low-velocity and low-resolution data, the helical magnetic field is likely a product of the wind-envelope interaction in the wind-driven outflows. The CO data reveal a PA of $\sim$ 30$\arcdeg$ deflection in the outflows. The variation in the CO polarization angle seems to correlate with the deflections. We speculate that the helical magnetic field contributes to $\sim$ 10$\arcdeg$ deflection of the outflows by means of Lorenz force.

We highlight distinct and systematic observational features of magnetic field morphologies in polarized submm dust continuum. We illustrate this with specific examples and show statistical trends from a sample of 50 star-forming regions.

Submillimeter dust polarization measurements of a sample of 50 star-forming regions, observed with the SMA and the CSO covering pc-scale clouds to mpc-scale cores, are analyzed in order to quantify the magnetic field importance. The magnetic field misalignment $\delta$ -- the local angle between magnetic field and dust emission gradient -- is found to be a prime observable, revealing distinct distributions for sources where the magnetic field is preferentially aligned with or perpendicular to the source minor axis. Source-averaged misalignment angles $\langle|\delta|\rangle$ fall into systematically different ranges, reflecting the different source-magnetic field configurations. Possible bimodal $\langle|\delta|\rangle$-distributions are found for the separate SMA and CSO samples. Combining both samples broadens the distribution with a wide maximum peak at small $\langle|\delta|\rangle$-values. Assuming the 50 sources to be representative, the prevailing source-magnetic field configuration is one that statistically prefers small magnetic field misalignments $|\delta|$. When interpreting $|\delta|$ together with an MHD force equation, as developed in the framework of the polarization-intensity gradient method, a sample-based scaling fits the magnetic field tension-to-gravity force ratio $\langle\Sigma_B\rangle$ versus $\langle|\delta|\rangle$ with $\langle\Sigma_B\rangle = 0.116 \cdot \exp(0.047\cdot \langle|\delta|\rangle)\pm 0.20$ (mean error), providing a way to estimate the relative importance of the magnetic field, only based on measurable field misalignments $|\delta|$. The force ratio $\Sigma_B$ discriminates systems that are collapsible on average ($\langle \Sigma_B\rangle <1$) from other molecular clouds where the magnetic field still provides enough resistance against gravitational collapse ($\langle \Sigma_B\rangle >1$) (abridged).

HH 211 is a young Class 0 protostellar system, with a flattened envelope, a possible rotating disk, and a collimated jet. We have mapped it with the Submillimeter Array in 341.6 GHz continuum and SiO J=8-7 at ~ 0.6 resolution. The continuum traces the thermal dust emission in the flattened envelope and the possible disk. Linear polarization is detected in the continuum in the flattened envelope. The field lines implied from the polarization have different orientations, but they are not incompatible with current gravitational collapse models, which predict different orientation depending on the region/distance. Also, we might have detected for the first time polarized SiO line emission in the jet due to the Goldreich-Kylafis effect. Observations at higher sensitivity are needed to determine the field morphology in the jet.

Massive stars ($M > 8$ \msun) typically form in parsec-scale molecular clumps that collapse and fragment, leading to the birth of a cluster of stellar objects. We investigate the role of magnetic fields in this process through dust polarization at 870 $\mu$m obtained with the Submillimeter Array (SMA). The SMA observations reveal polarization at scales of $\lsim$ 0.1 pc. The polarization pattern in these objects ranges from ordered hour-glass configurations to more chaotic distributions. By comparing the SMA data with the single dish data at parsec scales, we found that magnetic fields at dense core scales are either aligned within $40^\circ$ of or perpendicular to the parsec-scale magnetic fields. This finding indicates that magnetic fields play an important role during the collapse and fragmentation of massive molecular clumps and the formation of dense cores. We further compare magnetic fields in dense cores with the major axis of molecular outflows. Despite a limited number of outflows, we found that the outflow axis appears to be randomly oriented with respect to the magnetic field in the core. This result suggests that at the scale of accretion disks ($\lsim 10^3$ AU), angular momentum and dynamic interactions possibly due to close binary or multiple systems dominate over magnetic fields. With this unprecedentedly large sample massive clumps, we argue on a statistical basis that magnetic fields play an important role during the formation of dense cores at spatial scale of 0.01 - 0.1 pc in the context of massive star and cluster star formation.